范德华修正全解：从色散相互作用到方法差异与结果判据

说明：本文华算科技主要介绍范德华修正补充的色散相互作用、常用修正方法的物理差别，以及能量、结构和参照数据的判断条件。

范德华修正补充的色散能是什么？

两个有机分子结合时，总相互作用能含有静电、诱导、交换排斥和色散等贡献。材料计算中的范德华修正通常主要补充伦敦色散作用：不同片段的电子密度量子涨落通过库仑作用发生关联，使基态能量降低。没有永久偶极的闭壳层片段也存在这种吸引，零温下同样存在；它的来源区别于原子热振动。

精确的密度泛函理论包含色散，描述不足发生在常用的交换关联近似中。PBE 等半局域泛函用局部密度及其梯度构造能量，难以恢复相隔较远片段之间的涨落关联。半局域泛函在平衡接触附近可能已经包含部分吸引，却缺少正确的长程色散行为。混合泛函引入部分精确交换，也没有自动补齐这部分动态电子关联。

DMASQ 方酸染料二聚体既可错位堆积，也可形成氢键接触。图中的 a、b、c 采用含 D3 的方法优化，d 来自未加 D3 的优化。色散参与分子接触的能量竞争，同一分子的局部极小值可能随方法改变。若以各自优化的单体为参照，平行堆积与氢键构型的结合能差还包含静电吸引、排斥及单体形变的贡献。

对相距较远、尺寸有限的两个中性片段，在非迟滞条件下，领先色散项可写成 E_disp(R) ≈ −C₆/R⁶。C₆取决于片段的动态极化响应，R 是片段间距离；C₆ 为正时，该项给出吸引。原子对求和还可含 C₈/R⁸ 等高阶项。片段的尺寸与距离可比、电子云明显重叠或存在延展导电响应时，这个简单远距表达式的适用范围会缩小。

范德华力在分子间作用分类中还可包含永久偶极的取向作用和诱导作用，而DFT 中的色散修正有特定能量对象。常规自洽计算已经描述电荷分布产生的静电作用及一部分极化响应；添加色散时应与原有交换关联近似配合。一次修正的数值反映所选模型的色散贡献，整个分子结合能还含其他能量分量。

辨析成对修正、非局域泛函与多体色散

DFT-D 将色散项加入基础泛函的总能，E_DFT+D = E_DFT + E_disp。有限体系的领先成对项可示意为 E_disp = −s₆Σ_A f₆(R_AB)C_6,AB/R_AB⁶，s₆ 是与泛函配套的尺度因子。这里的原子对距离、色散系数与阻尼函数共同决定修正量；周期晶体还需计入周期镜像的贡献，修正随结构变化。

D2 主要采用按元素给定的系数，D3 用局部配位环境调整系数，并计入更高阶成对项；D4 又引入原子部分电荷对极化率的影响。常用 D4 电荷来自电荷平衡模型，和 Bader 分区电荷的定义不同。TS 方法则从电子密度分区获得有效原子体积，据此缩放自由原子的极化率与色散系数。上述方法对同一元素处在不同化学环境中的响应采用了不同近似。

短程阻尼限定远距公式向接触区的延伸，抑制无阻尼 R⁻⁶ 项在小距离处的发散，并调节与基础泛函已有相关能的衔接。D3 的零阻尼与 BJ 阻尼采用不同形式，接触附近的能量与导数会有差别；阻尼参数也按基础泛函配套。下方 PTCDA 势曲线显示，渐近势在吸附高度附近已经偏离短距估计，远距系数与接触区势能需要分别检验。

NTCDA、PTCDA、TTCDA 在 Au(111) 上的力测量表明，每碳原子的有效 C₃ 随分子尺寸增加，成对表面修正给出的尺寸趋势较弱，RPA 计算则呈现较强的尺寸依赖。C₃ 是分子—表面渐近势 V(z) ≈ −C₃/(z−z₀)³ 的系数，z₀ 为表面参考平面；它与原子对 C₆ 对应不同几何下的距离依赖。动态极化率面板区分长轴、短轴和面外响应，其虚频轴用于表示电子响应函数，区别于晶格振动中的虚频。

vdW-DF 和 VV10/rVV10 一类方法把非局域相关能写入泛函，示意式为 E_c^nl = ½∫∫ n(r)Φ(r,r′)n(r′)d³r d³r′。核函数 Φ 连接两个空间位置的电子密度，并依赖距离和局部密度信息。自洽实现中，该项可通过泛函导数影响电子势；最终结果由非局域相关项与配套交换、短程相关共同决定。

MBD 用相互耦合的量子德鲁德振子表示原子极化响应，通过耦合前后模式的零点能差描述多体色散。一个原子的涨落会受到其余原子的响应影响，屏蔽和集体响应可削弱或增强成对估计，方向随结构与环境改变。D3、D4 可包含三体项，但有限阶三体修正与 MBD 的集体响应处理仍有差别。

C₇₀—[8]-CPPA 复合物中的模型图把原子振子转换为集体模式。右上颜色表示振子的平均激发数，右下矩阵给出激发数之间的协方差，反映电子响应模型中的关联。这些振子模拟电子极化，图中位移及激发数没有直接表示晶格声子或原子得失电子数。非局域密度核与 MBD 的响应模型都处理色散，但采用的变量和近似不同。

怎样检验范德华修正对计算结果的影响？

在固定二聚体几何下，相互作用能可定义为 E_int = E_AB − E_A − E_B，单体采用复合物中的几何。复合物与两个单体应采用相同的色散处理，修正对相互作用能的贡献也由这三个能量相减得到。若单体改用各自优化结构，能量差还含单体形变代价，比较的物理过程随之改变。

DMASQ 的能量图横轴列出四种离散构型，纵轴为总相互作用能。不同基组和电子关联方法给出的点值及排序存在差别，图中的连线用于连接构型。色散分量与总相互作用能要分别比较：某种接触具有更负的色散能时，静电作用和交换排斥仍可能改变总能排序。

距离扫描中的总势能最低点给出平衡距离，曲线斜率决定受力；色散项也通过 F_disp = −∂E_disp/∂R 参与结构优化。固定几何的单点修正与重新优化后的能量回答不同问题。后者包含吸附高度、分子取向和内部结构的调整。常见 D3 在固定结构下主要补充能量与力，其对能带的影响常通过几何变化传递；非局域泛函的自洽密度响应则另有贡献。

SAPT 把相互作用能按静电、交换、诱导和色散等分量分解。DMASQ 的色散分量随基组和弥散函数改变，堆积构型的差别尤其突出。SAPT0、SCS-SAPT0 与 DFT-SAPT 也采用不同电子关联近似。这些变化属于参照方法的模型与基组误差，比较 D3 或 D4 的效果时，参照值本身的误差必须有依据。

实验吸附高度、解离能、晶体体积分别约束势能曲线的不同特征。零温电子能与实验热力学量之间的比较涉及零点能、温度、熵以及实际覆盖度；分子表面的渐近力系数也只约束相应距离范围。层状材料的层间距符合实验时，势阱深度仍需独立核对。若基础泛函存在电荷过度离域、磁态或局域电子描述偏差，补充色散后的总能仍保留相应误差。

对于能量相近的堆积构型，基组、周期模型尺寸与修正方案造成的相对能量变化要与构型间能差比较。原子受力收敛、周期镜像贡献和结构约束应与该比较一致；方法差别若与构型能差相当，当前计算仅能给出相近能量的候选构型，排序尚未分辨。