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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » NEB的计算出现负值，测试的其他能量都比初始的00和终止的10低
181/2返回列表12下一页
查看: 2782  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 085600 英一数二272求调剂 5+4 vida_a 2026-03-01 10/500 2026-03-02 12:33 by vida_a
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 9/450 2026-03-02 12:04 by 52hz~~
[考研] 材料工程专硕283求调剂 5+4 ，！？ 2026-03-02 6/300 2026-03-02 11:07 by 黑！在干嘛
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 7/350 2026-03-02 10:38 by lature00
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[基金申请] 本子写完了，给DS兄弟看了，得了92分 +3 Doma 2026-03-01 7/350 2026-03-02 00:00 by jnzsy
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 0856求调剂285 +10 吕仔龙 2026-02-28 10/500 2026-03-01 21:37 by 公瑾逍遥
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
信息提示
关闭
请填处理意见
关闭