身在長(zhǎng)沙

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 816.5
散金: 736
紅花: 8
帖子: 46
在線: 13.4小時(shí)
蟲號(hào): 2447869
注冊(cè): 2013-05-04
專業(yè): 金屬材料的微觀結(jié)構(gòu)

[求助] NEB的計(jì)算出現(xiàn)負(fù)值，測(cè)試的其他能量都比初始的00和終止的10低

問(wèn)題如下：
1.計(jì)算NEB的方法：首先利用nebmake.pl產(chǎn)生初始終態(tài)結(jié)構(gòu)，用xcrysden查看了路徑?jīng)]問(wèn)題，然后對(duì)初始終態(tài)結(jié)構(gòu)進(jìn)行幾何結(jié)構(gòu)優(yōu)化，接著采用優(yōu)化后的能量構(gòu)型進(jìn)行NEB的計(jì)算，按照VASP manual上說(shuō)明，在INCAR中設(shè)置如下：
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產(chǎn)生的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
初末位置的能量我多次驗(yàn)算都是最低的，不知道是什么原因，求大神們不吝指教我這個(gè)新手！謝謝啦

回復(fù)此樓

» 猜你喜歡

304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)
330求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

進(jìn)行過(guò)渡態(tài)計(jì)算時(shí)link died 已經(jīng)有10人回復(fù)

1樓 2013-11-04 12:04:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

liu-biao

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 4598.4
散金: 84
帖子: 464
在線: 247.2小時(shí)
蟲號(hào): 2102358
注冊(cè): 2012-11-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

求樓下解決我也遇到了大神們求支援

贊一下(1人)

回復(fù)此樓

2樓2013-11-04 12:16:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時(shí)
蟲號(hào): 1630291
注冊(cè): 2012-02-20
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
identation: 金幣+1, 謝謝交流 2013-11-04 13:53:17

你把03這個(gè)結(jié)構(gòu)拿出來(lái) 做一下結(jié)構(gòu)優(yōu)化，看看是否能夠回到00結(jié)構(gòu)
不能的話，就是你初始結(jié)構(gòu)沒(méi)找對(duì)，也就是說(shuō)沒(méi)有找到能量最低結(jié)構(gòu)

贊一下

回復(fù)此樓

3樓2013-11-04 12:43:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

身在長(zhǎng)沙

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 816.5
散金: 736
紅花: 8
帖子: 46
在線: 13.4小時(shí)
蟲號(hào): 2447869
注冊(cè): 2013-05-04
專業(yè): 金屬材料的微觀結(jié)構(gòu)

引用回帖:

3樓: Originally posted by emilyoyang at 2013-11-04 12:43:21
你把03這個(gè)結(jié)構(gòu)拿出來(lái) 做一下結(jié)構(gòu)優(yōu)化，看看是否能夠回到00結(jié)構(gòu)
不能的話，就是你初始結(jié)構(gòu)沒(méi)找對(duì)，也就是說(shuō)沒(méi)有找到能量最低結(jié)構(gòu)

謝謝你的回復(fù)，這個(gè)我做過(guò)了，弛豫之后又回到00結(jié)構(gòu)了。而且我做到這個(gè)構(gòu)型的能量最低位置與文獻(xiàn)上的是一致的。我懷疑是參數(shù)設(shè)置那里出了問(wèn)題，麻煩你再幫我看看。新手上路，請(qǐng)多多指教�。�

贊一下

回復(fù)此樓

4樓2013-11-04 13:10:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時(shí)
蟲號(hào): 1630291
注冊(cè): 2012-02-20
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★
franch: 金幣+2, 謝謝回帖交流 2013-11-04 16:22:20

參數(shù)沒(méi)什么大的問(wèn)題
NELM=200 默認(rèn)的60就行了  不用這么大
IOPT=7 最好用 1  我一般用1 3
IMAGES=9 這么多干嘛反應(yīng)路徑很長(zhǎng)嗎  9images的話計(jì)算量會(huì)很大
方便的話可以把你的initial和final POSCAR貼出來(lái)  我看看

贊一下

回復(fù)此樓

5樓2013-11-04 16:07:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zjh-111

銀蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 439.2
紅花: 1
帖子: 71
在線: 7.9小時(shí)
蟲號(hào): 1830673
注冊(cè): 2012-05-23
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
franch: 金幣+2, 謝謝回帖交流 2013-11-04 16:22:29

我的感覺(jué)是nebbarrier.pl的腳本的問(wèn)題，因?yàn)槲抑芭躈EB的時(shí)候就發(fā)現(xiàn)腳本各種錯(cuò)誤，我覺(jué)得只要你確定初終態(tài)是正確的，過(guò)渡態(tài)符合常理就可以了

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

身在長(zhǎng)沙

6樓2013-11-04 16:13:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

身在長(zhǎng)沙

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 816.5
散金: 736
紅花: 8
帖子: 46
在線: 13.4小時(shí)
蟲號(hào): 2447869
注冊(cè): 2013-05-04
專業(yè): 金屬材料的微觀結(jié)構(gòu)

引用回帖:

5樓: Originally posted by emilyoyang at 2013-11-04 16:07:14
參數(shù)沒(méi)什么大的問(wèn)題
NELM=200 默認(rèn)的60就行了  不用這么大
IOPT=7 最好用 1  我一般用1 3
IMAGES=9 這么多干嘛反應(yīng)路徑很長(zhǎng)嗎  9images的話計(jì)算量會(huì)很大
方便的話可以把你的initial和final POSCAR貼出來(lái)  我看 ...

謝謝你的回復(fù)，我的初末位置弛豫用的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
虛心求指導(dǎo)！

贊一下

回復(fù)此樓

7樓2013-11-04 16:37:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

身在長(zhǎng)沙

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 816.5
散金: 736
紅花: 8
帖子: 46
在線: 13.4小時(shí)
蟲號(hào): 2447869
注冊(cè): 2013-05-04
專業(yè): 金屬材料的微觀結(jié)構(gòu)

內(nèi)容已刪除

回復(fù)此樓

8樓2013-11-04 16:49:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zjh-111

銀蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 439.2
紅花: 1
帖子: 71
在線: 7.9小時(shí)
蟲號(hào): 1830673
注冊(cè): 2012-05-23
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

引用回帖:

8樓: Originally posted by 身在長(zhǎng)沙 at 2013-11-04 16:49:32
我的計(jì)算腳本都是師兄師姐用過(guò)留下來(lái)的，這讓我怎么搞��？還有這個(gè)腳本最新的你知道在哪里下載嗎？麻煩你了...

嗯，我都是自己寫，你可以上官網(wǎng)下個(gè)最新的看看，要么你那個(gè)腳本貼出來(lái)給我看看
一般情況我覺(jué)得問(wèn)題不大，很可能就是把負(fù)號(hào)當(dāng)正號(hào)就行了

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

9樓2013-11-04 18:40:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zjh-111

銀蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 439.2
紅花: 1
帖子: 71
在線: 7.9小時(shí)
蟲號(hào): 1830673
注冊(cè): 2012-05-23
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

內(nèi)容已刪除

贊一下

回復(fù)此樓

10樓2013-11-04 18:42:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	10/500	2026-03-18 20:41 by Wangjingyue
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321