计算化学公社

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

[新手求助] 求助利用TST原理计算反应速率常数的几个问题

[复制链接 Copy URL]

12

帖子

0

威望

81

eV
积分
93

Level 2 能力者

各位老师好,我有几个用sob老师编辑的TST Calculator计算反应速率常数的问题:
1.图1中红色剪头标记的V是什么意思?
2.图1中虚频输入的应该是绝对值吗,虚频是不是负数吗?
3.图1中Wigner和Skodje-Truhlar方法计算出来的透射系数κ不一样时,应该以哪个为准呢?
请各位老师,高手指教,谢谢。

202208052237356142..png (38.99 KB, 下载次数 Times of downloads: 29)

图1

图1

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

2#
发表于 Post on 2022-8-6 03:36:12 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“利用TST原理计算反应速率常数”改了,以后务必注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

3#
发表于 Post on 2022-8-6 03:39:38 | 只看该作者 Only view this author
1 好好看单元格的注释,又不是没写得超级无比明白


2 复习虚数和负数的定义。哪来的“虚频是不是负数吗”

3 好好看注释

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
kenhe + 3 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

8993

eV
积分
20741

Level 6 (一方通行)

4#
发表于 Post on 2022-8-6 04:34:04 | 只看该作者 Only view this author
虚频只是在量化软件里常常显示成负数而已,不代表它真的是负数。因为显示虚数要稍微多写个一两行的代码,显示成负数又不至于引起混淆,所以很多程序干脆就按负数显示了,仅此而已。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
kenhe + 2 谢谢

查看全部评分 View all ratings

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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

12

帖子

0

威望

81

eV
积分
93

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-8-6 20:32:08 | 只看该作者 Only view this author
sobereva 发表于 2022-8-6 03:39
1 好好看单元格的注释,又不是没写得超级无比明白

不好意思,Sob老师,我之前下载的计算表格里可能是因为我excel版本问题,没有显示注释,感谢您的回答。

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

6#
发表于 Post on 2023-4-10 21:25:11 | 只看该作者 Only view this author
wzkchem5 发表于 2022-8-6 04:34
虚频只是在量化软件里常常显示成负数而已,不代表它真的是负数。因为显示虚数要稍微多写个一两行的代码,显 ...

老师,这个表格中计算透射系数所需要的虚频是指TS过渡态中存在的唯一一项虚频吗
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

7#
发表于 Post on 2023-4-10 22:09:25 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2023-4-10 22:21 编辑
sobereva 发表于 2022-8-6 03:39
1 好好看单元格的注释,又不是没写得超级无比明白

sob老师,对于使用excel表格,TST理论计算k,我有几个关于透射系数相关的疑问想向您请教:
1.温度越低,隧道效应越显著,那针对常温以上的反应温度,即大于等于298.15K,还有必要计算透射系数吗,直接默认透射系数等于1是否可以
2.excel表格中透射系数的计算,所使用的虚频是指过渡态寻找中出现的唯一一个虚频数值吗3..excel表格中透射系数的计算所需要的温度与计算k的反应温度一致吗
4.使用完整的Skodje-Truhlar方法,正向势垒为过渡态减去反应物的内能,两者的内能都是0k下的吗
5.使用excel表格中完整的Skodje-Truhlar方法计算透射系数与KiSThelP使用Eckart方法相比,精度如何,使用excel表格中完整的Skodje-Truhlar方法计算透射系数校正k会被审稿人提问吗

目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

1万

帖子

0

威望

8993

eV
积分
20741

Level 6 (一方通行)

8#
发表于 Post on 2023-4-10 22:57:58 | 只看该作者 Only view this author
JCenter 发表于 2023-4-10 15:09
sob老师,对于使用excel表格,TST理论计算k,我有几个关于透射系数相关的疑问想向您请教:
1.温度越低, ...

1. 当且仅当已知你这个反应在室温下隧道效应已经可以忽略时可以。温度越低隧道效应越显著不假,但是注意这个“越”字。很多时候即便高于室温,隧穿仍然不可忽略,只不过室温时或低于室温时更不可忽略而已
2. 显然是的,既然虚频是唯一的,当然只能取那个唯一的虚频
3. 需要,否则计算出的做了透射系数校正的速率常数没有物理意义
4. 注释已经写了
5. 个人觉得Eckart方法用得多。能用Polyrate的SCT、LCT等方法最好,因为原理上最准,其次用Eckart,再其次用Skodje-Truhlar。如果我是审稿人,看到你用Skodje-Truhlar,会问你为什么不用Eckart,但是如果你能给出合理的解释,说明你的体系确实用不了Eckart,或者能说明对于你的体系Skodje-Truhlar比Eckart好,那我会放你过。假如充分了解隧穿计算的理论以后仍然找不出来这样的理由的话,那还是用Eckart或者更精确的方法比较好
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

9#
发表于 Post on 2023-4-10 23:09:21 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2023-4-10 23:10 编辑
wzkchem5 发表于 2023-4-10 22:57
1. 当且仅当已知你这个反应在室温下隧道效应已经可以忽略时可以。温度越低隧道效应越显著不假,但是注意 ...

好嘞,谢谢老师耐心指点。忘记补充最后一个疑问,利用excel表格计算k时,自由能垒的计算都是在高精度,比如双杂化泛函、高档3zeta基组下计算得到,而虚频的来源,振动分析是在一般泛函、一般2zeta基组下获得的,用虚频得到的透射系数,用于TST计算k的公式可以吗,虚频数需要重新将过渡态在单点能所使用的计算级别条件下进行振动分析,以获得相应级别下的虚频数呢,还是说虚频的来源和计算级别,与透射系数的数值精度、乃至k值有无关联呢
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

1万

帖子

0

威望

8993

eV
积分
20741

Level 6 (一方通行)

10#
发表于 Post on 2023-4-10 23:32:12 | 只看该作者 Only view this author
JCenter 发表于 2023-4-10 16:09
好嘞,谢谢老师耐心指点。忘记补充最后一个疑问,利用excel表格计算k时,自由能垒的计算都是在高精度,比 ...

透射系数在结构优化所用理论级别下计算得到。
如果只在单点能级别下重新做频率分析,会导致结构优化和频率分析理论级别不统一,结果无意义;如果在单点能级别下重新既做结构优化又做频率分析,那么结果可用,且比直接在原结构优化所用理论级别下计算要准,但那相当于你换了个更高的结构优化的级别,不改变“透射系数在结构优化所用理论级别下计算得到”这个规则本身。
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

11#
发表于 Post on 2023-4-10 23:45:12 | 只看该作者 Only view this author
wzkchem5 发表于 2023-4-10 23:32
透射系数在结构优化所用理论级别下计算得到。
如果只在单点能级别下重新做频率分析,会导致结构优化和频 ...

感谢老师,您说的很细致,也很有逻辑。作为新手的我,听得很明白。再次感谢老师的耐心、细心答复。学生受益匪浅。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

12#
发表于 Post on 2023-4-11 09:57:26 | 只看该作者 Only view this author
JCenter 发表于 2023-4-10 22:09
sob老师,对于使用excel表格,TST理论计算k,我有几个关于透射系数相关的疑问想向您请教:
1.温度越低, ...

DOI: 10.1002/anie.201708489里对SCT和Truhlar的公式算的透射系数做了对比,Truhlar的公式精度已经不错了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

13#
发表于 Post on 2023-4-11 10:18:48 | 只看该作者 Only view this author
sobereva 发表于 2023-4-11 09:57
DOI: 10.1002/anie.201708489里对SCT和Truhlar的公式算的透射系数做了对比,Truhlar的公式精度已经不错了

好嘞,谢谢老师!
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

103

帖子

0

威望

3829

eV
积分
3932

Level 5 (御坂)

14#
发表于 Post on 2023-4-11 14:11:16 | 只看该作者 Only view this author
wzkchem5 发表于 2023-4-10 23:32
透射系数在结构优化所用理论级别下计算得到。
如果只在单点能级别下重新做频率分析,会导致结构优化和频 ...

请问王老师,用Eckart方法计算透射系数时,如果代入的是高精度单点能“校正”过的能垒,得到的结果是否有意义?
「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

1万

帖子

0

威望

8993

eV
积分
20741

Level 6 (一方通行)

15#
发表于 Post on 2023-4-11 17:03:58 | 只看该作者 Only view this author
sobereva 发表于 2023-4-11 02:57
DOI: 10.1002/anie.201708489里对SCT和Truhlar的公式算的透射系数做了对比,Truhlar的公式精度已经不错了

这篇文章只做了heavy atom tunneling,不知道有没有氢的隧穿的benchmark文章?
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 18:19 , Processed in 0.163515 second(s), 24 queries , Gzip On.

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