VASP如何计算 CO₂RR *COOH/*CO 吸附能？

引言

CO₂RR 产 CO 路径看似只多了两个吸附态：*COOH 和 *CO。真正动手算时，麻烦通常出在三个地方：初始构型没有枚举，带电和自旋没有测试，气相参照前后不统一。*COOH 负责看 CO₂ 活化门槛，*CO 负责看 CO 脱附和位点中毒风险。本文就按这条线，把吸附能从建模、参照到复查讲清楚。

为什么 *COOH 和 *CO 要放在一条路径里看？

*COOH 到底卡住哪一步？

CO₂RR 产 CO 的常见写法是 CO₂ → *COOH → *CO → CO。这里的星号不是装饰，它表示中间体已经和表面活性位发生相互作用。*COOH 通常对应第一步质子–电子转移之后的吸附态，也就是 CO₂ 从线性分子被弯曲、被极化、开始进入反应路径的那一刻。

所以算 *COOH 时，不能只问“吸得牢不牢”。更准确的问题是：这个位点能不能把 CO2 拉成可反应的构型，同时又不把后续步骤压死。*COOH 形成自由能越低，通常说明 CO₂ 活化门槛越容易跨过，但低到过分也可能把路径锁在某个强吸附态上。

图1  NiM-HDAC 催化剂的 ECR 机理 DFT 计算，展示 *COOH、*H 竞争和 CO₂ 到 CO 的自由能台阶。

*CO 为什么不是越弱越好？

*CO 是产 CO 路径的后段节点。它吸得太弱，前面的 *COOH 断键后很难稳定成中间体；它吸得太强，CO 产物又不容易离开，表面位点会被 CO 占住。对于金属位、单原子位点和缺陷位点，这个平衡尤其敏感。

写吸附能时最怕只报一个 E_ads(*CO)，然后直接说选择性提高。更稳的说法，是把 *COOH 形成、*CO 稳定和 CO 脱附连在一起看。*COOH 决定入口是否顺，*CO 决定出口是否堵。入口和出口都合适，CO 路径才像一条能跑通的反应线。

图1里把 *COOH、*H 和 CO 生成自由能放在同一组比较中，正好说明这个思路：CO₂RR 不是孤立吸附能比赛，而是要同时躲开 HER 竞争和 CO 过强吸附。

*COOH 构型枚举应该从哪里开始？

初始构型为什么不能只放一个？

做 VASP 弛豫前，至少要把 *COOH 的连接方式列出来。常见初始构型包括 C 端靠近金属位、O 端靠近金属位、C/O 双齿配位、羟基 H 朝外或朝向表面、以及中间体沿不同晶向倾斜。对单原子位点，还要试一试 *COOH 是否会偏向邻近 N、C、O 或缺陷边缘。

式(1)  E_ads(*COOH) = E_slab+COOH – E_slab – E_COOH,ref

这些构型不是为了凑数量，而是为了避免把一个局部极小值当成全局趋势。*COOH 弛豫时很容易发生旋转、翻转、迁移，甚至从目标位点滑向邻位。构型枚举的目的，是先找到真正可比较的稳定吸附态，再谈吸附能高低。

图2中的 Fe-N₄ 位点示意、自由能和 *COOH 差分电荷可以一起看：结构图告诉你中间体坐在哪里，自由能告诉你路径是否顺，电荷分布再解释为什么这个吸附态能被稳定。

图2  Fe-N₄ 位点的优化结构、ECR 自由能图和 *COOH 吸附差分电荷分布。

*CO 构型为什么也要多试？

*CO 常被默认成 C 端顶位吸附，但这只是起点，不是结论。金属表面上可能出现 top、bridge、hollow 构型；二维单原子位点上也可能出现 CO 偏斜、C 端偏离金属中心、或 O 端被邻近原子弱相互作用牵住的情况。

真正要比较的是弛豫后的结构，而不是初始文件里摆得最漂亮的那一个。建议每个候选位点至少保存初始 POSCAR、CONTCAR、最终键长、吸附高度、C-O 键长和磁矩变化。*CO 的吸附构型会直接影响 CO 脱附判断，不能只拿一个随手放置的构型下结论。

*CO 的吸附能看起来差 0.1–0.2 eV 时，初始构型、k 点、偶极修正或自旋设置都可能改变排序。这个范围内不要把话写得太满，最好把“能量接近”写出来。

气相参照和自由能修正怎样统一？

Eads 的参照物怎么选才不乱？

吸附能公式看起来简单，真正麻烦的是参照态。E_slab+ads、E_slab 和气相分子必须来自同一套泛函、赝势、ENCUT、电子收敛标准；二维表面还要统一真空层、偶极修正和 slab 构型。否则 E_ads 的差别可能来自设置，而不是吸附本身。

式(2)  ΔG_*COOH = ΔE_*COOH + ΔZPE – TΔS – eU + ΔG_pH

CO₂、CO、H₂ 这类气相分子建议单独放在足够大的盒子里计算，使用Γ 点并检查分子构型是否合理。CO₂ 应接近线性，CO 和 H₂ 是闭壳层分子；如果开自旋后跑出奇怪磁矩，先不要急着拿来做参照。

CO₂RR 自由能通常会用计算氢电极处理 H⁺ + e^–。这时不要把所有中间体都硬拆成不存在的气相 *COOH 分子参照，而是要写清楚每一步采用的反应式。参照体系一旦前后混用，自由能台阶就失去可比性。

自由能修正要写到什么程度？

VASP 给的是电子总能，论文里讨论反应趋势时还要考虑 ZPE、熵、温度、电位和 pH。对吸附态而言，熵项通常比气相分子小很多；对 CO₂、CO、H₂ 这类气相物种，熵修正不能随便忽略。

如果只是做同一体系内部筛选，可以先用ΔE 看趋势；一旦要画 CO₂RR 自由能图，就要把ΔG 写完整。一个常见表达是ΔG_*COOH = ΔE_*COOH + ΔZPE – TΔS – eU + ΔG_pH。总能排序只能当初筛，自由能台阶才对应电化学路径判断。

图3中的 Sn-Bi 界面把不同路径的自由能和 PDOS 放在一起，提醒我们不要只盯某个中间体能量。路径竞争、电子态和产物脱附要同时看，CO 路径才不会被 HCOOH 或 HER 解释抢走。

图3  Sn-Bi 界面的 CO₂RR DFT 模型、CO/HCOOH 路径自由能和 PDOS 对比。

带电和自旋测试什么时候必须做？

哪些体系不能默认非自旋？

如果模型包含过渡金属、单原子位点、空位、边缘缺陷、自由基式中间体或未饱和配位环境，就不要默认 ISPIN = 1 直接跑完。*COOH 本身涉及质子–电子转移后的局域成键，某些位点还会出现局域磁矩变化；*CO 吸附也可能让金属 d 态重新分裂。

更稳妥的做法，是先用几个不同 MAGMOM 初值做测试，比较收敛后的总能、磁矩分布和关键键长。若最低能磁态不同，吸附能也要随之重算。自旋态改变的是电子占据，电子占据会反过来改变吸附强度，这不是后处理时随便补一句能解决的问题。

图4里的 Zn-N₄ 位点重构和自由能变化说明，局域配位对 CO₂RR 中间体稳定性影响很大。对于这类单原子或缺陷位点，磁性、配位畸变和吸附态必须一起复查。

图4  平面氯化诱导 Zn-N₄ 位点重构的 DFT 结构和 CO₂RR 自由能变化。

带电模型能不能直接拿来和中性模型比？

带电 slab 可以用来做电位、电场或电子注入的敏感性测试，但不要把不同总电荷的结果直接混在同一张吸附能表里。周期性边界条件下，带电体系会引入均匀背景电荷；如果 slab 不对称，还会叠加偶极和真空层依赖。

如果确实要比较带电结果，至少要固定同一电荷态、同一真空厚度、同一偶极修正和同一构型起点，并报告 NELECT 或总电荷设置。带电计算适合做趋势敏感性，不适合偷偷替代主线参照。

最后交付数据时，建议把每个吸附态的目录名写清楚：site、adsorbate、spin、charge、reference。比如 M1_COOH_Cbound_spin2_nelect-0.5 这种朴素命名，后期复查会少掉很多混乱。

表1  *COOH/*CO 构型枚举清单

表2  CO₂RR 吸附能复查表

本章要点总结

要点1：*COOH 主要用于判断 CO₂ 活化门槛，*CO 主要用于判断 CO 脱附和位点中毒风险。

要点2：构型枚举要覆盖连接方式、吸附位点、倾斜方向和弛豫后迁移，不能只取一个初始构型。

要点3：气相参照、slab 设置、ZPE/熵/电位修正必须前后一致，否则自由能台阶不可比。

要点4：过渡金属、单原子、空位和自由基式中间体要检查自旋；带电结果更适合做敏感性分析。