| 5 | 1/1 | 返回列表 |
| 查看: 1236 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
關(guān)于優(yōu)化收斂和頻率計(jì)算不收斂問題 已有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ì)算有問題呢?請(qǐng)高人指點(diǎn) |
版主 (著名寫手)
|
Sobereva在帖子《量子化學(xué)計(jì)算中幫助幾何優(yōu)化收斂的常用方法》中說: RFO、GDIIS、GEDIIS這三種優(yōu)化方法在走步時(shí)都需要Hessian矩陣(力常數(shù)矩陣),然而計(jì)算Hessian矩陣是很耗時(shí)的。因此,在默認(rèn)情況下,Gaussian會(huì)通過價(jià)力場的方法近似估計(jì)出初始的Hessian矩陣,在每一步優(yōu)化中只精確計(jì)算梯度,利用梯度對(duì)原先Hessian矩陣進(jìn)行修正不斷得到新的Hessian矩陣。所以,在默認(rèn)情況下,優(yōu)化從頭到尾使用的Hessian矩陣都只是近似的(這也是為什么在優(yōu)化出來的結(jié)構(gòu)上做Freq往往會(huì)顯示還沒收斂,因?yàn)镕req用的是精確的Hessian矩陣)。當(dāng)Hessian矩陣離精確值偏差較大,就會(huì)造成收斂緩慢,或者始終不收斂。為解決這個(gè)問題,可以用opt=calcfc,這會(huì)在優(yōu)化的第一步使用精確計(jì)算的Hessian矩陣,但是仍可能后續(xù)優(yōu)化過程中Hessian矩陣逐漸變得越來越不精確而依然收斂失敗。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ì)很大。 |
銀蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +6 | nhj_ 2026-03-25 | 6/300 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[有機(jī)交流] 考研調(diào)劑 +4 | watb 2026-03-26 | 4/200 |
|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-29 | 8/400 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|