24小时热门版块排行榜    

查看: 2095  |  回复: 11
当前主题已经存档。

sophia_999

金虫 (小有名气)

[交流] 【求助】急问关于IRC=restart

我用IRC=restart计算,还没从chk里读构型就死掉了,是什么原因呢?下面是我的输入文件,请高手们帮忙,谢谢啦!

%chk=TS2.chk
# ub3lyp/6-311g(d,p) irc=(Restart,calcfc,reverse,maxpoint=150,stepsize=1)
guess=read geom=check  optcyc=500 scfcyc=500  test

ts2

0 2
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangwc5020

银虫 (小有名气)

讨论

★ ★
yjcmwgk(金币+1,VIP+0): 12-7 21:24
sophia_999(金币+1,VIP+0):谢谢参与 12-11 12:21
请问你的出错信息是?能贴出来嘛?
还有,我看你的命令行里是calcfc,我们计算IRC时,用的是过渡态的频率文件的chk,且是用rcfc读取力常数。
如果是因为chk的问题,无法读取结构,可以把频率文件的坐标拿出来算IRC,命令行写成现在的你的样子,就合理了
2楼2009-12-07 21:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophia_999

金虫 (小有名气)

下面是结果文件.
之前的irc走了一些点就死掉了,我接着读的原来做irc的chk。

# ub3lyp/6-311g(d,p) irc=(restart,calcfc,reverse,maxpoint=200,stepsize
=1) optcyc=500 scfcyc=500 iop(5/13=1) nosymm guess=read geom=check tes
t
----------------------------------------------------------------------
1/6=500,10=4,14=-1,18=40,22=2,26=3,29=120002,35=1,39=1,42=200/15;
2/15=1,29=3,40=1/2;
3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,74=-5/1,2,3;
4/5=5,7=2,16=3/1;
5/5=2,7=500,13=1,38=5/2;
8/6=4,10=90,11=11/1;
11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10;
10/6=1,7=6,31=1/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7/10=1,18=20,25=1,30=1/1,2,3,16;
1/6=500,10=4,14=-1,22=2,39=1,42=200/15(3);
2/15=1,29=3/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99//99;
2/15=1,29=3/2;
3/5=4,6=6,7=101,11=2,16=1,25=1,30=1,74=-5/1,2,3;
4/5=5,7=2,16=3/1;
5/5=2,7=500,13=1,38=5/2;
7/30=1/1,2,3,16;
1/6=500,14=-1,22=2,39=1,42=200/15(-5);
2/15=1,29=3/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Berny reaction coordinate optimization.
Restoring state from the checkpoint file "R1-ts12-ir.chk".
Title:  2009 11 18
Route:  # ub3lyp/6-311g(d,p) irc=(calcfc,reverse,maxpoint=200,stepsize
=1) test nosymm guess=read geom=check optcyc=500 scfcyc=500 iop(5/13=1
)
Number of steps in this run= 500 maximum allowed number of steps= 100.
Error termination via Lnk1e in /home/gauss03/g03/l115.exe at Fri Nov 27 17:16:14 2009.
3楼2009-12-08 10:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uminj

木虫 (著名写手)

★
sophia_999(金币+1,VIP+0):谢谢参与 12-10 09:42
楼主能解释一下 optcyc=500 scfcyc=500  关键字是干什么的吗? 我计算的时候没用过,你的问题可能就出在这。
4楼2009-12-08 13:18:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)

★
sophia_999(金币+1,VIP+0):谢谢参与 12-10 09:42
Number of steps in this run= 500 maximum allowed number of steps= 100.
这个就是要你把optcyc改到100以下
5楼2009-12-10 00:44:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophia_999

金虫 (小有名气)

前面的部分是9999死掉了,我就想加大点循环。把optcyc这个关键词去掉或减小循环次数仍然还是死了。
6楼2009-12-10 09:42:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophia_999

金虫 (小有名气)

引用回帖:
Originally posted by zhangwc5020 at 2009-12-7 21:11:
请问你的出错信息是?能贴出来嘛?
还有,我看你的命令行里是calcfc,我们计算IRC时,用的是过渡态的频率文件的chk,且是用rcfc读取力常数。
如果是因为chk的问题,无法读取结构,可以把频率文件的坐标拿出来算I ...

用rcfc也死掉了
7楼2009-12-10 09:44:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★
sophia_999(金币+1,VIP+0):谢谢参与 12-11 12:17
引用回帖:
Originally posted by uminj at 2009-12-8 13:18:
楼主能解释一下 optcyc=500 scfcyc=500  关键字是干什么的吗? 我计算的时候没用过,你的问题可能就出在这。

人家是G98标准写法,呵呵。你是从03开始用的吧
用optcyc=500 scfcyc=500  的都是老前辈了,小卒我可不敢说什么
8楼2009-12-10 10:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uminj

木虫 (著名写手)

★ ★ ★ ★ ★
yjcmwgk(金币+2,VIP+0): 12-11 09:23
yjcmwgk(金币+2,VIP+0): 12-11 11:26
sophia_999(金币+1,VIP+0):谢谢参与 12-11 12:22
引用回帖:
Originally posted by yjcmwgk at 2009-12-10 10:11:

人家是G98标准写法,呵呵。你是从03开始用的吧
用optcyc=500 scfcyc=500  的都是老前辈了,小卒我可不敢说什么

98,03我都用过
我计算时只有
irc=(Restart,calcfc,reverse,maxpoint=150,stepsize=1)
只要过渡态优化好后(通过gview查看振动方式),就等于成功了,一般情况下,irc都不会有什么问题。另外感觉你的步长有点小。
9楼2009-12-11 07:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophia_999

金虫 (小有名气)

引用回帖:
Originally posted by uminj at 2009-12-11 07:55:


98,03我都用过
我计算时只有
irc=(Restart,calcfc,reverse,maxpoint=150,stepsize=1)
只要过渡态优化好后(通过gview查看振动方式),就等于成功了,一般情况下,irc都不会有什么问题。另外感觉你的步长 ...

看虚频振动模式是没有问题的。我现在需要步长小步数多的irc计算结果。

我想问一下,你做restart得到正常结束的结果了么?我做的仍是死掉了
10楼2009-12-11 12:19:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sophia_999 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 4/200 2026-10-26 03:04 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3xN9SualEUPt 2026-10-24 6/300 2026-10-26 02:42 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:38 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:25 by cLOojpOxkX0G
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:03 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 01:50 by cLOojpOxkX0G
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 3/150 2026-10-26 01:37 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:26 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-26 00:47 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:45 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:41 by cLOojpOxkX0G
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-26 00:34 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:32 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-26 00:28 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:19 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
信息提示
请填处理意见