计算化学公社

标题: 能量振荡无法收敛，已经很长时间了，实在没办法解决，各种办法都试过了，求助 [打印本页]

作者
Author: malong 时间: 2019-5-7 18:15
标题: 能量振荡无法收敛，已经很长时间了，实在没办法解决，各种办法都试过了，求助
本帖最后由 malong 于 2019-5-7 23:02 编辑

在寻找过渡态的时候，能量在两个结构之间振荡，能量最低的那个结构就是我要找的ts，但是由于振荡的原因无法收敛。通过在opt里面设置最大步长maxstep=2，能量依然振荡。这个能量优化图示用qst3的方法优化的

作者
Author: sobereva 时间: 2019-5-7 18:58
看轨迹变化，如果末端结构看着像TS，加上maxstep=3再试
如果末端结构都没TS的样子，重新调整初猜

以后传图时候每张图的含义都要进行解释，对问题进行确切的描述，别光发图

作者
Author: malong 时间: 2019-5-7 19:05
末端结构就是TS，我之前用maxstep=5还是振荡，我已经在用maxstep=2试了，图片描述的问题下次一定会注意的，非常感谢老师。

作者
Author: malong 时间: 2019-5-7 20:03
本帖最后由 malong 于 2019-5-7 22:27 编辑

sobereva 发表于 2019-5-7 18:58
看轨迹变化，如果末端结构看着像TS，加上maxstep=3再试
如果末端结构都没TS的样子，重新调整初猜

这个是我把qst3（帖子中的结构）振荡的最后一个结构提取出来然后用ts的方法继续寻找过渡态，第36个点是我设置优化步长maxstep=2从上次优化的结构继续restart，结果第37,38,39依然出现了振荡。
输入文件：%chk=ts-addition-enediyne.chk
%mem=120GB
%nprocshared=16
#p opt=(calcall,ts,maxstep=2,notrust,noeigentest,restart) freq b3lyp/6-31g(d) geom=check empiricaldispersion=gd3

作者
Author: alonewolfyang 时间: 2019-5-7 21:53

malong 发表于 2019-5-7 20:03
老师
这个是我把qst3（帖子中的结构）振荡的最后一个结构提取出来然后用ts的方法继续寻找过渡态，第36个 ...

你在搞笑吗，这坨图片谁认识

作者
Author: malong 时间: 2019-5-7 22:28

alonewolfyang 发表于 2019-5-7 21:53
你在搞笑吗，这坨图片谁认识

老师第一次发帖没有考虑全面望见谅，我把帖子重新整理了一遍，您看现在可以不

作者
Author: malong 时间: 2019-5-7 22:33
本帖最后由 malong 于 2019-5-7 22:49 编辑

alonewolfyang 发表于 2019-5-7 21:53
你在搞笑吗，这坨图片谁认识

我把帖子重新整理了，第一次发帖没有了解清楚规则，望见谅

作者
Author: Frank 时间: 2019-5-8 08:15

sobereva 发表于 2019-5-7 18:58
看轨迹变化，如果末端结构看着像TS，加上maxstep=3再试
如果末端结构都没TS的样子，重新调整初猜

借楼请教一下老师，如果用opt=(calcall,ts,noeigen)优化过渡态震荡不收敛，最后一帧结构接近过渡态结构，振动模式也对，还有必要尝试maxstep=3，notrust吗？还是说calcall优化不出来的用maxstep=3也不会优化出来？

作者
Author: sobereva 时间: 2019-5-8 09:28

Frank 发表于 2019-5-8 08:15
借楼请教一下老师，如果用opt=(calcall,ts,noeigen)优化过渡态震荡不收敛，最后一帧结构接近过渡态结构， ...

有必要
如果说calcall能解决80%概率的不收敛问题，加上maxstep可能能提升到90%~95%概率

作者
Author: sobereva 时间: 2019-5-8 09:30

malong 发表于 2019-5-7 20:03
这个是我把qst3（帖子中的结构）振荡的最后一个结构提取出来然后用ts的方法继续寻找过渡态，第36个点是我设 ...

我看不出你这是个什么反应。怎么碳上一个氢还跑到另一个氢上去了
当你要求优化TS，然而当前结构和实际TS偏离极大，都没正确TS的模样，反复震荡不收敛是正常的（哪怕收敛了，结果也没有意义）
PS：opt=TS的时候默认就用notrust，此时没必要写notrust

作者
Author: malong 时间: 2019-5-8 09:45

sobereva 发表于 2019-5-8 09:30
我看不出你这是个什么反应。怎么碳上一个氢还跑到另一个氢上去了
当你要求优化TS，然而当前结构和实际TS ...

我这个是三甲胺加氢的一个步骤，前面还有一个步骤是夺取附近一个碳原子的氢，主要是振荡的下面的点就是我要的ts，下面点的的振动模式也是正确的。

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