4. 吸附能计算的标准流程

吸附能是表面吸附第一性原理计算中最核心的物理量之一，它把“分子/原子是否愿意停留在表面、停留得有多牢”转化为一个可比较的数值。但吸附能并不是一次简单的能量相减：从构型枚举到参考态选择，每一步都会影响最终数值的可信度。本节按实际计算顺序，把这条链路拆开讲清楚，并给出可复现的操作要点与判断依据。

一、先明确吸附能定义与参考态

吸附能最常用的定义是：吸附体系总能量减去“洁净表面能量 + 孤立吸附质能量”。以符号表示，若吸附质为分子 A、表面为 S，则

E_ads = E(A/S) - E(S) - E(A)

约定：E_ads < 0 表示吸附放热、体系更稳定；
      E_ads > 0 表示吸附吸热、需要外界供能。

这里必须强调符号约定：不同文献可能定义成相反符号，阅读和汇报时必须先声明。本文统一采用“越负越稳定”的约定。参考态的选择同样关键：孤立吸附质能量应取其在同一计算设置下的气相（或自由态）能量，而不是实验值或另一套泛函的结果，否则误差会直接进入吸附能。若吸附质是 O₂、H₂ 等分子，还要注意自旋极化与基态构型；若吸附质是单个原子，则参考态为孤立原子能量，并需说明是否做了自旋处理。

二、吸附构型枚举：先穷举，再筛选

构型枚举的目标是覆盖所有物理上合理的吸附位点与取向，避免只算一个“看起来顺眼”的构型就下结论。对常见金属或氧化物表面，可按以下维度枚举：

• 吸附位点：顶位（top）、桥位（bridge）、空心位（hollow，fcc/hcp）、穴位、氧空位等；
• 分子取向：对多原子分子，考虑平行、垂直、倾斜以及不同原子朝向表面；
• 覆盖度：单吸附质、成对吸附、不同覆盖度（如 1/4 ML、1/2 ML）会改变侧向相互作用；
• 表面重构：某些体系吸附后表面原子会弛豫或重构，必要时把重构纳入候选。

枚举后不必全部做高精度弛豫。实用做法是：先用较粗的 k 点网格和较低截断能对所有候选做快速筛选，保留能量最低的若干构型，再用高精度参数复算。这样既控制机时，又不漏掉真正的基态构型。

三、初始结构搭建：把“物理图像”变成可算模型

初始结构搭建决定计算能否顺利收敛。以 3×3 的 fcc(111) 表面为例，典型 slab 模型包含 4 层原子，底部 2 层固定为体相位置，顶部 2 层允许弛豫；真空层厚度取 15 Å 以上，避免周期性镜像之间的人为相互作用。吸附质初始位置放在目标位点上方约 2.0 Å 处，并保证其与表面最近原子距离不小于 1.5 Å，防止初始重叠导致力发散。

搭建时还要注意：晶格常数应与体相优化结果一致；slab 的原子层数要足以让表面效应衰减；若研究带电吸附或极性表面，需考虑偶极修正。初始结构不是“随便摆一个”，而是把化学直觉（成键方向、配位数、位阻）编码进坐标。

四、结构弛豫：让原子找到自己的位置

结构弛豫是在势能面上寻找局部极小值。实际操作中，先固定底部若干层原子，放开表面层与吸附质，让它们自由弛豫。收敛判据通常同时看三个量：最大力小于 0.02 eV/Å、能量变化小于 1×10⁻⁵ eV、最大位移小于 1×10⁻³ Å。三者不必同时达到极严，但力判据最关键，因为它直接反映结构是否还在“动”。

弛豫策略上，建议分两步：先用较软的参数（如截断能 400 eV、k 点 3×3×1、力判据 0.05 eV/Å）快速得到近似构型，再用较硬参数（截断能 500 eV、k 点 5×5×1、力判据 0.02 eV/Å）精修。若弛豫过程中出现原子飞散、能量震荡或自洽不收敛，优先检查初始结构是否重叠、真空层是否足够、自旋设置是否合理，而不是盲目增加迭代步数。

五、能量收敛判据：不只是“算完了”

能量收敛包含两层含义：电子自洽循环收敛和离子弛豫收敛。电子步要保证总能变化小于设定阈值（如 1×10⁻⁵ eV），离子步要保证力和位移达标。对吸附能而言，还有一个容易被忽略的收敛维度——对计算参数的收敛测试：截断能、k 点网格、真空层厚度、slab 层数都会影响 E_ads。工程上通常要求吸附能随这些参数的变化小于 0.05 eV，才认为结果可用于比较。

六、参考态选择：最容易出错的一步

参考态选择直接决定吸附能的物理含义。常见参考态包括：孤立中性原子、气相分子（如 O₂、H₂、CO）、体相金属、以及某些反应中使用的“化学势参考”。选择哪一种，取决于你想回答什么问题。例如研究 O 原子吸附时，参考态可以是 1/2 E(O₂)，也可以是 E(O) 孤立原子，两者给出的数值不同，但都能在各自语境下使用，关键是全文一致并明确说明。

对于分子吸附，还要注意参考态的自旋与基态。O₂ 基态为三重态，若计算时误设为单重态，参考态能量偏高，吸附能会被系统性高估。类似地，H₂、CO、NO 等分子都有各自的基态电子结构，必须在参考态计算中正确设置。

七、一个可复现的示例流程

以 CO 在 fcc(111) 表面顶位吸附为例，示例参数如下（属于教学示例，非实测结果）：slab 为 3×3、4 层，底 2 层固定；真空层 15 Å；截断能 500 eV；k 点 5×5×1；力收敛 0.02 eV/Å。枚举顶位、桥位、fcc 空心位、hcp 空心位四种构型，粗筛后对最低两个构型精修。参考态为气相 CO 分子，放在 15 Å 立方盒子中，自旋极化开启。最终取能量最低构型，按 E_ads = E(CO/S) − E(S) − E(CO) 计算，并报告位点、取向、覆盖度及全部计算参数。

把上述流程固定下来，吸附能计算就从“一次能量相减”变成一条可检查、可复现、可比较的工程链路。构型枚举保证不漏掉基态，初始结构保证收敛顺利，弛豫与收敛判据保证数值可信，参考态选择保证物理含义清晰。四者缺一不可。