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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP Berry phase计算极化强度
191/2返回列表12下一页
查看: 5297  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 329求调剂 +8 想上学吖吖 2026-03-19 8/400 2026-03-20 21:07 by 无际的草原
[考研] 求调剂 +5 Mqqqqqq 2026-03-19 5/250 2026-03-20 20:46 by zhukairuo
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 299求调剂 +3 shxchem 2026-03-20 4/200 2026-03-20 19:23 by zhukairuo
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 304求调剂 +5 曼殊2266 2026-03-18 5/250 2026-03-20 09:00 by ZHANG0tao
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
信息提示
关闭
请填处理意见
关闭