范德华力模拟计算

我第一次做 范德华力模拟计算是算一个有机液体的密度，cutoff 只设了 10 Å 又没开长程色散校正，算出来密度比实验低 4%，导师一句话”你 LJ 尾巴没补”让我重来。那次让我记住，范德华力模拟计算里色散尾巴看着小，积起来却是体相性质的主宰项。

一、为什么 范德华力模拟计算先定 cutoff 与混规

LJ 势到 cutoff 处如果硬截断，会损失尾部能量和压强，必须用 dispcorr = EnerPres（GROMACS）或 LAMMPS 的 pair_modify tail yes 补长程色散。混规上，异种原子 ε/σ 默认用 Lorentz-Berthelot（算术平均 σ、几何平均 ε），但很多力场对关键相互作用给了专门的交叉参数。我做 范德华力模拟计算时，先确认体系有没有特殊混规，再定 cutoff，这两件事决定了密度和相分离能不能算对。

二、核心原理：色散为什么不能截断

LJ 的 r⁻⁶ 项积分到无穷贡献一个有限但不可忽略的尾部能量，对凝聚相压强尤其敏感。关掉 LRC，压强会系统性偏正，等温压缩率和相平衡点跟着偏。我用 NPT 算密度时，LRC 是必选项。关键认知：vdW 还决定了界面张力、吸附层厚度这些”弱但宏观”的量，参数错 5% 在界面上会被放大成 20% 的偏差。

三、关键技术要点：参数与极化

范德华力模拟计算最容易被低估的是极化效应，纯 LJ 对非极性体系够用，但对含卤素、芳香堆积或离子-π 的体系，各向异性色散（如 D3、IP 模型）才能描述准。我做过一组芳香堆积，普通 LJ 算出来的层间距比实验大 0.3 Å，换带各向异性色散的力场后回到合理值。力场选错，堆积和自组装方向全反。

四、实操流程：从参数到性质

我的标准动作：确认混规与交叉参数 → 开 LRC、设一致 cutoff → NPT 平衡取密度 → 验证相行为/界面 → 和实验对标。最容易被砍的是 LRC，但 范德华力模拟计算 的体相可信度全靠它，我默认保留。

补一个实算：某烷烃混合物，关 LRC 密度 0.78 g/cm³，开 LRC 后 0.81，实验 0.81。4% 的密度差足以让后续的扩散系数、黏度对标全偏。如果拿关 LRC 的轨迹去算输运，结论会系统性错。这之后我任何凝聚相模拟，LRC 是开工第一项，不再等密度不对才回头补。

五、常见踩坑：截断不一致与缩放

范德华力模拟计算常见错是 vdW cutoff 和静电 cutoff 设成不一样，导致能量基准漂移；或 1-4 相互作用的 vdW 缩放因子（如 0.5）在力场间不统一，骨架构象被悄悄改掉。另一个坑是粗粒化模型直接套全原子的 ε，尺度完全错。还有人忽略温度对色散的影响，高温下尾巴贡献更大。

六、复盘：弱作用决定强性质

回过头看，范德华力模拟计算真正难的不是打开一个开关，而是承认”最弱的那一项往往管着最宏观的性质”。我接项目默认先问是不是凝聚相/界面，再决定 LRC 和混规怎么配。把 cutoff、LRC、力场三件事做对，密度和相行为才敢拿去对标。被证明有用的，是愿意为 4% 的密度差返工半天的扎实。下次有人给我一个”密度偏低”的结果，我会先问——你的色散尾巴补了吗。