| 5 | 1/1 | 返回列表 |
| 查看: 1819 | 回復(fù): 8 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
1maomao金蟲 (初入文壇)
|
[求助]
用Gaussian03或09,MP2(full)/aug-cc-pvdz如何優(yōu)化HI?
|
||
|
用Gaussian優(yōu)化含I的分子總是出錯,請各位高人幫幫忙,看看怎么能優(yōu)化出來 gjf文件 # mp2(full)/gen opt(maxcycle=400) freq scf(maxcycle=400) Title Card Required 0 1 I -0.94029852 2.25373131 0.00000000 H -2.57029852 2.25373131 0.00000000 H 0 aug-cc-pvdz **** I 0 S 6 1.00 2.449790E+03 4.190000E-04 3.598080E+02 2.240000E-03 1.440580E+01 3.972230E-01 9.076320E+00 -9.322490E-01 2.088100E+00 9.371380E-01 1.034980E+00 3.920860E-01 S 6 1.00 2.449790E+03 1.750000E-04 3.598080E+02 1.057000E-03 1.440580E+01 1.690000E-01 9.076320E+00 -4.217930E-01 2.088100E+00 6.388640E-01 1.034980E+00 3.201150E-01 S 1 1.00 3.162840E-01 1.0000000 S 1 1.00 1.217190E-01 1.0000000 S 1 1.00 4.200000E-02 1.0000000 P 5 1.00 1.953010E+01 5.893400E-02 1.108820E+01 -2.309300E-01 2.715630E+00 6.648010E-01 1.204300E+00 4.506730E-01 3.399450E-01 2.898000E-02 P 5 1.00 1.953010E+01 -1.883600E-02 1.108820E+01 8.000600E-02 2.715630E+00 -3.066520E-01 1.204300E+00 -1.475940E-01 3.399450E-01 6.075060E-01 P 1 1.00 1.108810E-01 1.0000000 P 1 1.00 3.380000E-02 1.0000000 D 5 1.00 4.547650E+01 4.266000E-03 1.319280E+01 -1.362500E-02 4.227410E+00 3.097560E-01 1.942800E+00 5.097720E-01 8.397710E-01 2.974610E-01 D 1 1.00 3.000000E-01 1.0000000 D 1 1.00 1.191000E-01 1.0000000 **** out出錯信息 >>>>>>>>>> Convergence criterion not met. SCF Done: E(RHF) = -1631.30068077 A.U. after 401 cycles Convg = 0.4417D+00 -V/T = 9.3532 S**2 = 0.0000 Convergence failure -- run terminated. Error termination via Lnk1e in e:\G03W\l502.exe at Sun Jul 08 16:59:17 2012. Job cpu time: 0 days 0 hours 0 minutes 55.0 seconds. File lengths (MBytes): RWF= 11 Int= 0 D2E= 0 Chk= 4 Scr= 1 |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |

銅蟲 (初入文壇)
|
你這個是自洽場沒有收斂的問題,可以在本論壇搜一搜解決方法。 但剛看過類似的問題在如下的帖子里 http://www.gaoyang168.com/bbs/viewthread.php?tid=4696851&fpage=1 先可以在scf里加xqc等關(guān)鍵字試試,然后是再看看對碘原子的贗勢處理對不對,一般碘這樣的重原子容易出問題 |
|
SCF不能收斂,若你輸入正確則可能存在不穩(wěn)定態(tài),用ZMAT坐標或者目前的CART坐標改變一下結(jié)構(gòu)先跳到“阱”內(nèi)試一下,若還不行,若還不行 ,換MP2到其它更多包含關(guān)聯(lián)能方法嘗試一下,還不行就有可能是排斥態(tài)。 但從你的輸入看,似乎有問題,I元素你用的為aug-cc-pVDZ的ECP基組,你只輸入了價電子基組,對核電子的ECP沒有輸入,這樣不會收斂的,即使你做了上述的嘗試,把這段按照格式加在后邊應(yīng)該就沒問題了,注意空格和**** I 0 I-ECP 4 28 g-ul potential 1 2 1.00000000 0.00000000 s-ul potential 3 2 40.03337600 49.98964900 2 17.30057600 281.00655600 2 8.85172000 61.41673900 p-ul potential 4 2 15.72014100 67.41623900 2 15.20822200 134.80769600 2 8.29418600 14.56654800 2 7.75394900 28.96842200 d-ul potential 4 2 13.81775100 35.53875600 2 13.58780500 53.33975900 2 6.94763000 9.71646600 2 6.96009900 14.97750000 f-ul potential 4 2 18.52295000 -20.17661800 2 18.25103500 -26.08807700 2 7.55790100 -0.22043400 2 7.59740400 -0.22164600 |
金蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)279求調(diào)劑 +3 | 哈哈哈^_^ 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|
[考研] 367求調(diào)劑 +5 | 芋泥啵啵… 2026-03-28 | 5/250 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個楠 2026-03-25 | 5/250 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 壓國家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|