氢溢流不是单个 H 在表面平移，而是 H₂ 在金属上解离、H 跨越金属—载体边界、在载体表面或体相迁移并参与后续还原的连续过程。本文给出 VASP 界面模型、端点筛选和分段 CI-NEB 方法，并说明氧化物可还原性怎样改变主导路径。
一条 NEB 能覆盖全过程吗？
完整过程可拆为四类事件：H₂ 在金属团簇或界面解离；H* 从金属位转移到邻近载体 O 或金属位；H 在载体表面逐位迁移；迁移后的 H 与目标物种或晶格氧反应。氢溢流是一组相连的基元步骤，而不是一个跨越整块载体的单一跳跃。每一步都具有独立初态、终态和能垒。
若把金属解离、跨界面转移和长距离扩散塞进同一条插值路径，中间图像会同时改变 H—H、M—H 和 O—H 键，反应坐标难以收敛。先分别计算局部步骤，再用共同中间态连接能量剖面。共同中间态的组成、总电荷和结构优化口径必须一致。
图 1 比较多种氧化物负载金属体系，并把活性变化与氧空位形成能联系起来。它说明载体不只是几何背景：载体的可还原性、表面 O 配位和缺陷稳定性会改变 H 的接受位及后续反应。载体性质决定跨界面后的 H 以羟基、氢化物还是缺陷相关物种存在。
图1 不同氧化物负载 Ru/Ni 催化体系的反应行为及氧空位形成能关联。载体可还原性影响氢溢流。DOI: 10.1039/d2sc00871h。
最小模型要保留哪些结构？
界面模型从稳定氧化物晶面与金属团簇开始。载体 slab 的层数、真空层和固定层先做收敛测试；团簇尺寸要能提供 H₂ 解离位和与载体接触的边界位。方法案例可使用 Ru₅/TiO₂(110)、Ru₅/CeO₂(001) 或 Ru₅/WO₃(001)，但具体团簇构型仍需在目标体系中筛选。
团簇在载体上的锚定位至少比较多个平移、旋转和界面配位方式，并分别测试含氧空位与无空位表面。跨界面路径必须建立在已经确认稳定的金属—载体接触结构上。优化后记录金属—O 键、团簇形貌、Bader 电荷和局域磁矩，否则 NEB 中的主要能量变化可能来自团簇重排。
H₂ 解离既可能发生在金属顶部，也可能发生在金属—O 边界。分别建立分子吸附态、同裂后的两个 M—H、异裂后的 M—H/O—H 端点，并用 NEB 或 Dimer 比较。图2中的 XANES/EXAFS 信息反映不同载体上金属结构与合金化状态差异，可作为界面模型选择的实验约束。
图2 TiO₂、CeO₂、WO₃ 负载 Ru/Ni 体系的 XANES 与 EXAFS 结构信息。界面模型应与可观测配位环境相符。DOI: 10.1039/d2sc00871h。
H 从金属转移到哪里？
跨界面初态通常是边界金属上的 H*，终态可以是载体表面 O—H、载体金属—H或氧空位附近的 H。每个终态先单独优化，比较吸附能和局域键合。若 H 自动返回金属，说明该猜测不是稳定终态；若 H 诱导晶格 O 脱出，则过程已经包含载体还原，需要另设中间态。
原子映射保持 H 的身份连续，金属团簇和载体原子顺序不变。跨界面端点应分别对应清楚的 M—H 与 O—H 局域极小态。插值时处理周期边界，确保 H 沿实际界面距离移动；对于 M—H 断裂与 O—H 形成同步的路径，增加图像密度以解析键交换。
图 3 把 H₂ 解离、H 从 Ru₅ 转移到氧化物、表面迁移和还原 Ni 位点四步列在同一势能剖面中。这种分段方式允许比较究竟是边界跨越还是载体扩散控制整体速率，也能避免把多个键变换压缩成一个难以解释的最高点。
图3 Ru₅/氧化物体系的氢溢流基元步骤和 DFT 势能剖面。路径由 H₂ 解离、跨界面转移、载体迁移和还原步骤组成。DOI: 10.1039/d2sc00871h。
表面扩散端点怎样接起来？
载体表面迁移可按相邻 O 位、金属位或空位逐段建立。例如 O₁—H → O₂—H 是一个局部跳跃，O₂—H → O₃—H 是下一段；每一段终态也是下一段初态。对 TiO₂、CeO₂ 和 WO₃，应根据表面晶格和氢稳定形式分别选择路径，不能只沿几何最短直线。
CI-NEB 参数可从 5–7 个图像、`LCLIMB = .TRUE.`、`SPRING = -5` 开始，电子与离子收敛阈值按能垒精度测试。质子位置与电子局域位置共同决定可还原氧化物中的迁移态。若 H 迁移伴随小极化子转移，初末态还要设置一致且可追踪的局域电子状态；必要时用初始磁矩、占据矩阵或结构畸变引导电子局域化。
图 4 的 H₂/D₂ 原位红外结果显示温度和载体会改变表面氢交换行为。计算上可把同位素差异留在 ZPE 与频率修正中，电子总能面主要由相同原子势能决定。过渡态验证检查一个沿 H 跳跃方向的虚频，并对最高能结构两侧微扰。
图4 Ru/TiO₂、Ru/CeO₂ 与 Ru/WO₃ 的 H₂/D₂ 原位红外交换结果。载体与温度共同影响表面迁移行为。DOI: 10.1039/d2sc00871h。
多段能量剖面统一以 H₂ 气相加清洁界面为零点，也可以统一以解离吸附态为零点，但同一图中不能切换。每段活化能 Ea,i 取本段过渡态相对本段初态的差值，整体表观瓶颈还需结合前序中间态能量。最高局部能垒与最高绝对能量点是两个不同量，微观动力学使用哪一个取决于反应网络和覆盖状态。
Ea,i = E(TSi) – E(ISi) 式(1)
表面路径和体相路径怎样区分？
可还原氧化物能通过 O—H 形成、氧空位和电子局域稳定迁移氢，但不同材料的优选方向并不相同。比较 TiO₂、CeO₂ 和 WO₃ 时，模型记录至少包含晶面、缺陷状态、H 物种、迁移方向、每段能垒与终态反应能。不能只比较一个数值就归纳整个载体的溢流能力。
图 5 概括三类氧化物在不同温度下的表面与体相迁移路径。判断氢溢流方向需要同时比较表面跳跃、次表层进入和体相迁移。若表面扩散能垒低而体相进入能垒高，H 主要沿表面传播；若体相路径更有利，表面模型只能覆盖初始阶段。
图5 TiO₂、CeO₂ 和 WO₃ 在不同温度下的氢溢流路径示意。不同载体可能选择表面或体相迁移。DOI: 10.1039/d2sc00871h。
VASP 目录可设为 `00_interface`、`01_H2_ads`、`02_H2_split`、`03_metal_to_support`、`04_surface_hops`、`05_subsurface`、`06_target_reaction`。每段保存端点结构、NEB 图像、能量剖面、最大力、H 的键合对象和局域电荷。若载体含缺陷，再建立化学计量表面与缺陷表面的平行序列。
电子结构旁证怎样补充？
吸附能与能垒之外，还可用差分电荷密度、Bader 电荷、PDOS 和 O—H/M—H 频率确认氢物种。AIMD 可检查有限温度下界面 H 是否保持局域或发生快速迁移。静态 NEB 给出候选势能路径，实验约束和有限温度采样决定它是否对应实际主导过程。
当 H 从金属转移到载体 O 后，体系可能同时出现质子与还原电子。若电子仍局域在金属团簇，而质子已经形成 O—H，路径接近质子转移；若电子同步局域到载体金属位，则属于耦合的 H⁺/e⁻ 迁移。两种状态即使 H 的坐标接近，能量也可能不同。计算记录应同时保存 H 的键长、Bader 电荷、局域磁矩和带隙态分布。
缺陷浓度也会改变迁移网络。孤立氧空位可能稳定邻近 H，却也可能形成深能量陷阱，使后续扩散变慢；较高缺陷浓度则会改变表面电荷与可用 O 位数量。较稳定的 H 终态并不自动对应更快的氢溢流，连续路径仍取决于后续迁移能垒。以同一超胞分别计算无缺陷、单空位和远离界面的空位模型，可判断缺陷促进的是跨界面进入还是长距离传播。
温度比较需要把各段势垒转成统一的速率表达，并考虑吸附态覆盖度。跨界面转移很快而 H₂ 解离很慢时，金属供氢决定通量；H₂ 解离容易但载体表面跳跃较慢时，H 会积累在界面附近。由静态能垒确定候选控制步骤后，可再用过渡态理论或微观动力学计算温度依赖。
要点：氢溢流由 H₂ 解离、跨界面转移、载体迁移和后续反应等多个基元步骤组成。
要点：每段 NEB 都需要稳定端点、连续原子映射和统一能量零点。
要点：氧化物的可还原性、缺陷和电子局域化共同决定表面或体相迁移方向。
