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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有27人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有237人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复
半透性容器进行高湿或者低湿研究已经有0人回复
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习催化剂设计方向已经有0人回复
法国里昂大学招生CSC 资助博士生已经有0人回复
法国里昂大学接收CSC 资助访问学者已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）已经有1人回复
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材已经有0人回复

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

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

uminj

木虫 (著名写手)

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

★ ★ ★ ★ ★
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都不会有什么问题。另外感觉你的步长有点小。

赞一下(21人)

回复此楼

9楼2009-12-11 07:55:00

已阅回复此楼关注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 个回答