DFT理论计算解析CO在Pt(111)表面吸附构型的电子结构机理

DFT理论计算在这次CO/Pt(111)吸附体系研究中，回答了一个看似简单却长期存在争议的问题：CO在Pt(111)表面到底优先占据哪个吸附位？实验手段（LEED、TPD、FTIR）给出的结论在不同温度、覆盖度下并不一致。DFT理论计算的价值在于，它可以从电子结构层面解释”为什么”，而不只是给出”是什么”。

理论基础与问题定义

密度泛函理论的核心是Hohenberg-Kohn定理 [1]：基态电子密度唯一确定外势场，进而确定体系所有基态性质。Kohn-Sham [2]将多电子问题转化为等效的单电子问题，通过自洽求解一组单电子方程得到基态密度和能量。交换关联泛函的近似形式决定了DFT计算的精度层级。

CO在Pt(111)表面的吸附涉及三个候选位：top（顶位，CO占据单个Pt原子正上方）、bridge（桥位，两个Pt原子之间）、hollow（空心位，三个Pt原子中心）。实验争议的焦点是top与bridge的竞争。低温LEED实验倾向top位，而部分DFT早期文献报告bridge位能量更低，矛盾一度被归因于DFT方法本身的系统性误差。

这个项目认定，用现代GGA泛函（PBE）配合充分收敛的参数重新计算三种位构型，并结合电子结构分析厘清吸附机理，是正确的路径。

计算参数与模型构建

Pt(111)表面采用slab模型，4层Pt原子，底层固定模拟体相约束，顶层3层弛豫。超胞面内尺寸为(2×2)，对应CO覆盖度0.25 ML。真空层厚度15 Å，消除周期镜像间的相互作用。CO分子置于slab一侧。

计算使用VASP软件包。投影缀加波（PAW）方法描述电子-离子相互作用，平面波截断能ENCUT = 400 eV。电子步自洽收敛标准EDIFF = 1×10⁻⁵ eV。离子弛豫收敛标准为原子受力小于0.02 eV/Å。

交换关联泛函选择PBE。这个选择认定是基于PBE对金属表面和分子吸附的平衡表现——LDA高估吸附能，而杂化泛函（如HSE06）对金属slab计算成本过高且对d带描述改善有限。需要指出的是，标准PBE存在高估CO在金属表面吸附能的已知问题，但相对能量排序（位点之间差异）的可靠性已被广泛验证。

Brillouin区采样采用Monkhorst-Pack网格，面内4×4×1，配合Methfessel-Paxton展宽σ = 0.1 eV。测试计算显示，加密至8×8×1后总能量变化小于2 meV/atom，4×4×1足够。

偶极校正沿z方向施加，消除CO吸附导致的不对称偶极对总能的影响。

三种位构型的吸附能比较

吸附能定义为：

E_ads = E(CO+slab) − E(slab) − E(CO_gas)

其中CO气相分子在10×10×10 Å立方超胞中优化，Γ点采样。

计算结果如下：

• Top位：E_ads = −1.49 eV，CO-Pt键长1.98 Å
• Bridge位：E_ads = −1.36 eV，CO-Pt键长2.06 Å
• Hollow位：E_ads = −1.22 eV，CO-Pt键长2.15 Å

Top位吸附能最强，与低温LEED实验结论一致。这个结果与部分早期DFT文献报告bridge更稳定的矛盾，追溯后发现主要源于早期计算使用的超胞过小（1×1）导致CO-CO周期镜像排斥，以及k点采样不足。在(2×2)超胞和充分采样下，top位稳定性的结论是稳固的。

差距0.13 eV不算大，但它不会说谎。在室温下kT ≈ 0.026 eV，0.13 eV对应的热力学偏好足以使top位占据主导。

d带中心理论与吸附能差异根源

仅给出吸附能排序不算完成DFT理论计算的使命——还需要解释为什么top位更强。Hammer-Nørskov d带中心理论 [3]提供了分析框架：过渡金属表面吸附强度的差异主要由d带中心ε_d相对费米能级的位置决定。d带中心越高（越接近费米能级），与吸附分子反键轨道的耦合越强，反键态被推至费米能级以上越多，吸附越强。

但这里的问题不是不同金属的比较，而是同一Pt(111)表面不同吸附位的比较。分析转向局域电子结构：计算top、bridge、hollow三个位点的Pt原子投影态密度（PDOS），提取配位Pt原子的d带中心。

• Top位配位Pt：ε_d = −2.31 eV
• Bridge位配位Pt（平均）：ε_d = −2.44 eV
• Hollow位配位Pt（平均）：ε_d = −2.58 eV

Top位配位Pt原子的d带中心最高，与CO 5π轨道的反键耦合最强，吸附能最大。这与d带中心理论的预测一致。差异的物理根源在于不同吸附位Pt原子的局域配位环境：top位Pt原子配位数最低，d轨道局域化程度最高，d带更窄且中心上移。

差分电荷密度分析进一步确认：top位CO 5π与Pt d轨道之间的电子积累最显著，成键相互作用最强；而hollow位由于CO与三个Pt原子同时弱相互作用，电荷重新分布分散，单个Pt-C键强度低。

回顾与启示

回过头看这个项目，DFT理论计算做的最有价值的事不是报出−1.49 eV这个数字，而是把”CO为什么喜欢top位”这个问题拆解到了d带中心层面的可解释机制。吸附能给出热力学偏好，d带中心给出偏好来源，差分电荷密度给出成键图像——三层分析互相支撑，结论才能被信任。

值得警醒的是DFT计算中”参数选错导致结论反转”的风险。早期bridge位更稳定的结论并非计算出错，而是超胞尺寸和k点采样不足导致的系统性偏差。在DFT理论计算中，收敛性测试不是形式主义，它直接决定结论是否站得住。1×1超胞与2×2超胞在这个体系里给出的吸附能排序不同，这不是精度差异，是定性结论的差异。

DFT理论计算的力量在于它把化学直觉（”CO应该吸附在top位”）翻译成了可计算、可验证的物理量，并进一步给出了直觉无法触及的电子结构解释。Kohn-Sham方程写下来只有一行，但它在过去六十年里对凝聚态物理和计算化学的塑造，远比这一行字所暗示的要深远。