计算化学公社

标题: IRC的结果继续优化出现异常情况 [打印本页]

作者
Author: zc000 时间: 2018-6-28 09:09
标题: IRC的结果继续优化出现异常情况
本帖最后由 zc000 于 2018-6-28 14:41 编辑

老师，您好。我在做反应的时候，过渡态命令#p opt=(calcall,ts,noeigen,tight) geom=allcheck guess=read b3lyp/6-311g(d) scf=xqc。然后做IRC的命令是#p IRC=(calcall,maxpoint=30,stepsize=5,LQA) b3lyp/6-311g(d) guess=read geom=connectivity scf=xqc。出来的图也不错（一共56个点），IRC左边对应优化出来是我想要的反应物，右边从IRC的动图看也是我期望的生成物，但是对右边最后一个点进行优化（方法和基组没有变），得到一个莫名的结构（本来是一个原子慢慢远离分子，结果优化出来这个原子又跑回去了，生成了一个奇怪的结构）。这是什么原因导致的。不知道我是否把问题表达清楚了？

作者
Author: zjxitcc 时间: 2018-6-28 10:35
描述抽象，不如直接贴图

作者
Author: zc000 时间: 2018-6-28 14:55
本帖最后由 zc000 于 2018-7-4 09:29 编辑

zjxitcc 发表于 2018-6-28 10:35
描述抽象，不如直接贴图

IRC的左边优化出来是反应物，没有问题。从IRC最高点到最右边的变化趋势如下图1、2。即蓝色的原子逐渐远离，这也是我期望的结构。而当我拿最后一个点继续优化的时候，则出现图3的结构，该结构不是我期望的。

作者
Author: sobereva 时间: 2018-6-28 15:37
优化时候用opt(calcall,maxstep=5,notrust,maxcyc=100)以比较小的步长上限进行优化。如果到最后还是慢慢变成了最后一张图那样，说明在当前计算级别下，IRC右侧实际上就是非常长的，通过优化才真正走到了头

欢迎光临计算化学公社 (http://bbs.keinsci.com/)