CO2电解理论计算：自由能图怎么画、过电位怎么算

CO₂电还原（CO2RR）这个领域，实验文章数以万计，催化剂筛选可以说是寸土必争。Cu是唯一能在可观法拉第效率下把CO₂还原到C₂+产物的金属，但问题是Cu的选择性太差了——从CO、甲酸、甲烷到乙烯、乙醇，产物分布像撒胡椒面。

理论计算要回答的核心问题是：给一个催化剂表面，它能优先走哪条反应路径？对应的起始电位是多少？这两个问题的答案都藏在一张图上——从CO₂到终产物的自由能阶梯图。

CHE模型的物理假设

计算氢电极（Computational Hydrogen Electrode, CHE）模型是Nørskov课题组在2004年提出的，已经成为电催化理论计算的行业标准。核心简化是：质子-电子对（H⁺ + e⁻）的化学势在0V vs RHE条件下等价于½H₂的化学势：

μ(H++e−)=12μ(H2)−eUμ(H++e−)=21​μ(H2​)−eU

U是外加电势。这个假设把电化学环境中的、需要做隐式溶剂才能描述的质子迁移过程，转化成了一个可以纯DFT处理的自由能问题。代价是把溶剂效应、双电层效应、阳离子效应全压扁了——但作为筛选工具，CHE模型的精度已经足够了。

Cu(111)上的关键中间体

Cu(111)表面的CO₂还原路径网络相当复杂。为简化讨论，重点关注两条竞争路径：CO路径（CO₂→COOH*→CO*→CHO*→CH₂O*→…→CH₄）和甲酸路径（CO₂→HCOO*→HCOOH）。

每个中间体的吸附自由能需要算三样东西：DFT总能、零点振动能修正（ZPE）、振动熵修正（TS）。ZPE和TS都来自DFT做完后的频率计算。CO₂(g)的热力学数据从NIST标准表中取——气相的ZPE和熵不能用VASP算的，因为气相平移/转动自由度在DFT周期性计算中被截断了。

实际跑下来，COOH在Cu(111)上的吸附自由能是0.65 eV（相对于CO₂(g) + ½H₂），而HCOO是0.72 eV。0.07 eV的差距并不大——不到3k_BT——在数值精度范围内，但趋势是对的：COOH比HCOO略稳定。这意味着Cu(111)走CO路径的起始电位稍低于甲酸路径，与实验上Cu以CO为主要产物的现象定性一致。

从自由能图判断过电位

自由能图绘制完成后，找到决定电位步骤（PDS，Potential-Determining Step）——也就是具有最大正自由能变化的那一步。在Cu(111)上，对于CO路径，PDS通常是CO₂→COOH这一步（ΔG=0.65 eV），或者CO→CHO*（ΔG约0.7 eV，取决于CO*覆盖度）。

理论起始电位 = -ΔG_PDS / e = -0.7 V vs RHE。这意味着从热力学角度看，只有当阴极电位比RHE负0.7V以上时，CO₂→甲烷的整条反应路径才是自由能下坡的。

实验上Cu电极的CO₂还原到甲烷通常需要-0.8到-0.9 V vs RHE才有可观的电流密度。CHE模型给出的-0.7V与实验的起始电位相当接近，说明PDS的判断是正确的。中间的0.1-0.2V偏差可能来自动力学过电位（Tafel斜率效应）和溶剂效应。

CO vs C₂H₄的选择性分叉

Cu(111)上另一个关键问题是C₁和C₂路径的分叉。CO是共有的中间体——如果CO脱附，得到的就是CO产物；如果CO继续加氢到CHO，走的是C₁烃类路径；如果两个CO*偶联（CO-CO coupling），则进入C₂路径。

CO-CO偶联在低覆盖度下是不利的——两个CO之间的距离太大，无法有效耦合。高覆盖度（>0.25 ML）下，CO之间的排斥使它们在表面上聚集成岛状，偶联能垒降低到0.5 eV左右。这个覆盖度效应解释了为什么Cu需要在较负的电位下才有明显的乙烯和乙醇产物——需要先把表面CO*复盖度推上去。

这算出来有什么用

用CHE模型做CO₂电解理论计算，不能给你一个确切的”最佳催化剂”。它的价值在于：把几十种可能的中间体和反应路径在自由能图上排列出来，快速筛掉热力学不可行的路径，聚焦在过电位最小的候选路径上。一个10步的自由能计算也许需要几天机时，但做出来了，实验组就知道该把表征资源集中在哪几个反应中间体上。

Cu(111)只是一个起点。把这个框架搬到Cu的高指数面、Cu基合金、单原子催化剂上，就能做催化剂的高通量理论筛选——这才是DFT在电催化领域的真正力量所在。