计算化学公社

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

[Gaussian/gview] 求助:IRC优化的产物构象不同

[复制链接 Copy URL]

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

我看了一篇算燃烧动力学(低温氧化)的文章,算的是ROO---QOOH不同位点处分子内H转移,下图是作者给出的过渡态,对这两个过渡态俩端IRC再优化得到的反应物(R02-1)构象并不相同,能量也有区别,但是作者没有考虑构象变化的过程。我想问下这种情况能垒怎么确定呢?

4.jpg (271.17 KB, 下载次数 Times of downloads: 33)

4.jpg

2.jpg (205.34 KB, 下载次数 Times of downloads: 20)

2.jpg

SF`(N(F)8PF`_4RBQ~)H}4D.png (155.31 KB, 下载次数 Times of downloads: 26)

SF`(N(F)8PF`_4RBQ~)H}4D.png

E)1EY7H(1}FZ7LGYU9JK52P.png (101.35 KB, 下载次数 Times of downloads: 38)

E)1EY7H(1}FZ7LGYU9JK52P.png

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

2#
发表于 Post on 2021-5-15 14:52:52 | 只看该作者 Only view this author
DOI: 10.1021/ar1000956 参见这篇文章
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
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-5-15 20:40:19 | 只看该作者 Only view this author
wzkchem5 发表于 2021-5-15 14:52
DOI: 10.1021/ar1000956 参见这篇文章

您好,我在这篇文章里似乎没有找到关于这个的解释

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

4#
发表于 Post on 2021-5-15 20:50:09 | 只看该作者 Only view this author
baojiajing 发表于 2021-5-15 13:40
您好,我在这篇文章里似乎没有找到关于这个的解释

这篇文章的结论适用于任何反应,不仅适用于包含构象变化的反应,所以没专门提构象。
按照这篇文章的方法,找到TDI和TDTS,两个的自由能减一下就行了
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
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2021-5-15 21:34:22 | 只看该作者 Only view this author
wzkchem5 发表于 2021-5-15 20:50
这篇文章的结论适用于任何反应,不仅适用于包含构象变化的反应,所以没专门提构象。
按照这篇文章的方法 ...

你好,我看了一下,我觉得燃烧这块可能不是这么处理的,这么做耗时太高了

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

6#
发表于 Post on 2021-5-15 23:49:06 | 只看该作者 Only view this author
baojiajing 发表于 2021-5-15 14:34
你好,我看了一下,我觉得燃烧这块可能不是这么处理的,这么做耗时太高了

耗时不高啊,反应物、过渡态和产物的各个构象都找到了的话,很容易就做出来了啊
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
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2021-5-16 11:14:31 | 只看该作者 Only view this author
wzkchem5 发表于 2021-5-15 23:49
耗时不高啊,反应物、过渡态和产物的各个构象都找到了的话,很容易就做出来了啊

像我这样一个反应(这个图片是C4过渡态),C5位置处OO和C1,C3,C4,C7,C9,C8,C11,C12,C13发生H转移生成OOH,IRC跑了之后优化反应物,得到了9种构象,就是说反应物以九种不同构象发生反应,那应该怎么处理反应物呢,让他不是九种反应物而是一种

TS1.png (270.46 KB, 下载次数 Times of downloads: 17)

TS1.png

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

8#
发表于 Post on 2021-5-16 16:17:18 | 只看该作者 Only view this author
baojiajing 发表于 2021-5-16 04:14
像我这样一个反应(这个图片是C4过渡态),C5位置处OO和C1,C3,C4,C7,C9,C8,C11,C12,C13发生H转移生成OOH ...

取过渡态构象里自由能最低的,和反应物构象里自由能最低的相减就行了,哪怕这个过渡态构象不是这个反应物构象反应的过渡态。
这个方法有两个适用条件,第一是构象互变的自由能垒明显小于反应能垒,第二是你只研究这一步基元反应。
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
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2021-5-16 17:32:44 | 只看该作者 Only view this author
wzkchem5 发表于 2021-5-16 16:17
取过渡态构象里自由能最低的,和反应物构象里自由能最低的相减就行了,哪怕这个过渡态构象不是这个反应物 ...

我这个反应本来就有九种过渡态从一种反应物生成九种产物,反应物我直接取IRC两端自由能最低的那个就可以了吗?这个反应的能垒大概在20kcal左右,各种构象之间能量相差1kcal左右

1万

帖子

0

威望

9849

eV
积分
22083

Level 6 (一方通行)

10#
发表于 Post on 2021-5-16 17:53:15 | 只看该作者 Only view this author
baojiajing 发表于 2021-5-16 10:32
我这个反应本来就有九种过渡态从一种反应物生成九种产物,反应物我直接取IRC两端自由能最低的那个就可以 ...

我这样讲吧
假设你的反应物是A,有A1~A9这几个构象
过渡态有TSa~TSi这9种过渡态,分别反应得到Pa~Pi共9种产物
每个过渡态又有若干种构象,比如TSa有TSa1~TSa5这5种构象,TSb有TSb1~TSb7这7种构象
那么你算A到Pa的反应的自由能垒,应该选TSa1~TSa5里面最低的自由能,减去A1~A9里面最低的自由能。你要是算A到Pb,那就把TSa1~TSa5里面最低的自由能换成TSb1~TSb7里面最低的自由能
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
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

77

eV
积分
94

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2021-5-16 18:00:52 | 只看该作者 Only view this author
wzkchem5 发表于 2021-5-16 17:53
我这样讲吧
假设你的反应物是A,有A1~A9这几个构象
过渡态有TSa~TSi这9种过渡态,分别反应得到Pa~Pi共9 ...

我理解你的意思了,谢谢你

本版积分规则 Credits rule

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

GMT+8, 2026-2-17 04:01 , Processed in 0.168778 second(s), 23 queries , Gzip On.

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