| 17 | 1/1 | 返回列表 |
| 查看: 2804 | 回復(fù): 16 | ||||
[交流]
【求助】IRC計(jì)算問題
|
||||
|
我用TS方法一反應(yīng)的過渡態(tài),只有一個(gè)虛頻,看著虛頻的振動(dòng)方向可以用,想用IRC驗(yàn)證一下,要怎么做?我用了#p b3lyp/6-31++g(d,p) irc=(calcfc,stepsize=5,maxpiont=5) 計(jì)算,為什么一算就死掉了? [ Last edited by lijie110 on 2011-3-25 at 10:36 ] |
Bao_S_J |
» 搶金幣啦!回帖就可以得到:
+1/87
+1/81
+1/78
+1/71
+1/69
+1/68
+1/41
+1/35
+1/21
+1/9
+1/6
+1/5
+1/3
+1/2
+1/2
+1/1
+1/1
+1/1
+1/1
+1/1
|
如果你保存了計(jì)算TS的chk文件,你可以讀取chk來計(jì)算IRC,這樣可以節(jié)省機(jī)時(shí)。建議分開算前后路徑: copy TS的chk文件分別為forward.chk和reverse.chk,然后 irc=(rcfc,forward,stepsize=20,maxpoints=10),步長取大些,我一般設(shè)為20,前后我一般各計(jì)算15個(gè)點(diǎn)。將得到的最后一個(gè)點(diǎn)全優(yōu)化,看是否能回到反應(yīng)物和產(chǎn)物。如果能回去說明得到的過渡態(tài)是對的。 |
|
本帖內(nèi)容被屏蔽 |
|
這是輸出的文件。多謝! 一開始就死掉了 1/10=5,14=-1,18=40,22=1,26=3,29=120007,38=1,39=20,42=10/1,15; 2/29=3,40=1/2; 3/5=1,6=6,7=1111,11=2,16=1,25=1,30=1,74=-5/1,2,3; 4/5=1/1; 5/5=2,38=6/2; 6/7=2,8=2,9=2,10=2,28=1/1; 7//1,2,3,16; 1/14=-1,22=1,39=20,42=10/15(1); 99//99; 2/29=3/2; 3/5=1,6=6,7=1111,11=2,16=1,25=1,30=1,74=-5/1,2,3; 4/5=5,16=3/1; 5/5=2,38=5/2; 7//1,2,3,16; 1/14=-1,22=1,39=20,42=10/15(-5); 2/29=3/2; 6/7=2,8=2,9=2,10=2,19=2,28=1/1; 99/9=1/99; Leave Link 1 at Fri Mar 25 15:01:12 2011, MaxMem= 0 cpu: 1.0 (Enter d:\G03W\l101.exe) FileIO operation on non-existent file. FileIO: IOper= 2 IFilNo(1)=-20502 Len= 50 IPos= 20000 Q= 20905776 dumping /fiocom/, unit = 1 NFiles = 14 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 826880 FType=2 FMxFil=10000 Number 0 501 502 598 989 991 992 993 Base 385536 303104 310784 333312 304128 307200 306688 302592 End 826880 304104 332834 333314 306628 310481 306693 302692 End1 826880 304128 333312 333824 306688 310784 307200 303104 Wr Pntr 385536 303104 310784 333312 304128 307200 306688 302592 Rd Pntr 385536 303104 310784 333312 304128 310481 306693 302692 Length 441344 1000 22050 2 2500 3281 5 100 Number 994 995 996 997 998 999 Base 20480 302080 21504 22016 20992 333824 End 20510 302090 21604 302022 21042 385076 End1 20992 302592 22016 302080 21504 385536 Wr Pntr 20480 302080 21504 22016 20992 333824 Rd Pntr 20510 302090 21604 302022 21042 335076 Length 30 10 100 280006 50 51252 dumping /fiocom/, unit = 2 NFiles = 1 SizExt = 0 WInBlk = 512 defal = F LstWrd = 66048 FType=2 FMxFil=10000 Number 0 Base 20480 End 66048 End1 66048 Wr Pntr 20480 Rd Pntr 20480 Length 45568 dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 66048 FType=2 FMxFil=10000 Number 0 Base 20480 End 66048 End1 66048 Wr Pntr 20480 Rd Pntr 20480 Length 45568 Error termination in NtrErr: NtrErr Called from FileIO. [ Last edited by lijie110 on 2011-3-25 at 15:09 ] |
|
不是,應(yīng)該是你的chk文件的問題,無法從chk文件讀取坐標(biāo)信息以及力常數(shù)信息。你能將input文件貼出來嗎?我想應(yīng)該是你的輸入文件有問題。還有,如果你做IRC計(jì)算時(shí)使用與計(jì)算TS相同的chk文件名,那你的chk文件已經(jīng)被改寫。所以建議當(dāng)做IRC計(jì)算時(shí)要保留TS的chk文件并copy TS的chk為IRC的。這樣即使IRC計(jì)算出錯(cuò)也不用再重新計(jì)算力常數(shù)了。 |
| 17 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +5 | 為學(xué)666 2026-03-16 | 5/250 |
|
|
[考研] 301求調(diào)劑 +3 | yy要上岸呀 2026-03-17 | 3/150 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 341求調(diào)劑 +3 | 番茄頭--- 2026-03-10 | 3/150 |
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考研] 420求調(diào)劑 +4 | 莫向外求11 2026-03-10 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | 柯淮然 2026-03-10 | 8/400 |
|
|
[考研] 調(diào)劑 +5 | 呵唔哦豁 2026-03-10 | 5/250 |
|