本文将对近期发表在学术期刊*Journal of Catalysis*上的一项研究工作进行全面介绍。该研究由Shi-Qin Xiang, Shu-Ting Gao, Jun-Lin Shi, Wei Zhang, 以及 Liu-Bin Zhao等人完成，作者单位包括西南大学化学化工学院化学系与中国科学院重庆绿色智能技术研究院。论文发表于*Journal of Catalysis*第393卷（2021年），已于2020年11月25日在线发布。这是一项关于铜电极上电催化二氧化碳还原反应（CO2RR） 的原创性理论计算研究。
一、 研究的学术背景与目标
二氧化碳的电化学还原是一项极具前景的技术，有望将可再生的间歇性电能转化为可存储的燃料与高价值化学品，对于实现碳中和目标具有战略意义。在众多电催化剂中，铜（Cu）基材料因其独特的性质而备受关注：它是目前已知的唯一能够将二氧化碳高效还原为多种碳氢化合物（如甲烷、乙烯）和含氧化合物（如醇类）的金属，尽管这些多电子产物的法拉第效率较低且需要较高过电位。而对于双电子还原产物，如一氧化碳（CO）和甲酸（HCOOH），则可以在其他催化剂（如金、银）上以较低过电位和高选择性获得。
然而，理论模拟与实验结果之间尚未建立起理想的一致性。因此，深入探究二氧化碳电还原的反应机理，对于理解和调控催化活性与产物选择性至关重要。本研究的核心目标，正是针对CO2RR机理中两个关键且长期悬而未决的问题展开深入研究： 1. 竞争的反应路径机理：二氧化碳还原是一个涉及多电子与多质子的质子耦合电子转移（Proton-Coupled Electron Transfer, PCET） 过程。在此过程中，质子和电子的转移方式存在两种基本模式：协同质子电子转移（Concerted Proton-Electron Transfer, CPET），即质子和电子作为氢原子（H）同时一步转移；以及分步质子电子转移（Sequential Proton-Electron Transfer, SPET），即质子转移和电子转移先后分步进行。过去的研究通常简单采用CPET模型，但越来越多的证据表明，在实际的电极/电解质界面上，SPET路径也可能存在。两种路径的竞争关系以及其如何受到实验条件（如溶液pH值、电极电位）的影响尚不清楚。 2. 热力学控制与动力学控制的竞争：电催化反应的一个独特之处在于，施加的电极电位可以同时调控反应的驱动力（热力学）和反应速率（动力学）。通常，热力学控制的产物具有更早（更正）的起始电位，而动力学控制的产物则具有更低的反应能垒。在复杂的反应网络中，热力学有利的路径与动力学有利的路径可能并不一致，这直接决定了最终产物的分布。
为了系统揭示以上问题，研究团队旨在开发一个综合的微动力学模型（Micro-kinetic Model），将电化学速率理论与第一性原理密度泛函理论（Density Functional Theory, DFT） 模拟相结合，定量研究溶液pH值和电极电位如何影响铜表面CO2还原为CO和HCOOH的反应速率、反应路径和产物选择性，从而为理性设计电催化体系提供理论指导。
二、 研究的详细工作流程
本研究是一项理论计算工作，其工作流程高度依赖计算机模拟和理论建模，而非传统的实验操作。主要包含以下几个相互关联的阶段：
第一阶段：计算模型的建立与电子结构计算 本研究的计算对象是基于铜晶体的（111）晶面建模的表面模型。研究采用了周期性平板模型，使用4层原子厚的铜板，其中底部两层原子在几何优化时被固定以模拟体相。为了处理电化学反应环境，研究中引入了VASPsol隐式溶剂模型来模拟电解质溶液的影响，该模型基于线性化的Poisson-Boltzmann方程，考虑了静电力、空化能和弥散力等溶质-溶剂相互作用。所有的第一性原理计算均通过Vienna Ab Initio Simulation Package (VASP) 软件包完成，使用Projector Augmented Wave (PAW) 方法描述离子芯，交换关联泛函则采用广义梯度近似框架下的Perdew-Burke-Ernzerhof (PBE) 泛函。平面波截断能设为400 eV以保证能量收敛。
第二阶段：热力学数据获取与反应网络构建 这是研究的基础步骤，旨在获取所有关键中间体和反应步骤的吉布斯自由能。首先，对反应网络中涉及的所有关键物种（如反应物CO2，关键中间体CO2⁻, COOH, HCOO*，以及产物CO和HCOOH*）在铜（111）面上的最稳定吸附构型进行几何优化和频率计算。根据这些数据，可以计算每一步基元反应（包括CPET、SPET等）在不同电极电位和pH值下的反应自由能变。此阶段的关键是对电化学系统的自由能进行巨正则（Grand Canonical） 修正，以准确反映外电势和pH对体系能量的影响。基于这些热力学数据，研究者绘制了详细的反应网络图（对应原文Fig. 1），清晰展示了从CO2到CO和HCOOH的所有可能的平行与串联反应路径，包括CPET、电子转移（ET）、质子转移（PT）等各类基元步骤。
第三阶段：发展并应用微动力学模型 这是本研究方法学的核心创新部分。研究人员发展了一套结合电化学速率理论（主要是基于Marcus-Hush-Chidsey理论的扩展）和第一性原理计算的微动力学模型，用于计算每个基元反应步骤的电位和pH依赖的速率常数。 该模型的核心公式用于计算反应速率常数 ( k )： [ k = \frac{|v|^2}{\hbar} \sqrt{\frac{\pi}{kB T \kappa}} \left[ \exp\left( -\frac{E{a,f}}{kB T} \right) - \exp\left( -\frac{E{a,b}}{kB T} \right) \right] ] 其中，( v ) 是电子耦合强度，( \kappa ) 是重组能（Reorganization Energy），( E{a,f} ) 和 ( E_{a,b} ) 分别是正、逆向反应的活化能。 活化能的计算是本模型的关键。它并非固定值，而是与驱动自由能变（由DFT计算和施加的电极电位共同决定） 和重组能呈二次方关系。这种处理方式使得对称因子（活化能对电位的导数）随电位变化，与传统的Butler-Volmer方程中恒定的对称因子（通常为0.5）有本质区别。这能更真实地描述电化学反应动力学。模型中，所有必要的热力学和动力学参数（如重组能、活化能）均通过DFT计算获得，并整合到上述速率理论中。
第四阶段：系统分析计算与结果可视化 研究者利用所开发的模型，系统地计算了在不同pH值（从酸性到碱性）和不同电极电位（从起始电位到高过电位）下，所有竞争反应路径的速率常数。这些计算旨在回答核心问题：在给定的实验条件下，哪一种反应机理（CPET vs. SPET）占主导？哪一条反应路径是速率决定步骤？最终的产物选择性（CO vs. HCOOH）如何？研究通过绘制电位-速率常数（极化曲线）图（对应原文Fig. 4）和产物生成速率图（对应原文Fig. 5），直观地展示了这些复杂的竞争关系随外部条件的演化规律。此外，还通过对比在特定pH和电位下的势能曲线（Potential Energy Curves）（对应原文Fig. 2和Fig. 3），从活化能和反应热力学的角度详细解析了竞争机制的转变过程。
三、 研究的主要结果详述
研究结果深刻揭示了CO2在铜（111）面上还原为CO和HCOOH的复杂机制，主要发现如下：
1. 关键中间体的pKa值决定反应路径选择： 计算得到的两个关键中间体COOH*和HCOO*的pKa值（解离常数） 分别为2.01和4.66。这个数值是理解反应路径随pH变化的钥匙。在酸性（pH < pKa）条件下，质子供体充足，CO2分子倾向于通过协同质子电子转移（CPET）一步形成质子化的中间体COOH*或HCOO*。而在碱性（pH > pKa）条件下，质子浓度极低，CO2分子会先通过一个纯电子转移（ET）步骤形成CO2*⁻自由基阴离子，然后再发生质子化。因此，pKa值成为区分CPET机理和SPET机理的天然分界点。
2. 从热力学控制的CPET到动力学控制的SPET的机理转变： 通过分析不同pH值下的势能曲线，研究清晰地展示了一种“控制权”的转移。 * 低pH（酸性环境）： 在pH=0时，CO2还原为COOH*的热力学起始电位（-0.59 V vs. SHE）要远早于（更正于）其还原为CO2*⁻的电位（-0.71 V vs. SHE）。尽管CO2*⁻路径的活化能更低（动力学更有利），但热力学上的巨大优势使得COOH*路径在低过电位下优先发生。此时，反应遵循热力学控制的CPET机理。 * 高pH（碱性环境）： 当pH大于COOH*的pKa时，由于缺少质子，COOH*路径的热力学起始电位会向负方向移动（根据能斯特方程），变得比CO2*⁻路径更负。此时，虽然两者热力学上可能接近或逆转，但形成CO2*⁻的ET步骤具有更低的活化能。因此，反应路径由动力学优势主导，即优先通过SPET机理进行（先ET生成CO2*⁻，再PT生成COOH⁻*）。 HCOOH的生成路径也遵循完全相同的规律，其转变点在HCOO*的pKa（4.66）附近。这一发现澄清了关于PCET机理的争议，表明反应机制并非固定，而是强烈依赖于溶液的pH环境。
3. 过电位与pH共同调控最终产物选择性： 微动力学模型计算得出的产物生成速率图（对应原文Fig. 5）展示了CO和HCOOH选择性随电位和pH演变的复杂图像。 * HCOOH（甲酸）是热力学控制产物： 由于其生成的CPET路径（C4）具有最正（即所需过电位最低）的起始电位，因此在低过电位区域，无论pH如何，HCOOH通常是首先生成的主要产物。这源于HCOO*中间体与铜表面的结合能更强。 * CO（一氧化碳）是动力学控制产物： 随着过电位的增大（电极电位更负），其他竞争路径的速率开始提升。研究发现，即使在碱性条件下，CO通过SPET路径（E1→C3→P3）具有显著的动力学优势。因此，在中等至高过电位下，CO会逐渐取代HCOOH成为主要产物。 * 碱性环境利于提高HCOOH的法拉第效率： 在碱性条件下，水合氢离子浓度极低，极大地抑制了竞争反应——析氢反应（Hydrogen Evolution Reaction, HER）。同时，模型计算表明，在特定的“转折电位”下，HCOOH的生成速率随pH升高而增强。这意味着，在碱性环境中，通过精确控制电位，有可能实现高选择性的HCOOH生成。这与实验观察中碱性条件下某些催化剂（如银）的产物选择性从CO转向甲酸的现象一致。
4. 关于反应速控步与Tafel斜率的认识： 通过比较各步骤的速率常数，研究发现，无论是CO路径还是HCOOH路径，第一个电子转移步骤（E1, C1或C4）的速率常数远小于后续的第二个电子/质子转移步骤。因此，整个CO2RR过程的速控步（Rate-Determining Step, RDS） 是第一个电子转移阶段，具体是哪一步则取决于pH和电位。 此外，基于该模型的动力学分析还预测了反应的Tafel斜率。由于活化能与电位呈非线性关系（二次方依赖），Tafel斜率不是一个固定值，而是在102至112 mV/dec之间变化（在-0.8 V至-1.2 V vs. SHE电位区间内），这与许多实验观测到的Tafel行为相符合。
四、 研究的结论与价值
本研究的核心结论是：通过结合第一性原理计算与电化学微动力学模型，成功揭示了溶液pH值和电极电位对铜电极上CO2还原为CO和HCOOH的催化活性和选择性的精细调控机制。 * 科学价值： 1. 统一了对PCET机理的理解：明确了在CO2RR中，CPET和SPET机理并非互斥，而是相互竞争，其主导权取决于反应环境的pH值。低pH下为热力学控制的CPET，高pH下则转变为动力学控制的SPET。 2. 建立了热力学与动力学的关联框架：成功地将宏观的产物选择性（CO vs. HCOOH）与微观的热力学（起始电位）和动力学（活化能、重组能）控制因素联系起来，解释了为何HCOOH在低过电位占优（热力学控制），而CO在高过电位占优（动力学控制）。 3. 发展了一个普适性的理论工具：所开发的微动力学模型不仅限于CO2RR，其结合Marcus理论与DFT计算的方法，为研究其他复杂的电催化反应（如析氧反应（OER）、氧还原反应（ORR）、氮还原反应（NRR））的电位、pH依赖动力学提供了强大的理论框架。 * 应用价值： 研究结果为通过理性调控实验条件来优化催化剂性能提供了清晰的指导。例如，要高效生产甲酸，应倾向于使用碱性电解质并在较低过电位下运行；而要获得高产率的一氧化碳，则需要施加较高的过电位。同时，通过调控电极电位和溶液pH，可以有效抑制副反应（如HER），提高目标产物的法拉第效率。
五、 研究的亮点
六、 其他有价值的补充信息
本研究的支持信息中还提供了大量详细数据，如各中间体的吸附构型、所有基元反应在平衡电位下的重组能和活化能表格、不同pH下HER竞争反应的起始电位分析等。这些数据为读者深入理解模型的构建和结果的推导提供了充分的细节，增强了研究的透明度和可重复性。作者团队来自中国，研究得到了国家自然科学基金、重庆市人才计划、中央高校基本科研业务费等多个项目的支持，体现了中国科研工作者在能源催化基础理论研究领域的贡献。