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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC计算问题
171/1返回列表
查看: 1951  |  回复: 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英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬们，哪位大神知道E09有没有送啊？ +4 牛棕色左右 2024-05-15 6/300 2024-05-18 07:05 by Wing_spring
[教师之家] 今天才知道硕士考进大专比博士进本科工资高得多 +8 ZHONGWU_U 2024-05-14 9/450 2024-05-17 23:22 by 胡杨之殇
[基金申请] 申请基金代表性成果 +10 lancet0903 2024-05-17 10/500 2024-05-17 21:02 by lch2012
[教师之家] “直接受聘正高专业技术职务”怎么理解 +7 ZHONGWU_U 2024-05-17 8/400 2024-05-17 17:36 by jurkat.1640
[教师之家] 另一半是老师会加分哦 +11 acogw9004 2024-05-14 11/550 2024-05-17 11:18 by fysh2023
[教师之家] 南京林业大学老师首聘期考核不过，服药自杀 +30 babu2015 2024-05-12 35/1750 2024-05-17 10:55 by hunt_a_job
[有机交流] 无水硫酸钠除水 80+3 桃桃PXS 2024-05-13 5/250 2024-05-17 09:49 by FuMmm
[教师之家] 青椒的未来在哪里？ +14 sxhuabin 2024-05-15 14/700 2024-05-17 08:58 by Quakerbird
[考博] 申博 +3 己醛糖 2024-05-11 4/200 2024-05-17 07:33 by 化研迅
[论文投稿] 求期刊推荐 30+3 jxhx666666 2024-05-15 4/200 2024-05-16 22:47 by nono2009
[基金申请] 国自然青基把共同第一作者标成了唯一第一作者 +26 1137025804 2024-05-15 28/1400 2024-05-16 22:26 by 1135482512
[考博] 寻找高分子领域博导 +8 小政雅图 2024-05-14 11/550 2024-05-16 22:10 by 单调的胖子
[基金申请] 12个本子最多让给1a2b +16 地球e村长 2024-05-13 19/950 2024-05-16 17:48 by sars518
[考博] 急求博导收留 +9 8139173 2024-05-15 10/500 2024-05-16 13:06 by 北京炸酱面
[分析] 液相色谱基线上移 5+3 江上鸢飞 2024-05-13 3/150 2024-05-16 11:49 by gwmgyp
[考博] 深度学习申博 +4 酸菜鱼． 2024-05-15 4/200 2024-05-16 11:12 by 安塔瓦拉多
[硕博家园] 求助，请问我是否该转专业 +8 ygcbl 2024-05-12 10/500 2024-05-16 09:18 by 清晨和风一起
[基金申请] 2024工材口的国自然开始涵评了吗？ +3 Thanks1989 2024-05-13 3/150 2024-05-14 22:27 by weiyongscut
[论文投稿] 七个月了，还在selecting for review +3 g9522 2024-05-14 6/300 2024-05-14 19:00 by lizhengke06
[硕博家园] 哈工大硕博招生! +4 nailooo 2024-05-12 5/250 2024-05-13 22:04 by yuanjijoy
信息提示
关闭
请填处理意见
关闭