| 5 | 1/1 | 返回列表 |
| 查看: 1105 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
ww1987鐵桿木蟲 (小有名氣)
|
[交流]
【求助】計算錯誤 幫忙看看 已有5人參與
|
||
|
我計算C80(D5d)構(gòu)型的分子 計算方法用B3LYP/6-31G(d)水平 算到一半出現(xiàn)這種現(xiàn)象 請幫忙! >>>>>>>>>> Convergence criterion not met. SCF Done: E(RB+HF-LYP) = -3048.38575494 A.U. after 129 cycles Convg = 0.7361D-06 -V/T = 2.0094 S**2 = 0.0000 KE= 3.020073824779D+03 PE=-3.325628686430D+04 EE= 1.408728660716D+04 Convergence failure -- run terminated. Error termination via Lnk1e in C:\G03W\l502.exe at Sat Jul 03 19:20:23 2010. Job cpu time: 0 days 10 hours 15 minutes 11.0 seconds. File lengths (MBytes): RWF= 897 Int= 0 D2E= 0 Chk= 98 Scr= 1 |
鐵桿木蟲 (正式寫手)
|
網(wǎng)上的建議: 1. 檢查是否有初始文件錯誤 在命令行中加入 %kJob L301 or %kJob L302 如果通過則一般初始文件ok。 常見初級錯誤: a. 自旋多重度錯誤 b. 變量賦值為整數(shù) c. 變量沒有賦值或多重賦值 d. 鍵角小于等于0度,大于等于180度 e. 分子描述后面沒有空行 f. 二面角判斷錯誤,造成兩個原子距離過近 g. 分子描述一行內(nèi)兩次參考同一原子,或參考原子共線 2. SCF(自洽場)不收斂 則一般是L502錯誤 省卻情況做64個cycle迭代(G03缺省128 cycles) a. 修改坐標(biāo),使之合理 b. 改變初始猜 Guess=Huckel 或其他的,看Guess關(guān)鍵詞。 c. 增加疊代次數(shù) SCFCYC=N (對小分子作計算時最好不要增加,很可能結(jié)構(gòu)不合理) d. iop(5/13=1)這樣忽略不收斂,繼續(xù)往下做。 3. 分子對稱性改變 a. 修改坐標(biāo),強(qiáng)制高對稱性或放松對稱性 b. 給出精確的、對稱性確定的角度和二面角。 如CH4的角度給到109.47122 c. 放松對稱性判據(jù) Symm=loose d. 不做對稱性檢查 iop(2/16=1) (最好加這個選項) iop(2/16=2) 則保持新的對稱性來計算 4. Opt時收斂的問題 a. 修改坐標(biāo),使之合理 b. 增加疊代次數(shù)optcyc=N 5. 優(yōu)化過渡態(tài),若勢能面太平緩,則不好找到。 iop(1/8=10) 默認(rèn)30(下一個結(jié)構(gòu)和該結(jié)構(gòu)的差別0.3Å),可改成10。如果每一步都要用到小的步長,應(yīng)該加opt(notrustupdate) 6. 在CI(組態(tài))方法中 如QCISD(T),CCSD(T),CID方法中,省卻最大循環(huán)50,若出錯(L913錯誤) 解決方法:#P QCISD(maxcyc=N) 注:N≤512 7. 優(yōu)化過渡態(tài) opt=TS (給出過渡態(tài)) opt=qst2 (給出反應(yīng)物和產(chǎn)物) opt=qst3 (給出反應(yīng)物和產(chǎn)物和過渡態(tài)) a. 用G03時的出錯 opt=ts 必須加FC (force constant) 寫法:opt=(TS, calcFc) or opt=(TS,calchffc)計算HF力常數(shù),對QCISD,CCSD等方法用; or opt=(TS,modRedundant) (最好寫這個) b. 如果計算采用QCISD計算(不好計算FC) 則寫為QCISD opt=(TS, calcHFFC) (用HF計算FC) 8. 無法寫大的Scratch文件RWF a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1 b. 改變計算方法 MP2=Direct可以少占硬盤空間 c. 限制最大硬盤 maxdisk=N GB,****MB,有些系統(tǒng)寫2GB會出錯,可以寫2000MB 9. FOPT出錯 原因是變量數(shù)與分子自由度數(shù)不相等。 可用POPT 或直接用OPT 10. 優(yōu)化過渡態(tài)只能做一個STEP 原因是負(fù)本征數(shù)目不對 添加 iop(1/11)=1或者noeigentest |
木蟲 (小有名氣)

鐵桿木蟲 (正式寫手)
金蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +12 | zhubinhao 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +3 | Feisty。 2026-03-28 | 7/350 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考博] 26申博 +3 | 加油沖! 2026-03-26 | 3/150 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|