24小时热门版块排行榜

返回列表

gswylq

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2546.6
散金: 120
红花: 1
帖子: 782
在线: 730.1小时
虫号: 659377
注册: 2008-11-21
专业: 金属材料的凝固与结晶学

[求助] 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

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

请问LCALCPOL和LBERRY计算的极化有什么区别已经有55人回复
有关铁电材料极化的问题已经有2人回复
用vasp脚本算极化强度，脚本出错，急死了已经有6人回复
一个关于铁电薄膜极化翻转的问题已经有7人回复
vasp接口wannier90计算问题已经有24人回复
vasp自旋极化计算时，初始设置自旋能否不被优化已经有6人回复
请问vasp里面有Kubo-Greenwood公式的计算么已经有7人回复
利用VASP算压电极化率已经有10人回复
Vasp能计算电荷转移的具体数值么？已经有6人回复
求问，VASP计算born有效电荷后估算极化强度的方法已经有19人回复
压电陶瓷极化前后，剩余极化强度为什么差别这么大？已经有7人回复
berry phase 参考态已经有35人回复
berry phase！DIPOL? 已经有10人回复
多铁材料铁电极化（Berryphase方法）计算问题已经有6人回复
用berry phase计算Born effective charge时，最后一步的正负号问题已经有8人回复
【求助】vasp Berry phase计算自发极化已经有4人回复
【求助】VASP怎么考虑自旋极化？已经有14人回复
【求助】VASP计算出的电极化强度数据，怎么做后处理已经有15人回复
【求助成功】如何用VASP计算反铁磁已经有14人回复

1楼 2012-08-29 16:29:03

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

回帖支持 ( 显示支持度最高的前 50 名 )

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
散金: 118
红花: 76
帖子: 1362
在线: 543.8小时
虫号: 1315985
注册: 2011-06-05
专业: 无机非金属类电介质与电解

★
xueht987: 金币+1, 鼓励交流！ 2012-08-31 09:32:41

引用回帖:

11楼: Originally posted by gswylq at 2012-08-31 07:49:12
POTCAR采用PBE的，其中Fe的采用14个价电子的。具体如下。也不知道这是什么原因导致的？
O的如下：
PAW_PBE O_s 07Sep2000
6.00000000000000000
parameters from PSCTR are:
VRHFIN =O: s2p4
LEXCH ...

那你应该考虑的是GGA+U了，你的INCAR好像没有反映这一点
Bi PAW应该是用10个价电子那个。我计算中采用的是LDA+U，有比较小的bandgap。