说明：本文华算科技主要介绍范德华修正在第一性原理计算中的物理来源、能量写法、常见方法和结果判读边界。
范德华修正的物理来源是什么？
在 Kohn-Sham DFT 里，交换相关泛函决定电子密度怎样转成能量。LDA、GGA 和许多 meta-GGA 主要使用局域密度、密度梯度或动能密度来近似 Exc，对共价键、离子键、金属键一类短程成键描述较顺手。相隔较远的中性片段之间还存在一种吸引：瞬时电子密度涨落会诱导另一片段产生相关涨落，形成色散相互作用。范德华修正处理的正是这类长程相关吸引在常规半局域泛函中的缺口，该吸引在距离较远时仍会贡献能量。
色散能常用渐近形式写成 -C6/R6距离依赖吸引尾部，也可包含 C8/R8、C10/R10 等高阶项。这里的 R 是片段或原子之间的距离，C6 与极化率、局域化程度、化学环境有关。距离很短时，普通交换相关泛函已经包含部分相关能，色散项还需配合阻尼函数，避免同一段吸引被重复计入。修正项并非给所有能量统一加常数，它会改写势能曲线的形状，平衡距离、结合能和力都可能跟着变化。势能井的深度和最低点位置才是结构响应的来源。
图1. 苯和草酸晶体的晶胞体积-晶格能曲线显示，加入 D3 色散修正会移动能量最低点和势能井深度。DOI：10.1103/PhysRevB.94.115144
图 1 把苯和草酸晶体的晶胞体积放在横轴上，曲线最低点对应优化结构附近的平衡体积，箭头标出加入 D3 色散修正后能量最低点的移动方向。没有长程吸引时，势能井可能偏浅、平衡体积偏大；加入色散后，分子堆积会收紧，晶格能和体积一起变化。势能井深度或位置改变，结构优化得到的距离就会改变。
材料计算中常见的“加不加 vdW”争论，真正讨论的对象是能量面。吸附分子离表面太远时，短程电子云重叠小，半局域泛函容易漏掉吸引尾部；距离过短时，Pauli 排斥和化学成键主导，色散项只占一部分。范德华修正的合理读法，是看它怎样改变平衡距离附近的能量、力和相对构型，而不是把计算结果统一贴上“更准”标签。
色散参数怎样调控结合能？
最常见的写法把总能量拆成两部分：Etot=EKS+Edisp。EKS 来自所选泛函和电子密度自洽结果，Edisp 用额外模型补上长程色散。DFT-D2、DFT-D3、DFT-D4 属于这一路线。它们按原子对求和，再乘上阻尼函数和泛函相关的缩放因子。DFT-D 的计算位置很明确：色散项直接进入总能、力和应力，结构优化、吸附高度和晶胞参数都会读取这项贡献；应力项还会进入晶胞优化，层间距和分子晶体晶胞会随之响应。
阻尼函数控制短程区域怎样接上原泛函。没有阻尼，-C6/R6 在小距离处会给出过强吸引；阻尼太强，又会把应有的平衡距离吸引削弱。rVV10、VV10 和 vdW-DF 采用另一种写法，把电子密度两点之间的相关写入非局域项，并把非局域相关能记作 Ecnl，通过电子密度在空间中两点之间的核函数描述相关涨落。参数调节看似只改一个数，实际是在移动吸引尾部与短程阻尼的交接位置。
图2. r2SCAN+rVV10 中 b 参数会改变 Ar2 结合能误差和结合曲线。DOI：10.1103/PhysRevB.106.075422
图 2 把 b 参数和 Ar2 结合能放在一起。b 较小时，阻尼较弱，修正项较早进入短程区域；b 较大时，短程吸引被压低。曲线最低点的位置和深度共同决定平衡距离与结合能。势能曲线的斜率就是力，曲线形状改变后，离子步优化会沿着新的力场移动。这也是 vdW 修正影响结构优化，而不只影响最后单点能的原因。
图3. 苯二聚体 T 构型中，r2SCAN+rVV10 与 SCAN+rVV10 给出的结合能曲线和力曲线对参数、网格和参考数据表现出不同偏差。DOI：10.1103/PhysRevB.106.075422
苯二聚体比惰性气体多了分子形状、π 电子云和取向差异。图 3 中实线是结合能，虚线是力。低能构型附近，不同参数给出的能量最低点接近，但斥力壁和长程尾部仍有差别。对于表面吸附、孔道限域和分子晶体，分子姿态往往由几个很接近的构型竞争决定，小到几 kJ mol-1 的色散能差就足以改变构型排序。
各类色散修正该如何选用？
DFT-D 系列把色散写成原子对修正。D2 使用较固定的原子参数；D3 引入配位数，让 C6 系数随局域环境变化；D4 进一步把原子电荷信息纳入参数。Tkatchenko-Scheffler 方法从电子密度和有效体积估计极化率，MBD 把原子响应写成耦合振子，处理多体色散。vdW-DF、VV10、rVV10 直接在泛函中加入非局域相关能。这些方法的差别在于色散来源怎样建模、环境响应怎样进入、短程区域怎样阻尼，对应的参数、核函数和计算成本随之变化。
原子对模型计算量小，便于在大超胞吸附和分子晶体中使用；多体色散适合极化率较大、片段较多、屏蔽效应明显的体系；非局域 vdW 泛函把电子密度写进相关核，对分子、固体和表面之间的统一描述更自然。单一方法很难同时压低分子复合物、金属表面和层状晶体的全部误差。方法选择应跟研究对象绑定：弱分子复合物、层状晶体、金属表面吸附和柔性有机晶体对应的误差来源并不相同。
图4. S22、S66x8 和 L7 非共价作用集合中的平均绝对偏差显示，色散修正对不同尺寸复合物的误差影响不同。DOI：10.1103/PhysRevB.94.115144
图 4 把小二聚体 S22、距离扫描集合 S66x8 和更大的 L7 复合物放在同一张误差图里。SCAN-D3 在 S22 和 S66x8 上误差较低，到了 L7，大尺寸接触面、构象自由度和极化响应会放大方法差别。过强修正会把平衡距离压短或把结合能拉深。某个方法在小分子复合物上表现好，转到金属表面、孔道限域或强极化层状材料时仍要检查基准体系、结构误差和能量排序。
实际写计算方案时，方法名后面应标出母泛函。PBE-D3、PBE-D4、SCAN+rVV10、r2SCAN+rVV10、optB88-vdW 给出的修正对象和基础交换相关项不同。母泛函本身会影响键长、电荷密度、磁态和能带，色散项再改变弱相互作用部分。相同的 vdW 标签配上不同母泛函，得到的结构和吸附能不应混作同一套数据。两组数据合并前，应重新核对几何和能量基准。
vdW 修正会影响哪些电子性质？
层状材料的层间距、分子在表面的吸附高度、孔道内客体分子的姿态、有机晶体的堆积方式，都由短程排斥和长程吸引共同限定。弱相互作用主导的距离和相对能量，对 vdW 修正最敏感；强共价骨架内部的键长、面内晶格常数和局域配位通常变化较小，但也可能通过结构耦合出现间接变化。结构响应通常沿弱相互作用方向先出现，再影响派生图谱。
图5. L28 层状材料中层间结合能 Eb 的误差箱线图显示，不同 vdW 泛函给出的层间吸引强度偏差不同。DOI：10.1103/PhysRevB.106.075422
图 5 单独看 L28 层状材料的层间结合能 Eb。rev-vdW-DF2、SCAN+rVV10 和 r2SCAN+rVV10 的箱线位置和离散程度不同，说明层间吸引强度并不会因写上 vdW 标签就自动一致。层内通常是共价或离子-共价混合键，层间主要是色散、静电和少量轨道重叠。同一结构里不同方向的物理来源不同，vdW 修正对各向异性结构的影响会分方向出现。
层间距改变后，层间剪切模、呼吸模和低频声子分支会移动；吸附高度改变后，分子轨道与表面态的重叠、功函数和吸附态 DOS 也可能变化。电子结构位移若来自几何变化，应把结构效应和泛函本身的电子效应分开报告。同一材料可用固定几何单点计算和重新优化结构计算作对照，区分“色散项改变能量面”和“新结构改变电子态”。
图6. 石墨和 MoS2 的声子色散显示，加入 rVV10 后，沿层间方向的低频分支与实验数据的对应关系发生变化。DOI：10.1103/PhysRevB.106.075422
在石墨和 MoS2 这类层状体系里，低频分支对层间力常数敏感。图 6 中沿层间方向的 Γ-A 区域，加入 rVV10 后曲线更接近实验点。力常数矩阵来自能量对原子位移的二阶导数，弱相互作用势能面一旦改变，低频振动会给出响应。对于热导率、相变和层间滑移相关研究，声子图比单个晶格常数更能暴露这种影响。
图7. PPTA 有机晶体结构展示了柔性分子链和层状堆积中短程成键与弱相互作用并存的场景。DOI：10.1103/PhysRevB.106.075422
有机晶体、分子筛孔道和软界面还多了一层复杂度：分子内部共价骨架较硬，分子间堆积较软。PPTA 一类体系中，链间距离、π-π 堆积和氢键网络共同决定晶胞。若研究对象含有大分子、芳香环、层状堆积或疏水片段，未说明 vdW 处理方式会削弱结构和能量数据的可比性。报告中至少应写明母泛函、vdW 方法、阻尼形式或参数版本，以及结构是否在该方法下重新优化。
范德华修正适合回答弱相互作用能量面怎样变化。U 值对局域 d/f 电子的处理、PBE 带隙偏小、反应过渡态搜索不足和溶剂模型缺失，各自对应不同误差来源，不宜用 vdW 修正代替。吸附能、界面距离、分子取向、层间声子和柔性晶体堆积属于它最常改变的读数；价态、磁矩、带边位置和反应能垒则要看结构变化、电子相关、环境模型和反应路径设置各自贡献了多少。