析氧反应（Oxygen Evolution Reaction, OER）作为电解水制氢、金属-空气电池等能源转换技术中的关键阳极反应，涉及四个质子耦合电子转移（PCET）过程，具有复杂的动力学特征。在OER反应过程中，催化剂表面形成一系列含氧中间体，主要包括OH、O和OOH三种关键吸附物种。
这些中间体的吸附构型直接决定了反应的能垒和催化活性。DFT计算研究表明，OH与OOH中间体的吸附自由能之间存在稳定的标度关系：
ΔG（OOH） = ΔG（OH） + 3.2 ± 0.2 eV
这一关系适用于大多数金属氧化物和氢氧化物催化剂表面。该标度关系的可靠性与±0.2 eV的偏差反映了不同材料表面的分散性，而非计算误差本身。
吸附质演化机理（Adsorbate Evolution Mechanism, AEM）是目前最为广泛接受的OER反应机制。在AEM框架下，氧中间体经历吸附→去质子化→偶联→解吸的连续顺序过程，具体路径为：
OH → O → OOH → O₂
在碱性电解质中，AEM的四步反应路径分别为：
1.第一步：M + OH⁻ → M-OH + e⁻（*OH吸附）
2.第二步：M-OH + OH⁻ → M-O + H₂O + e⁻（*O形成）
3.第三步：M-O + OH⁻ → M-OOH + e⁻（*OOH形成）
4.第四步：M-OOH + OH⁻ → M + O₂ + H₂O + e⁻（O₂脱附）
其中，第三步（O到OOH的转化）通常是反应的速率决定步骤（RDS），具有最高的吉布斯自由能垒。
OH中间体的吸附构型
*OH中间体是最先吸附在催化剂表面的含氧物种。建模时通常在金属位点或桥位上放置OH基团，键长控制在1.85 Å左右。OH的氧原子与金属中心形成σ键，氢原子则朝向溶剂或表面空位。这种构型使得OH能够有效地从溶剂中接受电子，启动OER催化循环。
O中间体的吸附构型
O中间体是OH去质子化后形成的物种。从几何结构角度看，O仅由单个氧原子组成，通常以端基配位形式吸附在金属中心上，键长较短，通常在1.8-2.0 Å范围内。这种紧密的吸附构型导致O与金属中心之间形成较强的共价相互作用，但也使得后续*OOH的形成更加困难。
OOH中间体的吸附构型
*OOH中间体包含两个氧原子和一个氢原子，具有多个自由度，在构型搜索时需要尝试多种可能的结构。根据DFT计算结果，*OOH在催化剂表面主要有两种吸附构型：
单点吸附构型：整个*OOH基团作为一个整体吸附在单个金属原子上（如Bi原子）。
异位吸附构型：-OH基团吸附在一种金属原子上，而-O基团吸附在另一种金属原子上，形成双金属协同吸附结构。
*OOH的建模中，两个氧原子之间的O-O键长通常控制在1.4 Å左右，而O-H键长则保持在-0.97 Å左右。由于势能面复杂，单一构型优化往往无法找到全局最低点，需要通过多构型搜索来确定最稳定的吸附构型。
晶格氧机制（Lattice Oxygen Mechanism, LOM）是一种突破传统AEM限制的新兴机制。LOM的核心特征在于催化剂自身的晶格氧直接参与反应，与吸附的OH⁻或H₂O结合生成O₂，随后在催化剂表面形成氧空位，后续由电解质中的H₂O或OH⁻填补空位。
LOM机制的关键反应路径包括：
1.晶格氧活化（O_lat* + H⁺ + e⁻ → OH*）
2.O-O键形成（O_lat* + *OH → O₂ + V_O + H⁺ + e⁻）
3.空位再生（V_O + H₂O → OH* + H⁺ + e⁻）
这种机制绕过了AEM中*OOH形成和脱附的步骤，从而避免了线性标度关系的限制，理论上可以将过电位降低至远小于370 mV的水平。
LOM与AEM在中间体吸附构型上存在显著差异：
LOM路径中，由于晶格氧参与反应，中间体*OOH的形成步骤被绕过，取而代之的是晶格氧与吸附氢/羟基的直接耦合。
氧空位的形成是LOM的标志性特征，氧空位可作为潜在的活性位点，改变反应中间体的分布和能垒。
使用密度泛函理论（DFT）对OER中间体吸附构型进行优化是研究吸附能的关键方法。DFT计算通常采用以下策略：
1.模型构建：构建包含催化剂表面和吸附中间体的超胞模型，允许所有原子在三维空间弛豫。
2.收敛标准：原子受力收敛阈值通常设定为0.01 eV/Å，确保结构优化的准确性。
3.交换关联泛函：RPBE泛函被广泛用于吸附能计算，能够较好地描述过渡金属氧化物的表面吸附行为。
自由能台阶图是分析OER反应路径的重要工具。每个反应步骤的吉布斯自由能变化（ΔG）决定了反应的热力学难易程度。对于OER的四步反应：
ΔG₁：*OH的形成自由能
ΔG₂：*O的形成自由能
ΔG₃：*OOH的形成自由能
ΔG₄：O₂的脱附自由能
其中ΔG₃通常是最大值，决定了整个反应的能垒。
吸附自由能是评估OER催化活性的关键描述符。DFT计算揭示了OOH与OH的吸附能之间存在着稳定的线性关系：
ΔG（OOH） = 0.81ΔG（OH）+ 3.30 eV（R² = 0.93）
基于这一标度关系，OER最佳吸附自由能差为2.46 eV，对应最小热力学过电位。
火山图将催化活性与吸附能（如ΔG（O）- ΔG（OH）或ΔG（OOH） – ΔG（OH））关联起来，火山顶点对应最佳催化剂，其吸附能平衡了反应中间态的稳定性，使得各步骤的自由能差尽可能接近。
吸附质演化机理（AEM）中的中间体吸附构型：介绍了AEM机制的反应路径、各中间体的几何吸附构型
晶格氧机制（LOM）中的中间体吸附构型：介绍了LOM机制的反应特征、LOM与AEM的中间体吸附差异
DFT计算中的吸附构型优化方法：介绍了结构优化与几何参数、自由能计算与吸附能分析、吸附能标度关系与火山图
下一章将正式引入本次教程的核心—OER自由能台阶图。我们将从自由能台阶图的基本概念与物理意义、OER反应的四电子转移路径、自由能台阶图的计算方法与构建步骤方面详细介绍OER自由能台阶图。
