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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP Berry phase计算极化强度
191/2返回列表12下一页
查看: 5778  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +7 wangmingqi 2026-08-28 8/400 2026-08-28 18:42 by zhuzg0628
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭