为什么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₆双原子催化剂的金属选择——高通量DFTML客户的问题想设计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偶联、单原子合金、pH7或酸性介质中的CO₂RR首单计算资源费8折。发表奖励如计算部分协助发Nature子刊、JACS、Angew额外发放科研奖励金。如果你有计算方面的难题欢迎联系测试GO从电子结构到宏观性能我们可提供一站式多尺度模拟解决方案