| 24小時(shí)熱門(mén)版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 1217 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +4 | 吃吃吃才有意義 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|