表面能看的是把体相切成表面需要付出的代价。问题在于，氧化物、硫化物、氮化物或含缺陷表面并不一定有相同原子数；一个终止面多 O，另一个终止面少 O，缺陷表面又少一个 O 或多一个 H。若只比较 slab 总能，比较对象就不是同一批原子。统一化学势的目的，是把“多了什么、少了什么”先放到同一个热力学口径里。
这也是表面能公式比吸附能公式更容易出错的地方。吸附能常在同一个干净表面上加一个中间体；表面能则经常要跨终止面、跨缺陷数、跨氧/氢环境比较。只要 slab 组成不同，化学势就不是可选项，而是可比性的前提。
图1 Fe3O4(001) 不同表面结构示意。图中可见不同终止方式对应不同原子组成和表层配位。DOI: 10.1039/d5lf00022j。
以二元氧化物 MxOy 为例，一个 slab 可能恰好满足体相化学计量，也可能因为终止面选择而富 O 或缺 O。表面能中的 μM、μO 就是在说明：多出来的 M 或 O 来自哪里，少掉的 M 或 O 又还给了哪个热力学库。化学势不是给公式增加复杂度，而是在说明表面与环境交换原子的代价。
γ = (Eslab − nMμM − nOμO) / 2A 式(1)
式(1)适合上下两面等价、且形成两个相同表面的 slab。若上下两面不等价，分母不能机械写 2A，而要单独构造对称 slab、配对消元，或说明得到的是两面平均表面能。分母里的面积和表面数量，必须和模型几何一致。
μM 和 μO 也不能随便取孤立原子能。常见做法是让它们满足体相材料平衡条件，例如 xμM + yμO = μMxOy，再在富 O 和贫 O 边界之间扫描。这样每个终止面都被放到同一组环境变量下比较。
图2 Fe3O4(001) 表面能随氧化学势变化的相图示例。图中不同表面结构的相对高低会随 μO 改变。DOI: 10.1039/d5lf00022j。
不同终止面的总能差里混着三件事：表面切割代价、原子数差异、结构弛豫。若 A 终止面比 B 终止面多两个 O，A 的 Eslab 更低或更高都不能直接说明 A 更占优，因为你还没有扣掉两个 O 的化学势。终止面比较必须先扣组成差，再谈表面本身。
ΔγA−B = γA(μO, μH, …) − γB(μO, μH, …) 式(2)
式(2)强调比较发生在同一组 μ 条件下。若在富 O 条件下 A 更低，而在贫 O 条件下 B 更低，就不能只报一个“最优终止面”。更合适的说法是给出环境区间：当 μO 较高时富 O 终止面更有利；当 μO 降低时缺 O 或金属富集终止面可能出现。
对于极性表面，还要额外注意电荷补偿。某些终止面总能低，可能不是因为表面本身有利，而是 slab 通过重构、吸附、缺陷或电荷转移消除了极性发散。极性终止面比较时，组成口径和电荷补偿方式要同时说明。
图3 ZnGa2O4(111) 不同表面终止模型示例。图中展示终止面改变会带来表面原子种类和配位差异。DOI: 10.3390/s25020548。
缺陷表面看似只是在完整表面上去掉一个原子，但热力学意义已经变了。氧空位表面少了一个 O，氢化表面多了一个 H，羟基化表面可能多 H 又改变 O 的连接方式。缺陷表面和完整表面公平比较时，必须把缺失或新增原子的化学势写进同一条公式。
γdef = (Edef slab − Estoich reference + Σnremovedμi − Σnaddedμj) / A 式(3)
式(3)不是唯一写法，但它提醒我们：缺陷能不能作为工作表面，不能只看缺陷表面的吸附能。若去一个 O 后 *OH 吸附很强，却需要极低 μO 才能出现，这个缺陷在富 O 或水环境下未必占优。缺陷吸附强，不等于缺陷表面本身容易出现。
实际计算中，最好把完整终止面、缺陷终止面、羟基化终止面放在同一张 μ 图里，而不是分散给出三个单独能量。这样可以看出哪一种表面身份在目标环境下更有利，再把吸附能或自由能接上去。
图4 ZnGa2O4(111) 表面能与吸附相关能量示例。图中提示不同终止面需要放在统一热力学口径下比较。DOI: 10.3390/s25020548。
气相条件下，μO 常由 O2 的能量和温度、分压修正得到；μH 可由 H2 或 H2O/H2 平衡间接给出。电化学条件下，还要把 pH、电位和 RHE 参照纳入自由能口径。同一个表面，在真空、氧气、水汽和电解液中看到的是不同化学势窗口。
VASP 总能只给出电子能基准，温度和分压项通常来自热力学修正。若只是做趋势筛选，可以先用富 O/贫 O 两端边界；若要和实验气氛对应，就需要把 T、p 或 U、pH 写清楚。没有环境条件，表面能相图只能算模型比较，不能直接等同工作表面。
还要注意 O2 参照的系统误差。PBE 等泛函常对 O2 有偏差，氧化物表面能研究中常用实验形成焓校正、拟合化学势或相对比较来降低误差。关键不是追求一个绝对完美的 μO，而是同一篇计算里所有终止面采用同一套修正。
图5 羟基化 α-Fe2O3(0001) 表面的结构与热力学比较示例。图中说明 H/O 相关化学势会改变表面状态判断。DOI: 10.1021/acsomega.4c02113。
结果部分建议先写清楚三件事：比较了哪些终止面，允许交换哪些原子，化学势窗口怎样设定。然后再给出 γ 或 Δγ 的排序。表面能结论最好写成“在某个 μ 区间内谁更低”，而不是孤立说某个终止面最好。
如果后续还要接吸附能，应先确定目标条件下可能存在的表面，再在这些表面上放中间体。不要把一个在贫 O 条件下才有利的缺陷表面，拿去解释富 O 或强碱水环境下的吸附行为。先判表面身份，再判吸附强弱，这是表面能和反应能衔接的关键顺序。
图6 ACS Nano 中 In2O3(111) 表面随水化学势变化的表面相图。图中展示不同水/羟基覆盖结构的热力学稳定区间。DOI: 10.1021/acsnano.7b06387。
遇到多个表面能差距很小的终止面，不必强行指定唯一结构。可以保留两个候选，再用吸附物覆盖、动力学能垒或实验表征进一步区分。当能量差接近方法误差时，多候选表述比单一排名更可靠。
画表面相图时，还要注意横轴含义。横轴如果是 μO，就只能说明氧环境变化；若同时引入 H2O、H2 或电位，横轴应改写成对应的化学势组合。不要把不同横轴下得到的低能表面放进同一张排序里直接比较。
对于缺陷表面，建议把“形成缺陷”和“缺陷上吸附中间体”拆成两句结论。前一句回答这个缺陷在目标环境下是否容易存在，后一句回答缺陷一旦存在会怎样影响中间体。这两句话合在一起，才是缺陷表面参与反应的完整判断。
如果同一材料有多个终止面都可能出现，可以先按表面能筛掉明显不利的候选，再在剩余候选上做吸附和自由能。这样能避免把大量反应中间体放在不太可能出现的表面上。表面能不是最终活性指标，但它能决定哪些表面值得继续算。
还可以把表面能结果和实验暴露面保持适度对应。若计算显示某终止面只在极窄 μ 区间内较低，而实验长期暴露的是另一类表面，就要谨慎解释。表面能相图给的是热力学候选，不会自动等于样品中面积最大的暴露面。
对催化问题来说，最容易犯的错误是把低表面能直接写成高活性。低 γ 说明表面形成代价小，活性还要看中间体吸附和反应能垒。低表面能负责“能不能作为表面”，吸附自由能负责“能不能作为反应位”。
最后可以用一句话收束：不同终止面和缺陷表面不是比谁的 slab 总能低，而是比在同一化学势窗口下谁的单位面积代价更小；只有这个前提成立，后面的吸附能、自由能和活性位讨论才有公平边界。
要点一：不同终止面和缺陷表面原子数不同，必须用化学势统一组成差。
要点二：表面能分母要与真实表面数量和面积一致，非对称 slab 不能机械写 2A。
要点三：终止面排序应放在同一 μ 窗口下表达，不宜只比较 slab 总能。
要点四：缺陷表面是否可作为工作表面，要同时看缺陷形成代价和吸附行为。
要点五：表面能先确定表面身份，吸附能再判断中间体结合强弱。
