计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: 肖王木木。
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 求助无势垒反应计算速率常数的问题

[复制链接 Copy URL]

255

帖子

0

威望

1422

eV
积分
1677

Level 5 (御坂)

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

16#
发表于 Post on 2024-11-25 21:17:53 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2024-11-25 21:19 编辑
wzkchem5 发表于 2024-11-1 19:17
1. 可以证明任何双分子变成单分子的反应A+B -> C都有自由能垒。这是因为当A和B之间的距离趋于无穷远时, ...

老师,我看了变分过渡态的那篇综述。不知道我理解的是否正确:相比于经典的过渡态理论,广义的过渡态理论是在不考虑反应坐标这一自由度下,不再用势能面的鞍点作为过渡态,在某一温度下,从反应物到产物的分割面是在变化的。通过将其分割面做变分处理得到广义过渡态的方法即变分过渡态。其中正则变分过渡态方法CVT是做某一温度下,得到使k最小的过渡态-沿反应坐标做自由能计算。感觉Baik基于PMF的方法,很像cVT的思想。但是不知道Baik法和cVT的区别在哪,总觉得Baik法并不等于cVT。但是Baik法肯定属于广义过渡态的范畴。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

17#
发表于 Post on 2024-11-26 12:41:17 | 只看该作者 Only view this author
JCenter 发表于 2024-11-25 14:17
老师,我看了变分过渡态的那篇综述。不知道我理解的是否正确:相比于经典的过渡态理论,广义的过渡态理论 ...


一般CVT在IRC上面找变分过渡态,但是要想有IRC,能量势能面上必须有过渡态,才能跑IRC。Baik的方法用柔性扫描代替IRC,避开了“能量势能面上必须有过渡态”的要求
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

255

帖子

0

威望

1422

eV
积分
1677

Level 5 (御坂)

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

18#
发表于 Post on 2024-11-26 19:05:55 | 只看该作者 Only view this author
wzkchem5 发表于 2024-11-26 12:41

一般CVT在IRC上面找变分过渡态,但是要想有IRC,能量势能面上必须有过渡态,才能跑IRC。Baik的方法用 ...

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

42

帖子

0

威望

485

eV
积分
527

Level 4 (黑子)

19#
发表于 Post on 2026-1-26 16:15:54 | 只看该作者 Only view this author
wzkchem5 发表于 2020-11-15 14:57
这种反应要把自由能垒算准,除了用非常耗时的VRC-VTST以外,目前还没有公认的满意的解决方案。最常见的方法 ...

老师您好,我在计算一个催化过程,其中有一步中间体IM2转IM3整个过程ΔG=3kJ/mol,根据经验推测应该是一个势垒极低的反应。不知道是不是因为势能面太过平缓,用TS(Berry)算得不到过渡态结构,优化所得经常是没虚频的结构。
由于整个过程已知存在其他决速步,这一步能垒对算反应速率常数没啥影响,但整个机理缺失了这一环怕别人不认可。
请问如果不需要得到这一步准确能垒的话,可以做怎样的一个计算来让别人信服IM2转IM3的过程能垒很低,不影响整体动力学的推导。

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

20#
发表于 Post on 2026-1-27 09:20:32 | 只看该作者 Only view this author
ginlpein 发表于 2026-1-26 16:15
老师您好,我在计算一个催化过程,其中有一步中间体IM2转IM3整个过程ΔG=3kJ/mol,根据经验推测应该是一 ...

如果这一步是一个特别简单、well-studied的反应,比如质子转移、构象变化,那么往往可以不算。
如果你觉得审稿人问这一步的概率不可忽略,那么还是算一下,高斯收敛不了就用ORCA收敛(例如用NEB方法),再把结构放到高斯里再算即可
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-16 16:55 , Processed in 0.297485 second(s), 20 queries , Gzip On.

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