应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP Berry phase计算极化强度
191/2返回列表12下一页
查看: 5295  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-29 16:29:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

iamikaruk

木虫 (著名写手)
★
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。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

liutaifeng
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
12楼2012-08-31 08:18:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

iamikaruk

木虫 (著名写手)
★
xueht987: 金币+1, 鼓励交流！ 2012-08-30 09:27:39
你这个算出来是反铁磁么？铁磁似乎是金属
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
2楼2012-08-29 17:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gswylq

木虫 (正式写手)
我直接静态计算的，用反铁磁态，铁磁态也算了，还是金属
赞一下 回复此楼
3楼2012-08-29 17:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gswylq

木虫 (正式写手)
引用回帖:
2楼: Originally posted by iamikaruk at 2012-08-29 17:03:28
你这个算出来是反铁磁么？铁磁似乎是金属

我直接静态计算的，用反铁磁态，铁磁态也算了，还是金属，
原子内坐标根据的是文献的值（PHYSICAL REVIEW B 71, 014113）
赞一下 回复此楼
4楼2012-08-29 17:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

umklapp

银虫 (小有名气)
★
xueht987: 金币+1, 鼓励交流！ 2012-08-30 09:28:09
看不出 那个地方不对，要是完全按照文献的参数， 应该是能重复出来的， 况且又是牛人的paper;
你试试ISMEAR=0?(和文献不同) 我记得如果fermi level 上又多条能带的crossing， 可能会导致tetrahedron method的问题；
赞一下 回复此楼
5楼2012-08-30 00:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
xueht987: 金币+1, 鼓励交流！ 2012-08-30 09:28:17
弱弱地提醒楼主POTCAR要有两个Fe，即四个POTCAR依次合并在一起。。。。。。
赞一下 回复此楼
6楼2012-08-30 00:07:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★
xueht987: 金币+1, 鼓励交流！ 2012-08-31 09:32:13
引用回帖:
3楼: Originally posted by gswylq at 2012-08-29 17:07:10
我直接静态计算的，用反铁磁态，铁磁态也算了，还是金属

你给的spin是铁磁态啊。能够自洽收敛到反铁磁？
如果收敛到反铁磁态的话，你的磁矩是多少？我以前算的结果是3.8~4.0左右。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
7楼2012-08-30 08:28:16
已阅   回复此楼   关注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的回帖

gswylq

木虫 (正式写手)
引用回帖:
7楼: Originally posted by iamikaruk at 2012-08-30 08:28:16
你给的spin是铁磁态啊。能够自洽收敛到反铁磁？
如果收敛到反铁磁态的话，你的磁矩是多少？我以前算的结果是3.8~4.0左右。...

反铁态的R3C结构算了，Fe的磁矩是3.69， R-3C反铁磁态Fe的磁矩是3.499.可能是我加的U还不够吧。
赞一下 回复此楼
9楼2012-08-30 18:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
9楼: Originally posted by gswylq at 2012-08-30 18:41:59
反铁态的R3C结构算了，Fe的磁矩是3.69， R-3C反铁磁态Fe的磁矩是3.499.可能是我加的U还不够吧。...

U=4.0 eV足够了。
不应该有问题的啊，奇怪。你用的Fe PAW是哪一个？14个价电子还是16个价电子的？
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
10楼2012-08-30 21:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gswylq 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +5 怀瑾握瑜l 2026-03-20 5/250 2026-03-20 20:24 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 085700资源与环境308求调剂 +11 墨墨漠 2026-03-18 12/600 2026-03-20 19:43 by 丁丁*
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 材料学硕297已过四六级求调剂推荐 +6 adaie 2026-03-19 6/300 2026-03-20 14:12 by 星空星月
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂 +6 吃吃吃才有意义 2026-03-19 6/300 2026-03-20 10:47 by 尽舜尧1
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
关闭
请填处理意见
关闭