计算化学公社

标题: 过渡态优化SCF不收敛应如何调整？ [打印本页]

作者
Author: cheche710 时间: 2019-6-19 22:43
标题: 过渡态优化SCF不收敛应如何调整？
优化一个抽氢反应的过渡态，初猜将C-H键长拉长了30%左右，羟基自由基放在H原子附近，优化任务SCF不收敛关键词为：#p opt=(calcfc,ts,noeigen,cartesian) B3LYP 6-31g(d) scrf=smd empiricaldispersion=gd3bj
反应式：
(, 下载次数 Times of downloads: 44)

初猜结构：
(, 下载次数 Times of downloads: 42)
最后一帧结构：
(, 下载次数 Times of downloads: 52)
能量变化：
(, 下载次数 Times of downloads: 36)
请各位老师帮忙看看，感激不尽~

作者
Author: 小赵同学 时间: 2019-6-19 23:32
我觉得羟基应该摆在甲基旁边而不是氨基旁边

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

作者
Author: cheche710 时间: 2019-6-20 11:46

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

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

作者
Author: cheche710 时间: 2019-6-20 11:51

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

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

作者
Author: sobereva 时间: 2019-6-21 02:10

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（甚至直接砍掉）来节约时间，找到合适的过渡态后再用原本想用的基组进一步优化过渡态

作者
Author: cheche710 时间: 2019-7-1 22:04

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

谢谢老师的指点，在尝试之后还有2个问题想请教老师：
1.根据柔性扫描的结果，选择了第16个结构作为初猜，它的能量比初始结构低但是像过渡态的样子，是否合理？
(, 下载次数 Times of downloads: 29) (, 下载次数 Times of downloads: 35)
2.用初猜进行过渡态搜索，使用b3lyp/6-31g(d)优化得到的结构与反应物结构相似，使用m062x/def2SVP优化得到的结构与过渡态相似，是否应该选用m062x/def2SVP？
b3lyp/6-31g(d)结果：
(, 下载次数 Times of downloads: 31) (, 下载次数 Times of downloads: 32)
m062x/def2SVP结果:
(, 下载次数 Times of downloads: 29) (, 下载次数 Times of downloads: 36)

作者
Author: sobereva 时间: 2019-7-2 11:11

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

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

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

作者
Author: cheche710 时间: 2019-7-2 12:30
好的，非常感谢老师~

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