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

鐵桿木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|