氢吸附自由能计算

这张图是本项目做 Pt(111) 表面 HER 催化剂筛选时的标准输出：左边是三个典型吸附位点（Top、Bridge、Hollow），右边是它们对应的氢吸附自由能 ΔG_H*。对 HER 来说，Sabatier 原理已经讲得很明白——ΔG_H* 越接近 0，催化剂活性越高；负太多说明氢粘得太牢脱附困难，正太多则吸附都困难。本案例 Hollow 位 ΔG = −0.29 eV、Bridge 位 −0.15 eV、Top 位 +0.08 eV，一眼就能看出 Hollow 是最优活性位。

这张图在说什么

Hollow 位（三重空位）让 H 与三个 Pt 原子配位，成键最强，所以 ΔG 最负；Bridge 位两重配位，居中；Top 位只与一个 Pt 作用，最弱，ΔG 甚至转正。最有用的是把这三个位点和右边柱子连起来看：Hollow 绿柱扎进负区，说明氢自发吸附；Top 红柱翻到正区，说明这里基本不吸附。实际催化发生在配位环境复杂的缺陷或台阶位，但单晶面三个位点的排序是理解活性的起点。

自由能不是单点能

很多新手直接用 DFT 吸附能 E_ads 报结果，但催化活性看的是自由能 ΔG_H*。本项目计算时按下面公式把零点和熵补进去：

ΔG_H* = ΔE + ΔZPE − T·ΔS

– ΔE：DFT 吸附能（H* − slab − ½H₂）

– ΔZPE：零点能修正，通常 −0.05 ~ −0.15 eV

– −T·ΔS：H* 在表面近似无平动/转动，熵项主要扣掉气相 ½H₂ 的熵，室温下约 +0.20 eV

不加 ZPE 和熵修正，ΔG 能比 ΔE 偏出 0.2 eV 以上，催化剂排序可能直接翻过来。

VASP 计算设置速查

| 参数 | 取值 | 说明 |

|——|——|——|

| ENCUT | 450 eV（≥1.3×默认） | 截断能不足会低估吸附强度 |

| KPOINTS | 6×6×1 / 8×8×1 | 表面/slab 需加密面内 k 点 |

| ISMEAR | 0；SIGMA=0.05 | 金属表面用高斯展宽 |

| IVDW | 11 | DFT-D3 色散修正 |

| ISPIN | 2 | Pt 无强磁性但仍开自旋更稳 |

| ISIF | 2 | 弛豫离子，固定晶胞 |

位点怎么取

本项目取吸附位点的规则：

1. 先对原胞做结构优化，得到干净表面。
2. 在 Top、Bridge、Hollow 三种高配位点分别放 H，z 方向距最近金属原子 1.0–1.2 Å 初猜。
3. 每个位点从 2–3 个不同 z 高度松手弛豫，避免掉到同一个局部最小。
4. 取能量最低者作为该位点的 ΔE。

数值可信吗

– **ZPE 修正**：H 在 Hollow 位振动频率低，ZPE 小；Top 位频率高，ZPE 大。本项目对每种吸附构型单独算振动频率，不偷懒用固定值。

– **覆盖度**：本图是单个 H 的极限覆盖度。实际 HER 表面 H 覆盖可能 0.25–0.5 ML，相邻 H 的排斥会让 ΔG 变正，不能直接用单点结果外推高覆盖。

– **溶剂化**：水溶液中 H 和界面水有氢键网络，真空 ΔG 与水相可能差 0.05–0.10 eV。本项目对高精度需求会加隐式溶剂或 AIMD 采样。

这套结果能怎么用

把 Hollow 位 ΔG = −0.29 eV 放进火山图：Pt 位于火山峰左侧偏弱吸附区，但总体接近峰顶，解释了为什么 Pt 是 HER 基准催化剂。如果做合金（Pt-Ni、Pt-Co），目标就是把 ΔG 往 0 eV 调；若 ΔG 还更负，说明合金化方向错了。

本项目做完 Pt(111) 基线后，会把同一套流程批量套到台阶面、缺陷面、合金面上，比较谁的 ΔG 最接近 0，从而筛选出候选催化剂。

[氢吸附自由能计算](https://www.keyanxueshu.com/category/dft/) 是这类催化筛选的核心工具；完整计算流程与参数模板可访问科研学术网首页 https://www.keyanxueshu.com/ 获取。

局限

本图只覆盖清洁 Pt(111) 的极限覆盖度；真实电化学界面有溶剂、电场、动态覆盖度，真空 ΔG 只能做趋势判断，不能当活度定量。