dft计算表面能

晶体暴露什么晶面，直接决定它的催化活性、溶解速率甚至形貌——而表面能就是晶面”暴露成本”的量化指标。我这个项目客户在合成一种金属氧化物纳米晶，想控制它暴露高活性晶面来做催化剂载体，让我用 dft计算表面能 把各个 (hkl) 面的能量算出来，再用 Wulff 构型预测实际会长成什么形状。

计算方法本身是标准的 slab 模型。对每一个目标 (hkl) 晶面，我切出满足对称性的对称 slab（symmetric slab，上下表面相同，避免偶极矩和表面应力不对称带来的误差），厚度取 8–12 层、真空层 ≥ 15 Å。表面能公式 γ = (E_slab − N × E_bulk_per_atom) / (2A)，其中 E_slab 是 slab 总能、N 是 slab 原子数、E_bulk_per_atom 是体相单原子能、A 是表面积，乘 2 是因为上下两个表面。收敛上，截断能 500 eV、k 点沿表面两个周期方向取到 6×6 以上，力收敛 < 0.02 eV/Å。

我这个项目我算了 (100)、(110)、(111)、(211) 四个主要晶面，结果排成柱状图：低指数密排面 (111) 表面能最低（约 1.1 J/m²），(110) 次之，(100) 和 (211) 明显更高（约 1.6–1.9 J/m²）。判读逻辑很直接：表面能越低，这个晶面越”愿意”暴露，因为形成它的能量代价小。也就是说，热力学平衡下晶体倾向于被低表面能的晶面包裹。

Wulff 构型是这一步的输出。根据 Wulff 定理，晶体平衡形貌里每个晶面的法向距离正比于它的表面能 γ_(hkl)。我把算出的各面 γ 值按 Wulff 构造法则画出来：低 γ 的 (111) 面占据大部分外表面积，(100)/(211) 这些高 γ 面被”挤压”成小 facet 甚至消失。预测出的形貌是一个以 (111) 为主的截角多面体——这和客户电镜里看到的颗粒形状对得上。更有意思的是，客户想提高催化活性，而高活性往往来自高表面能的 (100)/(110) 面，所以我进一步算了”用表面活性剂封端低能面、迫使高能面暴露”的热力学代价，给了一个形貌调控的定量参考。

踩坑点：第一，slab 对称性。非对称 slab 会引入偶极矩，表面能算出来偏差明显，必须用能抵消偶极的对称 slab 或加偶极修正。第二，厚度收敛。厚度不够时上下表面相互作用会让 γ 偏低，我做了 8/10/12 层三档确认收敛。第三，表面重构。有些高能面（如 (100)）松弛后会发生表面原子重排，不弛豫直接算会高估表面能。第四，极性面。对离子晶体，某些 (hkl) 是极性面、带电，必须做非极化处理（如加一层补偿或取合适终止），否则表面能发散。

从我的工程经验看，dft计算表面能 的真正用处不是列一张表，而是把”形貌-晶面-性能”这条链路打通。很多合成以为靠试表面活性剂就能控形貌，其实底子是表面能排好的序。我给客户的建议是：想暴露高能活性面，得用封端剂把低能面先”盖住”，否则热力学自己会把高能面吞掉。更多表面能计算，见 [dft计算表面能](https://www.keyanxueshu.com/category/dft/)；需要定制，[dft计算表面能](https://www.keyanxueshu.com/) 上有入口。