第一性原理计算电化学

项目过程：RuO₂ 催化 OER 的过电位计算

客户拿到 RuO₂ 作为酸性电解水阳极催化剂，想知道它的理论过电位，和文献值对比一下，判断材料是否值得继续优化。我的任务是用第一性原理把四电子 OER 的自由能台阶算出来，找出决速步。

计算用 VASP，采用计算氢电极（CHE）模型。RuO₂ 取金红石结构，对 Ru 4d 电子加了 DFT+U 修正，U_eff=3.0 eV。平面波截断能 520 eV，K 点 4×4×6，真空层 15 Å。反应路径按 * → *OH → *O → *OOH → O₂(g) 四步处理，自由能校正包含 ZPE 和熵贡献，温度取 298 K。

电势分别取了 U=0 V、1.23 V 和 1.52 V 三个点。U=1.23 V 是热力学平衡电位，U=1.52 V 是根据自由能图读出的起始电位。

案例分析：自由能图直接锁定决速步

第一张案例图是三条自由能曲线。在 U=0 V 时，四步都是上坡，最大自由能变化出现在 *OOH→O₂ 这一步，约 1.5 eV；当电位升到 1.23 V 时，前三步基本被拉平，最后一步仍有约 0.29 eV 的上坡；到 1.52 V 时，所有台阶都变成下坡或平坡，整个反应自发进行。

第二张案例图把每一步的自由能变化单独画成柱状图。可以看到 *OOH→O₂ 的 ΔG 最大，是决速步（PDS）。从 U_eq=1.23 V 和 U_onset=1.52 V 得到过电位 η=U_onset−U_eq=0.29 V，这与 RuO₂ 作为优质 OER 催化剂的实验认知一致。

总结教训：CHE 模型的边界条件要交代清楚

第一，CHE 模型假设质子-电子对平衡于可逆氢电极，这个近似在强酸性介质里合理，但在强碱性或高 pH 体系里需要换成 RHE 标度并考虑 pH 校正。

第二，*OOH 的吸附能是误差最大的环节。不同泛函对这个中间体的描述差异可达 0.2–0.3 eV，做不同催化剂横向对比时必须用同一套泛函和参数。

第三，DFT+U 对 RuO₂ 是必需的。不加 U 时 Ru 3d 电子关联被低估，会导致 *O 吸附能过强，从而改变决速步判断。

项目经验：过电位不是越低越好，还要看稳定性

RuO₂ 的 η=0.29 V 很漂亮，但实际工业应用里更关心它在酸性 OER 条件下的溶解稳定性。这个项目算出来的热力学过电位只能回答“催化活性够不够好”，不能回答“能用多久”。如果客户要继续往下做，我建议把表面 Pourbaix 图也补上，把活性和稳定性一起评估。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：第一性原理计算电化学。