| 10 | 1/1 | 返回列表 |
| 查看: 3171 | 回復: 9 | |||||||||
beefly專家顧問 (職業(yè)作家)
地溝油冶煉專家
|
[交流]
Gaussian用冗余內(nèi)坐標優(yōu)化大分子結(jié)構(gòu)時的問題 已有7人參與
|
|
出錯信息: GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Eigenvalue 1920 is 9.98D-07 should be greater than 0.000001 Eigenvector: A602 A605 A526 A529 A590 1 -0.23218 0.21568 -0.17301 0.17222 -0.16746 A684 A687 D849 A603 D896 1 0.13759 -0.13443 -0.13232 0.11672 0.11547 NTrRot= -1 NTRed= 3567 NAtoms= 642 NSkip= 1647 IsLin=F Error in internal coordinate system. 問題分析: Gaussian默認使用冗余內(nèi)坐標做結(jié)構(gòu)優(yōu)化,計算中檢查使用的冗余坐標能否構(gòu)成3N-L(L=6或5)個獨立內(nèi)坐標空間,也就是檢查Wilson G矩陣是否有3N-L個非零特征值。由于產(chǎn)生冗余內(nèi)坐標的算法本身存在缺陷,隨著分子尺寸不斷增加(通常幾百個原子以上),有些重要的弱化學鍵(如氫鍵)會被丟掉,導致G矩陣的非零特征值個數(shù)小于3N-L,就會遇到上面的出錯信息。在Gaussian09.d中,由于重新定義了冗余內(nèi)坐標中的成鍵參數(shù),丟失的弱鍵可能更多,因此程序報錯的機會就更多一些。 解決方法: 1,自己寫個小程序,給定原子類型后,找出鍵長在一定范圍的所有的弱鍵。如果不包含在Gaussian產(chǎn)生的冗余內(nèi)坐標中,就把它用opt=ModRedundant加到輸入文件里。 2,更好的辦法是改用直角坐標做優(yōu)化,不僅不會遇到上述問題,還可以避免耗時的G矩陣對角化過程。 |
我學習計算的一些帖子 | 高斯細節(jié) | 經(jīng)驗 | Gaussian |

木蟲 (著名寫手)
禁蟲 (正式寫手)
物化程序猿

|
本帖內(nèi)容被屏蔽 |
木蟲 (著名寫手)

專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |

送紅花一朵 |
本帖內(nèi)容被屏蔽 |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +4 | 超級伊昂大王 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +6 | holy J 2026-03-21 | 6/300 |
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調(diào)劑 +5 | niko- 2026-03-22 | 5/250 |
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 考研化學308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 北科281學碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 一志愿華中科技大學071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|