24小时热门版块排行榜

>论坛更新日志 (3348)
>虫友互识 (626)
>导师招生 (271)
>文献求助 (195)
>考博 (177)
>休闲灌水 (177)
>论文投稿 (118)
>博后之家 (93)
>硕博家园 (82)
>考研 (80)
>教师之家 (47)
>招聘信息布告栏 (44)
>精细化工 (41)
>基金申请 (36)
>找工作 (27)
>公派出国 (27)

返回列表

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

身在长沙

新虫 (初入文坛)

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

回复此楼

» 猜你喜欢

版面费该交吗已经有9人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有13人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有8人回复
面上可以超过30页吧？已经有4人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有5人回复
什么是人一生最重要的？已经有4人回复
今年春晚有几个节目很不错，点赞！已经有12人回复

1楼 2013-11-04 12:04:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

身在长沙

新虫 (初入文坛)

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

引用回帖:

11楼: Originally posted by emilyoyang at 2013-11-04 21:23:21
贴出来 POSCAR看看 INCAR没啥问题...

谢谢你的关注！我的POSCAR是用计算出来的CONTCAR直接拷贝过来的具体如下：
W
1.00000000000000
   6.3299999237000000 0.0000000000000000 0.0000000000000000
   0.0000000000000000 6.3299999237000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 6.3299999237000000
  16 1
Direct
  0.0000000000000000 -0.0040008435774048 -0.0033687538950209
  0.0000000000000000 -0.0040008435774048  0.5033687538950209
  0.0000000000000000  0.5013936552370017 -0.0015767476815992
  0.0000000000000000  0.5013936552370017  0.5015767476815992
  0.5000000000000000  0.0008188547266570 -0.0006594841943691
  0.5000000000000000  0.0008188547266570  0.5006594841943691
  0.5000000000000000  0.5116831764744884 -0.0209389094477978
  0.5000000000000000  0.5116831764744884  0.5209389094477979
  0.2290610905522022  0.2383168235255117  0.2500000000000000
  0.2484232523184008  0.2486063447629983  0.7500000000000000
  0.2493405158056309  0.7491811452733428  0.2500000000000000
  0.2466312461049791  0.7540008435774049  0.7500000000000000
  0.7709389094477979  0.2383168235255117  0.2500000000000000
  0.7515767476815992  0.2486063447629983  0.7500000000000000
  0.7506594841943691  0.7491811452733428  0.2500000000000000
  0.7533687538950209  0.7540008435774049  0.7500000000000000
  0.5000000000000000  0.3750000000000000  0.2500000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00