24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5099  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gswylq

木虫 (正式写手)

[求助] VASP Berry phase计算极化强度

我计算BiFeO3的极化强度,选择R-3c为参考结构,但是我无论怎么加U,都是金属性质,因此不能计算Berry phase 极化,请高手指点一下,而相关的文献PHYSICAL REVIEW B 71, 014113却能够计算成功,我不知道在哪儿有问题?
INCAR如下
System = bifeo3
ISTART = 0
ICHARG = 2
ISMEAR = -5
SIGMA = 0.2
#NSW = 100
#IBRION = 2
#IALGO=38
#ISIF = 3
ALGO = Fast
#EDIFF = 1E-5; EDIFFG =-1E-3
ISPIN = 2
MAGMOM = 6*0 1*4 1*4 2*0
#GGA=PE
#VOSKOWN=1
PREC = Accurate
LPLANE = .TRUE.
NPAR = 1
# NSIM = 4
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 2 2 -1
LDAUU=0.0 4.0 4.0 0.0
LDAUJ=0.0 0.0 0.0 0.0
LREAL= Auto
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001
#RWIGS = 0.900 1.323 1.302 2.138
LORBIT = 10
#LELF = .TRUE.
#LREAL = .FALSE.
#LWAVE = .FALSE.
#LCHARG = .FALSE.
#NBANDS=81
#NGX=33
#NGY=33
#NGZ=33
POSCAR如下
feo3                                    
   1.00000000000000     
     4.4362798819116858    1.5524979352336064    2.7205063596444350
    -0.0341933132546832    4.6999640369110791    2.7205059734350789
    -0.0341939503958560   -0.0485663821309455    5.4303274268415738
   O    Fe   Fe   Bi
   6   1   1   2
Direct
  0.4169999994337562   0.0830000005662438   0.2500000000000000
  0.2500000000000000   0.4169999994337562   0.0830000005662438
  0.0830000005662438   0.2500000000000000   0.4169999994337562
-0.4169999994337562  -0.0830000005662438  -0.2500000000000000
-0.2500000000000000  -0.4169999994337562  -0.0830000005662438
-0.0830000005662438  -0.2500000000000000  -0.4169999994337562
  0.2500000000000000   0.2500000000000000   0.2500000000000000
-0.2500000000000000  -0.2500000000000000  -0.2500000000000000
  0.5000000000000000   0.5000000000000000   0.5000000000000000
  0.0000000000000000   0.0000000000000000  -0.0000000000000000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gswylq: 金币+5, 有帮助, 赶紧支持 2012-08-30 18:42:53
liliangfang: 金币+1, 谢谢交流 2012-08-31 07:40:39
金属态ISMEAR不设置为1吗?
未来就是现在
8楼2012-08-30 09:57:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-05-06 22:06:07
楼主,你INCAR 没有设置反铁磁啊,需要将两个Fe离子的磁矩设置成反铁磁的,也就是一正一负。算出来的能带是有带隙的
未来就是现在
16楼2013-03-19 14:06:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)

请问楼主,R-3c算极化的话也要设置反铁磁吗?多谢
未来就是现在
17楼2013-05-06 16:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gswylq 的主题更新
信息提示
请填处理意见