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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » P4mm相BiFeO3态密度
101/1返回列表
查看: 1459  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

豆豆王王

铜虫 (小有名气)

[求助] 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结构有问题吧，但又不知道该怎么改，大家帮帮忙呀，很急！！！！！！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-31 21:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

豆豆王王

铜虫 (小有名气)
Note that computations were performed with a 10-atom body-centered tetragonal unit cell, which was
doubled to accommodate the G-type antiferromagnetic ordering。 文献中确实有这样的话，难道是我理解错误？？
赞一下 回复此楼
5楼2013-04-01 10:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kollar
6楼2013-04-01 11:34:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhc198611

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流，， 2013-04-01 22:05:26
引用回帖:
5楼: Originally posted by 豆豆王王 at 2013-04-01 10:30:42
Note that computations were performed with a 10-atom body-centered tetragonal unit cell, which was
doubled to accommodate the G-type antiferromagnetic ordering。 文献中确实有这样的话，难道是我理解错 ...

三方结构R3c是可以用10个原子构建出G型反铁磁，四方结构P4mm是要用20个原子的。
赞一下 回复此楼
科研无止境
7楼2013-04-01 19:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆王王

铜虫 (小有名气)
那20个原子的原子位置是怎样的呢？有两个Fe原子沿【111】方向对称么？
赞一下 回复此楼
8楼2013-04-02 10:23:49
已阅   回复此楼   关注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的回帖

豆豆王王

铜虫 (小有名气)
引用回帖:
9楼: Originally posted by dhc198611 at 2013-04-02 18:36:02
对于四方结构，G型反铁磁在111方向上放铁磁，在111面内是铁磁。40个原子肯定能得到的，通过对角线的基矢变化，就能得到20个原子的模型。...

什么叫对角线的基矢变化？基础差的人表示一头雾水。
赞一下 回复此楼
10楼2013-04-03 11:03:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 豆豆王王 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭