VASP吸附能计算的完整流程与常见误差来源

吸附能公式背后的三个独立计算

VASP吸附能计算的标准公式是：

E_ads = E(surface+adsorbate) – E(surface) – E(adsorbate)

看似简单，但这三个能量值来自三个独立的VASP计算，每一个的参数设置都会影响最终结果。本项目中见过不少研究者三个计算用了不同的ENCUT或K点网格，导致吸附能误差超过0.5 eV——这在判断化学吸附还是物理吸附时是致命的。

核心原则：三个计算必须使用完全相同的ENCUT、K点密度（面内）、赝势、SIGMA和ISMEAR。唯一不同的是K点在真空方向的取值——孤立分子用Gamma点，表面用1×1×1（即仍然只取Gamma点，但slab模型的面内K点要与表面一致）。

Slab模型构建：厚度和真空层

吸附能计算的精度高度依赖slab模型的质量。以本项目近期完成的Pt(111)-O吸附体系为例：

Slab厚度的选择需要做收敛测试。本项目的方法是：固定其他参数，逐层增加slab厚度（3→4→5→6层），观察表面能的收敛情况。Pt(111)的测试结果显示，4层时表面能为1.32 J/m²，5层稳定在1.48 J/m²，6层为1.49 J/m²——说明5层已基本收敛。

超胞大小直接影响覆盖度。p(2×2)对应1/4 ML覆盖度，p(3×3)对应1/9 ML。覆盖度越高，吸附分子之间的排斥越强，吸附能越偏正（越弱）。如果要模拟低覆盖度极限，至少需要p(3×3)或更大的超胞。在[VASP/DFT计算服务](https://www.keyanxueshu.com/category/dft/)中，本项目通常推荐p(3×3)作为默认超胞。

吸附位点的系统搜索

对于fcc(111)表面，常见的吸附位点有四种：

• ontop：吸附质在表面原子的正上方
• bridge：两个表面原子中间
• fcc-hollow：三个表面原子围成的空位（无下层原子）
• hcp-hollow：三个表面原子围成的空位（正下方有下层原子）

本项目在每个吸附体系中都会做完整的位点搜索。具体做法：先用较粗的参数（ENCUT=400, EDIFFG=-0.05）对四个位点各做一轮快速优化，找到能量最低的位点后，再用精细参数（ENCUT=500, EDIFFG=-0.02）做最终优化。

一个容易被忽视的点：初始吸附高度对结果有影响。本项目取的经验值是：

孤立分子计算的技巧

孤立分子的能量计算看似最简单，实际最容易出错。关键点：

盒子大小：孤立分子放在足够大的立方盒子中，本项目通常用15Å×15Å×15Å。太小会有周期镜像相互作用，太大则浪费计算资源。可以通过比较10Å和15Å盒子的能量差来判断——如果差值超过1 meV，说明盒子不够大。

K点：Gamma点即可，无需更密。

注意分子构型：孤立分子的构型应该是自由优化后的，不能直接从实验晶体结构中提取。比如CO分子的C-O键长，气相中约1.128Å，如果直接用晶体中的1.15Å会导致分子的零点能偏高。

基组重叠误差（BSSE）校正

当吸附质与表面的距离较近时（化学吸附通常1-3Å），吸附质的波函数与表面原子的波函数会重叠，导致能量被人为降低。这就是基组重叠误差（BSSE），也叫Counterpoise Correction。

VASP中做BSSE校正的方法是引入鬼原子（ghost atom）：

1. 计算E(surface+adsorbate) —— 正常的吸附体系
2. 计算E(surface) —— 去掉吸附质，保留表面
3. 计算E(adsorbate) —— 去掉表面，保留吸附质
4. 计算E(adsorbate in ghost surface) —— 吸附质在表面位置但表面原子用鬼原子（无电子）
5. BSSE = E(adsorbate in ghost surface) – E(adsorbate)

校正后吸附能 = E_ads + BSSE（BSSE为正值，使吸附能变正，即吸附变弱）

对于弱吸附体系（物理吸附，vdW主导），BSSE可达0.2-0.5 eV，必须校正。对于强化学吸附（>2 eV），BSSE通常<0.05 eV，可以忽略。更多关于[分子模拟计算](https://www.keyanxueshu.com/)的误差分析方法，本项目在相关文章中有展开讨论。

vdW色散校正

标准DFT（PBE/LDA）无法正确描述范德华相互作用，这对物理吸附体系是致命的。VASP中常用的vdW校正方法：

本项目在Pt(111)-CO体系中对比了PBE和DFT-D3的结果：

• PBE：E_ads = -1.85 eV（过强）
• DFT-D3：E_ads = -1.52 eV
• 实验值：-1.43 ± 0.10 eV

DFT-D3将误差从0.42 eV降到0.09 eV，改善显著。

误差来源汇总

验证方法

计算完成后，本项目会用以下方式交叉验证：

6. 与实验值对比：查阅文献中TPD（程序升温脱附）或量热法给出的实验吸附能
7. 电荷转移分析：用Bader电荷分析计算吸附前后的电荷转移量，化学吸附通常有1-1.0 e的转移
8. 差分电荷密度：绘制ρ(complex) – ρ(surface) – ρ(adsorbate)，直观观察电荷重分布
9. 态密度分析：比较吸附前后表面d带中心和吸附质分子轨道的变化

以Pt(111)-O体系为例，Bader电荷分析显示O原子获得0.72 e，Pt失去等量电荷，与O的强电负性一致。差分电荷密度图在Pt-O键区域出现明显的电荷积聚，证实了化学吸附本质。这些分析在[科研学术网](https://www.keyanxueshu.com/)的其他文章中有详细案例展示。

VASP吸附能计算不是跑通三个SCF就完事的简单流程，从模型构建到误差校正，每一步都可能引入0.1 eV量级的偏差。系统性地控制这些误差源，才能得到可信的吸附能数值。