同一个晶体切出两个 slab，一个终止面总能低，另一个终止面总能高，能不能直接说前者更稳定？对化学组成完全相同、上下表面对称且面积一致的模型，有时可以初步比较；但只要终止面组成不同、表面吸附了 H/O/Cl/OH，或外界氧分压、电位、pH 发生变化，总能差就必须转化为表面巨正则势。它可以理解为从 slab 总能中扣除外界原子库贡献后得到的单位面积稳定性指标。slab 稳定性不是一个固定总能排序，而是随化学势窗口变化的表面相图问题。
总能里混入了什么？
slab 总能同时包含体相原子数、表面多余能量、吸附物贡献、缺陷贡献和弛豫能。若两个 slab 的原子数完全相同，直接比较总能差还有清晰意义；若一个终止面富 A、另一个终止面富 B，总能差会把多出来的 A 或 B 原子也算进去。此时必须引入化学势，把不同组成模型放到同一个原子库中比较。
γ = [Eslab – ΣᵢNᵢμᵢ] / A 式(1)
式(1)里 γ 是表面自由能或按具体模型定义的表面巨正则势，Eslab 是 slab 总能，Nᵢ 是第 i 类原子数，μᵢ 是对应化学势，A 是表面积。若上下两个表面等价，公式中常出现 2A；若上下表面不等价或采用非对称 slab，应按模型定义单独处理。这里的化学势用于把不同组成终止面换算到同一个外界原子库。
稳定性窗口从哪里来？
图1 给出 InP(100) 清洁表面随 In 化学势变化的相图。横轴不是随意参数，而是 In-rich 到 P-rich 条件下允许的化学势范围。不同重构线条的交点代表最稳定终止面的切换。也就是说，一个终止面可能在 In-rich 条件下稳定，在 P-rich 条件下被另一个重构取代。
这类图比单点表面能更适合指导 slab 选择。若研究的是富金属合成环境，就应优先考虑富金属窗口内稳定的终止面；若研究的是氧化、氯化或氢覆盖环境，就不能只取真空清洁表面的最低能结构。终止面稳定性应在目标化学势窗口内判断，而不是脱离环境给出全局排序。
图1 InP(100) 清洁表面随 In 化学势变化的稳定性相图。不同重构在不同 μIn 区间内成为稳定表面。DOI: 10.1021/acsami.4c20370。
结构枚举要覆盖哪些对象？
表面相图的可靠性取决于候选 slab 是否覆盖了可能出现的终止面和重构。图2 展示多种 InP(100) 表面重构，它们不是同一个 slab 稍微移动几个原子，而是具有不同表面原子数、配位方式和重构单元。若候选集合缺少某类低能重构，相图会错误地把次优结构判为稳定。
VASP 建模时，应先固定体相晶胞和晶面，再生成一组候选终止面：富金属、富非金属、缺原子、加原子、重构、吸附覆盖。每个候选都要用相同泛函、POTCAR、ENCUT、k 点密度、真空层和偶极处理。只算一个理想截断面，得到的是该模型的表面能；枚举多个终止面并加入化学势，才能判断哪些结构会进入稳定性窗口。
图2 InP(100) 的多种表面重构模型。该图说明终止面稳定性判断依赖候选 slab 集合，而不是单个理想截断面。DOI: 10.1021/acsami.4c20370。
化学势范围怎样限定？
化学势不能任意取值。对二元 AB 化合物，μA 与 μB 要满足体相 AB 的形成关系，同时还要避免析出单质 A、单质 B 或竞争相。常用做法是先选单质 A、单质 B 或稳定气相分子作为零点，写成 ΔμA、ΔμB；在不析出单质的条件下，Δμ 通常不能为正。
μA + μB = μAB,bulk 式(2)
式(2)只是最简单的二元化合物边界。若写成相对化学势，就是 ΔμA + ΔμB = ΔHf(AB)。A-rich 端通常取 ΔμA = 0，由式(2)推出 ΔμB；B-rich 端则反过来。这样得到的是主相本身允许存在的第一层范围。
mμA + nμB ≤ μAₘBₙ,bulk 式(3)
竞争相会继续收窄化学势窗口。若数据库或补充计算显示 A₂B、AB₂、氧化物、氢氧化物等相更容易形成，就要加入式(3)这类不等式；否则某个终止面看似稳定，实际对应的化学势点已经让体相分解。对多元氧化物、卤化物或含碳表面，这一步往往比单个 slab 总能更决定相图边界。
外界可变物种要单独说明参照。氧化物表面的 μO 可由 O₂、温度和氧分压换算，但 O₂ 参照常需要校正；含 H、Cl、OH 的表面，要说明使用 H₂、Cl₂、H₂O、HCl 或电化学氢电极等哪一种口径。若讨论电位和 pH，还应把 μH、μOH 与电化学条件关联起来，而不是把它们当成任意滑动参数。
H、Cl、O、OH 是新变量吗？
一旦表面吸附 H、Cl、O、OH 或水分子，表面稳定性就从一维化学势问题扩展为多变量问题。图3 给出 InP(100) 表面在 H 和 Cl 化学势下的稳定区切换。固定 μIn 后，H/Cl 覆盖会让原本稳定的清洁终止面被共吸附结构取代。
这对催化计算很重要。真空中最低表面能的 clean slab，不一定是反应环境中的工作表面。若 OER 条件下表面会羟基化，或者卤化物钙钛矿表面会被卤素/小分子覆盖，吸附自由能就应在对应覆盖态基底上重新计算。吸附物化学势可以改变稳定 slab 本身，从而改写后续吸附位点和反应路径。
图3 InP(100) 表面在 H 与 Cl 化学势下的稳定性相图。图中展示吸附物化学势如何改变稳定重构和覆盖态。DOI: 10.1021/acsami.4c20370。
电位和氧化环境怎样进入？
图4 展示金红石 IrO₂、RuO₂ 和 TiO₂ 不同低指数晶面的表面自由能随氧化环境变化。线条斜率与终止面中氧含量有关，折点通常对应最稳定终止面的切换。把这个逻辑移到电催化时，电位、pH 和氧化态会共同影响 *O、*OH、*OOH 或表面氧覆盖。
若研究 OER、ORR 或氧化物表面重构，清洁表面能只能作为入口。更接近工作条件的做法是构建含 O/OH/H₂O 的覆盖态，计算对应表面自由能，再决定哪个 slab 作为反应基底。相反，如果只是比较真空中不同晶面暴露趋势，加入电位项就会把问题复杂化。表面稳定性模型应服务于研究环境，而不是把所有变量一次性加入。
图4 IrO₂、RuO₂ 和 TiO₂ 低指数晶面的表面自由能随氧化环境变化。线条折点显示最稳定终止面会随氧化条件切换。DOI: 10.1002/cssc.202200015。
候选终止面太多怎么办？
多元碳化物、氮化物、高熵合金和缺陷氧化物会带来大量表面组成和原子排布。图5 展示 Fe-C 复杂表面在不同碳化学势下的稳定终止面比较，说明多元体系中稳定表面可能随 μC 连续切换。此时手工枚举少量终止面往往不够，需要结构筛选、主动学习或描述符辅助。
对常规 VASP 任务，可以先建立小规模候选集：低指数晶面、富金属/富非金属终止、一个或两个代表性缺陷、若干覆盖态。若候选空间很大，再考虑用 cluster expansion、机器学习势或主动学习筛选候选，最后用 DFT 精算稳定区域。复杂表面的关键不是穷尽所有模型，而是让候选集合覆盖可能改变稳定性的组成变量。
图5 Fe-C 多元表面稳定终止面随碳化学势变化的比较。该图说明复杂多元素表面需要在化学势维度上筛选稳定结构。DOI: 10.1016/j.xinn.2024.100571。
结果怎样落到 slab 选择？
稳定性窗口最终要回答一个实际问题：后续吸附、缺陷或反应计算应该用哪个 slab。若某终止面只在极窄化学势范围内稳定，且实验或工作条件很难达到，后续可以把它作为对比而不是主模型。若两个终止面在目标窗口内自由能接近，则应同时计算关键吸附态，因为真实表面可能存在混合暴露或局域切换。
同一个表面也可能在 clean 条件下稳定，却在 H/OH 覆盖下被另一表面取代。此时不能把 clean slab 上的吸附能作为工作条件结论，而应重建覆盖态 slab 并重新放置反应中间体。slab 稳定性窗口决定模型入口，吸附自由能和反应能垒是在这个入口之后才有可比性。
最小数据组有哪些？
每个候选 slab 至少需要记录晶面、终止面、表面积、上下表面是否等价、原子数 Nᵢ、总能 Eslab、对应 μᵢ 变量和是否含吸附物。若用非对称 slab，应说明偶极修正、真空层和平面平均静电势是否稳定。若表面含缺陷或吸附物，缺陷浓度和覆盖度也应写入模型身份。
计算顺序可以这样安排：先优化体相并确定允许化学势窗口；再生成候选 slab，完成层厚和真空层检查；随后计算每个候选的 Eslab；最后把Eslab – ΣNᵢμᵢ除以面积，画出随化学势变化的线图或相图。表面相图中的每条线都应能追溯到一个具体 slab 模型和一组清楚的原子计数。
哪些比较不能直接做？
不同面积、不同原子数、不同覆盖度、不同上下表面处理的 slab 不能只看总能。不同晶面之间若要比较暴露趋势，应统一面积归一和化学势口径；同一晶面的不同终止面若组成不同，也要进入化学势表达式；含吸附物与不含吸附物的表面，需要加入吸附物化学势或固定覆盖度口径。
最终表述应限定在目标窗口内。例如“在 O-rich 条件下终止面 A 的表面自由能最低”，“在给定 H/Cl 化学势组合下共吸附结构 B 稳定”，“在目标电位附近 clean slab 不是最低自由能表面”。这些判断比“某 slab 总能最低”更能服务后续 VASP 吸附与反应计算。
要点：不同组成的终止面不能只比较 slab 总能，必须用化学势把多余或缺少的原子纳入同一参照。
要点：稳定终止面会随 μA、μO、吸附物化学势、电位或 pH 等环境变量切换。
要点：表面相图依赖候选 slab 集合，遗漏重构、缺陷或覆盖态会改变稳定性判断。
要点：稳定性窗口决定后续吸附和反应计算的模型入口，不应与吸附能公式混用。
