范德华力理论计算 — vdW-DF泛函家族从optB86b到SCAN+rVV10的演进之路

为什么”标准”DFT算不准色散力

故事要从2015年说起。当时我们课题组在尝试计算一个有机半导体晶体（并五苯）的晶格参数。用PBE泛函优化得到的c轴晶格常数比实验值长了0.8 Å——这在晶体学上是个天文数字。问题出在哪里？PBE作为GGA泛函，对电子密度的非局域关联效应——也就是范德华力——几乎完全忽视。

在接下来的十年里，范德华力理论计算经历了一场静悄悄的革命。从半经验的DFT-D修正，到基于电子密度的非局域vdW泛函，再到meta-GGA级别的全自洽vdW处理，这个领域的发展速度令人惊叹。但与此同时，选择困难症也来了：面对optB86b-vdW、optPBE-vdW、rev-vdW-DF2、DFT-D3(BJ)、SCAN+rVV10等众多选项，哪一个是你的体系应该用的？

本文不是泛函综述，而是一次实战复盘的记录。我们将用一个并五苯晶体的计算案例作为主线，展示从”PBE算不准”到”找到合适的vdW方法”的全过程，以及在这个过程中对各种vdW-DF泛函的理解和误判。

第一代vdW-DF：Dion等开创的非局域关联方法

2004年Dion等人提出的vdW-DF是范德华力理论计算的里程碑。这个方法的精髓在于：它不像DFT-D那样在总能量上附加一个经验的C₆/R⁶项，而是直接在交换关联泛函的关联部分引入非局域贡献。具体来说，vdW-DF将关联能分解为局域部分（LDA关联）和非局域部分：E_c[n] = E_c^LDA[n] + E_c^nl[n]。其中E_c^nl通过一个核函数Φ(q₀, q₀’, |r-r’|)的双重积分来计算，依赖于局域等离激元频率q₀(r)。

这个方法在物理上很优美——色散力自然地产生于电子密度涨落之间的关联效应。但在实践中，原始的vdW-DF对结合能普遍高估，且交换部分选择revPBE也带来了一些问题。我们测试并五苯时，发现vdW-DF给出的层间结合能是-2.1 eV，而实验参考值约为-1.5 eV，高估了约40%。

不过，vdW-DF的框架是正确的，后续的改进正是围绕着交换泛函的选择和核函数的形式展开的。理解这个基础框架，是理解所有vdW-DF变体的前提。

optB86b-vdW：2010年代的”黄金标准”

optB86b-vdW由Klimeš等人在2010-2011年提出，核心改进是将vdW-DF中的revPBE交换泛函替换为优化过的optB86b。这个看似简单的更换产生了巨大影响：optB86b的交换增强因子在低密度区域的衰减行为更适合色散主导的体系，同时保持了GGA对共价键和离子键的合理描述。

在我们的并五苯测试中，optB86b-vdW将层间结合能修正到了-1.62 eV，与实验值的偏差缩小到了约8%，层间距也从PBE的3.8 Å改善到了3.5 Å（实验值3.4 Å）。这个水平对于晶体结构预测来说已经相当实用。

optB86b-vdW之所以成为”黄金标准”，还因为它对多种类型的相互作用（色散、氢键、静电）都能给出均衡的表现，不像某些方法专门优化色散力却在氢键上表现糟糕。当然，它不是完美的——在涉及过渡金属氧化物的表面吸附体系中，optB86b-vdW有低估吸附能的趋势，因为optB86b交换泛函对d电子的自相互作用误差修正不足。

rev-vdW-DF2和DFT-D3(BJ)：不同的哲学路线

rev-vdW-DF2代表了vdW-DF家族的另一条优化路线：通过改进核函数本身来提升精度，而不是交换泛函。Hamada在2014年重新参数化了vdW-DF2的核函数，使其在中等分离距离（3-5 Å）处的色散力描述更精确。这对有机分子晶体的适用性很好，但在层状材料（如石墨、MoS₂）中表现不如optB86b。

相比之下，DFT-D3(BJ)走了一条完全不同的路——纯经验的色散修正。Grimme的D3方法不涉及任何非局域关联的计算，而是在每个时间步后根据原子对的几何构型加上C₆/R⁶ + C₈/R⁸的色散能量项，其中系数来自预计算的原子极化率数据库。Becke-Johnson阻尼函数用于在短程区域平滑地关闭色散项。

DFT-D3(BJ)的优势是计算效率极高（几乎零额外开销）和对几乎所有泛函的即插即用兼容性。但它的根本局限在于：色散系数是原子对的函数，不依赖于化学环境——一个羰基碳和一个甲基碳的C₆系数是相同的，这显然不物理。

SCAN+rVV10：meta-GGA时代的vdW处理

SCAN+rVV10可能是目前范德华力理论计算中最令人兴奋的组合。SCAN泛函本身是2015年提出的meta-GGA泛函，满足所有17个已知的精确约束条件，在中程vdW相互作用上有比传统GGA更好的描述。但SCAN仍然缺少长程色散力，于是Peng等人将rVV10非局域关联核与SCAN交换泛函结合，形成了SCAN+rVV10。

这个组合的强大之处在于：SCAN的meta-GGA本质使其能比GGA更好地描述中等范围的电子关联（3-5 Å），而rVV10覆盖了长程（>5 Å）的色散尾部。两块拼图组合起来，实现了从共价键到色散力的无缝过渡。

我们对并五苯的测试中，SCAN+rVV10给出的层间结合能是-1.58 eV，层间距3.42 Å，比optB86b-vdW略好但差距不大。但SCAN+rVV10真正的优势在于体系多样性——同一个泛函可以同时准确描述水分子团簇的氢键网络（这是很多vdW-DF变体的痛点）和石墨的层间滑移能量面。如果你在做多相体系或多组分界面的模拟，SCAN+rVV10值得认真考虑。唯一需要注意的是计算成本：meta-GGA的动能密度依赖使得自洽收敛比GGA慢约2-3倍。

怎么选：一份基于体系特征的快速决策指南

经过多个项目的验证，我们总结了一套根据体系特征选择vdW方法的快速决策流程。如果主要关心中性有机分子晶体（如药物多晶型筛选），optB86b-vdW或optPBE-vdW是性价比最高的选择。如果体系包含过渡金属表面或氧化物界面，考虑rev-vdW-DF2或SCAN+rVV10——meta-GGA级别的交换对d电子的描述更好。如果做的是高通量筛选或需要快速计算数百个构型，DFT-D3(BJ)的效率优势无可替代。如果追求最高精度且预算允许，SCAN+rVV10是目前的顶点选择，但要确保有足够的计算资源。

最后一点建议：不要轻信任何一篇论文声称”XX方法最好”。vdW方法的表现高度依赖体系，我们每次接新项目时的第一件事就是用2-3个候选方法跑一个小规模的基准测试（比如你体系中最关键的几个构型），用高精度参考数据（如CCSD(T)/CBS或实验晶格参数）做标定。这2-3天的额外工作往往能避免后面几个月的错误方向。

——————————————————————————————

📐 科研计算代算服务

我们提供专业的第一性原理计算、分子动力学模拟、有限元仿真分析及材料表征测试服务。
团队由博士学历的计算化学、计算材料学专家组成，覆盖VASP、LAMMPS、GROMACS、ANSYS、COMSOL等主流工具链。

  🔗 立即咨询计算方案