晶面差异会怎样改变ΔGH*？

引言

同一种材料换一个晶面，ΔGH* 可能就会从接近 0 eV 变成偏正或偏负。原因不只是表面原子排得不一样，而是表面配位数、台阶密度、应变、电子态和吸附位点一起变了。晶面差异改变 ΔGH*，本质上是改变 H* 看到的局域原子环境。

所以，比较晶面不能只把 (111)、(100)、(110) 三个 slab 优化完，再直接排吸附能。要先确认每个晶面的 slab 厚度、真空层、覆盖度、k 点密度、固定层数和 H* 位点枚举是否一致。否则得到的差别里，可能既有晶面效应，也有模型设置带来的偏差。

晶面差异先改变哪些表面环境？

低指数晶面和高指数晶面差在哪里？

低指数晶面通常更平整，表面原子配位数较高；高指数晶面常带有台阶、边和 kink 位点，低配位原子更多。低配位位点往往更容易和 H* 成键，但不一定代表 HER 一定更好。如果 H* 绑得太牢，后续 H₂ 生成或脱附也会被拖慢。HER 需要的是适中氢吸附，而不是越强越好。

以金属催化剂为例，(111) 面常出现三重 hollow 位，(100) 面常有四重 hollow 或桥位，(110) 面更开放，顶位、桥位和沟槽位都可能参与吸附。晶面一变，H* 的候选位点也要跟着重新枚举，不能把一个晶面的最低能位点直接套到另一个晶面。

图1 HER/OER 中晶面调控、选择性封端和表面调制策略示意，说明暴露晶面会改变反应位点环境。DOI：10.1007/s40820-023-01024-6

为什么同一材料不同晶面不能只比 slab 总能？

slab 总能和原子数、层数、真空层、表面面积都有关，不能直接拿来判断哪个晶面更适合 HER。真正需要比较的是同一参照下的 H* 吸附自由能、表面能、关键位点稳定性和实验暴露比例。若某个高指数晶面 ΔGH* 很理想，但表面能很高、实验中难以大量暴露，也不能把它写成最现实的主导晶面。

晶面比较最好分成两层：第一层回答“这个晶面如果存在，H* 吸附强弱如何”；第二层回答“这个晶面是否容易形成并稳定暴露”。ΔGH* 说明位点活性，表面能和形貌证据说明这个位点有多可能出现。

图2 酸性与碱性 HER 中 Volmer、Heyrovsky 和 Tafel 步示意，提示不同晶面还可能改变水解离和氢脱附步骤。DOI：10.1007/s40820-023-01024-6

ΔGH* 为什么会随着晶面改变？

H* 吸附位点变了，自由能就会变吗？

会，而且经常是主要来源。H* 在顶位、桥位、三重 hollow 位和台阶边位上的成键数不同，金属—H 键长不同，电子态耦合也不同。平整晶面可能让 H* 稳定在 hollow 位，开放晶面可能让 H* 靠近低配位金属原子。晶面改变 ΔGH* 的第一条路径，就是改变最低能 H* 构型。

ΔGH* = ΔEH* + ΔZPE − TΔS

ΔEH* 来自 VASP 总能差，ΔZPE 和 TΔS 来自振动频率或常用热力学校正。晶面比较时，最容易出问题的是 ΔEH*，因为不同晶面上的 H* 初始位点、覆盖度和优化约束稍有变化，最低能构型就可能不同。

电子态也会让 ΔGH* 改变吗？

会。不同晶面的表面原子配位数不同，d 态宽度、d 带中心、费米能级附近态密度都会变。低配位金属原子常有更窄的 d 带和更强的局域态，可能增强 H* 吸附；合金晶面还会叠加应变效应和配体效应。

但 d 带中心不能代替 ΔGH*。它更适合解释趋势，比如某个晶面的 d 带中心上移，H* 吸附随之增强；若出现不一致，就要回到具体吸附位点、覆盖度和表面重构去看。描述符用于解释趋势，最终活性判断仍要回到自由能和构型。

图3 Pd3Au 不同晶面 DOS、Bader 电荷和表面模型示例，展示晶面会改变电子结构与局域电荷分布。DOI：10.3390/molecules28073169

VASP 里比较晶面时怎样做？

先让 slab 模型站在同一套标准上

不同晶面要尽量使用接近的 slab 厚度、真空层、固定层数和表面覆盖度。k 点密度要按表面尺寸换算，不能只复制同一个 KPOINTS。对于极性表面或上下不对称 slab，要考虑偶极修正；对于高指数晶面，要检查优化后台阶是否塌陷或发生过大重构。晶面比较的前提，是模型差异只来自晶面本身，而不是参数和边界条件。

Eads(H*) = Eslab+H − Eslab − 1/2EH₂

这个式子用于得到 H* 吸附能。后续再加零点能和熵修正，得到 ΔGH*。如果比较多个晶面，EH₂ 必须来自同一套计算设置；Eslab 必须对应各自未吸附 H 的同一晶面 slab；H 覆盖度也要保持一致。

H* 位点应该枚举到什么程度？

每个晶面至少要放顶位、桥位和 hollow 位；如果是台阶面，还要加 step edge、terrace、kink 或沟槽位。优化后记录初始位点、最终位点、H—金属键长、表面原子位移和相对能量。若某个初始位点总是迁移到另一个位点，这本身就是晶面信息，说明该晶面真正稳定的吸附中心并不是初始猜测的位置。

INCAR 参考示例可以这样写。晶面比较不靠参数堆得复杂，而靠对照组保持一致。

# ======== HER 晶面吸附优化参考示例 ========
SYSTEM = slab_111_H_relax
ENCUT  = 520
EDIFF  = 1E-5
EDIFFG = -0.03
ISMEAR = 1
SIGMA  = 0.10
IBRION = 2
NSW    = 120
ISIF  = 2
LDIPOL = .TRUE.
IDIPOL = 3

# ======== 静态单点参考示例 ========
SYSTEM = slab_111_H_static
ENCUT  = 520
EDIFF  = 1E-6
IBRION = -1
NSW    = 0
LORBIT = 11
LCHARG = .TRUE.

图4 Pt 与 SbPt 不同晶面/结构的 HER 自由能示例，说明晶面和表面组成会共同影响 ΔGH*。DOI：10.1007/s40820-023-01024-6

晶面结果怎样接到 HER 机理？

酸性和碱性 HER 要看的重点一样吗？

不完全一样。酸性 HER 中，ΔGH* 更直接地连接 H⁺ 还原和 H₂ 生成；碱性 HER 中，水解离和 *OH 脱附也会参与限制。某个晶面 ΔGH* 接近 0 eV，但如果水解离能力弱，在碱性条件下未必表现最好。晶面影响 HER，不只通过 H*，也可能通过 H₂O、*OH 和界面水结构。

因此，若体系用于碱性 HER，可以在 ΔGH* 之外补 H₂O 吸附、H—OH 断键构型、*OH 脱附自由能或水层氢键信息。若只讨论酸性 HER，重点可放在 H* 吸附和 Heyrovsky/Tafel 步；若实验电解液是碱性，直接用酸性 HER 逻辑会少一块关键证据。

图5 NiS2 不同晶面 H 吸附和 HER 性能示例，展示晶面暴露与 H* 吸附强度之间的联系。DOI：10.1007/s40820-023-01024-6

为什么 ΔGH* 接近 0 eV 也不能直接下最终判断？

ΔGH* 接近 0 eV 是很重要的指标，但它只覆盖氢吸附热力学。实际 HER 还受表面稳定性、电导、活性面积、晶面暴露比例、反应动力学和电解液环境影响。高指数晶面可能有漂亮的 ΔGH*，却因为表面能高、重构明显或暴露比例低，在实验中贡献有限。

更稳的表述是：某晶面上的目标位点提供更接近热中性的 H* 吸附，并可能降低 HER 中氢吸附相关步骤的热力学代价；若能结合表面能、形貌表征和稳定性数据，就可以进一步说明它是否可能成为主要活性晶面。ΔGH* 是判断晶面活性的核心证据之一，但不是完整机理的全部证据。

图6 TaC/GaN 体系中不同 Miller 晶面对 H 吸附与自由能路径的影响示例。DOI：10.1007/s40820-023-01024-6

晶面差异最后怎样写？

把晶面、位点和自由能三件事绑在一起

写晶面差异时，不建议只说“某晶面 ΔGH* 最接近 0 eV，所以活性最好”。更清楚的写法是：先说明该晶面暴露了哪类低配位或特殊配位位点，再说明 H* 最稳定地吸附在哪里，最后用 ΔGH*、PDOS 或差分电荷密度解释为什么这个位点更适合 HER。

如果多个晶面的 ΔGH* 差距很小，就不要强行拉开优劣。可以写成“多个晶面都可能提供接近热中性的 H* 吸附，差别更多来自晶面暴露比例、表面稳定性和水解离能力”。晶面结论要同时回答结构从哪里来、H* 吸在哪里、自由能为什么变。

本章要点总结

要点一：晶面差异会改变表面原子配位数、吸附位点、电子态和表面稳定性，进而影响 ΔGH*。

要点二：不同晶面比较时，slab 厚度、真空层、固定层数、k 点密度和 H* 覆盖度要尽量保持一致。

要点三：ΔGH* 接近 0 eV 说明 H* 吸附较适中，但不能替代表面能、形貌暴露和动力学证据。

要点四：碱性 HER 里还要关注 H₂O 解离和 *OH 脱附，不能只沿用酸性 HER 的 H* 逻辑。