过渡金属磷化物用于析氢反应（hydrogen evolution reaction, HER）时，P 空位常被看作调控活性的关键结构。真正进入 VASP 计算后，问题不只是删去一个 P 原子，而是要判断空位如何改变金属配位、H* 吸附、H₂O 解离和 *OH 位点再生。P 空位不是孤立几何缺陷，而是会同时牵动结构、电子态和反应端点的模型变量。
磷化物表面有哪些基本对象？
过渡金属磷化物通常由金属原子和 P 原子共同构成活性表面。Ni₂P、CoP、CoP₃、NiCoP 等晶相中，金属位负责提供 d 轨道参与吸附，P 位则会影响金属的电荷分布、局域配位和表面极性。HER 中常见的 H*、H₂O*、OH* 并不是随意放在表面任意位置，它们会响应表面的金属-P 配位关系。
因此，讨论 P 空位之前要先把“空位”定义清楚。表层 P 缺位、次表层 P 缺位、界面 P 缺位和高浓度 P 缺位对应的几何环境不同，暴露出的金属位数量也不同。P 空位的第一层作用，是把原本由 P 调节的金属配位环境改写成低配位金属位点。这一步如果没有讲清，后面得到的 H 吸附自由能就很难和具体活性中心对应。
HER 为什么不只看 H*？
HER 的酸性路径常围绕 H⁺ + e⁻ → H* 和 H₂ 生成展开，ΔGH* 是最常用的热力学描述符。碱性条件下，H 通常来自 H₂O，反应还要经历 H₂O 吸附、O-H 断裂和 *OH 处理。也就是说，同一个 P 空位既可能改善 H* 吸附，也可能改变 H₂O 的吸附取向和 *OH 的稳定性。
这也是金属磷化物缺陷模型容易被误读的地方。若只挑一个 H* 最低能构型，得到的只是热力学吸附趋势；若目标是解释碱性 HER，H₂O* 和 OH* 也必须进入同一组模型。P 空位能否提高 HER 活性，取决于 H* 形成、水解离和 *OH 后续处理是否形成同向改善。
图1 p-NiCoP/NCFs 的显微结构和元素分布。图中给出了磷化物纳米颗粒、晶格条纹和 Ni/Co/P 分布，可用于说明含 P 空位模型的材料对象。DOI: 10.1007/s40820-024-01471-9。
空位对应哪个表面？
图1 展示的 Ni-Co-P 磷化物具有纳米颗粒、碳载体和多元素分布，这类材料在计算中不能直接被压缩成一个孤立的理想晶胞。第一步应先确定研究对象：是 Ni₂P、CoP、CoP₃、NiCoP 还是负载在石墨烯或碳纤维上的磷化物界面。不同晶相中的 P 配位不同，删去不同位置的 P，会暴露不同数量的 Ni/Co 金属位。
对 VASP slab 来说，P 空位模型至少要保留三个对照：完整表面、单 P 空位表面和较高空位浓度表面。若研究 NiCoP(200)def/Gr，可在同一超胞里删除一个表面 P，再保持相同真空层、固定层、k 点采样、截断能、POTCAR 和色散修正。吸附能比较需要建立在同组设置上，否则数值差异会混入晶面、超胞和收敛口径的影响。
表征信号怎样进入建模？
图2 中的 XPS 和 EPR 信号说明，P 空位并不是单靠结构图命名的模型。P 2p 结合能变化、金属价态变化和 EPR 缺陷信号可以帮助判断实验材料中是否存在缺陷相关电子态。计算端则需要把这些信号转化为可操作模型：P 缺位、邻近金属低配位、局部氧化层或碳载体接触面。
一个常见对比情境是，实验显示 P 相关信号减弱，但 VASP 中只算完整 CoP 表面再放 H*。这种做法只能回答理想表面的 H 吸附，不能解释缺陷调控。另一个对比情境是，P 空位模型形成后邻近金属发生明显重构，此时应把重构后的最终结构作为吸附基底。P 空位建模要受晶面、缺陷位置和优化后局域重构共同约束。
图2 p-NiCoP/NCFs 的 XRD、XPS 和 EPR 表征。P 2p 与 EPR 信号用于识别缺陷和 P 空位相关特征。DOI: 10.1007/s40820-024-01471-9。
H* 位点矩阵怎样排？
P 空位通常会暴露低配位 Ni 或 Co 位，因此 H* 初始构型不能只放在 P 顶位。建议至少比较空位中心、空位邻近金属顶位、金属-金属桥位、远离空位的完整区域，以及界面金属位。图3 中红圈标出的 P 空位和 H 吸附自由能结果说明，缺陷表面与氧化或界面调控后的表面可能给出不同的 ΔGH*。
H 吸附能的核心是同一基底上吸附前后的能量差，其中 Eslab-H* 为吸附态总能，Eslab 为对应表面总能，EH₂ 为气相 H₂ 总能。若比较完整表面和 P 空位表面，Eslab 必须分别取各自优化后的基底，不能用完整表面能量替代缺陷表面能量。
Eads(H*) = Eslab-H* – Eslab – 1/2EH₂ 式(1)
ΔGH* 接近 0 eV 够吗？
HER 中常把 ΔGH* 接近 0 eV 作为热力学参照，但 P 空位体系要多走一步。自由能校正需要包含 ΔZPE、TΔS、电位项和 pH 项；其中 ΔZPE 和 ΔS 来自频率计算或同一口径文献校正，ΔGU 与 ΔGpH 取决于电位和 pH 标尺。对一组位点比较，ZPE 和熵修正可以先用统一近似筛选，再对代表构型做频率任务。
如果 P 空位让 H* 吸附过强，后续 H₂ 形成可能受到限制；如果 H* 吸附过弱，Volmer 后形成的 H* 覆盖度不足。图3 同时给出差分电荷密度分析，提示 H* 自由能变化需要和电子密度重排一起解释。这里更值得标出的不是某个最低值，而是位点几何、金属 d 态参与和吸附自由能三项证据必须相互支撑。
ΔGH* = ΔEH* + ΔZPE – TΔS + ΔGU + ΔGpH 式(2)
图3 NiCoP(200)def/Gr 的 P 空位模型、H 吸附自由能和差分电荷分析。红圈标出了 P 空位位置，可对应 VASP 缺陷 slab 的位点设定。DOI: 10.1007/s40820-024-01471-9。
碱性 HER 多了哪一步？
酸性 HER 可以通过计算氢电极（computational hydrogen electrode, CHE）把 H⁺ + e⁻ 连接到 1/2H₂；碱性 HER 中，H 来源通常要经过 H₂O 吸附和 O-H 断裂。P 空位邻近金属位可能增强 H₂O 吸附，也可能稳定 *OH。两种效应方向不同，因此只看 ΔGH* 会遗漏 Volmer 步的动力学限制。
图4 的碱性 HER 计算给出了 H* 自由能、H₂O 吸附能和水解离反应能的并列比较。迁移到 P 空位主题时，可以把 H₂O 放在空位邻近金属位，让 O 指向低配位金属，H 指向邻近 P 或金属桥位，再构建 H₂O* → H* + OH* 的端点。若空位只改善 H* 吸附而不降低水解离能垒，碱性 HER 仍可能受 Volmer 步限制。
NEB 端点怎样设置？
水解离路径建议先做端点优化：初态为 H₂O* 吸附在 P 空位邻近金属位，终态为 H* 和 *OH 分别吸附在相邻位点。若终态中 *OH 自动迁移到远离空位的位点，说明单一空位中心可能不能同时容纳 H* 与 *OH。NEB 插值时应保持原子映射一致，并检查中间图像是否出现非目标的分子翻转或表面重构。
过渡态确认需要频率计算。若虚频对应 O-H 断裂方向，说明该路径可用于描述水解离；若虚频主要是水分子转动或表面原子摆动，则需要调整端点或路径分段处理。水解离计算要确认的是 O-H 断裂的反应坐标和能垒，不是给 H₂O 单独寻找一个最低吸附态。
图4 Co@CoFe-P 相关模型在碱性 HER 中的吸附构型、H* 自由能、H₂O 吸附能和水解离反应能。该图说明 H₂O 活化需要与 H* 热力学同时判断。DOI: 10.1039/d1ee02249k。
*OH 会改变位点可用性吗？
碱性条件下，H₂O 解离后留下的 *OH 不能被忽略。若 P 空位邻近金属位同时稳定 H* 和 *OH，就要比较共吸附构型：H* + OH* 在邻位、H* 在空位中心而 OH* 在远端、以及 OH* 覆盖后再次吸附 H₂O 的结构。覆盖度不同会改变表面偶极和电子态，不能直接与孤立 H* 结果混合排序。
图5 的 HER 性能结果提供了实验侧参照，但计算不能把高活性简单归因于一个空位。较完整的闭环应包含 P 空位形成能、H* 自由能、H₂O 解离能垒、*OH 稳定性、PDOS/差分电荷密度和必要的结构稳定性。热力学吸附、动力学水解离和反应后位点再生三者同时成立时，P 空位对 HER 的解释才更完整。
VASP 任务怎样排？
建议把计算对象放回 HER 物理问题中。先优化完整表面和 P 空位表面，用同一套晶面、超胞、真空层、固定层和电子收敛口径确认缺陷基底；随后放入 H*、H₂O* 和 OH*，比较它们在空位邻近金属位、桥位和远离空位区域的构型差异。优化完成后，再选择代表结构做静态计算，输出 DOSCAR、CHGCAR、AECCAR 等结果，用于 PDOS、Bader 电荷和差分电荷密度分析。
结构优化阶段可先采用 PBE 泛函、PAW 赝势、自旋极化和按体系收敛后的截断能；力收敛阈值可用 -0.02 eV/Å 作为起始测试值，电子收敛可从 10⁻⁵ eV 量级开始检查。含 Ni、Co、Fe 的磷化物还要比较不同 MAGMOM 初值，必要时测试 DFT+U 对局域磁矩、PDOS 和吸附自由能的影响。
图5 p-NiCoP/NCFs@CC 的 HER/OER 极化曲线、Tafel 曲线和稳定性结果。性能图为缺陷工程和反应活性关联提供实验参照。DOI: 10.1007/s40820-024-01471-9。
要点：P 空位首先改变磷化物表面的金属配位和局域电子态，H* 位点应围绕空位邻近金属位展开比较。
要点：ΔGH* 接近 0 eV 只是 HER 热力学筛选口径，不能单独覆盖碱性 HER 的 H₂O 活化。
要点：水解离需要独立端点、NEB 路径和虚频确认，尤其要检查 *OH 是否改变位点可用性。
要点：P 空位解释活性提升时，应同时给出空位形成、H* 吸附、水解离能垒和电子结构旁证。
