机器学习原子间势揭示α-Fe一般晶界氢脆起源

机器学习原子间势揭示α-Fe一般晶界氢脆起源

1  研究概述

氢脆是制约高强结构材料在绿氢制备、输运与利用全链条中应用的关键瓶颈，尤其沿一般晶界的氢致开裂因其缺乏晶格对称性而更难预测与防范。传统经验势（如EAM）在描述含氢晶界的变形断裂行为时存在显著精度不足，而密度泛函理论（DFT）虽精度高却受限于计算规模，无法直接模拟含氢晶界的大尺度力学响应。机器学习原子间势（MLIP）为弥合这一精度-尺度的鸿沟提供了新途径，但现有Fe-H MLIP的训练集多局限于对称倾转晶界，难以准确描述一般晶界的氢偏聚与力学行为。

针对上述问题，Kazuma Ito、Takashi Otaki、Tatsuya Yokoi、Katsutoshi Hyodo等人在“Communications Materials"期刊上发表了题为“Machine Learning Interatomic Potential Reveals Hydrogen Embrittlement Origins at General Grain Boundaries in α-Iron"的研究论文。该研究通过并行学习（concurrent learning）策略，系统构建了覆盖广义层错、表面、对称倾转晶界、位错及空位团簇等所有与氢致脆化相关原子环境的Fe-H二元MLIP训练集（图1）。基于Behler-Parrinello神经网络势（BNNP）框架，最终实现了对α-Fe中一般晶界氢偏聚、<111>/2全位错发射、变形孪生及裂纹萌生扩展的DFT精度描述。大规模分子动力学模拟揭示了氢在不同变形机制晶界中的差异化影响，为高强合金抗氢脆设计提供了原子尺度的力学判据。

图1  Fe-H二元体系MLIP训练数据集构建流程：(a) DFT初始数据集；(b) 四个并行深度势（DP）训练；(c) 多条件分子动力学候选构型采集；(d) 候选构型DFT标记与迭代

2  MLIP构建方法与精度验证

2.1  并行学习策略与训练集

训练集构建经历169轮并行学习迭代，最终包含36,438个构型（约155万局部原子环境）。训练集系统涵盖了：纯α-Fe与应变Fe中的H溶解能、H-空位团簇、相邻四面体间隙位的双氢对、广义层错能、七种对称倾转晶界的氢偏聚能、以及纳米多晶中一般晶界的结构与偏聚评估等。采用节点并行训练与修正的n2p2超参数（截断半径从6.5 Å 优化至6.0 Å），在保证精度的前提下显著提升计算效率。

2.2  基本物理性质与缺陷-氢相互作用精度

图2展示了构建MLIP对Fe-H体系关键物理量DFT结果的复现能力。氢在四面体间隙位（TIS）的溶解能可准确匹配单轴应变与静水压条件的变化（图2a,b）；螺位错核心对氢的偏聚能预测与DFT高度一致（图2c），尽管训练集中未明确包含螺位错-氢构型，验证了MLIP优异的可迁移性；广义层错能曲线与氢-空位团簇相互作用能同样达到DFT精度（图2d,e）。全训练集总能量与原子力的均方根误差分别为0.13 meV/atom和35.8 meV/Å。

图2  构建MLIP的基本物性精度验证：(a-b) TIS氢溶解能；(c) 螺位错氢偏聚能；(d) (112)广义层错能（氢位于TIS）；(e) 空位-氢团簇相互作用能；(f) 折射率表面分离能；(g)能量的键长依赖性；(h) 四面体位置原子相互作用能；

2.3  对称倾转晶界氢偏聚精度

图3展示了七种对称倾转晶界（Σ3-Σ11）的氢偏聚能预测结果。MLIP整体准确复现了DFT的偏聚能趋势与相对次序，仅在高Σ值晶界的极低偏聚能区出现轻微系统偏差，但偏差幅度在可接受范围内。这表明MLIP能够精确描述氢在不同晶界能量/结构下的偏聚行为差异。

图3  MLIP对七种对称倾转晶界氢偏聚能的预测精度：(a-g) 分别为Σ3(112)、Σ3(111)、Σ5(210)、Σ5(310)、Σ9(114)、Σ9(221)、Σ11(113)晶界；(h) DFT和原子间势在七个晶界位点的偏析能对比。

2.4  一般晶界能量与氢偏聚精度

针对最难描述的一般晶界，研究从含8晶粒的纳米多晶（28.8 nm）中提取计算单元，评估MLIP对一般晶界氢偏聚能的预测能力（图4）。MLIP预测的平均晶界能与实验外推值（图4a）基本吻合，远优于EAM势。一般晶界氢偏聚能的RMSE仅为0.042–0.103 eV（图4b,c），表明即使在训练集中未出现相关构型，MLIP仍对一般晶界氢偏聚展现出优异的外推精度。这一结果为后续大规模MD模拟提供了必要的置信度。

图4  MLIP对一般晶界氢偏聚的预测精度：(a) 平均晶界能与模型尺寸关系；(b) 纳米多晶中一般晶界氢偏聚能的评估方法；(c-e) MLIP与DFT氢偏聚能对比

3  一般晶界氢偏聚行为

基于构建的MLIP，研究建立了六个随机取向的双晶一般晶界模型（约30万原子），考察不同体氢浓度条件下的晶界氢偏聚状态（图5）。Voronoi分析表明，氢优先偏聚于晶界附近体积膨胀的八面体间隙位（图5a,b），晶界氢占据数与体浓度呈近似正比关系（图5c）。在中等与高体浓度条件下，一般晶界的氢面密度分别达到1.4和6.7 atoms/nm²，与Sato等人在高纯Fe中阴极充氢实验测得的2.3–2.9 atoms/nm²处于同一数量级，验证了模拟结果的合理性。

图5  一般晶界氢偏聚状态：(a) 双晶模型示例（GGB 1）；(b) 原子能量着色图；(c) 体氢浓度与晶界氢浓度的关系；(d) GGB1晶界附近的氢浓度分布；(e-f) 氢在GGB1晶界中心的一个原子层（2.5 nm）区域的Voronoi体积和配位数的概率密度；(g) 体积氢浓度为0.0044 at.%时，GGB1晶界中心附近的原子结构

4  氢对一般晶界变形断裂行为的影响

4.1  两种变形模式

单轴拉伸MD模拟揭示了六个一般晶界模型的两种本征变形机制（图6a）：GGB 1–3由晶界处的<111>/2全位错发射主导塑性；GGB 4–6由晶界附近连续发射<111>/6不全位错形成变形孪晶主导塑性。前者的最大应力（11.2–13.0 GPa）通常略高于后者（9.3–9.8 GPa）。与Σ3(111)（15.4 GPa）和Σ3(112)（13.3 GPa）晶界相比，低Σ值晶界的最大应力显著更高，体现了晶界有序度对应力响应的重要影响。

图6  氢偏聚对一般晶界变形与断裂行为的影响：(a) 清洁晶界应力-应变曲线；(b-c) GGB 1和GGB 3受氢影响的应力-应变曲线；(d-f) 氢对最大应力与断裂应变的影响；(g) GGB 3晶界迁移与位错发射；(h) GGB 4变形孪晶形成；(i) GGB 4-6断裂应变随氢浓度的变化；(j) GGB 4最大应力附近应变处的原子结构；(k) GGB4晶界变形孪晶的形核机制

4.2  氢对位错主导型晶界的影响：抑制位错发射、促进晶界断裂

对于位错主导塑性的GGB 1–3模型（图6b–g），氢偏聚呈现一致趋势：增大最大应力（抑制位错发射的宏观表现），并降低晶界断裂应变。以GGB 1为例，中等氢浓度下最大应力显著增大，断裂应变显著降低（图6b）；位错分析（DXA）表明，位错发射位置不受氢影响，但晶界内位错线伸入晶内即触发发射，氢在位错线附近偏聚提升了位错发射的能垒（图6e）。在GGB 2的中等氢浓度条件下，最大应力后并未立即开裂，而是继续以位错发射维持变形，约在14%应变达到第二应力峰后才发生断裂。在GGB 3的低/中氢浓度条件下，晶界从约6%应变起发生显著迁移，由迁移形成的突出部位发射位错，但氢偏聚使应力单调下降并最终断裂。高氢浓度时晶界迁移被抑制，弹性应变能无法释放，导致在位错发射之前即发生低应力脆断。

4.3  氢对孪生主导型晶界的影响：弱化显著降低

GGB 4–6模型由变形孪生负责塑性，其氢敏感性显著弱于位错主导型晶界（图6h,i）。氢偏聚几乎不改变最大应力与断裂应变，即使晶界氢浓度与GGB 1–3相当（甚至更高）。变形孪生通过相邻(112)晶面上连续发射<111>/6不全位错形成（图6j,k），该机制对氢的敏感度远低于全位错发射。这一发现与实验报道的氢添加改善CoCrFeMnNi高熵合金力学性能的现象具有内在一致性——孪生主导变形的材料对氢脆的敏感性天然较低。

5  拉伸过程的外推精度验证

为确认大变形条件下的预测可靠性，研究基于局部原子环境计算了拉伸过程中所有约30万原子的外推等级（extrapolation grade）（图7）。对位错主导的GGB 1和孪生主导的GGB 4，无论氢浓度与应变条件，绝大部分原子始终处于内插（interpolation）区域（外推等级<1），仅有极少量原子接近精确外推边界（图7a–d）。这表明MLIP在拉伸过程中未遭遇超出训练数据范围的新原子环境，所有模拟结果均具有可靠的DFT精度。

图7  拉伸过程MLIP外推精度验证：(a-b) GGB 1和GGB 4全部原子的外推等级随应变的变化；(c-d) 12%应变时GGB 1和GGB 4的原子结构与外推等级着色图；(e-f) GGB 4单轴拉伸试验中原子结构中所有原子平均和最大外推级的应变依赖性

6  讨论与工程意义

6.1  氢抑制位错发射的原子机制

模拟结果直接支持了实验上关于氢偏聚于晶界可抑制位错发射的推测。氢占据晶界附近的膨胀八面体间隙位后，提高了位错线从晶界向晶内突出所需应力，这与HELP机制（氢降低位错运动阻力促进滑移）看似矛盾，实则反映了氢在不同应力/浓度阶段的差异化影响：低浓度促进位错运动，高浓度偏聚反而抑制位错发射，迫使晶界在低应变下发生脆性断裂。这一发现为理解HEDE与HELP在高强合金晶界氢脆中的协同竞争提供了新视角。

6.2  对高强合金抗氢脆设计的启示

本研究从原子尺度揭示了晶界本征变形机制是决定氢脆敏感性的核心参数：位错主导变形的晶界为高敏感性，变形孪生主导的晶界为低敏感性。这为抗氢脆微结构设计指明了方向：（1）通过晶界工程增加特殊晶界（低Σ值）比例；（2）利用合金化调控层错能，促进变形孪生而非位错滑移；（3）引入纳米析出相阻碍氢向一般晶界的扩散和偏聚。此外，团队正基于本研究训练数据构建精度-效率更优的新一代MLIP，未来结合加速分子动力学方法，有望在更接近工程应变率的条件下系统探索晶界特征与氢脆敏感性的构效关系。

7  总结与展望

本研究通过并行学习策略成功构建了Fe-H二元体系的BNNP型MLIP，实现了一般晶界氢偏聚、位错发射、变形孪生及裂纹萌生扩展的DFT精度描述。主要成果可归纳为三点：（1）证明了系统覆盖氢脆相关原子环境的训练集设计策略可使MLIP对训练集外构型（螺位错-氢、一般晶界-氢）保持优异的外推精度；（2）从原子尺度揭示了一般晶界的两种本征变形模式——全位错发射与变形孪生，以及氢对两者的差异化影响；（3）氢通过偏聚于晶界膨胀八面体位、提升位错发射能垒的方式抑制全位错发射并促进晶界断裂，而对变形孪生的影响极其有限。

本研究为高强铁基合金的抗氢脆晶界工程提供了原子尺度力学判据。未来工作将聚焦于：（1）研发精度-效率更优的MLIP以覆盖更广泛的应变率与晶界特征空间；（2）引入加速分子动力学方法实现更接近工程实际的低应变率模拟；（3）将MLIP构建方法学推广至Fe-C-H等多元素体系，为碳中性背景下高强结构合金的理性设计提供高通量计算工具。

需注意的是，当前分析受限于高应变率条件（MD固有局限）和有限数量的晶界模型，且晶界迁移等因素在某些条件下（如GGB 3高氢浓度）亦发挥额外作用。实验验证方面，Sugiyama等人的观测为本研究的位错发射抑制结论提供了间接支持，但变形孪生对氢不敏感的预言有待进一步实验检验。

原文阅读：Machine learning interatomic potential reveals hydrogen embrittlement origins at general grain boundaries in α-iron | Communications Materials