计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8304|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 环加成反应过渡态求助

[复制链接 Copy URL]

26

帖子

0

威望

193

eV
积分
219

Level 3 能力者

本帖最后由 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键,结果也是失败。

途径三:接着途径二最后一步结构,在固定C4-C35键的同时,柔性扫描C3-C27键(NO2-scan2.log),发现在2.27A处有最高点,如下图。对该点固定键长的单点做频率分析,振动模式正确。放开优化找过渡态,失败(NO2-scan2-TS.log),结果仍旧是两个底物分子跑开了。
途径四:基于前三种失败的反馈结果,猜想反应是协同过程,且过渡态键长较含-CN的过渡态键长短一些,于是尝试在更短的固定键长的优化结构基础上(固定优化的单点振动频率正确),放开优化,失败(NO2-fix-TS.log),结果仍旧是两个底物分子跑开了。

求助各位老师,先谢谢了!



Find-TS.rar

7.77 MB, 下载次数 Times of downloads: 181

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

2#
发表于 Post on 2021-8-5 19:41:23 | 只看该作者 Only view this author
在途径三的最高点附近,用更小的步长做扫描,然后优化过渡态,优化时步长设小一些,且算Hessian算得频繁一些
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

193

eV
积分
219

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-8-6 09:36:06 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-5 19:41
在途径三的最高点附近,用更小的步长做扫描,然后优化过渡态,优化时步长设小一些,且算Hessian算得频繁一 ...

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

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

4#
发表于 Post on 2021-8-6 14:34:40 | 只看该作者 Only view this author
ananbaby2309 发表于 2021-8-6 02:36
非常感谢回复!已经试了,还是不行,两个底物又跑开了。

是扫描都没有成功,还是扫描成功了但是优化没成功?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

193

eV
积分
219

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-8-6 15:50:33 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-6 14:34
是扫描都没有成功,还是扫描成功了但是优化没成功?

扫描成功了,优化不成功

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

6#
发表于 Post on 2021-8-6 16:36:00 | 只看该作者 Only view this author
ananbaby2309 发表于 2021-8-6 08:50
扫描成功了,优化不成功

如果这次扫描也有最高点,那过渡态应该是存在的,可以尝试进一步减小步长,并且每一步结构优化都算Hessian
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

655

帖子

1

威望

5314

eV
积分
5989

Level 6 (一方通行)

7#
发表于 Post on 2021-8-7 19:32:44 | 只看该作者 Only view this author
怀疑没戏啊,我拿 M062X/6-31G* 扫了一个简化的模型体系,二维势能面扫描,扫到很远了也没看到有到过渡态的趋势,再远我感觉没意义了

Find-TS.zip

1.44 MB, 下载次数 Times of downloads: 135

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

26

帖子

0

威望

193

eV
积分
219

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-8-13 20:09:39 | 只看该作者 Only view this author
北大-陶豫 发表于 2021-8-7 19:32
怀疑没戏啊,我拿 M062X/6-31G* 扫了一个简化的模型体系,二维势能面扫描,扫到很远了也没看到有到过渡态的 ...

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

26

帖子

0

威望

193

eV
积分
219

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-8-13 20:10:10 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-6 16:36
如果这次扫描也有最高点,那过渡态应该是存在的,可以尝试进一步减小步长,并且每一步结构优化都算Hessia ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 13:26 , Processed in 0.172580 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list