应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC计算问题
北京石油化工学院2026年研究生招生接收调剂公告
171/1返回列表
查看: 2844  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lijie110

新虫 (小有名气)

[交流] 【求助】IRC计算问题

我用TS方法一反应的过渡态，只有一个虚频，看着虚频的振动方向可以用，想用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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-24 17:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

caipingl

木虫 (正式写手)
★
lijie110(金币+1):谢谢参与
lijie110(金币+10): 谢谢回复 送上金币 2011-03-30 15:42:04
如果你保存了计算TS的chk文件，你可以读取chk来计算IRC，这样可以节省机时。建议分开算前后路径：
copy TS的chk文件分别为forward.chk和reverse.chk，然后 irc=(rcfc,forward,stepsize=20,maxpoints=10)，步长取大些，我一般设为20，前后我一般各计算15个点。将得到的最后一个点全优化，看是否能回到反应物和产物。如果能回去说明得到的过渡态是对的。
赞一下(2人) 回复此楼
7楼2011-03-25 12:03:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangshash

铜虫 (小有名气)
引用回帖:
Originally posted by boylc789 at 2011-04-05 13:36:12:
iop，怎么加？

iop（5/13=1，1/11=1）  不加容易出错
赞一下(1人) 回复此楼
16楼2011-04-07 09:50:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dx357

铜虫 (初入文坛)
★
lijie110(金币+1):谢谢参与
引用回帖:
Originally posted by lijie110 at 2011-03-24 17:33:22:
我用TS方法一反应的过渡态，只有一个虚频，看着虚频的振动方向可以用，想用IRC验证一下，要怎么做？我用了#p b3lyp/6-31++g(d,p) ir=(calcfc,stepsize=5,maxpiont=5) 计算，为什么一算就死掉了？

关键词是irc=(.......)吧，是不是把irc，c弄丢了。呵呵
赞一下(1人) 回复此楼
2楼2011-03-25 10:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowerge

禁虫 (小有名气)
★
lijie110(金币+1):谢谢参与
本帖内容被屏蔽
3楼2011-03-25 10:17:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie110

新虫 (小有名气)
引用回帖:
Originally posted by flowerge at 2011-03-25 10:17:13:
关键词是irc

你给点步长很小，可以多做几个点，可以更好的看出趋势

我把stepsize,maxpiont都加大了，试了几种情况还是不行，能帮忙分析一下原因吗？
赞一下 回复此楼
4楼2011-03-25 10:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie110

新虫 (小有名气)
引用回帖:
Originally posted by dx357 at 2011-03-25 10:12:48:
关键词是irc=(.......)吧，是不是把irc，c弄丢了。呵呵

是我打错了，计算时写的是IRC，是什么原因呢？帮忙分析下...十分感谢！
赞一下 回复此楼
5楼2011-03-25 10:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

panyang811

新虫 (初入文坛)
★
lijie110(金币+1):谢谢参与
lijie110(金币+10): 送上金币 谢谢回复 2011-03-30 15:41:39
应该是maxpoint，另外可以加上guess=read geom=allcheck
赞一下(1人) 回复此楼
6楼2011-03-25 11:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caipingl

木虫 (正式写手)
你应该把出错信息贴出来，这样大家才好帮你分析原因啊。
赞一下 回复此楼
8楼2011-03-25 12:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie110

新虫 (小有名气)
引用回帖:
Originally posted by caipingl at 2011-03-25 12:04:44:
你应该把出错信息贴出来，这样大家才好帮你分析原因啊。

这是输出的文件。多谢！！
一开始就死掉了

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 ]
赞一下 回复此楼
9楼2011-03-25 15:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caipingl

木虫 (正式写手)
你读取chk文件？“FileIO operation on non-existent file.”显示读取文件不存在。
赞一下 回复此楼
10楼2011-03-26 07:26:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luntan_tjnu

木虫 (小有名气)
★
lijie110(金币+1):谢谢参与
forward或reverse 忘了写了吧
赞一下(1人) 回复此楼
11楼2011-03-26 16:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caipingl

木虫 (正式写手)
lijie110(金币+10): 多谢，给你金币表示感谢 2011-03-30 15:42:55
引用回帖:
Originally posted by lijie110 at 2011-03-25 15:03:31:
这是输出的文件。多谢！！
一开始就死掉了

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,3 ...

不是，应该是你的chk文件的问题，无法从chk文件读取坐标信息以及力常数信息。你能将input文件贴出来吗？我想应该是你的输入文件有问题。还有，如果你做IRC计算时使用与计算TS相同的chk文件名，那你的chk文件已经被改写。所以建议当做IRC计算时要保留TS的chk文件并copy TS的chk为IRC的。这样即使IRC计算出错也不用再重新计算力常数了。
赞一下 回复此楼
12楼2011-03-28 22:43:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)
★
lijie110(金币+1):谢谢参与
引用回帖:
Originally posted by lijie110 at 2011-03-25 15:03:31:
这是输出的文件。多谢！！
一开始就死掉了

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,3 ...

这蛋疼的错误信息...一看就知道是压根没读取到chk文件...
赞一下(1人) 回复此楼
13楼2011-03-31 16:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangshash

铜虫 (小有名气)
★
lijie110(金币+1):谢谢参与
貌似错的很多
没读CHECK文件，iop也没加，reverse，forward也没加。
先看书啦
赞一下(1人) 回复此楼
14楼2011-03-31 18:53:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boylc789

木虫 (著名写手)
★
lijie110(金币+1):谢谢参与
引用回帖:
Originally posted by liangshash at 2011-03-31 18:53:35:
貌似错的很多
没读CHECK文件，iop也没加，reverse，forward也没加。
先看书啦

iop，怎么加？
回复此楼
15楼2011-04-05 13:36:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
boylc78917楼
2011-04-07 21:58   回复  
引用回帖:
Originally posted by liangshash at 2011-04-07 09:50:10: iop（5/13=1，1/11=1） 不加容易出错

3q
相关版块跳转 我要订阅楼主 lijie110 的主题更新
171/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 青春裁为三截 2026-03-29 4/200 2026-03-29 15:01 by 唐沐儿
[考研] 291求调剂 +5 Y-cap 2026-03-29 6/300 2026-03-29 13:18 by mumin1990
[考研] 一志愿郑州大学，080500学硕，总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 各位老师好，我的一志愿为北京科技大学085601材料专硕 +8 Koxui 2026-03-28 8/400 2026-03-29 09:50 by laoshidan
[考研] 0703化学 +11 妮妮ninicgb 2026-03-27 11/550 2026-03-29 06:45 by 544594351
[考研] 数一英一271专硕（085401）求调剂，可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 346求调剂 一志愿070303有机化学 +3 萝卜炖青菜 2026-03-28 3/150 2026-03-28 14:11 by 唐沐儿
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 086000调剂 +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[考研] 一志愿211院校 344分 东北农业大学生物学学硕，求调剂 +5 丶风雪夜归人丶 2026-03-26 8/400 2026-03-27 19:22 by 丶风雪夜归人丶
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
[考研] 中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培 +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by 猫咪猫咪呀
[考研] 机械学硕310分，数一英一，一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 一志愿哈工大，085400，320，求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
信息提示
关闭
请填处理意见
关闭