1. 辐照损伤分子动力学模拟工程导论

核反应堆结构材料在服役期间持续承受中子辐照，高能中子与晶格原子发生弹性或非弹性碰撞，将原子撞离平衡位置，形成离位原子和空位。这些初级缺陷在皮秒尺度内迁移、复合、聚集，最终演化成位错环、空洞、堆垛层错四面体等微观缺陷团簇，宏观上表现为硬化、脆化、肿胀和蠕变加速。理解这一从埃到毫米、从飞秒到年的跨尺度过程，是核材料工程的核心难题，也是本课程要建立的工程能力。

课程定位：从物理模型到工程判断

本课程不是分子动力学（MD）软件操作手册，也不是纯理论固体物理课。它面向的是"用MD解决核材料工程问题"的完整能力链：把工程失效现象翻译成可计算的物理问题，选择合适的势函数和边界条件，设计可复现的模拟方案，从轨迹中提取有物理意义的缺陷统计量，并判断结果是否可信、能否外推到工程条件。学完本课程，你应当能够独立完成一次从问题定义到结果判读的辐照损伤MD模拟任务，并清楚它的适用边界。

辐照损伤多尺度模拟框架

单一尺度的模拟无法覆盖辐照损伤的全貌，工程上采用多尺度接力框架。下图给出各尺度方法的分工与信息传递关系。

图中信息流的方向是：DFT 计算缺陷形成能和迁移势垒，用于拟合或验证 MD 势函数；MD 模拟级联损伤，输出缺陷类型、数量、空间分布和团簇尺寸，作为 KMC 的初始条件和反应参数；KMC 在扩散时间尺度上模拟微结构演化，为连续介质模型提供肿胀率、硬化系数等本构关系。实验数据（离子辐照、中子辐照、TEM 观察）则反向校验各尺度模型的可靠性。MD 处于链条中枢，向上承接电子尺度参数，向下输出介观演化规律。

分子动力学在核材料工程中的角色

MD 在核材料工程中承担三类具体任务。第一类是机制揭示：例如模拟 10 keV 初级撞出原子（PKA）在 α-Fe 中产生的级联损伤，观察热峰持续时间、离位原子峰值数量、缺陷复合比例，回答"辐照初期缺陷如何产生和湮灭"。第二类是参数提取：通过 NEB 方法计算空位和间隙原子的迁移能垒，为 KMC 提供输入；通过分子静力学计算缺陷团簇的结合能，判断其热稳定性。第三类是性能预测：模拟辐照后晶格的热导率变化、裂纹扩展行为、氢氦泡形核过程，为工程寿命评估提供原子尺度依据。

需要清醒认识 MD 的局限：时间尺度通常限于纳秒到微秒，空间尺度限于百万原子以内，势函数的精度直接决定结果可信度。因此 MD 结果不能直接等同于工程数据，必须经过尺度外推和实验验证。

典型工程问题与模拟目标

以反应堆压力容器钢（RPV 钢）中子辐照脆化为例，工程问题是：随着辐照剂量增加，材料韧脆转变温度升高，威胁压力容器完整性。MD 模拟目标是：模拟富铜团簇和 Mn-Ni-Si 团簇在级联损伤区的形核过程，统计团簇尺寸分布和数密度，判断其与位错的交互作用强度。模拟中需要设置的关键参数包括：PKA 能量（典型 5~20 keV）、辐照温度（约 290°C，即 563 K）、势函数（如 Fe-Cu-Ni 三元势）、边界条件（三维周期性）、系综（NVT 或 NVE）、时间步长（1~2 fs）。

模拟结果的判读依据包括：级联损伤产生的 Frenkel 对数量是否随 PKA 能量线性增长；缺陷团簇中空位型和间隙型比例是否合理；富铜团簇的临界形核尺寸是否与实验观察一致。若模拟得到的团簇数密度比实验高一个数量级，通常提示势函数中 Cu-Cu 结合能偏强或模拟时间不足，需要调整势函数或延长退火时间。

本章建立了课程的整体坐标系：MD 是辐照损伤多尺度模拟框架中的原子尺度核心环节，向上依赖 DFT 参数，向下为 KMC 和连续介质模型提供缺陷演化数据。后续所有章节的模拟任务，都可以在这个框架中找到位置和意义。