| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 4262 | 回復(fù): 12 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
cenwanglai榮譽版主 (知名作家)
|
[交流]
反應(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)該是走大運了。 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)了?梢杂嬎闾擃l驗證之。如果出現(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)主題 |
» 搶金幣啦!回帖就可以得到:
+2/162
+1/84
+1/79
+1/35
+1/34
+1/34
+1/31
+1/19
+1/18
+1/17
+1/17
+1/17
+1/16
+1/16
+1/13
+1/11
+1/6
+1/5
+1/1
+1/1
榮譽版主 (知名作家)
榮譽版主 (知名作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么! 2026-03-15 | 21/1050 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 8/400 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|