求助在溶剂模型下离子的溶解自由能计算解释溶解性的问题

cyr

      各位老师好！想请教一下，已知在溶剂B中，离子之间结合生成沉淀，而在溶剂A中，是溶液状态。我如果通过计算该离子体系在溶剂B中的溶解自由能小于在溶剂A中的溶解自由能，是否能说明问题。
计算使用的关键词与结果如下：
几何优化：A：#p opt b3lyp/6-311+g(d,p) scrf=(smd,solvent=water)     B： #p opt b3lyp/6-311+g(d,p) scrf=(smd,solvent=generic,read)
能量计算：A：#p M062X/def2TZVP scrf=(smd,solvent=water)              B：#p M062X/def2TZVP scrf=(smd,solvent=generic,read)
结果：

	E	ΔE/hartree	E/(kcal/mol)
A	-763.947	-0.113	-70.908
B	-763.936	-0.10194	-63.9687
vaccum	-763.834

wzkchem5

如果沉淀不包含溶剂分子（如结晶水等），那么可以

cyr

wzkchem5 发表于 2023-10-23 22:00
如果沉淀不包含溶剂分子（如结晶水等），那么可以

好的，谢谢老师！

cyr

wzkchem5 发表于 2023-10-23 22:00
如果沉淀不包含溶剂分子（如结晶水等），那么可以

不好意思，老师，我想问一下如果沉淀中有溶剂分子，那是要把溶剂分子加进去再进行计算，还是说这种方法就不可行了。感谢老师的解答！

wzkchem5

cyr 发表于 2023-10-25 12:53
不好意思，老师，我想问一下如果沉淀中有溶剂分子，那是要把溶剂分子加进去再进行计算，还是说这种方法就 ...

那样需要想办法算沉淀的自由能，例如用支持周期性边界条件的程序做计算。
溶解自由能指的是溶质从气相到溶液相的自由能变化，即便溶质在常温常压下不是气体也是如此。所以你这个方法只有在A溶剂、B溶剂里生成的沉淀的组成、晶型完全一样的时候才适用，因为此时溶质从固相到气相的自由能变化在比较的时候刚好抵消掉了。但如果沉淀晶格里带了溶剂，就会导致A、B溶剂里沉淀组成不一样，此时固相到气相的自由能变化抵消不掉。

cyr

wzkchem5 发表于 2023-10-26 15:40
那样需要想办法算沉淀的自由能，例如用支持周期性边界条件的程序做计算。
溶解自由能指的是溶质从气相到 ...

好的，谢谢老师！