过渡态优化能量震荡求助

星月无光 · 发表于 Post on 2021-2-24 23:27:52

本帖最后由星月无光于 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（报错），这些都试过了，还是发生震荡，请教用什么方法可以解决啊？

hebrewsnabla · 发表于 Post on 2021-2-24 23:36:49

本帖最后由 hebrewsnabla 于 2021-2-24 23:38 编辑

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

北大-陶豫 · 发表于 Post on 2021-2-25 10:00:54

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

星月无光 · 发表于 Post on 2021-2-25 13:27:17

本帖最后由星月无光于 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多，但是就是能力震荡无法收敛

北大-陶豫 · 发表于 Post on 2021-2-25 14:22:52

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

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

星月无光 · 发表于 Post on 2021-2-25 14:58:30

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

哪里多出两个氢啊

北大-陶豫 · 发表于 Post on 2021-2-25 15:46:50

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

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

北大-陶豫 · 发表于 Post on 2021-2-25 15:55:02

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

这两个氢哪来的？
如果把这两个氢去掉，你的产物就是个不稳定的卡宾，不太可能稳定存在，所以找不到过渡态一点也不奇怪。

北大-陶豫 · 发表于 Post on 2021-2-25 16:09:33

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

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

星月无光 · 发表于 Post on 2021-2-25 17:36:52

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

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

北大-陶豫 · 发表于 Post on 2021-2-25 19:22:43

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

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

DoubeeTwT · 发表于 Post on 2021-2-26 18:48:13

能不能把输入的gjf文件的头传上来，有可能是电荷和自旋多重度设定错了

北大-陶豫 · 发表于 Post on 2021-2-27 09:47:35

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

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

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

[Gaussian/gview] 过渡态优化能量震荡求助

浏览过的版块