表面能怎么算？

引言

比较不同晶面时，最容易犯的错是把 slab 总能直接拿来排队。slab 原子数不同、层数不同、表面积不同，甚至上下表面终止也不同，总能当然不能直接比较。表面能γ 比较的是单位面积上切出表面所付出的能量代价。本文围绕γ 的定义、建模和复查，讲清楚为什么低总能不等于低表面能。

表面能为什么要除以面积？

只看 slab 总能会错在哪里？

slab 总能随原子数增加而降低或升高，本质上先反映体系大小。一个 6 层 slab 和一个 9 层 slab 即使来自同一晶面，总能也不能直接比较，更不用说不同晶面往往有不同表面积和原子密度。

表面能把体相部分扣掉，再除以两个表面的总面积。γ 的核心是把“切面成本”归一到单位面积，而不是比较整个模型谁的总能更低。

式(1) γ = (E_slab – Nμ_bulk)/(2A)

公式里的 2A 为什么不是永远都对？

对上下表面对称、终止一致的 slab，常用 γ = (E_slab – Nμ_bulk)/(2A)。这里 2A 表示两个等价表面。如果 slab 上下终止不同，或者一面被钝化，2A 就不再是简单答案。

图1把表面能纳入稳定性筛选，说明它服务的是可形成性和暴露倾向。只有在表面数目、终止方式和面积定义一致时，表面能比较才有物理意义。

图1 催化剂筛选中表面能、稳定性和活性之间的综合比较。DOI: 10.1038/s41467-025-57949-6。

体相参照该怎样取？

μ_bulk 为什么要用同一套计算设置？

体相参照μ_bulk 通常来自充分收敛的 bulk 单胞能量除以原子数。泛函、POTCAR、ENCUT、U 值和磁性设置要与 slab 保持一致，否则 E_slab – Nμ_bulk 会把计算设置差异误当成表面成本。

复杂氧化物、氮化物或含缺陷表面还会牵涉化学势范围。此时不能只拿一个总能公式硬套，而要说明 O-rich、O-poor 或金属富集条件。表面能的参照态一旦变了，不同晶面的排序也可能跟着变。

式(2) μ_bulk = E_bulk/N_bulk

层数收敛为什么会影响γ？

图2显示不同氧化物表面弛豫后的结构差异。表面层数太少时，中间层还没有恢复体相环境，扣除 Nμ_bulk 后剩下的不只是表面成本，还混入了 slab 太薄造成的整体畸变。

实际操作中建议测试 4、5、6、7 层或相近序列，看 γ、关键键长和中间层电荷是否收敛。表面能不是建一个 slab 就能报数，而要先证明层数和真空层足够稳定。

图2 Al₂O₃ 不同表面结构及 DFT 弛豫后稳定性示意。DOI: 10.1038/s41467-026-73690-0。

不同晶面应该怎样放在同一把尺子上？

面积 A 应该从哪里读？

A 来自 slab 面内两个晶格矢量叉乘的面积，不能用屏幕上看起来的大小估计。旋转超胞、重构表面或大超胞比较时，尤其要确认 A 是否对应同一表面周期。

图3的手性表面模型提醒我们，晶面比较常常伴随复杂终止和大超胞。只要超胞大小不同，就必须回到单位面积γ，而不是比较 E_slab 的绝对值。

式(3) A = |a × b|

图3 手性陶瓷表面 slab 与吸附模型示意。DOI: 10.1038/s41467-024-54414-8。

重构和吸附物会不会改变表面能定义？

表面重构会改变 slab 的局域结构，但只要组成不变，仍可按同一化学势计算 γ。若表面多了 O、H、OH 或其他吸附物，就要进入表面自由能或吸附修正，不能继续把它当成干净晶面表面能。

含环境物种时，更稳的是写γ(T,p) 或 ΔG_surf，把外部化学势列出来。干净表面能和吸附覆盖下的表面自由能不是同一个指标。

式(4) ΔG_surf = γA + Σn_iμ_i

VASP里怎样做？

目录和数据列要怎样安排？

建议每个晶面单独建目录，例如 facet_001、facet_011、facet_111。每个目录下再按层数和真空层分子目录，最后统一提取 E_slab、N、A、μ_bulk、γ、最大残余力和偶极修正状态。

不要把不同层数的 CONTCAR 混着继续算，也不要把 bulk 的 INCAR 偷换成更粗的设置。表面能表最重要的不是列很多晶面，而是保证每一列来自同一套参照和收敛标准。

层数收敛最好不要只看γ 的最后一位。若表层键长还在明显变化，中间层没有恢复体相配位，或者真空电势平台随层数漂移，即使γ 数字看着差不多，也要谨慎。对金属表面可以更关注能量和结构，对极性氧化物则要额外看电势和电荷重排。

记录时建议把固定层数也写进表格。固定底层过多会压住表面弛豫，固定太少又可能让整个 slab 漂移。不同晶面若固定规则不一致，γ 的差异就不再只是晶面本身造成的。

式(5) γ₁ – γ₂ = [(E₁ – N₁μ_bulk)/2A₁] – [(E₂ – N₂μ_bulk)/2A₂]

晶面对比应该怎样读？

图4把 TiO₂ 低指数晶面和水吸附行为放在一起，说明晶面稳定性和反应行为是两层问题。低表面能晶面更容易暴露，但未必就是某个反应的最高活性面。

写论文时可以先用γ 说明哪些晶面更可能稳定存在，再用吸附能、反应自由能或能垒解释活性差异。晶面总能、表面能和反应活性要分层讨论，不能用一个总能数字包办。

还有一种常见情况是 slab 弛豫后出现轻微重构。只要结构没有明显塌陷，重构本身不是错误；它可能正是该晶面降低表面能的方式。真正需要警惕的是不同晶面采用了不同弛豫限制，例如有的固定底两层，有的固定底三层，这会让γ 的差异混入人为约束。

若表面含极性终止，最好补充说明是否使用对称 slab、是否做偶极修正、上下两面的终止是否一致。对极性氧化物而言，表面能往往不是一个孤立数字，而是和化学势窗口、表面重构以及补偿机制一起出现。

先说在给定化学势和收敛设置下哪个晶面的γ 更低，再说该晶面是否更可能暴露；随后另起一层讨论吸附或反应。这样即使后续补了 Wulff 构型、覆盖度或溶剂效应，也不会把表面稳定性和反应活性搅在一起。

图4 TiO₂ 低指数晶面及水吸附、解离行为对比。DOI: 10.1038/s41467-023-41865-8。

本章要点总结

要点1：不同 slab 总能不能直接比较，因为原子数、层数、面积和终止方式可能不同。

要点2：表面能γ 是单位面积切出表面的能量代价，常用公式只适用于上下表面对称的干净 slab。

要点3：体相参照、层数、真空层、表面终止和化学势条件都会影响γ 的数值。

要点4：低表面能说明晶面更可能稳定暴露，不等于该晶面对所有反应都更活泼。