 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4018 | 回復(fù): 20 | ||||||||
zhangfan2192金蟲(chóng) (著名寫(xiě)手)
|
[交流]
【其他】vasp neb 計(jì)算討論 已有11人參與
|
|||||||
|
希望各位用vasp的蟲(chóng)友能夠討論下neb計(jì)算中的注意事項(xiàng)： 不知道為什么，我在用neb計(jì)算過(guò)度態(tài)的時(shí)候總是感覺(jué)計(jì)算結(jié)果不是很理想，很難找到自己需要的結(jié)果，比如說(shuō)計(jì)算的8個(gè)image的能量是單調(diào)上升或下降的，而不是拋物線形狀，不知道大家在計(jì)算的時(shí)候遇到這樣的情況沒(méi)，如果遇到的話大家都是怎么解決的？相信現(xiàn)在用vasp軟件的不在少數(shù)，歡迎大家積極討論建議！ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 VASP | first principle | 第一性原理和電化學(xué) | vasp計(jì)算總結(jié) |
| 計(jì)算資料收集byFrank |
» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
-
341求調(diào)劑(一志愿湖南大學(xué)070300)
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- vasp計(jì)算過(guò)渡態(tài)問(wèn)題
已經(jīng)有13人回復(fù)
- 新手請(qǐng)教用NEB方法計(jì)算energy barrier.
已經(jīng)有15人回復(fù)
- vasp用neb計(jì)算遷移勢(shì)壘
已經(jīng)有11人回復(fù)
- vasp NEB 如何選取原胞大小 ？
已經(jīng)有10人回復(fù)
- vasp 過(guò)渡態(tài)計(jì)算 線性插值方法。
已經(jīng)有8人回復(fù)
- 【求助】vasp計(jì)算過(guò)渡態(tài)結(jié)果出錯(cuò)
已經(jīng)有17人回復(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ù)
- 【求助】NEB算過(guò)渡態(tài)出現(xiàn)致命錯(cuò)誤
已經(jīng)有5人回復(fù)
- 【求助成功】怎么樣提高NEB計(jì)算速度?
已經(jīng)有19人回復(fù)
- 【討論】NEB尋找過(guò)渡態(tài)專題討論
已經(jīng)有49人回復(fù)
uuv2010
榮譽(yù)版主 (職業(yè)作家)
- 1ST強(qiáng)帖: 6
- 應(yīng)助: 160 (高中生)
- 貴賓: 3.14
- 金幣: 3681.7
- 散金: 5802
- 紅花: 79
- 沙發(fā): 35
- 帖子: 4560
- 在線: 916小時(shí)
- 蟲(chóng)號(hào): 966216
- 注冊(cè): 2010-03-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
★ ★ ★
gavinliu7390(金幣+2):謝謝交流！ 2010-07-06 11:25:33
xiaohunhun(金幣+1):謝謝 2010-07-06 12:34:51
zhangfan2192(金幣+5): 2010-07-06 14:08:52
gavinliu7390(金幣+2):謝謝交流！ 2010-07-06 11:25:33
xiaohunhun(金幣+1):謝謝 2010-07-06 12:34:51
zhangfan2192(金幣+5): 2010-07-06 14:08:52
做NEB要保證初態(tài)和末態(tài)處于能量極小值點(diǎn)（局域的或是全局的），也就是初態(tài)和末態(tài)都是經(jīng)過(guò)優(yōu)化之后的，在這個(gè)基礎(chǔ)上再做中間的imagines可能就可以避免單調(diào)上升和下降了。個(gè)人觀點(diǎn)，歡迎指正。 [ Last edited by uuv2010 on 2010-7-6 at 10:23 ] |
3樓2010-07-06 10:20:14
lfhuang
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 39 (小學(xué)生)
- 金幣: 3576.1
- 散金: 2141
- 紅花: 37
- 帖子: 1432
- 在線: 681.4小時(shí)
- 蟲(chóng)號(hào): 505824
- 注冊(cè): 2008-02-17
- 性別: GG
- 專業(yè): 金屬材料表面科學(xué)與工程
★
lvjian8596(金幣+1):感謝交流！ 2010-07-07 12:58:39
zhangfan2192(金幣+5):謝謝 2010-07-07 20:13:51
lvjian8596(金幣+1):感謝交流！ 2010-07-07 12:58:39
zhangfan2192(金幣+5):謝謝 2010-07-07 20:13:51
|
為什么LCLIMB=.False.呢？Climbing-Image NEB 能夠幫助找到路徑中的能量最高點(diǎn)。 還有，如果路徑較長(zhǎng)，可能還需要增加images 數(shù)量來(lái)保證路徑收斂正確。 再則，如果樓主大概知道反應(yīng)勢(shì)壘，并且很小的話，每一個(gè)image能量的計(jì)算就需要特別精細(xì)了，因?yàn)橐稽c(diǎn)誤差都會(huì)導(dǎo)致結(jié)果和預(yù)期很不一樣。 |
8樓2010-07-07 12:21:50
tangxiaowei
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 8.8
- 紅花: 2
- 帖子: 221
- 在線: 142.2小時(shí)
- 蟲(chóng)號(hào): 888937
- 注冊(cè): 2009-10-30
- 專業(yè): 電化學(xué)
18樓2010-11-20 21:29:59
stractor
金蟲(chóng) (著名寫(xiě)手)
- 1ST強(qiáng)帖: 6
- 應(yīng)助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時(shí)
- 蟲(chóng)號(hào): 394087
- 注冊(cè): 2007-06-06
- 專業(yè): 半導(dǎo)體材料
2樓2010-07-05 22:16:51
zhangfan2192
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 1007.3
- 散金: 62
- 紅花: 3
- 帖子: 1170
- 在線: 171.3小時(shí)
- 蟲(chóng)號(hào): 648132
- 注冊(cè): 2008-11-07
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2010-07-06 14:10:44
fanchen021
鐵桿木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 7994.9
- 散金: 1844
- 紅花: 4
- 帖子: 776
- 在線: 569小時(shí)
- 蟲(chóng)號(hào): 773188
- 注冊(cè): 2009-05-18
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2010-07-07 10:30:40
zhangfan2192
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 1007.3
- 散金: 62
- 紅花: 3
- 帖子: 1170
- 在線: 171.3小時(shí)
- 蟲(chóng)號(hào): 648132
- 注冊(cè): 2008-11-07
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
|
INCAR文件如下： SYSTEM = Sul-Fe Startparameter for this Run: NWRITE = 2 verbosity write-flag LPETIM= F write-flag & timer # ISTART = 0 0-new 1-cont 2-samecut job ISPIN = 2 1-no, 2-yes INIWAV = 1 0-jellium WF; 1-random numbers (recomm.) Electronic Relaxation ! ENCUT = 400.00 #Ecut (eV) PREC = Medium * #precission: #Low/Medium/High NELM = 200 max. # of electronic SC steps NELMDL = 5 # of non-consistent steps at the beggining:>0 or <0 EDIFF = 1.0e-4 stopping-criterion for electronic upd. VOSKOWN = 1 ISMEAR = 1 SIGMA = 0.1 ! AMIX = 0.2 ! BMIX = 0.0001 Ionic relaxation EDIFFG = 1.0e-3 stopping-conditions for ionic relax.loop (EDIFF*10) NSW = 200 * # of steps for ionic update (def:0) IBRION = 2 * ionic relax: -1:no move; 0-MD;1-Newton;2-CG;3-damped;4-? ISIF = 2 * calc.stress:0-no;1-tr;2-7 yes; 3-change vol.4-change shape ISYM = 0 symmetry:0-OFF, 1-ON,2-ON+memory conserving for CD POTIM = 0.3 if IBRION=1,2,3:scaling constant forces;IBRION=0,time st. LPLANE = .TRUE. NPAR = 8 IALGO = 48 * algorithm: use only 8 (CG) or 48 (RMM-DIIS) LDIAG = .TRUE. sub-space diagonalisation GGA = 91 * xc-type: PB, PW (86), LM, 91 (Perdew-Wang 91) LWAVE = .F. * write WAVECAR LCHARG = .T. * write CHGCAR and CHG LVTOT = .F. * write the local potential LOCPOT LELF = .F. create ELFCAR file LORBIT = .F. create PROOUT LREAL = Auto ! NEB IMAGES = 8 ! Number of images in the band, excluding endpoints SPRING = -5.0 eV/A2 ! Spring force between images; negative value turns on nudging LCLIMB = .False. ! Flag to turn on the climbing image algorithm |
6樓2010-07-07 10:40:56
7樓2010-07-07 12:07:50
9樓2010-07-12 09:39:04
10樓2010-07-12 09:43:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhangfan2192 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
