



OER 机理判断最容易混在一起的地方,是把所有含氧中间体都放进 *OH、*O、*OOH 三步自由能里,然后默认它就是吸附质演化机制。对于含氧化物、羟基氧化物和高价过渡金属位点,晶格氧可能参与 O-O 成键,氧空位和金属-氧能带关系也会改变路径。AEM 和 LOM 的差别不只是图上有没有晶格氧,而是氧原子来源、电子空穴位置和闭合循环是否不同。后续建模要先把这个差别写进反应式。
这里讨论的是 VASP/DFT 能真正落地的 OER 机理区分:AEM 与 LOM 的反应式要分开建,自由能、氧空位形成能、PDOS、Bader 电荷、磁矩和同位素实验旁证要对应到各自位置。若只比较一个过电位数字,容易把两个机制压成同一条曲线;若把晶格氧参与直接写成结论,又会缺少能量与电子结构证据。
机制区分的最小计算包应同时包含反应式、能量闭合、结构继承和电子结构旁证。




AEM 先回答什么?
吸附质演化机制(adsorbate evolution mechanism, AEM)把表面活性位看作反应中心,典型路径为 * → *OH → *O → *OOH → *。在这个框架里,O-O 键通常由吸附态 *O 与水或 OH⁻ 相关的外来氧形成,表面晶格氧主要提供配位环境和电子结构背景。AEM 的核心计算量是每一步吸附态自由能差,而不是单独某个中间体总能。因此,它对参照态和吸附位一致性很敏感。
AEM 计算要固定同一表面模型、同一活性位和同一参照态。干净 slab、*OH、*O、*OOH 使用相同固定层、k 点、赝势和磁构型;含电位修正时,常用计算氢电极近似把 H⁺+e⁻ 与 1/2 H₂ 相连。比较 NiFe-LDH、FeOOH 或过渡金属氧化物时,不同终止面和不同金属位不能混成一条自由能序列。
LOM 多算了哪一层?
晶格氧机制(lattice oxygen mechanism, LOM)把部分晶格氧作为反应物处理,路径中会出现表面晶格氧参与 O-O 成键、氧空位形成和后续补氧。单晶格氧机制只让一个晶格氧进入产氧循环,双晶格氧机制则可能包含两个邻近晶格氧的耦合。图1 给出了 AEM、sLOM 和 dLOM 的示意边界,适合先用来校准反应对象。

图1 AEM、单晶格氧机制和双晶格氧机制的示意。图中把吸附氧与晶格氧区分开,有助于判断 O-O 成键对象来自表面吸附态还是晶格氧。 DOI: 10.1002/advs.202303321。
LOM 不能只在 *OOH 构型旁边标一个晶格 O。它至少需要三类计算:第一,晶格氧被氧化或形成氧空穴的电子结构证据;第二,O-O 成键后是否留下稳定氧空位或过氧型中间体;第三,空位被 H₂O/OH⁻ 补回后能否闭合循环。如果氧空位形成能过高、补氧步骤不可行,LOM 路径即使某一步自由能较低,也不能单独作为完整机制。
氧空位形成能和补氧自由能决定 LOM 是否只是局部成键片段,还是可以形成完整产氧循环。




AEM 台阶图怎样写?
AEM 的四电子路径可以写成 ΔG₁=*OH−*,ΔG₂=*O−*OH,ΔG₃=*OOH−*O,ΔG₄=O₂−*OOH。常用表达为 ΔG(U,pH)=ΔE+ΔZPE−TΔS−neU+kBT ln10×pH,其中 ΔE 来自 VASP 总能差,ZPE 和熵需要按气相分子与吸附态分别处理。参数不要和公式挤在一句话里,先把自由能口径说清楚,再给可复查的输入示例。
INCAR 参考示例可这样排,具体数值只用于说明 OER 吸附态静态计算口径:
/// INCAR 参考示例:OER 吸附态静态计算
SYSTEM = oer_aem_static
ISTART = 0
ICHARG = 2
ENCUT = 520
EDIFF = 1E-6
ISMEAR = 0
SIGMA = 0.05
ISPIN = 2
MAGMOM = 按过渡金属价态和初始磁态设置
LORBIT = 11
LASPH = .TRUE.
LDAU = .TRUE. # 仅在过渡金属强关联体系中按校准使用
LWAVE = .TRUE.
LCHARG = .TRUE.
这组输入适合放在结构优化后的静态单点或电子结构复查阶段。ENCUT 和 k 点要分别做收敛,U 值要来自文献、线性响应或同一项目内的校准,MAGMOM 要至少测试几组合理初值;若体系不使用 DFT+U,可删去 LDAU 相关设置。真正可比较的是同一套泛函、U 值、k 点和磁构型下得到的一组吸附态能量。
同一催化剂上,AEM 的限制步来自最大 ΔG_i 对 1.23 eV 的偏离。若 *OH→*O 过高,说明脱氢生成 *O 的热力学代价较大;若 *O→*OOH 过高,说明 O-O 成键前驱体不够有利。这里的判断仍然在吸附态路径内完成,不能因为 *OOH 台阶下降就自动说明晶格氧已经参与。
LOM 台阶图怎样补齐?
LOM 的自由能序列需要把晶格氧位点和空位态写进反应式。含晶格氧参与的 O-O 初态和终态先单独优化,随后计算氧空位表面 E(slab−O)、补水或补 OH⁻ 后的再生结构,并把缺陷态表面与完整表面放在同一能量基准下比较。LOM 的台阶图必须包含空位形成和补氧闭合,否则只是一段局部成键过程。这些结构要分别保存,不能在后处理阶段临时混用。
图2 的氧缺陷表征说明,同一 NiFe-LDH/FeOOH 体系中,氧缺陷含量会随合成条件变化。建模时不能把有缺陷与无缺陷表面混成一个参照;AEM 在完整表面上比较 *OH、*O、*OOH,LOM 则加入含氧空位或氧空穴的表面,并对比有无缺陷时的自由能差异。

图2 NiFe-LDH/FeOOH 中氧缺陷和价态变化的表征。图中 XPS、EPR 和结构示意可作为建立有缺陷表面模型的实验旁证。 DOI: 10.1002/advs.202303321。




氧 2p 和金属 d 带怎样读?
LOM 常与氧 2p 态靠近费米能级、金属 d 态和氧 p 态强杂化、氧空穴更易形成有关。VASP 静态计算可用 LORBIT=11 输出 PDOS,并在同一费米能级对齐口径下比较完整表面、缺陷表面和吸附态表面。若氧 2p 态在 EF 附近贡献明显增加,只能说明晶格氧更可能参与氧化还原,仍需和结构与自由能闭环一起判断。
氧 2p 态靠近费米能级是 LOM 倾向的电子结构线索,不是独立机制判据。
以 NiFeOOH 表面为例,若完整表面的 AEM 曲线中 *O→*OOH 台阶最高,而邻近氧空位模型中晶格氧参与的 O-O 终态更稳定,且补氧步骤自由能没有显著升高,LOM 才具备竞争路径的热力学基础。相反,PDOS 即使出现氧 2p 态上移,只要氧空位形成能和补氧台阶仍然较大,结果更适合作为电子结构倾向,而不是机制切换证据。
同位素和 DFT 怎样互相约束?
图3 同时给出 ¹⁸O 同位素质谱、AEM/LOM 自由能、过电位和 PDOS。¹⁸O 信号回答晶格氧是否进入 O₂,自由能图回答对应路径能否闭合,PDOS 说明氧 2p 与金属 d 态是否提供电子结构支撑。同位素信号可以证明晶格氧参与产氧,DFT 则需要说明哪条晶格氧路径能量上可闭合。

图3 氧同位素、AEM/LOM 自由能、过电位和 PDOS 的联合分析。图中展示了缺陷引入后不同路径的自由能差异和电子结构变化。 DOI: 10.1002/advs.202303321。
结果归档时,把四类输出分开列出:AEM 中 *OH、*O、*OOH 的相对自由能;LOM 中晶格氧成键、氧空位和补氧步骤的自由能;完整与缺陷表面的 O 2p/金属 d PDOS 对比;Bader 电荷和局域磁矩变化。若只保留一张台阶图,后续很难判断差异来自机制本身、缺陷浓度,还是来自表面终止面变化。




模型矩阵怎样设计?
模型矩阵从四组结构展开:完整表面、含一个氧空位表面、AEM 吸附态表面、LOM 晶格氧成键/空位再生表面。每组模型都应记录层数、固定层、真空层、U 值、磁构型和吸附位点。若体系是 LDH 或羟基氧化物,还要明确 H 覆盖度和脱氢程度,因为这些变量会直接改变金属价态与氧空穴分布。
若完整表面上 Fe 位 *OOH 更稳定,而缺陷表面上邻近晶格氧成键终态更稳定,完整表面的 *OOH 能量不能和缺陷表面的晶格氧终态直接拼接为同一条路径。更清楚的比较方式是分别画 AEM(完整)、AEM(缺陷)、LOM(缺陷) 三条曲线,再比较限制步与闭合循环。

图4 NiFe-LDH/FeOOH 氧缺陷工程与 OER 性能关联。图中把缺陷调控、电子结构和反应性能放在同一体系内对照。 DOI: 10.1002/advs.202303321。
哪些输出最值得留下?
结构优化后先检查 O-O 键长、金属-氧键长和氧空位附近配位恢复;电子结构再看 PDOS、差分电荷密度、Bader 电荷、局域磁矩和必要的 COHP/ICOHP。AEM 更依赖吸附态自由能矩阵,LOM 更依赖晶格氧活化、空位形成和再生步骤。判断机制时,能量闭合、氧来源证据和电子结构旁证三者缺一项,结论都会失去完整性。
图5 展示了同一体系的 OER 极化曲线、Tafel 斜率和稳定性信息。同一材料体系中的性能差异要求 DFT 模型保持可追溯的结构关系,不能只在彼此无关的模型之间比较能量。VASP 模型中,不同缺陷浓度、不同 Fe/Ni 邻位环境应作为平行结构分别计算,而不是让一个理想表面承载所有实验特征。

图5 NiFe-LDH/FeOOH 的 OER 电化学性能对比。图中性能差异为缺陷调控与 AEM/LOM 路径比较提供实验背景。 DOI: 10.1002/advs.202303321。




要点:AEM 主要比较 *OH、*O、*OOH 吸附态自由能,LOM 需要把晶格氧参与、氧空位形成和补氧闭合纳入同一能量基准。
要点:PDOS 中氧 2p 态靠近费米能级只能提供倾向,不能单独替代自由能闭合和结构证据。
要点:有缺陷表面和完整表面应分别建立路径,不能把不同参照下的中间体能量拼成一条台阶图。
要点:较完整的 VASP 判断应同时保留 AEM/LOM 自由能、氧空位形成能、补氧步骤、PDOS、Bader 电荷和局域磁矩。
