Cu 表面 CO₂RR 最绕人的地方之一，是 *CO 的身份会变。它可以是产 CO 路径的终点，也可以是生成 CH₄、CH₃OH 或 C₂ 产物的中间体；在覆盖度过高或后续加氢太慢时，它又会变成阻塞位点的物种。判断 *CO 是中间体还是毒化物，不能只看 *CO 吸附能一个数，还要同时看覆盖度、位点、后续加氢和 C-C 偶联。
图1 展示了 CO 和 OH 覆盖度如何影响 Cu(111) 与 Cu(100) 的表面稳定性。对 VASP 建模来说，这说明 *CO 不是孤立存在的单点吸附问题；在真实反应电位下，表面可能被 *CO、*OH 或其他物种共同重塑。若模型只算一个 *CO@Cu(111)，很难回答它到底是有用中间体还是堵位物种。
图1 DOI: 10.1038/s41467-024-45538-y
第一步不是给 *CO 贴标签，而是先问它处在什么路径位置。若目标产物是 CO，*CO 更接近产物前体；若目标是 CH₄、CH₃OH 或 C₂ 产物，*CO 就是继续加氢或偶联的入口；若它在关键位点上停留太久，又可能变成阻塞物。同一个 *CO，在不同覆盖度、电位和晶面上可能换身份。
图2 给出了 CO 在 Cu(111) 不同位点的吸附构型、吸附能、Bader 电荷和键长。它提醒我们，top、bridge、fcc 和 hcp 位点不能混在一起比较。最低能位点只能说明“这个构型站得住”，不能自动推出“这条路径走得通”。
图2 DOI: 10.1021/acs.jpcc.4c07418
这时再写吸附能才有意义。它是身份判断的入口指标，不是最终结论。若 *CO 吸附太弱，可能还没来得及加氢就脱附；若吸附太强，可能让后续步骤或位点更新变慢。
式(1) Eads(*CO) = E*CO/Cu – ECu – ECO(g)
这个表达式是 *CO 吸附能的常用写法。吸附能越负通常表示 *CO 越稳，但这只说明它留在表面的倾向；是否继续反应，还要看 *CO → *CHO、*CO → *COH 或 *CO + *CO → *OCCO 等后续步骤。
低覆盖度模型适合看单个 *CO 与 Cu 位点的本征相互作用，高覆盖度模型才接近电解过程中 CO 富集的表面。覆盖度升高后，*CO 之间会出现排斥、偶极耦合和位点竞争；表面也可能从 Cu(111) 偏向 Cu(100) 或台阶富集。如果只用一个孤立 *CO 判断毒化，结论很容易过度简化。
图3 展示了双 CO 吸附和 OCCO 二聚构型。对 Cu 表面 CO₂RR 来说，这类构型非常重要：如果两个 *CO 靠近后能够形成 *OCCO，*CO 就是 C₂ 路径入口；如果两个 *CO 只是互相排斥并抬高后续能垒，它更像覆盖层阻塞物。
图3 DOI: 10.1021/acs.jpcc.4c07418
覆盖度建议写成模型条件，而不是写成一句笼统描述。比如 3×3 Cu(111) 上 1 个 *CO、2×2 Cu(100) 上 1 个 *CO、同一超胞里 2 个相邻 *CO，这三类模型回答的问题不同：前者偏向单分子吸附，中间情况适合比较晶面，后者才开始触及偶联或阻塞。
式(2) θCO = N*CO / Nsite
这里 θCO 是 CO 覆盖度，N*CO 是表面 *CO 数量，Nsite 是可用活性位数量。比较不同覆盖度时，超胞大小、位点定义和每个 *CO 的平均吸附能要写清楚。
*CO 后续常见分叉有三类。第一类是 *CO 加氢到 *CHO，常被用来讨论 CH₄ 或 CH₃OH 路径；第二类是 *CO 加氢到 *COH，在部分模型中可能更接近 C-C 偶联前的活化；第三类是两个 *CO 偶联成 *OCCO，通向 C₂ 产物。*CO 能否继续走下去，要看这些后续台阶，而不是只看 *CO 本身稳不稳。
图4 展示了原位光谱中 *CO atop 和 *CO bridge 信号随条件变化的趋势。计算上可以对应两类构型：线性吸附的 *CO 和桥式吸附的 *CO。若实验显示桥式 *CO 积累而产物分布变差，VASP 中就应检查桥式 *CO 的覆盖度、迁移和后续加氢是否卡住。
图4 DOI: 10.1021/jacs.5c02954
如果只想先做一个最小判断，可以从 *CO → *CHO 这一步开始。它不覆盖所有路径，却能快速区分“*CO 留在表面后能否进入深度还原”。
式(3) ΔG*CO→*CHO = G*CHO – G*CO – 1/2GH₂
这个自由能差表示 *CO 加氢到 *CHO 的难度。若这一步很高，而 *CO 脱附也不容易，*CO 就可能在热力学上积累；若这一步较低，*CO 更像深度还原中间体。实际电化学条件下还要加入电位项和振动校正。
在计算语境里，不建议只因为 *CO 吸附能很负就说毒化。更稳的判断需要三个条件：它占据关键 Cu 位点，它的脱附或后续转化都不容易，它还会抬高其他关键中间体的形成自由能。若 *CO 稳定但能快速偶联或加氢，它反而是有价值的中间体。
图5 给出了 CO₂RR 后续 C-C 偶联和加氢步骤的 DFT 结果示例。它说明 *CO 后面的路径并不止一条：有些表面会把 *CO 引向乙烯，有些会引向乙醇，还有些会停在 CO。写结论时应说明你比较的是哪条分支，而不是笼统说 *CO 好或不好。
图5 DOI: 10.1002/anie.202520291
脱附自由能可以作为反证之一：如果 *CO 不容易脱附，后续加氢和偶联也不顺，它才更接近阻塞物。若 *CO 脱附困难但 *OCCO 或 *CHO 台阶很低，它仍可能只是反应链上的高覆盖度中间体。
式(4) ΔGdes(*CO) = GCO(g) + G* – G*CO
这个量可以用来衡量 *CO 脱附倾向。ΔGdes 越大，*CO 越不容易离开表面。但是否构成阻塞，还要和 *CO 后续加氢、*H 竞争吸附和 *CO 覆盖度一起判断。
建议把自己的 VASP 结果按同一套条件放在一起。最小列可以包括：晶面、位点、覆盖度、Eads(*CO)、ΔGdes(*CO)、ΔG*CO→*CHO、ΔG*CO→*COH、ΔG*OCCO 和 ΔG*H。这些列不是越多越好，而是要对应你的目标产物。
如果目标产物是 CO，判断重点放在 *CO 脱附和 *COOH 形成之间的平衡。*CO 太弱会导致 *COOH 难以稳定或 CO 过早离开，*CO 太强又会造成脱附困难。此时重点比较 ΔG*COOH、Eads(*CO) 和 ΔGdes(*CO)，不要只拿最低 *CO 吸附能排序。
如果目标是 C₂ 产物，覆盖度和相邻 *CO 构型要提前进入矩阵。低覆盖度 *CO 可能只说明单分子吸附稳定，高覆盖度下两个 *CO 的距离、取向和偶极排斥才决定 *OCCO 是否可行。建议至少保留同一 Cu 面上的低、中、高覆盖度模型，并固定相同表面层数、K 点和色散修正。覆盖度变化本身就是反应条件的一部分。
图6 给出了不同电流密度下 CO₂/CO 还原产物分布的例子。它适合放在这一节，是因为计算矩阵最后要回到选择性：如果某组模型预测 *CO 更容易进入 C₂ 路径，就应能和产物分布、原位谱峰或覆盖度变化形成互相印证，而不是只停在单个吸附能数值上。
图6 DOI: 10.1038/s41467-024-45538-y
竞争物种也要留一个入口。CO₂RR 表面往往同时存在 *H、*COOH、*OCHO、*OH 或水分子，*CO 是否阻塞位点，要看它是否挤掉了更关键的反应入口。若 *CO 占位后 *H 或 *COOH 形成自由能显著升高，它更接近阻塞物；若 *CO 与邻近 *H 共同降低 *CHO 台阶，它仍可能是可用中间体。
写结论时可以采用条件句：在低覆盖度 Cu(111) 上，*CO 更像可脱附产物；在 Cu(100) 或台阶位高覆盖度下，*CO 可能进入偶联路径；若后续 *CHO、*COH 或 *OCCO 台阶都偏高，则表现为位点阻塞。这样的表达比简单说“*CO 有利”或“*CO 毒化”更接近计算数据本身。
Cu 表面 *CO 的身份取决于吸附强度、覆盖度和后续反应能力。低覆盖度下的 *CO 吸附能只能回答它是否愿意留在 Cu 上；高覆盖度、双 CO 偶联、*CHO/*COH 加氢、CO 脱附和 *H 竞争才决定它是中间体还是阻塞物。稳妥写法是把 *CO 相关台阶拆开，用同一套参照比较，而不是用一个吸附能给全部路径下结论。
