| 查看: 5874 | 回复: 18 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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 |
回复此楼» 猜你喜欢
请问社科规划办官网上的“项目查询模块”与“项目申报模块”的账号密码一样吗？
已经有3人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有9人回复
-
40+中年人失业后，还能找到好工作吗？
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
极限拉扯，如何抉择
已经有15人回复
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
★
xueht987: 金币+1, 鼓励交流！ 2012-08-31 09:32:41
xueht987: 金币+1, 鼓励交流！ 2012-08-31 09:32:41
|
那你应该考虑的是GGA+U了，你的INCAR好像没有反映这一点 Bi PAW应该是用10个价电子那个。我计算中采用的是LDA+U，有比较小的bandgap。 |
liutaifeng
12楼2012-08-31 08:18:53
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
2楼2012-08-29 17:03:28
3楼2012-08-29 17:07:10
4楼2012-08-29 17:09:22