高斯计算过渡金属离子溶解自由能该使用什么溶剂模型？

取不起名字

各位大神好，我目前已经在IEFPCM溶剂模型下对钼离子配合水分子和氢氧根离子的结构进行了优化，在计算溶解自由能时遇到了问题如下：
1. 由于钼离子已经与水分子进行了配位，那么是否意味着我在隐式溶剂模型下计算溶解自由能就可以了?还是说我需要在杂化溶剂模型下计算？
2. 如果在隐式溶剂模型下计算就可以，我计算的体系涵盖了电荷数-2~+2的配合物结构，那使用uESE更好还是SMD？如果需要用杂化溶剂模型，那么可以提供一点文献参考吗？
希望各位大神可以指点迷津，谢谢大家！

sobereva

如果算的就是配合物的溶解自由能，隐式溶剂模型一般就够了。建议uESE

ljh123

uESE不是很稳定。建议看一下你的体系算这些量需不需要对比。金属配合物体系可能涉及到好几种带电状态，我试过有一些体系带个+1、+2之类的算出来比SMD准，但是中性下又远不如SMD。

取不起名字

ljh123 发表于 2025-2-25 16:15
uESE不是很稳定。建议看一下你的体系算这些量需不需要对比。金属配合物体系可能涉及到好几种带电状态，我试 ...

谢谢您的解答。我想算质子化过程的吉布斯自由能变，算的体系带的电荷数是从-2~+2的，包含了中性的配合物，但理论方法和基组需要统一的。

取不起名字

sobereva 发表于 2025-2-25 11:58
如果算的就是配合物的溶解自由能，隐式溶剂模型一般就够了。建议uESE

谢谢老师的解答。我还有一个问题：我的目的是计算溶质在溶剂中的自由能，需要气相自由能+溶解自由能，您提到这两个部分得分别尽量算准确。那么我可以在原先用IEFPCM溶剂模型下优化的结构算高精度单点能，在B3LYP/def2TZVP下优化的结构算溶解自由能吗？还是说必须统一用一个结构，我看了您的博文后感觉用B3LYP/def2TZVP对过渡金属配合物的优化不太合适，因此对于气相自由能的计算有些迷惑。