计算化学公社

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

[新手求助] 寻找转动过渡态问题

[复制链接 Copy URL]

139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

结构在附件里:我想找1结构到2结构的右边那个面旋转的过渡态,qst2方法找不到,直接将平面旋转一半之后再做TS也找不到,试了柔性扫描但是找完的过渡态振动频率只有-30.61,而且还掺杂了其他的振动,刚性扫描也试了,还是找不到过渡态,老师们有什么好的方法吗

求助.doc

142 KB, 下载次数 Times of downloads: 46

21

帖子

0

威望

1313

eV
积分
1334

Level 4 (黑子)

2#
发表于 Post on 2018-5-9 10:42:48 | 只看该作者 Only view this author
我之前用DFT找过类似的转动过渡态,情况跟你差不多。一位老师告诉我,DFT做转动不是很合适,用动力学更合适。

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

3#
发表于 Post on 2018-5-9 10:47:03 | 只看该作者 Only view this author
尽量直接贴图,不要放到文档里
刚性扫描对此问题没有任何用处
对你当前找的过渡态做IRC,能连通两侧就说明找对了
DFT做转动的过渡态非常适合,没有任何问题。当前问题牵扯不到动力学

注意恰当写标题,给你改了
http://bbs.keinsci.com/thread-9348-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

4#
发表于 Post on 2018-5-9 15:47:27 | 只看该作者 Only view this author
jgslinjin 发表于 2018-5-9 10:42
我之前用DFT找过类似的转动过渡态,情况跟你差不多。一位老师告诉我,DFT做转动不是很合适,用动力学更合适 ...

如果它说的转动与楼主例子中的转动是一回事儿,那不敢苟同这个意见

21

帖子

0

威望

1313

eV
积分
1334

Level 4 (黑子)

5#
发表于 Post on 2018-5-9 17:20:50 | 只看该作者 Only view this author
我之前做的是苯乙烯环氧化,端基碳连接氧后生成自由基,可能发生C-C键旋转,如下图。当初还scan过二面角,也出现了峰,后来就没往后做了。







139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2018-5-9 19:52:44 | 只看该作者 Only view this author
sobereva 发表于 2018-5-9 10:47
尽量直接贴图,不要放到文档里
刚性扫描对此问题没有任何用处
对你当前找的过渡态做IRC,能连通两侧就说 ...

老师做了IRC但是只IRC曲线上只有三个点,两端不是反应物和产物,IRC是不正确的

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

7#
发表于 Post on 2018-5-10 03:22:00 | 只看该作者 Only view this author
大慧慧 发表于 2018-5-9 19:52
老师做了IRC但是只IRC曲线上只有三个点,两端不是反应物和产物,IRC是不正确的

刚一跑就停的问题仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

8#
发表于 Post on 2018-5-11 22:24:16 | 只看该作者 Only view this author
本帖最后由 mutron 于 2018-5-11 22:30 编辑

题主碰到的估计是IRC每一边都只跑一个点,有可能是碰到势能面比较平缓的情况,可以尝试把stepsize设置大一点解决问题。用PM7扫了一下:

    stepsize设为15的话,只有三个点,两边都是跑一个点就停了

    可以看到两边的第一点受力已经非常小,估计程序自动判断为终点停了

    stepsize设置为20(或者更大点也可以)就比较正常了

139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2018-5-24 21:27:50 | 只看该作者 Only view this author
mutron 发表于 2018-5-11 22:24
题主碰到的估计是IRC每一边都只跑一个点,有可能是碰到势能面比较平缓的情况,可以尝试把stepsize设置大一 ...

您的意思是说我将步长设的大一点,还是我要换成PM7的方法重新跑一下呀

139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2018-5-24 21:29:13 | 只看该作者 Only view this author
本帖最后由 大慧慧 于 2018-5-24 21:30 编辑
mutron 发表于 2018-5-11 22:24
题主碰到的估计是IRC每一边都只跑一个点,有可能是碰到势能面比较平缓的情况,可以尝试把stepsize设置大一 ...

也就是说这个面旋转的过渡态是对的,是不。irc是可以跑出来的对吗?我现在就是不太确定我的这个过渡态找的对不对,它的振动频率只有-30点多

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

11#
发表于 Post on 2018-5-25 19:26:41 | 只看该作者 Only view this author
本帖最后由 mutron 于 2018-5-27 19:41 编辑
大慧慧 发表于 2018-5-24 21:27
您的意思是说我将步长设的大一点,还是我要换成PM7的方法重新跑一下呀

步长设大一点就行(逐渐加大到可以连续跑出点),PM7只是省时间看看问题出在哪里而已,最终还是需要你用原来的方法跑一遍。留意到你给的结构里面居然是双自由基,我8楼跑的是桥头碳旁边的碳原子不是自由基(更合理)、氧原子和氧原子在同一边的双线态,虚频为-53;B3LYP/6-31G*下是-57:
     



用PM7算一遍你图3中的三线态结构,虚频为-35:
           只是IRC更短了

139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2018-5-26 21:10:29 | 只看该作者 Only view this author
mutron 发表于 2018-5-25 19:26
步长设大一点就行(逐渐加大到可以连续跑出点),PM7只是省时间看看问题出在哪里而已,最终还是需要你用 ...

您的意思是我把右边的按个面反向旋转到一个位置在跑过渡态比较合理一点对吗

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

13#
发表于 Post on 2018-5-27 19:40:41 | 只看该作者 Only view this author
大慧慧 发表于 2018-5-26 21:10
您的意思是我把右边的按个面反向旋转到一个位置在跑过渡态比较合理一点对吗

我觉得两个都是,毕竟两个环不对称

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 03:02 , Processed in 0.172886 second(s), 23 queries , Gzip On.

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