Vacancy Diffusion Across FeCrAl Alloy Composition Space for Accident-Tolerant Fuel Cladding

技术摘要：FeCrAl 合金成分空间内的空位扩散研究，用于事故容错燃料包 cladding

问题陈述
铁-铬-铝 (FeCrAl) 合金因其相比于锆合金具有更优异的高温抗氧化性能，成为轻水反应堆 (LWR) 中事故容错燃料 (ATF) 包 cladding 的领先候选材料。然而，这些包 cladding 在反应堆内的寿命受辐照诱导微观结构演化的支配，其中点缺陷传输——特别是空位扩散——是主要的机制。虽然保护性氧化层的形成已得到充分表征，但 FeCrAl 在辐照下的基础缺陷物理学仍不为人所知。现有的理论处理通常依赖于理想化或成分简化的环境，未能考虑到高浓度固溶体合金中固有的复杂局部化学波动。这些局部环境显著改变了空位结合、迁移路径以及诸如团簇和捕获等集体行为。目前迫切需要厘清各组成元素（Fe、Cr、Al）对整体空位传输的个体贡献，以预测在热与辐照复合载荷下的长期材料演化、相稳定性及蠕变。

方法论
作者采用了一种成分解析的动力学蒙特卡洛 (KMC) 框架，来模拟在固定 5 at.% Al 浓度下，跨越系统性 Fe/Cr 成分范围（从富 Fe 到富 Cr）的空位扩散。该方法论通过以下步骤进行：

1. 势函数与数据生成： 利用拟合至密度泛函理论 (DFT) 数据库的嵌入原子法 (EAM) 势函数来评估原子间相互作用。使用“跳跃装饰”（Hop-Decorate）工作流生成空位迁移数据，其中将 BCC Fe 参考结构根据目标全局成分随机装饰上 Fe、Cr 和 Al 原子。通过爬升图像能垒搜索法 (CI-NEB) 计算全球成分为 $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$、$\text{Fe}_{47.5}\text{Cr}_{47.5}\text{Al}_5$ 和 $\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$ 的最小能量路径。
2. 代理模型： 为了避免在 KMC 过程中进行高昂的即时 NEB 计算，研究基于生成的能垒数据库训练了线性代理模型。每个跳跃的局部化学环境通过描述符进行编码，这些描述符统计了鞍点和迁移原子周围固定径向壳层内的 Fe、Cr 和 Al 原子数量。针对每种迁移物种（Fe、Cr、Al）分别拟合了独立的线性回归模型，以捕捉物种特异性的能垒依赖性。
3. KMC 模拟： 训练好的代理模型驱动在 500 K 至 1500 K 温度范围内的驻留时间 KMC 模拟。模拟追踪单个空位的轨迹超过 $10^5$ 步，以计算均方位移 (MSD) 并提取扩散系数 ($D$)、激活能 ($E_a$) 和指前因子 ($D_0$)。

关键结果

• 迁移能垒与局部环境： 研究证实，迁移能垒对局部化学环境高度敏感。虽然纯 Fe 的能垒约为 0.63 eV，Cr 约为 0.95 eV，Al 约为 0.69 eV；但在合金中，能垒分布显著展宽，尤其是 Cr 和 Al，表明局部化学不对称性调节了迁移路径。
• 铝的作用： 尽管 Al 的浓度很低 (5 at.%)，但它对迁移能垒产生了不成比例的影响。在所有三个代理模型中，第一最近邻 (1NN) 鞍点壳层中存在 Al 原子会显著提高能量能垒（对于 Fe 为 +0.048 eV，对于 Cr 为 +0.073 eV，对于 Al 为 +0.151 eV）。
• 扩散率的成分依赖性：
  • 富 Fe 到等原子比： 在 $\text{Fe}_{80}\text{Cr}_{15}\text{Al}_5$ 到 $\text{Fe}_{32.5}\text{Cr}_{62.5}\text{Al}_5$ 的成分范围内，空位扩散率对全局成分并不敏感。激活能收敛于 0.70–0.84 eV 的范围，表明在此类机制中 Fe 迁移主导了传输。
  • 富 Cr 机制： 富 Cr 成分 ($\text{Fe}_{15}\text{Cr}_{80}\text{Al}_5$) 表现出显著升高的激活能，达到 1.135 eV。因此，在较低温度下（例如在 ~570 K 时，富 Fe 合金约为 ~$10^{-13}$ m$^2$ s$^{-1}$，而富 Cr 合金仅为 ~$10^{-17}$ m$^2$ s$^{-1}$），其宏观扩散率比富 Fe 合金低几个数量级。
• 渗透与“迟缓”扩散： 结果表明存在一个渗透阈值，即随着 Cr 含量的增加，低能垒路径（由 Fe 跳跃主导）的连接网络发生崩溃。在富 Cr 机制中，空位迁移由高能垒 Cr 跳跃主导，导致“迟缓扩散”。作者观察到，在化学粗糙的景观中，空位轨迹涉及从捕获口袋中的逃逸，导致在达到扩散机制之前出现非线性的均方位移行为。

意义与主张
本文声称，FeCrAl 合金的宏观扩散率不能通过简单的纯端元能垒的成分加权平均来预测；相反，它是由局部能量景观的异质性和低能垒迁移路径的具体连通性决定的。

• 设计启示： 研究确定了一个明显的成分机制（富 Cr 机制），在该机制下空位迁移率受到抑制。这种抑制意味着在辐照下形成的富 Cr $\alpha'$ 相中，空位迁移性较低，这可能降低空位-间隙子复合的速率并改变缺陷积累动力学。
• 事故容错： 在温度升高至 1200 K 以上的失水事故 (LOCA) 条件下，扩散率的成分依赖性变得非常显著。作者认为，具有高激活能和低扩散率的富 Cr 成分，可能通过提供有利的空位动力学特性来补充其已知的抗氧化性能，从而提升事故容错包 cladding 的性能。
• 方法论贡献： 本工作证明了通过以物种解析的线性代理模型（基于 NEB 数据库参数化）为基础的降阶动力学蒙特卡洛框架，可以有效地绘制复杂成分空间内的缺陷传输图谱。这种方法在原子过程与用于燃料性能计算的介观模型之间建立了机械联系，超越了合金设计的经验试错法。

作者注明的局限性
作者明确指出了一些局限性：模型将晶格视为静态且冻结的，忽略了局部弛豫和随时间变化的辐照诱导成分偏析（segregation）。线性代理模型无法捕捉方向性或对称性破缺效应，也无法捕捉能量景观中的磁性贡献。此外，本研究仅关注空位传输，并未涉及间隙子团簇传输或位错环的辐照诱导转变，而这些因素也可能影响缺陷汇密度。