DFT计算吸附能：d带中心理论与吸附构型搜索实战

DFT计算吸附能时，项目组在Pt(111)/CO体系上吃过一次亏：初始只算了top位吸附，结果吸附能约-1.5 eV，与实验值（约-1.35至-1.55 eV）偏差不大就交了差。后来审稿人要求补充其他位点，补算后发现PBE给出bridge位的吸附能比top位更负——约低0.1-0.2 eV。CO在Pt(111)上优先吸附在top位是实验共识，DFT却给出bridge位更稳定。问题出在泛函：PBE对CO的描述有系统性偏差，这个问题在文献中被称为”CO/Pt puzzle”。

吸附能的定义与计算

吸附能的定义是：

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

负值表示放热吸附（稳定），正值表示吸热（不稳定）。这个公式看似简单，但每一项都有坑。

E(surface)的计算：需要构建slab模型。Pt(111)表面用4层slab，上2层弛豫、下2层固定。真空层15 Å（避免周期镜像间相互作用）。K-mesh用5×5×1（面内），Γ-centered。

E(adsorbate)的计算：CO分子在10×10×10 Å的立方盒子中优化，Γ点采样。CO的键长优化到约1.144 Å（实验值1.128 Å，PBE系统性偏长约1.4%）。

E(surface+adsorbate)的计算：CO放在slab上方，上2层Pt+CO一起弛豫。这里有一个微妙的点——slab的层数。项目组测试过3层、4层、5层slab对吸附能的影响：

3层到4层变化约0.1 eV，4层到5层变化约0.02 eV——4层已基本收敛。用3层slab会导致吸附能低估约0.1 eV。

d带中心理论：预判吸附位点

在搜索吸附构型之前，项目组先用d带中心理论做预判。Hammer-Nørskov模型指出：吸附能与金属表面d带的中心位置ε_d近似线性相关——d带中心越靠近费米能级，吸附越强。

通过VASP计算Pt(111)表面不同位点（top、bridge、hollow）的投影态密度（PDOS），提取d带中心：

根据d带中心理论，top位应该吸附最强。但PBE计算给出的bridge位吸附能更负（比top低约0.1-0.2 eV），与理论预判矛盾。这就是”CO/Pt puzzle”的核心：PBE高估了CO的5σ轨道在bridge位的反馈键贡献。

解决方案是使用包含范德华校正的泛函或rVV10校正。用optPBE-vdW或rVV10重新计算后，top位和bridge位的吸附能差距缩小到约0.03 eV以内，在DFT精度范围内可以认为两者竞争吸附。这与实验观察到的top位优先（但bridge位也有可观布居）一致。

构型搜索：6种初始构型的系统筛选

项目组对每种吸附分子执行系统构型搜索。以CO在Pt(111)上为例，6种初始构型：

1. top位垂直吸附（C端朝下）：最常见构型
2. bridge位垂直吸附：C在两个Pt之间
3. hollow(fcc)位垂直吸附：C在三层Pt原子空位上方
4. hollow(hcp)位垂直吸附：C在三层Pt原子空位上方（另一种堆垛）
5. top位倾斜吸附：CO轴线偏离法线30°
6. bridge位倾斜吸附：CO轴线偏离法线20°

每种构型独立优化，取能量最低者为最稳定构型。倾斜构型的意义在于：某些分子（如O2、NO）在金属表面倾向于倾斜或平躺吸附，如果只算垂直构型会错过最稳定状态。

项目组在O2/Pt(111)体系上验证过这个策略：只算垂直top位得到E_ads约-0.85 eV，补充平躺bridge位后得到E_ads约-1.2 eV——差了约0.4 eV，构型搜索不足会导致吸附能严重低估。

覆盖度效应：孤立吸附 vs 高覆盖

上述计算都是1/4 ML覆盖度（2×2超胞中1个吸附分子）。实际催化条件下覆盖度可能更高，吸附分子间的排斥会显著影响吸附能。

项目组测试了1/4 ML、1/2 ML和1 ML覆盖度下CO/Pt(111)的吸附能：

从1/4到1 ML，每增加一个CO分子，吸附能减弱约0.5 eV——这是CO之间的偶极排斥和Pauli排斥导致的。如果只算低覆盖度的吸附能来推断高覆盖度下的催化活性，会高估催化剂的活性。

BSSE校正：是否必要

基组重叠误差（BSSE）在平面波DFT中通常很小，因为平面波在空间中是正交的。但在计算小分子在表面上的弱吸附（如CH4在石墨烯上的物理吸附，E_ads~0.1 eV）时，BSSE仍然不可忽略。

项目组用counterpoise correction方法估算BSSE：在CO/Pt(111)体系中BSSE约0.03 eV——可以忽略。但在CH4/graphene体系中BSSE可达0.08 eV——占吸附能的80%，必须校正。

反思：吸附能不是催化活性的唯一指标

DFT吸附能是热力学量，描述的是吸附的平衡态。但催化反应涉及多个中间体的吸附-脱附竞争——Sabatier原理指出催化活性在”吸附不太强也不太弱”时最优。项目组在实际催化筛选中，不只看单一吸附能，而是构建反应能图（reaction energy diagram），追踪从反应物到产物的每一步能量变化。

d带中心理论提供了一个快速的预筛选指标，但它只是零级近似。精确的吸附能计算仍然需要系统构型搜索+合理的泛函选择+覆盖度校正。更多表面催化计算的实战经验，可以参考VASP/第一性原理栏目，或访问科研学术网首页。