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

Views: 4930|回复 Reply: 17

[新手求助]

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

[复制链接 Copy URL]

kenhe

kenhe

当前离线

UID39372

在线时间6 小时

注册时间2022-6-28

最后登录2023-11-23

好友0

12

帖子0

威望81

eV

Level 2 能力者, 积分 credits 93, 距离下一级还需 57 积分 credits

积分93

Level 2 能力者

发消息 Send PM

电梯直达 Go to

楼主

发表于 Post on 2022-8-5 22:41:19

|

只看该作者 Only view this author

|只看大图 View big picture

|倒序浏览 Reverse view

|阅读模式 Reading model

各位老师好，我有几个用sob老师编辑的TST Calculator计算反应速率常数的问题：

1.图1中红色剪头标记的V是什么意思？

2.图1中虚频输入的应该是绝对值吗，虚频是不是负数吗？

3.图1中Wigner和Skodje-Truhlar方法计算出来的透射系数κ不一样时，应该以哪个为准呢？

请各位老师，高手指教，谢谢。

202208052237356142..png

(38.99 KB, 下载次数 Times of downloads: 33)

下载附件 Download

2022-8-5 22:37 上传 Uploaded

图1

收藏 Add to favorites4

回复 Reply

举报 Report

sobereva

sobereva

当前在线

UID1

在线时间6122 小时

注册时间2014-10-7

最后登录2025-8-18

好友67

6万

帖子99

威望5万

eV

积分120205

管理员

公社社长

发消息 Send PM

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

回复 Reply

举报 Report

sobereva

sobereva

当前在线

UID1

在线时间6122 小时

注册时间2014-10-7

最后登录2025-8-18

好友67

6万

帖子99

威望5万

eV

积分120205

管理员

公社社长

发消息 Send PM

3#

发表于 Post on 2022-8-6 03:39:38

|

只看该作者 Only view this author

1 好好看单元格的注释，又不是没写得超级无比明白

1.png (16.39 KB, 下载次数 Times of downloads: 36)

下载附件 Download

2022-8-6 03:39 上传 Uploaded

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

3 好好看注释

2.png (13.38 KB, 下载次数 Times of downloads: 37)

下载附件 Download

2022-8-6 03:39 上传 Uploaded

评分 Rate

参与人数Participants 1eV +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

回复 Reply

举报 Report

wzkchem5

wzkchem5

当前离线

UID21811

在线时间1126 小时

注册时间2020-8-7

最后登录2025-8-17

好友75

1万

帖子0

威望9004

eV

积分20764

Level 6 (一方通行)

发消息 Send PM

4#

发表于 Post on 2022-8-6 04:34:04

|

只看该作者 Only view this author

虚频只是在量化软件里常常显示成负数而已，不代表它真的是负数。因为显示虚数要稍微多写个一两行的代码，显示成负数又不至于引起混淆，所以很多程序干脆就按负数显示了，仅此而已。

评分 Rate

参与人数Participants 1eV +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

本团队长期招收研究生，有意者可私信联系

回复 Reply

举报 Report

kenhe

kenhe

当前离线

UID39372

在线时间6 小时

注册时间2022-6-28

最后登录2023-11-23

好友0

12

帖子0

威望81

eV

Level 2 能力者, 积分 credits 93, 距离下一级还需 57 积分 credits

积分93

Level 2 能力者

发消息 Send PM

5#

楼主 Author|

发表于 Post on 2022-8-6 20:32:08

|

只看该作者 Only view this author

sobereva 发表于 2022-8-6 03:39

1 好好看单元格的注释，又不是没写得超级无比明白

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

回复 Reply

举报 Report

JCenter

JCenter

当前离线

UID29779

在线时间73 小时

注册时间2021-7-27

最后登录2025-8-15

好友7

236

帖子0

威望1229

eV

Level 4 (黑子), 积分 credits 1465, 距离下一级还需 35 积分 credits

积分1465

Level 4 (黑子)

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

发消息 Send PM

6#

发表于 Post on 2023-4-10 21:25:11

|

只看该作者 Only view this author

wzkchem5 发表于 2022-8-6 04:34

虚频只是在量化软件里常常显示成负数而已，不代表它真的是负数。因为显示虚数要稍微多写个一两行的代码，显 ...

老师，这个表格中计算透射系数所需要的虚频是指TS过渡态中存在的唯一一项虚频吗

目前专攻：

基于DFT、MD模拟的自由能计算

基于过渡态理论的自由能垒和反应速率常数计算

基于MD模拟的交联聚合物计算

回复 Reply

举报 Report

JCenter

JCenter

当前离线

UID29779

在线时间73 小时

注册时间2021-7-27

最后登录2025-8-15

好友7

236

帖子0

威望1229

eV

Level 4 (黑子), 积分 credits 1465, 距离下一级还需 35 积分 credits

积分1465

Level 4 (黑子)

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

发消息 Send PM

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模拟的交联聚合物计算

回复 Reply

举报 Report

wzkchem5

wzkchem5

当前离线

UID21811

在线时间1126 小时

注册时间2020-8-7

最后登录2025-8-17

好友75

1万

帖子0

威望9004

eV

积分20764

Level 6 (一方通行)

发消息 Send PM

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

本团队长期招收研究生，有意者可私信联系

回复 Reply

举报 Report

JCenter

JCenter

当前离线

UID29779

在线时间73 小时

注册时间2021-7-27

最后登录2025-8-15

好友7

236

帖子0

威望1229

eV

Level 4 (黑子), 积分 credits 1465, 距离下一级还需 35 积分 credits

积分1465

Level 4 (黑子)

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

发消息 Send PM

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模拟的交联聚合物计算

回复 Reply

举报 Report

wzkchem5

wzkchem5

当前离线

UID21811

在线时间1126 小时

注册时间2020-8-7

最后登录2025-8-17

好友75

1万

帖子0

威望9004

eV

积分20764

Level 6 (一方通行)

发消息 Send PM

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

本团队长期招收研究生，有意者可私信联系

回复 Reply

举报 Report

JCenter

JCenter

当前离线

UID29779

在线时间73 小时

注册时间2021-7-27

最后登录2025-8-15

好友7

236

帖子0

威望1229

eV

Level 4 (黑子), 积分 credits 1465, 距离下一级还需 35 积分 credits

积分1465

Level 4 (黑子)

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

发消息 Send PM

11#

发表于 Post on 2023-4-10 23:45:12

|

只看该作者 Only view this author

wzkchem5 发表于 2023-4-10 23:32

透射系数在结构优化所用理论级别下计算得到。

如果只在单点能级别下重新做频率分析，会导致结构优化和频 ...

感谢老师，您说的很细致，也很有逻辑。作为新手的我，听得很明白。再次感谢老师的耐心、细心答复。学生受益匪浅。

目前专攻：

基于DFT、MD模拟的自由能计算

基于过渡态理论的自由能垒和反应速率常数计算

基于MD模拟的交联聚合物计算

回复 Reply

举报 Report

sobereva

sobereva

当前在线

UID1

在线时间6122 小时

注册时间2014-10-7

最后登录2025-8-18

好友67

6万

帖子99

威望5万

eV

积分120205

管理员

公社社长

发消息 Send PM

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

回复 Reply

举报 Report

JCenter

JCenter

当前离线

UID29779

在线时间73 小时

注册时间2021-7-27

最后登录2025-8-15

好友7

236

帖子0

威望1229

eV

Level 4 (黑子), 积分 credits 1465, 距离下一级还需 35 积分 credits

积分1465

Level 4 (黑子)

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

发消息 Send PM

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模拟的交联聚合物计算

回复 Reply

举报 Report

LittlePupil

LittlePupil

当前离线

UID20457

在线时间377 小时

注册时间2020-6-10

最后登录2025-8-17

好友1

103

帖子0

威望3837

eV

Level 5 (御坂), 积分 credits 3940, 距离下一级还需 60 积分 credits

积分3940

Level 5 (御坂)

发消息 Send PM

14#

发表于 Post on 2023-4-11 14:11:16

|

只看该作者 Only view this author

wzkchem5 发表于 2023-4-10 23:32

透射系数在结构优化所用理论级别下计算得到。

如果只在单点能级别下重新做频率分析，会导致结构优化和频 ...

请问王老师，用Eckart方法计算透射系数时，如果代入的是高精度单点能“校正”过的能垒，得到的结果是否有意义？

「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

回复 Reply

举报 Report

wzkchem5

wzkchem5

当前离线

UID21811

在线时间1126 小时

注册时间2020-8-7

最后登录2025-8-17

好友75

1万

帖子0

威望9004

eV

积分20764

Level 6 (一方通行)

发消息 Send PM

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

本团队长期招收研究生，有意者可私信联系

回复 Reply

举报 Report