24小时热门版块排行榜    

查看: 1611  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见