| 查看: 3359 | 回复: 17 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
NEB的计算出现负值，测试的其他能量都比初始的00和终止的10低
|
||
|
问题如下： 1.计算NEB的方法：首先利用nebmake.pl产生初始终态结构，用xcrysden查看了路径没问题，然后对初始终态结构进行几何结构优化，接着采用优化后的能量构型进行NEB的计算，按照VASP manual上说明，在INCAR中设置如下： ISTART = 0; ICHARG = 2 ENCUT = 400 ISMEAR = 2 LREAL=Auto NELM=200 NSW=1000 EDIFFG=-0.02 LWAVE = F LCHARG = F ICHAIN=0 IMAGES=9 SPRING=-5 LCLIMB=.TRUE. IBRION=3 POTIM=0 IOPT=7 算完之后用nebbarrier.pl产生的neb.dat如下： 0 0.000000 0.000000 0.000000 0 1 0.151183 -2.145471 -0.242258 1 2 0.302705 -2.106921 -0.362263 2 3 0.452469 -2.051050 -0.403422 3 4 0.603020 -2.012050 -0.282203 4 5 0.753502 -1.992645 0.000084 5 6 0.904120 -2.012146 0.253874 6 7 1.054939 -2.051164 0.407455 7 8 1.204729 -2.107090 0.352265 8 9 1.356227 -2.145477 0.236853 9 10 1.507559 0.000000 0.000000 10 初末位置的能量我多次验算都是最低的，不知道是什么原因，求大神们不吝指教我这个新手！谢谢啦 |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有12人回复
-
2026国自然函评费到账
已经有3人回复
-
科研孤儿太难了
已经有17人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有19人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 进行过渡态计算时link died
已经有10人回复
4楼2013-11-04 13:10:11
liu-biao
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 4598.4
- 散金: 84
- 帖子: 464
- 在线: 247.2小时
- 虫号: 2102358
- 注册: 2012-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-11-04 12:16:40
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
3楼2013-11-04 12:43:21
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
5楼2013-11-04 16:07:14
40