4. Lennard-Jones势的选型判断

Lennard-Jones（LJ）势是分子动力学中使用最广的对势模型，形式简单、计算便宜，对惰性气体和简单分子的范德华相互作用描述得相当好。但“能用”和“该用”是两回事。本节要解决的核心问题是：拿到一个体系时，如何判断LJ势是否合适，参数从哪里来，截断半径怎么定，长程修正要不要加，以及什么情况下必须果断放弃LJ。

一、LJ势的数学形式与两个参数的物理含义

LJ势的标准形式为：

V(r) = 4ε [ (σ/r)^12 - (σ/r)^6 ]

其中 r 是原子间距离，ε 是势阱深度，代表相互作用的强度；σ 是势能为零时的距离，近似反映原子的有效直径。第一项 (σ/r)^12 是排斥项，来自泡利不相容原理导致的电子云重叠排斥；第二项 (σ/r)^6 是吸引项，来自瞬时偶极—诱导偶极的色散力。注意：LJ势只有两个可调参数，这意味着它把所有范德华相互作用的复杂性压缩成两个数，这是它便宜的原因，也是它局限的根源。

在LAMMPS中，LJ势的典型写法是：

pair_style lj/cut 10.0
pair_coeff 1 1 0.0103 3.40
pair_coeff 1 2 0.0080 3.20
pair_coeff 2 2 0.0060 3.00

这里 10.0 是全局截断半径（单位Å），每行的两个数分别是 ε（eV）和 σ（Å）。不同原子对需要分别指定，混合规则默认用几何平均求 ε、算术平均求 σ。

二、参数来源：三条路，各有适用场景

LJ参数不是随便填的，来源主要有三类。

第一类：实验数据拟合。最经典的是惰性气体。例如氩（Ar）的 ε ≈ 0.0104 eV，σ ≈ 3.40 Å，这组参数来自第二维里系数和黏度数据的拟合，对液态氩的密度、扩散系数复现得很好。氖、氪、氙也有类似的高质量参数。这类体系的LJ势是“物理正确”的，可以放心用。

第二类：从头算或高精度参考数据拟合。对于小分子或离子，可以用CCSD(T)等量子化学方法算出二聚体相互作用能曲线，再拟合出 ε 和 σ。这类参数比实验拟合更灵活，能处理实验难以测量的体系，但拟合质量取决于量子化学计算的精度和采样点的覆盖范围。

第三类：力场数据库中的现成参数。OPLS、TraPPE、COMPASS等力场都内置了LJ参数。使用时要看清参数是为哪种体系、哪种截断方案拟合的。例如TraPPE的参数通常配合1.4 nm截断和长程修正使用，如果你改成2.5 nm截断又不加修正，结果会偏。

三、截断半径：不是越大越好，但太小一定出问题

LJ势的吸引项按 r^-6 衰减，理论上要到无穷远才为零。实际计算必须截断，截断半径 rc 的选择直接影响计算量和精度。

经验规则是 rc ≥ 2.5σ。对氩来说，σ = 3.40 Å，2.5σ = 8.5 Å，所以常用 10 Å 或 12 Å。如果 rc 太小，比如只取 2.0σ，截断处的势能还有明显数值，原子在截断边界会感受到“突然消失”的力，导致能量不守恒、密度偏低。

截断方式也有讲究。最简单的 truncation 是直接把 rc 以外的势能设为零，这会在 rc 处产生势能跳跃，力不连续。更常用的是 shift 方法，把整个势能曲线平移，使 V(rc) = 0，能量连续但力仍不连续。最平滑的是 force-switch 或 potential-switch，在 rc 附近的一个窗口内把力和势能逐渐过渡到零，适合对能量守恒要求高的体系。

# 平滑截断示例：从9.0 Å开始过渡，10.0 Å处完全截断
pair_style lj/cut 10.0
pair_modify shift yes

四、长程修正：截断丢掉的能量要补回来

截断不仅丢掉了 rc 以外的相互作用，还改变了体系的内压和能量。长程修正（long-range correction, LRC）就是假设 rc 以外体系是均匀的，用积分把平均贡献补回来。

对能量，修正项为：

U_LRC = (8πNρ/3) ε σ^3 [ (1/3)(σ/rc)^9 - (σ/rc)^3 ]

对压力也有类似修正。在LAMMPS中，只要在 pair_style 后加 pair_modify tail yes，就会自动加上tail修正。注意：tail修正假设 rc 以外是均匀流体，对气液界面、团簇、表面体系不适用，因为这些区域在 rc 尺度上并不均匀。

五、什么时候必须放弃LJ

LJ势的适用边界很清晰，以下情况必须换模型：

• 有方向性相互作用：水、醇、酰胺等体系存在氢键，氢键有明确的方向性和饱和性，LJ势是各向同性的球对称模型，无法描述。必须用带点电荷和氢键项的模型，如TIP3P、SPC/E或更复杂的极化力场。
• 有静电相互作用：离子、带电基团、极性分子之间的库仑力远大于色散力，LJ势完全缺失这一项。必须用 coulomb 加 LJ 的组合，并配合PPPM等长程静电求解器。
• 需要描述化学键的形成与断裂：LJ势是物理吸附模型，不能描述共价键的断裂和重组。反应体系必须用ReaxFF、AIREBO等反应力场，或走QM/MM路线。
• 需要高精度能量和结构：金属体系用LJ势会得到错误的配位数和熔点，必须用EAM、MEAM等嵌入原子方法。半导体和碳材料需要Stillinger-Weber、Tersoff或机器学习势。
• 需要极化效应：强电场、界面、离子通道等场景中，电子云会重新分布，固定电荷模型不够，需要Drude振子、 fluctuating charge或极化力场。

下面这张图把LJ势的选型判断整理成一条决策路径，从体系特征出发，逐步判断是否继续使用LJ，以及需要补充或替换什么。

回到实践层面，一个典型的判断流程是：先看体系有没有氢键、静电、反应、金属键——有就放弃纯LJ；没有，再确认参数来源是否与力场一致；然后按 2.5σ 规则定截断半径，均匀体系加tail修正，非均匀体系不加；最后跑一个短程NVE检查能量守恒，如果总能量漂移超过 10^-4 量级，优先怀疑截断设置而不是时间步长。

LJ势不是“低级”模型，它是范德华相互作用的正确一阶近似。问题从来不在LJ本身，而在于把它用在了它不该出现的地方。选型判断的核心，就是先问体系里有没有LJ描述不了的物理，再问参数和截断是否自洽，最后用能量守恒和结构量做一次快速验证。