1. 势函数选型总论

在分子动力学（MD）中，原子势函数（interatomic potential）是唯一决定“原子之间如何相互作用”的物理模型。它把原子坐标映射为势能，再通过势能对坐标的负梯度给出每个原子所受的力：F_i = −∇_iU(r₁, r₂, …, r_N)。积分器随后用这些力更新速度与位置，统计力学再把轨迹转化为温度、压力、扩散系数、弹性常数等宏观量。因此，势函数不是“一个可随意替换的数值参数”，而是整个模拟的物理地基：地基选错，后续再长的采样、再精细的后处理都无法恢复正确的物理。

势函数选型本质上是一次“物理问题 → 模型能力”的匹配。不同势函数在键合描述、长程作用、电荷转移、磁性、激发态等方面能力差异极大。例如，Lennard-Jones 只能描述范德华与简单密堆积；嵌入原子法（EAM）适合金属；Tersoff/Brenner 类键序势适合共价体系；ReaxFF 允许键断裂与形成；而机器学习势（如基于 DFT 数据的 NEP、MTP、DP）在精度与体系适应性上更强，但依赖训练集覆盖度。选型的核心不是“哪个势最好”，而是“哪个势的近似假设与我的科学问题相容”。

选型失误的后果往往具有隐蔽性。第一类是物理机制缺失：用固定电荷的离子势去模拟质子转移，体系根本不会发生反应，得到的“反应路径”只是原子被强行推过势垒的假象。第二类是参数外推失效：用 300 K 下拟合的势去模拟 2000 K 的熔化，势函数未覆盖的高能构型会导致密度、熔点、潜热系统性偏差。第三类是量纲与单位混用：LAMMPS 的 metal 单位下能量为 eV、距离为 Å，而某些势文件以 kcal/mol 或 nm 给出参数，若未换算，力常数会偏差数十倍，表现为体系瞬间“爆炸”或冻结。第四类是截断与长程处理不当：对离子晶体忽略库仑长程求和，或对高分子把截断半径设得小于链段相关长度，会引入虚假的表面与密度振荡。

因此，选型必须与验证构成闭环，而不是“选完就跑”。一个可操作的闭环包含四个环节：需求定义 → 候选筛选 → 基准验证 → 生产模拟与回归。需求定义阶段要写清目标性质、温度压力区间、体系尺寸与时间尺度；候选筛选阶段查阅文献与势数据库（如 NIST Interatomic Potentials Repository、OpenKIM），优先选择已有独立验证的势；基准验证阶段用小型体系做“可证伪”的测试；生产模拟阶段持续监控守恒量与异常信号，并在条件变化时回到验证环节。

基准验证应至少覆盖三类指标。其一是结构量：晶格常数、径向分布函数第一峰位置、配位数。其二是能量与力学量：内聚能、弹性常数、空位形成能、表面能。其三是动力学量：均方根位移给出的扩散系数、声子色散或热导。以面心立方铜为例，若用 EAM 势，可先做 0 K 弛豫得到晶格常数 a；实验值约 3.615 Å，合格势通常落在 3.60–3.63 Å。再做 300 K、NPT 系综下 20 ps 的平衡，检查密度是否在 8.9–9.0 g/cm³、势能是否收敛、温度是否稳定在 300±5 K。若这些基础量偏差超过 2%–3%，就不应直接进入生产模拟。

验证不是一次性的。当模拟条件改变（温度、压力、界面、缺陷浓度）或体系规模扩大时，势函数的适用边界可能被跨越。此时应回到闭环，补充针对新条件的验证，例如高温下的熔点、界面体系中的粘附功、缺陷体系中的迁移势垒。只有把“选型—验证—生产—再验证”作为持续过程，才能保证模拟结论可信。