电化学材料dft计算

电解水制氢、金属空气电池这些方向，阳极的析氧反应（OER）是最慢的一步，决定了整个器件的效率天花板。我这个项目客户在筛 OER 催化剂，拿了几组过渡金属氧化物/氢氧化物来找”哪一个过电位最低”。我做的 电化学材料dft计算 核心，就是把 OER 的四步自由能台阶图画出来，用计算直接给出过电位排序。

OER 在碱性条件下走吸附演化机制（AEM）：* + OH⁻ → *OH + e⁻ → *O + H₂O + e⁻ → *OOH + H₂O + e⁻ → * + O₂ + H₂O + e⁻。每一步都涉及一个质子-电子对转移，我用计算氢电极（CHE）模型把每一步的自由能算出来。每个吸附中间体（*OH、*O、*OOH）都在催化剂表面优化，自由能用 G = E_DFT + ZPE − TS 修正，其中 ZPE 和熵从振动频率拿。

判读的关键是”过电位”。在理想情况下，四步平均每个台阶应该是 1.23 eV（对应水的理论分解电压），这样过电位为 0。但实际催化剂因为 *OH、*O、*OOH 之间的标度关系（scaling relation），三个吸附能不是独立的，导致至少有一个台阶会高于 1.23 eV，多出来的部分就是过电位 η。我画了四步台阶图，找到那个”最高台阶”，η = 最高台阶 − 1.23 eV。我这个项目里催化剂 A 的最高台阶是 1.75 eV，η≈0.52 V；催化剂 B 只有 1.45 eV，η≈0.22 V——B 的活性明显更好，和客户的 LSV 起始电位趋势一致。

更进一步的洞察来自火山图。我把一系列候选催化剂的 *OH 吸附能（作为描述符）和过电位画成火山曲线：吸附太弱（左翼）时 *OOH 形成困难，吸附太强（右翼）时 *O 脱附困难，中间有个火山顶。客户最初最看好的那个高载量样品，其实落在右翼（吸附过强），解释了它过电位反而高的反常。我用这个图给客户一个清晰的选材判据：找 *OH 吸附能在火山顶附近的体系。

踩坑记录：第一，标度关系的近似。我用 *OOH 和 *OH 的吸附能差约 3.2 eV 这个经验标度来约束，但某些体系偏离，必须逐个算 *OOH 验证。第二，溶剂和电场。气相算的台阶和真实电化学界面差很大，我用了隐式溶剂 + 外加电场（或 CHE 的 pH 修正）来逼近。第三，Hubbard U。含 3d 金属的氧化物，O 的 p 带和金属 d 带强关联，不加 U 带隙和吸附能都错，U 值我从文献+线性响应法标定。第四，表面终止和覆盖度。同样的材料不同表面晶面 OER 活性差几倍，我挑了实验上暴露的主要晶面算，并考虑了含氧中间体的覆盖度效应。

从我的工程经验看，电化学材料dft计算 最怕的是”拿体相能带去判断催化活性”——催化只看表面吸附,体相带隙几乎不相关。正确的做法是盯死表面吸附自由能这个描述符。我给客户的建议是优先推进火山顶附近的 B 类配方，并通过应变或掺杂把 A 类往火山顶推。更多 OER/ORR 计算案例，见 [电化学材料dft计算](https://www.keyanxueshu.com/category/dft/)；需要定制，[电化学材料dft计算](https://www.keyanxueshu.com/) 上有入口。