24小时热门版块排行榜    

查看: 1686  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhangzihang

金虫 (初入文坛)

[求助] 求助:已经算了过渡态,想从TS出发做IRC分析去找反应复合物和产物复合物 已有1人参与

乙醚+OH的脱氢反应,已经通过TS的方法找到了过渡态,然后就想从过渡态做IRC分析找打对应的反应复合物(RC)和产物复合物(PC),
指定的maxpoints=100,stepsize=10,可是它算到第87步的时候就算不下去了,报错如下:

Cartesian Forces:  Max     0.116797105 RMS     0.024200336
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Error in corrector energy =          -0.0531605844
Magnitude of corrector gradient =     5.2275905845
Magnitude of analytic gradient =      0.1728249652
Magnitude of difference =             5.2448124915
Angle between gradients (degrees)=   94.7773
Pt 87 Step number  20 out of a maximum of  20
CORRECTOR INTEGRATION CONVERGENCE:
   Recorrection delta-x convergence threshold:    0.010000
   Delta-x Convergence NOT Met
Maximum number of corrector steps exceded.
Error termination via Lnk1e in D:\G09W\l123.exe at Sat May 16 22:35:27 2015.
Job cpu time:  0 days  0 hours 54 minutes 32.0 seconds.
File lengths (MBytes):  RWF=     74 Int=      0 D2E=      0 Chk=      5 Scr=      1

我把能算好的最后一个结构提取出来,对它进行结构的优化,希望能够优化到RC,可是每次都是失败的,提示的错误是:


GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
NTrRot=    -1 NTRed=    91 NAtoms=    17 NSkip=    46 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in D:\G09W\l103.exe at Sat May 16 22:59:08 2015.
Job cpu time:  0 days  0 hours 13 minutes 46.0 seconds.
File lengths (MBytes):  RWF=     77 Int=      0 D2E=      0 Chk=      3 Scr=      1

请问各位高手,我该怎么做才能找到我想要找的RC和PC呢?万分感谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

wy_001049

银虫 (初入文坛)

引用回帖:
3楼: Originally posted by zhangzihang at 2015-05-18 15:22:58
您好,我就是拿最后一个构型做了普通的opt,但是算的时候总是报错:
Error in internal coordinate system.
Error termination via Lnk1e in D:\G09W\l103.exe at Sat May 16 22:59:08 2015.

还有请问iop是什 ...

L103出错不用怕。你稍稍调一下构型,稍稍哈,键长键角你随意(例如Fe-O键1.6998,你可以给调成1.6999),在OPT,就可以了。
4楼2015-05-18 15:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

gauss98

禁虫 (正式写手)

感谢参与,应助指数 +1
本帖内容被屏蔽

2楼2015-05-18 12:39:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzihang

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by gauss98 at 2015-05-18 12:39:49
算到87步基本上已经能够确定反应物和产物的构型了

最后一个构型做普通的opt几何优化就能得到反应物或者产物

断掉的原因大多数都是 走到最后分子复合物作用力太弱,导致势能面平坦而不收敛,可以用一个iop直接 ...

您好,我就是拿最后一个构型做了普通的opt,但是算的时候总是报错:
Error in internal coordinate system.
Error termination via Lnk1e in D:\G09W\l103.exe at Sat May 16 22:59:08 2015.

还有请问iop是什么东西?我在高斯09的说明书上也没有看到它的具体说明,请问我可以在哪里看到呢?谢谢
3楼2015-05-18 15:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzihang

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by wy_001049 at 2015-05-18 15:46:30
L103出错不用怕。你稍稍调一下构型,稍稍哈,键长键角你随意(例如Fe-O键1.6998,你可以给调成1.6999),在OPT,就可以了。...

谢谢你啦
5楼2015-05-19 21:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +8 呼呼呼。。。。 2026-03-17 8/400 2026-03-18 09:07 by 无际的草原
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 10/500 2026-03-18 08:50 by xujiaoszu
[考研] 312求调剂 +5 陌宸希 2026-03-16 6/300 2026-03-18 08:00 by Iveryant
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 303求调剂 +3 睿08 2026-03-17 5/250 2026-03-17 22:01 by 睿08
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见