计算离子型晶体表面能的思路是否正确？

27@790624905

请教各位老师，在用CASTEP计算晶体表面能时，图中所给公式的思路计算是否正确（即切表面后建立6*6*6的体相晶体，结构优化计算表面能Eslab；计算通用胞中的单个晶胞的体相能Ebulk；按照6*6*6的体相晶体给出a*b得到表面积；依据公式计算得出表面能）。目前出现的问题是：单个晶胞结构优化出来的Final energy（Ebulk）和切表面后建立6*6*6的体相晶体结构优化得到的表面能（Eslab）差值为2eV。请问这个问题具体出在哪里？

丁越

E_bulk不用你构建超胞计算，直接计算原胞就行，这里的E_bulk指的是每个原子的能量（即晶胞能量/晶胞原子数）。然后N指的是slab所含的原子数。
另外你的slab构建必须要足够厚，并且需要将slab的中间部分固定住以体现体相环境。

27@790624905

丁越 发表于 2025-12-7 10:11
E_bulk不用你构建超胞计算，直接计算原胞就行，这里的E_bulk指的是每个原子的能量（即晶胞能量/晶胞原子数 ...

感谢了，经过你这个提点确实能够算出来的！

27@790624905

丁越 发表于 2025-12-7 10:11
E_bulk不用你构建超胞计算，直接计算原胞就行，这里的E_bulk指的是每个原子的能量（即晶胞能量/晶胞原子数 ...

还有一个问题想请教，我按这个思路一共计算了6层晶胞，计算出来的表面能比文献中的给的低了一个量级，单位之间的换算是没有问题的，这个可能是什么原因呢？

丁越

27@790624905 发表于 2025-12-7 19:38
还有一个问题想请教，我按这个思路一共计算了6层晶胞，计算出来的表面能比文献中的给的低了一个量级，单 ...

没具体数据没法说

27@790624905

你好，再请教一下！参考的文献计算参数和我实际计算参数对比，计算过程的实际情况如附件，请指导一下是哪里存在问题导致实际计算氟化镁表面能（0.87 J/cm2）低于文献表面能（1.5 J/cm2）？