求助：计算锕系配合物吉布斯自由能变中溶剂模型选择的问题

wangzh

老师们好，对于含有锕系离子配合物形成过程的吉布斯自由能变的计算，所用的泛函和基组一样的情况下，以此为例Th(H2O)4+aq + Lorg + 4NO3-aq → ThL(NO3)4org + H2Oaq，发现用IEFPCM模型优化+频率，SMD模型用更高精度单点计算的吉布斯自由能变与IEFPCM模型优化+频率，在IEFPCM模型用更高精度单点计算得到的吉布斯自由能变差距很大，前者算出来是正值，后者算出来是负值。
1.请问老师哪个模型合理些？或者对于含有锕系的体系，比较推荐那种溶剂模型？
2.为什么不同的溶剂模型，在结构，基组，泛函相同的情况下会得到如此大的误差？
3.这个反应过程中，发现两个溶剂模型对Th(H2O)4+aq的吉布斯自由能影响比较大，这是为什么？
感谢老师们的解惑！

zrm1234

吉布斯自由能变涉及的修正项是不是应该与优化+频率的计算级别严格一致？SDD是指SMD吗？如果你在SMD模型下用高精度单点并修正获得吉布斯自由能，你是否应该在优化+频率的时候也用SMD模型来获得低精度的修正项？

wangzh

zrm1234 发表于 2024-10-27 22:48
吉布斯自由能变涉及的修正项是不是应该与优化+频率的计算级别严格一致？SDD是指SMD吗？如果你在SMD模型下用 ...

是SMD模型，我已经修改过来了，不好意思。优化+频率不一定要用SMD模型，我的金属配合物结构在SMD模型下优化不出来。