1. 项目概述从原子碰撞到材料失效我们如何用AI看清辐照损伤在核反应堆的堆芯、航天器的外壁或是未来聚变装置的面向等离子体部件材料无时无刻不在承受着高能粒子的“轰击”。这些粒子比如中子或离子携带着巨大的能量撞入材料内部就像一颗高速子弹射入一片紧密排列的台球阵中。最初的撞击初级撞出原子PKA会引发一连串的原子级“车祸现场”——级联碰撞。这个过程在皮秒万亿分之一秒量级内产生大量脱离晶格位置的原子间隙原子和留下的空位即弗伦克尔对。这些点缺陷并非静止它们会迁移、聚集形成位错环、空洞等缺陷团簇最终导致材料肿胀、硬化、脆化性能急剧下降。理解并预测这一微观损伤过程是设计下一代抗辐照材料的关键。然而实验手段难以实时捕捉原子尺度的瞬态过程而理论计算的第一性原理方法又受限于计算量无法模拟包含数百万原子、持续上百皮秒的级联碰撞事件。这时分子动力学模拟成为了连接微观机制与宏观性能的桥梁。但这座桥的稳固性完全取决于其基石——描述原子之间如何相互作用的“势函数”。传统的经验势函数如嵌入原子法势其函数形式固定在描述复杂多组分合金如包含钼、铌、钽、钒、钨五种元素的高熵合金时往往力不从心难以准确预测缺陷形成能、表面能等关键性质导致模拟结果与真实物理图像存在偏差。近年来机器学习势函数的出现为这一困境带来了转机。它不再依赖预设的物理公式而是像一个极度聪明的“学徒”通过“学习”大量第一性原理计算的高精度数据自己总结出原子间相互作用的复杂规律。这种数据驱动的方法使得我们能够以前所未有的精度和效率模拟复杂材料中的辐照损伤过程。与此同时如何量化损伤也至关重要。长期以来NRT-dpa模型是国际标准但它简单地将所有原子视为等同忽略了不同元素在质量、尺寸上的差异对于高熵合金这类“化学复杂”的材料其预测偏差可能很大。因此我们面临的核心问题是能否结合高精度的机器学习势函数与更物理的损伤量化模型深入揭示像MoNbTaVW这样的难熔高熵合金其卓越抗辐照性能背后的原子机制这项工作正是对此的探索。我们将训练一个针对Mo-Nb-Ta-V-W体系的机器学习势函数进行从1 keV到150 keV能区的系统级联模拟并发展一个考虑元素特性的扩展arc-dpa模型来量化损伤。最终目标是看清化学复杂性如何通过“亚级联分裂”等机制巧妙地抑制有害缺陷团簇的形成从而赋予材料更强的抗损伤能力。无论你是材料模拟领域的研究者还是对先进核材料设计感兴趣的工程师本文将带你深入这一交叉前沿理解从算法到物理图像的全链条逻辑。2. 方法论基石构建可靠的模拟与量化框架要开展可靠的辐照损伤模拟研究必须搭建一个坚实的计算框架。这包括两个核心部分一是能够准确描述MoNbTaVW高熵合金中复杂原子间相互作用的势函数二是能够合理量化这种多组分材料中辐照损伤的物理模型。两者缺一不可前者决定了我们看到的“电影”是否逼真后者决定了我们如何科学地“统计”电影中的关键情节。2.1 为复杂合金“量身定制”机器学习势函数传统经验势函数在模拟高熵合金时的主要瓶颈在于“泛化能力”不足。一个针对纯钨优化的EAM势很难准确预测钒原子掺入后形成的间隙原子构型及其能量。机器学习势函数则采取了“海量数据灵活模型”的策略。2.1.1 训练数据的构建覆盖所有可能场景我们采用的训练数据集并非随意生成而是经过精心设计旨在覆盖材料可能经历的所有物理状态和缺陷类型确保势函数具有广泛的适用性。数据集主要包含以下几类结构纯金属基础数据包括每种元素Mo, Nb, Ta, V, W的体心立方晶体、高温bcc相、包含空位和自间隙原子的缺陷结构以及表面和液态结构。这确保了势函数能准确描述每种元素的单质特性。合金成分空间采样从二元到五元的所有等原子比随机合金结构。这是理解高熵合金“鸡尾酒效应”的关键迫使模型学习不同元素随机混合时的相互作用。有序合金结构包含从二元到五元的有序合金。这有助于模型捕捉可能的短程有序或局部偏聚现象。高熵合金中的缺陷专注于MoNbTaVW成分构建包含最多五个空位或自间隙原子的缺陷超胞。特别地为了精确描述级联碰撞中原子极度接近时的强排斥作用数据集还包含了在HEA晶体中随机插入一个间隙原子、且使其与近邻原子距离非常近但又不至于过近导致不物理的结构以捕获短程排斥力。液态结构涵盖了从二元到五元的所有等原子比组合的液态。这对于模拟级联过程中产生的局部熔化区热峰至关重要。整个数据集包含21,672个原子构型和334,858个原子并通过第一性原理计算获得了每个构型的能量、原子受力和维里应力。如此全面的数据确保了训练出的势函数能够可靠地模拟材料的弹性、热学、缺陷性质以及表面能、熔化过程和液相结构。2.1.2 模型选择与训练神经进化势与ZBL修正我们选用神经进化势作为机器学习势函数的框架。NEP模型以其高精度和优异的计算效率著称其“NEP4”版本在描述多组分体系方面表现尤其出色。描述符与神经网络模型的输入是每个原子周围的局部化学环境通过一组径向和角向描述符来数学表征。我们设置的截断半径为5埃这意味着模型只“看到”每个原子周围5埃内的邻居。径向和角向部分各使用7个径向函数每个函数由9个切比雪夫多项式基函数线性组合而成其系数是可训练的。这些描述符信息随后输入到一个含80个神经元的单隐藏层神经网络最终输出该原子的能量贡献。损失函数训练目标是让NEP预测的能量、力和维里应力尽可能接近第一性原理计算值。损失函数是这三项误差的加权和权重分别设为1, 1, 0.1。这意味着我们最看重力和能量的准确性同时也适当考虑应力。短程排斥势的融合在级联模拟的初始阶段PKA能量极高原子间距离可能小于1埃此时电子云重叠相互作用以强库仑排斥为主。NEP主要从平衡态附近的数据学习对这类极端短程相互作用的描述可能不够精确。为此我们采用了NEP-ZBL方案将经典的Ziegler-Biersack-Littmark屏蔽库仑势与NEP势平滑耦合。总势能表示为U_i U_i^NEP 1/2 * Σ_{j≠i} U^ZBL(r_ij)。在1-2埃的过渡区间NEP被训练去拟合DFT与ZBL势的差值从而保证从短程强排斥到中长程相互作用的平滑过渡。最终训练出的模型其能量、力和维里应力的预测均方根误差分别为4.37 meV/原子、145 meV/埃和47.8 meV/原子与DFT结果高度吻合。更重要的是在单块RTX 4090 GPU上该模型达到了每秒1千万原子步的计算速度与传统的EAM势相当使得模拟包含数千万原子的大型级联过程成为可能。2.2 从“平均原子”到“个性原子”扩展损伤模型量化辐照损伤的经典模型是NRT-dpa其公式N_d 0.8 * T_d / (2 * E_d)隐含了一个关键假设所有原子都是相同的。它用一个平均的离位阈值能量E_d来代表整个材料。这对于纯金属或简单合金尚可接受但对于MoNbTaVW这种五种元素随机分布的高熵合金这个假设就过于粗糙了。钒原子小而轻钨原子大而重它们被撞离位的难易程度E_d截然不同。因此我们提出了针对多组分材料的扩展损伤模型。其核心思想是在一次碰撞事件中考虑碰撞双方原子的类型差异并引入一个有效离位阈值能量E_multi。2.2.1 模型推导的核心逻辑推导从基本的Kinchin-Pease模型开始。考虑一次碰撞入射原子类型i将能量T传递给靶原子类型j。产生的离位原子数ω_ij(E)取决于入射原子剩余能量(E-T)产生的离位加上靶原子获得能量T后可能产生的离位。通过对所有可能的能量转移T进行积分平均并考虑不同类型原子碰撞的概率我们得到了适用于多组分体系的总离位数表达式。其中最关键的是靶原子选择概率p_i的引入p_i (c_i * r_i^2) / Σ_j (c_j * r_j^2)这里c_i是原子浓度r_i是原子半径。p_i正比于该种原子的丰度浓度和它的“靶子面积”正比于半径平方。这意味着在随机碰撞中一个原子被选为靶原子的概率既取决于它有多少浓度也取决于它“个头”有多大截面。基于此我们定义了有效离位阈值E_multiE_multi Σ_i p_i ( E_d^i Σ_j p_j E_d^j / Λ_ij )其中E_d^i是第i种原子在多组分材料中的离位阈值Λ_ij 4M_i M_j / (M_i M_j)^2是质量因子表征从i到j的最大能量传递效率。这个公式的物理意义很直观E_multi是各种原子作为靶原子时其离位难易程度E_d^i以及它被其他类型原子撞击时能量接收效率通过Λ_ij和p_j加权的综合体现。2.2.2 扩展的arc-dpa模型将E_multi代入到经典的arc-dpa模型中我们就得到了扩展版本N_d, arc-dpa(T_d) 0.8 * (0.8 * T_d / E_multi) * ξ_arcdpa(T_d), 当 T_d E_multi/0.8其中ξ_arcdpa(T_d) 1 - c_arcdpa * (E_multi/0.8)^b_arcdpa * T_d^b_arcdpa c_arcdpa。b_arcdpa和c_arcdpa是材料相关的参数需要通过分子动力学模拟或实验来确定。ξ函数描述了在高能区由于热峰效应导致的缺陷复合使得稳定存活缺陷数偏离线性增长的趋势。这个扩展模型通过E_multi和拟合的参数b,c将化学复杂性元素种类、浓度、质量、尺寸对损伤产生效率的影响纳入了量化体系。注意这个扩展模型基于原子随机均匀分布的假设。在实际高熵合金中可能存在短程有序或局部化学起伏这可能会对碰撞概率产生影响。本工作作为初步探索采用了随机固溶体假设更复杂的化学有序效应值得在未来研究中深入探讨。3. 势函数验证与基本物性确保模拟的可靠性在将机器学习势函数投入大规模的级联模拟之前我们必须像检验一把新尺子的刻度一样对其预测材料基本物性的能力进行全面的“校准”。只有当它在各种测试中都与已知的实验数据或高精度第一性原理计算结果吻合良好我们才能信任它模拟复杂非平衡过程如辐照级联的结果。3.1 静态物性从晶格常数到熔点我们首先计算了纯金属及MoNbTaVW高熵合金的一系列基本性质并与文献中的实验值或DFT计算结果进行对比。3.1.1 晶格常数与弹性性质对于五种纯金属V, Nb, Mo, Ta, WNEP模型预测的晶格常数和体弹模量与DFT数据高度一致。例如钨的晶格常数预测值为3.186埃与DFT值3.185埃几乎相同其体弹模量预测为309 GPa与DFT的303 GPa也非常接近。对于MoNbTaVW高熵合金我们通过对50个包含2000原子的弛豫超胞取平均得到其平衡晶格常数为3.195埃体弹模量为219 GPa。这些值与之前的DFT研究结果吻合表明NEP模型准确地捕捉了合金的基态结构和弹性响应。一个有趣的趋势是钨和钼是这些元素中最“硬”的体弹模量最高当它们与较“软”的钒、铌、钽形成合金后整体的弹性刚度会有所下降。3.1.2 形成能与混合能单空位形成能是衡量材料中点缺陷稳定性的关键参数。在高熵合金中由于化学环境的复杂性每个空位周围的元素排列都不同因此其形成能有一个分布。我们计算了HEA中50个不同位置的单空位形成能。DFT计算的平均值为3.33 eV而NEP模型的预测平均值为3.29 eV分布范围也与DFT结果高度重合。这表明模型能准确描述空位在复杂化学环境中的形成能涨落。对于自间隙原子我们在1000个不同的HEA超胞中随机插入一个原子弛豫后观察其稳定的构型。结果显示含钒的哑铃状间隙原子构型尤其是V-V哑铃占据主导地位。这与之前的DFT和MD研究结论一致归因于钒原子尺寸最小倾向于形成更短的键从而在能量上更有利。3.1.3 熔点预测我们采用固液共存法估算了HEA的熔点。在零压条件下构建一个包含13500个原子的双相体系一半固体bcc相一半液相在熔点附近进行模拟。NEP模型预测的MoNbTaVW的熔点约为2840 K。需要指出的是对于合金固相线和液相线之间存在一个温度区间我们的模拟只给出了固液共存的一个温度点可以近似看作合金的熔化温度。这个值处于各组分纯金属熔点之间是合理的。3.2 动态门槛离位阈值能量离位阈值能量是辐射损伤研究中最基础的参数它定义了一个原子被永久撞离其晶格位置所需的最小动能。我们利用NEP模型在300 K温度下对纯金属和HEA进行了大规模的TDE计算。3.2.1 纯金属的TDE各向异性对于纯金属我们沿500个随机方向计算了TDE并绘制了极射投影图。结果显示钨的TDE最高平均约109 eV且其值分布范围广对晶向非常敏感在密排方向如111需要更高的能量才能离位。钼次之平均约99 eV。钒的TDE最低平均约54 eV最容易离位。 这个趋势与纯金属中间隙原子形成能的趋势一致形成能越高原子结合越牢越难被撞离位。3.2.2 高熵合金的TDE“均匀化”在MoNbTaVW高熵合金中我们为每种元素作为PKA计算了1000个随机方向的TDE。结果呈现出与纯金属截然不同的特征分布更集中五种元素PKA的TDE值分布都更接近高斯分布角度依赖性显著减弱。这是因为高熵合金中严重的晶格畸变和化学无序使得每个原子所处的局部环境都各不相同平均化了沿不同晶向的离位难度差异。TDE与原子质量的逆相关计算得到的平均TDE值为V (58 eV), Nb (53 eV), Mo (54 eV), Ta (45 eV), W (45 eV)。一个明显的趋势是较轻的PKAV, Nb反而需要更高的能量才能造成稳定离位。这与直觉似乎相悖因为轻原子更容易被加速。其深层因在于在HEA中一个轻原子如V被撞离位后它更容易与周围较重的原子发生非弹性碰撞而损失能量或者其形成的间隙原子构型如V-V哑铃迁移能力较强导致其在冷却阶段更容易与空位复合从而无法形成“稳定”的弗克尔对。统计也发现在5000次随机反冲事件中最终形成稳定缺陷的PKA85%是钒而较重的钽和钨分别只占0.6%和2.4%。这再次印证了V原子在缺陷产生过程中的活跃角色。基于这些元素特定的TDE值我们利用扩展模型公式计算了MoNbTaVW HEA的有效离位阈值能量E_multi 108 eV。这个值高于大多数纯组元的TDE表明从整体上看HEA中的原子“更难”被离位这初步提示了其潜在的抗辐照特性。4. 级联模拟全景缺陷产生、演化与量化有了经过验证的可靠势函数我们便可以深入辐照损伤的核心舞台——原子碰撞级联。我们系统模拟了PKA能量从1 keV到150 keV的碰撞过程对比了纯钨和MoNbTaVW高熵合金在缺陷产生、存活和团簇化行为上的差异。4.1 缺陷产量与扩展arc-dpa模型拟合我们统计了级联模拟最终阶段约100 ps后存活的稳定弗伦克尔对数量。图5展示了MoNbTaVW HEA和纯W在不同损伤能量T_d下的FP数量。4.1.1 与经典模型的对比对于纯钨我们采用了文献中报道的arc-dpa参数b -0.56,c 0.12。在低能区T_d 20 keV我们的MD模拟结果与arc-dpa模型预测吻合良好。但在高能区MD模拟产生的稳定FP数量显著高于模型预测。这种偏差是预期的因为arc-dpa模型基于纯金属数据建立且在高能区级联分裂、热峰效应等复杂过程使得缺陷产生效率偏离简单标度律。4.1.2 HEA的损伤特征与模型参数在整个能量范围内HEA中存活的FP数量始终高于纯钨。这意味着在单次碰撞事件中HEA产生了更多的点缺陷。这似乎与“抗辐照”的直观印象矛盾。然而抗辐照性能不仅取决于缺陷的“产量”更取决于缺陷的“命运”——它们是保持孤立、易于迁移复合还是聚集形成稳定的、有害的大团簇为了量化HEA的损伤我们用MD模拟数据拟合了扩展arc-dpa模型的参数。得到b -0.73,c 0.21。作为对比如果使用原始的、不考虑元素差异的模型即用平均TDE代替E_multi进行拟合会得到c 0.19。c值代表了在极高能量下缺陷存活分数趋近的常数。扩展模型给出了更高的c值表明它预测在高能区HEA会产生更多的损伤这与我们的MD结果更一致。这凸显了考虑化学复杂性对于准确量化多组分材料损伤的重要性。4.2 缺陷团簇化行为抑制与促进的博弈抗辐照性能的关键在于抑制缺陷特别是间隙原子团簇的聚集和长大因为它们是位错环的前驱体会导致材料硬化。我们分析了空位和间隙原子团簇的比例及尺寸。4.2.1 团簇分数与尺寸图6揭示了HEA与纯W在缺陷团簇化行为上的根本差异间隙原子团簇在纯钨中随着PKA能量增加间隙原子团簇的比例急剧上升在150 keV时接近90%。这意味着绝大多数间隙原子都聚集在了一起。相反在HEA中间隙原子的团簇比例始终低于25%且不随能量显著增加。空位团簇HEA中空位的团簇比例始终高于纯钨。同时HEA中最大空位团簇的尺寸也大于纯钨而最大间隙原子团簇的尺寸则远小于纯钨。4.2.2 现象解读亚级联分裂的主导作用这些数据指向一个核心机制在高熵合金中辐照级联倾向于“分裂”成多个更小的、空间分离的亚级联。这种“亚级联分裂”效应是HEA抗辐照性能提升的关键。对间隙原子亚级联分裂将产生的间隙原子分散到多个孤立的小区域极大地减少了它们相遇、聚集形成大团簇的机会。即使在小区域内形成了团簇其尺寸也有限。对空位空位迁移能力通常远低于间隙原子。亚级联分裂产生的多个小热峰区域在冷却过程中其中心区域更容易因原子快速填充而形成空位富集区从而促进了空位团簇的形成。但这些空位团簇相对孤立且后续与分散的间隙原子复合概率较低。这种“促进空位团簇、抑制间隙原子团簇”的现象与在WTaCrV等其它难熔HEA中的模拟观察一致。它改变了缺陷演化的格局有害的、可移动的间隙原子团簇被有效抑制而相对静止的空位团簇被保留。后者虽然也是缺陷但其对材料力学性能如硬化的直接影响通常小于位错环。4.3 PKA元素类型的影响谁触发了分裂为了探究化学复杂性如何影响级联形态我们在150 keV能量下针对HEA中五种不同的PKA元素Mo, Nb, Ta, V, W各进行了10次模拟共计50次以获取统计规律。4.3.1 亚级联分裂的触发概率统计分析给出了清晰的结论V和Nb作为PKA时10次模拟全部观察到了显著的亚级联分裂。Mo作为PKA时9次分裂1次未分裂。W和Ta作为PKA时约有一半的模拟发生了分裂。而在相同能量下纯钨的自我辐照模拟中完全没有观察到亚级联分裂。文献指出纯钨中自离子辐照产生亚级联分裂的阈值能量在160 keV附近。我们的结果表明HEA的化学无序显著降低了级联分裂的阈值能量。4.3.2 分裂与否的后果图8和图9清晰地展示了分裂与未分裂级联的迥异结局发生亚级联分裂的案例级联早期就分裂成多个分支峰值损伤时形成多个分散的损伤区。最终大多数点缺陷呈孤立或小团簇状分散间隙原子团簇化比例很低。未发生分裂的案例级联保持一个紧凑的、高能量密度的“热峰”。这会形成一个大的熔融区在冷却凝固后容易产生大型的间隙原子和空位团簇尺寸可超过100个缺陷。虽然这种大热峰在冷却过程中也会促进点缺陷复合导致最终FP总数可能比分裂案例少但一旦形成大团簇它们就非常稳定难以通过后续热迁移消除对材料性能危害极大。表III的数据量化了这种差异在未分裂的案例中间隙原子团簇的比例和最大团簇尺寸都显著增加。特别地当Ta作为PKA且未分裂时产生了超过200个间隙原子的大团簇。4.3.3 物理机制分析为什么V和Nb更容易触发分裂这可以从原子质量和局部环境两方面理解质量效应碰撞截面V和Nb是体系中最轻的两个元素。作为PKA它们具有更高的初始速度。在与周围原子尤其是重原子如W、Ta碰撞时根据碰撞动力学轻粒子更容易发生大角度散射将能量更有效地传递给近邻原子自身则迅速偏转方向。这有利于能量在初始阶段就向多个方向传播而不是沿一个方向深入从而诱发级联分支和分裂。化学与晶格畸变HEA中严重的晶格畸变和化学无序使得原子间的碰撞截面不再是均匀的。PKA在运动路径上会遇到质量、尺寸、键合强度各不相同的原子其运动轨迹更容易被随机扰动和分散进一步促进了级联的分裂。热峰寿命未分裂的紧凑级联会形成一个集中、高温的热峰其寿命较长见图9中W和Ta作为PKA的曲线热峰衰减更慢。较长的热峰寿命虽然增加了点缺陷在热峰内的复合机会降低最终FP总数但也为缺陷在高温下的迁移和聚集提供了时间窗口有利于大团簇的形成。而分裂的级联产生多个小热峰冷却更快缺陷被“冻结”在分散状态。5. 实操、心得与未来展望5.1 分子动力学模拟实操要点与避坑指南基于此项研究进行类似的高熵合金辐照模拟时有几个关键实操点需要特别注意5.1.1 势函数的选择验证是生命线切勿直接套用纯金属势这是最常见的错误。即使有Mo、Nb、Ta、V、W各自的EAM势简单通过混合规则构建的合金势对于HEA往往误差巨大尤其在缺陷能量、表面能等方面。机器学习势是目前的最优解但必须确保其训练数据集覆盖了合金成分空间、各类缺陷、液态及高压短程排斥等关键区域。务必进行系统性验证在跑级联之前必须像本文第三部分那样对势函数进行“体检”。至少包括晶格常数、弹性常数、空位/间隙形成能分布、熔点或至少是熔化趋势。与可靠的DFT或实验数据对比。对于辐照研究离位阈值能量的验证尤为重要它直接关系到损伤产量的基准。关注ZBL势的平滑过渡高能级联模拟必须考虑短程排斥势。确保NEP-ZBL或类似方案在1-2埃的过渡区间内平滑衔接。检查 dimer curve二聚体能量-距离曲线是否平滑、无异常震荡。5.1.2 模拟设置中的细节决定成败PKA方向选择为避免沟道效应channeling导致PKA长距离运动、低估损伤应选择高指数方向如135。本文即采用此方向。在实际研究中若想全面评估可能需要采样多个随机方向。系统尺寸与边界条件盒子必须足够大以确保级联产生的热峰和压力波不会受到周期性镜像的干扰。一个经验法则是盒子边长至少是级联损伤区域最大尺寸的2-3倍。边界区域通常3-5个晶格常数厚度应采用热浴控温以耗散级联产生的热量模拟无限大基体的热沉作用。电子阻止本领对于能量高于10 eV的原子必须考虑电子能量损失电子阻止。本文使用了SRIM数据。忽略此项会导致高能PKA射程过长热峰温度过高影响缺陷产量和分布。统计显著性由于级联过程的随机性单次模拟毫无意义。必须在每个能量点、每种PKA类型下进行足够多次通常≥10次本文达50次的独立模拟取统计平均。误差棒标准差是评估结果可靠性的重要指标。5.1.3 缺陷识别与分析技巧Wigner-Seitz 方法识别点缺陷空位和间隙原子最稳健的方法是Wigner-Seitz元胞法。它基于理想晶格位置不依赖于固定的距离截断对于畸变严重的级联区域更可靠。团簇分析截断半径的选择将点缺陷归类为团簇时截断半径的选择需要小心。通常空位团簇的截断半径设在第二和第三近邻距离之间间隙原子团簇的截断半径设在第三和第四近邻之间。需要测试不同截断值对团簇统计结果的影响确保结论的稳健性。可视化是关键使用OVITO等工具可视化级联的峰值损伤状态和最终缺陷构型。直观观察亚级联分裂、熔融区形态、缺陷分布是理解物理机制不可替代的一步。不要只依赖数字统计。5.2 研究的心得体会与延伸思考完成这项研究后我对于利用机器学习势和分子动力学模拟研究复杂材料辐照行为有了几点深刻的体会“化学无序”是一把双刃剑但利大于弊传统观点认为高熵合金中严重的晶格畸变会增加散射可能产生更多点缺陷。我们的模拟证实了这一点FP产量更高。但更重要的发现是这种无序性通过促进亚级联分裂改变了缺陷演化的路径将可能形成的大害大间隙原子团簇转化为多个小害分散的小缺陷或空位团簇。这种“以空间换稳定”的策略是HEA抗辐照的核心机制之一。轻元素PKA的角色被低估了在考虑辐照粒子种类时人们往往关注重离子如自离子、中子撞出的重反冲核。但我们的工作表明在HEA中轻元素作为PKA可能更易触发级联分裂从而间接提升了抗辐照性能。在实际辐照环境中如聚变堆包含D、T、He等多种粒子轻粒子的这种效应需要被纳入考量。损伤模型需要“个性化”NRT-dpa或标准arc-dpa模型用于评估HEA的辐照损伤剂量如dpa值可能会产生显著偏差。本文提出的扩展模型提供了一个考虑元素特性的框架。虽然其参数仍需通过模拟或实验确定但这指明了未来更精确量化多元合金损伤剂量的方向。机器学习势开启了高通量筛选的大门传统势函数开发周期长一个合金体系一个势。而机器学习势一旦建立了可靠的训练流程和数据集可以相对快速地扩展到新的成分。结合自动化模拟流程未来有望对庞大的高熵合金成分空间进行高通量的辐照损伤初步筛选加速抗辐照材料的发现。5.3 未来工作展望本研究聚焦于初级损伤约100 ps内。材料的最终性能退化取决于缺陷在长时间尺度秒、小时、年下的演化。未来的工作可以沿着以下几个方向深入长时间尺度模拟结合动力学蒙特卡洛或加速分子动力学方法研究间隙原子/空位团簇在HEA中的迁移、融合、湮灭过程以及它们与溶质原子的相互作用。化学有序的影响本文假设完全随机固溶体。实际HEA中可能存在短程有序甚至中程有序。研究这种有序性如何影响碰撞概率、PKA路径和亚级联分裂将非常有趣。多粒子辐照与高温效应模拟氢、氦等气体原子的共注入行为以及高温如服役温度下缺陷产生与退火率的竞争。与宏观性能关联将模拟得到的缺陷类型、密度、尺寸分布作为输入参数通过离散位错动力学或晶体塑性有限元等方法预测材料的辐照硬化、脆化等宏观力学性能变化实现从原子到宏观的跨尺度关联。这项工作就像用一台超高分辨率的原子摄像机拍摄下了高熵合金在辐照冲击下的“第一现场”。我们看到化学无序带来的并非混乱而是一种精妙的“耗散”艺术它将毁灭性的集中能量冲击化解为无数可管理的微小扰动。这不仅是材料科学中的一个具体机制更是一种关于如何在复杂体系中实现鲁棒性和抗损伤能力的深刻启示。机器学习势函数作为这架摄像机的核心镜头正让我们看得越来越清晰也必将引领我们在设计面向极端环境的新材料道路上走得更远。