为什么同一材料切成不同晶面后 γ 会不同？

晶面稳定性差异来自表层原子的局域环境。密排面保留更多近邻，断键少，表面弛豫幅度通常较小；开放面、高指数面和极性面会露出更多低配位原子、台阶或带净偶极的层状结构。

γ 反映的是这些表面状态相对体相多出的单位面积能量；同一种 MnS，(100)、(010)、(111)-S、(311)-S 对应的外层配位和化学计量都不同。

Miller 指数只给出切面方向，真正进入计算的是最外层元素、面内面积、上下两侧是否等价、表层弛豫后的结构。

例如同一 ZB-MnS 晶面可以有 Mn 终止或 S 终止；某些终止在 Mn-rich 条件下低能，换到 S-rich 条件后排序会改变。把“(111) 面”写成一个固定 γ，容易漏掉终止和环境变量。

在极性或含过渡金属 d 电子的表面上，泛函误差会被低配位和局域态放大。RS-MnS 的非极性 (010)/(100) 面在 r2SCAN 与 HSE06 之间差异较小；

S 终止的极性面含有 S 3p dangling-bond states，并与次表层 Mn 3d 态耦合，r2SCAN 会低估这些表面的 γ。加入 U = 2.7 eV 后，(111)-S 与 (131)-S 更接近 HSE06 基准。

图2. RS-MnS 多个晶面的 r2SCAN、HSE06 和 r2SCAN+U 表面能对比。

由断键数走到电子局域后，不同晶面稳定性不能只归因于断键数。断键数、表面配位和原子面密度给出结构起点，电子局域、磁态、极性偶极和弛豫方式决定最终能量。

对过渡金属硫化物、氧化物和含强关联 d 态材料，泛函和 U 值会参与定义可比较的 γ，而不只是影响小数点后的误差。

图3. ZB-MnS 对称与非对称 slab 模型及不同终止的稳定区间。

ZB-MnS 的终止比较把问题压缩到对称 slab、非对称 slab 和终止选择上：(100)、(111)、(210)、(221) 等晶面需要在 S 终止与 Mn 终止之间选择。

Δγ 大于 0 区域代表 S 终止更低能，Δγ 小于 0 区域代表 Mn 终止更低能；(221) 和 (210) 的稳定终止会在 ΔμS = −1.62 eV 附近切换。

用表面能判断稳定晶面时要绑定哪些条件？

用 γ 比较晶面时，先把能量基准锁住：同一体相相、同一磁态、同一泛函、同一 U 值或无 U 设置，相近 k 点密度、相近截断能、足够厚的 slab 和足够宽的真空层。

对极性 slab，需说明偶极修正、对称化方案或 wedge/相对能方法。未收敛的 slab 厚度和不一致的 bulk 参考会直接变成虚假的晶面差异。

清洁真空表面得到的是 0 K 静态结构的 γ，它不能自动外推到带配体、带溶剂或带电位的实验表面；合成溶剂、气氛、pH、电位和吸附覆盖物会改变表面自由能。

若实验里看到某高能面大面积暴露，它可能来自动力学生长、配体选择性吸附、重构或尺寸效应，不一定推翻清洁表面能排序。

把 γ 用于形貌、吸附或催化问题时，应先区分它能回答的对象：表面能回答某晶面在给定条件下是否容易出现在平衡形貌里；吸附能、反应自由能、NEB 能垒、PDOS、Bader 电荷和差分电荷密度才用于解释位点反应。

一个低 γ 晶面可能平整且稳定，但低配位活性位少；一个高指数面可能面积小，却在边角或台阶处提供更强吸附。

图8. RS-MnS 纳米晶 TEM 形貌和表观表面能拟合。

RS-MnS 的 TEM 与量热结果给出实验侧参照条件：颗粒从小尺寸、不规则形貌逐步变成立方块，表观表面能拟合值 1.15 ± 0.38 J·m-2 高于理论平衡立方体的 {100} 值。

原因不只来自 DFT 误差，还包括高指数面暴露、表面积估计、油胺包覆、非理想表面构型和小颗粒中的边角占比。

写表面能结果时，最少应同时列出材料相、(hkl)、终止、化学势范围、slab 层数、真空层、弛豫层、泛函、U 值、单位和面积定义。

γ 的单位也要统一，常见 J·m-2 与 eV·Å-2 之间相差 16.02 倍。

只有这些条件对齐后，低 γ 才能解释平衡形貌里的大面积暴露，高 γ 才能指向台阶、极性、重构或环境稳定化带来的额外能量代价。

表面能计算的有效结论应停在具体模型上：RS-MnS 中 (100)/(010) 低能面支持纳米立方体；ZB-MnS 的 {110} 与 {111}-S 竞争随 ΔμS 改变形貌；

WZ-MnS 的 polar facet 需要用相对表面能处理端面稳定。离开这些模型条件，单个 γ 数值不能直接替代实验形貌、工作态表面或催化活性判断。