24小时热门版块排行榜

返回列表

gooselin

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1900.1
帖子: 79
在线: 87.2小时
虫号: 3809417
注册: 2015-04-14
专业: 凝聚态物性 II ：电子结构

[求助] 急！VASP计算铁电极化时提示体系非绝缘体已有1人参与

最近在用VASP计算铁电极化，用了Berry Phase 的方法计算的，我分别计算了加入自旋轨道耦合（SOC）和没有加SOC的极化，但是加入SOC计算时总是在最后一步提示体系不是绝缘体，我的体系能带结构计算过的，绝对是绝缘的。所以不知道哪里出了问题。我的计算步骤是
1：加上磁性的结构优化；2：加上磁性的静态计算得到CHGCAR；3：读取上一步的CHGCAR分别计算加SOC与不加SOC的极化

错误提示：
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CALC_POLARIZATION: ERROR, system is not insulating, skipping...
CALC_POLARIZATION: ERROR, system is not insulating, skipping...
CHECK_OCCUPATIONS: ERROR, system is not insulating.

加入SOC的计算铁电极化的INCAR：
ISTART=0
ICHARG=11
ENCUT=600
ISMEAR=-5
#SIGMA=0.05
#NSW=100
NELM=200
#NELMIN=2
#IBRION=2
#POTIM=0.1
#ISIF=3
ISPIN=2
#MAGMOM=0 0 4 4 4 -4 -4 -4 5 -5 5 -5 5 -5 5 -5 24*0
MAGMOM=6*0 0 0 4 0 0 4 0 0 4 0 0 -4 0 0 -4 0 0 -4 0 0 5 0 0 -5 0 0 5 0 0 -5 0 0 5 0 0 -5 0 0 5 0 0 -5 72*0
LSORBIT=.TRUE.
LMAXMIX=6
SAXIS=1 0 0
NBANDS=432
ISYM=0
GGA_COMPAT=.FALSE.
LORBMOM=.TRUE.
NGX=62
NGY=62
NGZ=62
NCORE=1
EDIFF=1E-6
#EDIFFG=-0.03
PREC=Accurate
LREAL=Auto
LASPH=.TRUE.
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 2 2 -1
LDAUU=0 5.5 4 0
LDAUJ=0 1.5 1 0
LCALCPOL=.TRUE.
DIPOL=0.5 0.5 0.5
LWAVE=.FALSE.
LCHARG=.FALSE.

回复此楼