计算化学公社

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

作者
Author: DHT 时间: 2025-6-4 23:20
标题: 求助掺杂稀土离子体系结构优化，c/a跟文献不一样是什么问题导致的
刚开始接触掺杂相关的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（跟其他文章接近）】

欢迎光临计算化学公社 (http://bbs.keinsci.com/)