24小时热门版块排行榜    

查看: 1833  |  回复: 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的回帖

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

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

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

04nylxb

木虫 (正式写手)

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

用ismear=0后的DOS还是在远高于EF的地方不回零。
想请教下,Ni表面吸附乙醇这样的体系,ISMEAR方法不同是否会导致错误的结果?
集中精力发文章
6楼2012-06-05 13:47:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)


xueht987: 金币+1, 谢谢回帖,鼓励交流! 2012-06-06 09:03:25
引用回帖:
6楼: Originally posted by 04nylxb at 2012-06-05 13:47:11
用ismear=0后的DOS还是在远高于EF的地方不回零。
想请教下,Ni表面吸附乙醇这样的体系,ISMEAR方法不同是否会导致错误的结果?...

道理上讲只要收敛了,SIGMA也足够小,结果应该都是一样的。但我之前有过“换个ISMEAR就没有空带dos“的情况,所以推荐你试一试抱歉呢,没解决问题
楼主是关心体系的什么性质呢?一般那些dos都没什么影响的:根本没电子往那去
7楼2012-06-05 21:58:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04nylxb

木虫 (正式写手)

引用回帖:
7楼: Originally posted by mousekingadv at 2012-06-05 21:58:52
道理上讲只要收敛了,SIGMA也足够小,结果应该都是一样的。但我之前有过“换个ISMEAR就没有空带dos“的情况,所以推荐你试一试抱歉呢,没解决问题
楼主是关心体系的什么性质呢?一般那些dos都没什么影响的: ...

收到,非常感谢,哈哈。
我是想作下表面吸附的成键情况
集中精力发文章
8楼2012-06-06 09:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 04nylxb 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见