24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1715  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求生物学学硕调剂——364分 +3 云朵遛弯指南 2026-04-04 3/150 2026-04-04 15:39 by dongzh2009
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 工科求调剂 +10 11ggg 2026-04-03 10/500 2026-04-03 22:35 by syjjj0321
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 282求调剂 不挑专业 求收留 +7 Yam. 2026-03-30 8/400 2026-04-03 14:12 by zhangdingwa
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 298分 070300求调剂 +11 zwen03 2026-04-02 11/550 2026-04-03 10:19 by macy2011
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 282求调剂 +3 aaa车辆 2026-04-02 3/150 2026-04-02 21:55 by zllcz
[考研] 环境科学与工程334分求调剂 +7 王一一依依 2026-03-30 9/450 2026-04-02 21:15 by 1104338198
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +9 真名有冰 2026-03-29 10/500 2026-03-31 18:01 by xhai2011
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
信息提示
请填处理意见