计算化学公社

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

[新手求助] Gaussian限制性优化结果有振动方向正确的虚频,但TS后虚频方向不对

[复制链接 Copy URL]

14

帖子

0

威望

93

eV
积分
107

Level 2 能力者

各位老师好,目前在研究如下两个分子间H转移的过渡态,32H从33N转移到31O
因为过渡态比较难找,所以先做了限制性优化固定32、33;31、32的键长。(gjf见附件)
得到的优化最后一帧有一个-1464的虚频,对应振动方向为初猜的反应方向。于是保存这一帧、改关键词为opt=(ts,noeigen,calcfc)。结果中只有一个-31的虚频,且振动方向不是反应方向(如图)
请问各位老师,出现这种情况有可能是什么原因?如何解决以寻找过渡态?






in.gjf

3.5 KB, 下载次数 Times of downloads: 13

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

2#
发表于 Post on 2024-3-6 02:33:16 | 只看该作者 Only view this author
做32-31距离缩短的柔性扫描,如果有极大点,当做过渡态初猜
详谈使用Gaussian做势能面扫描
http://sobereva.com/474http://bbs.keinsci.com/thread-12660-1-1.html

为了节约柔性扫描时间,可以把烃链截短。转移的H应当带上极化函数
北京科音自然科学研究中心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!

174

帖子

0

威望

2865

eV
积分
3039

Level 5 (御坂)

3#
发表于 Post on 2024-3-6 13:33:33 | 只看该作者 Only view this author
你的做法总体上是没问题的,这种状况其实也挺常见的。解决方案按照sob的来,问题也基本能解决。但是我有一个不太明白的地方,你做这个反应不用隐式溶剂模型吗?我记得有些反应在溶剂下能找到过渡态,但是在真空条件下就是找不到过渡态的。

14

帖子

0

威望

93

eV
积分
107

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2024-3-7 13:17:15 | 只看该作者 Only view this author
GoldenBaby 发表于 2024-3-6 13:33
你的做法总体上是没问题的,这种状况其实也挺常见的。解决方案按照sob的来,问题也基本能解决。但是我有一 ...

谢谢老师,正在按照Sob老师的方法尝试计算。
在跑过渡态的时候都是带上了隐式溶剂模型的(文献中为二氯甲烷

14

帖子

0

威望

93

eV
积分
107

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2024-3-8 10:13:22 | 只看该作者 Only view this author
本帖最后由 jsczfx 于 2024-3-8 11:12 编辑
sobereva 发表于 2024-3-6 02:33
做32-31距离缩短的柔性扫描,如果有极大点,当做过渡态初猜
详谈使用Gaussian做势能面扫描
http://sobere ...

SOB老师好,尝试了柔性扫描,使用极大点作为过渡态初猜、加上隐式溶剂(二氯甲烷)、基组改为6-31**。但是跑opt=(calcfc,ts,noeigen),不收敛且最后震荡的几帧H跑回N上了(反应物),请问是应该尝试小步长看看吗?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

6#
发表于 Post on 2024-3-9 05:39:51 | 只看该作者 Only view this author
jsczfx 发表于 2024-3-8 10:13
SOB老师好,尝试了柔性扫描,使用极大点作为过渡态初猜、加上隐式溶剂(二氯甲烷)、基组改为6-31**。但 ...

好好写基组名

用recalc=3结合maxstep=10再试
北京科音自然科学研究中心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!

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 08:43 , Processed in 0.189264 second(s), 23 queries , Gzip On.

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