手机版

氢吸附自由能计算

发布时间：2026-08-28 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

析氢反应（Hydrogen Evolution Reaction，HER）是电解水制氢的核心阴极反应，也是研究电催化机理的模型反应。其总反应 2H⁺ + 2e⁻ → H₂ 可拆分为多个基元步骤，而吸附氢中间体 H 的结合强度是决定催化活性的关键描述符。根据 Nørskov 提出的计算氢电极（CHE）模型，析氢反应的自由能图可由氢吸附自由能 ΔG_H 刻画：当 ΔG_H* ≈ 0 时，H 的吸附与脱附均接近热中性，催化活性最高；ΔG_H 过正（吸附太弱）或过负（吸附太强）都会使反应受限。

经典的 HER 火山图将实验交换电流密度 log(i₀) 对 ΔG_H* 作图，呈火山形曲线，铂（Pt）位于火山峰顶，而 MoS₂ 边缘位、Au、Ag 等材料位于两翼。本项目通过 DFT 计算多种金属与二维材料的 ΔG_H*，构建火山图并定位催化活性最优的材料。

二、计算方法与理论背景

2.1 氢吸附自由能

ΔG_H* = ΔE_H* + ΔZPE − TΔS，其中 ΔE_H* 为氢原子在表面的吸附能（相对 ½H₂ 的参考），ΔZPE 与 ΔS 为吸附态与气相氢气的零点能和熵差。在标准条件下（T = 298 K，pH = 0，U = 0 V），吸附氢的 ΔG_H* 表达式简化为 ΔG_H* = ΔE_H* + 0.24 eV（经验修正项，包含 ZPE 与熵变）。

2.2 吸附能计算

采用周期性 slab 模型，氢原子吸附在表面最稳定位点（如 fcc 空位、顶位）。吸附能 E_ads = E(slab+H) − E(slab) − ½E(H₂)。计算需保证 slab 厚度、真空层与 k 点收敛。

2.3 火山图

火山图以”吸附强度”为横轴、催化活性（交换电流密度）为纵轴。其物理本质是 Sabatier 原理：吸附太弱，质子还原为 *H 的步骤（Volmer 步）受限；吸附太强，H 脱附产氢（Heyrovsky/Tafel 步）受限。火山峰对应 ΔG_H ≈ 0 的理想催化剂。

三、计算设置与模型构建

软件：VASP；泛函：GGA-PBE，PAW 赝势，截断能 400–500 eV。
表面模型：fcc(111) 4×4 超胞，4 层 slab，真空层 ≥ 15 Å，底层固定。
k 点：Monkhorst–Pack 网格，标准精度。
研究体系：Pt(111)、Ni、Co、Cu、Au、Ag(111) 与 MoS₂ 边缘位。
收敛判据：原子力 < 0.02 eV/Å，能量 < 1×10⁻⁵ eV。

四、结果分析与案例图解读

图 2 的左面板给出 HER 火山图。灰色虚线为火山曲线两翼，各金属的数据点（彩色圆点）分布在火山两侧：Pt(111) 的 ΔG_H* = −0.09 eV，非常接近零点，位于火山峰顶附近，对应最高活性；MoS₂ 边缘位 ΔG_H* = +0.08 eV 同样接近峰顶，解释了 MoS₂ 作为廉价析氢催化剂的活性来源；Au（+0.30 eV）、Ag（+0.38 eV）因吸附过弱落在右翼，活性较低；Ni（−0.27 eV）、Co（−0.34 eV）、Cu（−0.23 eV）因吸附过强落在左翼。火山图直观揭示了”适度吸附”的催化原则，是论文中”催化剂活性描述符”章节的核心图。

图 2 的右面板给出 Pt(111) 表面在 U = 0 V 时的氢吸附自由能阶梯图（红色圆点线）。路径为 H⁺ + e⁻ → H → ½H₂(g)：第一步 ΔG_H = −0.09 eV（轻微放热，吸附容易），第二步由 *H 脱附产氢需克服 +0.09 eV（轻微吸热）。两步骤的自由能差很小，说明 Pt 表面 H 吸附接近热中性，这正是 Pt 催化活性最优的电子结构根源。该自由能图将”火山峰”的宏观描述落实到具体的基元步骤能量学。

两面板结合，火山图建立”吸附强度-活性”的全局规律，自由能图提供 Pt 的微观证据，共同回答”为什么 Pt 是最好的 HER 催化剂”这一核心问题，并为替代催化剂筛选提供描述符判据。