| 2 | 1/1 | 返回列表 |
| 查看: 834 | 回復(fù): 1 | |||
wqg1806木蟲 (著名寫手)
科研農(nóng)民工
|
[交流]
【求助】盼高手指點(diǎn),NEB計(jì)算
|
|
我在計(jì)算過(guò)渡態(tài)的時(shí)候插入了4個(gè)點(diǎn),計(jì)算收斂之后處理數(shù)據(jù)的時(shí)候發(fā)現(xiàn)找到的不是能量的鞍點(diǎn)。在這種情況下應(yīng)該怎么辦?是不是繼續(xù)優(yōu)化,找到鞍點(diǎn)?或者提高計(jì)算精度,去尋找鞍點(diǎn)?我把計(jì)算的精度設(shè)置以及得到的數(shù)據(jù)貼一下:SYSTEM =Nb O ENCUT =400 Ionic relaxation EDIFF = 1E-5 ! low accuracy EDIFFG =-0.01 !force #PREC = high #NELM = 100 ! maximum electronic steps #NELMIN = 4 ! minimum electronic steps #NELMDL = -8 ! update of charge for 3 steps (for continuation: 0) #IBRION = -1 ! use CG algorithm for ions #IBRION = 2 ! use CG algorithm for ions NSW =1550 ! No. of steps for ion #NSW =0 ! No. of steps for ion # new added lines: AMIX = 0.05 #AMIX = 0.1 BMIX = 0.0001 ! almost zero, but 0 will crash some versions AMIX_MAG = 0.3 #AMIX_MAG = 0.5 BMIX_MAG = 0.0001 ! almost zero, but 0 will crash some versions #LPARD=.TRUE. #NBMOD=-3 #EINT=-1 0 #ISIF=3 #ISYM=1 POTIM = 0.3 #add for NEB IMAGES=4 SPRING=-5.0 LCLIMB=.TRUE. LTANGENTOLD=.FALSE. LDNEB=.FALSE. IBRION=3 #IALGO = 48 ! RMM-DIIS algorithm for electrons ALGO = Fast LREAL =.TRUE. ! evaluate projection operators in real space #NPAR = 1 NSIM = 4 ! blocked algorithm update, four bands a time GGA = PE ! PBE DOS related values LVTOT = .TRUE. 計(jì)算收斂之后處理的數(shù)據(jù):3 1.171351 -0.650848 1.218286 3.05 1.1930247 -0.677556873901334 1.24581005302295 3.1 1.2146984 -0.704827288335876 1.27010368993822 3.15 1.2363721 -0.732589228234031 1.29116691074581 3.2 1.2580458 -0.760772678526208 1.30899971544572 3.25 1.2797195 -0.789307624142813 1.32360210403796 3.3 1.3013932 -0.818124050014252 1.33497407652251 3.35 1.3230669 -0.847151941070934 1.34311563289939 3.4 1.3447406 -0.876321282243264 1.34802677316859 3.45 1.3664143 -0.90556205846165 1.34970749733011 3.5 1.388088 -0.9348042546565 1.34815780538395 3.55 1.4097617 -0.963977855758219 1.34337769733011 3.6 1.4314354 -0.993012846697216 1.33536717316859 3.65 1.4531091 -1.0218392124039 1.32412623289939 3.7 1.4747828 -1.05038693780867 1.30965487652251 3.75 1.4964565 -1.07858600784194 1.29195310403796 3.8 1.5181302 -1.10636640743411 1.27102091544572 3.85 1.5398039 -1.1336581215156 1.24685831074581 3.9 1.5614776 -1.1603911350168 1.21946528993822 3.95 1.5831513 -1.18649543286814 1.18884185302295 4 1.604825 -1.211901 1.154988 所以,3不是鞍點(diǎn)?CNEB不是可以自動(dòng)找到能量鞍點(diǎn)的嗎? |

鐵桿木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒(méi)有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|