24小时热门版块排行榜    

查看: 1876  |  回复: 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的回帖
信息提示
请填处理意见