水分解反应（H₂O decomposition）是电催化、光催化和热催化等领域的核心反应过程，包括两个基本半反应：阴极的析氢反应（HER）和阳极的析氧反应（OER）。在水的分解过程中，H₂O分子需要经历O-H键的断裂和H-H键、O-O键的重新形成，这是一个涉及多个电子转移和化学键重组的复杂过程。
水分解路径包含以下关键组成部分：
反应物：H₂O（气态或吸附态）
中间体：H₂O*（吸附水）、OH*（羟基吸附物）、H*（吸附氢原子）、O*（吸附氧原子）
产物：1/2 H₂（氢分子）、1/4 O₂（氧分子）
过渡态：反应路径上的能量最高点
在实际催化过程中，H₂O解离通常经历物理吸附、O-H键断裂及进一步解离三个步骤，关键中间体包括H₂O*、OH和H，其键长和解离能垒决定了反应活性。
密度泛函理论（Density Functional Theory, DFT）是计算催化反应路径和能垒的核心方法。DFT通过近似求解薛定谔方程来预测材料的性质，其中Kohn-Sham方程是常见的实现方法 。在电催化领域，常用的交换关联泛函包括广义梯度近似（GGA），其中Perdew-Wang 1991（PW91）GGA是最流行的泛函之一 。
势能面描述了分子能量随几何结构变化的关系，是理解反应机理的基础工具。反应坐标λ表示参与反应的所有原子的位置，为了精确描述势能，必须使用多维坐标系，这将在多维空间中产生一个势能面 。
势能面与反应坐标的关系揭示了结构、稳定性、反应性和选择性与能量之间的联系。反应坐标图显示反应物、产物和过渡态的结构，势能面还揭示了反应的热力学和动力学信息，如反应的放热或吸热性质，以及反应速率与活化能的关系。
HER和OER反应的吉布斯自由能计算采用CHE（计算氢电极）模型，其基本公式为：
其中：
ΔE：通过DFT计算得到的电子能量差
ΔEZPE：零点振动能修正
TΔS：温度与熵变的乘积修正项
在298.15 K条件下，HER自由能计算中的综合修正值约为0.24 eV 。
过渡态是指分子从反应物转变为产物过程中能量最高的中间状态，处于反应物向生成物转化的势能面鞍点位置。过渡态决定了反应的活化能（activation energy），因此是理解反应速率和选择性的重要因素 。
Nudged Elastic Band（NEB）方法是寻找两个稳定状态之间的最小能量路径和过渡态的常用方法 。NEB通过在两个初始状态之间创建等间距的图像并应用恢复力来保持图像间距，从而高效地确定过渡态的位置 。
CI-NEB（Climbing Image NEB）是NEB方法的改进版本，通过移除最高能量图像的“弹簧“来捕获过渡态，将一个图像放置在鞍点位置，从而提供更准确的能量壁垒估计 。
NEB与CI-NEB方法对比：
|
方法 |
特点 |
精度 |
|
常规NEB |
低分辨率图像，插值估计 |
活化能偏低 |
|
CI-NEB |
一个图像在鞍点，精确搜索 |
能量壁垒准确 |
过渡态搜索的标准流程包括：
初始路径构建：基于反应物和产物的几何结构生成初始路径
NEB计算：使用NEB方法寻找最小能量路径
CI-NEB优化：通过攀爬图像方法精确定位过渡态
频率计算：验证过渡态结构（应存在一个虚频）
IRC计算：确认过渡态连接正确的反应物和产物
在VASP计算中，过渡态计算通常先使用CI-NEB方法初步搜索，能量收敛精度为10-5 eV，力的收敛精度为0.05 eV/Å 。然后使用改进二聚体法（Improved Dimer Method, IDM）对CI-NEB搜索到的过渡态结构进行精确计算，力的收敛精度为0.05 eV/Å，能量的收敛标准为10-5 eV 。
H2O分解反应的基本路径与化学背景：介绍了水分解反应概述、反应路径的构成要素
DFT计算方法与理论基础：介绍了密度泛函理论的基本框架、势能面与反应坐标、吉布斯自由能计算
过渡态搜索方法：介绍了过渡态的定义与特征、NEB与CI-NEB方法、过渡态搜索流程
下一章将正式引入本次教程的核心—OER反应路径与活性位点。我们将从OER反应路径、催化剂活性位点、反应中间体方面详细介绍。
