第4章常用力场介绍

4. 常用力场介绍：Lennard-Jones势、EAM势、ReaxFF反应力场、CHARMM/AMBER生物力场

做分子动力学模拟，选力场就像选工具。你想想看，拧螺丝用扳手，切木头用锯子，力场也是一样的道理。不同的材料体系，得用不同的力场来描述。

我个人习惯把力场分成两大类：经典力场和反应力场。经典力场里，原子之间的连接关系是固定的，不能断键也不能成键。反应力场则允许化学键断裂和形成，适合模拟燃烧、催化这类过程。

今天咱们就聊聊四个最常用的力场。嗯，我当年刚入行时，就是被这几个力场搞得晕头转向，后来踩了不少坑才慢慢摸清楚。

核心观点：没有万能的力场，只有最合适的力场。选力场之前，先想清楚你要模拟什么物理过程。

LJ势，说白了就是描述两个中性原子之间"近了排斥、远了吸引"的相互作用。它的数学形式很简单，就两项：一项是排斥项（r⁻¹²），一项是吸引项（r⁻⁶）。

公式长这样：

U(r) = 4ε[(σ/r)¹² - (σ/r)⁶]

其中ε是势阱深度，σ是原子直径。这两个参数决定了相互作用的强度和范围。

我记得刚做模拟那会儿，用LJ势模拟氩气的液气相变。参数调了半天，结果发现临界温度跟实验值差了20%。后来一查，原来是截断半径设得太小了。嗯，这里要注意：LJ势的截断半径一般要取2.5σ以上，不然能量不连续，模拟会出问题。

我的经验：LJ势虽然简单，但用来模拟惰性气体、有机分子间的范德华力，效果相当不错。不过别用它来模拟金属，那会闹笑话的。

EAM势，全称是Embedded Atom Method，嵌入原子法。它跟LJ势最大的区别在于：EAM考虑了原子所处的局部环境。

你想想看，一个原子在金属表面和金属内部，受到的电子云密度是不一样的。EAM势正好能描述这种差异。它的能量由两部分组成：

公式看起来有点复杂：

E_total = Σᵢ Fᵢ(ρᵢ) + ½ Σᵢⱼ φᵢⱼ(rᵢⱼ)

其中F是嵌入函数，ρ是电子密度，φ是对势。

我曾经用EAM势模拟铜纳米线的拉伸过程。结果发现，用不同的EAM参数文件，得到的屈服强度能差出30%。后来我才明白，EAM参数是针对特定温度范围拟合的，不能随便拿来用。

避坑指南：我曾经吃过一次亏——用室温下拟合的EAM参数去模拟高温扩散，结果扩散系数完全不对。记住：EAM参数有温度适用范围，别越界使用。

前面说的LJ和EAM，原子之间的键是固定的。但ReaxFF不一样，它允许化学键断裂和形成。说白了，它用键级（bond order）来描述原子间的连接关系。

ReaxFF的能量项特别多，包括：

我建议初学者别被这么多项吓到。你只需要知道：ReaxFF能模拟燃烧、爆炸、催化这类有化学反应的过程。它的计算量比经典力场大一个数量级，但能换来化学反应的细节。

举个例子，模拟甲烷燃烧：

# ReaxFF 输入文件片段（LAMMPS格式）
pair_style      reaxff NULL
pair_coeff      * * ffield.reax.cho H C O

这里ffield.reax.cho就是力场参数文件，包含了C、H、O三种元素的相互作用参数。

注意：ReaxFF的参数文件是经过大量量子化学计算拟合出来的。不同版本的参数文件，结果可能差异很大。我习惯用之前先跑几个小体系验证一下。

这两个力场，做生物模拟的人应该不陌生。CHARMM和AMBER都是专门为蛋白质、DNA、脂质这类生物大分子设计的。

它们的能量形式很相似，都包含：

能量项	描述	CHARMM	AMBER
键伸缩	化学键的伸缩振动	K_b(b-b₀)²	K_r(r-r₀)²
角弯曲	键角的弯曲	K_θ(θ-θ₀)²	K_θ(θ-θ₀)²
二面角	绕单键的旋转	K_φ[1+cos(nφ-δ)]	K_φ[1+cos(nφ-δ)]
非键	范德华+静电	LJ + Coulomb	LJ + Coulomb

你可能会问，这两个力场有什么区别？说实话，现在它们越来越像了。早期CHARMM更侧重蛋白质，AMBER更侧重核酸。但现在两者都能处理生物大分子体系。

我个人习惯是：做蛋白质模拟用CHARMM36，做DNA模拟用AMBER ff14SB。当然，这只是个人偏好，没有绝对的对错。

小技巧：用CHARMM或AMBER时，水模型的选择很关键。TIP3P是最常用的，但如果你想模拟蛋白质折叠，我建议用TIP4P-Ew，精度更高。

说了这么多，到底怎么选？我整理了一个简单的对照表：

最后说一句：力场参数不是万能的。我见过有人拿着LJ势去模拟金属，结果能量发散；也见过有人用AMBER模拟高温下的蛋白质，结果结构全乱套。选力场之前，先问问自己：我要模拟什么？这个力场能描述我关心的物理过程吗？

嗯，今天就聊到这儿。力场这东西，用多了自然就有感觉了。

公众号:蓝海资料掘金营，微信deep3321