应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » NEB的计算出现负值，测试的其他能量都比初始的00和终止的10低
181/2返回列表12下一页
查看: 2649  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

身在长沙

新虫 (初入文坛)

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-04 12:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liu-biao

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
求楼下解决  我也遇到了 大神们 求支援
赞一下(1人) 回复此楼
2楼2013-11-04 12:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
identation: 金币+1, 谢谢交流 2013-11-04 13:53:17
你把03这个结构拿出来 做一下结构优化，看看是否能够回到00结构
不能的话，就是你初始结构没找对，也就是说没有找到能量最低结构
赞一下 回复此楼
3楼2013-11-04 12:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

身在长沙

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by emilyoyang at 2013-11-04 12:43:21
你把03这个结构拿出来 做一下结构优化，看看是否能够回到00结构
不能的话，就是你初始结构没找对，也就是说没有找到能量最低结构

谢谢你的回复，这个我做过了，弛豫之后又回到00结构了。而且我做到这个构型的能量最低位置与文献上的是一致的。我怀疑是参数设置那里出了问题，麻烦你再帮我看看。新手上路，请多多指教啊！
赞一下 回复此楼
4楼2013-11-04 13:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-04 16:22:20
参数没什么大的问题
NELM=200 默认的60就行了  不用这么大
IOPT=7 最好用 1  我一般用1 3
IMAGES=9 这么多干嘛 反应路径很长吗  9images的话计算量会很大
方便的话可以把你的initial和final POSCAR贴出来  我看看
赞一下 回复此楼
5楼2013-11-04 16:07:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjh-111

银虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流 2013-11-04 16:22:29
我的感觉是nebbarrier.pl的脚本的问题，因为我之前跑NEB的时候就发现脚本各种错误，我觉得只要你确定初终态是正确的，过渡态符合常理就可以了

[ 发自小木虫客户端 ]
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

身在长沙
6楼2013-11-04 16:13:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

身在长沙

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by emilyoyang at 2013-11-04 16:07:14
参数没什么大的问题
NELM=200 默认的60就行了  不用这么大
IOPT=7 最好用 1  我一般用1 3
IMAGES=9 这么多干嘛 反应路径很长吗  9images的话计算量会很大
方便的话可以把你的initial和final POSCAR贴出来  我看 ...

谢谢你的回复，我的初末位置弛豫用的INCAR都是如下的
SYSTEM =W-bcc
ENCUT = 400
ISTART = 0 ; ICHARG = 2
ISMEAR = 2
NSW = 100 ; IBRION =2
ISIF = 2
EDIFF = 1E-04 ; EDIFFG = -1E-02
LREAL= Auto
PREC = Accurate
虚心求指导！
赞一下 回复此楼
7楼2013-11-04 16:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

身在长沙

新虫 (初入文坛)
内容已删除
回复此楼
8楼2013-11-04 16:49:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjh-111

银虫 (小有名气)

【答案】应助回帖

引用回帖:
8楼: Originally posted by 身在长沙 at 2013-11-04 16:49:32
我的计算脚本都是师兄师姐用过留下来的，这让我怎么搞啊？还有这个脚本最新的你知道在哪里下载吗？  麻烦你了...

嗯，我都是自己写，你可以上官网下个最新的看看，要么你那个脚本贴出来给我看看
一般情况我觉得问题不大，很可能就是把负号当正号就行了

[ 发自小木虫客户端 ]
赞一下 回复此楼
9楼2013-11-04 18:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjh-111

银虫 (小有名气)

【答案】应助回帖

内容已删除
赞一下 回复此楼
10楼2013-11-04 18:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 身在长沙 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭