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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP Berry phase计算极化强度
191/2返回列表12下一页
查看: 5273  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭