24小时热门版块排行榜    

查看: 1924  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

豆豆王王

铜虫 (小有名气)

[求助] P4mm相BiFeO3态密度

小女纸最近在计算P4mm相BiFeO3态密度,但怎么算都是金属性质,没有带隙,而
Epitaxial BiFeO3 Multiferroic Thin Film Heterostructures
DOI: 10.1126/science.1080615;   Science 299, 1719 (2003)却能算出来,我把我的输入文件贴出来还望大家指正!!!!!
INCAR
SYSTEM = BFO
ENCUT=500
PREC=A
ISMEAR =0
LORBIT=11
ISTART=1

ICHARG=11

ISPIN=2
MAGMOM=2*0 4 -4 6*0

LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 2 2 -1
LDAUU= 0 2 2 0
LDAUJ=0 0 0 0
LDAUPRINT=2


POSCAR
BFO                                    
   3.93500000000000     
     0.9557393861722605    0.0000000000000000    0.0000000000000000
     0.0000000000000000    0.9557393861722605    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.9661910873973605
   Bi   Fe   Fe   O
     2     1     1     6
Direct
  0.0000000000000000  0.0000000000000000 -0.0000000096756567
  0.0000000000000000  0.0000000000000000  0.5000000096756592
  0.5000000000000000  0.5000000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000  0.7500000000000000
  0.5000000000000000  0.0000000000000000  0.2500000000000000
  0.0000000000000000  0.5000000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000 -0.0000000227334085
  0.5000000000000000  0.0000000000000000  0.7500000000000000
  0.0000000000000000  0.5000000000000000  0.7500000000000000
  0.5000000000000000  0.5000000000000000  0.5000000227334159

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00


KPOINTS
Monkhorst Pack
0
Monkhorst Pack
12 12 6
0 0 0
我自己估计可能是POSCAR结构有问题吧,但又不知道该怎么改,大家帮帮忙呀,很急!!!!!!!!
回复此楼

» 猜你喜欢

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

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

dhc198611

木虫 (小有名气)

【答案】应助回帖

★
fzx2008: 金币+1, 谢谢回帖! 2013-04-02 22:17:09
引用回帖:
8楼: Originally posted by 豆豆王王 at 2013-04-02 10:23:49
那20个原子的原子位置是怎样的呢?有两个Fe原子沿【111】方向对称么?

对于四方结构,G型反铁磁在111方向上放铁磁,在111面内是铁磁。40个原子肯定能得到的,通过对角线的基矢变化,就能得到20个原子的模型。
科研无止境
9楼2013-04-02 18:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

identation

金虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-04-01 22:05:03
Fe的U可以大一点,一般可以用4eV,或者3.8eV
``Revisiting properties of ferroelectric and multiferroic thin  films under tensile strain from  first principles,''  Yurong Yang, Wei Ren, Massimiliano Stengel, X. H. Yan and L. Bellaiche, Physical Review Letters 109, 057602 (2012).
``Epitaxial short-period PbTiO3/BiFeO3 superlattices studied by first-principles calculations,''  Yurong Yang, Massimiliano Stengel, Wei Ren, X. H. Yan and L. Bellaiche, Physical Review B 86, 144114 (2012).
l另外,P4mm的BFO应该用20个原子的晶胞,磁性按照G型反铁磁给,也就是ab平面和c方向上Fe的磁矩都应该是反平行,而你在INCAR中似乎没有指定磁矩
ISPIN=2
MAGMOM=。。。。
2楼2013-04-01 02:38:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆王王

铜虫 (小有名气)

其实U=4我也算过,也是没有带隙。
P4mm为什么用20个原子呢?10个就够了吧,Epitaxial BiFeO3 Multiferroic Thin Film Heterostructures
DOI: 10.1126/science.1080615;   Science 299, 1719 (2003)中用的就是10个。而且磁矩我设置了呀
ISPIN=2
MAGMOM=2*0 4 -4 6*0
3楼2013-04-01 08:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-04-01 22:05:17
引用回帖:
3楼: Originally posted by 豆豆王王 at 2013-04-01 08:45:10
其实U=4我也算过,也是没有带隙。
P4mm为什么用20个原子呢?10个就够了吧,Epitaxial BiFeO3 Multiferroic Thin Film Heterostructures
DOI: 10.1126/science.1080615;   Science 299, 1719 (2003)中用的就是10个 ...

10个原子无法造出G型反铁磁,你给的磁性应该是C型的 你要保证在ab平面内体系也是磁矩反平行的
4楼2013-04-01 09:12:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:46 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:38 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:27 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 FChKX8meu1ds 2026-11-09 8/400 2026-11-14 22:04 by Equinoxhua
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-14 21:58 by Equinoxhua
[考博] 27博士 +6 锦鲤送福至 2026-11-10 6/300 2026-11-14 21:23 by xupeii
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:00 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 FChKX8meu1ds 2026-11-09 12/600 2026-11-12 09:43 by HajghKOMynmz
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 06:03 by 9fOTqj074tIg
信息提示
请填处理意见