为什么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计算服务优势
维度 | 测试Go |
路径覆盖 | 不只算CO₂→*COOH→*CO→CO这一条标准路径，自动生成6种产物12条分支路径，包括C-C偶联和HER竞争 |
溶剂处理 | 显式水层+隐式溶剂双验证，溶剂化修正直接嵌入计算流程 |
过渡态搜索 | CINEB+二聚体验证+初猜路径自动生成，收敛率超过85% |
电位依赖 | CHE模型快速修正+外加电场显式方法，两个都支持 |
数据交付 | 完整分析报告：参数测试、收敛性证明、与实验对比建议，图表可直接用于论文 |
复杂体系 | C-C偶联、单原子合金、双电层模型这些难算的东西，我们有专门的处理方案 |
CO₂RR计算专项支持计划
免费可行性评估：你提供催化剂体系和实验现象（产物分布、关键电位），我们24小时内给计算方案和预期能解决的科学问题。
合作送顶刊图表模板：自由能台阶图、NEB路径图、PDOS组合图等，按JACS/Angew格式出图，省排版时间。
疑难体系8折：涉及C-C偶联、单原子合金、pH>7或酸性介质中的CO₂RR，首单计算资源费8折。
发表奖励：如计算部分协助发Nature子刊、JACS、Angew，额外发放科研奖励金。
如果你有计算方面的难题
欢迎联系测试GO
从电子结构到宏观性能
我们可提供一站式多尺度模拟解决方案