| 查看: 3446 | 回复: 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 初末位置的能量我多次验算都是最低的，不知道是什么原因，求大神们不吝指教我这个新手！谢谢啦 |
回复此楼» 猜你喜欢
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有6人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
liu-biao
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 4598.4
- 散金: 84
- 帖子: 464
- 在线: 247.2小时
- 虫号: 2102358
- 注册: 2012-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
16楼2013-11-12 18:31:51
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
4楼2013-11-04 13:10:11
