求指导显式溶剂化自由能计算方法

aipacino

各位老师好，我目前得到了一个机器学习势能模型，想计算Li离子在各种溶剂（以EC为例）中的溶剂化自由能，卢老师在这篇帖子（谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 - 量子化学 (Quantum Chemistry) - 计算化学公社）中的解释我仍不完全理解，希望得到进一步指导：

1. 我想计算Li离子在EC中的溶剂化自由能的话，是不是先确认Li和EC的配位构型（通过MD得到径向分布函数），然后根据这个配位结构搭建团簇结构做构型搜索，找到最稳定的构型，最后根据Li离子在气相的自由能，EC在EC环境下的自由能，与EC+Li在EC环境下的自由能作差得到？请问这个自由能只能通过Gaussian等软件的Gibbs Free Energy值得到吗？我的模型只能计算势能是不是就用不上？另外，需要考虑第二溶剂化层吗？


2.我还想绘制自由能-(Li-O)距离曲线，由于机器学习势模型固定，不能额外添加势约束，可不可以固定一系列的Li-O距离，然后做限制性分子动力学模拟，最后积分得到自由能变？

问题较多，感谢老师们解惑