过渡金属氧化物表面吸附能计算中的范德华修正——如何判断你的体系是否需要加VDW

一个做催化计算的同行把INCAR发过来，说他用标准PBE算出来的甲醇在CeO₂(111)面的吸附能只有0.15eV，但实验文献报道的TPD脱附温度换算下来吸附能应该在0.6-0.8eV之间。差了四五倍，问我们有没有遇到过类似的情况。

遇到过。不止一次。而且问题不在INCAR设置上——他所有的收敛标准、k点密度、截断能都挑不出毛病。问题出在他没加范德华修正。

什么时候必须加VDW修正

对于过渡金属氧化物表面的分子吸附，这个问题的判断依据其实很明确。标准GGA泛函（PBE、PW91）对色散力的描述几乎是零——它不是”算不准”，是”根本没算”。电子密度梯度只能捕捉到电荷密度重叠区间的交换-关联效应，而伦敦色散力来自瞬时偶极-偶极相互作用，跟电子密度梯度没直接关系。

所以判断要不要加VDW修正，不看体系大小，看吸附质和基底之间的物理吸附成分有多重。两个信号可以帮你判断：

第一个，吸附质有没有闭合电子壳层。甲醇是闭壳层分子，没有未配对电子与表面形成强共价键，主要的吸附驱动力是氢键和范德华力。闭壳层分子在氧化物表面的吸附，不加VDW修正基本就是underbinding。

第二个，标准PBE算出来的吸附能和实验值的差距是不是系统性的。如果差了两三倍甚至更多，不是0.1-0.2eV的量级偏差，那大概率不是电子结构方法的精度问题，是少算了物理吸附分量。

这个项目里，吸附质甲醇是闭壳层，表面是CeO₂(111)（离子型氧化物，表面极化也会贡献额外的静电-色散耦合），两个信号都到了——必须加VDW。

三种方案的实际表现

加VDW容易，选哪个加就没那么简单了。这个项目里我们并行跑了三种方案，用同一套计算参数（3×3×3 CeO₂超胞，截断能500eV，k点2×2×1，能量收敛1×10⁻⁶ eV），只换VDW方法，对比每个方案给出来的吸附能。

DFT-D3（零阻尼）：Grimme的经典方案，在总能上加一个C₆·R⁻⁶的阻尼函数对。跑出来吸附能0.72eV，跟实验值吻合得很好。但有一个隐患——零阻尼在短程会过度修正，R太小的原子对会被C₆/R⁶项拉出非物理的吸引力。CeO₂体系里Ce-O键本来就短（约2.34Å），如果吸附质靠近表面太近，零阻尼的短程过修正会出问题。

DFT-D3(BJ)：用Becke-Johnson阻尼替代零阻尼，在短程做衰减。吸附能0.61eV，比零阻尼低了约0.1eV。这0.1eV的差距来自BJ阻尼对近距离Ce-O对的色散项做了衰减——甲醇的氧原子距离表面Ce原子最近处约2.8Å，这个距离上BJ阻尼已经开始起作用了。

vdW-DF2：非局域vdW泛函，色散能从电子密度直接算出来，不是后加的修正项。吸附能0.58eV。比DFT-D3(BJ)再低一点，但趋势一致。vdW-DF2的优势是物理图像更干净——色散能是自洽计算的一部分，不存在”加多少”的选择问题。代价是计算量比DFT-D3大了约40%。

三组数据摆在一起：DFT-D3=0.72eV，DFT-D3(BJ)=0.61eV，vdW-DF2=0.58eV。实验参考值在0.6-0.8eV之间。DFT-D3(BJ)和vdW-DF2都在合理的误差带内，DFT-D3虽然最接近实验值的上界，但那0.72eV里可能混入了短程过修正的贡献。

最终选择和检查

最终选的是DFT-D3(BJ)。理由有三个：一，比vdW-DF2便宜40%的计算量，后续还要跑过渡态搜索和频率分析，计算量不能翻倍；二，BJ阻尼在物理上更合理，短程衰减避免了零阻尼的过修正风险；三，0.61eV落在实验范围的偏中位置，没有系统性地偏高或偏低。

选定方案之后还做了一轮检查——把甲醇替换成水分子（极性更强、分子更小），看DFT-D3(BJ)给的水的吸附能是否也合理。水在CeO₂(111)面的吸附能用DFT-D3(BJ)算出来是0.52eV，跟文献报道的0.45-0.60eV范围一致。这意味着方案不仅对甲醇成立，对同类小分子吸附也自洽。

这个项目的结论可以推广：过渡金属氧化物表面的闭壳层分子吸附，不加VDW修正的误差不是百分之几，是成倍的。选DFT-D3(BJ)作为默认方案是一个性价比和物理合理性的平衡点，但如果吸附质含有重原子（碘、金属有机配体等），C₆系数的精度可能不够，VDW-DF2会是更安全的选择。

每一个涉及表面吸附的项目，第一件事不是拿起INCAR就跑，而是问自己：这个吸附，物理吸附占多少？如果占得不少，VDW修正不是可选项，是必选项。