计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 谈谈隐式溶剂模型下溶解自由能和体系自由能的计算

  [复制链接 Copy URL]

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

136#
发表于 Post on 2021-7-15 22:32:20 | 只看该作者 Only view this author
guoxh 发表于 2021-7-15 13:39
sob老师:

您好,我是一个学习理论计算的硕士研究生,可以请教您一下关于分子溶解自由能计算的一些细节 ...

看你贴的链接的2(6)节就知道了,溶剂化熵不需要单独考虑。
但是审稿人说的第一配位层是否要加,这个需要结合你的体系来考虑,如果溶质和溶剂分子之间有比较强的氢键、静电相互作用或配位作用,那么需要用显式溶剂模型考虑第一配位层
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

137#
 楼主 Author| 发表于 Post on 2021-7-16 15:52:13 | 只看该作者 Only view this author
guoxh 发表于 2021-7-15 20:39
sob老师:

您好,我是一个学习理论计算的硕士研究生,可以请教您一下关于分子溶解自由能计算的一些细节 ...

在文章里强调清楚你的溶解自由能的计算过程,并且引用Truhlar的SMD的原文。
并且说明,溶解自由能里面的熵效应在SMD模型拟合非极性部分参数的时候已经等效体现在参数里面了,所以直接用SMD模型按照作者的方式(即我博文里的方式)去计算溶解自由能,就已经能恰当地把熵效应等效表现出来了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

15

帖子

0

威望

151

eV
积分
166

Level 3 能力者

138#
发表于 Post on 2021-7-16 16:16:25 | 只看该作者 Only view this author
sobereva 发表于 2021-7-16 15:52
在文章里强调清楚你的溶解自由能的计算过程,并且引用Truhlar的SMD的原文。
并且说明,溶解自由能里面的 ...

好哒,谢谢老师

15

帖子

0

威望

151

eV
积分
166

Level 3 能力者

139#
发表于 Post on 2021-7-16 16:25:28 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-15 22:32
看你贴的链接的2(6)节就知道了,溶剂化熵不需要单独考虑。
但是审稿人说的第一配位层是否要加,这个需要 ...

好哒,谢谢

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

140#
发表于 Post on 2021-8-5 10:02:01 | 只看该作者 Only view this author
sob老师,我想请教一下,如果我想要研究的是单个离子(比如钠离子),在水溶剂中逐步与水分子形成团簇Na(H2O)n+过程中每一步的自由能变化,那么在第一步(Na+)+H2O=Na(H2O)+这个反应的自由能计算过程中,是否能对Na+采用您这篇帖子里介绍的方法结合scrf(smd,solvent=water)计算其溶解自由能呢,如果不行的话想请问老师怎么计算(Na+)+H2O=Na(H2O)+这个反应在水溶液中的自由能变化呢

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

141#
 楼主 Author| 发表于 Post on 2021-8-6 08:06:02 | 只看该作者 Only view this author
Qing 发表于 2021-8-5 10:02
sob老师,我想请教一下,如果我想要研究的是单个离子(比如钠离子),在水溶剂中逐步与水分子形成团簇Na(H2 ...

这实际上只是一个完全虚构的过程。对于当前目的,考察每一步的自由能变化,倒是可以
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

142#
发表于 Post on 2021-8-6 08:26:17 | 只看该作者 Only view this author
sobereva 发表于 2021-8-6 08:06
这实际上只是一个完全虚构的过程。对于当前目的,考察每一步的自由能变化,倒是可以

好的,谢谢sob老师

269

帖子

0

威望

1758

eV
积分
2027

Level 5 (御坂)

143#
发表于 Post on 2022-1-19 09:34:05 | 只看该作者 Only view this author
请教各位:卢老师说“计算某分子在乙醇中的溶解自由能,就用比如# M052X/6-31G* scrf(SMD,solvent=ethanol)计算得到的单点能,减去# M052X/6-31G*计算得到的单点能”,那么用ORCA是不是可以这样计算:用DLPNO-CCSD(T) + SMD得到的单点能(优化时用SMD,计算DLPNO-CCSD(T)时也用SMD),减去不用SMD得到的DLPNO-CCSD(T)单点能(优化时不用SMD,计算DLPNO-CCSD(T)时也不用SMD),就得到溶解自由能?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

144#
 楼主 Author| 发表于 Post on 2022-1-19 09:55:17 | 只看该作者 Only view this author
dnlx 发表于 2022-1-19 09:34
请教各位:卢老师说“计算某分子在乙醇中的溶解自由能,就用比如# M052X/6-31G* scrf(SMD,solvent=ethanol) ...

不建议这样
文中明确说了SMD在M05-2X/6-31G*下计算溶解自由能精度是最好的,结合DLPNO-CCSD(T)来用精度如何没法估计,被质疑都不好回答

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
dnlx + 5 懂了,谢谢社长!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

40

帖子

0

威望

565

eV
积分
605

Level 4 (黑子)

145#
发表于 Post on 2022-2-23 20:18:00 | 只看该作者 Only view this author
老师您好,我研究的是一个电子包合物,文献中“需要通过PCM 模型采用介电常数ε = 3.2 来模拟周围的冰环境”来看该条件下它的稳定情况,想问下老师这里应该如何设置这个冰环境?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

146#
 楼主 Author| 发表于 Post on 2022-2-24 11:26:56 | 只看该作者 Only view this author
Lyy1 发表于 2022-2-23 20:18
老师您好,我研究的是一个电子包合物,文献中“需要通过PCM 模型采用介电常数ε = 3.2 来模拟周围的冰环境 ...

自定义溶剂
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

34

帖子

0

威望

482

eV
积分
516

Level 4 (黑子)

147#
发表于 Post on 2022-3-1 17:02:43 | 只看该作者 Only view this author
本帖最后由 作业抄 于 2022-3-1 17:06 编辑

老师您好,我所研究的方向为钠离子电池电解液方向,比较常用到的就是Na+的脱溶剂化能。您在博文里所提到的物质的溶解自由能可以理解为我的Na盐溶解到电解液溶剂中并与周围的溶剂分子形成溶剂化结构所需要的能量,而脱溶剂化过程是Na+与周围的溶剂化分子分离。我想问的是,这里的溶解自由能与脱溶剂化过程所需要的能量是否是互为相反数,谢谢老师

34

帖子

0

威望

482

eV
积分
516

Level 4 (黑子)

148#
发表于 Post on 2022-3-1 17:10:39 | 只看该作者 Only view this author
本帖最后由 作业抄 于 2022-3-1 17:12 编辑

请问老师,利用隐形模型是否可以算出不同温度下的溶解自由能,博文中的两种说法是否是矛盾的

202203011710105516..png (11.29 KB, 下载次数 Times of downloads: 74)

202203011710105516..png

202203011710349561..png (16.19 KB, 下载次数 Times of downloads: 75)

202203011710349561..png

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

149#
发表于 Post on 2022-3-1 17:19:11 | 只看该作者 Only view this author
作业抄 发表于 2022-3-1 10:10
请问老师,利用隐形模型是否可以算出不同温度下的溶解自由能,博文中的两种说法是否是矛盾的

不矛盾。第一张图是说得到精确的不同温度下的溶解自由能没戏,第二张图是说,在假设溶解自由能受温度的影响可忽略不计的情况下,可以近似地得到总自由能。两张图说的一个是溶解自由能,一个是总自由能,而且一个是不能精确计算、一个是可以近似计算。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

34

帖子

0

威望

482

eV
积分
516

Level 4 (黑子)

150#
发表于 Post on 2022-3-1 17:28:17 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-1 17:19
不矛盾。第一张图是说得到精确的不同温度下的溶解自由能没戏,第二张图是说,在假设溶解自由能受温度的影 ...

谢谢老师,我明白了

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 09:36 , Processed in 0.217954 second(s), 26 queries , Gzip On.

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