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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC计算问题
北京石油化工学院2026年研究生招生接收调剂公告
171/1返回列表
查看: 2877  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +3 一样YWY 2026-04-01 3/150 2026-04-01 10:49 by Jaylen.
[考研] 265求调剂 +6 yelck 2026-04-01 7/350 2026-04-01 09:53 by yelck
[考研] 273求调剂 +14 李芷新1 2026-03-31 14/700 2026-04-01 09:32 by fangnagu
[考研] 一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业 +6 zju51 2026-03-31 8/400 2026-04-01 09:29 by JourneyLucky
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 0856调剂 +7 曲听筠 2026-03-30 7/350 2026-04-01 08:51 by JourneyLucky
[考研] 359求调剂 +8 邓邓邓书书 2026-03-25 8/400 2026-04-01 08:20 by herarysara
[考研] 一志愿华东理工大学，080500学硕，317分，求调剂 +9 s1145 2026-03-31 9/450 2026-03-31 18:47 by JourneyLucky
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[考研] 279求调剂 +12 j的立方 2026-03-29 12/600 2026-03-30 20:30 by dick_runner
[考研] 304材料求调剂 +5 钟llll 2026-03-26 5/250 2026-03-30 19:37 by 源_2020
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 压国家一区线，求导师收留，有恩必谢！ +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[有机交流] 高温高压反应求助 10+4 chibby 2026-03-25 4/200 2026-03-27 21:08 by BT20230424
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 机械学硕310分，数一英一，一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
信息提示
关闭
请填处理意见
关闭