平均吸附能和逐步吸附能有什么区别？

引言

同一块表面依次放入两个、三个或更多吸附物时，一个平均值无法直接回答“新加入的这一个是否还能稳定吸附”。本文讨论平均吸附能、逐步吸附能和差分吸附能的定义，说明怎样从 VASP 总能序列识别覆盖度效应、横向相互作用与构型重排。

两个能量分别回答什么问题

平均值包含了哪些贡献？

设同一 slab 上有 n 个相同吸附物 X，E(nX*) 是优化后体系总能，E(*) 是干净表面总能，E(X) 是统一参照态能量。平均吸附能描述整个 nX* 构型的总体稳定性。它把从 0 到 n 个吸附物的总能变化除以 n，适合比较不同覆盖态的综合吸附强弱。

E_avg(n) = [E(nX*) – E(*) – nE(X)] / n　　式(1)

逐步吸附能只比较相邻两个覆盖态，用 E(nX*) 减去 E((n-1)X*) 和一个参照态。它回答第 n 个 X 加入现有覆盖层后的能量变化。逐步吸附能保留了吸附顺序和当前表面环境，因此更适合判断新增吸附是否仍有利，以及哪一步开始出现明显排斥或位点耗尽。

E_step(n) = E(nX*) – E((n-1)X*) – E(X)　　式(2)

图 1 直接比较平均吸附能与差分吸附能随覆盖度的变化。两条曲线斜率不同，说明同一批 VASP 总能可以导出两个物理问题：一条面向整体覆盖态，另一条面向最后加入的吸附物。若只保留 E_avg，第 n 步突然变弱的现象可能被前面较强的吸附贡献稀释。

两种定义的符号约定也要统一。有些工作把放热吸附写成负值，也有工作把结合强度写成正值；比较前先从公式确认符号。本文采用吸附放热为负的口径，因此 E_step(n) 越负，表示第 n 个吸附物加入后的热力学稳定化越强。

图1 CO/Pt(111) 的平均吸附能、差分吸附能及覆盖度相关脱附结果。两种能量对应不同的统计与动力学含义。DOI: 10.1021/acs.jctc.2c01005。

第n个吸附物怎样单独计价

相邻覆盖态必须满足什么关系？

式(2)成立的前提是 nX* 与 (n-1)X* 属于同一模型序列：slab 尺寸、固定层、泛函、赝势、k 点、ENCUT、电子展宽和能量读取口径都保持一致。第 n 个吸附物的参照态也不能中途更换。例如 H* 全部用 1/2E(H₂)，CO* 全部用 E(CO)，这样相邻能量差才只包含新增吸附与结构响应。

构型继承同样关键。计算第 3 个 H 时，应从最低能 2H* 覆盖态继续加入候选位点，同时保留若干接近的 2H* 构型作为分支。逐步能量只有在相邻覆盖态具有可追踪的结构关系时，才接近“第 n 个吸附物”的贡献。若 2H* 取分散构型、3H* 却从团簇构型起算，E_step(3) 会同时包含新增 H 与原有两个 H 的大幅重排。

图 2 给出 CO—CO 距离相关的吸附物相互作用势，并用 Pt(111) 位点示意说明邻居范围。短程排斥、表面介导相互作用和偶极作用都会改变新增吸附物的能量。若第 n 个 CO 被迫靠近已有 CO，逐步吸附能变弱可能来自横向排斥，并不等于单个活性位的本征成键能力突然下降。

若同一 n 下近邻和分散构型的能量差随 n 增大，应把构型差记入吸附物相互作用分析。只比较各覆盖度最低能结构时，这部分信息会被最低值筛选压缩，因此候选构型能量仍要保留。

图2 CO 吸附物间相互作用随距离的变化及 Pt(111) 邻近位点定义。逐步吸附能会响应局域排布。DOI: 10.1021/acs.jctc.2c01005。

构型序列为什么决定结果

同一覆盖度为什么不止一个能量？

同一 n 值往往对应许多位点组合。覆盖度只规定数量比例，不能唯一确定原子排布。2×2 金属表面放入两个 CO，可以形成近邻、对角分散、顶位—桥位混合等构型；二维材料上的 H 还可能分布在缺陷、边缘与基面。每个候选构型都要充分弛豫，并确认吸附物身份没有变化。

对热力学最低能序列，可以在每个 n 下选择最低能构型，再计算 E_avg 与 E_step。对动力学问题，还要保留亚稳构型，因为实验覆盖层可能受扩散能垒限制。图3中的吸附能分布说明，同一覆盖度不是单一数值，而是一组局域排布的稳定性分布；最低值、均值与分布宽度分别承载不同信息。

当新增吸附触发表面重构时，能量变化需要单独解释。以柔性金属表面为例，近邻吸附可能增加吸附物排斥，却通过表面畸变获得额外稳定；两项相互抵消后，平均能变化很小，而逐步能仍能暴露当前加入步骤的变化。构型、表面形变和吸附物相互作用应在同一优化结构上共同判断。

图3 不同 CO 数量和排布形成的平均吸附能分布。分布宽度反映同一覆盖度下的构型差异。DOI: 10.1021/acs.jctc.2c01005。

平均值会掩盖哪些变化

什么时候应优先看逐步能量？

若前两个吸附物结合很强，第三个吸附物结合很弱，三者平均后仍可能得到负值。此时 E_avg(3) 只能说明 3X* 整体相对参照态稳定，不能推出第三个 X 的加入仍然容易。判断覆盖层能否继续增长、某一位点何时饱和或脱附从哪一覆盖度开始时，应优先读取 E_step(n)。

平均值也可能掩盖位点身份切换。低覆盖度 H* 占据缺陷位，高覆盖度新增 H* 进入基面位；E_avg 混合了两类位点，而 E_step 能标出切换发生在哪一步。平均能用于评价覆盖态，逐步能用于识别吸附序列中的转折点。对 CO₂RR 的 *CO 覆盖层，新增 *CO 的差分能还能与 *CO 脱附、加氢和 C—C 偶联竞争相联系。

图 4 比较多个覆盖度的程序升温脱附曲线，表明脱附行为对覆盖度和吸附能分布敏感。微观动力学需要的是与当前覆盖状态相匹配的化学势或差分能，而不是把所有覆盖度都替换成一个平均常数。若只讨论静态构型稳定性，平均吸附能仍然有用，但要同时注明 n、表面积与位点排布。

还有一种情况是 E_avg 变化平缓，而 E_step 在某个 n 发生跳变。先检查该覆盖度是否出现位点切换、吸附物成对、表面原子抬升或磁态改变；若结构保持同类，跳变才主要反映吸附物间相互作用。能量转折必须回到对应优化结构确认物理来源。

将 E_step 与吸附物间最短距离同时作图，可以区分近邻排斥和位点耗尽。若距离缩短时逐步能持续变弱，横向排斥是合理解释；若距离近似不变但能量突变，则应继续核对表面重构、电子占据或新位点类型。

图4 不同 CO 覆盖度的程序升温脱附曲线。覆盖度相关动力学需要差分能或能量分布。DOI: 10.1021/acs.jctc.2c01005。

VASP算例怎样组织

总能序列怎样保持一致？

算例从一个统一母体开始：`00_clean` 保存干净 slab，`01_1X` 到 `06_6X` 分别保存不同 n 的候选构型。每个 n 建立位点组合，完成结构优化后用一致参数做静态单点。金属表面可测试合适的 ISMEAR 与 SIGMA，半导体表面要核对带隙和占据；所有能量都取最后一个完成的静态步对应的同类能量项。

参数可采用 PBE/PAW 作为起点，ENCUT、k 点、slab 厚度与真空层均需按体系收敛；同一序列保持固定层和偶极修正设置一致。结果记录包含 n、覆盖度 θ、构型编号、E(nX*)、E_avg(n)、E_step(n)、关键间距和表面位移。若涉及自由能，再为每一步加入统一的 ΔZPE 与 -TΔS 口径。

数值核对可利用恒等关系：nE_avg(n) 等于前 n 个逐步吸附能之和。平均吸附能与逐步吸附能应由同一组总能同时生成。若两边不相等，通常来自参照能、构型编号或能量列混用；统一生成也能让单位、符号和参照态一次闭合。

nE_avg(n) = Σ_i=1ⁿE_step(i)　　式(3)

若吸附物之间不完全等价，还应把“加入顺序”明确写入构型编号。例如先占据缺陷位再占据基面位，与先占据基面位再占据缺陷位，最终可能弛豫到同一覆盖态，也可能保留不同局域极小点。前者若收敛到同一结构，只保留最低能序列；后者若保持不同结构，则分别计算逐步能，并把它们视为两条吸附路径。这样可以判断差异来自位点顺序还是覆盖度本身。

自由能处理中，E_step 还要加入第 n 步对应的振动和熵修正。对 H*，新增一个 H 的参照常取 1/2G(H₂)；对电化学吸附，还可能加入 -eU 与 pH 项。所有修正都应在相邻覆盖态之间做差，不能把 n 个吸附物的总修正直接当成第 n 步修正。逐步自由能的每一项都来自两个相邻覆盖态之差，这一口径能直接接入覆盖度相关的微观动力学。

若吸附物会发生缔合或解离，n 的计数单位也要保持统一。两个 H* 形成 H₂ 后，最终结构已经不再是单纯的 2H* 覆盖态；两个 CO 形成 *OCCO 后，也应转入反应路径数据，而不继续作为普通 2CO* 平均吸附能。结构身份一旦变化，就为新物种单独定义能量参照。

本章要点总结

要点：平均吸附能评价整个覆盖态的平均稳定性，逐步吸附能评价第 n 个吸附物加入后的变化。

要点：相邻覆盖态必须使用同一 slab、同一参照态和可追踪的构型序列。

要点：覆盖层增长、位点饱和和微观动力学更依赖逐步或差分能量。