计算化学公社

标题: 溶解能是正的怎么处理？ [打印本页]

作者
Author: gauss98 时间: 2022-2-17 23:11
标题: 溶解能是正的怎么处理？
计算发现一个物种溶解能是正的，比真空能量高一个多千卡，有什么物理含义？是表示难溶解么
这种情况在反应曲线中是不是应该以真空状态的能量来计算？
谢谢

作者
Author: wzkchem5 时间: 2022-2-17 23:44
说明当这个物种的分压（蒸汽压）等于1atm时，它在溶液里的平衡浓度小于1mol/L。
这种情况下按真空算和按溶解态算都行，但是当按真空算时，应该明确说明对这个物种用的是1atm理想气体的参考态，以防止别人误认为你用的是1mol/L溶解态作为参考态。

作者
Author: gauss98 时间: 2022-2-18 00:04
谢谢
高斯里有溶剂效应关键词时给出的自由能就是1mol/l么？
我看社长关于做溶剂效应的博文里专门说要溶解能加上1.89千卡来转换气相跟1mol/l浓度？

作者
Author: wzkchem5 时间: 2022-2-18 04:53

gauss98 发表于 2022-2-17 17:04
谢谢
高斯里有溶剂效应关键词时给出的自由能就是1mol/l么？
我看社长关于做溶剂效应的博文里专门说要溶解 ...

程序给的是1atm的（还是1bar？忘了），需要自己换算

作者
Author: sobereva 时间: 2022-2-18 07:58
应当说明具体怎么算的。下文说的用隐式溶剂模型算对应的1 M的溶解自由能
使用Multiwfn预测晶体密度、蒸发焓、沸点、溶解自由能等性质
http://sobereva.com/337（http://bbs.keinsci.com/thread-3571-1-1.html）
溶解自由能为正体现溶质和溶剂亲和性极差
一般实际是什么环境下发生反应就在什么环境下算，也省得额外解释。
注意分清楚kcal、cal、kcal/mol

作者
Author: gauss98 时间: 2022-2-18 10:00
本帖最后由 gauss98 于 2022-2-18 10:33 编辑

sobereva 发表于 2022-2-18 07:58
应当说明具体怎么算的。下文说的用隐式溶剂模型算对应的1 M的溶解自由能
使用Multiwfn预测晶体密度、蒸发 ...

谢谢社长。
我算的就是“巨懒”的#p B3LYP/def2tzvp Em=GD3BJ OPT scrf=(smd,solvent=toluene) freq temperature=373.15 nosym
找到
Sum of electronic and thermal Free Energies=
得到的自由能
跟社长推荐的溶解态几何结构下的真空电子能+溶解态的热力学校正值+M052X 溶解态电子能-M052X电子能
二者相差0.05~0.3kcal/mol。
我疑问的是，这个自由能是 1M 的自由能还是 1atm 的自由能？

另外，试用shermo 2.3.4版本，根据说明设置 settings.ini conc=1.0M 应该有1atm——>1M的自由能转换值
我在输出里始终没有找到这一项。（执着于这个是要用于concvar的气液体系的速度模拟）

作者
Author: sobereva 时间: 2022-2-18 16:31

gauss98 发表于 2022-2-18 10:00
谢谢社长。
我算的就是“巨懒”的#p B3LYP/def2tzvp Em=GD3BJ OPT scrf=(smd,solvent=toluene) freq ...

严格来说既不是1M也不是1atm
本质上可以当成的B3LYP-D3(BJ)/def2-TZVP的1atm气相自由能加上B3LYP-D3(BJ)/def2-TZVP结合SMD算的1M的溶解自由能。
虽然这在逻辑上比较怪，但inconsistency在求差时很大程度会抵消，而且这种inconsistency相对于其它各种各样的误差来源并不显著

欢迎光临计算化学公社 (http://bbs.keinsci.com/)