计算化学公社

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

[新手求助] 吉布斯自由能计算

[复制链接 Copy URL]

125

帖子

1

威望

1047

eV
积分
1192

Level 4 (黑子)

本帖最后由 fankui1990129 于 2015-6-11 15:20 编辑

本人新手。如图所示,我想计算这个反应在水溶液中的吉布斯自由能(副产物是水)

搭建模型的是我将A+B放在一起进行opt+freq计算,同样将C+H2O放在一起进行opt+freq计算。在不加溶剂模型的时候一切正常,加完之后就不收敛了。昨天在论坛进行了咨询,尝试之后效果不大。

我现在想问一下,能不能将A,B,C和水分别进行opt+freq计算。如果这样操作得出的结果是完全不可以信呢?还是存在误差呢?
(因为我是以有机合成实验为主,我们在实验的过程中发现了一个现象,希望用计算化学说明。所以,出现误差也不是很重要,只要能和实验结果定性符合就行,不要定量的完全拟合。)

untitled.gif (1.46 KB, 下载次数 Times of downloads: 114)

untitled.gif

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

2#
发表于 Post on 2015-6-11 14:04:47 来自手机 | 只看该作者 Only view this author
不就是应该分别计算吗 放到一起计算应该不合理吧

125

帖子

1

威望

1047

eV
积分
1192

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2015-6-11 14:55:09 | 只看该作者 Only view this author
本帖最后由 fankui1990129 于 2015-6-11 15:20 编辑
lastzealot 发表于 2015-6-11 14:04
不就是应该分别计算吗 放到一起计算应该不合理吧

谢谢,也就是说原料之间的弱作用忽略不计了么?

因为这次我们必须考虑产物和水之间的弱作用,这是我们研究的重点,所以产物和水最好能要放一起。后来转念一想,水又是溶剂,所以加入溶剂模型的时候就应该会有影响了

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

4#
发表于 Post on 2015-6-11 18:55:15 | 只看该作者 Only view this author
实际上并不需要在溶剂模型下优化结构。

在气相下优化A+B复合物并做振动分析,得到气相下自由能,然后在这个优化好的结构下计算溶剂化自由能(带着溶剂模型和不带情况下的单点能的差值),将二者相加,即是溶剂下的自由能。对于H2O+C也这么计算。最后将反应物和产物求差即可。

一些人认为在气相下而非溶剂模型下优化和做振动分析并不是近似,反倒是必要的,如J. Phys. Chem. A 2010, 114, 13442–13444
北京科音自然科学研究中心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

125

帖子

1

威望

1047

eV
积分
1192

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-6-11 21:17:11 | 只看该作者 Only view this author
sobereva 发表于 2015-6-11 18:55
实际上并不需要在溶剂模型下优化结构。

在气相下优化A+B复合物并做振动分析,得到气相下自由能,然后在 ...

非常感谢!初学者,懂得不多,自学很难。最近多亏了sob老师!!

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

6#
发表于 Post on 2015-6-11 21:38:40 | 只看该作者 Only view this author
sobereva 发表于 2015-6-11 18:55
实际上并不需要在溶剂模型下优化结构。

在气相下优化A+B复合物并做振动分析,得到气相下自由能,然后在 ...

请问初始结构怎么设计?两个反应物的距离、构象如何设计?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

7#
发表于 Post on 2015-6-11 21:59:11 | 只看该作者 Only view this author
lastzealot 发表于 2015-6-11 21:38
请问初始结构怎么设计?两个反应物的距离、构象如何设计?

摆一起,优化即可。初始结构不合适经过优化后也就自然得到实际二聚体构型了。

但很可能有多种二聚体构型,如果自己拿不准哪种二聚体构型最稳定就必须摆各种初猜结构分别优化,或者利用molclus来自动处理了(使用molclus程序做分子团簇构型搜索和分子构象搜索http://bbs.keinsci.com/forum.php?mod=viewthread&tid=577
北京科音自然科学研究中心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

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

8#
发表于 Post on 2015-6-11 22:45:07 | 只看该作者 Only view this author
sobereva 发表于 2015-6-11 21:59
摆一起,优化即可。初始结构不合适经过优化后也就自然得到实际二聚体构型了。

但很可能有多种二聚体构 ...

谢谢老师

21

帖子

0

威望

378

eV
积分
399

Level 3 能力者

9#
发表于 Post on 2017-5-3 17:23:21 | 只看该作者 Only view this author
sobereva 发表于 2015-6-11 18:55
实际上并不需要在溶剂模型下优化结构。

在气相下优化A+B复合物并做振动分析,得到气相下自由能,然后在 ...

sob老师,那做优化和做单点的方法和基组要保持一致,对么?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

10#
发表于 Post on 2017-5-4 11:17:56 | 只看该作者 Only view this author
qxx 发表于 2017-5-3 17:23
sob老师,那做优化和做单点的方法和基组要保持一致,对么?

完全不需要一致
北京科音自然科学研究中心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

21

帖子

0

威望

378

eV
积分
399

Level 3 能力者

11#
发表于 Post on 2017-5-4 14:47:08 | 只看该作者 Only view this author
sobereva 发表于 2017-5-4 11:17
完全不需要一致

谢谢老师!

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 05:27 , Processed in 0.171950 second(s), 23 queries , Gzip On.

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