DFT 计算表面能

本项目以某晶体为模型，用 DFT 计算了 (100)、(110)、(111)、(211) 四个低指数晶面的表面能，判断哪个晶面在热力学上最稳定。这张图用柱状图给出各晶面表面能，最低者即为最易暴露的稳定晶面。

这张图一眼能看出什么

四个柱子高低不同，比的就是哪个晶面”更省能量”。柱越矮，说明切开晶体长出这个面越容易——也就是热力学上越稳定。本案例里 (111) 面最矮（1.38 J/m²），(211) 最高（2.05）。所以不用算太多，看图就能先猜：这个晶体平衡时最可能露出来的就是 (111) 面。

计算逻辑

表面能定义为：

γ = (E_{slab} − N·E_{bulk}) / (2A)

其中 E_{slab} 是 slab 模型总能量，E_{bulk} 是单个原子的体相能量，N 是 slab 中原子数，A 是单面表面积，因子 2 是因为 slab 有上下两个表面。

本项目流程：

1. 优化体相结构，得到 E_{bulk}。
2. 构建各晶面的 slab 模型，真空层 ≥ 15 Å，保证上下表面不相互作用。
3. 固定底部若干层原子，弛豫表层，计算 E_{slab}。
4. 代入公式得到各晶面 γ。

关键参数

• Slab 厚度：≥ 5 层原子，否则中心未恢复体相环境会引入误差。
• 真空层：≥ 15 Å，避免周期性镜像相互作用。
• 底层固定：通常固定底部 2 层，让表层自由弛豫。
• k 网格：对低对称 slab 要适当加密，避免积分误差。

Slab 模型怎么切（以 (111) 面为例）

表面能的可靠性 80% 取决于 slab 构建。本项目切 (111) 面的实操步骤：

5. 取对称切法：沿 (111) 法向切出 slab，使上下两个表面化学环境相同，避免人为偶极矩。
6. 层数收敛测试：分别算 5、7、9 层，比较表面能变化，本案例 7 层后波动 < 0.02 J/m²，定为最终厚度。
7. 真空层定向：真空层放在 c 轴（法向），厚度 ≥ 15 Å，让上下表面互不干扰。
8. 偶极修正：若上下表面不等价（极性面），必须开 `LDIPOL = .TRUE.` 并设 `IDIPOL` 方向，否则表面能会被虚假偶极拉偏。

这个案例能借鉴什么

• 表面能的绝对值对 slab 厚度敏感，比较不同晶面时用相同厚度才公平。
• 极性表面需要额外处理（如表面钝化、对称 slab），否则会出现虚假偶极矩。
• 最稳定晶面不一定就是催化活性最好的晶面，活性位点往往出现在台阶、棱或高指数面上。

表面能能用来干什么

算完四个晶面 γ 值，本项目把结果推进到两个实际应用：

• Wulff 构型预测：以各晶面 γ 为权重，用 Wulff 构造法预测平衡晶形。本案例 (111) 最低，预测多面体以 (111) 面为优势暴露面，与实验中该晶体形貌吻合。
• 形貌演化判断：若某高 γ 晶面在实际样品中大量出现，说明存在动力学控制或吸附稳定，反过来提示你该调工艺而非质疑计算。

局限

表面能是 0 K 热力学量，未考虑温度、气氛、吸附物种与动力学因素。实际催化或生长环境中，高表面能晶面可能因动力学稳定或吸附稳定而存在。更多 DFT 表面与界面计算可参考本站的第一性原理计算栏目。整体服务见科研学术网首页。