氢析出反应（HER）的催化活性评价核心在于准确计算反应过程中各步骤的吉布斯自由能变化（ΔG）。根据密度泛函理论（DFT）结合计算氢电极（CHE）模型，自由能计算的完整表达式为：
其中，ΔE 是通过DFT计算得到的电子能量差，ΔE_ZPE 代表零点振动能修正，TΔS 为温度与熵变的乘积修正项。该公式将电子结构计算与热力学修正相结合，是电催化自由能台阶图构建的基础。
在实际计算中，ΔE通常通过吸附能公式获得：
式中，EH@surface 表示氢原子吸附在催化剂表面的总能量，Esurface为清洁表面的能量，EH2为气态氢分子的总能量。
计算氢电极（Computational Hydrogen Electrode, CHE）模型是Nørskov团队在2004年提出的理论框架，其核心思想是将标准氢电极（SHE）条件下的质子–电子对（H⁺ + e⁻）的能量定义为：
在298.15 K和1 atm标准条件下，这一关系使得电化学反应的电极电位可以转化为化学势，从而实现DFT计算与实验电化学条件的对接。根据CHE模型，在U=0 V（相对于标准氢电极）条件下，理想催化剂每一步反应的自由能变化均为0 eV，这为催化剂的活性评价提供了统一的热力学基准。
在HER反应中，CHE模型使得ΔGH*（氢吸附自由能）的计算公式简化为：
其中0.24 eV是ZPE修正与熵修正的综合贡献值 。
零点振动能（Zero-Point Vibrational Energy, ZPE）是量子力学的基本概念，源于海森堡不确定性原理。即使在绝对零度下，分子仍然具有残余振动能量，其值为：
其中，ħ为约化普朗克常数，ω_i 为第i个振动模式的频率，3N为分子的振动自由度（N为原子数）。
在HER催化中，ZPE修正主要涉及吸附氢原子（H*）与气态氢分子（H₂）之间的零点能差异。由于吸附态的H*受到表面束缚，其振动频率显著低于气态H₂分子，导致ZPE值出现差异，这一差异通常在0.1-0.3 eV范围内 。
在VASP等主流DFT软件中，ZPE的计算通过频率分析实现。具体计算步骤包括：
1.结构优化后频率计算：在优化后的几何结构基础上，设置IBRION=5（有限差分法计算Hessian矩阵）和NSW=0（固定晶胞参数），计算所有振动模式的本征频率 。
2.积分法频率计算：设置NFREE=2采用中心差分方法，可更精确地获取振动频率信息 。
3.ZPE求和：将所有振动频率代入公式求和，得到总ZPE值。
典型计算示例中，H*吸附在Pt(111)表面时，其ZPE值约为0.10-0.15 eV，而气态H₂的ZPE值约为0.27 eV，两者的差值即为ZPE修正项 。
ZPE修正虽然数值不大（0.1-0.3 eV量级），但对于催化剂活性的精确评价至关重要。文献统计显示，忽略ZPE修正可能导致ΔG_H*误差达到0.2 eV以上，进而使催化剂活性评级出现偏差 。
以Pt催化剂为例，不考虑ZPE时计算的ΔE_H约为-0.35 eV，加入ZPE修正（+0.14 eV）和熵修正（+0.09 eV）后，最终ΔG_H为-0.12 eV，与实验测得的-0.09 eV高度吻合 。
在HER反应过程中，熵变主要来自两部分：
1.气态氢分子的平动、转动和振动熵：气态H₂分子具有显著的平动自由度（3个）、转动自由度（2个）和振动自由度（1个），在298 K时的标准熵值约为130.6 J/(mol·K)，对应TS项约为0.41 eV 。
2.吸附态氢原子的振动熵：吸附态H*被束缚在催化剂表面，自由度减少，振动熵显著低于气态分子。通常近似认为吸附态的熵贡献接近0，或采用计算得到的振动频率进行精确修正 。
根据Nørskov团队的标准处理方法，在298.15 K条件下，HER自由能计算中的熵变修正为：
这一数值来源于：气态H₂的TS值为0.41 eV，吸附态H*的TS近似为0.17 eV，两者差值为0.24 eV 。
具体数据表格展示了不同体系中熵变修正的详细数值：
|
体系 |
ZPE修正 (eV) |
TΔS修正 (eV) |
综合修正 (eV) |
|
Pt(111)-H* |
0.14 |
0.10 |
0.24 |
|
MoS₂边缘-H* |
0.12 |
0.08 |
0.20 |
|
高熵合金-H* |
0.16 |
0.11 |
0.27 |
熵变修正具有显著的温度依赖性，其关系可表示为：
当温度从298 K升高到350 K时，TΔS值相应增加约17%。在实际工况下（如高温电解池），需要考虑这一温度效应对自由能计算的影响 。
HER自由能计算的理论框架与公式：吉布斯自由能计算的基本公式、计算氢电极模型（CHE）与标准氢电极
零点振动能（ZPE）的物理意义与计算方法：零点振动能的物理本质、ZPE的计算方法、ZPE对催化剂筛选的影响
熵变（TΔS）修正的来源与数值：熵变的物理来源、TΔS的具体数值、温度依赖性
下一章将正式引入本次教程的核心—吉布斯自由能计算。我们将从HER反应吉布斯自由能计算的基本理论、HER反应路径与自由能台阶图、修正项的物理意义与计算方法方面详细介绍吉布斯自由能计算。
