2. 辐照损伤物理基础

辐照损伤的分子动力学模拟，本质上是在原子尺度上重现“一个高能粒子撞进晶体后发生了什么”。要读懂模拟输出，必须先建立一条从入射粒子到宏观性能退化的物理链条：级联碰撞产生点缺陷，点缺陷聚集形成团簇，团簇与气泡进一步改变材料的力学与尺寸行为。本节按这条链条梳理关键概念与可量化的判据。

级联碰撞：从一次离位到热峰

当入射粒子（中子、离子或高能反冲原子）把足够能量传给晶格原子时，该原子离开格点成为离位原子，并继续撞击其他原子，形成分支状的碰撞级联。级联可分为三个阶段：弹道阶段（皮秒以内，原子以直线或近直线运动，产生大量离位）、热峰阶段（局部区域原子剧烈无序化，温度等效可达数千开尔文，持续约1–10 ps）、退火/复合阶段（热峰冷却，绝大多数离位原子回到格点，只有少部分以稳定缺陷形式残留）。

分子动力学模拟级联时，典型做法是：先对完美晶体进行能量最小化和一定时间的弛豫，再给初级撞出原子（PKA）赋予初始动能和方向，随后在微正则或恒温系综下演化。PKA能量常用1–50 keV量级，模拟盒子从几万到几百万原子不等。判断级联是否“成功”的关键不是看瞬时离位数，而是看退火结束后残留的稳定缺陷数。

离位阈能与NRT模型

离位阈能（E_d）是把一个原子从格点搬到远离格点、不再自发复合所需的最小能量。对金属，E_d通常在20–40 eV；对半导体如Si，方向依赖性很强，常用值约15–20 eV。E_d不是单一常数，它随晶向、温度、表面距离变化，因此模拟中应明确所采用的E_d定义与判据。

工程上估算离位损伤量常用 NRT模型（Norgett–Robinson–Torrens）。其核心是：先由能量分配得到损伤能 E_damage，再除以2E_d得到离位数：

N_d = κ · E_damage / (2 · E_d)
κ ≈ 0.8   （NRT标准效率因子）
E_damage = E_PKA − E_loss（电子能损等）

例如，一个10 keV的PKA在Fe中，若E_d=40 eV、E_damage≈8 keV，则NRT估算N_d≈0.8×8000/(2×40)=80个离位原子。注意：NRT给出的是“离位事件数”，不是最终稳定缺陷数。MD模拟通常显示，稳定Frenkel对数量只有NRT值的10%–30%，这个比值称为存活率，是连接理论估算与模拟结果的重要参数。

点缺陷与缺陷团簇

级联退火后残留的缺陷主要有两类：空位（vacancy）和自间隙原子（SIA）。单个空位和单个SIA是点缺陷；当它们聚集，就形成缺陷团簇。SIA团簇通常呈小尺寸、可移动的位错环（如1/2<111>或<100>环），空位团簇则多为三维空位簇或堆垛层错四面体。

在MD中识别缺陷常用两种方法：一是Wigner–Seitz缺陷分析，把每个格点归属到最近原子，统计空位与SIA；二是共同邻域分析（CNA）或配位数，识别非晶或位错结构。输出量包括：缺陷总数、团簇尺寸分布、团簇类型、空间分布。这些量直接决定后续硬化、肿胀的预测。

辐照硬化与脆化

缺陷团簇对位错运动构成障碍。位错线在滑移面上遇到不可穿透的团簇时，需要额外应力才能绕过或切过，宏观表现为屈服强度升高，即辐照硬化。硬化量常用分散障碍强化模型估算：Δσ ∝ √(N·d)，其中N是障碍数密度，d是障碍平均尺寸。MD中可通过计算位错与团簇的相互作用能、临界剪切应力来量化这一效应。

硬化的同时往往伴随脆化：材料强度上升但塑性下降，韧脆转变温度（DBTT）升高。对反应堆压力容器钢，辐照后DBTT可升高数十到上百摄氏度。MD层面，脆化表现为位错形核困难、裂纹扩展阻力降低、孪生或解理倾向增强。模拟时可通过拉伸/剪切加载，观察应力–应变曲线中屈服后软化、裂纹萌生位置与扩展路径。

气泡与肿胀的物理图像

当材料中存在氦、氢等不溶性气体时，气体原子与空位结合形成气泡。气泡内压随气体原子数增加而升高，当内压超过基体屈服强度时，气泡周围发射位错环，气泡继续长大。气泡与空位团簇共同导致材料体积膨胀，即肿胀。肿胀率通常用 ΔV/V 表示，典型快堆材料在几十dpa下可达百分之几到十几。

MD模拟气泡时，常用方法是在空位簇中插入He原子，然后弛豫并升温，观察气泡是否稳定、是否发射位错环、内压与体积的关系。关键参数包括：He/V比、气泡半径、内压、发射位错环的临界条件。这些结果可与速率理论或团簇动力学模型耦合，预测长期肿胀行为。