计算化学公社

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

[ORCA] 溶解自由能修正问题

[复制链接 Copy URL]

108

帖子

0

威望

286

eV
积分
394

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
请教各位老师,我计算了两个反应 1: A+B — A-B 以及 2: A+C— A-C 的活化能, A-B和A-C 分别是两个反应的过渡态结构,实验结果是反应1的反应温度为60℃,反应2的反应温度为 0℃左右,用B3LYP SVP 方法计算得到的反应活化能 ,气相优化后得到的自由能能垒的结果是 反应1比反应2高1Kcal/mol , 然后用 M062X 6-31G*计算了溶解自由能,溶解自由能算完修正后,溶剂下自由能能垒反应1比反应2低6Kcal/mol .我用的能垒计算方法:结构优化后 分别对反应物A、B 以及过渡态A-B用! M062X 6-31G*计算单点能,再用! M062X 6-31G* % cosmo smd true  solvent "THF" end 计算反应物A、B 以及过渡态A-B 溶剂下单点能 ,用{溶剂下单点能-气相单点能+1.89Kcal/mol } 分别加到各自气相freq算出的自由能上,再用过渡态的自由能-反应物自由能的和 得到能垒。
过渡态结构的自由能 我也是把优化后过渡态的坐标拿出来,直接用 M062X 6-31G*计算气相和溶剂下的单点。
本来想用TZVP 但是TZVP算的特别慢  想用SVP先试一下。
我的溶剂模型计算方法有没有弄错什么??为什么我在气相中反应1的能垒高,溶剂中却反过来了?
跟以前请教大家的这个帖子里算的是一样的物质http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3311&extra=





6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

2#
发表于 Post on 2016-6-21 18:40:26 | 只看该作者 Only view this author
为了清晰,我再说一下计算流程,以A + B -> A-B为例

1 用m062X/def2-SVP优化A···B复合物结构,并做振动分析得到气相下自由能
2 用m062X/def2-SVP优化A-B过渡态结构,并做振动分析得到气相下自由能
3 用m052X/6-31G*对1的结构计算气相单点,用m052X/6-31G*结合SMD计算溶液下单点。两个单点求差得到复合物结构的溶解自由能
4 用m052X/6-31G*对2的结构计算气相单点,用m052X/6-31G*结合SMD计算溶液下单点。两个单点求差得到过渡态的溶解自由能
将3、4得到的溶解自由能都加上1.89kcal/mol,再加到1、2得到的气相下自由能,就得到了复合物和过渡态的溶液下的自由能。然后求差得到势垒。

如果结果还是和实验冲突,改用def2TZVP计算复合物和过渡态的单点,与m062X/def2-SVP的自由能热校正量相加得到气相自由能,或者索性用PWPB95/def2-TZVP计算单点。
如果还是不好,就不是计算级别或者计算流程的问题了,可能过渡态不对、计算模型有问题,或者没考虑一些其它因素之类。
北京科音自然科学研究中心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

108

帖子

0

威望

286

eV
积分
394

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2016-6-22 10:09:21 | 只看该作者 Only view this author
sobereva 发表于 2016-6-21 18:40
为了清晰,我再说一下计算流程,以A + B -> A-B为例

1 用m062X/def2-SVP优化A···B复合物结构,并做振 ...

谢谢sob老师,之前我算的时候没有计算“A···B复合物”,用的是A和B各自的结构,计算能垒时的自由能用的是A和B各自的自由能相加起来。再请教sob老师,怎么优化“A···B复合物”?是不是把A 和B 的坐标放在一起,然后用结构优化的命令让它自己去寻找A···B复合物的结构??

6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

4#
发表于 Post on 2016-6-22 13:23:25 | 只看该作者 Only view this author
dingliping1108 发表于 2016-6-22 10:09
谢谢sob老师,之前我算的时候没有计算“A···B复合物”,用的是A和B各自的结构,计算能垒时的自由能用 ...

是。
或者你找到过渡态之后,走IRC,对IRC在反应物那一端的结构进行优化。
北京科音自然科学研究中心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

108

帖子

0

威望

286

eV
积分
394

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-6-22 20:33:59 | 只看该作者 Only view this author
sobereva 发表于 2016-6-22 13:23
是。
或者你找到过渡态之后,走IRC,对IRC在反应物那一端的结构进行优化。

谢谢sob老师,orca里面没有IRC,我把两个反应物放到一起进行基态优化,但是发现两个反应物 越跑离得越远,有什么更好的办法寻找复合物结构吗?

6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

6#
发表于 Post on 2016-6-23 02:06:27 | 只看该作者 Only view this author
dingliping1108 发表于 2016-6-22 20:33
谢谢sob老师,orca里面没有IRC,我把两个反应物放到一起进行基态优化,但是发现两个反应物 越跑离得越远 ...

说明初猜结构太差
要凭化学直觉
北京科音自然科学研究中心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

108

帖子

0

威望

286

eV
积分
394

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2016-6-23 08:52:39 | 只看该作者 Only view this author
本帖最后由 dingliping1108 于 2016-6-23 09:50 编辑
sobereva 发表于 2016-6-23 02:06
说明初猜结构太差
要凭化学直觉

请教sob老师,如果复合物的结构实在找不到,是不是也可以用自由两个反应物自由能之和来代替复合物的自由能。

6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

8#
发表于 Post on 2016-6-24 03:41:21 | 只看该作者 Only view this author
dingliping1108 发表于 2016-6-23 08:52
请教sob老师,如果复合物的结构实在找不到,是不是也可以用自由两个反应物自由能之和来代替复合物的自由 ...

复合物结构肯定是存在的。即便不做IRC,也可以按照虚频的振动矢量方式将结构移动一点然后再优化。如果死活找不到复合物结构就只能按你说的用反应物单体加和做近似了
北京科音自然科学研究中心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

108

帖子

0

威望

286

eV
积分
394

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2016-6-24 11:01:29 | 只看该作者 Only view this author
sobereva 发表于 2016-6-24 03:41
复合物结构肯定是存在的。即便不做IRC,也可以按照虚频的振动矢量方式将结构移动一点然后再优化。如果死 ...

谢谢sob老师

36

帖子

0

威望

1770

eV
积分
1807

Level 5 (御坂)

10#
发表于 Post on 2016-9-22 20:48:07 | 只看该作者 Only view this author
sob老师,如果一个反应是 aA +bB -> cC+dD,那么算该反应的自由能还是分别优化A和B形成的复合物,C和D形成的复合物的自由能,然后再求差值吗? 这种情况下不考虑反应系数吗?还是需要分别单独计算A、B、C和D的自由能,然后再根据各自的反应系数算自由能呢?这点我一直存在疑惑,请教sob老师!!

6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

11#
发表于 Post on 2016-9-22 21:15:40 | 只看该作者 Only view this author
花开半夏 发表于 2016-9-22 20:48
sob老师,如果一个反应是 aA +bB -> cC+dD,那么算该反应的自由能还是分别优化A和B形成的复合物,C和D形成 ...

看你计算目的。
如果是要和实验对比,用每个单体单独优化的结构,单独计算单体自由能
如果从理论上研究复合物之间的变化,就用复合物优化的结构,计算复合物自由能
北京科音自然科学研究中心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

36

帖子

0

威望

1770

eV
积分
1807

Level 5 (御坂)

12#
发表于 Post on 2016-10-3 15:28:57 | 只看该作者 Only view this author
sob老师,谢谢您的解答。最近我一直在用热力学组合的方法进行气相下吉布斯自由能的计算,又遇到了下面几个问题:
1. 我算的反应式里面是包含Pb和I 原子的,用G4一直算不出来。我查了相关的文献,发现 The Gaussian-4 theory (G4 theory) for the calculation of energies of compounds containing first-(Li–F), second-(Na–Cl), and third-row main group (K, Ca, and Ga–Kr) atoms is presented. 对于涉及到周期表中第五周期及之后的元素,是否就不能再使用热力学组合的方法了?
2.我用G4方法计算了H2和O2复合物的自由能  G4 Free Energy= -151.425752 , H2O的自由能为G4 Free Energy= -76.414891 ,这样的话生成水的自由能该如何计算呢?怎样将反应系数考虑进去呢?


6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

13#
发表于 Post on 2016-10-3 22:35:08 | 只看该作者 Only view this author
花开半夏 发表于 2016-10-3 15:28
sob老师,谢谢您的解答。最近我一直在用热力学组合的方法进行气相下吉布斯自由能的计算,又遇到了下面几个 ...

1 对
2 2(H2)+O2=2(H2O)
算复合物的自由能没有意义。直接算水、H2、O2的自由能,按照上式一算就出来了,就是基本的热力学知识
北京科音自然科学研究中心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

36

帖子

0

威望

1770

eV
积分
1807

Level 5 (御坂)

14#
发表于 Post on 2016-10-5 19:56:16 | 只看该作者 Only view this author
多谢sob老师指点。sob老师,我在算一个含有Pb原子和I原子的反应的吉布斯自由能,体系在25个原子左右,
我用赝势基组 b3lyp/genecp

H C N 0
6-311g**
****
Pb I 0
Lanl2DZ
****

Pb I 0
Lanl2DZ

进行优化,然后用b2plyp/def2tzvp 进行单点能的计算,这样选择可以吗?另外,对于混合基组是没有校正因子的,那么该如何解决呢?

6万

帖子

99

威望

6万

eV
积分
125138

管理员

公社社长

15#
发表于 Post on 2016-10-5 21:04:00 | 只看该作者 Only view this author
花开半夏 发表于 2016-10-5 19:56
多谢sob老师指点。sob老师,我在算一个含有Pb原子和I原子的反应的吉布斯自由能,体系在25个原子左右,
我 ...

可以。
校正因子就用大多数原子用的基组的校正因子。
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

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

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