| 2 | 1/1 | 返回列表 |
| 查看: 835 | 回復(fù): 1 | |||
wqg1806木蟲 (著名寫手)
科研農(nóng)民工
|
[交流]
【求助】盼高手指點(diǎn),NEB計算
|
|
我在計算過渡態(tài)的時候插入了4個點(diǎn),計算收斂之后處理數(shù)據(jù)的時候發(fā)現(xiàn)找到的不是能量的鞍點(diǎn)。在這種情況下應(yīng)該怎么辦?是不是繼續(xù)優(yōu)化,找到鞍點(diǎn)?或者提高計算精度,去尋找鞍點(diǎn)?我把計算的精度設(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. 計算收斂之后處理的數(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不是可以自動找到能量鞍點(diǎn)的嗎? |

鐵桿木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 362求調(diào)劑 +14 | 西南交材料專碩3 2026-03-31 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 085600,材料與化工321分求調(diào)劑 +11 | 大饞小子 2026-03-28 | 11/550 |
|
|
[考研] 環(huán)境285分,過六級,求調(diào)劑 +6 | xhr12 2026-04-02 | 6/300 |
|
|
[考研] 286分調(diào)劑 +20 | Faune 2026-03-30 | 22/1100 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-04-02 | 4/200 |
|
|
[考研] 262求調(diào)劑 +5 | 勵志一定發(fā)文章 2026-04-02 | 6/300 |
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +18 | 深V宿舍吧 2026-03-30 | 19/950 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +14 | c297914 2026-03-28 | 14/700 |
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
[考研] 材料專碩調(diào)劑 +17 | 椰椰。 2026-03-29 | 17/850 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
[考研] 調(diào)劑 +5 | 好好讀書。 2026-03-28 | 7/350 |
|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-01 | 9/450 |
|
|
[考研] 土木304求調(diào)劑 +3 | 兔突突突, 2026-03-31 | 3/150 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 求佛 +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|