分子动力学模拟力场

力场是分子动力学的灵魂，决定了整个模拟”算的是什么物理”。我这个项目我不是跑一个具体反应，而是系统梳理并验证了一套全原子力场，目标让客户真正理解 分子动力学模拟力场 里每一项势能项的物理意义和参数来源，免得以后调参像盲人摸象。

一个典型全原子力场包含四个部分。第一，键伸缩（bond）：谐振子势 U = ½k_b(r − r₀)²，描述成键原子间的距离涨落，k_b 是力常数、r₀ 是平衡键长。第二，键角弯曲（angle）：同样谐振子形式 U = ½k_θ(θ − θ₀)²。第三，二面角扭转（dihedral）：周期势 U = k_φ[1 + cos(nφ − γ)]，n 是周期数、γ 是相位，控制分子围绕化学键的旋转自由度。第四，非键相互作用：Lennard-Jones 势描述范德华（吸引+排斥），Coulomb 势描述静电，这是力场里最耗算时也最容易出问题的部分。

参数来源是力场的命门。我梳理了常用力场：OPLS-AA、AMBER、CHARMM 对小分子、蛋白、脂质都有成熟参数；金属/无机材料则要换 ReaxFF 或 MEAM 这类反应力场。我这个项目验证时跑了一个小分子在真空中的 NVT 模拟，看能量守恒和温度漂移。结果发现：1 fs 步长下能量漂移很小（< 0.1%），但换成 2 fs 后漂移明显变大——原因是含氢键的振动频率高，2 fs 步长会让积分失真。解决方法是加 SHAKE/RATTLE 约束把氢相关的键长冻住，之后 2 fs 也能稳。这个细节很多新手会忽略，导致模拟”跑着跑着就炸了”。

另一个关键坑是 1-4 非键相互作用的 scaling。很多力场对相隔 3 个键（1-4）的原子，其 LJ 和静电会乘一个系数（常见 0.5 或 0.833），写输入文件时必须和力场文档对上，否则体系能量会系统偏移、构象偏掉。我特意做了一个对照：用错 scaling 时，小分子的二面角分布明显畸变。

验证清单我给客户列了四项：①能量守恒（NVE 下总能量是否常数）；②温度漂移（NVT 下是否围绕设定值波动）；③结构合理性（RMSD 是否稳定、键长分布在合理范围）；④与实验对照（如密度、扩散系数误差 < 10%）。只有四项全过，力场才算”可用”。

踩坑点：第一，力场-水模型匹配。AMBER 力场要配 TIP3P 水、CHARMM 配 TIP3P 也有差异，混用会破坏溶剂化结构。第二，参数外推。通用力场对含氟化合物、离子液体等特殊体系不准，需要专门参数化，我这个项目的小分子恰好含氟，我用了针对性的 RESP 电荷。第三，长程静电。PME 的截断半径和网格设置不当会引入静电误差，我检查了力的收敛。第四，约束与积分器。不同约束算法（SHAKE vs LINCS）对步长容忍度不同，须匹配。

从我的工程经验看，分子动力学模拟力场 最忌”拿来就用”。我给客户的建议是：选力场先看体系类型，验证四步一个不能少，含氢体系务必加约束。更多力场选择与验证，见 [分子动力学模拟力场](https://www.keyanxueshu.com/category/md/)；需要力场参数化，[分子动力学模拟力场](https://www.keyanxueshu.com/) 上有服务入口。