24小時熱門版塊排行榜

返回列表

【獎勵】本帖被評價28次，作者wxcl8826增加金幣 25.05 個

wxcl8826

木蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1835.8
帖子: 269
在線: 26.3小時
蟲號: 700198

[資源] 【原創(chuàng)】GAUSSION計算常見錯誤及解決方案

初學Gaussian03常見出錯分析
最初級錯誤
1. 自旋多重度錯誤
2. 變量賦值為整數
3. 變量沒有賦值
4. 鍵角小于等于0度，大于等于180度
5. 分子描述后面沒有空行
6. 二面角判斷錯誤，造成兩個原子距離過近
7. 分子描述一行內兩次參考同一原子,或參考原子共線
運行出錯
1. 自洽場不收斂  SCF
            a. 修改坐標，使之合理
            b. 改變初始猜  Guess
            c. 增加疊代次數 SCFCYC=N
            d. iop（5/13=1）
2.  分子對稱性改變
            a. 修改坐標，強制高對稱性或放松對稱性
            b. 給出精確的、對稱性確定的角度和二面角
            c. 放松對稱性判據  Symm=loose
            d. 不做對稱性檢查 iop(2/16=1)
3. 無法寫大的Scratch文件RWF
            a. 劈裂RWF文件
%rwf=loc1,size1,loc2,size2,……..,locN,-1
            b. 改變計算方法 MP2=Direct可以少占硬盤空間
            c. 限制最大硬盤 maxdisk=N GB
4. FOPT出錯原因是變量數與分子自由度數不相等。可用POPT 或直接用OPT
5. 優(yōu)化過渡態(tài)只能做一個STEP 原因是負本征數目不對添加 iop（1/11）=1
6. 組態(tài)相互作用計算中相關能疊代次數不夠，增加疊代次數
QCISD（Maxcyc=N）
Default.Rou設置
• 在Scratch文件夾中的Default.Rou文件中設置G03程序運行的省缺參數：
•    -M-  200MW
• -P- 4
• -#- MaxDisk=10GB
• -#- SCF=Conventional or Direct
• -#- MP2=NoDirect or  Direct
• -#- OPTCYC=200
• -#- SCFCYC=200
• -#- IOPs 設置如iop（2/16=1）
Default.Rou設置中的沖突
• Default route:  MaxDisk=2GB SCF=Direct MP2=Direct OPTCYC=200 SCFcyc=100 iop(2/16=1) iop(5/13=1)
•    ------------------
•    # ccsd/6-31G** opt
•    ------------------
•    L903/L905 and L906 can only do MP2.
問題在于，ＭP2=Direct�。　�  去掉這個設置，CCSD的作業(yè)就能進行了。
因此，建議在Default設置中只設置，內存，最大硬盤，等項
振動分析的關鍵詞
振動分析的選項
拉曼光譜                Raman
非諧性                   Anharmonic          勢能函數不是簡諧的
內轉動，受阻轉子          Hinderedrotor       力常數過小的轉動
振-轉耦合                VibRot
振動圓二色性             VCD                      用于光活性分子
同位素效應，溫度，壓力  ReadIsotopes  不影響頻率數值
算法選項：
解析頻率                Analytic
數值頻率                Numeric
二階導數數值       Enonly
FREQ的缺省設置是：有解析方式就計算解析頻率，不計算Raman性質
MP3，MP4，CISD，CCSD，QCISD等方法沒有解析頻率，只能算數值頻率。
一般各種理論方法下，計算的頻率和實驗測得頻率之間有系統(tǒng)誤差，需要用校正因子校正（Scale Factor）
各種方法計算頻率的Scale因子
方法：          約化因子    約化因子
                     （頻率）    （ZPE）
HF/3-21G                0.9085                0.9409
HF/6-31G(d)             0.8929                0.9135
MP2(Full)/6-31G(d) 0.9427                0.9646
MP2(FC)/6-31G(d)    0.9434                0.9676
BLYP/6-31G(d)          0.9940                1.0119
B3LYP/6-31G(d)       0.9613                0.9804
SVWN/6-31G(d)       0.9833                1.0079
Freq=ReadIsotopes輸入
# 。。。Freq=ReadIsotopes 。。。
Title
0 2
H
O 1 0.98
298.15 1.5 0.8929    溫度，壓力約化因子
H 2       用質量數是2（程序讀取精確值）的同位素
O 18       用質量數是18（程序讀取精確值）的同位素
            空行
Freq=（ReadFC，ReadIsotopes）
同一方法基組下計算不同同位素的力常數矩陣（Hessian）相同，因此從前面頻率計算的chechpoint文件里讀取力常數，計算同位素效應的頻率和熱力學性質，可以節(jié)省許多時間。
注：同種元素的不同同位素，幾何結構是相同的，電子結構也是相同的。
也可以用Freqcheck.exe工具做不同溫度，壓力，同位素的熱力學分析。
數值頻率問題
• 數值頻率計算不需要很大內存和硬盤空間，但是精度較差，時間也長。
Freq=(Numeric,(Step=N))
Step用于指定微分步長，缺省為0.001*5 Angstrom
Freq=Enonly  計算解析一階導數。數值二階導數

振動分析中的問題
• 振動分析原則上只對穩(wěn)定點有意義
• 振動分析一定要在結構優(yōu)化相同的方法和基組上進行。
• 理論計算的振動頻率和試驗結果有系統(tǒng)誤差，需要乘以Scale因子校正
• 考慮相關能方法的解析頻率計算需要大的硬盤和內存
Warning -- explicit consideration of 3 degrees of freedom as vibrations may cause significant error
• 如果出現這個警告，一般來說是出現兩個方面的問題：
      1. 優(yōu)化的結構不夠精確，需要用更嚴格收斂判據優(yōu)化。
      2. 有3個振動模式實際上應該按內轉動處理。包括：自由轉子模型，受阻轉子模型。
內轉動模式對熱力學性質計算有影響，特別是低頻模式
穩(wěn)定點優(yōu)化中出現小的虛頻問題
1。優(yōu)化的結構還不是局域極小值點，需要用更嚴格的收斂判據優(yōu)化。OPT（tight  or verytight）
2。對應轉動模式的虛頻，可以用稍微修改轉動模式對應的二面角，使之合理。
3。用OPT=CalcAll來計算
結構優(yōu)化和頻率計算中的常見錯誤
1. 自洽場不收斂  SCF Convergence Failure
      a. 修改坐標，使之合理
      b. 改變初始猜  Guess
      c. 增加疊代次數 SCFCYC=N
      d.先用較松的SCF判據優(yōu)化結束后再用正常SCF判據優(yōu)化
      e.不管不收斂，繼續(xù)計算，iop（5/13=1）
      f. 嘗試用SCF=qc算法
2.  分子對稱性改變a. 修改坐標，強制高對稱性或放松對稱性b. 給出精確的、對稱性確定的角度和二面角如對Td對稱性，二面角值給為，109.47122 度c. 放松對稱性判據  Symm=loose d. 不做對稱性檢查 NoSymm或iop(2/16=1)
3. FOPT出錯  原因是變量數與分子自由度數不相等。可用POPT 或直接用OPT
4. 超出最大優(yōu)化圈數Step Exceed
• 檢查分子結構是否合理，如合理可增加循環(huán)圈數，OPTCYC=N
• 讀取Checkpoint文件最后幾何，繼續(xù)優(yōu)化。
• 選取上次優(yōu)化過程中最好的幾何結構繼續(xù)優(yōu)化。 GEOM=(step=n)
5. 無法寫大的Scratch文件RWF
a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1 b. 改變計算方法
c. 限制最大硬盤 maxdisk=N GB 6. 組態(tài)相互作用計算中相關能疊代次數不夠，增加疊代次數 QCISD(Maxcyc=N) CCSD (Maxcyc=N)等。

[ Last edited by fegg7502 on 2009-12-7 at 22:35 ]

[考研] 352求調劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 324求調劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 一志愿070300浙大化學358分，求調劑！ +3	酥酥魚.. 2026-03-21	3/150	2026-03-22 11:31 by 楊楊楊紫
[基金申請] 學校已經提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 265求調劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 22 350 本科985求調劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 機械專碩299求調劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 初始318分求調劑（有工作經驗） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 332求調劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 華東師范大學-071000生物學-293分-求調劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 353求調劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 招收調劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 一志愿中國海洋大學，生物學，301分，求調劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 材料專碩306英一數二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 考研求調劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 277調劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 考研調劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 333求調劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

24小時熱門版塊排行榜

[資源] 【原創(chuàng)】GAUSSION計算常見錯誤及解決方案

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★ 三星級,支持鼓勵

★★★★★ 五星級,優(yōu)秀推薦

★ 一星級,一般

★★★★★ 五星級,優(yōu)秀推薦

★★★ 三星級,支持鼓勵

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★ 三星級,支持鼓勵

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦

★★★★★ 五星級,優(yōu)秀推薦