求助掺杂稀土离子体系结构优化，c/a跟文献不一样是什么问题导致的

DHT

刚开始接触掺杂相关的vasp计算，有很多基础的问题没有明白。
文献做的是BTO,PTO,KNO在A位/B位掺杂稀土离子相关内容，超胞是2*2*2，稀土离子替换一个BTO的Ba，文中有关赝式的描述如图1，所以我选了Ce_3进行后边的复现，就是一个Ce替换BTO中的一个Ba，但是结构优化过后的结果跟文章中差别很大，文章结果2图，可以看到c/a接近1，但是我的还是1.02左右，而且从2图也能知道他开始用的也是四方相，但我的没有优化到立方附近。
INCAR 如下
SYSTEM=BTO         ##正常结构优化               
##        start parameter for this run in/out
ISTART=0         
ICHARG=2                     
LREAL=Auto
ADDGRID= .TRUE.
LCHARG=.F                  
LWAVE=.F                  
NPAR=4                           
ISYM=0                  
                        
##        ionic relaxation               
IBRION =2                    
ISIF=3                           
NSW =60                          
EDIFFG=-0.01                           
##POTIM=0.001
##        electronic relaxtion         电子弛豫
ENCUT=500                          
NELM=100                    
EDIFF=1E-6                  
PREC=Accurate                          
ISMEAR=0
SIGMA=0.05         
##GGA=PS        

##        加U
LDAU = .TRUE.          #（打开+U设置，控制计算中是否考虑在位库伦较正相）
LDAUTYPE = 2          #（+U的类型，1/2/4；2为默认值，U值由LDAUU减去LDAUJ确定）
LDAUL = -1 -1 -1 3          # (控制具体轨道上加U；-1：不加U；1-p轨道；2-d轨道；3-f轨道)
LDAUU = 0 0 0 5.5         #（O不加U,Ce的f轨道上+U，U值取5）
LDAUJ = 0 0 0 0.5           
LMAXMIX = 6          #（对+U日系，设置线性混合参数，f轨道设为6，d轨道设为4）         

为什么结果会差这么多呢，希望有明白的人能指点一下，谢谢     

另外：关于铈赝式的选择，如果选择了Ce_3还需要加U计算吗 【我前边不加U算过DOS，只有0.4，加上了有1.7（跟其他文章接近） 】