liyuanhe211 发表于 2023-2-8 18:04 谢谢老师解答！确实一般而言自由基裂解反应速率会远高于自由基生成的反应，在反应体系中并不是决速步。 |
alanzhou 发表于 2023-2-8 17:44 我的意思是，如果你研究的是某个实际问题，你现在作为底物的自由基物种在反应中一定是通过其他步骤生成的，这些步骤的速率（尤其是从闭壳层分子生成自由基物种的速率、或者自由基转移的速率）可能远低于你当前帖子中所研究的这个步骤，当前步骤不是决速的，因此，把这个步骤的速率搞清楚对于与实验反应速率比较可能没有多少贡献。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 谢谢 |
liyuanhe211 发表于 2023-2-8 11:39 老师您好！反应物就是IRC曲线最右端，单分子的醇自由基（从-OH基团起第二个C位点失去一个H原子形成的），然后基元反应是醇燃烧过程分子自由基的beta裂解反应，C-O键断裂，生成产物烯烃与OH自由基。通过旋转不同的二面角，而后优化获得所有的分子构型，选择所有构型中能量最低的用于速率常数计算。其中优化得到的反应物的能量明显低于过渡态，产物能量却高于过渡态。现在经过分析这个反应可能不存在理论意义上的过渡态，无法用经典TST理论计算速率常数。请问老师您看是不是需要用其他的方法计算？是VTST吗？ |
alanzhou 发表于 2023-2-6 08:32 那你文件中的反应物是怎么生成的？ |
suntree 发表于 2023-2-5 17:34 Hi @suntree I was using Chemcraft to generate these two images. Just search for it by using the keyword "Chemcraft". I am afraid that I am not allowed to give a link to a commercial program here. Note that the CPK style was used in this case. Cheers |
chands 发表于 2023-2-6 00:14 谢谢老师解答！ |
本帖最后由 alanzhou 于 2023-2-6 08:35 编辑 liyuanhe211 发表于 2023-2-6 00:05 谢谢老师解答！现在我的课题是计算这个基元反应在气相条件下的反应速率，因为产物能量就高于过渡态，原本是想先计算得到一个能量低于过渡态的中间体，然后用TST理论计算反应速率，现在看来好像无法解决。能请问老师，在产物能量高于过渡态的情况下，用什么方法可以得到自由基裂解反应的速率常数吗？ |
dnlx 发表于 2023-2-5 23:48 谢谢老师解答！ |
alanzhou 发表于 2023-2-5 23:11 很多自由基反应都是无势垒反应。本来你的醇类自由基就不稳定，生成的烯烃和OH自由基相对稳定，这个反应肯定是thermodynamically favoured。IRC是基于单点能（电子能）的。如dnlx所说，如果IRC跑出来的（电子）能垒很小，算上热力学校正量，可能TS的自由能反而比反应物（产物）的自由能低，就没有过渡态了。如liyuanhe所说，就算有（自由）能垒，（自由）能垒太低也没什么意义，因为这一步未必是决速步。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 |
alanzhou 发表于 2023-2-5 23:15 可不可靠看你的需求。应该考虑某种计算机别的一般误差，总是要在考虑误差范围的基础上考虑问题。考虑到你优化级别的误差，与单点级别相差不到1 kcal/mol很正常。如果你需要研究的课题是所优化得到的TS和底物这两个结构谁能量低，那优化级别就是不可靠的；但如果你要研究的课题是这个反应在一般溶液条件下的反应速率，那优化级别在这一步中的误差对研究就没什么显著影响，因为这个反应的反应速率肯定是其他步骤控制的（或者是扩散控制的），1 kcal/mol的活化能不管是有还是没有、对一般溶液反应的动力学研究没啥实际意义。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 |
alanzhou 发表于 2023-2-5 23:11 个人猜测是，这是个熵增反应，所以自由能垒低于电子能垒。依据电子能垒判断，存在过渡态；依据自由能垒，或许没有。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 |
liyuanhe211 发表于 2023-2-5 21:27 谢谢老师的分析！请问老师您是认为这个优化得到的中间体与过渡态能量差距较小，从而确实可能在不同计算级别下能量差距不同，甚至翻转吗？那这样的话是不是用低级别泛函计算得到的过渡态、中间体结构，甚至反应路径的分析也不是十分的可靠呢？ |
chands 发表于 2023-2-5 19:05 谢谢老师！请问您是如何分析反应可能没有过渡态的呢？这个几何结构的频率计算结果确实只有一个虚频，跑IRC也确实连接了反应物（另一端虽然连接了中间体，但是确实是产物生成的方向）。 |
|
IRC上的过渡态和底物之间一共只有4 kJ/mol的活化能，换了级别之后顺序翻过来也不很奇怪。 一个与此无关的问题是过渡态转一下CCCC二面角估计会更优。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-22 04:12 , Processed in 0.184405 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites