在 MoS₂、CoS₂、NiSₓ 和多金属硫化物中，S 空位常被写成提升析氢反应（hydrogen evolution reaction, HER）活性的关键因素。真正进入计算时，问题并不是“有空位就更好”，而是空位把哪个原子暴露出来、H* 更适合落在哪类位点、缺陷浓度是否改变相邻位点，以及该结构在反应电位下是否还能保持。S 空位的核心作用，是把原本配位饱和或电子态远离费米能级的表面，改造成低配位金属位与缺陷诱导态共同参与的吸附环境。
图1 MoS₂ 基面 S 空位浓度与 HER 自由能的关系。图中比较了不同空位浓度下 H* 吸附自由能接近热中性区域的变化。DOI: 10.1038/ncomms15113。
空位只是少了一个 S 吗？
在完整硫化物表面，S 原子通常位于最外层，部分金属原子被 S 配位包围。对 MoS₂ 这类层状硫化物而言，基面在理想结构中常表现为较弱的 H* 吸附，真正活跃的常是边缘位或缺陷位。去掉一个表面 S 后，邻近 Mo、Co、Ni 等金属原子的配位数降低，局部键长和电子态重新分布，H* 不再只面对一个平整的 S 终止面，而会感受到一个金属暴露与缺陷电荷共同形成的局域环境。
图1 中可以看到，S 空位浓度变化会让 ΔGH* 从偏离 0 eV 的区域逐渐靠近热中性。这里不能把空位理解成“多一个吸附坑位”这么简单。空位改变 HER 位点的第一层含义，是把吸附中心从完整 S 终止面转移到空位邻近的低配位金属原子及其桥连区域。因此建模时至少要比较空位中心、邻近金属顶位、金属-硫桥位和远离空位的完整基面位。
缺陷态为什么会影响 H*？
S 空位还会改变费米能级附近的投影态密度（projected density of states, PDOS）。完整 2H-MoS₂ 基面的 S-p 态和 Mo-d 态对 H-s 轨道的耦合较弱，H* 吸附常偏弱；引入 S 空位后，邻近金属 d 态更容易在费米能级附近出现，H-s 与金属 d 态的杂化增强，吸附能随之改变。这个判断需要用 PDOS、差分电荷密度和结构键长一起确认，不能只靠一个 ΔGH* 数值给出电子结构结论。
对不同硫化物，缺陷态的方向也不完全一致。某些体系中 S 空位让 H* 吸附更接近 0 eV；另一些体系中，空位邻近金属位可能吸附过强，使后续 H₂ 形成或 Heyrovsky 步骤变得不利。判断 S 空位是否真正改善 HER，必须同时看 H* 热力学、电子态来源和反应环境下的结构稳定性。
图2 MoS₂ 表面 S 相关物种和脱硫路径的自由能分析。图中说明 S 空位生成与电位条件、表面物种稳定性有关。DOI: 10.1038/ncomms15113。
先算空位形成能吗？
如果只比较完整表面和含 S 空位表面的 ΔGH*，容易忽略一个前提：含空位模型是否在给定硫化学势、电位和反应气氛下可达到。空位形成能可写为式(1)，其中 E(slab – VS) 为含 S 空位表面能量，Eslab 为完整表面能量，μS 为 S 的化学势。μS 不能随意取一个固定数，常需要结合 S-rich、S-poor 或实验处理条件设置边界。
Evac = E(slab – VS) + μS – Eslab 式(1)
在 VASP 中，完整 slab 和含空位 slab 应保持相同的晶面、层数、真空层、固定层、泛函、POTCAR 版本、ENCUT、KPOINTS 和色散修正设置。若研究二维 MoS₂ 基面，可先用 4×4 或更大超胞放置单个 S 空位，再用更高空位浓度模型做趋势对照。只有当空位形成代价和 HER 吸附改善同时被量化时，才能说明缺陷位点具有计算上的可行性。
H* 应该放在哪些位置？
HER 计算通常先比较 H* 在完整基面、空位中心、空位邻近金属顶位、邻近 S 顶位、桥位和边缘位的吸附。吸附能可按式(2)计算，其中 Eslab-H* 是吸附态总能，Eslab 是对应缺陷表面总能，EH₂ 是气相 H₂ 总能。对同一组位点比较时，所有模型必须使用相同超胞和相同电子设置，否则吸附能差异会混入数值误差。
Eads(H*) = Eslab-H* – Eslab – 1/2EH₂ 式(2)
自由能校正进一步写成式(3)。ΔZPE 与 ΔS 来自频率计算或可靠文献修正，ΔGU 表示电位项，pH 项需按所用参照标尺说明。对酸性 HER，H⁺ + e⁻ 可通过计算氢电极（computational hydrogen electrode, CHE）连接到 1/2H₂；对碱性 HER，还应显式考虑 H₂O 解离和 *OH 稳定性。S 空位只让 ΔGH* 接近 0 eV，并不自动说明 Volmer 步、Heyrovsky 步或 Tafel 步都同步优化。
ΔGH* = ΔEH* + ΔZPE – TΔS + ΔGU + ΔGpH 式(3)
图3 富 S 空位 MoS₂ 中电子态调节与 HER 路径示意。图中对比了不同 MoS₂ 相和空位结构对 H₂O/H* 转化的影响。DOI: 10.1038/s41467-022-33636-8。
结构优化怎样设置？
第一步是建立完整表面和含 S 空位表面。对层状硫化物，slab 厚度和真空层先用完整表面收敛，再在相同超胞中删去表面 S。结构优化建议使用自旋极化，过渡金属硫化物还要给出合理的初始 MAGMOM；EDIFFG 可按力收敛目标设置为示例值 -0.02 eV/Å，具体阈值需按体系收敛测试确认。若体系存在明显层间作用，vdW 修正也应在完整表面、缺陷表面和吸附态中保持一致。
第二步是构建 H* 位点矩阵。每个初始构型优化后，要核对 H-M 键长、H-S 距离、空位邻近金属位是否发生大幅重构，以及 H* 是否迁移到非预期位置。对于高空位浓度模型，还要确认两个空位之间是否产生明显相互作用。同一篇对比中，完整基面位、空位邻近位和边缘位应分组讨论，不能把不同几何环境直接排成一个单一活性序列。
电子结构怎样辅助判断？
完成吸附自由能后，可以对代表位点进行静态单点，输出 DOSCAR、CHGCAR、AECCAR 和差分电荷密度所需结果。PDOS 应重点比较 H-s、邻近金属 d 态和 S-p 态在费米能级附近的重叠；差分电荷密度则用于观察 H* 与空位邻近金属之间的电子密度重排。Bader 电荷可作为定量辅助，但净电荷变化很小并不代表没有局域成键变化。
图3 展示的富 S 空位体系说明，缺陷引入后的电子态可改变 H₂O 到 H₂ 的路径效率。对计算结果而言，较合理的证据组合是：ΔGH* 接近 0 eV、H-M 键长处于合理范围、PDOS 显示 H-s 与金属 d 态存在相互作用、差分电荷密度显示界面电子密度重排。电子结构图应服务于解释为什么某个位点吸附更合适，而不是替代自由能或能垒判断。
图4 MoS₂ 缺陷结构中 H* 吸附位点、PDOS 与 HER 自由能示意。图中给出了缺陷调节后 H* 位点和自由能路径的差异。DOI: 10.1038/s41467-022-29929-7。
碱性环境还要补什么？
在酸性 HER 中，ΔGH* 常是第一判断量；在碱性 HER 中，水解离通常会进入关键步骤。此时仅有 H* 吸附自由能还不够，还需要比较 H₂O 在空位邻近位的吸附取向、O-H 键活化、*OH 脱附或迁移，以及可能的 H₂O → H* + OH⁻ Volmer 步能垒。S 空位若增强 H₂O 吸附但同时使 *OH 过度稳定，反而可能让位点被 *OH 覆盖。
实际建模可先做三组结构：空位表面 + H*、空位表面 + H₂O、空位表面 + *OH。若 H₂O 吸附态中 O 原子靠近空位邻近金属，且 O-H 键明显拉长，可进一步设置 NEB 端点计算水解离能垒。碱性 HER 中，S 空位是否有效取决于 H₂O 活化、H* 形成和 *OH 去除三者之间的平衡。
空位会不会在反应中继续变化？
硫化物表面在电化学条件下可能发生进一步脱硫、补硫、氧化、羟基化或相转变。图2 已经提示 S 相关物种的稳定性与电位有关，因此含 S 空位模型最好不要只作为固定几何使用。更可靠的做法是把空位形成能、H* 自由能、*OH 或 O* 吸附、表面重构能和必要的 AIMD 稳定性结果组合起来，判断该位点在目标条件下是否仍能代表反应表面。
如果论文目标是解释实验活性提升，还应把空位浓度与表征证据相连，例如 XPS 中 S/M 比例变化、EPR 或 Raman 缺陷信号、原位谱学中的 S 相关物种变化。计算端则要说明模型浓度与实验缺陷浓度不必一一相同，但趋势应能解释位点类型、吸附强度和反应路径。S 空位不是单独的活性标签，而是一个同时改变几何配位、局域电子态和表面稳定性的结构变量。
要点：S 空位会暴露低配位金属位，并引入缺陷相关电子态，使 H* 吸附位置和吸附强度发生变化。
要点：HER 位点比较应同时包含完整基面、空位邻近位、边缘位和不同空位浓度，不能只挑一个最优结构。
要点：ΔGH* 需要结合空位形成能、PDOS、差分电荷密度和结构稳定性共同判断。
要点：碱性 HER 还要检查 H₂O 活化和 *OH 稳定性，单独的 H* 热力学并不能覆盖全部反应限制。
