| 7 | 1/1 | 返回列表 |
| 查看: 1104 | 回復(fù): 6 | |||
ww1987鐵桿木蟲 (小有名氣)
|
[交流]
【求助】計(jì)算錯(cuò)誤 幫忙看看 已有5人參與
|
|
我計(jì)算C80(D5d)構(gòu)型的分子 計(jì)算方法用B3LYP/6-31G(d)水平 算到一半出現(xiàn)這種現(xiàn)象 請(qǐng)幫忙! >>>>>>>>>> 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 |
木蟲 (小有名氣)

至尊木蟲 (職業(yè)作家)
鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
|
網(wǎng)上的建議: 1. 檢查是否有初始文件錯(cuò)誤 在命令行中加入 %kJob L301 or %kJob L302 如果通過則一般初始文件ok。 常見初級(jí)錯(cuò)誤: a. 自旋多重度錯(cuò)誤 b. 變量賦值為整數(shù) c. 變量沒有賦值或多重賦值 d. 鍵角小于等于0度,大于等于180度 e. 分子描述后面沒有空行 f. 二面角判斷錯(cuò)誤,造成兩個(gè)原子距離過近 g. 分子描述一行內(nèi)兩次參考同一原子,或參考原子共線 2. SCF(自洽場(chǎng))不收斂 則一般是L502錯(cuò)誤 省卻情況做64個(gè)cycle迭代(G03缺省128 cycles) a. 修改坐標(biāo),使之合理 b. 改變初始猜 Guess=Huckel 或其他的,看Guess關(guān)鍵詞。 c. 增加疊代次數(shù) SCFCYC=N (對(duì)小分子作計(jì)算時(shí)最好不要增加,很可能結(jié)構(gòu)不合理) d. iop(5/13=1)這樣忽略不收斂,繼續(xù)往下做。 3. 分子對(duì)稱性改變 a. 修改坐標(biāo),強(qiáng)制高對(duì)稱性或放松對(duì)稱性 b. 給出精確的、對(duì)稱性確定的角度和二面角。 如CH4的角度給到109.47122 c. 放松對(duì)稱性判據(jù) Symm=loose d. 不做對(duì)稱性檢查 iop(2/16=1) (最好加這個(gè)選項(xiàng)) iop(2/16=2) 則保持新的對(duì)稱性來計(jì)算 4. Opt時(shí)收斂的問題 a. 修改坐標(biāo),使之合理 b. 增加疊代次數(shù)optcyc=N 5. 優(yōu)化過渡態(tài),若勢(shì)能面太平緩,則不好找到。 iop(1/8=10) 默認(rèn)30(下一個(gè)結(jié)構(gòu)和該結(jié)構(gòu)的差別0.3Å),可改成10。如果每一步都要用到小的步長(zhǎng),應(yīng)該加opt(notrustupdate) 6. 在CI(組態(tài))方法中 如QCISD(T),CCSD(T),CID方法中,省卻最大循環(huán)50,若出錯(cuò)(L913錯(cuò)誤) 解決方法:#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時(shí)的出錯(cuò) opt=ts 必須加FC (force constant) 寫法:opt=(TS, calcFc) or opt=(TS,calchffc)計(jì)算HF力常數(shù),對(duì)QCISD,CCSD等方法用; or opt=(TS,modRedundant) (最好寫這個(gè)) b. 如果計(jì)算采用QCISD計(jì)算(不好計(jì)算FC) 則寫為QCISD opt=(TS, calcHFFC) (用HF計(jì)算FC) 8. 無法寫大的Scratch文件RWF a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1 b. 改變計(jì)算方法 MP2=Direct可以少占硬盤空間 c. 限制最大硬盤 maxdisk=N GB,****MB,有些系統(tǒng)寫2GB會(huì)出錯(cuò),可以寫2000MB 9. FOPT出錯(cuò) 原因是變量數(shù)與分子自由度數(shù)不相等。 可用POPT 或直接用OPT 10. 優(yōu)化過渡態(tài)只能做一個(gè)STEP 原因是負(fù)本征數(shù)目不對(duì) 添加 iop(1/11)=1或者noeigentest |
新蟲 (小有名氣)
金蟲 (正式寫手)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|