计算化学公社

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

[新手求助] 求助过渡态吉布斯自由能比初态更低

[复制链接 Copy URL]

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
一个H转移过程,我用M052X/6-31G*找出过渡态,再走了IRC曲线连接了初末态。然后对IRC得到的初末态进行优化,算能量。得到的从初态到过渡态居然放热-14.56kcal/mol.(M052X/6-31G*优化,M052X/DEF2-TZVP算能量)
然后我用M052X/TZVP进行的计算得到初态到过渡态放热-4.49kcal/mol.(M052X/tzvp优化,M052X/DEF2-TZVP算能量)
是因为涉及到H的转移,所以用6-31g*没有考虑H的极化不合理吗?但是用tzvp算出来的吉布斯自由能过渡态的还是比初态更低。优化都收敛没有虚频,IRC出来的曲线都是正常的,过渡态最高分别往左右下降

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

2#
发表于 Post on 2021-7-30 16:22:26 | 只看该作者 Only view this author
本帖最后由 wzkchem5 于 2021-7-30 09:23 编辑

检查一下是不是单点能计算没有收敛到稳定波函数。另外M052X早就过时了,只有在计算SMD溶解自由能这一个场合有一定理由用,其他所有情况下都不如M062X。
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
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-7-30 20:11:21 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-30 16:22
检查一下是不是单点能计算没有收敛到稳定波函数。另外M052X早就过时了,只有在计算SMD溶解自由能这一个场合 ...

我做了初态的stable,显示波函数stable了,所以是我泛函的原因吗?我做的是一个有机体系

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

4#
发表于 Post on 2021-7-30 20:28:06 | 只看该作者 Only view this author
Pamela 发表于 2021-7-30 13:11
我做了初态的stable,显示波函数stable了,所以是我泛函的原因吗?我做的是一个有机体系

是不是初态是两个分子?如果是的话,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
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-7-30 20:59:36 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-30 20:28
是不是初态是两个分子?如果是的话,IRC看似走到了类似初态的结构,实际可能是走到反应复合物。这时应该 ...

对啊 我就是把IRC的跑出来的初态整体优化算能量,并不是把两个单体的能量叠加的,我那算的就是是反应复合物得能量吧? 还有我用M062X跑了一个IRC跑了一周都没跑完,体系有七十几个原子

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

6#
发表于 Post on 2021-7-30 22:19:40 | 只看该作者 Only view this author
Pamela 发表于 2021-7-30 13:59
对啊 我就是把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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

7#
发表于 Post on 2021-8-1 06:11:58 | 只看该作者 Only view this author
此情况即便优化也必须给氢,至少是迁移的那个氢,带极化
但凡问这类问题,都必须说清楚过渡态相对于反应物的电子能量差和自由能差各是多少,由此搞清楚是电子能量能垒层面就有问题还是自由能的热力学校正量变化所导致的
北京科音自然科学研究中心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

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-8-2 17:01:09 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-30 22:19
在结构优化的理论级别下,初态到过渡态能量变化多少?

是看第一步优化出来的Sum of electronic and thermal Free Energies吗?用6-31g*优化的初态到过渡态能量差是-13多kcal/mol,tzvp优化出来的初态到过渡态是-6kcal/mol的能量。

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

9#
发表于 Post on 2021-8-2 17:09:00 | 只看该作者 Only view this author
Pamela 发表于 2021-8-2 10:01
是看第一步优化出来的Sum of electronic and thermal Free Energies吗?用6-31g*优化的初态到过渡态能量 ...

我说的是电子能量。如果在优化的理论级别下,电子能仍然是过渡态比反应物低,那么唯一的解释是这个过渡态不是和反应物直接相连的,建议仔细检查IRC到底是不是跑到了反应物,还是跑到了和反应物相似但是和反应物不等价的另一个结构。
如果IRC跑完了,且跑到了正确的反应物,那么IRC反应物侧的能量应该等于优化的理论级别下反应物的电子能,IRC过渡态能量等于优化的理论级别下过渡态的电子能。又因为IRC的基本原理,IRC过渡态能量一定高于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
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2021-8-2 17:12:34 | 只看该作者 Only view this author
sobereva 发表于 2021-8-1 06:11
此情况即便优化也必须给氢,至少是迁移的那个氢,带极化
但凡问这类问题,都必须说清楚过渡态相对于反应物 ...

好的,谢谢sob老师
我用的是溶剂模型,电子能量就是读取第二歩单点能能量吧。用M052X/6-31G*优化,M052X/DET2TZVP算单点能,初态到过渡态电子能量差是-7.34kcal/mol。用M052X/tzvp优化,M052X/DET2TZVP算单点能,初态到过渡态电子能量差是0.11kcal/mol,所以tzvp优化的是由于自由能热力学矫正导致的。请问sob老师该怎么解决呢?

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2021-8-2 17:22:40 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-2 17:09
我说的是电子能量。如果在优化的理论级别下,电子能仍然是过渡态比反应物低,那么唯一的解释是这个过渡态 ...

您的意思是我可以用优化的级别 比如M052X/6-31G*算一下反应物和过渡态的单点能进行对比是吗?

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

12#
发表于 Post on 2021-8-2 19:12:42 | 只看该作者 Only view this author
Pamela 发表于 2021-8-2 10:22
您的意思是我可以用优化的级别 比如M052X/6-31G*算一下反应物和过渡态的单点能进行对比是吗?

对,不过不用单独算单点能,从结构优化的输出文件里读取就行了
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

13#
发表于 Post on 2021-8-2 20:09:31 | 只看该作者 Only view this author
Pamela 发表于 2021-8-2 17:12
好的,谢谢sob老师
我用的是溶剂模型,电子能量就是读取第二歩单点能能量吧。用M052X/6-31G*优化,M052X ...

如果电子能垒为轻微正值,自由能垒却为负,说明没有真正意义的过渡态,是个无势垒过程

评分 Rate

参与人数
Participants 2
eV +10 收起 理由
Reason
Xx_DerEnTh + 5 找到答案
超限制抱怨 + 5 对我有用

查看全部评分 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

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2021-8-2 21:13:37 | 只看该作者 Only view this author
sobereva 发表于 2021-8-2 20:09
如果电子能垒为轻微正值,自由能垒却为负,说明没有真正意义的过渡态,是个无势垒过程

谢谢sob老师,所以用6-31g*优化这个体系还是误差太大了哦。sob老师我现在用M052X这个泛函进行优化和算单点能以及溶解自由能,我的是有机体系,用M062X我的体系跑得太慢了。请问老师我用M052X合理吗?

4290

帖子

4

威望

9558

eV
积分
13928

Level 6 (一方通行)

MOKIT开发者

15#
发表于 Post on 2021-8-2 21:18:44 | 只看该作者 Only view this author
Pamela 发表于 2021-8-2 21:13
谢谢sob老师,所以用6-31g*优化这个体系还是误差太大了哦。sob老师我现在用M052X这个泛函进行优化和算单 ...

M062X与M052X耗时没什么区别。
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 03:36 , Processed in 0.335424 second(s), 21 queries , Gzip On.

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