| 17 | 1/1 | 返回列表 |
| 查看: 2813 | 回復(fù): 16 | ||||
[交流]
【求助】IRC計算問題
|
||||
|
我用TS方法一反應(yīng)的過渡態(tài),只有一個虛頻,看著虛頻的振動方向可以用,想用IRC驗證一下,要怎么做?我用了#p b3lyp/6-31++g(d,p) irc=(calcfc,stepsize=5,maxpiont=5) 計算,為什么一算就死掉了? [ Last edited by lijie110 on 2011-3-25 at 10:36 ] |
Bao_S_J |
» 搶金幣啦!回帖就可以得到:
+5/355
+1/165
+1/81
+1/80
+1/36
+1/32
+1/31
+1/31
+1/30
+1/26
+2/24
+1/17
+1/16
+1/7
+1/7
+1/5
+1/4
+1/4
+1/2
+1/1
|
如果你保存了計算TS的chk文件,你可以讀取chk來計算IRC,這樣可以節(jié)省機(jī)時。建議分開算前后路徑: copy TS的chk文件分別為forward.chk和reverse.chk,然后 irc=(rcfc,forward,stepsize=20,maxpoints=10),步長取大些,我一般設(shè)為20,前后我一般各計算15個點(diǎn)。將得到的最后一個點(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 ] |
| 17 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +7 | 煙汐憶海 2026-03-16 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 5/250 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +11 | 柴郡貓_ 2026-03-12 | 11/550 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|