过渡态优化SCF不收敛应如何调整？

cheche710 · 发表于 Post on 2019-6-19 22:43:48

优化一个抽氢反应的过渡态，初猜将C-H键长拉长了30%左右，羟基自由基放在H原子附近，优化任务SCF不收敛关键词为：#p opt=(calcfc,ts,noeigen,cartesian) B3LYP 6-31g(d) scrf=smd empiricaldispersion=gd3bj
反应式：

初猜结构：

最后一帧结构：

能量变化：

请各位老师帮忙看看，感激不尽~

小赵同学 · 发表于 Post on 2019-6-19 23:32:42

我觉得羟基应该摆在甲基旁边而不是氨基旁边

sobereva · 发表于 Post on 2019-6-20 03:55:07

问题不够确切，你说是SCF不收敛，但你的描述并没充分体现出是遇到了SCF不收敛，更像是几何优化达到步数上限而报错结束
你的关键词没什么好改的，反复尝试不同初猜结构。当前最终结构都没过渡态的苗头，继续跑下去也白搭

cheche710 · 发表于 Post on 2019-6-20 11:46:30

小赵同学发表于 2019-6-19 23:32
我觉得羟基应该摆在甲基旁边而不是氨基旁边

请问为什么应该摆在甲基旁边呢？（想考察的是在氨基部分发生的抽氢反应，所以放在了氨基旁边）

cheche710 · 发表于 Post on 2019-6-20 11:51:40

sobereva 发表于 2019-6-20 03:55
问题不够确切，你说是SCF不收敛，但你的描述并没充分体现出是遇到了SCF不收敛，更像是几何优化达到步数上限 ...

抱歉老师，确实描述的不充分。请问老师初猜结构应该如何改变呢，初学过渡态搜索没有经验，麻烦老师指点一下，谢谢老师

sobereva · 发表于 Post on 2019-6-21 02:10:08

cheche710 发表于 2019-6-20 11:51
抱歉老师，确实描述的不充分。请问老师初猜结构应该如何改变呢，初学过渡态搜索没有经验，麻烦老师指点一 ...

反复调节OH的位置、N-H的键长进行尝试
实在搞不出来，可以借助柔性扫描H...O距离帮助确定初猜，看
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

另外，为了省时间，试过渡态的阶段可以把苯环右边的那些结构都用很小的基组比如3-21G（甚至直接砍掉）来节约时间，找到合适的过渡态后再用原本想用的基组进一步优化过渡态

cheche710 · 发表于 Post on 2019-7-1 22:04:13

sobereva 发表于 2019-6-21 02:10
反复调节OH的位置、N-H的键长进行尝试
实在搞不出来，可以借助柔性扫描H...O距离帮助确定初猜，看
详谈 ...

谢谢老师的指点，在尝试之后还有2个问题想请教老师：
1.根据柔性扫描的结果，选择了第16个结构作为初猜，它的能量比初始结构低但是像过渡态的样子，是否合理？

2.用初猜进行过渡态搜索，使用b3lyp/6-31g(d)优化得到的结构与反应物结构相似，使用m062x/def2SVP优化得到的结构与过渡态相似，是否应该选用m062x/def2SVP？
b3lyp/6-31g(d)结果：

m062x/def2SVP结果:

sobereva · 发表于 Post on 2019-7-2 11:11:06

cheche710 发表于 2019-7-1 22:04
谢谢老师的指点，在尝试之后还有2个问题想请教老师：
1.根据柔性扫描的结果，选择了第16个结构作为初猜 ...

你既然用了m062x/def2SVP就没必要再用b3lyp/6-31g(d)，前者对这类问题肯定比后者强

对前者优化出的看似是过渡态的结构跑IRC，走通了过渡态就对了，没走成功就还是不对。

cheche710 · 发表于 Post on 2019-7-2 12:30:13

好的，非常感谢老师~

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

[新手求助] 过渡态优化SCF不收敛应如何调整？

浏览过的版块