4. 原子间势函数工程选择

在辐照损伤分子动力学模拟中，原子间势函数是决定模拟可信度的第一性输入。它定义了原子在给定相对位置下的势能，进而决定力、缺陷形成能、迁移势垒、级联演化路径乃至最终损伤形貌。同一材料、同一辐照条件，换一套势函数，结果可能从“形成致密位错环”变成“缺陷大量复合”，因此势函数选择不是软件参数问题，而是物理建模问题。本章围绕嵌入原子方法（EAM）、Finnis-Sinclair势、Tersoff势、机器学习势，以及势函数验证与迁移性评估、核材料常用势函数库展开。

一、从对势到多体势：为什么金属需要EAM

最简单的对势只依赖两原子距离，无法正确描述金属的弹性常数与空位形成能之间的Cauchy关系。EAM把体系能量写成嵌入能加对势排斥项：

E_total = Σ_i F_i(ρ_i) + 1/2 Σ_{i≠j} φ_ij(r_ij)
ρ_i = Σ_{j≠i} f_j(r_ij)

其中ρ_i是i原子处的电子密度，F是嵌入能，φ是短程排斥对势。嵌入能项使能量对局部配位环境敏感，因此EAM能同时给出合理的空位形成能、表面能和弹性常数。对辐照损伤而言，EAM的关键优势是能较稳定地描述Frenkel对、间隙原子和空位团簇的相对稳定性。工程上选择EAM势时，应优先检查其是否针对缺陷性质拟合，而非只看晶格常数和结合能。

二、Finnis-Sinclair势：EAM的紧束缚近似来源

Finnis-Sinclair势可视为EAM的一种早期形式，其嵌入能取平方根形式，来源于紧束缚近似下d带电子密度的二阶矩。它形式简洁、计算高效，在bcc金属（如Fe、W、Mo）中应用广泛。FS势的典型特点是能较好再现bcc金属的弹性各向异性，但对<111>间隙原子的迁移势垒和位错芯结构往往需要额外修正。实际使用中，若模拟对象是bcc铁中的级联损伤，应确认势函数是否给出了合理的<110>哑铃间隙构型与迁移能，否则级联后期的缺陷团簇统计会失真。

三、Tersoff势：共价键材料的定向成键描述

对于Si、SiC、GaN等共价材料，EAM的电子密度近似不再适用。Tersoff势引入键序概念，使键强依赖局部配位和键角：

E = Σ_i Σ_{j>i} f_C(r_ij)[f_R(r_ij) + b_ij f_A(r_ij)]
b_ij = (1 + β^n ζ_ij^n)^(-1/2n)

其中ζ_ij由周围原子对键角的贡献累加。Tersoff势能描述sp3杂化、键角畸变和缺陷重构，适合SiC、Si等辐照损伤模拟。其局限是参数拟合范围有限，对高能级联中的极端短程碰撞可能不准确，通常需要与ZBL短程修正拼接。

四、机器学习势：精度与迁移性的新平衡

机器学习势（MLP）用高维描述符或图神经网络拟合第一性原理数据，能在保持近DFT精度的同时达到经典势的计算速度。常见路线包括Behler-Parrinello对称函数、GAP/SOAP、以及基于图神经网络的NequIP、MACE等。对核材料而言，MLP的吸引力在于可同时覆盖多种缺陷构型和合金组分，但其外推能力依赖训练集覆盖度。若训练集只含平衡态和少量缺陷，遇到级联热峰中的高能原子时可能给出非物理能量。因此MLP必须做迁移性评估，不能默认“精度高就通用”。

五、势函数验证：不能只看晶格常数

拿到一套势函数后，至少应做四类验证。第一，基本性质：晶格常数a₀、结合能E_c、弹性常数C₁₁、C₁₂、C₄₄，与实验或DFT偏差一般应控制在几个百分点内。第二，点缺陷性质：空位形成能E_v、间隙形成能E_i、迁移势垒E_m。以bcc Fe为例，实验空位形成能约1.7～2.0 eV，<110>哑铃间隙迁移能约0.3 eV；若某EAM势给出E_m超过0.6 eV，则级联在室温下会过度冻结，缺陷复合率被高估。第三，缺陷团簇：位错环形成能、堆垛层错能、氦-空位结合能。第四，热力学：熔点、热膨胀系数、声子谱。验证时建议固定同一初始构型和积分步长，只换势函数，做对照计算。

六、迁移性评估：势函数能走多远

迁移性指势函数在训练或拟合范围之外是否仍给出物理合理结果。评估方法包括：温度迁移性——从300 K外推到1000 K以上，检查是否出现非物理熔化或相变；组分迁移性——从纯金属外推到合金，检查混合焓和偏析能；缺陷迁移性——从单空位外推到空位团簇、位错环和层错；能量迁移性——从热平衡外推到级联热峰中的高能碰撞。对MLP，还应做不确定性量化，识别描述符空间中的外推区域。若模拟对象是聚变堆第一壁钨在14 MeV中子下的级联，必须确认势函数在PKA能量达到数十keV时仍与ZBL短程势平滑衔接。

七、核材料常用势函数库

核材料领域常用的势函数库包括：Fe-He EAM势（如Juslin、Bonny等版本），用于bcc铁中氦泡和级联损伤；W-Re Finnis-Sinclair势，用于钨基面向等离子体材料；SiC Tersoff势，用于包壳和共价陶瓷；U-Zr、U-Mo等金属燃料的EAM或MLP；UO₂的Buckingham、Rigid Ion或MLP势。选择时需注意版本差异：同一元素不同作者拟合的势，缺陷性质可能相差20%以上。建议优先选用在核材料辐照损伤文献中被反复验证、且公开了完整参数和验证数据的势。

总结而言，势函数工程选择的核心逻辑是：先判断键合类型确定势函数族，再用缺陷性质而非仅弹性常数做验证，最后通过温度、组分、缺陷和能量四个维度的迁移性评估界定适用范围。只有把势函数当作有明确适用边界的物理模型，而不是黑箱参数，辐照损伤模拟结果才具备可解释性和可复现性。