第4章 | 2026新版VASP吸附与催化教程:吉布斯自由能计算！

HER反应吉布斯自由能计算的基本理论

吉布斯自由能计算的基本公式

氢析出反应（HER）的催化活性评价核心在于准确计算反应过程中各步骤的吉布斯自由能变化（ΔG）。根据密度泛函理论（DFT）结合计算氢电极（CHE）模型，自由能计算的完整表达式为：

其中，ΔE 是通过DFT计算得到的电子能量差，ΔE_ZPE代表零点振动能修正，TΔS 为温度与熵变的乘积修正项。该公式将电子结构计算与热力学修正相结合，是电催化自由能台阶图构建的基础 。

在实际计算中，ΔE通常通过吸附能公式获得：

式中，E_H@surface 表示氢原子吸附在催化剂表面的总能量，E_surface为清洁表面的能量，E_H2为气态氢分子的总能量。

计算氢电极模型（CHE）的核心原理

计算氢电极（Computational Hydrogen Electrode, CHE）模型是Nørskov团队在2004年提出的理论框架，其核心思想是将标准氢电极（SHE）条件下的质子–电子对（H⁺ + e⁻）的能量定义为：

在298.15 K和1 atm标准条件下，这一关系使得电化学反应的电极电位可以转化为化学势，从而实现DFT计算与实验电化学条件的对接。根据CHE模型，在U=0 V（相对于标准氢电极）条件下，理想催化剂每一步反应的自由能变化均为0 eV，这为催化剂的活性评价提供了统一的热力学基准。

在HER反应中，CHE模型使得ΔG_H*（氢吸附自由能）的计算公式简化为：

其中0.24 eV是ZPE修正与熵修正的综合贡献值 。

HER反应路径与自由能台阶图

HER反应的基本步骤

HER在酸性介质中的反应路径包括三个基本步骤：

Volmer步骤（吸附步骤）：H⁺ + e⁻ + * → H*

Heyrovsky步骤（电化学脱附）：H* + H⁺ + e⁻ → H₂ + *

Tafel步骤（化学脱附）：2H* → H₂ + 2*

不同的电催化剂HER决速步骤不同，相应机理也不同，一般可分为Volmer-Heyrovsky机理或Volmer-Tafel机理 。

决速步的Tafel斜率识别

Tafel斜率是评价析氢反应动力学及反应机理的重要参数，可用于判定电极表面析氢反应的决速步骤。若Volmer反应为决速步骤，则Tafel斜率约为120 mV·dec⁻¹；若斜率约为30 mV·dec⁻¹，表明Tafel步骤是决速反应，HER通过Volmer-Tafel机制进行；如果Tafel斜率约40 mV·dec⁻¹，说明Heyrovsky过程为决速步骤，反应机制为Volmer-Heyrovsky机制 。

例如，NiCu/NF电极的Tafel斜率为56.9 mV/dec，表明Volmer-Heyrovsky反应为决速步骤 。商业化Pt/C催化剂的Tafel斜率在30 mV/dec及以下，表明Tafel过程为决速步骤 。

决速步与理论过电位的关系

决速步对应的最大自由能变化（ΔGmax）与理论过电位（η）存在直接关联：

该关系表明，理论过电位直接取决于反应路径中自由能升高最大的步骤，这为催化剂活性优化提供了明确的目标。

修正项的物理意义与计算方法

零点振动能（ZPE）的修正

零点振动能（Zero-Point Vibrational Energy, ZPE）来源于量子力学的基本原理，即使是在绝对零度下，分子仍然具有残余振动能量。在HER催化中，ZPE修正主要涉及吸附氢原子（H*）与气态氢分子（H₂）之间的零点能差异 。

在VASP等主流DFT软件中，ZPE的计算通过频率分析实现：

结构优化后频率计算：在优化后的几何结构基础上，计算所有振动模式的本征频率

ZPE求和：将所有振动频率代入公式求和，得到总ZPE值

典型计算示例中，H*吸附在Pt(111)表面时，其ZPE值约为0.10-0.15 eV，而气态H₂的ZPE值约为0.27 eV，两者的差值即为ZPE修正项 。

熵变（TΔS）的修正

在HER反应过程中，熵变主要来自两部分：

1.气态氢分子的平动、转动和振动熵：气态H₂分子具有显著的平动自由度、转动自由度和振动自由度，在298 K时的标准熵值约为130.6 J/(mol·K)

2.吸附态氢原子的振动熵：吸附态H*被束缚在催化剂表面，自由度减少，振动熵显著低于气态分子

根据Nørskov团队的标准处理方法，在298.15 K条件下，HER自由能计算中的熵变修正为-0.24 eV 。

修正的综合影响

ZPE和熵修正虽然数值量级不大（0.1-0.3 eV量级），但对于催化剂活性的精确评价至关重要。文献统计显示，忽略ZPE和熵修正可能导致ΔG_H*误差达到0.2 eV以上，进而使催化剂活性评级出现偏差 。

以Pt催化剂为例，不考虑ZPE和熵修正时计算的ΔE_H约为-0.35 eV，加入ZPE修正（+0.14 eV）和熵修正（+0.10 eV）后，最终ΔG_H为-0.11 eV，与实验测得的-0.09 eV高度吻合 。

本章要点总结

吉布斯自由能计算知识点梳理

HER反应吉布斯自由能计算的基本理论：吉布斯自由能计算的基本公式、计算氢电极模型（CHE）的核心原理

HER反应路径与自由能台阶图：HER反应的基本步骤、决速步的Tafel斜率识别、决速步与理论过电位的关系

修正项的物理意义与计算方法：零点振动能（ZPE）的修正、熵变（TΔS）的修正、修正的综合影响

下一步学习建议

下一章将正式引入本次教程的核心—HER台阶图。我们将从HER自由能台阶图的定义与理论基础、基本公式与修正项、反应路径与自由能台阶图方面详细介绍HER台阶图。