范德华力理论计算

我第一次做 范德华力理论计算是算石墨层间结合能，用裸 PBE 跑出来层间作用只有 0.02 eV，文献值却是 0.05 eV 量级，还差了一倍。后来加上 DFT-D2，数字立刻回到合理范围。那次让我记住，范德华力理论计算里色散不是锦上添花，是主导项，关掉它你算的根本不是同一个物理。

一、为什么 范德华力理论计算先做泛函选型

标准 LDA/GGA（如 PBE）在长程色散上天生缺失，对层状、堆积、主客体这类靠 vdW 稳定的体系系统性低估。解决路线有三类：加经验色散校正（DFT-D2/D3/D4）、用 vdW-DF 系列非局域泛函、或用杂化+DFT-D。我默认小体系用 PBE-D3(BJ)，层状材料用 optB88-vdW 或 rVV10，因为不同方案的层间距离和结合能差能到 20%。做 范德华力理论计算时，泛函选型比 k 点更该先定。

二、核心原理：色散项怎么加进去

DFT-D 是在 DFT 总能上额外加一个 −C₆R⁻⁶ 形式的经验项，D3 还引入三体 AXB 项、用可调阻尼函数避免近程双计。vdW-DF 则是把非局域关联能直接写进泛函，对层间距离更敏感。我用 D3 时习惯看它给出的 pairwise 贡献，确认色散项量级和物理预期一致。关键认知：色散校正的 C₆ 系数是元素相关的，主族和过渡金属不能混用同一套，否则层间力算飞。

三、关键技术要点：基准对标不能少

范德华力理论计算最怕”盲加 D3 就交差”。我每接一个新体系，先用一个已知基准（比如 S22 复合物、层状材料实验层间距）校验所选方案的误差。曾经用 D3 算一个金属有机框架的客分子结合，误差 0.1 eV，换成 rVV10 后降到 0.03 eV——原因是骨架里有大量色散主导的弱作用，D3 的经验参数覆盖不到。泛函选错，后面所有结合能都建在沙子上。

四、实操流程：从建模到结合能

我的标准动作：选基准体系校验泛函 → 结构弛豫（色散开着的条件下）→ 分层扫描层间距或结合曲线 → 取能量最低点算结合能 → 和实验或高阶方法（如 CCSD(T) 小体系）对标。最容易被砍的是基准校验，但 范德华力理论计算 的可信度全靠它，我默认保留。

补一个主客体识别的实算：杯芳烃包结一个客分子，PBE 算结合能 0.08 eV，PBE-D3 升到 0.21 eV，和实验量热 0.19 eV 接近。如果只用裸 PBE，你会得出”这个主客体几乎不结合”的错误结论，进而误判分离选择性。这之后我做弱作用体系，色散方案是开工第一项，不再等结果不对才回头补。

五、常见踩坑：近程双计与阻尼

范德华力理论计算最常见的坑是 DFT-D 的阻尼函数在近程和 DFT 自身的交换关联重叠，导致短程能量双计、键长缩短。我见过有人用 D2 不调阻尼，氢键体系算出来的几何明显偏紧。另一个坑是 vdW-DF 和 DFT-D 混用，两套色散加两次，结果毫无意义。还有人忽略范德华对频率的影响，算出来的振动模式和实验对不上。

六、复盘：弱作用才是很多体系的命门

回过头看，范德华力理论计算真正难的不是加个开关，而是选对能描述你体系相互作用层次的方案。我接项目默认先判断”这个体系靠不靠 vdW 稳定”，再决定走 D3 还是 vdW-DF。把基准校验、泛函选型两件事做对，结合能才敢拿去解释选择性。被证明有用的，是愿意为一个层间距多扫五个点的扎实。下次有人给我一个”结合很弱”的结论，我会先问——你的色散开没开。