版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sophia_999

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 2791.9
散金: 10
帖子: 187
在线: 75.4小时
虫号: 369787
注册: 2007-05-13
专业: 理论和计算化学

[交流] 【求助】急问关于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

回复此楼

» 猜你喜欢

请问如何鉴别TPU和普通PU材料已经有4人回复
汕头大学化学招调剂，相关内容如下。已经有4人回复
物理化学论文润色/翻译怎么收费? 已经有126人回复
求助火焰封管的时候管子炸了已经有1人回复
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？已经有2人回复
七嗪类物质合成求助已经有0人回复
半透性容器进行高湿或者低湿研究已经有0人回复
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材已经有0人回复

1楼 2009-12-07 11:34:43

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

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★
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 的都是老前辈了，小卒我可不敢说什么

赞一下(2人)

回复此楼

8楼2009-12-10 10:11:06

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

查看全部 12 个回答

zhangwc5020

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 335.8
散金: 28
帖子: 198
在线: 39.4小时
虫号: 683182
注册: 2008-12-28
专业: 理论和计算化学

讨论

★ ★
yjcmwgk(金币+1,VIP+0): 12-7 21:24
sophia_999(金币+1,VIP+0):谢谢参与 12-11 12:21

请问你的出错信息是？能贴出来嘛？
还有，我看你的命令行里是calcfc,我们计算IRC时，用的是过渡态的频率文件的chk，且是用rcfc读取力常数。
如果是因为chk的问题，无法读取结构，可以把频率文件的坐标拿出来算IRC，命令行写成现在的你的样子，就合理了

赞一下(11人)

回复此楼

2楼2009-12-07 21:11:40

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

sophia_999

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 2791.9
散金: 10
帖子: 187
在线: 75.4小时
虫号: 369787
注册: 2007-05-13
专业: 理论和计算化学

下面是结果文件.
之前的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

木虫 (著名写手)

应助: 7 (幼儿园)
金币: 6315.5
红花: 1
帖子: 1848
在线: 488小时
虫号: 886601
注册: 2009-10-28
专业: 橡胶及弹性体

★
sophia_999(金币+1,VIP+0):谢谢参与 12-10 09:42

楼主能解释一下 optcyc=500 scfcyc=500 关键字是干什么的吗？我计算的时候没用过，你的问题可能就出在这。

赞一下(2人)

回复此楼

4楼2009-12-08 13:18:11

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

查看全部 12 个回答