吉布斯自由能计算

一、项目背景与科学问题

电催化二氧化碳还原反应（CO₂RR）是将温室气体 CO₂ 转化为高附加值燃料与化学品的关键技术路径，其产物包括 CO、甲酸、甲醇及多碳产物等。反应的热力学驱动力与动力学瓶颈可由各基元步骤的吉布斯自由能变化（ΔG）描述。计算氢电极（CHE）模型由 Nørskov 课题组提出，将质子-电子对（H⁺ + e⁻）的化学势与标准氢电极相关联，使得每个基元步骤的自由能变化可仅由 DFT 总能量与零点能/熵修正计算得到，极大降低了电化学模拟的复杂度。

本项目选取 Au(111) 表面作为模型催化剂，沿 CO₂ → *COOH → *CO → CO(g) 路径计算各基元步骤的 ΔG，确定反应的电位决定步骤（PDS）与极限电位，为理解 Au 电极选择性产 CO 的机理提供热力学依据。

二、计算方法与理论背景

2.1 计算氢电极（CHE）模型

在 CHE 模型中，对于步骤 *+CO₂ + H⁺ + e⁻ → *COOH，质子-电子对自由能通过公式 μ(H⁺+e⁻) = ½μ(H₂) − eU 计算，其中 U 为外加电位。因此各步骤的 ΔG 在 U = 0 V vs RHE 时仅取决于吸附物种的 DFT 能量差，并可通过 −eU 项线性移动到任意电位。

2.2 自由能修正

ΔG = ΔE_DFT + ΔZPE − TΔS，其中 ZPE 为零点能修正，由吸附物种振动频率计算；S 为熵，气相物种取标准熵值，吸附态取振动熵近似。本计算温度取 298.15 K。

2.3 吸附能计算

吸附能定义为 E_ads = E(slab+ads) − E(slab) − E(gas)，采用 4×4 Au(111) 超胞、4 层金属 slab，固定底层 2 层。k 点与截断能设置经收敛测试确定，保证吸附能误差小于 0.05 eV。

三、计算设置与模型构建

• 表面：Au(111) 4×4 超胞，真空层 ≥ 15 Å。
• 泛函：GGA-PBE，PAW 赝势，截断能 400 eV。
• 吸附位点：顶位（top）、桥位（bridge）、fcc 与 hcp 空位均作考察，取最稳定构型。
• 电位扫描：U = 0 V 与 U = −0.5 V vs RHE。

四、结果分析与案例图解读

图 3 的左面板为 U = 0 V（灰色圆点线）与 U = −0.5 V（红色方块线）两种电位下的自由能阶梯图，横坐标为反应坐标（CO₂ → *COOH → *CO → CO），纵坐标为 ΔG。在 U = 0 V 时，CO₂ 活化生成 *COOH 的自由能上升最大，ΔG(*COOH) = 0.58 eV，是电位决定步骤；随后 *COOH 还原为 *CO 的 ΔG = 0.06 eV，接近热中性；最后 *CO 脱附为 CO(g) 的 ΔG = −0.30 eV，为放热步骤。整个路径的 ΔG_max = 0.58 eV，对应极限电位 U_L = −0.58 V，与实验报道 Au 电极上 CO 生成起始电位（约 −0.5 V）相符。

当电位移动到 U = −0.5 V 时，含质子化步骤的自由能被拉低 0.5 eV，此时 *COOH 生成步骤的 ΔG 降至 0.08 eV，接近热中性，反应速率显著提升。该结果解释了 Au 电极在中等过电位下即可高效产 CO 的实验现象。

图 3 的右面板对比了 Ag、Au、Cu、Zn、Pd 五种金属上 *COOH 的自由能（文献值）。Au（0.58 eV）与 Ag（0.45 eV）的 ΔG(*COOH) 相对较低，有利于 CO 路径；而 Cu 的 ΔG(*COOH) 较高（0.68 eV），其限速步骤特征不同，这与 Cu 上多碳产物选择性相关联。该对比图可直接支撑论文中”不同金属催化活性比较”的论述。

五、结论与展望

本项目采用 DFT + CHE 模型计算了 Au(111) 上 CO₂RR 关键路径的吉布斯自由能分布，识别出 *COOH 生成为电位决定步骤（ΔG = 0.58 eV），极限电位约 −0.58 V，与实验趋势一致。案例图以自由能阶梯图直观展示电位对反应路径的调制，辅以多金属对比图突出 Au 的选择性优势。

后续可拓展：(1) 引入显式溶剂与电场效应，修正吸附自由能；(2) 计算 *CO 覆盖度效应与竞争析氢路径，完善机理；(3) 结合微动力学模型，定量预测电流密度与法拉第效率。