溶解自由能的计算结果与实验对比问题

Eliza-NZH · 发表于 Post on 2018-4-23 16:14:20

老师好，我这边发现我的计算的溶解自由能与实验结果是相反的，我这边计算的一个有机小分子在不同溶剂里面的溶解度，实验上，常温下，它在二氯甲烷（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。以上为我的问题和计算过程，请老师指导我问题出在哪里，谢谢！

Graphite · 发表于 Post on 2018-4-23 16:43:27

请教一下, 请问你SMD模型是怎么实现混合溶剂输入的呢?

Eliza-NZH · 发表于 Post on 2018-4-23 17:58:40

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

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

ggdh · 发表于 Post on 2018-4-23 18:30:44

1没有opt么
2没有算纯EC DMC中的溶解自由能么
3把DCM EC DMC EC/DMC 的参数都列出来看看

Eliza-NZH · 发表于 Post on 2018-4-23 18:52:18

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

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

ggdh · 发表于 Post on 2018-4-23 19:48:49

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

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

mutron · 发表于 Post on 2018-4-23 21:20:53

本帖最后由 mutron 于 2018-4-24 00:11 编辑

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

Eliza-NZH · 发表于 Post on 2018-4-23 21:33:41

本帖最后由 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。

liyuanhe211 · 发表于 Post on 2018-4-24 02:22:17

SMD那么多自定义参数哪儿来的？感觉计算过程的问题

刘士晗 · 发表于 Post on 2018-12-26 22:49:29

老师您好，我想问一下溶解自由能越负越容易溶解，那我怎么界定他不溶呢，如果是正值就说明他不溶吗，还是说他高于某个数值溶解度就会变差。

sobereva · 发表于 Post on 2018-12-27 07:23:36

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

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

刘士晗 · 发表于 Post on 2018-12-27 11:37:06

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

谢谢sober老师！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 溶解自由能的计算结果与实验对比问题