24小时热门版块排行榜    

查看: 1685  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 zhangzihang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +9 小熊joy 2026-03-14 9/450 2026-03-18 04:33 by anny19840123
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-16 10/500 2026-03-18 00:01 by 星空星月
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 299求调剂 +4 △小透明* 2026-03-17 4/200 2026-03-17 20:09 by peike
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 274求调剂0856材料化工 +13 z2839474511 2026-03-11 14/700 2026-03-17 16:51 by share_joy
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见