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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC计算问题
171/1返回列表
查看: 2585  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭