计算化学公社

标题: 溶解自由能的计算结果与实验对比问题 [打印本页]

作者
Author: Eliza-NZH 时间: 2018-4-23 16:14
标题: 溶解自由能的计算结果与实验对比问题
老师好，我这边发现我的计算的溶解自由能与实验结果是相反的，我这边计算的一个有机小分子在不同溶剂里面的溶解度，实验上，常温下，它在二氯甲烷（DCM）里面是不溶的，在EC/DMC里面的溶解度比DCM里面好，那为什么我的计算结果在DCM里面的溶解自由能是-20Kcal/mol, 而在EC/DMC里面计算的溶解自由能是-12 Kcal/mol, 根据越负溶解度越好原理，应该DCM溶解度大于EC.DMC。这与实验结果相反。
以下是我的溶解自由能计算过程，输入文件命令行：（1）#p m062x/def2tzvp scrf=(smd,solvent=dichloromethane) ，这条命令计算出来的总能为E1；然后去掉溶剂：（2）：#p m062x/def2tzvp 计算的总能为E2。则溶解自由能G=E1-E2。以上为我的问题和计算过程，请老师指导我问题出在哪里，谢谢！

作者
Author: Graphite 时间: 2018-4-23 16:43
请教一下, 请问你SMD模型是怎么实现混合溶剂输入的呢?

作者
Author: Eliza-NZH 时间: 2018-4-23 17:58

Graphite 发表于 2018-4-23 16:43
请教一下, 请问你SMD模型是怎么实现混合溶剂输入的呢?

根据溶剂体积比设置相关参数，也就是混合溶剂模型。

作者
Author: ggdh 时间: 2018-4-23 18:30
1没有opt么
2没有算纯EC DMC中的溶解自由能么
3把DCM EC DMC EC/DMC 的参数都列出来看看

作者
Author: Eliza-NZH 时间: 2018-4-23 18:52

ggdh 发表于 2018-4-23 18:30
1没有opt么
2没有算纯EC DMC中的溶解自由能么
3把DCM EC DMC EC/DMC 的参数都列出来看看

老师，结构是事先优化好的，在纯EC DMC 里面的溶解自由能我没算呢。

作者
Author: ggdh 时间: 2018-4-23 19:48

Eliza-NZH 发表于 2018-4-23 18:52
老师，结构是事先优化好的，在纯EC DMC 里面的溶解自由能我没算呢。

结构不是应该分别不同的溶剂中优化么？
EC DMC最好都算一下，毕竟有选择性溶剂化的问题。

作者
Author: mutron 时间: 2018-4-23 21:20
本帖最后由 mutron 于 2018-4-24 00:11 编辑

EC/DMC是常用的锂电池电解液，因为EC的偶极距和介电常数都比较大，不少文章用杂化溶剂模型来计算Li离子的溶解自由能。如果你的有机小分子有活泼氢，是否需要考虑有机小分子溶解的时候形成氢键

作者
Author: Eliza-NZH 时间: 2018-4-23 21:33
本帖最后由 Eliza-NZH 于 2018-4-23 22:19 编辑

ggdh 发表于 2018-4-23 19:48
结构不是应该分别不同的溶剂中优化么？
EC DMC最好都算一下，毕竟有选择性溶剂化的问题。

您意思是我应该这样，#p m062x/def2tzvp scrf=(smd,solvent=A) opt ，计算出来的总能为E1？A 为EC或者DMC。

作者
Author: liyuanhe211 时间: 2018-4-24 02:22
SMD那么多自定义参数哪儿来的？感觉计算过程的问题

作者
Author: 刘士晗 时间: 2018-12-26 22:49
老师您好，我想问一下溶解自由能越负越容易溶解，那我怎么界定他不溶呢，如果是正值就说明他不溶吗，还是说他高于某个数值溶解度就会变差。

作者
Author: sobereva 时间: 2018-12-27 07:23

刘士晗发表于 2018-12-26 22:49
老师您好，我想问一下溶解自由能越负越容易溶解，那我怎么界定他不溶呢，如果是正值就说明他不溶吗，还是 ...

溶解自由能不是决定是否溶解的唯一因素。最理想、严格的考察方法是做分子动力学

作者
Author: 刘士晗 时间: 2018-12-27 11:37

sobereva 发表于 2018-12-27 07:23
溶解自由能不是决定是否溶解的唯一因素。最理想、严格的考察方法是做分子动力学

谢谢sober老师！

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