吸附能是表面催化、电化学界面和气体传感中最基础的物理量之一，VASP是吸附能计算的主流工具，新手常因参数设置不规范、模型构建不合理，得到与文献偏差较大的结果。这篇文章将为大家整理VASP计算吸附能的标准操作与关键注意事项。

 

一、吸附能的定义与计算原理

吸附能描述吸附质分子与基底表面结合前后的能量变化，反映吸附作用的强弱与自发性。VASP 中采用基态电子总能进行计算，标准公式为：

E_ads = E_{slab +mol} −E_slab−E_mol

其中E_{slab +mol}为吸附后复合体系的总能量，E_slab为清洁基底表面的总能量，E_mol为气相孤立吸附质分子的总能量，三项均取自VASP优化完成后的自洽总能（OUTCAR中energy(sigma->0)对应数值）。

计算结果为负值表示吸附放热，数值越负吸附作用越强；通常大于-0.5 eV为物理吸附（范德华作用为主），小于-1eV为化学吸附（形成化学键）。

 

二、VASP 吸附能计算标准流程

完整计算需完成清洁表面、孤立分子、吸附复合体系三类体系的运算，全程保持计算参数完全统一，具体分为三步：

1.模型构建：保证合理性是前提

截取目标晶面构建周期性 slab 模型，密排金属体系通常堆叠 3~5 层原子，氧化物、二维材料需增至 5~7 层，确保表面原子弛豫充分；固定底部 1~2 层原子模拟体相晶格约束，避免优化时整体结构形变。垂直表面方向设置 15~20Å 真空层，若吸附质为大尺寸有机分子，需增至 20Å 以上，消除周期镜像的相互作用干扰。

孤立吸附质分子置于 15Å×15Å×15Å 的立方体真空盒中，避免周期性边界带来的分子间作用。吸附初始构型需覆盖表面高对称位点（顶位、桥位、fcc 穴位、hcp 穴位），同时调整分子吸附角度，避免遗漏全局最优构型。

 

2.参数配置：统一标准是核心

所有体系必须使用同一套赝势与泛函，INCAR 关键参数统一设置：截断能ENCUT需高于赝势推荐值的 1.2~1.5 倍，过渡金属体系通常设为 400~500eV，确保基组收敛；电子自洽阈值EDIFF=1E-5，离子优化采用共轭梯度法（IBRION=2），力收敛阈值EDIFFG=-0.02 eV/Å；统一开启 DFT-D3 (BJ) 色散校正（IVDW=12），保证弱相互作用计算精度。

K 点网格保持一致：表面体系根据超胞大小设置 2×2×1 或 3×3×1 的 Gamma 中心网格，孤立分子统一使用 1×1×1 Gamma 点，兼顾精度与计算效率。

 

3.分步计算与能量提取

建议先对多个吸附构型做粗优化快速筛选，再对能量最低的优势构型做高精度全弛豫精算，节省算力。所有体系收敛后，从 OUTCAR 文件中提取energy(sigma->0)对应的零展宽外推总能，而非直接读取末步能量，降低数值误差。磁性体系需额外确认磁矩收敛稳定，避免自旋态跳变导致结果失真。

 

三、新手注意避坑

1.参数一致性是底线

三个体系的泛函、赝势、截断能、k 点密度必须完全一致，任何一项改动都会导致能量差失去物理意义。

2.色散修正不能省略

对于物理吸附、有机分子体系，常规 GGA 泛函会严重低估吸附强度，需添加 DFT-D3 (BJ) 校正（VASP 中设置IVDW=12）。

3.警惕假收敛

不能仅以离子步数判断收敛，需确认受力与能量变化均满足阈值；过渡金属体系建议适当增加电子迭代步数。

4.偶极修正补误差

对于极性表面、不对称吸附体系，需开启LDIPOL=.TRUE.修正表面偶极矩，消除周期性带来的能量误差。

 

总结来说，VASP 吸附能计算是入门级任务却处处是细节，规范的模型与参数设置是结果可靠的前提。掌握以上要点，就能快速得到可重复、可对比的计算结果。

 

【科学计算 优选至简】 

至简科技提供专业的第一性原理、分子动力学、量子化学、机器学习、有限元仿真、相图计算等代算服务。 

✅50➕全职工程师 500➕合作高校企业
✅部分合作项目己发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊