说明：本文华算科技主要介绍表面能在 DFT slab 模型中的计算口径、不同晶面 γ 值产生差异的结构原因，以及用 γ 判断晶面稳定性时需要绑定的终止、化学势、泛函和形貌条件。

表面能 γ 的计算对象不是一个孤立原子，也不是把晶面名称代入公式后得到的标签值。它比较的是含外露面的 slab 总能和同数量体相结构单元的参考能，差值再按外露面积归一。

对上下两个表面等价的模型，常见写法为 γ = (E_slab − N E_bulk) / 2A，其中 A 是单面面积，分母 2 来自上下两个相同表面。

这条式子能用在清洁、对称、化学计量不变的 slab 上。只要上下两侧终止不同，或者表面少了某种元素、吸附了 O、S、H、OH 等物种，能量差就不能简单平均。

此时需要引入元素化学势 μ_i，把 E_slab 与 Σn_iμ_i 比较；在 MnS 这类体系里，硫化学势 Δμ_S 会直接进入极性终止面的 γ。

对称 slab 求单个表面，非对称 slab 往往只能先得到两个表面的和。若晶体沿极性方向切开，两个外表面分别露出不同元素或不同配位层，γ_top 和 γ_bottom 被 E_slab 捆在一起。

要拆出某一个极性面的表面能，常需要把 slab 结果与 wedge 模型、固定层/弛豫层组合使用。

图1 中的 RS-MnS(100) 对称模型，两个表面相同，E_slab 与体相参考相减后除以 2A 即可。WZ-MnS(0001)/(0001) 这类非对称模型只能给出两侧表面能之和；

ZB-MnS 与 WZ-MnS 的 polar facet 需要借助 wedge 方程，把表面面积、固定层、弛豫面和化学计量差逐项分开。

表面能计算的第一层核对不是数值大小，而是模型是否对应同一个物理量：

E_bulk 是否来自同一泛函和同一磁态，slab 中间层是否恢复体相配位，真空层是否隔开周期镜像，固定层是否改变表面弛豫。

任何一项不一致，γ 的差值都会混入模型误差。

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

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

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

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

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

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

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

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

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

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

多组分材料的表面能常写成 γ(Δμ)。在 Mn-rich 条件下，含 Mn 较多或少 S 的表面可能更低能；在 S-rich 条件下，S 终止极性面会被化学势项拉低。

稳定晶面排序随 Δμ 变化，来源于比较对象从一个环境窗口移动到另一个环境窗口。

RS-MnS 具有较高对称性，(100) 与 (010) 面在 r²SCAN+U 下长期保持最低 γ。

即使把 Miller 指数最高到 3 的高指数面纳入，Wulff 构型仍以纳米立方体为主；只有接近 S-rich 端点时，S 终止的 {311} 面才略微截去立方体边角。

这里的稳定性指给定硫化学势窗口内的平衡外露倾向。(100)/(010) 低能面面积大，Wulff 形貌中自然成为主外露面；

(131)-S 或 (311)-S 即便随 Δμ_S 降低，也要低到足以进入 Wulff 包络，才会在形貌里占据可见面积。

从 Mn-rich 移到 S-rich 的窗口内，ZB-MnS 的低能面集合会直接替换：{110} 面在 Mn-rich 时控制菱形十二面体；

Δμ_S 升高后，{111}-S 面的 γ 持续下降，S-rich 端出现由 4 个三角 {111} 面和 12 个三角 {110} 面围成的 16 面体。这里改变的不是同一形貌上的颜色，而是低能面集合本身发生替换。

WZ-MnS 沿 [0001] 方向缺少反演和镜面对称，某些 polar facet 的单独 γ 难以直接拆出，只能用 Δγ 描述某个端面相对参考面的高低。

终止层外层是 Mn²、Mn³ 还是 S¹/S²，会改变能量曲线斜率，也会改变纳米棒端部是否截断。

WZ-MnS 的 Wulff 序列显示，[0001] 顶部形貌对 Δμ_S 不太敏感，而 [0001] 端部会随硫化学势出现截断差异。

同一材料、同一晶相、不同极性方向都可能给出不同稳定面，这就是晶面稳定性不能脱离终止和环境讨论的原因。

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

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

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

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

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

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

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