计算化学公社

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

[综合交流] 求助:关于活化能垒比较问题

[复制链接 Copy URL]

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

本帖最后由 beijudw 于 2022-6-23 22:51 编辑

各位老师好,我在使用gaussian计算有机胺吸收CO2以及仲胺两性离子去质子化的活化能垒,得到了不同的一系列活化能垒的比较图,比如仲胺与仲胺+CO2的去质子化之间比较。叔胺与叔胺+CO2+H2O之间比较(胺与胺之间只差1~2个甲基),这样我想请问一下各位老师,可以通过什么方法从理论上去解释:为什么某个胺的能垒更低? 通过multifwn里面的计算静电势图?或者IRI方法图形化?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

2#
发表于 Post on 2022-6-16 14:18:35 | 只看该作者 Only view this author
主要就是静电、位阻和氢键,静电项用静电势看,位阻和氢键用IRI以及类似的其他方法看
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
本团队长期招收研究生,有意者可私信联系

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-6-16 15:22:59 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-16 14:18
主要就是静电、位阻和氢键,静电项用静电势看,位阻和氢键用IRI以及类似的其他方法看

好的,谢谢老师

16

帖子

0

威望

1805

eV
积分
1821

Level 5 (御坂)

4#
发表于 Post on 2022-6-17 09:35:54 | 只看该作者 Only view this author
可以考虑用Distortion-interaction分析方法,把过渡态的自由能分解为不同片段扭曲造成的能量上升以及片段之间相互作用造成的能量降低,来解释同类过渡态中能垒的差异。可参考http://bbs.keinsci.com/thread-23990-1-1.html

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

5#
发表于 Post on 2022-6-17 16:40:37 | 只看该作者 Only view this author
手打100遍Gaussian、Multiwfn记住怎么拼
北京科音自然科学研究中心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

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-6-23 22:50:37 | 只看该作者 Only view this author
baisl 发表于 2022-6-17 09:35
可以考虑用Distortion-interaction分析方法,把过渡态的自由能分解为不同片段扭曲造成的能量上升以及片段之 ...

好的,谢谢老师

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-6-23 22:51:16 | 只看该作者 Only view this author
sobereva 发表于 2022-6-17 16:40
手打100遍Gaussian、Multiwfn记住怎么拼

对不起对不起,问的比较急,打错了,后面会注意的

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-6-27 11:18:53 | 只看该作者 Only view this author
本帖最后由 beijudw 于 2022-6-27 11:28 编辑
wzkchem5 发表于 2022-6-16 14:18
主要就是静电、位阻和氢键,静电项用静电势看,位阻和氢键用IRI以及类似的其他方法看

老师,您好,我在观察对应结构的静电势,发现自由能垒低的胺,反而N孤对电子方向的静电势更高了,可是按道理来说,这里是质子化反应,不应该静电势更负的,能垒更低吗?想请教一下老师该怎么去解释,下图是初态复合物里面的胺结构静电势图:左边为质子化自由能垒更低的

YX@$]Y{6K(X$QU{4(AGS0LH.png (286.99 KB, 下载次数 Times of downloads: 30)

YX@$]Y{6K(X$QU{4(AGS0LH.png

16

帖子

0

威望

1805

eV
积分
1821

Level 5 (御坂)

9#
发表于 Post on 2022-6-27 13:35:05 | 只看该作者 Only view this author
本帖最后由 baisl 于 2022-6-27 13:36 编辑
beijudw 发表于 2022-6-27 11:18
老师,您好,我在观察对应结构的静电势,发现自由能垒低的胺,反而N孤对电子方向的静电势更高了,可是按 ...

这是在你选取的电子密度等值面上的静电势极小值,不一定是整个空间上的极小值。应当考虑静电势的全局极小点,可参考社长的这篇博客。http://sobereva.com/493

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2022-6-27 13:47:15 | 只看该作者 Only view this author
baisl 发表于 2022-6-27 13:35
这是在你选取的电子密度等值面上的静电势极小值,不一定是整个空间上的极小值。应当考虑静电势的全局极小 ...

但是我计算反应的位点是N的孤对电子方向,用来与H2O+CO2发生质子化反应形成AH+ + HCO3- ,不应该去看N的孤对电子附近的静电势来说明吗?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

11#
发表于 Post on 2022-6-27 15:13:42 | 只看该作者 Only view this author
beijudw 发表于 2022-6-27 04:18
老师,您好,我在观察对应结构的静电势,发现自由能垒低的胺,反而N孤对电子方向的静电势更高了,可是按 ...

说明可能有位阻、色散的因素影响,不一定是纯静电问题
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
本团队长期招收研究生,有意者可私信联系

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-6-27 15:18:22 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-27 15:13
说明可能有位阻、色散的因素影响,不一定是纯静电问题

位阻的话怎么看呢,色散影响就是绘制IRI图像,然后计算氢键强度吗?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

13#
发表于 Post on 2022-6-27 15:22:44 | 只看该作者 Only view this author
beijudw 发表于 2022-6-27 08:18
位阻的话怎么看呢,色散影响就是绘制IRI图像,然后计算氢键强度吗?

对,就是IRI、IGMH这类方法,具体什么情况适合用什么方法看http://sobereva.com/621
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
本团队长期招收研究生,有意者可私信联系

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2022-6-27 15:23:47 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-27 15:22
对,就是IRI、IGMH这类方法,具体什么情况适合用什么方法看http://sobereva.com/621

好的 谢谢老师

48

帖子

0

威望

245

eV
积分
293

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2022-6-27 23:04:42 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-27 15:22
对,就是IRI、IGMH这类方法,具体什么情况适合用什么方法看http://sobereva.com/621

老师,您好,对于IRI分析中,是通过对比断键,成键,还有分子间作用力的键能大小来解释能垒高低嘛?比如断键键能更低,成键键能更高,然后还有一定的范德华吸引力,所以能垒更低这样嘛?但是如果有高有低,是不是就比较分子间键能的总和来解释能垒?

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 23:42 , Processed in 1.044720 second(s), 23 queries , Gzip On.

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