氢气泡在纳米孔洞内的形成是导致金属材料氢致脆化的关键机制，但其原子尺度起源长期未被揭示→本研究开发了首个能够准确捕捉氢陷阱和纳米孔洞中H₂形成过程的机器学习原子间势（NEP-WH），在保持DFT计算精度的同时实现了百万原子级模拟→大规模模拟揭示了独特的氢聚集路径：平面氢团簇在{100}晶面附近形核生长，高浓度区域形成HCP结构→单轴拉伸实验表明，这些氢聚集体抑制位错发射，促进{100}裂纹扩展，最终导致脆性断裂→这一发现为理解氢致损伤提供了原子级视角，为未来聚变堆钨部件的设计提供了关键指导。

背景定位

氢致退化是材料科学领域的百年难题，曾导致无数结构金属的过早失效。钨（W）因其在极端环境下的优异性能，被认为是未来聚变反应堆面向等离子体部件的首选材料。然而，实验观察到钨在氢富集环境中会出现大量氢气泡和失效现象。令人困惑的是，尽管实验证据充分，氢气泡形核与生长的原子级过程却始终是个"黑箱"——现有表征技术难以直接观测纳米尺度气泡的内部结构和压力状态。这个领域长期被两类方法所困扰：第一性原理计算精度高但只能处理几百个原子，经验势计算效率高却无法可靠描述氢-缺陷相互作用。

追问升级

以往研究存在一个关键瓶颈：没有一种原子间势能同时满足"高精度"和"大尺度"的需求。传统经验势（如EAM、BOP）在预测氢陷阱行为时与DFT结果偏差巨大，更无法复现实验观测到的GPa级气泡内压。2023年发表的DP-WH模型虽然有所改进，但其训练集未包含多空位-多氢原子的复合构型，预测氢气泡动力学的能力存疑。本文的核心追问是：能否开发一个既能达到DFT精度、又能模拟百万原子体系的机器学习势，从而直接"看见"氢气泡从形核到破裂的全过程？

方法优势

本研究采用了神经演化势（NEP）框架，这是一种将神经网络与进化算法相结合的机器学习方法。通过主动学习策略——即利用模型集成的不确定性估计自动探索构型空间——研究团队构建了包含57,400个结构、827万个原子的训练数据集。关键创新在于：训练数据涵盖了从完美晶体到缺陷结构、从熔融态到受应变态的广泛构型，特别纳入了纳米孔洞与氢的复合构型。最终的NEP-WH模型在单块GPU上就能实现450万原子的模拟，计算速度比DP-WH快两个数量级，仅比经验势慢几倍——这是该领域首次实现精度与效率的兼得。

核心发现

最令人震惊的发现是氢致脆化的全新机制：与传统认知中氢"促进位错运动"不同，本研究发现氢气泡周围的氢聚集体会*抑制*位错发射。具体而言，氢原子在{100}晶面形成平面团簇（planar hydrogen clusters），这些团簇像"预制裂纹"一样在拉伸载荷下直接引发脆性 cleavage fracture，裂纹沿氢团簇平面扩展，完全绕过位错机制。氢浓度越高，脆性越显著——从无氢时的韧性断裂（应变更0.12）转变为高氢浓度时的脆性断裂（应变仅0.06）。这解释了实验观测到的钨中{100}面优先开裂现象。

研究闭环

从问题揭示到机制阐明，本研究形成了完整链条：开发高精度势函数→模拟氢气泡形核生长→发现{100}平面氢团簇→揭示其抑制位错、促进裂纹的机制→与实验观测的{100} blisters直接关联。这项工作不仅解决了该领域长期存在的精度-效率矛盾，更建立了从原子尺度模拟到宏观失效行为的桥梁，为聚变堆材料的长寿命设计提供了理论依据。

研究亮点与数据支撑

核心发现1：NEP-WH模型达到DFT精度且计算效率极高

• 证据支撑
  ：在57,400个结构上训练后，能量预测RMSE仅4 meV/原子，力预测RMSE仅122 meV/Å，覆盖了纯W和W-H系统的晶格、缺陷、表面、晶界等性质
• 图表解读
  ：

- (a) 训练数据集的Sketch-map可视化：每个点代表一个原子构型投影到2D主成分空间，蓝色/白色大球为BCC/其他结构的W原子，橙色小球为H原子，左半边为纯W构型（含完美、缺陷、熔融、间隙原子等），右半边为W-H二元系统和H₂分子

- (b) 能量预测对比：NEP-WH预测值vs DFT参考值，RMSE=4 meV/原子

- (c) 力预测对比：NEP-WH预测值vs DFT参考值，RMSE=122 meV/Å

- (d) 计算速度对比：在单块NVIDIA V100 GPU上，NEP-WH达到4.5×10⁶原子·步/秒，比DP-WH快约100倍，仅比EAM-Wang慢几倍

• 结论总结
  ：NEP-WH模型首次实现了W-H系统模拟中DFT精度与百万原子级规模的兼得，为后续氢气泡动力学模拟奠定了基础。

核心发现2：NEP-WH是首个能准确描述纳米孔洞中氢陷阱和H₂形成的势函数

• 证据支撑
  ：在V₄Hn（4空位+ n个氢原子）体系中，NEP-WH准确复现了DFT计算的氢陷阱能，而DP-WH在H₂分子形成后出现显著偏差，其他经验势在所有阶段均表现糟糕
• 图表解读
  ：

- (a) NEP-WH预测 vs DFT：氢陷阱能随H原子数变化，黑色圆圈为DFT参考，NEP-WH（红色）几乎完美重合

- (b) DP-WH预测 vs DFT：在n>10后出现明显偏离，无法正确描述H₂形成

- (c) BOP-Juslin：整体偏差大

- (d) BOP-Li：整体偏差大

- (e) EAM-Wang：在H₂形成区域严重偏离

- (f) EAM-Mason：表现最差

- 注：空心符号表示H吸附在孔洞表面，实心符号表示H₂分子在孔洞内形成

• 结论总结
  ：这是目前唯一能可靠模拟氢气泡形核与演化的W-H势函数，解决了该领域长期以来的关键技术瓶颈。

核心发现3：NEP-WH在高压下准确描述H₂分子结构

• 证据支撑
  ：在300 K、5 GPa条件下，NEP-WH准确复现了DFT-MD计算的H-H径向分布函数（RDF）的两个峰——分别对应分子内H-H距离和分子间H-H距离
• 图表解读
  ：

- (a) BOP模型：两个峰均无法正确重现

- (b) EAM-Wang：第一个峰（分子内距离）偏离，第二个峰（分子间距离）缺失

- (c) EAM-Mason：同样无法重现分子间距离

- (d) DP-WH：第一个峰准确但第二个峰缺失

- (e) NEP-WH：两个峰均与DFT参考（黑色）完美吻合

• 结论总结
  ：NEP-WH是唯一能够可靠模拟气泡内GPa级H₂气压的势函数，为后续气泡压力计算提供了保障。

核心发现4：氢气泡压力与孔洞半径呈反比关系，验证Young-Laplace定理

• 证据支撑
  ：对半径0.75-5 nm的纳米孔洞进行氢溶解模拟，平衡气泡压力分别为18.73 GPa、17.58 GPa、16.95 GPa、16.01 GPa，与γ/R（γ为表面能）呈完美线性关系
• 图表解读
  ：

- (a) 气泡压力随模拟时间变化：对于0.75 nm、1 nm、1.25 nm、1.5 nm四种孔洞，H₂分子快速填充后压力达到平衡，插入图显示压力与γ/R的线性关系

- (b) 1 nm孔洞在0.109 ns的原子构型：H₂分子占据孔洞核心，H原子吸附在表面

- (c) 1 nm孔洞在1.417 ns的构型：显示{100}面附近的平面氢团簇开始形成

- (d) 1 nm孔洞在5.450 ns的构型：平面氢团簇扩展，交叉处形成HCP结构（橙色）

- 注：(c)(d)中粉色连线表示氢原子簇，H原子按[100]方向坐标着色（蓝到红），W原子颜色：CNA方法识别为BCC灰色、FCC绿色、HCP橙色、其他白色

• 结论总结
  ：首次在原子模拟中验证了宏观的Young-Laplace方程，为氢气泡热力学提供了定量数据支撑。

核心发现5：氢浓度驱动韧脆转变，强度随氢含量下降

• 证据支撑
  ：单轴拉伸模拟显示，无氢时极限强度27.46 GPa、应变~0.12（含位错发射的韧性断裂）；氢浓度1.45 at.%时强度降至24.73 GPa；2.62 at.%时强度21.44 GPa；4.96 at.%时强度仅17.73 GPa、应变~0.06（脆性断裂）
• 图表解读
  ：

- 应力-应变曲线：紫色实线（0 at.%氢）最高，蓝色虚线（1.45 at.%）、橙色点划线（2.62 at.%）、粉色点线（4.96 at.%）依次降低，高氢浓度下极限应变显著减小

• 结论总结
  ：氢气泡导致的材料软化效应触发了典型的韧脆转变，直接揭示了氢致失效的力学本质。

核心发现6：氢抑制位错发射，促进{100}裂纹扩展

• 证据支撑
  ：无氢时位错从孔洞表面发射并产生塑性变形；有氢时平面氢团簇引发应力集中，直接导致沿{100}面的裂纹形核与扩展，位错活动被完全抑制
• 图表解读
  ：

- (a)(b) 0 at.%氢：孔洞在应变下发射1/2⟨111⟩位错环（橙色线），裂纹穿过孔洞中心，断裂面不规则

- (c)(d) 1.45 at.%氢：拉伸过程中氢原子迁移至应力集中区形成(100)平面团簇，裂纹对称扩展

- (e)(f) 2.62 at.%氢：预存在的{100}平面团簇引导裂纹沿团簇面扩展，产生不对称裂纹

- (g)(h) 4.96 at.%氢：HCP相（橙色）与平面团簇共存，位错发射后裂纹继续沿新形成的团簇扩展，形成高度不对称的裂纹

- 视图方向：[010]方向，[001]×[100]面，5 Å厚的切片

- W原子颜色：CNA识别结果；粉色球为H原子；橙色线为1/2⟨111⟩位错，红线为⟨100⟩位错

• 结论总结
  ：氢团簇"劫持"了断裂路径——从位错主导的韧性模式转变为氢团簇引导的脆性cleavage模式，这是氢致脆化的原子级机制。

核心发现7：高氢浓度下应力分布与裂纹扩展特征

• 证据支撑
  ：von Mises应力分析显示，平面氢团簇和HCP区域产生显著的局部应力集中，加速裂纹萌生与扩展
• 图表解读
  ：

- (a)(b) 2.62 at.%氢：裂纹萌生前后的应力分布，{100}平面团簇区域应力最高

- (c)(d) 4.96 at.%氢：同时含平面团簇和HCP区域，HCP相耗尽后发射⟨100⟩位错（红线），应力重新分布

- 视图：[010]方向，5 Å切片，W原子按应力着色，H原子粉色

• 结论总结
  ：氢聚集体的空间分布直接决定了裂纹扩展的路径和形态，解释了实验观测到的{100}面优先开裂现象。