溶解自由能的温度依赖性如何计算

New-boy

各位老师好，看了Sobereva老师在谈谈隐式溶剂模型下溶解自由能和体系自由能的计算博文，隐式溶剂化模型下无法求解不同温度下的溶解自由能，请问有什么办法计算不同温度下的溶解自由能呢。
我在算不同温度下有机胺共轭酸的酸度系数，按照热力学循环框架计算的。
温度相关性只考虑了气相自由能，是采用Shermo软件扫描了包含振动分析的输出文件获得的。
溶解自由能是在M052X/6-31G*级别SMD隐式溶剂化模型下的自由能减去相同计算级别下的气相自由能获得的，不同温度下都采用了这一溶解自由能，即标况下的溶解自由能。
这样算出来的pKa值随温度的变化和实验值差距较大。目前分析是没有考虑溶解自由能的温度相关性。
请问使用uESE溶剂模型是否可以求解不同温度下的溶解自由能，或者有什么其他的方法可以计算。

sobereva

考虑温度计算溶解自由能的主流做法是跑分子动力学
uESE同样不行