求助铁电性、极化强度问题

Ava

我是研究材料的铁电性的，想搞出双势阱图，就先拿简单的KNbO3试，方法是这样的，先找一个中心对称的结构，把它当成顺电相，poscar文件如下：K1 Nb1 O3
1.0
        4.0574688911         0.0000000000         0.0000000000
        0.0000000000         4.0574688911         0.0000000000
        0.0000000000         0.0000000000         4.0574688911
    K   Nb    O
    1    1    3
Direct
     0.000000000         0.000000000         0.000000000
     0.500000000         0.500000000         0.500000000
     0.000000000         0.500000000         0.500000000
     0.500000000         0.000000000         0.500000000
     0.500000000         0.500000000         0.000000000

然后，我把Nb原子在Z方向的0.5改成从 0.40 0.42  ..... 0.58 0.60 其他的坐标不动，得到这么多poscar文件，把这些当成铁电相结构，然后再优化，和一开始的顺电相一起算出极化强度，然后把铁电相能量减去顺电相能量，得到的值和极化强度作图，没有得到文献中的形状，不知是哪的问题，请各位老师，师兄解答一下，特别感谢