CsPbBr3上下表面不对称如何计算表面能

zyd

各位老师好我想咨询的问题是关于上下表面不对称的表面能的计算（体块以及表面的截断能等参数都一致，k点相差不大）
（1）对体块结构CsPbBr3体块进行了优化，体块的原子数为5。得到能量Ebluk

（2）我切体块（001）方向并且厚度为3，原子数为15，构建真空层（20埃米）后，固定最下一层进行优化没得到能量Eslab，但是上下两层终止层不对称（一层是Pb-Br,一层是Cs-Br）,
所以公式Y=(Eslab-nEbluk)/2A,计算出来相当离谱

（3）若是保持上下表面对称的话。只能截取半数个厚度，并且切出来的表面原子数不等于体块的整数倍，如图所示（厚度2.5.原子数12，不是整数倍）所以,这个公式也不能用

所以想请问老师们，对于这样的体系怎么算能呢，
还有层数收敛性测试还能用什么方法呢

zyd

有相关的文章的话，也请各位老师介绍一下