计算化学公社

标题: 过渡态优化能量震荡求助 [打印本页]

作者
Author: 星月无光 时间: 2021-2-24 23:27
标题: 过渡态优化能量震荡求助
本帖最后由星月无光于 2021-2-25 12:04 编辑

大家好
在做过渡态优化的时候，得到的结构是我想要的过渡态结构，有且只有一个虚频（-300左右），原子移动方向也是朝着反应物与生成物方向。
但是优化曲线发生能量震荡，并且不收敛，四个优化参数都是NO

      Item             Value    Threshold  Converged?
Maximum Force             0.014024    0.000450    NO
RMS    Force                0.002260    0.000300    NO
Maximum Displacement    0.169431    0.001800    NO
RMS    Displacement       0.063329    0.001200    NO

G16版本，过渡态参数设置：# opt=(calcall,ts,noeigen) b3lyp/6-31g(d)  scf=qc
过渡态优化结果在图片上

scf=qc，scf=xqc，scf=NoDIIS（报错），这些都试过了，还是发生震荡，请教用什么方法可以解决啊？

作者
Author: hebrewsnabla 时间: 2021-2-24 23:36
本帖最后由 hebrewsnabla 于 2021-2-24 23:38 编辑

你这是结构优化不收敛，用了一堆scf的关键词当然没用。要对症下药，看这个http://sobereva.com/164

作者
Author: 北大-陶豫 时间: 2021-2-25 10:00
这是啥反应？α酮酸的羧基氢转移到α羰基上吗？或许你应该先扫描一下看看反应是不是真的有过渡态。溶剂化等有没有设对。
以后问这种问题应该把输入文件输出文件贴上来。
另外给自己的脚本打个广告：
【分享】监督高斯结构优化进度的脚本
http://bbs.keinsci.com/thread-21118-1-1.html

作者
Author: 星月无光 时间: 2021-2-25 13:27
本帖最后由星月无光于 2021-2-25 13:28 编辑

是一个羟基氢转移到羰基氢上，分子裂解为二氧化碳分子和相应的小分子
优化的方法我试过
1 对于DFT，特别是明尼苏达系列泛函（如M062X、M06、M06L），首先加上int=ultrafine再试（从G16开始此为默认，不再需要尝试此项）
2 如果体系小，尝试opt=calcall。如果体系大算得慢，此关键词会令计算量进一步激增，实在没辙再考虑这个，或者尝试opt=recalc=3或5
3 尝试opt=gdiis/gediis/rfo/caresian
4 尝试opt(gdiis,maxstep=x,notrust)，x为3~5
这些都试过了，还是发生震荡，最后一步得到的结果就是过渡态结构，并且虚频300多，但是就是能力震荡无法收敛

作者
Author: 北大-陶豫 时间: 2021-2-25 14:22

星月无光发表于 2021-2-25 13:27
是一个羟基氢转移到羰基氢上，分子裂解为二氧化碳分子和相应的小分子
优化的方法我试过
1 对于DFT，特别 ...

所以你是平白无故多出来两个氢？这怎么可能……

作者
Author: 星月无光 时间: 2021-2-25 14:58

北大-陶豫发表于 2021-2-25 14:22
所以你是平白无故多出来两个氢？这怎么可能……

哪里多出两个氢啊

作者
Author: 北大-陶豫 时间: 2021-2-25 15:46

星月无光发表于 2021-2-25 13:27
是一个羟基氢转移到羰基氢上，分子裂解为二氧化碳分子和相应的小分子
优化的方法我试过
1 对于DFT，特别 ...

你自己数数反应前后的氢数。。。

作者
Author: 北大-陶豫 时间: 2021-2-25 15:55

星月无光发表于 2021-2-25 14:58
哪里多出两个氢啊

(, 下载次数 Times of downloads: 68)
这两个氢哪来的？
如果把这两个氢去掉，你的产物就是个不稳定的卡宾，不太可能稳定存在，所以找不到过渡态一点也不奇怪。

作者
Author: 北大-陶豫 时间: 2021-2-25 16:09

北大-陶豫发表于 2021-2-25 15:55
这两个氢哪来的？
如果把这两个氢去掉，你的产物就是个不稳定的卡宾，不太可能稳定存在，所以找不到过 ...

当然也可能是脱羧之后羧基α位的羰基变成醛基，这样原子数也是合理的。但这个反应不太可能是直接分子内氢转移和碳碳键断裂一步完成的，更可能是一种能极性翻转的物质催化完成的（可以回顾一下生化里维生素B1催化丙酮酸脱羧生成乙醛的内容），至少也是负离子脱羧机理，绝对不可能像你想的一步完成那么简单。

作者
Author: 星月无光 时间: 2021-2-25 17:36

北大-陶豫发表于 2021-2-25 15:55
这两个氢哪来的？
如果把这两个氢去掉，你的产物就是个不稳定的卡宾，不太可能稳定存在，所以找不到过 ...

一语惊醒梦中人啊，太感谢了，我居然没有发现！

作者
Author: 北大-陶豫 时间: 2021-2-25 19:22

星月无光发表于 2021-2-25 17:36
一语惊醒梦中人啊，太感谢了，我居然没有发现！

想想也很简单，从 -CO-COOH 的羧基拿走一个二氧化碳，剩下的显然是醛基，而不是羟甲基。

作者
Author: DoubeeTwT 时间: 2021-2-26 18:48
能不能把输入的gjf文件的头传上来，有可能是电荷和自旋多重度设定错了

作者
Author: 北大-陶豫 时间: 2021-2-27 09:47

DoubeeTwT 发表于 2021-2-26 18:48
能不能把输入的gjf文件的头传上来，有可能是电荷和自旋多重度设定错了

他模型就建错了，给一个不存在的反应找过渡态……

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