吸附能计算分析：CO 在 Pt(111) 表面的 DFT 研究

一、案例背景与科学问题

CO 在过渡金属表面的吸附是多相催化、燃料电池阳极抗毒化、汽车尾气催化净化等领域的基础问题：

• Pt 是性能优异的催化材料，但其表面容易被 CO 强吸附而中毒失活；
• 吸附能的大小决定了吸附物种在表面的稳定性、迁移性与反应活性；
• 准确计算吸附能是后续催化反应机理研究（如 CO 氧化、水煤气变换）的第一步。
• 

本案例要回答的科学问题：

1. CO 在 Pt(111) 表面更倾向吸附在哪个位点（顶位/桥位/穴位）？
2. 各吸附位点的吸附能差异有多大？
3. 吸附过程中电荷如何重新分布、成键机制是什么？

二、计算模型与 DFT 方法

2.1 表面模型（Slab 模型）

• 基底：切取 Pt(111) 表面，采用 2×2 超胞，共 4 层原子；
• 底层 2 层固定（模拟体相），顶层 2 层弛豫（允许结构优化）；
• 层间加入 15 Å 真空层，消除周期性镜像相互作用；
• 采用偶极校正消除表面非对称带来的虚假电场。

2.2 交换关联泛函与参数

数值均为示意参数，实际计算需通过收敛性测试确定。

2.3 吸附能定义

吸附能定义为吸附体系总能量与裸表面、自由分子能量之差：

$$ E_{ads} = E_{slab+CO} – E_{slab} – E_{CO(g)} $$

• $E_{slab+CO}$：CO 吸附后体系总能量；
• $E_{slab}$：裸 Pt(111) 表面能量；
• $E_{CO(g)}$：气相 CO 分子能量（置于相同尺寸盒子中计算，做自旋与盒子尺寸修正）。

约定：吸附能为负值表示吸附放热、热力学有利；数值越负，吸附越强。

三、吸附构型与位点分析

在 Pt(111) 密排面上存在三类高对称吸附位点（见案例图 a 面板）：

对每种位点分别构建初始构型（C 端朝下的垂直构型），进行结构弛豫，得到稳定吸附构型。

• 弛豫后 CO 分子始终保持 C 端向下（C 与表面成键），这与实验观测一致；
• 吸附高度 d_ads ≈ 2.0 Å（C 原子到最上层 Pt 面的距离）。

四、吸附能计算结果与位点稳定性

结构弛豫收敛后，代入吸附能公式得到各吸附位点的能量（见案例图 b 面板）：

结果分析：

1. 全部位点吸附能均为负值 → CO 在 Pt(111) 表面吸附是自发放热过程；
2. 稳定性排序：hollow > bridge > top，穴位吸附能比顶位低约 0.6 eV，说明 CO 更倾向占据穴位；
3. 该趋势与文献中 CO/Pt(111) 体系”穴位最稳、顶位次稳”的结论定性一致（实际 DFT 结果因泛函不同，顶位/穴位相对稳定性存在争议，是 GGA 泛函的已知局限，可通过杂化泛函或实验校正）；
4. 吸附能差异（0.20.6 eV）对应了 CO 在表面的迁移势垒量级，影响表面扩散与催化活性。

物理解释：位点配位数越高，CO 的 2π* 反键轨道与表面 d 带的杂化越强，形成更稳定的吸附键——这正对应下文的电荷密度分析。

五、差分电荷密度分析（吸附成键机制）

为揭示吸附成键本质，计算差分电荷密度：

$$ \Delta\rho = \rho_{slab+CO} – \rho_{slab} – \rho_{CO} $$

其中 $\rho_{slab}$、$\rho_{CO}$ 为吸附前各自独立体系的电荷密度（保持与吸附体系相同的原子坐标）。

案例图 c 面板展示了沿表面法线方向的平面平均差分电荷分布 $\Delta\rho(z)$：

• 界面区（C–Pt 之间，z≈1.5 Å）出现明显的电荷积累（正值峰）：表明 C 与表面 Pt 之间形成共价键，电子在成键区域富集；
• CO 分子内部（z≈2.1 Å）出现电荷耗尽（负值谷）：CO 的 5σ 轨道向金属 d 带捐赠电子，同时金属 d 电子回授到 CO 的 2π* 反键轨道，即经典的 σ-捐赠 / π-回授机制；
• O 原子周围再次出现电荷积累（z≈2.65 Å）：回授电子主要分布在 CO 的 π* 轨道上；
• 表层 Pt 原子下方出现轻微电荷耗尽：表面金属态参与成键后电子重排。

这一 “积累—耗尽—再积累” 的振荡特征，是化学吸附（区别于弱物理吸附）的典型指纹。

六、结论

本案例通过 DFT 计算完整分析了 CO 在 Pt(111) 表面的吸附行为：

1. 吸附能：CO 在穴位吸附最稳定（E_ads = −2.21 eV），排序为 hollow > bridge > top，均为放热吸附；
2. 吸附构型：CO 以 C 端垂直吸附，吸附高度约 2.0 Å；
3. 成键机制：差分电荷密度揭示 C–Pt 界面电荷积累 + CO 内部电荷耗尽的 σ-捐赠/π-回授特征，证实为化学吸附；
4. 方法学意义：吸附能计算是催化机理研究的基石——准确的吸附能是构建微观动力学模型、预测反应路径与活性的前提。

附：方法流程速览

构建 Pt(111) slab 模型 → 收敛性测试（k点/截断能）
→ 优化裸表面 → 构建 3 种位点吸附构型
→ 结构弛豫 → 静态自洽
→ 代入 E_ads = E_slab+CO − E_slab − E_CO
→ 差分电荷密度分析 → 结论与机理讨论