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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC计算问题
171/1返回列表
查看: 1956  |  回复: 16
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 lijie110 的 12 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会

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返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 12/600 2024-06-01 22:36 by 嘟噜嘟1
[基金申请] 化学B02口青基 代表作都是什么水平的？向大佬求助 +8 arthas_007 2024-06-01 10/500 2024-06-01 22:05 by photocatalysis
[基金申请] 化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？ +13 salmon95 2024-05-30 22/1100 2024-06-01 21:22 by salmon95
[硕博家园] 各位同学能否分享一下实验室的学生劳务发放标准呀？ +18 ma3252788 2024-05-30 18/900 2024-06-01 17:02 by 落雨风霜
[硕博家园] 论大家对6070后普通教授导师的看法 +6 SNaiL1995 2024-05-28 10/500 2024-06-01 16:55 by 篮球易飞
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[基金申请] 数理的人才答辩通知发了吗 +7 yzy3327 2024-05-30 7/350 2024-06-01 09:15 by 撒哈拉西部
[考博] 申博求助?本硕双非一篇三区和四区去985工科非天坑专业有没有希望？ +4 LYS1200 2024-05-29 6/300 2024-06-01 08:50 by 学术渣渣神
[考博] 求生物学的博导收留 +4 申博成功功成博� 2024-06-01 4/200 2024-06-01 08:01 by 起名字哎
[教师之家] 选择 +6 459582015 2024-05-31 6/300 2024-06-01 07:32 by otani
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[有机交流] 苯酚和炔溴反应成醚 6+3 TONGMEIMEI 2024-05-26 12/600 2024-05-31 14:41 by TONGMEIMEI
[材料综合] 真空封石英管 北京 +4 dessha 2024-05-29 5/250 2024-05-30 16:40 by mpdfwxgui
[论文投稿] 审稿专家比较坚定的让补充实验，但实在没法补充实验，修回还有希望吗？ (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[论文投稿] 有没有老师需要发表论文 +4 金老师论文助理- 2024-05-29 4/200 2024-05-29 16:51 by liuyupu132
[基金申请] 信息学部函评结束了吗？ +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
[基金申请] 感觉自然基金限制通过比例就是有点扯，学学B口，化学学部，不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
信息提示
关闭
请填处理意见
关闭