 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2105 | 回復(fù): 6 | |||
cenwanglai榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
[交流]
【求助成功】NEB計(jì)算的一個stdout輸出怎么解讀? 已有3人參與
|
|
放假前一天把vststool Code和vasp4.6一起編譯了.在兩個端點(diǎn)中間插入3個images.下面是中間02文件夾中的stdout,OSZICAR及INCAR.希望得到幫助. 問題: (1)stdout中第一步后能量收斂就很快(很快到1E-7),說明什么?后面的quench提示是什么意思? (2)OSZICAR中電子步中能量是-400多,結(jié)果出來F=-1400多.這個是輸出錯誤吧?stdout中的離子步能量也是這樣的.是不vststool與vasp編譯的問題?OUTCAR中每離子步的能量與電子步收斂后的能量是相符的. stdout: vasp.4.6.28 25Jul05 complex 02/POSCAR found : 2 types and 55 ions LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro 00/POSCAR found : 2 types and 55 ions 04/POSCAR found : 2 types and 55 ions POSCAR, INCAR and KPOINTS ok, starting setup WARNING: wrap around errors must be expected FFT: planning ... 9 reading WAVECAR entering main loop N E dE d eps ncg rms rms(c) DAV: 1 0.554736476659E+04 0.55474E+04 -0.18260E+05 3200 0.138E+03 DAV: 2 0.154389023593E+04 -0.40035E+04 -0.37247E+04 3860 0.348E+02 .............................................................................................................. RMM: 31 -0.467372515298E+03 -0.41006E-03 -0.26808E-04 2111 0.242E-02 0.431E-02 RMM: 32 -0.467372569584E+03 -0.54286E-04 -0.10334E-04 1827 0.215E-02 1 F= -.14034855E+04 E0= -.14032665E+04 d E =-.140349E+04 quench: g(F)= 0.000E+00 g(S)= 0.000E+00 dE (1.order)= 0.000E+00 bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.467372620179E+03 -0.46737E+03 -0.14539E-04 2198 0.171E-02 0.109E-02 RMM: 2 -0.467372627169E+03 0.69903E-05 -0.48155E-05 1624 0.143E-02 0.168E-02 RMM: 3 -0.467372628565E+03 -0.13963E-05 -0.16158E-05 1509 0.936E-03 0.156E-02 RMM: 4 -0.467372628187E+03 0.37851E-06 -0.53661E-06 1176 0.382E-03 0.389E-03 RMM: 5 -0.467372628897E+03 -0.71051E-06 -0.92475E-07 750 0.114E-03 2 F= -.14034856E+04 E0= -.14032666E+04 d E =-.140349E+04 quench: g(F)= 0.000E+00 g(S)= 0.000E+00 dE (1.order)= 0.000E+00 bond charge predicted OSZICAR: N E dE d eps ncg rms rms(c) DAV: 1 0.554736476659E+04 0.55474E+04 -0.18260E+05 3200 0.138E+03 DAV: 2 0.154389023593E+04 -0.40035E+04 -0.37247E+04 3860 0.348E+02 DAV: 3 -0.316744761608E+03 -0.18606E+04 -0.18065E+04 4256 0.194E+02 ............................................................................................. RMM: 30 -0.467372105235E+03 -0.17295E-02 -0.16898E-03 2655 0.635E-02 0.330E-02 RMM: 31 -0.467372515298E+03 -0.41006E-03 -0.26808E-04 2111 0.242E-02 0.431E-02 RMM: 32 -0.467372569584E+03 -0.54286E-04 -0.10334E-04 1827 0.215E-02 1 F= -.14034855E+04 E0= -.14032665E+04 d E =-.140349E+04 N E dE d eps ncg rms rms(c) DAV: 1 -0.467372620179E+03 -0.46737E+03 -0.14539E-04 2198 0.171E-02 0.109E-02 RMM: 2 -0.467372627169E+03 -0.69903E-05 -0.48155E-05 1624 0.143E-02 0.168E-02 RMM: 3 -0.467372628565E+03 -0.13963E-05 -0.16158E-05 1509 0.936E-03 0.156E-02 RMM: 4 -0.467372628187E+03 0.37851E-06 -0.53661E-06 1176 0.382E-03 0.389E-03 RMM: 5 -0.467372628897E+03 -0.71051E-06 -0.92475E-07 750 0.114E-03 2 F= -.14034856E+04 E0= -.14032666E+04 d E =-.140349E+04 N E dE d eps ncg rms rms(c) DAV: 1 -0.467372629463E+03 -0.46737E+03 -0.17888E-06 2112 0.391E-03 0.181E-03 RMM: 2 -0.467372629846E+03 -0.38340E-06 -0.27115E-09 505 0.434E-06 0.261E-03 RMM: 3 -0.467372629831E+03 0.14523E-07 0.00000E+00 504 0.000E+00 0.131E-03 RMM: 4 -0.467372629813E+03 0.18692E-07 -0.28324E-09 505 0.402E-06 0.181E-03 RMM: 5 -0.467372629420E+03 0.39265E-06 -0.48393E-07 635 0.531E-04 3 F= -.14034856E+04 E0= -.14032665E+04 d E =-.140349E+04 INCAR: ENCUT = 480 ISTART= 0 ;ICHARG= 2 NPAR = 2 NELMIN = 5 EDIFF = 1E-4 EDIFFG = -0.05 LREAL = Auto NSW = 60 ALGO = Fast LPLANE = .TRUE. NSIM = 4 GGA = 91 POTIM = 0.0 ICHAIN = 0 LCLIMB = .TRUE. LGLOBAL = .TRUE. MAXMOVE = 0.2 INVCURV = 0.01 IMAGES = 3 SPRING = -5 IBRION = 3 ISMEAR=0 VOSKOWN=1 #PREC=Accurate [ Last edited by cenwanglai on 2010-7-1 at 20:45 ] |
回復(fù)此樓» 猜你喜歡
食品專碩 一志愿雙一流 328
已經(jīng)有4人回復(fù)
-
材料調(diào)劑
已經(jīng)有6人回復(fù)
-
300分，材料，求調(diào)劑，英一數(shù)二
已經(jīng)有5人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有13人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有6人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有6人回復(fù)
-
086003食品工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有16人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 新手請教用NEB方法計(jì)算energy barrier.
已經(jīng)有15人回復(fù)
- vasp用neb計(jì)算遷移勢壘
已經(jīng)有11人回復(fù)
- vasp NEB 如何選取原胞大小 ？
已經(jīng)有10人回復(fù)
- 【求助】請教CNEB
已經(jīng)有15人回復(fù)
- 【求助】NEB計(jì)算outcar中無各個鏡像的能量輸出
已經(jīng)有6人回復(fù)
- 【求助】關(guān)于CI-NEB
已經(jīng)有21人回復(fù)
- 【文獻(xiàn)學(xué)習(xí)】JACS上一篇純DFT計(jì)算的文獻(xiàn)： CI-NEB計(jì)算，Leuwdin布局分析，態(tài)密度分析
已經(jīng)有30人回復(fù)
- 【求助】關(guān)于vasp NEB
已經(jīng)有8人回復(fù)
- 【其他】vasp neb 計(jì)算討論
已經(jīng)有20人回復(fù)
- 【求助成功】怎么樣提高NEB計(jì)算速度?
已經(jīng)有19人回復(fù)
- 【求助/交流】請教N(yùn)EB的pmeI的酶切體系
已經(jīng)有7人回復(fù)
- 【討論】NEB尋找過渡態(tài)專題討論
已經(jīng)有49人回復(fù)
- 【求助】vasp的NEB計(jì)算,輸出需要哪些參數(shù)?
已經(jīng)有7人回復(fù)
2樓2010-06-17 12:31:31
cenwanglai
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
- 1ST強(qiáng)帖: 5
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 仿真模擬
|
請問OUTCAR文件中: (1)POSITION TOTAL-FORCE 這個是每個原子的坐標(biāo)及受力.對于非端點(diǎn)的image,所考察的移動原子在不同image中受力是不受EDIFFG=-xE-y(如-1E-2)約束的吧?根據(jù)NEB的原理,垂直與chain方向受力應(yīng)該->0,但是沿著chain方向是受力的.請問EDIFFG對哪些力有約束呢? (2)TANGET CHAIN-FORCE 跟(1)有對應(yīng)的六個數(shù),表示什么意思呢?分別表示切向受力和垂直與CHAIN方向受力? (3)前兩著之間有如下信息: energy of chain is (eV) 0.0000 for this image 0.000000 tangential force (eV/A) -2.644444 left and right image 0.6379 0.6379 A 這行說明中的文字部分是什么意思?數(shù)字表示什么物理意義? |
3樓2010-06-17 14:35:54
★ ★
hedaors(金幣+2):謝謝回復(fù) 2010-06-17 23:21:13
cenwanglai(金幣+6):呵呵,謝謝回復(fù)! 2010-06-18 08:42:12
cenwanglai(金幣+1):呵呵,有空請看看下面一條. 2010-06-18 09:52:40
hedaors(金幣+2):謝謝回復(fù) 2010-06-17 23:21:13
cenwanglai(金幣+6):呵呵,謝謝回復(fù)! 2010-06-18 08:42:12
cenwanglai(金幣+1):呵呵,有空請看看下面一條. 2010-06-18 09:52:40
|
(1)EDIFFG應(yīng)該是垂直方向上的力的收斂條件 (2)這個部分我算的輸出和你的有點(diǎn)不一樣 可能是版本更新了吧 前三個數(shù)應(yīng)該是TANGENT方向 CHAIN-FORCE不知道具體指什么 (3) 我也沒看到這部分輸出 energy of chain 不知道什么東西 而且是0 有點(diǎn)奇怪 tangential force就是切向力 left and right image是和旁邊image的距離 這個距離應(yīng)該差不多等距 不過climbing image兩邊是不等距離的 另外, 你可以在網(wǎng)頁上下載vtstscripts里面有nebef.pl這個script 可以直接輸出每個image上的力和能量 這個力就是收斂標(biāo)準(zhǔn)所用的力。 里面還有很多和neb相關(guān)的script |
4樓2010-06-17 23:18:36
cenwanglai
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
- 1ST強(qiáng)帖: 5
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 仿真模擬
5樓2010-06-18 09:51:58
6樓2010-06-18 10:48:59
劉小偉110
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1635.7
- 散金: 51
- 帖子: 256
- 在線: 120.5小時
- 蟲號: 2096954
- 注冊: 2012-10-30
- 專業(yè): 理論和計(jì)算化學(xué)
7樓2014-04-04 09:03:08
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cenwanglai 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +5 | @taotao 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工調(diào)劑 +8 | A-哆啦Z夢 2026-03-23 | 13/650 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6 | 小羊睡著了? 2026-03-23 | 9/450 |
|
|
[考研] 340求調(diào)劑 +3 | 話梅糖111 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +4 | 。，。，。，。i 2026-03-23 | 4/200 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +3 | 我愛學(xué)電池 2026-03-23 | 3/150 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|