| 查看: 2876 | 回复: 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 初末位置的能量我多次验算都是最低的，不知道是什么原因，求大神们不吝指教我这个新手！谢谢啦 |
回复此楼» 猜你喜欢
366求调剂
已经有9人回复
-
材料工程085601，270求调剂
已经有37人回复
-
279学硕食品专业求调剂院校
已经有18人回复
-
290调剂生物0860
已经有31人回复
-
一志愿085802 323分求调剂
已经有13人回复
-
277求调剂
已经有23人回复
-
322求调剂，08工科
已经有4人回复
-
材料工程281还有调剂机会吗
已经有30人回复
-
化学070300 求调剂
已经有16人回复
-
化工学硕294分，求导师收留
已经有12人回复
【答案】应助回帖
|
我实在不是什么高手  ，才跑NEB没多久。。。nebbarrier.pl 里面能量有关部分是 $energy = `grep 'energy w' $directories[$i]/OUTCAR|tail -1`; # $dist = `grep 'NEB: distance' $directories[$i]/OUTCAR|tail -1`; $force = `grep 'NEB: projections' $directories[$i]/OUTCAR|tail -1`; if ($ssneb_flag =~ /T/){ $force = `grep 'NEBCELL: projections' $directories[$i]/OUTCAR|tail -1`; }; $energy =~ s/\s+$//g; @energy = split(/\s+/,$energy); $energy = $energy[@energy-1]; if ($i==0) { $energy0 = $energy; } $energy -= $energy0; 取的是OUTCAR里面energy(sigma->0) = 之后的数值，然后跟00文件夹下的OUTCAR里面的能量相减得到的值。两个脚本都没什么问题。你可以进各个文件夹用 grep energy OUTCAR 搜一搜，看看数据，然后用它的方法 grep 'energy w' OUTCAR，应该就大概知道是什么问题了 |
15楼2013-11-06 10:04:28
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
