4. 原子尺度实战

当高能粒子（中子、离子或电子）进入金属晶格时，会在皮秒到纳秒的时间尺度、纳米到数十纳米的空间尺度上引发一连串原子级事件：级联碰撞、离位损伤、热峰形成与淬灭、缺陷团簇的迁移与演化。分子动力学（MD）正是捕捉这一过程最直接的工具——它把原子当作经典粒子，用经验势函数描述原子间相互作用，通过数值积分牛顿方程追踪每个原子的轨迹，从而在原子尺度上"看见"辐照损伤的完整演化。

一、物理图像：从级联碰撞到热峰

一次辐照事件通常经历几个阶段。初始阶段，初级撞出原子（PKA）获得 keV 量级动能，在晶格中引发级联碰撞，产生大量弗仑克尔对（空位-间隙原子对）。级联核心区域原子剧烈扰动，局部动能极高，等效温度可达数千开尔文，形成所谓"热峰"。热峰在几皮秒内通过向周围晶格导热而淬灭，留下一个由空位团簇和间隙原子团簇组成的缺陷核心。随后，缺陷在长时间尺度上迁移、复合或聚集，最终形成位错环、空洞等可观测的损伤组织。

MD 模拟的关键在于：势函数要能正确描述远离平衡态的原子相互作用，时间步长要足够小以捕捉高能碰撞，而模拟盒子要足够大以避免边界效应。常用势函数包括嵌入原子法（EAM）及其改进版本，对 Fe、W、Cu 等金属已有较成熟的参数化。

二、模拟流程与关键参数

以纯铁（BCC Fe）为例，一个典型的辐照损伤 MD 模拟可按以下步骤组织。首先构建含约 10⁵～10⁶ 个原子的 BCC 超胞，尺寸约 20×20×20 nm³，三个方向均采用周期性边界条件。在 NPT 或 NVT 系综下于 300 K 弛豫 20 ps，使体系达到热平衡。随后在盒子中心选取一个原子作为 PKA，赋予其 1～10 keV 的动能，方向通常取低指数晶向（如 <135> 或 <100>），以研究不同晶向的级联差异。

积分采用 Velocity-Verlet 算法，时间步长取 0.1～1 fs。高能碰撞阶段用 0.1 fs 保证能量守恒，热峰淬灭后可放宽到 1 fs 以提高效率。模拟总时长通常为 20～50 ps，足以覆盖级联产生、热峰淬灭和初级缺陷复合。温度控制采用 Nose-Hoover 热浴，但需注意：在级联核心区域应避免强耦合热浴，否则会人为抽走能量、抑制热峰。

# LAMMPS 输入片段（示例参数，非实测结果）
units           metal
atom_style      atomic
boundary        p p p
pair_style      eam/alloy
pair_coeff      * * Fe_Mendelev.eam.alloy Fe
timestep        0.001          # 1 fs，级联阶段可改为 0.0001
velocity        all create 300.0 12345 mom yes rot yes
fix             1 all nvt temp 300.0 300.0 0.1
run             20000          # 20 ps 弛豫
# 选取中心原子作为 PKA，赋予 5 keV 动能
variable        pka equal 5000
# ...（此处省略 PKA 初始化与方向设置）
run             30000          # 30 ps 级联演化

三、缺陷识别与结果分析

模拟结束后，核心任务是识别缺陷。最常用的方法是 Wigner-Seitz 胞分析：将参考完美晶格的每个原子位置作为 WS 胞中心，统计每个胞内原子数。原子数为 0 的胞对应空位，原子数大于 1 的胞对应间隙原子或挤列子。通过统计空位和间隙的数量、空间分布及团簇尺寸，可定量描述损伤程度。

另一个重要指标是缺陷团簇的尺寸分布。级联结束后，空位往往聚集成三维空位团簇，间隙原子则形成位错环或间隙团簇。团簇分析通常采用连通性判据：若两个缺陷原子间距小于某一截断半径（如 0.3a₀），则归为同一团簇。统计团簇中缺陷原子数，即可得到尺寸分布直方图。

热峰效应可通过动能或温度的空间分布来观察。将模拟盒子划分为若干薄层，统计每层原子的平均动能，可得到热峰沿深度方向的衰减曲线。典型结果显示：热峰中心温度在级联初期迅速升高，随后在几皮秒内指数衰减，最终回到环境温度。热峰持续时间与 PKA 能量、材料热导率密切相关。

四、从轨迹到物理结论

MD 模拟的价值不仅在于"看到"缺陷，更在于建立工艺参数与损伤结果之间的定量关系。例如，改变 PKA 能量，统计残余弗仑克尔对数量，可拟合出损伤效率随能量的变化曲线；改变温度，可研究热峰淬灭速率对缺陷复合的影响；改变晶向，可比较不同方向的级联形貌差异。这些关系为更高尺度的速率理论或团簇动力学模型提供输入参数。

需要注意的是，MD 的时间尺度（皮秒到纳秒）和空间尺度（纳米级）远小于实际辐照条件，因此模拟结果通常用于揭示机制、提取参数，而非直接预测宏观寿命。将 MD 得到的缺陷产生率、迁移能垒等参数传递给动力学蒙特卡洛（KMC）或团簇动力学（CD）模型，才能实现从原子尺度到介观尺度的跨尺度衔接。

五、实践中的常见问题与判断

MD 辐照模拟有几个容易出错的环节。第一，势函数选择不当会导致级联形貌失真，例如某些 EAM 势在高温下会高估或低估热峰温度，建议先用已知实验或第一性原理数据做验证。第二，时间步长过大可能引起能量不守恒，表现为总能量漂移，可通过监控总能量曲线判断。第三，边界条件若采用非周期性，级联产生的缺陷可能被表面吸收，导致残余缺陷数偏低。第四，热浴耦合过强会人为抑制热峰，建议在级联阶段使用弱耦合或仅对边界层控温。

结果分析时，除了统计缺陷数量，还应关注缺陷的空间关联性。例如，空位团簇与间隙团簇是否空间分离，是否形成位错环结构。这些信息可通过可视化工具（如 OVITO）结合共同近邻分析（CNA）或配位分析获得。最终，将 MD 得到的缺陷产生率、团簇尺寸分布等参数整理成表格，作为 KMC 或 CD 模型的输入，完成从原子尺度到介观尺度的参数传递。