24小时热门版块排行榜    

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

w474964260

新虫 (初入文坛)

[求助] NEB计算扩散能垒出现鞍点过高,求大神帮忙分析!已有1人参与

用NEB方法计算H扩散的能垒,计算出的结果与文献中的曲线相似,但是鞍点(峰值为0.21ev)与文献差距很大(文献的峰值为0.074ev),不知道是哪里出问题,困惑很久。计算过程大概,第一:优化初始态和终态,第二:用nebmake.pl分配若干image点,第三:VASP跑起。INCAR的参数截图已经贴出来,另外贴出nebef文件mep.eps图,初学VASP,希望不吝赐教。也为后面学NEB的同学,立基础贴,相互讨论
另外还有2个问题:
1.如何验证频率,确定过渡态是否正确,看过很多贴,只知道需要验证频率,但还没有具体学会
2.H扩散的位置在初始态的POSCAR已经给出,若是几何优化后,不就会改变初始的位置吗,与自己的设想位置不符,还是在优化的时候,要设定好ISIF参数,不让H原子在优化时跑动?NEB计算扩散能垒出现鞍点过高,求大神帮忙分析!
INCAR.png


NEB计算扩散能垒出现鞍点过高,求大神帮忙分析!-1
文献上的峰值.png
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : mep.eps
  • 2015-06-25 16:52:45, 17.78 K

» 猜你喜欢

永远不要埋没自己
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑然49

金虫 (小有名气)

引用回帖:
13楼: Originally posted by only_大师 at 2016-09-09 13:50:16
请问一下楼主,mep.eps是怎么生成的啊,我的为什么没有啊

计算结束后,使用命令nebbarrier.pl和nebspline.pl对计算结果进行处理,会出现neb.dat、spline.dat和mep.eps等文件
smile
14楼2016-09-09 16:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

xxxxxxxfc

金虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流。 2015-06-29 20:04:29
不清楚,你这个体系是什么,至少对我来说0.074ev是不是太低啦,完全可以忽略不计,相反0.21ev看着还正常些。
2楼2015-06-26 09:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w474964260

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by xxxxxxxfc at 2015-06-26 09:00:50
不清楚,你这个体系是什么,至少对我来说0.074ev是不是太低啦,完全可以忽略不计,相反0.21ev看着还正常些。

非常感谢回复,这是V体系的氢渗透,对于能垒的是否为正常值,我还不会确定,经过查文献大都在0.07~0.15ev之间,而这里的0.074ev确实是文献上给出的,但我根据文献计算,得不到这个值。
永远不要埋没自己
3楼2015-06-26 10:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
w474964260: 金币+5, 有帮助 2015-06-29 11:21:17
引用回帖:
3楼: Originally posted by w474964260 at 2015-06-26 10:34:38
非常感谢回复,这是V体系的氢渗透,对于能垒的是否为正常值,我还不会确定,经过查文献大都在0.07~0.15ev之间,而这里的0.074ev确实是文献上给出的,但我根据文献计算,得不到这个值。...

至少我的理解,我这儿一般的体系,势垒一般在0.2到1点几,才正常。
4楼2015-06-26 11:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见