VASP吸附能计算：从表面建模到能量收敛的完整经验总结

VASP吸附能计算在催化、表面科学和新能源材料研究中出现的频率相当高，但真正跑出可靠结果的，往往需要在建模阶段就把坑踩清楚。

吸附能的定义本身不复杂：E_ads = E(表面+吸附物) – E(表面) – E(吸附物)。问题出在三个单独计算之间的一致性——同一套截断能、同一个自旋设置、同一种泛函，缺一不可。不少计算结果出现几十到上百meV的偏差，往往不是因为物理模型选错了，而是三个计算的参数没有对齐。

表面模型：层数、真空层与固定方式

表面模型的构建直接决定计算的物理合理性。常规做法是取4-6层slab，底部2-3层固定模拟体相约束，顶部层充分弛豫。真空层厚度通常设置在15Å以上，目的是消除周期镜像之间的虚假相互作用。

这里有一个容易被忽略的点：真空层厚度的收敛测试很少有人做，但它对吸附能的影响并不小。特别是极性表面或者有偶极修正需求的体系，真空层不够会让能量出现系统性偏移。建议在正式计算前先跑一组真空层从12Å到20Å的对比，确认结果趋于收敛后再固定下来。

k点设置与ENCUT的取舍

VASP吸附能计算中，k点密度往往是精度与效率矛盾最突出的地方。表面体系的布里渊区在z方向已经被压缩，通常只取1个k点就够;但xy平面方向的k点仍需要根据超胞大小认真收敛。

ENCUT的选取建议参考POTCAR中各元素ENMAX值的1.3倍，宁可偏高不要偏低。对于含有第一行过渡金属的体系，ENCUT不足带来的误差会比真空层不够更难排查，因为它会均匀地影响所有能量项，看起来结果很”稳定”，但绝对值已经偏离。

自旋极化与DFT+U的处理

含3d过渡金属的表面体系，ISPIN=2是标配。如果涉及Fe、Co、Ni等强关联体系，还需要考虑DFT+U修正。U值的选取是个老问题，文献中同一元素的U值范围往往跨度很大，建议优先参照与研究体系最接近的实验数据来校准，而不是直接套用文献中随机挑选的U值。

还有一个细节：孤立吸附物分子的参考能量计算，也必须开自旋。O2分子的三重态是典型案例——如果用单重态去算参考能量，吸附能会出现0.5-1 eV量级的误差，而这个误差不会在任何收敛测试中暴露出来。

VASP吸附能计算的结果解读与误差来源总结

拿到吸附能数值之后，建议首先检查吸附构型是否合理——键长、键角与文献比对，异常值往往指向局部弛豫陷入了局域极值。其次检查磁矩：吸附前后表面磁矩变化过大，通常意味着自旋态在吸附过程中发生了翻转，此时直接比较能量会引入额外误差。

VASP吸附能计算的精度上限，很大程度上取决于泛函的选择和参数的一致性。PBE在大多数金属表面体系上的表现是可接受的，但对于弱相互作用主导的物理吸附，van der Waals修正(如DFT-D3)不能省。把这几个环节都对齐之后，VASP吸附能计算的结果才值得拿去和实验数据比较。