氢吸附自由能计算

氢吸附自由能（ΔG_H*）是电催化析氢（HER）、氢氧化与储氢材料评价的核心描述符，直接决定活性位点对氢的中间体结合是否”恰到好处”。本项目计算 ΔG_H* 时，严格区分气相、真空表面与带电界面的处理，因为三者给出的值可能相差 0.3 eV 以上。下文围绕氢吸附自由能的计算协议与判读展开。

为什么是自由能而不是能量

吸附过程在室温、液相、电位下发生，熵与溶剂贡献不可忽略。本项目计算的是 ΔG_H* = ΔG(gas H₂→surface H*) − ½ G(H₂)，而非单纯电子能量差。忽略零点能（ZPE）与熵项，单一吸附构象的能量差会偏离真实自由能 0.1–0.2 eV，足以改变对”火山图”位置的判断。

计算协议（以金属表面 HER 为例）

1. 表面模型：用足够厚 slab（≥ 4 层，底层固定），建 (2×2) 或更大 supercell 降低覆盖度效应。
2. 吸附构象：枚举顶位、桥位、空位（fcc/hcp），各做结构优化找最稳位点。
3. 能量分解：E_ads = E(slab+H) − E(slab) − ½ E(H₂)。本项目用 H₂ 在 PBE 下能量作基准。
4. 自由能修正：加 ZPE 与熵差（气相 H₂ 与吸附 H 的振动/平动熵），用 ΔG = ΔE + ΔZPE − TΔS。
5. 电位修正：对电化学界面加 U（外加电位）与参比修正，本项目用计算氢电极（CHE）模型把电位折进自由能。

方法要点

• 泛函：PBE 普遍偏低，本项目加 Hubbard U（对 d 带金属）或用 RPBE/r2SCAN 改善吸附能；色散（D3/D4）对弱吸附位点有帮助。
• 覆盖度：高覆盖度下 H-H 排斥显著，本项目按目标条件选覆盖度，避免低估稳定化。
• 溶剂：显式水层（AIMD）比隐式更准但贵；本项目对关键位点补做几 ps AIMD 确认溶剂化效应。

判读：火山图与活性

ΔG_H* ≈ 0 eV 的位点理论活性最高。本项目算得某合金表面桥位 ΔG_H* = −0.08 eV，接近理想，与实验过电位低吻合；而纯金属母相 ΔG_H* = −0.45 eV，过强吸附锁氢，活性差——解释了合金化的增益来源。

坦诚的局限

CHE 模型假设 H* 与 H₂ 平衡、忽略动力学势垒；真实过电位还受 Tafel 步骤与传质影响。本项目把 ΔG_H* 定位为”热力学活性描述符”而非完整动力学预测。氢吸附与电催化计算可进一步查本站的第一性原理计算栏目。更多科研计算服务见科研学术网首页。