CeO2(100)晶面建模与氧空位形成能计算

Jihang

我用VASP计算CeO2不同晶面的氧空位形成能，用MS切的晶面，文献上(100)面的氧空位形成能是高于(110)的，可是我怎么计算都是(100)晶面的氧空位形成能更低，不知道是不是建模方面的问题，请各位老师指点一下。

U加了5，开了偶极校正，开了自旋，K点网络用的的1 1 1，我试了多切几层，也试了扩大超胞，怎么算都不行，不知道是哪里的问题了。

我切出来的(100)面的Ce-O是晶胞的对角方向，文献上看起来是我这个结构旋转45°的情况，请问这两种建模对结果会不会有影响呀？还有我应该怎么调整我切的晶面的方向呢？还是说有其他需要注意的因素？