VASP 吸附能计算

一、项目背景与科学问题

吸附能是表面科学、多相催化与气体传感研究的核心热力学量，它定量描述了分子与固体表面的结合强度。吸附能的大小直接决定催化反应的决速步骤、吸附质的表面覆盖度以及传感器的响应灵敏度。对同一表面，吸附质在不同吸附位点（顶位、桥位、穴位）的吸附能差异，决定了吸附质在表面的优先落位与扩散路径，进而影响反应路径的选择。

第一性原理 DFT 计算可从原子层面精确预测吸附能。本项目以一氧化碳 CO 在 Pt(111) 表面的吸附为模型体系，利用 VASP 计算 CO 在顶位（atop）、桥位（bridge）、面心立方穴位（fcc）与六方密排穴位（hcp）四种高对称位点的吸附能，分析最稳定吸附构型与吸附位点偏好。

二、计算方法与理论背景

2.1 吸附能的定义

吸附能定义为 E_ads = E_{CO/slab} − E_{slab} − E_{CO}，其中 E_{CO/slab} 为吸附体系总能量，E_{slab} 为裸表面总能，E_{CO} 为气相 CO 分子能量。E_ads 为负值表示放热吸附，绝对值越大结合越强。该定义已隐式包含基组叠加误差（BSSE）与零点能的处理约定，实际计算中需保持三套体系使用相同的截断能与 k 点设置以保证误差抵消。

2.2 表面模型

Pt(111) 面为面心立方金属最密排表面，采用周期性平板模型（slab）描述。为避免相邻平板间的非物理相互作用，需在表面法向添加约 15 Å 真空层。为准确描述吸附质与表面的相互作用并减小镜像作用，通常采用 2×2 或 3×3 超胞，并固定底部两层原子模拟体相环境、放开顶层原子与吸附质弛豫。

2.3 位点辨识

Pt(111) 表面具有四种高对称吸附位点：顶位（atop，位于表面 Pt 原子上方）、桥位（bridge，两个相邻 Pt 原子之间）、fcc 穴位（下表层无 Pt 原子）与 hcp 穴位（下表层有 Pt 原子）。CO 在这些位点的吸附强度差异来自金属 d 轨道与 CO 前沿轨道（5σ、2π*）的重叠程度。

三、计算设置与模型构建

• 软件：VASP；泛函：GGA-PBE，PAW 赝势。
• 模型：Pt(111) 3×3 超胞，4 层原子平板，真空层 15 Å。
• 平面波截断能：400 eV；k 点网格采用 Monkhorst–Pack 方案，标准精度。
• 结构优化：底层两层固定，顶层与 CO 弛豫；能量收敛 1×10⁻⁵ eV，力收敛 0.02 eV/Å。
• 参考能量：气相 CO 在相同盒子中计算（仅 Gamma 点，偶极修正）。

四、结果分析与案例图解读

图 1 的左面板给出 CO 在 Pt(111) 四种位点的吸附能柱状图：顶位 −1.75 eV、桥位 −1.58 eV、fcc 位 −1.42 eV、hcp 位 −1.38 eV。计算结果表明 CO 最偏好顶位吸附，其吸附能绝对值最大，其次为桥位，两种穴位能量相近。该结果与低温实验（CO 在 Pt(111) 上主要占据顶位）及高精度计算文献一致。顶位吸附的优先性源于 CO 的 5σ 轨道与 Pt 原子 d 轨道沿轴向的强相互作用，而多配位穴位虽提供更多金属原子，但因轨道重叠方式改变反而结合较弱。四种位点的能量差异（约 0.4 eV）意味着在低温下 CO 会定域在顶位，而高温下可在位点间迁移，为理解 CO 在 Pt 表面的扩散与覆盖度依赖提供了定量基础。

图 1 的右面板给出 Pt(111) 表面的俯视示意图，用不同颜色与形状标记四种吸附位点：红色圆点表示顶位（正对表面 Pt 原子）、橙色方块表示桥位（两原子连线中点）、绿色三角表示 fcc 穴位、蓝色菱形表示 hcp 穴位。该示意图将抽象的位点概念可视化，配合左面板的能量数据，直观展示了”能量最低位点”与”几何位置”的对应关系，使读者能快速建立吸附位点空间结构与吸附能之间的关联。

两面板结合，左面板从能量角度筛选出最稳定吸附位点，右面板从几何角度定位位点空间分布，共同完成 CO/Pt(111) 吸附能的系统评估，为后续过渡态搜索、反应路径分析与催化活性预测提供输入结构。

五、结论与展望

本项目利用 VASP 计算了 CO 在 Pt(111) 表面四种位点的吸附能，确认顶位为最稳定吸附位（−1.75 eV），吸附能排序为 atop > bridge > fcc ≈ hcp，与实验和文献一致。案例图以吸附能柱状图与表面位点示意图双面板呈现。

后续可拓展：(1) 计算 CO 的吸附振动频率与实验红外对照，验证吸附构型；(2) 构建 CO 吸附的差分电荷密度，分析成键机制；(3) 以吸附能为基准，计算 CO 氧化反应的过渡态与活化能，构建完整催化循环。