| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 4274 | 回復(fù): 12 | |||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||||
cenwanglai榮譽(yù)版主 (知名作家)
|
[交流]
反應(yīng)路徑計(jì)算(NEB&CNEB)
|
||||
|
最近的一篇文章是關(guān)于反應(yīng)路徑的,算了幾條NEB,算了更多的CNEB. 文章投出前,趁熱寫(xiě)出來(lái),跟大家交流一下,歡迎批評(píng)指正。 (1)收斂速度 收斂速度有兩個(gè)標(biāo)準(zhǔn):開(kāi)始計(jì)算到最后得到結(jié)果的時(shí)間;和開(kāi)始計(jì)算到得到結(jié)果經(jīng)歷的總離子步數(shù)。我更關(guān)注第一個(gè)。 就我的金屬氧化物體系來(lái)看,IOPT = 7的fire算法30步內(nèi)有明顯的收斂效果,100步出現(xiàn)定性結(jié)果。覺(jué)得比其他方法都快。 (2)測(cè)試計(jì)算 反應(yīng)路徑或者過(guò)渡態(tài)計(jì)算失敗的機(jī)會(huì)很大。給出兩個(gè)端點(diǎn),直接拉一條線(xiàn),一個(gè)run到終點(diǎn)就收斂,這樣應(yīng)該是走大運(yùn)了。 NEB或者CNEB開(kāi)始計(jì)算的時(shí)候多次出錯(cuò),不是一件壞事情。這可能意味著我們開(kāi)始的設(shè)想是錯(cuò)誤的。克服這些困難,就離真相更近了,研究的意義就體現(xiàn)出來(lái)。我覺(jué)得有幾點(diǎn)需要認(rèn)真考慮: a. 兩個(gè)端點(diǎn)是不是穩(wěn)定的。有可能精度不夠,或者該處“過(guò)于平坦”,收斂了,但是不是該領(lǐng)域的全局最小點(diǎn)。比如NEB中第二點(diǎn)或者第三點(diǎn)能量更低,但是他們與起點(diǎn)的幾何構(gòu)型很接近,那就果斷停掉NEB,從第二點(diǎn)或者第三點(diǎn)結(jié)構(gòu)優(yōu)化,找新的起始點(diǎn)。 b. 計(jì)算方法導(dǎo)致的出錯(cuò)。有種情況,路徑初設(shè)施合理的,但是計(jì)算方法或者流程上不合理,導(dǎo)致計(jì)算過(guò)程中出現(xiàn)錯(cuò)誤。這種情況最容易在CNEB計(jì)算中發(fā)生。一開(kāi)始就高精度加CLIMBING,導(dǎo)致整個(gè)路徑嚴(yán)重偏離合理的范圍。這個(gè)時(shí)候,我不清楚,堅(jiān)持算下去會(huì)不會(huì)恢復(fù)到合理范圍。我會(huì)選擇停掉計(jì)算,擦掉眼淚和失去的時(shí)間,從頭再來(lái)。避免這種錯(cuò)誤出現(xiàn)的最好辦法,是多觀(guān)察中間結(jié)果,先從低精度快速收斂到一個(gè)較好的狀態(tài),保存中間結(jié)果慢慢推進(jìn)。 (3)結(jié)果分析 按照CNEB收斂的,最高點(diǎn)應(yīng)該就是過(guò)渡態(tài)了。可以計(jì)算虛頻驗(yàn)證之。如果出現(xiàn)多過(guò)虛頻,那最可能的,是收斂精度不夠。一般我的EDIFFG=-0.03,不會(huì)超過(guò)兩個(gè)虛頻。 如果只有兩三個(gè)虛頻,分析一下各個(gè)頻率對(duì)應(yīng)的振動(dòng)。如果是反應(yīng)路勁上只有一個(gè)虛頻,那應(yīng)該就是可以了。如果還有其他虛頻,應(yīng)該也是極小的。 有一種說(shuō)法,不知對(duì)不對(duì):如果虛頻對(duì)應(yīng)的是整體的平移,就忽略之。 關(guān)于零點(diǎn)校正,我不確定是不是一定需要。我見(jiàn)過(guò)沒(méi)有校正就發(fā)表的。也有不校正,校正,以及校正并考慮范德華力的過(guò)渡態(tài)計(jì)算。討論的時(shí)候到底看什么,我并沒(méi)有十足的把握。 我想,計(jì)算結(jié)果給出的能壘也只能是定性的,絕對(duì)值意義并不大。一般而言,一個(gè)體系里,關(guān)注點(diǎn)放在能量相對(duì)大小,得出定性結(jié)論就可以。 第一性原理過(guò)渡態(tài)能壘的絕對(duì)值到底有多大意義呢?希望各位發(fā)表意見(jiàn)。 [ Last edited by cenwanglai on 2012-7-26 at 19:35 ] |
計(jì)算化學(xué)-輔助工具 | VASP相關(guān)主題 |
» 搶金幣啦!回帖就可以得到:
+3/1778
+2/138
+1/81
+1/80
+1/79
+2/56
+1/44
+1/40
+1/33
+1/27
+1/20
+1/15
+5/15
+1/9
+1/4
+1/4
+1/3
+1/3
+1/3
+1/1
金蟲(chóng) (正式寫(xiě)手)
鐵桿木蟲(chóng) (著名寫(xiě)手)
榮譽(yù)版主 (知名作家)
木蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +5 | vv迷 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21 | 4/200 |
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤(rùn)采123 2026-03-18 | 11/550 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線(xiàn)與星海 2026-03-19 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|