| 10 | 1/1 | 返回列表 |
| 查看: 3174 | 回復(fù): 9 | |||||||||
beefly專家顧問 (職業(yè)作家)
地溝油冶煉專家
|
[交流]
Gaussian用冗余內(nèi)坐標(biāo)優(yōu)化大分子結(jié)構(gòu)時(shí)的問題 已有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默認(rèn)使用冗余內(nèi)坐標(biāo)做結(jié)構(gòu)優(yōu)化,計(jì)算中檢查使用的冗余坐標(biāo)能否構(gòu)成3N-L(L=6或5)個(gè)獨(dú)立內(nèi)坐標(biāo)空間,也就是檢查Wilson G矩陣是否有3N-L個(gè)非零特征值。由于產(chǎn)生冗余內(nèi)坐標(biāo)的算法本身存在缺陷,隨著分子尺寸不斷增加(通常幾百個(gè)原子以上),有些重要的弱化學(xué)鍵(如氫鍵)會被丟掉,導(dǎo)致G矩陣的非零特征值個(gè)數(shù)小于3N-L,就會遇到上面的出錯信息。在Gaussian09.d中,由于重新定義了冗余內(nèi)坐標(biāo)中的成鍵參數(shù),丟失的弱鍵可能更多,因此程序報(bào)錯的機(jī)會就更多一些。 解決方法: 1,自己寫個(gè)小程序,給定原子類型后,找出鍵長在一定范圍的所有的弱鍵。如果不包含在Gaussian產(chǎn)生的冗余內(nèi)坐標(biāo)中,就把它用opt=ModRedundant加到輸入文件里。 2,更好的辦法是改用直角坐標(biāo)做優(yōu)化,不僅不會遇到上述問題,還可以避免耗時(shí)的G矩陣對角化過程。 |
我學(xué)習(xí)計(jì)算的一些帖子 | 高斯細(xì)節(jié) | 經(jīng)驗(yàn) | Gaussian |

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

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

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

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

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 359求調(diào)劑 +3 | 王了個(gè)楠 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 081700 調(diào)劑 267分 +10 | 迷人的哈哈 2026-03-23 | 10/500 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 271求調(diào)劑 +4 | 生如夏花… 2026-03-22 | 4/200 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+4
|
糖葫蘆888ll 2026-03-22 | 9/450 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|