24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1106  |  回复: 8

calos818

木虫 (著名写手)

新人菜鸟

[交流] 【求助】dos计算为导体,和实验不符合,问题出在哪

磁性材料,U值取自别人的计算文献,计算和实验都为绝缘体,采用和文献同样的设置,我的dos算下来是导体,且改变U值还是导体,似乎U值不影响。大致设置如下,麻烦大家帮我分析,可能的出错原因,非常感谢!
ISIF=2
NSW=0
NELM=200
IBRION=-1
LDAU      = .TRUE.
LDAUTYPE  = 2
LDAUL     = -1 2 -1
LDAUU     = 0 5.0 0
LDAUJ     = 0 1.0 0
LDAUPRINT = 2
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

★ ★
zzy870720z(金币+2):谢谢提示 2010-08-12 18:46:22
引用回帖:
Originally posted by calos818 at 2010-08-12 14:53:15:
磁性材料,U值取自别人的计算文献,计算和实验都为绝缘体,采用和文献同样的设置,我的dos算下来是导体,且改变U值还是导体,似乎U值不影响。大致设置如下,麻烦大家帮我分析,可能的出错原因,非常感谢!
ISIF= ...

ISMEAR=-5      ?
好像对于绝缘体,vasp说明上提到用ISMEAR>1可能产生错误,而它的默认值就是1,因此不设置的话也许问题就出在这。
2楼2010-08-12 17:55:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

引用回帖:
Originally posted by nkleof at 2010-08-12 17:55:40:

ISMEAR=-5      ?
好像对于绝缘体,vasp说明上提到用ISMEAR>1可能产生错误,而它的默认值就是1,因此不设置的话也许问题就出在这。

您好,我确实用的-5的ISMEAR,采用这种方法 可能的出错原因是什么,麻烦您指点下
3楼2010-08-12 19:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)


calos818(金币+1): 2010-08-12 19:41:20
zxzj05(金币+1):奖励一下 2010-08-12 20:19:58
引用回帖:
Originally posted by calos818 at 2010-08-12 19:12:38:

您好,我确实用的-5的ISMEAR,采用这种方法 可能的出错原因是什么,麻烦您指点下

那就不知道了,我看你的帖子里没有写出ISMEAR的设置,以为你用的是默认值,才回了那个贴。既然你设置了ISMEAR,那看来就不是了。
不过我觉得你应该把计算设置列得更详细一些,方便大家找问题。
4楼2010-08-12 19:26:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

5楼2010-08-13 18:03:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)


calos818(金币+1): 2010-08-14 07:49:25
zzy870720z(金币+1):谢谢提示 2010-08-14 08:36:54
DFT的错误,低估带隙...著名的例子参见Fe.
6楼2010-08-14 07:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

引用回帖:
Originally posted by condensed at 2010-08-14 07:11:14:
DFT的错误,低估带隙...著名的例子参见Fe.

您好,别的计算和实验文献坐下来确实是绝缘体,这点是肯定的,问题是我的j计算不对l可能原因是什么,请教下您
7楼2010-08-14 07:49:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
zzy870720z(金币+2):谢谢提供思路 2010-08-14 08:37:12
这样多半是你的设置有误,从晶体结构开始到INCAR逐个检查。或者画出被加U的那个元素在不同U下的分波态密度图,看它随U的变化。
8楼2010-08-14 08:04:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


zzy870720z(金币+1):谢谢侯老师补充 2010-08-14 08:37:28
Fe单质没有能隙。Fe好像是LDA不能正确算出Fe的基态,需GGA。

[ Last edited by valenhou001 on 2010-8-14 at 08:09 ]
9楼2010-08-14 08:07:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 calos818 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见