| 5 | 1/1 | 返回列表 |
| 查看: 4301 | 回復(fù): 12 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
cenwanglai榮譽(yù)版主 (知名作家)
|
[交流]
反應(yīng)路徑計算(NEB&CNEB)
|
||||
|
最近的一篇文章是關(guān)于反應(yīng)路徑的,算了幾條NEB,算了更多的CNEB. 文章投出前,趁熱寫出來,跟大家交流一下,歡迎批評指正。 (1)收斂速度 收斂速度有兩個標(biāo)準(zhǔn):開始計算到最后得到結(jié)果的時間;和開始計算到得到結(jié)果經(jīng)歷的總離子步數(shù)。我更關(guān)注第一個。 就我的金屬氧化物體系來看,IOPT = 7的fire算法30步內(nèi)有明顯的收斂效果,100步出現(xiàn)定性結(jié)果。覺得比其他方法都快。 (2)測試計算 反應(yīng)路徑或者過渡態(tài)計算失敗的機(jī)會很大。給出兩個端點,直接拉一條線,一個run到終點就收斂,這樣應(yīng)該是走大運(yùn)了。 NEB或者CNEB開始計算的時候多次出錯,不是一件壞事情。這可能意味著我們開始的設(shè)想是錯誤的?朔@些困難,就離真相更近了,研究的意義就體現(xiàn)出來。我覺得有幾點需要認(rèn)真考慮: a. 兩個端點是不是穩(wěn)定的。有可能精度不夠,或者該處“過于平坦”,收斂了,但是不是該領(lǐng)域的全局最小點。比如NEB中第二點或者第三點能量更低,但是他們與起點的幾何構(gòu)型很接近,那就果斷停掉NEB,從第二點或者第三點結(jié)構(gòu)優(yōu)化,找新的起始點。 b. 計算方法導(dǎo)致的出錯。有種情況,路徑初設(shè)施合理的,但是計算方法或者流程上不合理,導(dǎo)致計算過程中出現(xiàn)錯誤。這種情況最容易在CNEB計算中發(fā)生。一開始就高精度加CLIMBING,導(dǎo)致整個路徑嚴(yán)重偏離合理的范圍。這個時候,我不清楚,堅持算下去會不會恢復(fù)到合理范圍。我會選擇停掉計算,擦掉眼淚和失去的時間,從頭再來。避免這種錯誤出現(xiàn)的最好辦法,是多觀察中間結(jié)果,先從低精度快速收斂到一個較好的狀態(tài),保存中間結(jié)果慢慢推進(jìn)。 (3)結(jié)果分析 按照CNEB收斂的,最高點應(yīng)該就是過渡態(tài)了。可以計算虛頻驗證之。如果出現(xiàn)多過虛頻,那最可能的,是收斂精度不夠。一般我的EDIFFG=-0.03,不會超過兩個虛頻。 如果只有兩三個虛頻,分析一下各個頻率對應(yīng)的振動。如果是反應(yīng)路勁上只有一個虛頻,那應(yīng)該就是可以了。如果還有其他虛頻,應(yīng)該也是極小的。 有一種說法,不知對不對:如果虛頻對應(yīng)的是整體的平移,就忽略之。 關(guān)于零點校正,我不確定是不是一定需要。我見過沒有校正就發(fā)表的。也有不校正,校正,以及校正并考慮范德華力的過渡態(tài)計算。討論的時候到底看什么,我并沒有十足的把握。 我想,計算結(jié)果給出的能壘也只能是定性的,絕對值意義并不大。一般而言,一個體系里,關(guān)注點放在能量相對大小,得出定性結(jié)論就可以。 第一性原理過渡態(tài)能壘的絕對值到底有多大意義呢?希望各位發(fā)表意見。 [ Last edited by cenwanglai on 2012-7-26 at 19:35 ] |
計算化學(xué)-輔助工具 | VASP相關(guān)主題 |
» 搶金幣啦!回帖就可以得到:
+5/840
+1/188
+2/126
+1/70
+1/40
+2/36
+1/35
+3/22
+1/16
+1/16
+1/9
+1/9
+1/8
+1/6
+1/6
+1/6
+1/4
+1/3
+1/2
+1/1
榮譽(yù)版主 (知名作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑:085600材料與化工,考材科基,總分319 +16 | 678lucky 2026-03-31 | 20/1000 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +7 | Doleres 2026-03-31 | 7/350 |
|
|
[考研] 086000生物與醫(yī)藥 初試274求調(diào)劑 +4 | 小叮當(dāng)來了 2026-03-30 | 4/200 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +12 | 材料學(xué)257求調(diào)劑 2026-03-28 | 13/650 |
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +15 | Micky11223 2026-03-25 | 16/800 |
|
|
[考研] 土木304求調(diào)劑 +4 | 兔突突突, 2026-03-31 | 5/250 |
|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我!! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|