第9章 | 2026新版VASP吸附与催化教程：HER反应模型构建全解 —— 表面、界面、单原子位点模型构建对比

HER的DFT建模首先要回答：当前模型解释的是本征吸附、界面反应，还是单原子位点的局域结构？表面 slab、含水界面与 M–N_x 位点并非彼此替代，而是服务于不同问题。本文按反应步选择模型，并给出可执行的 VASP 计算路径。

图1  Au@Ru 核壳纳米颗粒上的碱性 HER 反应示意与局域水环境。图中给出了 Ru 表面位点及 H/O 物种的界面分布。DOI: 10.1038/s41467-023-41030-1。

三类模型的边界

裸表面先给出可比较的吸附基线

表面模型用周期性slab比较晶面、缺陷、掺杂和吸附位点，是筛选 H* 结合趋势的第一步。它应在一致的超胞、层数、覆盖度和参考态下比较 ΔG_H*，用于判断位点是否存在过强或过弱吸附。

式(1)  ΔG_H* = ΔE_H* + ΔE_ZPE − TΔS

ΔG_H* 接近零不等于实际电极一定最快。真空 slab 缺少水分子取向、离子屏蔽、双电层与电位响应，因此它更适合回答“本征吸附趋势”，而不是独立给出完整电化学机理。

slab 的横向尺寸、层数和真空层都需收敛测试；非对称 slab 还应说明偶极修正。若不同材料采用不同超胞或覆盖度，比较结果可能来自模型差异而非材料差异。

界面模型把反应环境放回结构中

若研究碱性 HER 的水解离、Heyrovsky 步的质子/电子转移，或活性随 pH、电位明显变化，就应加入界面描述。显式水提供局域氢键与取向，隐式溶剂近似介电屏蔽；加入水层不等于已经实现恒电位。

图2  Ru 表面碱性 HER 的原位拉曼与 DFT 结果。d 图比较了不同 Ru 表面/氧化层模型上的水与吸附物构型。DOI: 10.1038/s41467-023-41030-1。

界面模型的关键难点是构型采样：不同水层初始构型可能给出不同氢键网络和势垒。实用策略是先用表面模型筛选，再对代表性位点进行界面验证。

界面模型还应区分“结构环境”与“热力学参考”。显式水可以改善局域氢键的描述，却不能替代自由能修正；隐式溶剂可以提供连续介电响应，也不会自动给出真实的水解离路径。若文章报告的是ΔG_H*，需要交代 CHE 参考、温度与 ZPE/熵修正；若报告的是势垒，则需交代初态、终态、路径搜索方法以及水分子的参与方式。把这些条件写清楚，读者才能判断不同计算之间能否并列。

单原子位点需要把局域配位当作变量

金属原子不是一个固定的圆点

单原子模型的关键不是金属元素名称，而是配位数、配位原子、局域电荷、自旋态、载体缺陷与轴向吸附物。理想M–N₄ 构型不应被默认视为电化学条件下唯一活性结构。

单原子位点至少比较两类信息：不同配位/自旋态的稳定性，以及这些结构上的 H* 吸附、关键反应步和电荷变化。若介质会占位，还需加入竞争吸附结构。

对M–N₄/C 位点，应检查金属是否偏离平面、H* 或 OH* 是否诱发重构、不同自旋态能量排序，以及邻近缺陷是否改变电子转移。未做这些比较时，ΔGH* 只代表一个设定构型。

图3  氮掺杂石墨烯上 Pt 单原子与 Pt 团簇的显微表征及形成示意。该图说明载体、配位与团簇化均会影响单原子位点。DOI: 10.1038/ncomms13638。

按反应步升级模型，避免一次性堆复杂度

Volmer、Heyrovsky 与 Tafel 的建模重点不同

酸性Volmer步可先从 H* 吸附入手；碱性条件下水解离更敏感。Heyrovsky 步涉及溶剂化质子/电子与界面构型，通常需要界面验证；Tafel 步则重点检查相邻 H* 距离与覆盖度。

在VASP实操中，可先完成超胞大小、层数、真空层、k 点、自旋设置与固定层策略的收敛测试，再比较吸附构型。进入界面阶段后，需明确水层数量、初始构型、是否采用隐式溶剂、是否处理偶极修正及其适用条件。进入过渡态或能垒计算前，还要确认初态和终态处于同一模型与同一参考体系；否则能垒或自由能差即使数值漂亮，也缺少可比性。

判断顺序很简单：先明确要筛选材料还是解释机理；再判断决速步是否涉及水、离子或电位；最后检查活性中心是否会重构。这三问决定是否由 slab 升级到界面与单原子配位搜索。

图4  不同 Pt 单原子配位模型的自由能路径与 PDOS 分析。图中展示局域配位对 H* 中间体和电子结构的影响。DOI: 10.1038/s41467-020-14848-2。

把模型落实为 VASP 的计算流程

先统一模型与能量参考

三类模型都应先定义同一套参考体系：清洁表面（或清洁位点）、吸附态、气相H₂，以及在需要时加入水、OH* 或带电界面的参考状态。先由总能量得到吸附能，再加入零点振动能与熵修正；比较时需统一超胞、覆盖度、参考态和电位定义。

式(2)  E_ads(H*) = E(*H) − E(*) − 1/2 E(H₂)

式(3)  ΔG_H* = ΔE_ads + ΔE_ZPE − TΔS

若讨论电位或 pH，计算氢电极（CHE）可把质子—电子对写入自由能表达式；常见写法为 ΔG(U,pH) = ΔG(0,0) − neU + kBT ln(10)×pH。这里的 n、U 与 pH 必须和具体反应步对应。公式给出的是计算框架，不能代替对溶剂、水层和界面构型的实际处理。

VASP 与辅助软件

VASP 用于几何优化、静态能量、频率与 NEB 势垒计算；POSCAR 定义 slab、吸附物或 M–N_x 位点，INCAR 控制电子与离子迭代，KPOINTS 与 POTCAR 决定采样和赝势。VASPKIT、ASE 或 pymatgen 可用于批量生成结构、检查输入和提取能量，VESTA 用于人工核对晶面、吸附位点与结构重构。最终仍以 VASP 输出文件与收敛测试作为判断依据。

建议按六步推进计算

确定体系与反应步：明确酸/碱性环境、关注 Volmer、Heyrovsky 还是 Tafel 步，以及需要比较的位点。

建立清洁模型：对 slab 检查层数、真空层、横向超胞和固定层策略；对单原子位点同时准备不同配位与自旋初猜。

完成收敛测试：ENCUT、k 点、超胞尺寸和层数没有通用最优值；文中参数只能作为示例起点，需按体系验证。

优化初态与吸附态：分别优化清洁表面、*H、*OH、*H₂O 等必要构型，并记录总能量与最终磁矩。

进行静态与热力学校正：在收敛结构上用更严格电子收敛进行单点计算；对关键吸附态做频率分析，获得ΔE_ZPE 与熵项。

计算势垒并回到机理：若决速步涉及水解离或氢气形成，可用 CI-NEB 搜索路径；初态、终态和中间图像必须处于同一模型与同一参考条件。

参数需要写清适用边界

作为起始示例，结构优化可取 EDIFF=1×10⁻⁵ eV、EDIFFG≈−0.02 eV/Å，并以 IBRION=2、ISIF=2 进行固定晶胞下的离子弛豫；静态能量计算可设 NSW=0。这些参数仅是起点；磁性单原子、强偶极表面、显式水层或弱相互作用体系仍需检查自旋、偶极修正、色散和电子展宽。

建模时最容易混淆的四件事

把不同表面、不同覆盖度、不同溶剂状态下的ΔG 直接排在一起比较。应先统一参考态、超胞和电位定义。

把“单原子结构优化收敛”当成“电化学条件下稳定”。应额外检查配位重构、竞争吸附和可能的迁移路径。

只因实验在碱性介质中测试，就把所有结论交给真空 slab。若水解离或 OH* 占位参与决速步，界面影响需要进入验证。

为追求复杂而直接使用巨大水层模型，却没有对多个水构型或关键参数做检查。模型复杂度不能替代可验证性。

本章要点总结

表面slab用于建立可比较的吸附与位点筛选基线。

界面模型用于补充水、离子、电场和质子/电子转移等环境因素。

单原子位点的核心是局域配位、自旋和载体协同，不是只放置一个金属原子。

模型升级应由反应步与待验证假设驱动，而不是由模型大小驱动。

下一步学习建议

下一章可以进一步讨论 H* 吸附自由能如何从 VASP 总能量走到 ΔG_H*：需要哪些结构能量、ZPE 与熵修正如何进入公式，以及不同覆盖度和参考态为什么会改变结果解释。把这一步算清楚，才能把“模型选择”落实到可比的自由能数据上。

#华算科技朱老师#VASP吸附与催化计算教程 #HER反应模型#表面模型#电化学界面#单原子催化剂#DFT计算#VASP计算