计算化学公社

 找回密码 Forget password
 注册 Register

能量振荡无法收敛,已经很长时间了,实在没办法解决,各种办法都试过了,求助

查看数: 11304 | 评论数: 10 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-5-7 18:15

正文摘要:

本帖最后由 malong 于 2019-5-7 23:02 编辑 在寻找过渡态的时候,能量在两个结构之间振荡,能量最低的那个结构就是我要找的ts,但是由于振荡的原因无法收敛。通过在opt里面设置最大步长maxstep=2,能量依然振荡。 ...

回复 Reply

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

我这个是三甲胺加氢的一个步骤,前面还有一个步骤是夺取附近一个碳原子的氢,主要是振荡的下面的点就是我要的ts,下面点的的振动模式也是正确的。
sobereva 发表于 Post on 2019-5-8 09:30:49
malong 发表于 2019-5-7 20:03
这个是我把qst3(帖子中的结构)振荡的最后一个结构提取出来然后用ts的方法继续寻找过渡态,第36个点是我设 ...

我看不出你这是个什么反应。怎么碳上一个氢还跑到另一个氢上去了
当你要求优化TS,然而当前结构和实际TS偏离极大,都没正确TS的模样,反复震荡不收敛是正常的(哪怕收敛了,结果也没有意义)
PS:opt=TS的时候默认就用notrust,此时没必要写notrust
sobereva 发表于 Post on 2019-5-8 09:28:35
Frank 发表于 2019-5-8 08:15
借楼请教一下老师,如果用opt=(calcall,ts,noeigen)优化过渡态震荡不收敛,最后一帧结构接近过渡态结构, ...

有必要
如果说calcall能解决80%概率的不收敛问题,加上maxstep可能能提升到90%~95%概率
Frank 发表于 Post on 2019-5-8 08:15:50
sobereva 发表于 2019-5-7 18:58
看轨迹变化,如果末端结构看着像TS,加上maxstep=3再试
如果末端结构都没TS的样子,重新调整初猜

借楼请教一下老师,如果用opt=(calcall,ts,noeigen)优化过渡态震荡不收敛,最后一帧结构接近过渡态结构,振动模式也对,还有必要尝试maxstep=3,notrust吗?还是说calcall优化不出来的用maxstep=3也不会优化出来?
malong 发表于 Post on 2019-5-7 22:33:48
本帖最后由 malong 于 2019-5-7 22:49 编辑
alonewolfyang 发表于 2019-5-7 21:53
你在搞笑吗,这坨图片谁认识

我把帖子重新整理了,第一次发帖没有了解清楚规则,望见谅
malong 发表于 Post on 2019-5-7 22:28:47
alonewolfyang 发表于 2019-5-7 21:53
你在搞笑吗,这坨图片谁认识

老师第一次发帖没有考虑全面望见谅,我把帖子重新整理了一遍,您看现在可以不
alonewolfyang 发表于 Post on 2019-5-7 21:53:59
malong 发表于 2019-5-7 20:03
老师
这个是我把qst3(帖子中的结构)振荡的最后一个结构提取出来然后用ts的方法继续寻找过渡态,第36个 ...

你在搞笑吗,这坨图片谁认识
malong 发表于 Post on 2019-5-7 19:05:04
末端结构就是TS,我之前用maxstep=5还是振荡,我已经在用maxstep=2试了,图片描述的问题下次一定会注意的,非常感谢老师。
sobereva 发表于 Post on 2019-5-7 18:58:48
看轨迹变化,如果末端结构看着像TS,加上maxstep=3再试
如果末端结构都没TS的样子,重新调整初猜

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

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 00:33 , Processed in 0.178364 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list