版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

emilyoyang

木蟲 (正式寫手)

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

引用回帖:

7樓: Originally posted by 身在長沙 at 2013-11-04 16:37:21
謝謝你的回復(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
...

貼出來 POSCAR看看 INCAR沒啥問題

贊一下

回復(fù)此樓

11樓2013-11-04 21:23:21

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

身在長沙

新蟲 (初入文壇)

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

引用回帖:

11樓: Originally posted by emilyoyang at 2013-11-04 21:23:21
貼出來 POSCAR看看 INCAR沒啥問題...

謝謝你的關(guān)注！我的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

贊一下

回復(fù)此樓

12樓2013-11-05 18:38:29

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

身在長沙

新蟲 (初入文壇)

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

引用回帖:

9樓: Originally posted by zjh-111 at 2013-11-04 18:40:42
嗯，我都是自己寫，你可以上官網(wǎng)下個最新的看看，要么你那個腳本貼出來給我看看
一般情況我覺得問題不大，很可能就是把負(fù)號當(dāng)正號就行了
...

謝謝你的持續(xù)關(guān)注，感覺您是個高手，麻煩你指導(dǎo)指導(dǎo)。謝謝我的腳本太大，我傳上來。麻煩你了。

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : nebbarrier.pl

2013-11-05 18:50:05, 3.56 K

附件 2 : nebmake.pl

2013-11-05 18:50:25, 5.2 K

13樓2013-11-05 18:50:04

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

身在長沙

新蟲 (初入文壇)

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

送紅花一朵

引用回帖:

6樓: Originally posted by zjh-111 at 2013-11-04 16:13:14
我的感覺是nebbarrier.pl的腳本的問題，因為我之前跑NEB的時候就發(fā)現(xiàn)腳本各種錯誤，我覺得只要你確定初終態(tài)是正確的，過渡態(tài)符合常理就可以了

我的腳本都穿上來了，麻煩大神一點時間幫我看看，不勝感激！

贊一下

回復(fù)此樓

14樓2013-11-05 18:51:58

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

zjh-111

銀蟲 (小有名氣)

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

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

引用回帖:

13樓: Originally posted by 身在長沙 at 2013-11-05 18:50:04
謝謝你的持續(xù)關(guān)注，感覺您是個高手，麻煩你指導(dǎo)指導(dǎo)。謝謝我的腳本太大，我傳上來。麻煩你了。...

我實在不是什么高手

，才跑NEB沒多久。。。
nebbarrier.pl 里面能量有關(guān)部分是

$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) = 之后的數(shù)值，然后跟00文件夾下的OUTCAR里面的能量相減得到的值。兩個腳本都沒什么問題。你可以進(jìn)各個文件夾用 grep energy OUTCAR 搜一搜，看看數(shù)據(jù)，然后用它的方法 grep 'energy  w' OUTCAR，應(yīng)該就大概知道是什么問題了

贊一下

回復(fù)此樓

15樓2013-11-06 10:04:28

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

liu-biao

木蟲 (正式寫手)

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
身在長沙: 金幣+85, ★★★★★最佳答案, 我最后發(fā)現(xiàn)確實是版本問題。換成5.2的就可以了。。。算了，金幣就給你吧 2013-11-12 18:35:19

版本問題你換個版本就行了

贊一下

回復(fù)此樓

16樓2013-11-12 18:31:51

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

身在長沙

新蟲 (初入文壇)

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

謝謝。我的問題也解決了換成5.2版本就行了。

贊一下

回復(fù)此樓

17樓2013-11-12 18:33:39

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

AV兄

木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 654.8
散金: 350
紅花: 9
帖子: 319
在線: 598.1小時
蟲號: 2129181
注冊: 2012-11-15
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

17樓: Originally posted by 身在長沙 at 2013-11-12 18:33:39
謝謝。我的問題也解決了換成5.2版本就行了。

是VASP換成5.2嗎？可是我的VASP本來就是5.2啊我算出來的過渡態(tài)能壘也比初態(tài)和末態(tài)低，怎么解決？

贊一下

回復(fù)此樓

歡迎光臨我的 AI 社區(qū) http://keepdl.org

18樓2015-10-10 21:00:35

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主身在長沙的主題更新

返回列表

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

[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 314求調(diào)劑 +8	無懈可擊的巨人 2026-03-12	8/400	2026-03-18 14:50 by haxia
[考研] 化工學(xué)碩306求調(diào)劑 +10	42838695 2026-03-12	10/500	2026-03-18 14:42 by haxia
[考研] 085601材料工程專碩求調(diào)劑 +6	慕寒mio 2026-03-16	6/300	2026-03-18 14:26 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant