手机版

电催化模拟计算

发布时间：2026-09-02 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

析氧反应（OER）是电解水制氢、金属-空气电池与可再生燃料合成中的瓶颈半反应。OER 涉及四电子四质子转移，动力学缓慢、过电位高，严重制约能源转换效率。贵金属氧化物 RuO₂、IrO₂ 虽活性优异但成本高昂，因此通过理论计算理解 OER 的决速步骤、量化过电位、指导催化剂设计，是电催化研究的核心任务。

本项目的科学问题是：RuO₂(110) 表面 OER 四步反应路径的自由能变化如何？理论过电位是多少？不同催化剂（RuO₂、IrO₂、NiFe-LDH、Pt）的过电位差异如何？计算结果能否解释实验上的活性趋势？

二、计算方法与理论背景

2.1 OER 四步机理

在酸性条件下，OER 经典四步机理为：*+H₂O → *OH + H⁺ + e⁻；*OH → *O + H⁺ + e⁻；*O + H₂O → *OOH + H⁺ + e⁻；*OOH → * + O₂ + H⁺ + e⁻。每一步的吉布斯自由能变化 ΔG_i 由 DFT 能量、零点能、熵变以及电位相关项组成：ΔG = ΔE + ΔZPE − TΔS + eU。在 U = 0 V 时各步 ΔG 之和为 4.92 eV（对应 1.23 V × 4e⁻）。

2.2 过电位定义

理论过电位 η = max(ΔG_i) / e − 1.23 V，其中 max(ΔG_i) 为 U = 0 时各步自由能的最大值，即热力学决速步。理想催化剂各步 ΔG 均为 1.23 eV，η = 0；实际催化剂因中间体结合能不成比例，η 通常为 0.2–0.5 V。OER 的火山型活性曲线正是由 *OH、*O、*OOH 之间的标度关系决定。

2.3 计算氢电极（CHE）模型

计算氢电极模型将 H⁺ + e⁻ 的化学势与 1/2H₂ 在标准电极电位下的化学势关联，从而无需显式模拟溶液中的质子，即可得到每步自由能，是电催化 DFT 计算的标准方法。

三、计算设置与模型构建

软件与方法：VASP，PAW 赝势，PBE+U（Ru 4d、Ir 5d，U=2.5–3 eV）+ 隐式溶剂（溶剂化修正）。
表面模型：RuO₂(110) 四层原子层板，2×1 表面超胞，真空层 15 Å，底部固定。
k 点与截断：4×4×1 网格，平面波截断 500 eV。
中间体构型：*OH、*O、*OOH 分别在表面 Ru 位优化，考虑氢键网络。
修正：零点能与熵变用谐波近似计算；对 *OH/*OOH 施加经验性溶剂化修正（约 0.1–0.2 eV）；比较 RuO₂、IrO₂、NiFe-LDH 与 Pt(111) 的过电位。

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

图 4 的左面板（a）给出 RuO₂(110) 表面 OER 在 U = 0 V 与 U = 1.23 V 下的四步自由能图。U = 0 时各步自由能台阶为：*OH 1.23 eV、*O 累计 2.66 eV、*OOH 累计 4.19 eV、释放 O₂ 累计 4.92 eV，对应每步 ΔG 为 1.23、1.43、1.53、0.73 eV。最大步为第三步（*O → *OOH），ΔG = 1.53 eV，因此理论过电位 η = 1.53 − 1.23 = 0.30 V。在 U = 1.23 V 时各步均变为下坡（红色曲线），说明在该电位下反应在热力学上自发进行，但动力学上仍需克服活化能。图中标注的 η = 0.30 V 对应 OER 决速步为 *OOH 的生成，这与 RuO₂ 表面 *O 结合过强、*OOH 结合偏弱导致的标度关系瓶颈一致。

图 4 的右面板（b）给出四种催化剂的过电位对比：RuO₂ 为 0.30 V、IrO₂ 为 0.34 V、NiFe-LDH 为 0.26 V、Pt(111) 为 0.62 V。Pt(111) 过电位最高，与其对含氧中间体结合过弱、*OH 生成步成为决速步的经典结论一致；NiFe-LDH 在碱性条件下活性最高，过电位仅 0.26 V；RuO₂ 与 IrO₂ 的过电位接近，均属于高效的 OER 催化剂。该活性趋势与实验上”RuO₂ ≈ IrO₂ < NiFe-LDH ≪ Pt”的观测一致，验证了计算方法的可靠性。

两面板结合，左面板揭示了 OER 决速步与过电位的定量来源，右面板提供了不同催化剂的横向对比，共同完成电催化模拟计算的核心分析，为催化剂筛选提供了理论依据。