计算化学公社

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

[新手求助] 关于求反应速率常数的疑问

[复制链接 Copy URL]

113

帖子

0

威望

1822

eV
积分
1935

Level 5 (御坂)

本帖最后由 409407227 于 2019-6-12 19:41 编辑

读到一篇求无过渡态反应的速率常数的计算,文章中这么来求反应速率常数:如下所述:For barrierless reactions, a variational canonical transition state theory (VTST) approach was adopted to compute the rate constants. VTST was employed in the following manner. After the ‘relaxed-scan’ calculations on bond lengths (R), a vibrational frequency (υ) calculation was performed to obtain the frequencies of the normal vibrational modes perpendicular to the reaction path. At each R, ‘trial’ canonical TST rate constants (kCaTST) were computed according to formula, based on υ and potential energy obtained from quantum chemistry calculations.其公式和运算过程如下图所述。

卢老师您好,我的疑问有两点:1是这篇文章求垂直于反应路径的振动频率的作用是什么,是指的虚频吗?  
                                              2就是他是通过计算柔性扫描所有点的自由能,然后找出了自由能垒吗?如果自由能也是单调升高的话,是不是就不能用KiSHELP求反应速率常数了?

TIM图片20190612193011.png (14.51 KB, 下载次数 Times of downloads: 37)

TIM图片20190612193011.png

TIM图片20190612193034.png (35.95 KB, 下载次数 Times of downloads: 29)

TIM图片20190612193034.png

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

2#
发表于 Post on 2019-6-12 20:04:18 | 只看该作者 Only view this author
1 即计算振动频率不考虑虚频模式(方向顺着IRC)。本身Gaussian输出的自由能也是忽略了虚频的贡献
2 是。对
北京科音自然科学研究中心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!

113

帖子

0

威望

1822

eV
积分
1935

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2019-6-12 20:24:55 | 只看该作者 Only view this author
sobereva 发表于 2019-6-12 20:04
1 即计算振动频率不考虑虚频模式(方向顺着IRC)。本身Gaussian输出的自由能也是忽略了虚频的贡献
2 是。 ...

谢谢卢老师的解答

103

帖子

0

威望

780

eV
积分
883

Level 4 (黑子)

4#
发表于 Post on 2019-7-19 15:14:31 | 只看该作者 Only view this author
sobereva 发表于 2019-6-12 20:04
1 即计算振动频率不考虑虚频模式(方向顺着IRC)。本身Gaussian输出的自由能也是忽略了虚频的贡献
2 是。 ...

老师:1、那如果是共价键直接断裂的,柔性扫描后发现没有自由能最大点,能量是一直增大的,请问这种情况怎么计算反应速率呢?不能用KiShelp计算速率,请问怎么办呢?
2、我之前发过帖子问过简单分子C2H4断键的反应,对称破缺态下柔性扫描了C2H4的C-C键断裂的过程,发现当C-C原子之间的距离到一定程度后,自由能变化上下起伏,一直搞不懂怎么回事。之前您说过如果每个点都收敛到了真正的基态波函数,能量应该也是连续的,但是我算的一直不连续,很是困扰。

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

5#
发表于 Post on 2019-7-19 23:46:33 | 只看该作者 Only view this author
长颈鹿先森 发表于 2019-7-19 15:14
老师:1、那如果是共价键直接断裂的,柔性扫描后发现没有自由能最大点,能量是一直增大的,请问这种情况 ...

1 扫描到能量平坦之后,把那个位置的自由能减去平衡状态的自由能作为自由能垒
2 轨迹、电子能量、振动频率和热力学校正量分别去考察搞清楚原因
北京科音自然科学研究中心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 02:35 , Processed in 0.225077 second(s), 23 queries , Gzip On.

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