IRC的结果继续优化出现异常情况

zc000 · 发表于 Post on 2018-6-28 09:09:50

本帖最后由 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的动图看也是我期望的生成物，但是对右边最后一个点进行优化（方法和基组没有变），得到一个莫名的结构（本来是一个原子慢慢远离分子，结果优化出来这个原子又跑回去了，生成了一个奇怪的结构）。这是什么原因导致的。不知道我是否把问题表达清楚了？

zjxitcc · 发表于 Post on 2018-6-28 10:35:09

描述抽象，不如直接贴图

zc000 · 发表于 Post on 2018-6-28 14:55:32

本帖最后由 zc000 于 2018-7-4 09:29 编辑

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

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

sobereva · 发表于 Post on 2018-6-28 15:37:31

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] IRC的结果继续优化出现异常情况

浏览过的版块