gaussian计算Li+的溶剂化能的几何优化和计算公式的问题

Dreamer1

老师好，我在用gromacs跑完分子动力学后，搜索得到最有可能存在的溶剂化结构 Li-2TFSI-2DMF。Li是阳离子，TFSI是阴离子，DMF是溶剂分子。现在我想通过Gaussian计算这个溶剂化结构的溶剂化能。有以下三个问题：
(1)计算公式：E(溶剂化能）=E(Li-2TFSI-2DMF)-2 E(TFSI)-2 E(DMF)-E(Li)。这样的公式正确吗？
(2)看过谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 - 思想家公社的门口：量子化学·分子模拟·二次元这条帖子后，想问，(1)中的每一个E 是不是得参照帖子中的“4 溶剂环境下溶质自由能的计算实例”的部分计算出来，然后再用（1）公式计算得到？
(3)昨天我尝试用“#P B3LYP/6-311G** em=GD3BJ opt freq geom=connectivity”  对 “Li-2TFSI-2DMF”整个结构进行几何优化，发现优化后的成键都变了，甚至原来的分子的键都断开了，这样的优化对吗？然后想请教老师，几何优化该用什么基组和泛函合适呢？是不是还得加上SMD隐式溶剂模型进行优化？

下面附上了溶剂化结构的图片。

sobereva

不正确。这把TFSI、DMF之间的相互作用能也算进去了。而且也不是体现把Li+插入到原本体系中的自由能能量变化
没有很简单的办法算准Li+在这种环境中的溶解自由能