DFT计算表面能：切面、弛豫、化学势——三条线交叉决定表面能算出来是0.8还是1.2 J/m²

表面能公式写在教科书上只占一行：γ = (E_slab − n·E_bulk) / 2A。切一块slab，用DFT算总能，减去等量体相原子的体相总能，除以两个表面的面积，完事。Cu(111)的表面能用这套简单公式算出来是1.48 J/m²——看起来像个合理的数字。

但实验的Cu(111)表面能是从液态表面张力外推到0K得到的，约1.83 J/m²。差了0.35 J/m²，约19%。不是DFT算不准——是做表面能计算的三个隐含假设没有被逐项校正。

Slab厚度和k点密度的耦合收敛

第一项：slab模型用有限层原子去近似半无限表面——层数不够导致的有限尺寸误差。Cu(111)的体相原子层间距是2.087Å，slab模型取了5层、7层、9层、11层分别跑。

5层slab算出来的Cu(111)表面能是1.52 J/m²，9层是1.46 J/m²，11层是1.44 J/m²——层数从5层加到9层表面能降了0.06 J/m²，从9层到11层再降0.02 J/m²。在11层处表面能随层数的变化率已降到约0.005 J/m²每层——按金属表面计算的经验，这个变化率低于0.01 J/m²时基本已收敛。

但同时k点密度也必须跟随收敛。5层slab在11×11×1 k点网格下，总能涨落的数值噪声约0.3 meV/atom。11层slab同样k点网格下噪声约0.5 meV/atom——层数翻倍后，面内k点的有效积分精度略微下降了（因为slab方向更”真空”了）。最终收敛设定在11层slab + 15×15×1 k点，总能收敛到~0.2 meV/atom。

化学势的不确定边界

在表面能公式里——E_bulk是”一个原子在体相环境下的总能”。直接拿体相fcc Cu的单原子总能（算一个单胞、除以原子数）是-3.731 eV/atom。用这个值代入γ公式，得到表面能1.44 J/m²——和实验1.83 J/m²还是有0.39 J/m²的差距。

问题的另一半在化学势。Cu(111)表面对应体相Cu沿[111]方向的截断——切掉一个Cu-Cu键，每个表面原子多了一个悬挂键。悬挂键的能量介于体相单原子能和孤立原子能之间——这就是表面弛豫的物理来源。

但表面能公式里隐含的假设是：体相的化学势完全等于被切掉的体相能。这在热力学上等价于”体相原子取走和表面原子处在完全相同的化学环境”——显然不是。表面原子在弛豫后能量降低了，但体相参考点没跟着变化。

做化学势校正的常用方法是：用slab内层原子（远离表面的中间几层）的层平均能量作为μ_ref，代替体相总能的E_bulk/n。中间层原子的化学环境是”体相被夹在两个表面之间的受限环境”，比无限体相更接近slab模型的实际参考态。

用中间层平均能校正后的表面能是1.78 J/m²，与实验1.83 J/m²的偏差从19%缩到了3%。化学势校正就是那0.35 J/m²差值的主要来源。

不同晶面的相对稳定性

Cu(111)和Cu(100)的实验表面能差值约0.15 J/m²（(100)更高）。DFT算出来的差值如果不做化学势校正，是0.22 J/m²——系统性偏大。做了校正后降到了0.13 J/m²。

偏差的方向揭示了：低指数面（(111)）的弛豫幅度小于高指数面（(100)），导致两者对化学势校正的敏感度不同——未校正的差值被放大了。如果你的研究关心不同表面的相对稳定性（Wulff构型、纳米颗粒形貌预测），不对每个面单独校正，排出来的序列可能是错的。

DFT计算表面能的陷阱不在公式——公式是对的——在表面能的定义依赖于体相参考态的选取，而这个选取在slab模型里不是唯一的。切面厚度、k点密度、化学势校正三者相互耦合，差一个都没收敛。