24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

身在长沙

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 816.5
散金: 736
红花: 8
帖子: 46
在线: 13.4小时
虫号: 2447869
注册: 2013-05-04
专业: 金属材料的微观结构

[求助] 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
初末位置的能量我多次验算都是最低的，不知道是什么原因，求大神们不吝指教我这个新手！谢谢啦

回复此楼