24小时热门版块排行榜    

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

04nylxb

木虫 (正式写手)

[求助] 请教表面吸附DOS计算中Ef很大于0后仍然不回0

请教大家,我用vasp5.2.12计算DOS,计算普通的金属表面(如Ni111)吸附乙醇这样的体系,出来的DOS(请见图),在能量远大于费米能的地方,如E-Ef=5的地方,自旋向上和向下的DOS都没有回到0,请问是体系的原因还是我INCAR设置的原因?在能量高于费米能5eV或更高是否还会有电子存在?非常感谢
我的INCAR如下:
PREC=Accurate
ISTART=0 ; ICHARG=11
ENCUT=400
ISMEAR=-5
SIGMA=0.05
ALGO=F
LORBIT=11
EMIN=-15
EMAX=15
NEDOS=4500
ISPIN=2
ISYM=0
EDIFF=1E-4;
LPLANE=.TRUE.
NPAR=4
NSIM=4
LREAL=Auto

DOS
回复此楼

» 猜你喜欢

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

集中精力发文章
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xueht987: 金币+2, 鼓励交流,谢谢 2012-05-31 10:03:41
04nylxb: 金币+5, ★★★很有帮助, 非常感谢,就是能量远高于费米能级的地方也有这么多的电子,当时觉得是不是设置错了 2012-05-31 13:38:23
Ni 是有磁性的,你设置 ISPIN=2 是正确的。从你的DOS上可以看出体系有磁性的,因为自旋向上、向下的电子态密度是不对称分布的。
还有,DOS 计算的是电子在不同能级上可能的占据状态。具体的占据与否要看体系的费米能级,因为费米能级表示体系电子的最高占据能级。也就是说,基态情况下费米能级以下的态是被占据的,以上的态是没有电子占据的。
3楼2012-05-31 08:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

likai0106

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
xueht987: 金币+1, 谢谢回帖,鼓励交流! 2012-05-31 10:03:23
04nylxb: 金币+3, ★★★很有帮助, 非常感谢 2012-05-31 13:36:24
这个正常吧,本身Ni属于磁性原子,如果都对称了也没有磁性了。态穿过费米面说明Ni表面有弟子丢失,这个正常,催化么表面肯定有电子去辅助反应的分解。
2楼2012-05-31 07:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
franch: 金币+2, 谢谢回帖交流 2012-06-04 14:07:06
04nylxb: 金币+2, ★★★很有帮助, 收到,非常感谢,我试下,感觉很有可能 2012-06-04 22:33:11
04nylxb: 金币+3, 有帮助, 试了下,用ismear=0还是不归零……谢谢啊 2012-06-05 13:43:37
我觉得这个和ISMEAR有关系(确实很奇怪),我以前用ISMEAR=-5就有,用0就没有了。楼主可以试一下
4楼2012-06-03 17:22:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04nylxb

木虫 (正式写手)

引用回帖:
4楼: Originally posted by mousekingadv at 2012-06-03 17:22:03
我觉得这个和ISMEAR有关系(确实很奇怪),我以前用ISMEAR=-5就有,用0就没有了。楼主可以试一下

非常感谢啊,我试试,到时候再补加金币
集中精力发文章
5楼2012-06-04 22:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见