钨PFM氢氦辐照理论模拟

你大概想不到，人类离聚变能最近的那个零件——托卡马克偏滤器面向等离子体的那层钨——它的理论摸拟正处在一个奇特的十字路口。一边是机器学习势被捧上天：速度提升百万倍，精度逼近第一性原理。一边是五份从不同角度切入的报告交叉审视后发现：模拟内部算得热闹，但与实验现实脱节得令人不安。

这个路口离 ITER 的偏滤器不到十年。

钨（W）被选为面向等离子体材料（PFM）不是因为别的——熔点 3422°C，高-Z 元素溅射率低，热导率在室温下约 170 W/mK。聚变堆偏滤器表面要承受 10 MW/m² 的热流和 14.1 MeV 中子的轰击，以及同时涌入的氢（H）、氘（D）、氚（T）和氦（He）等离子体。H 会滞留、会起泡。He 会在近表面形成纳米气泡，撑出蓬松的绒毛状结构——"fuzz"。两者一起作用时，更会产生协同效应，让损伤叠加。

理论摸拟的任务是从原子尺度理解这些过程，然后预测零件能用多久。这个任务有多难？DFT 计算几百个原子，时间尺度是飞秒。实验关心的场景是几个月到几年。中间差十二个数量级。

*He 会跑。跑的还很快。*

这一点是五份独立报告最大的共识之一。He 在完美钨晶格中从一个四面体间隙跳到相邻四面体的势垒只有 0.06 eV——DFT 给出 0.06 eV，分子动力学给出 63 meV，本质上是同一个数。这意味着在偏滤器的工作温度（600-1200 K）下，He 原子每皮秒跳跃几万亿次，几乎不受晶格阻碍。它会在晶格中飞快地游荡，直到撞上一个缺陷（空位、晶界、位错）停下来——然后系统性地聚集。

问题是，它停下来以后发生了什么。

*H 会停下。但停下来以后又被别的物种改变了。*

H（以及它的同位素 D 和 T）的行为更复杂。H 在完美钨中的迁移势垒约 0.30 eV，也是快的。但它一旦遇到空位，结合能高达 1.23 eV——空位就像一口深井，一个空位最多可以装下 12 个 H 原子。在偏滤器环境下，H 的滞留总量可以达到纯体扩散预测的 28 倍（FESTIM 模型，一种有限元连续模拟工具的预测）。

H 和 He 在一起时，事情变得更奇怪。He 先进去占据空位，形成 He-V 团簇。H 进不去核心，就在外壳聚集——形成一种"核-壳"结构，He 在内、H 在外。ASIPP（中科院等离子体物理所）李晓春团队的经验势模拟和 UKAEA 的实验观察都证实了这一点。

但故事还没完。

当级联碰撞（14.1 MeV 中子在晶格中撞出的连锁反应）发生时，预存的 H 或 He 把损伤放大了 5-15 倍（Abu-Shams 等人 2021 年的发现）。Jiang 等人在 2025 年《Journal of Nuclear Materials》上发表了目前最全面的级联模拟：PKA（初级撞出原子）能量从以往的 50 keV 扩展到 100 keV。他们发现 He 显著增加 Frenkel 缺陷对的产生，机理有两个：一是 He 占了空位，减少间隙-空位的复合；二是 He 改变了阈位移能，让钨原子更容易被撞走。而 H 的影响相比之下微妙得多——它使空位团簇尺寸减小，不影响位错环类型。

*然就在这时，我们发现了一个"抵消"效应。*

Hatton 等人（2025, LANL/UT Knoxville）在 MD 模拟中看到：纯 He 气泡会产生 glissile（可动的）1/2⟨111⟩ 位错环，这些环滑移到表面沉积钨吸附原子，是 fuzz 生长的前驱。但当 He:H 比例达到约 1:1 时，位错环类型发生了切换——从可动环变成了 sessile（不可动）的 ⟨100⟩ 环。表面沉积的钨吸附原子从纯 He 时的 161 个暴降至几乎为零。也就是说，H 在一定程度上*抑制*了 He 造成的表面形貌改变。

这个"抵消效应"的微观机制有两层：中低 H 浓度下，H 通过 W-H 结合抑制已有自间隙原子形成位错环；高 H 浓度下，H 增加 Frenkel 对形成的能量势垒，从源头抑制。

报告之间另一个张力点埋得更深。四份独立报告从不同路径确认了"协同增强滞留"——大连理工倪维元/刘东平组的速率理论模型（Nucl. Fusion 2024）、FESTIM 连续模拟 28 倍增强、USTC 张媛媛的团簇动力学 IRadMat-TDS 模型（Nucl. Fusion 2024）——都指向同一个方向：H 和 He 在一起时，滞留比单独存在时多得多。大连理工的模型还被 VTT Finland（芬兰国家技术研究中心）在其欧洲聚变项目中使用。

然而报告 E 列出的一串系统偏差让人醒过来。Anderl 复合系数——多年来被广泛使用的氢表面复合系数——被发现是错的。用它做的模拟与实验严重不符。W-H 经典 Tersoff 势根本无法再现氢分子的复合与脱附，因为三体作用不足以描述这个过程。纳米气泡中的平衡压力——Young-Laplace 方程预测 4.40 GPa，MD 得出 2.08 GPa，偏差超过两倍。He 气泡的状态方程（最常用的是 Mills-Liebenberg-Bronson 自由气体模型）从 933 K 起系统性地高估压力，温度越高偏差越大。

模拟内部自洽，但与现实脱节。

*这里插进来一个更根本的争议：机器学习势到底可不可信？*

报告 A 像推介会一样列出了 MLIP 的发迹史：Deep Potential, NEP, SNAP, tabGAP, NNAP, MTP——五种以上的 ML 势框架在 2025-2026 年爆发。北科大包等人的 NEP-WH 势能量 RMSE 4 meV/atom，计算速度达到 4.5×10^6 步/秒（在 V100 GPU 上），是现有 W-H 二元机器学习势中最精确的之一。它首次捕捉到了纳米空洞中的 H₂ 分子形成过程，揭示了 H 沿 {100} 面形成 HCP 结构进而驱动裂纹扩展的图景。Sandia 的 SNAP 势已经可以在 Perlmutter 超算上部署。

但报告 E 泼了一盆冷水：MLP 的三重争议——精度与效率不可兼得、可迁移性不可靠、W-He 二元开发尤其困难——被一个接一个摆出来。MLIP 本质上是插值工具。当模拟条件超出训练集覆盖时（高温、高缺陷密度、复杂应力场），漂亮数字会毫无预警地崩溃。不同势函数在同一位错环相对稳定性上的预测甚至可以定性相反——EAM 倾向产生 ⟨111⟩ 和 ⟨100⟩ 位错，而更精确的 MLP 预测占主导的可能是 C15 环状结构。

更令人惊讶的是，在这个 MLIP 大辩论中，有一个重要的方法论进展几乎没有其他报告触及。日本 QST 的 Saito 等人 2026 年开发了三阶段 3D-CNN 流水线来预测迁移势垒，平均绝对误差 0.124 eV，每势垒推断只要 2.7 ms，比传统 NEB 方法快约 23000 倍——实现了 O(1) 复杂度。但在当前的五份报告全景中，这项突破就像一颗在无人区炸开的照明弹。

*还有一件事，叫核量子效应。*

合肥工大孙飞团队 2025 年首次用路径积分分子动力学（PIMD）研究了钨空位中的 H 同位素交换。他们发现在 500 K 以下，核量子效应显著改变 H/D/T 的行为——也就是经典模拟中完全被忽略的原子核零点振动效应，对同位素差别有实质性影响。ITER 每年要处理公斤级的氚滞留，这个发现在五份报告中只出现在报告 A 里，完全被其他报告遗漏。

中国团队在这张图景中的位置有点复杂。大连理工的速率理论模型在国际上局部领先——倪维元和刘东平 2024 年 Nucl. Fusion 文章是国际上首次定量揭示 H-He 协同增强高温滞留机制的速率理论模型。USTC 潘必才/王炯荣的紧束缚法可以处理电子激发态效应，这是国内唯一、国际上与法国 CEA 并行的方向。北科大包 2025 年做出了国内首篇 W-H 机器学习势论文。ASIPP 李晓春的三元键序势被 UKAEA 作为欧洲聚变材料模拟的参考势之一。

但人才密度是个硬缺口。报告 D 识别出中国核心团队约 5-7 个，总人数 30-40 人——对比美国 DOE PSI-SciDAC 和欧洲 EUROfusion 的百人级规模，差距至少三倍。ML 势覆盖率只有 W-H 二元（国际已到 W-H-He 三元）。高通量计算近乎空白。嬗变产物 Re/Os 与 H/He、温度梯度、应力梯度等多物理场全耦合模拟——中国发表了零篇。

综合五份报告，最有价值的判断可能不是某一个具体的物理发现。fuzz 形成的速率限制步骤是什么？是钨表面迁移还是 He 扩散？不知道。晶界对 He 到底是"通道"还是"陷阱"？不知道。V3（三空位）介导的空位聚集机制比传统单空位扩散快一千倍——自适应 KMC 算出来了——但在其他模拟框架里没有体现。每个团队都在优化自己那块拼图的清晰度，整幅画布却没人画。

2026 年 7 月，我们同时拥有三项互相矛盾的事实：MLIP 在飞速进步，但经典势在不同体系上的 >50% 差异从未被系统解决；H+He 协同滞留被多路径独立确认，但几乎所有理论预测都建立在未被实验充分校准的模型上；中国有局部领先（速率理论、紧束缚、三元势），但人才密度和计算基础设施差距巨大。

偏滤器那层钨还是会被装上 ITER。问题是，我们算得够准了吗？

---

参考文献

[1] Ni W, Niu C, Liu D et al. Temperature-dependent bubble growth under synergistic interactions of hydrogen and helium in tungsten. *Nuclear Fusion*, 2024, Vol. 64, 056039.

[2] Jiang et al. Impacts of helium and hydrogen on the defect evolution in tungsten under high-energy cascades: A molecular dynamics study. *Journal of Nuclear Materials*, 2025, Vol. 611.

[3] Hatton P, Perez D, Uberuaga B. Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten. *Journal of Nuclear Materials*, 2025, Vol. 605, 155555.

[4] Bao et al. NEP-WH: a machine learning potential for tungsten-hydrogen systems. *npj Computational Materials*, 2026, Vol. 12, 108. arXiv:2508.20350.

[5] Li X C et al. Analytical W-H-He bond-order potential for plasma-facing materials. *Journal of Nuclear Materials*, 2025.

[6] Tirumala, Mason et al. An empirical potential to simulate helium and hydrogen in highly irradiated tungsten. *J. Phys.: Condens. Matter*, 2025, Vol. 38. arXiv:2509.22012.

[7] Kaur, Mason et al. A general multi-occupancy trapping site model for hydrogen in tungsten. *Physical Review Materials*, 2025, Vol. 9, 125404.

[8] Kong et al. Temperature-dependent trap mutation criterion for He bubbles in tungsten. *Materials*, 2025, Vol. 18(15), 3518.

[9] You et al. Transmutation Re/Os reduce hydrogen retention in tungsten. *Nuclear Fusion*, 2026, Vol. 66, 056001.

[10] Sun F et al. Nuclear quantum effects in hydrogen isotope exchange in tungsten vacancies: A PIMD study. *Nuclear Fusion*, 2025.

[11] Zhang Y et al. IRadMat-TDS model for hydrogen retention in tungsten. *Nuclear Fusion*, 2024.

[12] Xu Y C et al. Insight into the loop-punching mechanism for He bubble growth in tungsten. *Computational Materials Science*, 2023.

[13] Xu Y C et al. Synergistic influence of dislocations and helium clusters on hydrogen atoms in tungsten. *Journal of Nuclear Materials*, 2024, Vol. 596, 155110.

[14] Dasgupta D, Maroudas D, Wirth B D. Surface morphological response of plasma-facing tungsten: Effects of surface crystallographic orientation. *Acta Materialia*, 2024.

[15] Lasa A et al. Xolotl: A multi-scale framework for plasma-facing materials. *Journal of Nuclear Materials*, 2024, Vol. 594, 155011.

[16] Saito et al. 3D-CNN accelerated transition state search for hydrogen in tungsten. QST Japan, 2026.

[17] Cui et al. Molecular dynamics simulations of cascade overlap with void/helium bubble in tungsten. *Nuclear Materials and Energy*, 2022.

[18] Wang et al. Possible approaches for simulating the formation of fuzz structure on tungsten surface under helium irradiation. *Computational Materials Science*, 2024.

[19] Byggmästar et al. Machine learning interatomic potential for tungsten: training and validation. *Acta Materialia*, 2025/2026.

[20] Fedrigucci et al. High-throughput DFT screening of plasma-facing material candidates. *PRX Energy*, 2024, Vol. 3, 043002.

[21] Abu-Shams et al. Displacement cascade evolution in tungsten with pre-existing helium and hydrogen clusters. *International Journal of Materials Research*, 2021.

[22] Li Y G, Zeng Z et al. Multiscale simulation platform IMSM for PFM. theory.issp.ac.cn/IM3D.

[23] Wang J R, Pan B C et al. Non-adiabatic tight-binding study of electron excitation effects on H-defect interactions in W. 2023-2026.