



分子吸附的 Eads 同时受几何高度、表面极化、分子取向和色散相互作用影响。vdW 修正改变的不只是吸附能大小,还可能改变吸附高度、位点排序和后续反应端点。对芳香分子、卤代分子、层状材料和多离子储能体系,D3、D4 或非局域 vdW 泛函的选择本身就属于计算口径。




分子吸附缺了哪类相互作用?
常规 GGA 泛函对长程色散描述有限。苯、卤代苯、CO₂、N₂ 和多硫化锂这类物种靠近表面时,除了局域成键和电荷转移,还存在由瞬时偶极相关产生的吸引作用。若忽略这部分能量,吸附高度可能偏大,弱吸附物的 Eads 也可能偏弱。
以苯在 Cu(111)、Ag(111) 和 Au(111) 上的高度与吸附能对照为例,同一个吸附物在不同金属上的色散贡献并不相同,因为表面电子响应、吸附距离和分子极化率都会变化。若只报告 PBE 吸附能而不说明色散处理,芳香分子和大分子体系的可比性会明显下降。

图1 苯在 Cu(111)、Ag(111) 和 Au(111) 上的平衡高度与吸附能随泛函变化。该图用于说明色散修正会同时改变几何和能量。 DOI: 10.1002/cphc.202400865。
D3、D4 和 vdW-DF 有何差别?
D3 与 D4 属于经验或半经验色散修正,通常作为总能的附加项加入;D4 会进一步考虑电荷相关极化率。optB86b-vdW、vdW-DF、rev-vdW-DF2 等非局域泛函则在交换-相关形式中处理非局域相关。它们的计算代价、适用体系和吸附高度表现并不完全一致。
D3 的优势是设置简单、计算代价小,适合大批量吸附构型筛选;D4 对电荷环境更敏感,适合离子性或电荷转移较明显的体系作为对照;非局域 vdW 泛函在层状材料、二维界面和芳香分子吸附中更常用于检查平衡距离。三类方法的选择要对应目标量:吸附能、层间距和反应端点对色散项的敏感位置并不相同。
在 VASP 中,D3/D4 更适合快速筛选和与既有文献口径对齐;非局域 vdW 泛函更适合对吸附高度和弱相互作用更敏感的体系。选哪种方案应由研究对象和对照口径决定。若旧数据全部采用 PBE-D3,新数据突然改用 optB86b-vdW,吸附能差异不能直接解释为材料本征差异。




吸附高度为什么会变?
色散吸引会让分子更靠近表面,但靠近后又会增强 Pauli 排斥和轨道重叠。最终吸附高度来自两类作用的平衡。卤代苯体系中,Br 或 I 原子极化率较大,分子平面、卤素端高度和金属表面距离都可能对所选泛函敏感。
图2 显示溴苯和碘苯的高度分量随泛函变化。吸附高度变化会继续传递到键长、PDOS 和反应端点。若优化结构来自一种 vdW 方案,而静态能量来自另一种方案,几何和能量口径就被拆开了。更合适的做法是先确定主计算方案,再对代表性结构做交叉对照。

图2 溴苯和碘苯在贵金属表面的吸附高度对比。该图说明重原子取代会放大色散与吸附构型耦合。 DOI: 10.1002/cphc.202400865。
位点排序会不会改变?
对弱吸附体系,不同位点的能量差可能很小,色散项的变化足以改变排序。对化学吸附体系,局域成键可能占主导,但色散仍会通过吸附高度和分子取向影响能量。图3 中卤代苯吸附能随泛函变化,说明同一分子在不同金属表面上可能出现不同敏感性。
若顶位、桥位和空位构型的 Eads 排序随 vdW 方案改变,应先检查构型是否都充分弛豫、表面层数和真空层是否一致、k 点和 SIGMA 是否收敛。位点排序只有在同一 vdW 方案下才有直接比较意义;跨方案结果更适合用于敏感性评估。

图3 溴苯和碘苯在 Cu、Ag、Au 表面的吸附能随泛函变化。该图用于说明不同 vdW 方案可改变吸附能幅度。 DOI: 10.1002/cphc.202400865。




端点和过渡态要统一吗?
反应能垒来自初态、过渡态和终态之间的能量差。若初态和终态采用 PBE-D3 优化,而 NEB 或过渡态单点采用无 vdW 的 PBE,势能剖面会混入方法差异。对于分子脱氢、脱卤、CO₂ 活化或多硫化物分解,这个影响可能直接改变限速步判断。
图4 的脱氢势能剖面把构型和能量路径放在一起。反应路径中的每个端点都应使用同一色散口径。若为了节省计算只在终态做 vdW 单点,应把结果定义为近似对照,而不宜与全程 vdW 优化路径混作同一结论。

图4 苯在 Cu、Ag、Au 表面的脱氢反应势能剖面。该图说明吸附构型和色散处理会进入反应路径能量。 DOI: 10.1002/cphc.202400865。
自由能校正要不要重算?
vdW 修正改变结构后,振动频率也可能改变,尤其是低频平动、转动受限模式。对强吸附小分子,ZPE 差异可能不大;对大分子、层状吸附或柔性物种,低频模式对熵项更敏感。若只比较电子吸附能,可以说明未包含振动贡献;若讨论 ΔG,则要定义 ZPE 与熵校正的来源。
当吸附高度变化较大时,低频振动模式会最先受到影响。苯环平行表面的滑移、卤代苯的倾斜和多硫化物链段的摆动,都可能在不同 vdW 方案下形成不同频率分布。若目标是自由能台阶而不是单点吸附能,代表性端点的频率校正就不能完全省略。
先用主 vdW 方案得到优化结构和静态能量,再对关键构型计算频率。若两种 vdW 方案给出相反排序,需要比较吸附高度、关键键长、差分电荷密度和 PDOS,而不是只比较总能。vdW 敏感性本身也是结果边界,说明该吸附体系受弱相互作用控制较强。




分子和层状材料如何选?
芳香分子、有机配体、CO₂RR 多碳中间体、Li-S 多硫化物和 MXene/石墨烯界面常需要色散修正。若目标是与多数催化文献比较,PBE-D3 或 PBE-D4 常见且便于复现;若目标是吸附高度、层间距或界面结合能,optB86b-vdW 等非局域方案可作为主方案或关键对照。
对金属表面小分子,色散项未必支配成键,但仍可能改变弱吸附端点的初始高度。对氧化物或半导体表面,极性位点和氢键会与 vdW 共同作用。选择方案时,应把吸附物大小、表面极化程度和是否存在层间相互作用一起考虑。
图5 的 V₂N 储离子算例直接比较 DFT-D3、vdW-DF 和 optPBE-vdW 的吸附能随覆盖度变化。它说明覆盖度、离子种类和 vdW 方案会同时影响吸附能曲线。当 Li、Na、K、Mg 或 Zn 等离子在二维表面逐步覆盖时,方案差异不只来自单个离子吸附,还会来自离子间排斥和表面弛豫。

图5 V₂N 单层上不同离子覆盖度吸附能及 DFT-D3、vdW-DF、optPBE-vdW 对照。该图用于说明储离子体系也要统一色散口径。 DOI: 10.1021/acsomega.2c00936。
VASP 输入怎样保持清楚?
D3 常用 IVDW 参数调用,D4 需确认当前 VASP 版本和编译支持;非局域 vdW 泛函通常要设置对应 GGA/AGGAC/LUSE_VDW 或相关参数,具体以本机 VASP 版本文档为准。所有参数应写成示例,并注明需按体系测试 ENCUT、k 点和真空层。
若研究对象是 CO₂、N₂、苯环或多硫化锂,可用同一套初构型分别完成无 vdW 与主 vdW 方案优化,再比较吸附高度、关键键长和 Eads。如果无 vdW 结构明显远离表面,而加 vdW 后出现合理接触距离,说明色散项对几何端点有决定性影响。
若目标是比较不同材料,不宜让每个材料使用不同 vdW 方案。PBE-D3 可以作为主结果口径,少数代表性材料再用 optB86b-vdW 对照;大批量筛选保持同一色散处理后,关键排序是否依赖修正方法也能单独判断。
vdW 修正的核心不是追求某个固定数值,而是保持同一组比较的物理口径一致。同一篇计算中,清洁表面、吸附态、孤立分子、NEB 端点和频率校正应尽量采用同一方案;确需跨方案比较时,应把它定义为敏感性对照,并说明哪些排序保持一致、哪些排序依赖修正方法。




要点:vdW 修正会同时影响吸附能、吸附高度、位点排序和反应端点。
要点:D3/D4 便于筛选和文献对齐,非局域 vdW 泛函更适合检查高度和弱相互作用。
要点:反应路径中的初态、过渡态和终态应尽量使用同一色散口径。
要点:跨 vdW 方案的结果应作为敏感性对照,而不是直接当作材料差异。
