| 查看: 4016 | 回復: 15 | |||
[交流]
Gaussian計算的隨機性 已有8人參與
|
|
同一個版本Gaussian軟件,在三臺不同的服務器上進行優(yōu)化計算,同一個輸入文件,前27步的能量完全一樣,從28步后能量就開始不一樣了,難道Gaussian計算也有隨機性,請教各位高手是否也遇到類似的問題。 20 -2679.23298795 -2679.23298795 -2679.23298795 21 -2679.23319064 -2679.23319065 -2679.23319065 22 -2679.23339199 -2679.23339199 -2679.23339199 23 -2679.23355888 -2679.23355888 -2679.23355889 24 -2679.23367177 -2679.23367177 -2679.23367177 25 -2679.23374824 -2679.23374824 -2679.23374825 26 -2679.23381226 -2679.23381226 -2679.23381226 27 -2679.23384277 -2679.23384277 -2679.23384277 28 -2679.23388423 -2679.23388425 -2679.23388425 29 -2679.23390839 -2679.23390840 -2679.23390840 30 -2679.23399369 -2679.23399347 -2679.23399351 31 -2679.23397296 -2679.23417440 -2679.23413908 32 -2679.23400632 -2679.23419135 -2679.23414835 33 -2679.23413415 -2679.23422701 -2679.23419652 |
鐵桿木蟲 (職業(yè)作家)
站在計算化學入門的門檻上

木蟲 (著名寫手)
瞌睡木蟲

鐵桿木蟲 (職業(yè)作家)
站在計算化學入門的門檻上

木蟲 (著名寫手)
瞌睡木蟲
|
我的感覺是類似蝴蝶效應,就是在某一步計算出現(xiàn)了一個看似微不足道的差別,但這種差別不斷累積下來最后導致了能量的不盡一致。 不過樓主如果能繼續(xù)優(yōu)化完,可能最終的優(yōu)化結構和總能量也不會有太大差別。 我加的那個關鍵詞是讓計算積分時的誤差盡量小,看看能不能通過減小一點這方面的誤差來保持能量的相同。 另外,總能量的差別還不到0.0001 a.u.,不算很大的差別吧? [ Last edited by hairan on 2011-5-5 at 18:26 ] |

木蟲 (著名寫手)
未來國家凍涼

| 謝謝以上各位的講解,事實上,在一臺服務器上能夠優(yōu)化出最終的結構,需要130步,盡管有一個虛頻。但在另一臺服務器上,優(yōu)化到41步就停止了,結果出現(xiàn)很多的NaN等錯誤的信息,最后是Problem with the distance matrix.就停止優(yōu)化了。這才引起我的注意。我用的Gaussian版本是G03 E.01,以前沒有注意過這種情況。只是想測試一下幾臺服務器的運行速度,才發(fā)現(xiàn)這種情況,感覺不應該出現(xiàn)。 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
可能的原因非常多。例如根據內存大小的不同自動使用不同積分儲存方式(在寫入/讀取時引入數(shù)值差異),并行計算時數(shù)據累積運算順序存在的一定隨機性等等。實際上前27步也并非相同,只是在有限的數(shù)值輸出精度中沒顯示出來罷了,隨著累積運算的進行差異被逐漸放大到可感知的地步。 有興趣的話建議你看看J Comput Aided Mol Des (2008) 22:39–51。雖然是討論數(shù)值誤差對分子力學優(yōu)化影響的,但對理解這個問題也有益。 |
榮譽版主 (著名寫手)
阿超
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調劑 +7 | RuiFairyrui 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 316求調劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 藥學105500求調劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 330一志愿中國海洋大學 化學工程 085602 有讀博意愿 求調劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 070300化學求調劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 考研調劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 考研調劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 求調劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 321求調劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
[考研] 340求調劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 081700 調劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 一志愿哈工大,085400,320,求調劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 0854AI CV方向招收調劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 一志愿國科過程所081700,274求調劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|