计算化学公社

标题: Gaussian计算氧酞菁铁催化丙烯环氧化的反应过渡态不收敛问题求助 [打印本页]

作者
Author: delltF 时间: 2021-7-20 17:42
标题: Gaussian计算氧酞菁铁催化丙烯环氧化的反应过渡态不收敛问题求助
各位老师好：
我在计算氧酞菁铁催化丙烯环氧化的反应，在以B3LYP/6-31G*/Lanl2dz级别寻找三重态下第二步过渡态（即碳-碳-氧闭环过程）时，计算文件如下，碰到了问题。
遇到的问题：(1)发现过渡态TS2已经优化到一个能量相对恒定的结构，但始终不收敛，且能量, 受力和位移都反复小幅波动，(2)之后我尝试了诸如缩小步长，改变优化方法（gdiis），尝试calcall，微调初始结构，降低收敛标准等方法发现均没有起效。(3)为了验证，我还选择优化途中的一个形式接近但未收敛结构单独计算了频率发现存在两个虚频（-497与-23），且较大虚频的振动方向正是目标虚频，
求助：鉴于这个情况，我想不出什么方法来帮助其收敛了，因此寻求各位的帮助，希望可以给出一些帮助收敛的方法。

输入文件：
(, 下载次数 Times of downloads: 16)
计算的反应和初始结构：
(, 下载次数 Times of downloads: 36) (, 下载次数 Times of downloads: 26)
优化过程中的能量，受力，位移变化：
(, 下载次数 Times of downloads: 35) (, 下载次数 Times of downloads: 32) (, 下载次数 Times of downloads: 40)

作者
Author: sobereva 时间: 2021-7-21 02:50
初猜质量太差，O和C之间距离也太远了。做个C-O距离的柔性扫描帮助判断合适的过渡态初猜。如果实在没经验，把闭环前后的结构都优化出来，用QST2找过渡态

作者
Author: delltF 时间: 2021-7-21 10:07

sobereva 发表于 2021-7-21 02:50
初猜质量太差，O和C之间距离也太远了。做个C-O距离的柔性扫描帮助判断合适的过渡态初猜。如果实在没经验， ...

多谢老师解答，我尝试一下。

作者
Author: 喵星大佬 时间: 2021-7-21 12:48
这体系用B3LYP不太好吧。。。。Fe上到sdd，可以先做个柔性扫描看看

作者
Author: delltF 时间: 2021-7-21 15:01

喵星大佬发表于 2021-7-21 12:48
这体系用B3LYP不太好吧。。。。Fe上到sdd，可以先做个柔性扫描看看

前置计算的级别都用了这个，为了保持一致也不是很好改，不过确实不太合适。
之前的初猜确实有点差，我现在先尝试着用QST2找过渡态，之后还不行的话就做柔性扫描试试。

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