计算化学公社

标题: 环加成反应过渡态求助 [打印本页]

作者
Author: ananbaby2309 时间: 2021-8-5 18:11
标题: 环加成反应过渡态求助
本帖最后由 ananbaby2309 于 2021-8-5 18:15 编辑

各位老师大家好！我在找环加成反应的过渡态过程中，遇到了以下问题：当取代基为-CN时，反应能顺利找到过渡态（文件CN-TS.log），当取代基为NO2时，怎么都找不到过渡态。我采取了以下办法：
途径一：在CN-TS.log结构基础上，将-CN换为-NO2，固定（C4-C32, C3-C27键）优化，频率分析发现振动模式正确，然后放开找过渡态，结果两个底物跑开了，过渡态寻找失败（NO2-TS1.log)。也尝试过iop(1/8=5)，以及calcall的条件，用m062x/def2svp，结果均是两个底物分子跑开了。
途径二：尝试从产物结构出发，柔性扫描拉长C4-C35键（NO2-scan1.log），如下图，并对最后一个固定优化单点做频率分析，发现振动模式为正确模式，便放开优化找过渡态，失败(NO2-scan1-TS.log)。同时，也尝试过扫描拉长C3-C27键，结果也是失败。
(, 下载次数 Times of downloads: 23)
途径三：接着途径二最后一步结构，在固定C4-C35键的同时，柔性扫描C3-C27键（NO2-scan2.log），发现在2.27A处有最高点，如下图。对该点固定键长的单点做频率分析，振动模式正确。放开优化找过渡态，失败(NO2-scan2-TS.log)，结果仍旧是两个底物分子跑开了。 (, 下载次数 Times of downloads: 41)
途径四：基于前三种失败的反馈结果，猜想反应是协同过程，且过渡态键长较含-CN的过渡态键长短一些，于是尝试在更短的固定键长的优化结构基础上（固定优化的单点振动频率正确），放开优化，失败(NO2-fix-TS.log)，结果仍旧是两个底物分子跑开了。

求助各位老师，先谢谢了！

作者
Author: wzkchem5 时间: 2021-8-5 19:41
在途径三的最高点附近，用更小的步长做扫描，然后优化过渡态，优化时步长设小一些，且算Hessian算得频繁一些

作者
Author: ananbaby2309 时间: 2021-8-6 09:36

wzkchem5 发表于 2021-8-5 19:41
在途径三的最高点附近，用更小的步长做扫描，然后优化过渡态，优化时步长设小一些，且算Hessian算得频繁一 ...

非常感谢回复！已经试了，还是不行，两个底物又跑开了。

作者
Author: wzkchem5 时间: 2021-8-6 14:34

ananbaby2309 发表于 2021-8-6 02:36
非常感谢回复！已经试了，还是不行，两个底物又跑开了。

是扫描都没有成功，还是扫描成功了但是优化没成功？

作者
Author: ananbaby2309 时间: 2021-8-6 15:50

wzkchem5 发表于 2021-8-6 14:34
是扫描都没有成功，还是扫描成功了但是优化没成功？

扫描成功了，优化不成功

作者
Author: wzkchem5 时间: 2021-8-6 16:36

ananbaby2309 发表于 2021-8-6 08:50
扫描成功了，优化不成功

如果这次扫描也有最高点，那过渡态应该是存在的，可以尝试进一步减小步长，并且每一步结构优化都算Hessian

作者
Author: 北大-陶豫 时间: 2021-8-7 19:32
怀疑没戏啊，我拿 M062X/6-31G* 扫了一个简化的模型体系，二维势能面扫描，扫到很远了也没看到有到过渡态的趋势，再远我感觉没意义了

作者
Author: ananbaby2309 时间: 2021-8-13 20:09

北大-陶豫发表于 2021-8-7 19:32
怀疑没戏啊，我拿 M062X/6-31G* 扫了一个简化的模型体系，二维势能面扫描，扫到很远了也没看到有到过渡态的 ...

非常感谢！我这几天也想了其他办法，也都失败了

作者
Author: ananbaby2309 时间: 2021-8-13 20:10

wzkchem5 发表于 2021-8-6 16:36
如果这次扫描也有最高点，那过渡态应该是存在的，可以尝试进一步减小步长，并且每一步结构优化都算Hessia ...

Calcall 加了，也没用……。非常感谢建议和帮助！

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