CO₂RR计算的“疑难解答手册”--80+课题组验证的实战方案分享

为什么CO₂RR这么需要计算？

CO₂电还原是碳中和方向里发文章最密集的催化领域之一，Angew、JACS、Nature Catalysis上几乎每期都有。但做实验的人经常卡在同一个问题上：产物分布测得出来，原因说不清楚。

►电化学测试能告诉你CO、HCOOH、CH₄各占多少，但解释不了"为什么这个催化剂偏偏走CO路径，不走甲酸路径"。

►合成了一堆单原子、合金、氧化物衍生材料，活性位点的真实结构和电子态还是猜的。

►审稿人基本必问"请补充DFT计算解释选择性"，但自己搭集群、算NEB、做溶剂化修正，门槛不低。

测试GO在这一方向已服务

80+课题组算过CO₂RR，

其中12篇发在Nature子刊、JACS或Angew上。

下面按实际遇到的问题来讲我们怎么做的、

算出来的东西能解决什么。

三个常见计算难题，我们怎么处理

1

选择性说不清——CO还是甲酸还是多碳？

现象：同一个催化剂，电位变了或pH调了一下，产物就从CO变成甲酸，甚至冒出乙烯。实验只能看到结果。

我们的做法：

◎把竞争路径全部算一遍。比如*COOH路径和*OCHO路径同时跑，C1路径和C-C偶联路径同时跑，给每一步的ΔG和极限电位（UL）差值列出来。

◎用*CO吸附能和*H吸附能画二维火山图，看催化剂落在"选择性地图"的哪个区域，天然偏向CO、甲酸还是析氢。

◎加恒电势修正（CHE模型+外加电场显式模型），算不同电位下的自由能变化，结果可以直接和LSV曲线、产物分布数据对标。

交付物：多路径自由能台阶对比图、极限电位差值表、选择性热力学分析——讨论部分可以直接用。

2

活性位点太多筛不过来——单原子/合金表面构型爆炸

现象：M-N-C单原子催化剂的配位环境可能是M-N₄、M-N₃C₁、M-N₂C₂……上百种构型，不可能一个一个合成验证；双原子催化剂更是组合爆炸。

我们的做法：

◎与客户共建了3000+种金属-配位环境的结构数据库，用d带中心、配位数、Bader电荷、CO吸附能等描述符加上高斯过程回归做快速筛选，一周内可以从200个候选构型缩到5个。

◎初筛结果再用DFT精算，迭代训练模型，降低计算成本。

◎案例：帮一个985团队筛Fe-Co双原子位点，从78种构型里预测出FeCo-N₆为CO₂→CO最优活性中心，后续合成验证法拉第效率96%，发在Adv. Mater.上。

交付物：构型-吸附能热图、ML特征重要性排序、推荐位点的PDOS/COHP分析。

3

反应路径里的隐藏高能垒找不到——尤其是C-C偶联

现象：C₂+产物（乙烯、乙醇）是CO₂RR的圣杯，但算起来很难：*CO二聚能垒高、中间体多、过渡态搜索容易失败；甲酸路径的*OCHO构型特殊，常规优化经常跑偏。

我们的做法：

◎C-C偶联（*CO+*CO→*OCCO）用Climbing Image NEB做过渡态搜索，配合初猜路径自动生成和二聚体方法验证，保证找到真实能垒而不是局部极小。

◎长路径反应（CH₄、C₂H₅OH）里质子-电子转移步骤经常有水分子参与，我们用VASPsol隐式溶剂+AIMD显式水分子快照双验证，消除气相模型的误差。

◎把DFT算出的各基元步骤能垒输入微观动力学模型（kMC或稳态近似），理论预测产物法拉第效率随电位变化的曲线，和实验极化曲线对标。

交付物：C-C偶联NEB能垒图（含虚频验证）、不同电位下理论产物分布图、pH影响模拟。

三层递进方案：从基础到顶刊

我们不按"算一个吸附能"这种碎片化方式收费，而是围绕论文设计了三层方案。

针对C₂+产物、单原子合金、金属氧化物界面这些更复杂的体系，我们提供显式-隐式混合溶剂模型+恒电势双电层模型，让计算结果更接近真实电化学环境。

服务客户案例

Cu₃N纳米晶催化CO₂→CH₄，审稿人质疑HER竞争

客户的问题：Cu₃N在-0.8V vs RHE下CH₄法拉第效率约55%，HER效率35%。审稿人问：CH₄路径能垒是否真的低于HER？请补DFT。

我们做了什么：

◎建Cu₃N(111)表面模型，算*CO→*CHO→*CH₂O→ *OCH₃→*CH₄完整路径。

◎同时算竞争反应*H→*H₂的Tafel/Heyrovsky步能垒。

◎发现Cu₃N表面*CO加氢为*CHO的能垒0.45eV，而*H结合能垒0.78eV，热力学上解释了CH₄优于HER。

◎补充*CHO的NEB过渡态虚频验证，消除审稿人对过渡态真实性的疑虑。

FeM-N₆双原子催化剂的金属选择——高通量DFT+ML

客户的问题：想设计Fe基双原子催化剂（FeM-N₆，M= Mn, Co, Ni, Cu, Zn等）做CO₂→CO，不知道哪种第二金属最好。

我们做了什么：

◎构建36种FeM-N₆构型（不同金属、不同配位环境）。

◎高通量DFT计算*COOH吸附能作为关键描述符。

◎训练随机森林模型，发现金属d电子数和电负性是主要描述符，预测FeCo-N₆和FeNi-N₆活性最高。

◎额外算*CO脱附能，确认产物CO容易释放、不会毒化位点。

Cu(111) vs Cu(100)乙烯选择性差异——

覆盖度和溶剂的隐藏效应

客户的问题：实验发现Cu(100)乙烯选择性显著高于Cu(111)，但常规DFT算*CO二聚能垒两者差不多，解释不了。

我们做了什么：

◎引入覆盖度效应：在高*CO覆盖度（θ>0.33）下，表面吸附物之间的排斥作用改变了二聚过渡态的几何构型。

◎加显式溶剂化模型（3层水分子+VASPsol），发现Cu(100)表面水层协助稳定*OCCOH中间体，C-C偶联能垒降低了0.18eV。

◎建微观动力学模型，预测Cu(100)在-0.9V时C₂产物选择性比Cu(111)高2.3倍，和实验吻合。

测试GO计算服务优势

CO₂RR计算专项支持计划

免费可行性评估：你提供催化剂体系和实验现象（产物分布、关键电位），我们24小时内给计算方案和预期能解决的科学问题。

合作送顶刊图表模板：自由能台阶图、NEB路径图、PDOS组合图等，按JACS/Angew格式出图，省排版时间。

疑难体系8折：涉及C-C偶联、单原子合金、pH>7或酸性介质中的CO₂RR，首单计算资源费8折。

发表奖励：如计算部分协助发Nature子刊、JACS、Angew，额外发放科研奖励金。

如果你有计算方面的难题

欢迎联系测试GO

从电子结构到宏观性能

我们可提供一站式多尺度模拟解决方案