氧空位形成能不是一个只由缺陷结构决定的固定数。氧从表面离开以后进入哪种氧库，周围是富氧还是贫氧，温度和 pO₂ 如何设定，都会改变形成能的热力学含义。氧空位形成能本质上是缺陷结构与氧化学势共同定义的反应能，因此同一个 VASP 总能差不能脱离气氛边界解释。
被移走的氧到哪里去？
最常见写法是从完整表面移走一个 O，得到含 VO 的表面，同时把被移走的 O 放入氧库。若氧库取 1/2 O₂，形成能可写成 Ef(VO) = Eslab-VO + μO – Eslab。这里的 μO 不是任意常数，而是由气氛和材料稳定区间限定。
在目录设计上，完整表面和缺陷表面应使用同一个 POSCAR 祖结构编号，便于追踪被移走的氧位。若一个表面有桥氧、端位氧和羟基氧，三类氧位要分别命名，否则后续很难判断能量差对应哪一种缺陷。
Ef(VO) = Eslab-VO + μO – Eslab 式(1)
式(1)中 Eslab 和 Eslab-VO 应来自同一表面、同一层厚、同一固定方式和同一电子设置。若完整表面与缺陷表面优化到不同磁态或不同重构形态，形成能会同时包含缺陷与结构态变化，需单独说明。
图1 对比了多种金属氧化物的氧空位形成自由能，并展示不同参照会给出不同数值。参照态决定被移走氧原子的热力学身份。热催化常用 O₂ 和温压修正，电化学环境还可能用 H₂O/H₂ 或 CHE 口径表达。
图1 多种金属氧化物在标准条件及不同参照下的氧空位形成自由能。该图用于说明氧空位形成能依赖参照反应式和环境项。 DOI: 10.1021/acsomega.6c03286。
O-rich/O-poor 怎样进入？
O-rich 条件通常把 μO 取接近 1/2 EO₂ 加温压校正的上限；O-poor 条件则受金属析出或低价氧化物形成边界限制。对于 ABO₃ 钙钛矿或尖晶石氧化物，A、B 金属化学势也要满足母相稳定条件，不能只单独移动 μO。
如果只做表面趋势比较，可以先给出两个边界值，再在表格中标出每个位点形成能随 μO 的线性移动。若要和实验气氛对应，则需要把温度、氧分压和相稳定区间一起列出。
μO(T,p) = 1/2EO₂ + ΔμO(T,p°) + 1/2kBT ln(pO₂/p°) 式(2)
式(2)说明温度升高或 pO₂ 降低时，μO 会下降，氧空位形成能随之改变。若论文只报告一个 0 K 总能差，却把它解释为高温还原气氛下的空位稳定性，热力学条件就没有对齐。
图6 中不同界面空位构型随氧化学势改变而移动，适合读出 O-rich/O-poor 的实际作用。贫氧条件降低的是氧库化学势项，不是自动改变所有缺陷结构的相对总能。不同位点之间的排序仍要看缺陷 slab 的结构能和局域弛豫。
图6 ZrO₂ 与 In₂O₃ 界面中不同氧空位构型随氧化学势变化的形成能。该图对应 O-rich/O-poor 条件下空位稳定性的读图方式。 DOI: 10.1038/ncomms11892。
表面空位有什么额外变量？
表面氧空位比体相空位多出终止面、配位数、羟基化、吸附物覆盖和偶极修正等变量。桥氧、顶层晶格氧、羟基氧和次表层氧的局域环境不同，移走一个氧以后周围金属价态、磁矩和键长弛豫也不同。
若同一氧化物有多个终止面，应先分别建立完整 slab，再在每个终止面内移走等价或非等价氧位。图4 显示表面氧空位形成自由能与体相氧化物形成自由能存在趋势关联，但不能把体相形成能直接替代表面位点能。
图4 多种氧化物表面氧空位形成自由能与体相氧化物形成自由能的关联。该图用于说明表面空位和体相稳定性不能直接混用。 DOI: 10.1002/aenm.201401082。
表面和次表层怎样比较？
次表层氧空位会引起更接近体相的局域弛豫，表面氧空位则会直接改变吸附位和界面电荷。图5 对比了表面与次表层氧空位形成自由能，说明空位深度本身就是一个计算变量。
如果表面 VO 更低，缺陷可能直接参与 O₂、H₂O、CO 或 *OH 的吸附；如果次表层 VO 更低，则缺陷对催化的影响可能通过电子结构或晶格氧迁移间接体现。空位深度与表面终止面应在同一组模型中分开比较。
对 CeO₂、TiO₂、Co₃O₄ 或 NiFeOOH 一类可变价体系，还要检查缺陷电子是否局域在邻近金属上。若 Bader 电荷、局域磁矩和 PDOS 都显示还原态形成，空位形成能才和局域价态变化对应起来。
图5 次表层氧空位形成自由能与表面氧空位形成自由能的对照。该图用于比较不同氧位点形成空位的热力学差异。 DOI: 10.1002/aenm.201401082。
电位项进入哪里？
在 OER 或 ORR 条件下，表面移走氧不一定写成产生 1/2 O₂。若空位形成伴随 H₂O、OH⁻、H⁺ 和 e⁻，反应式需要按照 CHE 或碱性参照重新配平。电位 U 和 pH 项进入反应自由能，而不是简单改动 VASP 总能。
同一个氧位在酸性和碱性写法下可能对应不同的粒子交换项。比较不同材料时，应先统一电极电位标尺，再统一 H₂O、H₂ 或 OH⁻ 参照，最后再讨论表面总能差。
ΔG(U,pH) = ΔEDFT + ΔZPE – TΔS – n e U + m kBT ln(10) pH 式(3)
式(3)中 n 是转移电子数，m 是质子或 OH⁻ 参照带来的 pH 系数。若比较不同电位下的 VO 稳定性，应明确每个反应式中的电子数和质子数。图2 显示 URHE 改变后，金属氧化物表面空位形成自由能会整体移动，并可能改变材料之间的相对趋势。
若没有配平反应式，直接把所有空位形成能都加同一个电位修正，会让不同氧位、羟基位和晶格氧位的比较失去化学意义。在这种比较中，电化学空位形成能要先写反应式，再代入电位和 pH 项，再讨论缺陷稳定性。
图2 金属氧化物在不同 URHE 下的氧空位形成自由能。该图对应电位项进入表面空位稳定性的方式。 DOI: 10.1021/acsomega.6c03286。
羟基位和晶格氧能等同吗？
NiOOH、CoOOH 或 FeOOH 表面常同时存在晶格氧和羟基氧。移走晶格 O 与移走 OH 对应不同缺陷，前者可能形成 VO，后者可能形成羟基空位或脱氢态。图3 中氧空位和羟基空位被分开标出，正是因为二者的反应式和自由能项不同。
若 OER 工作态已经羟基化或脱氢，完整表面的定义要先固定。以 NiFeOOH 为例，含 OH、O 或 OOH 的覆盖态不能和干净氧化物表面直接混用；否则空位形成能会同时包含脱氢、补水和吸附覆盖变化。
羟基空位不能简单当作晶格氧空位。VASP 目录中应分开计算完整羟基化表面、去 H 表面、去 OH 表面和去晶格 O 表面，并用对应的 H₂O/H₂/O₂ 参照写出自由能。
图3 3d 过渡金属羟基氧化物表面中氧空位与羟基空位模型及形成自由能。该图用于区分晶格氧位和羟基位缺陷。 DOI: 10.1021/acsomega.6c03286。
哪些结构先算？
第一组结构是完整 slab，包括不同终止面和不同羟基化程度；第二组结构是在每个非等价 O 位点生成 VO；第三组结构是表面、次表层和界面空位；第四组结构是缺陷后吸附 O₂、H₂O、*OH 或 CO 的端点。这样可以把空位稳定性和后续反应活性分开。
参数上，氧化物通常需要检查 ISPIN、MAGMOM、LASPH 和 DFT+U。U 值、初始磁矩和空位邻近金属的价态会影响缺陷电子局域。所有参数只是参考起点，实际数值应按体系、赝势和目标性质做收敛测试。
含空位 slab 的优化还要观察邻近金属是否发生过大的位移。若移走顶层 O 后，相邻金属进入新的配位环境，这个结果代表缺陷弛豫后的稳定结构；若整个表面层发生明显翻转或吸附覆盖同时改变，形成能应和原先的单空位问题分开归类。缺陷身份由移走的氧位和优化后的局域结构共同决定。
若完整表面是绝缘体而缺陷表面出现金属性，ISMEAR、SIGMA 和能量读取口径也要统一。最终总能应来自最后一个完成的离子步或静态步对应的能量项，不能混用不同输出项。
对含水表面，还要把吸附水、羟基化和空位三者分开设定，避免把水解离后的结构误归为单个 VO。
位点、气氛和电子态怎样放在一起？
同一张缺陷表中可保留终止面、氧位编号、配位数、空位深度、Eslab、Eslab-VO、μO、O-rich/O-poor 边界、Ef(VO)、邻近金属磁矩、Bader 电荷、PDOS 缺陷态和吸附后自由能。
若某个表面空位形成能较低，但缺陷态深埋带隙且 O₂ 或 H₂O 吸附很弱，它不一定直接提升目标反应；若形成能较高但一旦形成就显著改变 *OH 或 *OOH 自由能，则它更像反应条件下少量活性缺陷。
热催化和电催化结果也要分开读。热催化更关注温度、pO₂ 和还原气氛下 VO 是否热力学可达；电催化更关注 U、pH 和表面覆盖物共同定义的工作态。若同一材料在贫氧高温下容易形成 VO，并不等价于在 OER 电位下也能维持同一种缺陷表面。
最终报告中，氧空位形成能应和气氛、位点和后续吸附端点一起报告，只给一个空位总能差，无法说明该缺陷是在富氧、贫氧、高温、还原气氛还是电化学电位下更有热力学依据。
要点：氧空位形成能由缺陷 slab 总能和氧化学势共同决定。
要点：O-rich/O-poor 条件改变 μO，不同位点排序仍取决于缺陷结构能和局域弛豫。
要点：表面氧、次表层氧、羟基氧和界面氧应分开建模并写出各自反应式。
要点：电化学条件下需先配平 H₂O、H⁺/OH⁻ 和 e⁻，再加入 U 与 pH 修正。
