3. 势函数选择与评估

在晶界与位错的分子动力学模拟中，势函数是唯一决定原子间相互作用的核心输入。它既定义了体系的能量景观，也决定了位错核心结构、层错能、晶界能和迁移势垒等关键量。因此，势函数选择不是“随便挑一个能跑”的步骤，而是决定模拟结论是否可信的第一道关口。本节围绕 EAM、MEAM、ADP 三类常用势的适用边界，以及一套可操作的验证与标定流程展开。

三类势函数的能力边界

嵌入原子法（EAM）把体系能量写成嵌入能加对势的形式，嵌入能由局域电子密度决定。它的优势是计算量小、对 fcc 金属（Cu、Al、Ni、Ag、Au）的弹性常数、层错能、空位形成能拟合较好，因此长期是晶界与位错模拟的主力。但 EAM 的电子密度是球对称的，无法描述方向性键合，对 bcc 金属（Fe、W、Mo）的位错核心（尤其是 screw 位错的孪晶/滑移竞争）和晶界偏聚往往偏差明显。

MEAM（Modified EAM）在 EAM 基础上引入角度依赖的电子密度修正，能更好地描述共价性和方向性键合，适合 Si、bcc 过渡金属以及含 C、H 等轻元素的体系。代价是参数更多、拟合更复杂，且不同版本（MEAM、2NN-MEAM、MEAM/C）对同一元素的预测可能差异较大。ADP（Angular-Dependent Potential）则在 EAM 中显式加入偶极张量项，保留 EAM 的高效性同时改善堆垛层错能、孪晶能和位错核心的角向描述，对 hcp 金属（Mg、Ti、Zr）和部分 fcc 体系表现突出。

验证与标定的实操流程

拿到一个候选势文件后，不要直接投入大体系模拟。建议按以下顺序做小体系验证，每一步都给出可判断的数值依据。

第一步，检查基本物理量。用 0 K 弛豫后的完美晶格计算晶格常数 a₀、内聚能 E_coh 和弹性常数。以 fcc Cu 为例，EAM 势通常给出 a₀≈3.61 Å、E_coh≈-3.54 eV/atom、C44≈75 GPa。若 a₀ 偏差超过 1%，后续所有缺陷能都会系统性偏移，应优先排查势文件与元素映射是否匹配。

第二步，计算层错能。沿 [111] 方向施加刚性位移，构造内禀层错（ISF）和孪晶层错（TSF），得到 γ_isf 与 γ_tsf。这两个量直接决定位错分解宽度和孪晶倾向。例如 fcc Cu 的 γ_isf 参考值约 45 mJ/m²，若某 EAM 势给出 25 mJ/m²，则位错会过度分解，模拟出的孪晶行为会明显偏软。

第三步，计算晶界能。以 Σ5(310)[001] 对称倾斜晶界为例，构建双晶模型，固定一侧原子、弛豫另一侧，得到 γ_GB。fcc Cu 的该晶界能参考值约 0.7～0.9 J/m²。若计算值偏离 20% 以上，说明势对晶界结构的描述不可靠，需谨慎用于偏聚或断裂模拟。

第四步，检查位错核心。用差分位移图（differential displacement map）观察 screw 位错核心是否分解、分解宽度是否合理。bcc Fe 的 screw 位错核心在 EAM 下常呈非简并的紧凑结构，而 DFT 预测其呈简并或部分分解，这是 EAM 用于 bcc 时最典型的失效点。

第五步，热稳定性与熔点。用 NPT 系综升温，观察体系在熔点附近是否出现异常膨胀或原子重叠。若熔点偏离实验值超过 15%，说明势的高温行为不可靠，涉及高温晶界迁移的模拟需重新评估。

常见误区与判断依据

一个常见误区是“用文献里最常被引用的势就一定安全”。实际上，同一元素的不同势版本可能针对不同性质拟合：有的侧重弹性常数，有的侧重缺陷能，有的侧重液相结构。选择时应先明确本节模拟最关心的输出量，再反向筛选在该量上有验证数据的势。

另一个误区是忽略截断半径与近邻数的一致性。EAM 的电子密度截断若设置过小，会导致晶界处原子配位数变化时能量描述失真；MEAM 的 2NN 修正若未启用，bcc 的剪切行为会明显偏硬。这些细节应在势文件说明中确认，而不是默认“能跑就行”。

最后，验证不是一次性动作。当模拟条件从 0 K 静态计算切换到有限温度、从单晶切换到多晶、从纯金属切换到合金时，都应重新检查关键量是否仍在可接受范围内。只有经过验证与标定的势函数，才能支撑晶界与位错行为的可靠结论。