CO2RR计算：电催化还原路径与选择性预测的量化分析

CO2RR计算在电催化CO₂还原反应研究中是核心手段，但多步反应路径复杂、电位依赖建模困难、产物选择性预测偏差大，如何建立可靠计算模型输出反应动力学数据指导催化剂设计，是电催化计算领域的关键技术挑战。

项目背景是一个铜基电催化剂的CO₂还原反应机理研究。客户合成了Cu-Zn合金催化剂，实验测得CO₂还原的主要产物是CH₄和C₂H₄，但选择性随电位变化，需要用CO2RR计算解析选择性变化的微观机理。CO2RR计算在电催化领域是研究反应路径和选择性的标准手段，但电位效应的建模和质子耦合电子转移（PCET）的处理是两个核心技术难点。

催化剂表面建模是第一步。Cu-Zn合金的(211)面是阶梯面，催化活性高于平坦(111)面。构建了4层slab模型，Cu和Zn原子按1:1比例交替排列在阶梯位点。底层2层固定，顶层2层弛豫。真空层15Å。CO2RR计算中阶梯面的台阶位是CO₂吸附的活性位点——平坦面上CO₂吸附能太弱（-0.18 eV），阶梯面上增强到-0.45 eV，因为台阶位的配位不饱和Ni原子提供了更强的电子注入能力。

电位建模是CO2RR计算中最具争议的技术环节。电催化反应在施加电位条件下发生，但DFT计算是零电位基态计算。用了计算氢电极（CHE）方法来引入电位效应：在标准条件下，化学势μ(H⁺+e⁻)=½μ(H₂)，外加电位U时每个电子的能量修正为-eU。这意味着在不同电位U下，涉及n个电子转移的反应步骤自由能修正-n·eU。CHE方法的优点是不需要显式建模电解质和电极电位，但缺点是无法描述电解质中的离子效应和双电层结构。CO2RR计算中CHE方法是目前最实用的近似方案，更精确的恒电位方法（如Grand Canonical DFT）还在发展中。

反应路径包含8个关键步骤。CO₂吸附后依次经历：CO₂→COOH→CO→CHO→CH₂O→CH₃O→CH₄（甲烷路径）或CO₂→CO→CO-CO→OCCO→CH₃CH₂→C₂H₄（乙烯路径）。两条路径在CO中间体处分叉——CO是甲烷路径和乙烯路径的共同起点，选择哪条路径取决于CO在表面的覆盖度和二聚化概率。CO2RR计算中这个分叉点是解释选择性的关键。

活化能垒计算揭示了电位对选择性的影响。在U=0 V（零电位）时：甲烷路径第一步CO→CHO活化能0.78 eV，乙烯路径第一步CO二聚化活化能0.92 eV——甲烷路径能垒更低，CO倾向于继续加氢生成CH₄。但在U=-0.8 V时：甲烷路径活化能修正后为0.78-0.8=接近0 eV（几乎无垒），乙烯路径0.92-1.6=也大幅降低但需两步电子转移，甲烷路径在低电位下更有利。然而在U=-1.2 V时，CO覆盖度升高导致CO二聚化概率增大，乙烯路径的实际选择性反而上升。CO2RR计算需要同时考虑热力学能垒和动力学覆盖度效应才能准确预测选择性。

产物选择性预测的量化方法用了微动力学模型。将CO2RR计算得到的各步骤活化能和吸附能输入微动力学方程，用稳态近似求解各中间体覆盖度和各产物生成速率。在U=-0.6 V时甲烷选择性72%、乙烯28%；U=-1.0 V时甲烷45%、乙烯55%；U=-1.2 V时甲烷30%、乙烯70%。实验数据：U=-0.6 V时甲烷68%、乙烯32%；U=-1.0 V时甲烷42%、乙烯58%；U=-1.2 V时甲烷25%、乙烯75%。CO2RR计算预测的选择性趋势与实验完全吻合，绝对值偏差5-8个百分点——这个精度在电催化领域已经属于高质量预测。

电荷分析解释了Cu-Zn合金优于纯Cu的原因。Zn的d电子比Cu更局域化，在阶梯位上Zn提供了更强的电子注入能力——CO₂在Cu-Zn阶梯位获得0.42个电子（纯Cu阶梯位只有0.28个），更多的电荷注入意味着C-O键被更充分活化，降低了后续加氢步骤的能垒。CO2RR计算通过Bader电荷和d带中心分析把”合金为什么更好”从经验判断上升到电子结构层面的定量解释。

计算资源方面，8步反应路径×5个NEB image=40个并行计算，加上3个电位点的自由能修正，总计算时间约120小时（32核）。CO2RR计算的多电位分析虽然增加了计算量，但电位依赖性数据对指导实验优化电极条件具有直接参考价值。

https://www.keyanxueshu.com/