| 5 | 1/1 | 返回列表 |
| 查看: 4023 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
Gaussian計(jì)算的隨機(jī)性 已有8人參與
|
|||
|
同一個(gè)版本Gaussian軟件,在三臺不同的服務(wù)器上進(jìn)行優(yōu)化計(jì)算,同一個(gè)輸入文件,前27步的能量完全一樣,從28步后能量就開始不一樣了,難道Gaussian計(jì)算也有隨機(jī)性,請教各位高手是否也遇到類似的問題。 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 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
可能的原因非常多。例如根據(jù)內(nèi)存大小的不同自動使用不同積分儲存方式(在寫入/讀取時(shí)引入數(shù)值差異),并行計(jì)算時(shí)數(shù)據(jù)累積運(yùn)算順序存在的一定隨機(jī)性等等。實(shí)際上前27步也并非相同,只是在有限的數(shù)值輸出精度中沒顯示出來罷了,隨著累積運(yùn)算的進(jìn)行差異被逐漸放大到可感知的地步。 有興趣的話建議你看看J Comput Aided Mol Des (2008) 22:39–51。雖然是討論數(shù)值誤差對分子力學(xué)優(yōu)化影響的,但對理解這個(gè)問題也有益。 |
木蟲 (著名寫手)
瞌睡木蟲

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

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求化學(xué)調(diào)劑 +5 | wulanna 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +10 | leaflight 2026-03-24 | 10/500 |
|
|
[考研] 283求調(diào)劑(080500) +7 | A child 2026-03-27 | 7/350 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +7 | WW.' 2026-03-23 | 9/450 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經(jīng)歷,求調(diào)劑! +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|