公式： 参考资料：https://en.wikipedia.org/wiki/Surface_energy

其中，

• $E_{slab}$，DFT计算表面得到的总能量（或总键能）。在BAND模块中，生成Slab模型完成计算后，logfile尾部有总键能
• $N$，表面单胞的原子个数
• $E_{bulk}$，体材料中每个原子的能量（或键能）。在BAND中，Bulk模型（也就是三维周期体系）计算完毕后，logfile尾部有总键能，平均到每个原子即可。
• A，表面积（注意有上下两个表面），因此乘以2。