24小時熱門版塊排行榜

返回列表

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

wxcl8826

木蟲 (小有名氣)

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

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

初學(xué)Gaussian03常見出錯分析
最初級錯誤
1. 自旋多重度錯誤
2. 變量賦值為整數(shù)
3. 變量沒有賦值
4. 鍵角小于等于0度，大于等于180度
5. 分子描述后面沒有空行
6. 二面角判斷錯誤，造成兩個原子距離過近
7. 分子描述一行內(nèi)兩次參考同一原子,或參考原子共線
運行出錯
1. 自洽場不收斂  SCF
            a. 修改坐標，使之合理
            b. 改變初始猜  Guess
            c. 增加疊代次數(shù) SCFCYC=N
            d. iop（5/13=1）
2.  分子對稱性改變
            a. 修改坐標，強制高對稱性或放松對稱性
            b. 給出精確的、對稱性確定的角度和二面角
            c. 放松對稱性判據(jù)  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出錯原因是變量數(shù)與分子自由度數(shù)不相等。可用POPT 或直接用OPT
5. 優(yōu)化過渡態(tài)只能做一個STEP 原因是負本征數(shù)目不對添加 iop（1/11）=1
6. 組態(tài)相互作用計算中相關(guān)能疊代次數(shù)不夠，增加疊代次數(shù)
QCISD（Maxcyc=N）
Default.Rou設(shè)置
• 在Scratch文件夾中的Default.Rou文件中設(shè)置G03程序運行的省缺參數(shù)：
•    -M-  200MW
• -P- 4
• -#- MaxDisk=10GB
• -#- SCF=Conventional or Direct
• -#- MP2=NoDirect or  Direct
• -#- OPTCYC=200
• -#- SCFCYC=200
• -#- IOPs 設(shè)置如iop（2/16=1）
Default.Rou設(shè)置中的沖突
• 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　！　  去掉這個設(shè)置，CCSD的作業(yè)就能進行了。
因此，建議在Default設(shè)置中只設(shè)置，內(nèi)存，最大硬盤，等項
振動分析的關(guān)鍵詞
振動分析的選項
拉曼光譜                Raman
非諧性                   Anharmonic          勢能函數(shù)不是簡諧的
內(nèi)轉(zhuǎn)動，受阻轉(zhuǎn)子          Hinderedrotor       力常數(shù)過小的轉(zhuǎn)動
振-轉(zhuǎn)耦合                VibRot
振動圓二色性             VCD                      用于光活性分子
同位素效應(yīng)，溫度，壓力  ReadIsotopes  不影響頻率數(shù)值
算法選項：
解析頻率                Analytic
數(shù)值頻率                Numeric
二階導(dǎo)數(shù)數(shù)值       Enonly
FREQ的缺省設(shè)置是：有解析方式就計算解析頻率，不計算Raman性質(zhì)
MP3，MP4，CISD，CCSD，QCISD等方法沒有解析頻率，只能算數(shù)值頻率。
一般各種理論方法下，計算的頻率和實驗測得頻率之間有系統(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       用質(zhì)量數(shù)是2（程序讀取精確值）的同位素
O 18       用質(zhì)量數(shù)是18（程序讀取精確值）的同位素
            空行
Freq=（ReadFC，ReadIsotopes）
同一方法基組下計算不同同位素的力常數(shù)矩陣（Hessian）相同，因此從前面頻率計算的chechpoint文件里讀取力常數(shù)，計算同位素效應(yīng)的頻率和熱力學(xué)性質(zhì)，可以節(jié)省許多時間。
注：同種元素的不同同位素，幾何結(jié)構(gòu)是相同的，電子結(jié)構(gòu)也是相同的。
也可以用Freqcheck.exe工具做不同溫度，壓力，同位素的熱力學(xué)分析。
數(shù)值頻率問題
• 數(shù)值頻率計算不需要很大內(nèi)存和硬盤空間，但是精度較差，時間也長。
Freq=(Numeric,(Step=N))
Step用于指定微分步長，缺省為0.001*5 Angstrom
Freq=Enonly  計算解析一階導(dǎo)數(shù)。數(shù)值二階導(dǎo)數(shù)

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

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

[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 一志愿華東理工大學(xué)081700，初試分數(shù)271 +4	kotoko_ik 2026-03-23	5/250	2026-03-23 23:23 by 呆呆師姐
[考研] 材料調(diào)劑 +3	匹克i 2026-03-23	3/150	2026-03-23 23:18 by peike
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 考研化學(xué)308分求調(diào)劑 +7	你好明天你好 2026-03-23	8/400	2026-03-23 18:39 by macy2011
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴

24小時熱門版塊排行榜

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

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

» 猜你喜歡

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

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

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

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

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

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

★★★ 三星級,支持鼓勵

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

★ 一星級,一般

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

★★★ 三星級,支持鼓勵

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

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

★★★ 三星級,支持鼓勵

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

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

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

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