双原子催化剂做 CO₂RR 时，很多误判来自“位点”没有定义清楚。两个金属原子可能分别承担 CO₂ 活化、C 端吸附、O 端稳定和质子转移环境调节；同一个 *COOH 或 *OCHO，也可能是 端位、桥位或偏向单金属的构型。先把中间体吸附方式拆成可比较的构型集合，比直接讨论最终产物更可靠。
两个金属原子提供了什么？
M₁-M₂/N-C 或 M₁M₂-N₆ 这类模型中，两个金属原子可以形成相邻配位中心。一个金属可能更容易与 C 端相互作用，另一个金属可能更稳定 O 端或调节局域电荷；当 CO₂ 弯曲、*COOH 形成或 *OCHO 桥连时，双金属几何会直接改变中间体的稳定性。
图1 展示的 Zn₁Sn₁/SNC 属于双原子位点与掺杂碳载体耦合的材料对象。与单金属 M-N₄ 不同，金属组合、金属间距、配位 N 数和邻近掺杂环境 都会改变吸附结果。双原子位点不能只写成“一个金属中心 + 一个吸附物”。
为什么 *COOH 与 *OCHO 要一起看？
*COOH 通常连接 CO 生成路径，*OCHO 通常连接甲酸盐或甲酸路径。二者的区别不只是 H 加到 C 还是 O 上，也包括中间体通过 C 端、O 端或 O-O/C-O 桥连方式与双金属位点相互作用。若只优化一个 *COOH 构型，就无法判断 CO 路径和甲酸路径的早期选择性。
一个对比情境是，M₁ 端位 *COOH 比桥位 *COOH 稳定，但 *OCHO 在 M₁-M₂ 桥位更稳定。另一个情境是，*COOH 与 *OCHO 能量接近，但后续 *CO 脱附或 HCOO* 转化差异明显。CO₂RR 选择性要从首个质子化中间体和后续脱附/转化步骤共同判断。
图1 Zn₁Sn₁/SNC 双原子催化剂的结构设计、合成示意和显微表征。图中给出了双金属位点与 N/S 掺杂碳载体的材料对象。DOI: 10.1038/s41467-
025-57573-4。
载体模型先定哪些参数？
氮掺杂石墨烯或多孔碳模型通常需要足够大的平面超胞，避免相邻双原子位点周期镜像过近。可从 5×5 或更大石墨烯超胞起步，移除相邻 C 原子形成双空位，再用 N 原子替换边缘 C，最后锚定 M₁ 和 M₂。若论文对象含 S、P 或 B 等邻近掺杂，也要把掺杂位置纳入模型变量。
图2 的原子级结构和配位环境提醒我们，双原子模型不能只看俯视图。载体尺寸、配位原子和双金属间距 决定了 *COOH/*OCHO 构型是否具有可比性。应记录 M₁-M₂ 距离、M-N 键长、金属是否突出碳平面、载体是否弯曲，以及吸附前后的位点重构。
金属组合怎样进入自旋和电荷？
Fe、Co、Ni、Cu、Zn、Sn 等金属在 N-C 载体中的局域电荷和磁矩差异很大。结构优化前应设置自旋极化，并给过渡金属合理 MAGMOM 初值；优化后要检查局域磁矩是否连续，避免把不同磁态下的能量混为同一构型比较。对含 Zn、Sn 这类组合，还要关注 p 轨道与金属 d 态或载体 p 态的耦合。
若比较 Zn-Sn、Fe-Ni 和 Co-Cu 三类双原子位点，不能把金属组合变化、配位环境变化和吸附构型变化同时混在一个结论里。更合适的做法是在同一载体和同一配位数下换金属，或在同一金属组合下改变 N/S 掺杂位置，逐步拆分变量。
图2 Zn₁Sn₁/SNC 的原子级结构和配位环境表征。该图可对应 VASP 中 Zn-Sn 双原子位点、N 配位和掺杂邻域的模型搭建。DOI: 10.1038/s41467-025-57573-4。
*COOH 有哪些端位和桥位？
*COOH 至少应比较三类初始构型：C 端连接 M₁、C 端连接 M₂，以及 C/O 或 O/H 参与的桥位构型。优化后要记录 C-M、O-M、C-O、O-H 键长和 CO₂ 弯曲程度。若 *COOH 优化后转成 *CO + *OH 或迁移到单金属位，说明该初始构型不稳定，应单独记录为路径变化而不是简单删除。
吸附能比较需要让 Eslab-X 和 Eslab 对应同一个双原子基底。对 *COOH 来说，气相参照常与 CO₂、H₂ 和 CHE 口径连接；若只比较相对构型，可先用同一参照筛选，再进入自由能修正。桥位 *COOH 的判据不是名称，而是优化后是否仍由两个金属共同稳定关键键合片段。
Eads(*X) = Eslab-X – Eslab – EX,ref 式(1)
*OCHO 怎样避免和 *COOH 混淆？
*OCHO 是 O 端主导的吸附物，常见构型包括单 O 端位、双 O 桥位和倾斜桥位。它和 *COOH 的元素组成相同，但氢的位置、C-H/O-H 键关系和与金属连接的原子不同。检查 CONTCAR 时应确认 H 是否连在 O 上或 C 上，不能只凭初始 POSCAR 命名。
图3 的产物选择性结果说明，CO₂RR 中 CO、甲酸盐和 H₂ 竞争会同时存在。放到 VASP 中，*COOH、*OCHO 和 H* 都应在同一双原子位点上比较。若 *OCHO 通过双 O 桥位稳定，而 H* 在同一金属上也很稳定，就要进一步判断 HER 竞争是否会占据该位点。
图3 Zn₁Sn₁/SNC 的 CO₂RR 电化学性能和产物选择性。性能结果提示双原子位点需要与 *COOH/*OCHO 路径选择联系分析。DOI: 10.1038/s41467-025-57573-4。
首个质子化步骤怎样算？
CO₂RR 的首个质子化常比较 CO₂ + H⁺ + e⁻ → *COOH 与 CO₂ + H⁺ + e⁻ → *OCHO。CHE 近似下，H⁺ + e⁻ 可由 1/2H₂ 表示；自由能差还需加入零点能、熵修正、电位项和必要的溶剂修正。对于极性较强的 *COOH 和 *OCHO，隐式溶剂或显式水层可能改变相对稳定性。
自由能修正可用于同一电位下比较两个路径。若 ΔG*COOH 低于 ΔG*OCHO，只能说明 首个质子化更偏向 CO 路径入口；后续 *CO 脱附、*CHO 形成或 HCOO* 转化仍需继续判断。首步自由能偏好不能单独等同于最终产物选择性。
ΔG*COOH/*OCHO = ΔE + ΔZPE – TΔS – eU + ΔGsolv 式(2)
电子结构怎样支撑路径判断？
自由能只能告诉我们哪个构型在给定参照下更稳定，不能直接说明稳定性来自 C 端成键、O 端配位还是双金属协同。对 *COOH/*OCHO 比较来说，PDOS 应重点看吸附物 C/O 2p 态与金属 d 态或 Sn/Zn p 态之间的重叠，差分电荷密度则用于判断电子密度是否集中在 C-O、O-M 或 C-M 键附近。
若桥位 *OCHO 的 O-M 相互作用强，但 C 端活化不足，路径可能偏向甲酸盐相关物种；若 *COOH 的 C-M 相互作用适中且 *CO 不被过度稳定，CO 路径更容易延续。图4 和图5 中的 DFT 结果可以作为读谱参照，但正文判断应落回 p-d 轨道耦合、中间体键合方式和自由能入口 三者是否一致。
图4 Zn₁Sn₁/SNC 上 CO₂RR 中间体构型、自由能路径和电子结构分析。图中展示了 *HCOO、*COOH 等路径相关物种的计算证据。DOI: 10.1038/s41467-025-57573-4。
图5 双原子位点中 p-d 轨道耦合与 CO₂RR 机制关系。该图可用于理解金属组合、掺杂环境和中间体吸附方式之间的关联。DOI: 10.1038/s41467-025-57573-4。
桥位构型能直接和端位比吗？
桥位和端位可以比较，但前提是它们来自同一个双原子基底、同一个电荷和磁态口径、同一个自由能修正方案。若桥位构型引发双金属位点重构，而端位构型没有重构，能量差就同时包含吸附方式和位点重构两部分，需要在正文中分开解释。
有些构型在优化后会改变原本的分类。初始设置为桥位 *COOH 时，如果优化后 C 端只连接 M₁，另一个 O 原子远离 M₂，它就不应继续命名为桥位 *COOH；*OCHO 优化后若 H 转移到载体 N 或邻近水分子，对应的是质子转移事件，需要重新定义端点和路径。
结果矩阵应留下哪些量？
最终结果建议保留吸附物名称、初始构型、优化后构型、关键键长、ΔG、局域磁矩、Bader 电荷、PDOS 特征和差分电荷密度描述。若包含 HER 竞争，还应加入 H* 的 ΔGH*。这些量能帮助区分“*COOH 更稳定”“*OCHO 更稳定”和“二者能量接近但后续步骤不同”三种情形。
对氮掺杂石墨烯负载双原子模型来说，VASP 任务可按基底优化、CO₂ 预吸附、*COOH/*OCHO 构型矩阵、H* 竞争、自由能修正和电子结构分析推进。双原子 CO₂RR 的可靠比较来自同一模型族内的构型矩阵，而不是来自两个孤立优化结果。
要点：双原子位点需要同时定义金属组合、配位环境、金属间距和载体掺杂位置。
要点：*COOH 与 *OCHO 要按端位、桥位和偏单金属构型建立矩阵，优化后再重新命名。
要点：首个质子化自由能只能说明路径入口偏好，最终选择性还要看后续脱附、转化和 HER 竞争。
要点：PDOS、差分电荷和 Bader 电荷应服务于吸附方式解释，而不是孤立展示。
