| 24小時(shí)熱門(mén)版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 1218 | 回復(fù): 3 | ||
左岸冰新蟲(chóng) (初入文壇)
|
[求助]
關(guān)于優(yōu)化收斂和頻率計(jì)算不收斂問(wèn)題 已有2人參與
|
|
我剛剛接觸高斯,最近做了一個(gè)25個(gè)原子的化合物,結(jié)構(gòu)優(yōu)化收斂是4個(gè)yes,如下: Maximum Force 0.000428 0.000450 YES RMS Force 0.000103 0.000300 YES Maximum Displacement 0.000127 0.001800 YES RMS Displacement 0.000032 0.001200 YES Predicted change in Energy=-2.977102D-07 Optimization completed. 但優(yōu)化好的結(jié)構(gòu)計(jì)算頻率只有一個(gè)yes,如下: Maximum Force 0.000914 0.000450 NO RMS Force 0.000294 0.000300 YES Maximum Displacement 0.041760 0.001800 NO RMS Displacement 0.008580 0.001200 NO Predicted change in Energy=-4.117907D-05 這個(gè)正常嗎?是優(yōu)化的結(jié)構(gòu)不對(duì)還是頻率計(jì)算有問(wèn)題呢?請(qǐng)高人指點(diǎn) |
銀蟲(chóng) (正式寫(xiě)手)
新蟲(chóng) (初入文壇)
版主 (著名寫(xiě)手)
|
Sobereva在帖子《量子化學(xué)計(jì)算中幫助幾何優(yōu)化收斂的常用方法》中說(shuō): RFO、GDIIS、GEDIIS這三種優(yōu)化方法在走步時(shí)都需要Hessian矩陣(力常數(shù)矩陣),然而計(jì)算Hessian矩陣是很耗時(shí)的。因此,在默認(rèn)情況下,Gaussian會(huì)通過(guò)價(jià)力場(chǎng)的方法近似估計(jì)出初始的Hessian矩陣,在每一步優(yōu)化中只精確計(jì)算梯度,利用梯度對(duì)原先Hessian矩陣進(jìn)行修正不斷得到新的Hessian矩陣。所以,在默認(rèn)情況下,優(yōu)化從頭到尾使用的Hessian矩陣都只是近似的(這也是為什么在優(yōu)化出來(lái)的結(jié)構(gòu)上做Freq往往會(huì)顯示還沒(méi)收斂,因?yàn)镕req用的是精確的Hessian矩陣)。當(dāng)Hessian矩陣離精確值偏差較大,就會(huì)造成收斂緩慢,或者始終不收斂。為解決這個(gè)問(wèn)題,可以用opt=calcfc,這會(huì)在優(yōu)化的第一步使用精確計(jì)算的Hessian矩陣,但是仍可能后續(xù)優(yōu)化過(guò)程中Hessian矩陣逐漸變得越來(lái)越不精確而依然收斂失敗。opt=calcall則不僅在第一步,在后續(xù)的每一步中也都精確計(jì)算Hessian矩陣,這使得很多優(yōu)化失敗的情況都能得到解決,優(yōu)化所需步數(shù)通常也會(huì)減少很多,而且能夠保證最終優(yōu)化結(jié)果準(zhǔn)確(因?yàn)樽罱K判斷是否收斂時(shí)是基于精確的Hessian矩陣所得結(jié)果),但代價(jià)是每一步計(jì)算量會(huì)很大。 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 0817調(diào)劑 +3 | 沒(méi)有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 085601專(zhuān)碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|