溶液反应自由能变计算中，各物种是否需要使用同一方法处理？

Ender_Wang

各位老师好！最近我在计算如下反应的自由能变：
Y(aq) + B(OH)4-(aq) → YB-(aq) + nH2O
Y(aq) + B(OH)3(aq) → YB(aq) + nH2O
其中Y为有机分子，YB、YB-为分子与硼酸、硼酸根离子的配合物。根据帖子《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》（sobereva.com/327）和《谈谈弥散函数和“月份”基组》（sobereva.com/119）的内容，我对阴离子YB-(aq)的结构优化后使用了M06-2X/ma-def2-TZVPP基组，再放到Shermo计算G(YB-, g)，然后用M05-2X/6-31G*计算了溶解自由能，算出了它的液相自由能G(YB-, aq)。但是，对于中性分子Y(aq)、YB(aq)和其他小分子，是否都需要加上弥散函数进行计算，还是应使用不加弥散函数的版本？即：计算反应自由能是每种反应物/生成物分别算准能量，还是必须统一用相同的方法？当对比不同的Y对硼酸的配位能力时，不同的Y之间是否也需要相同的方法？

RichardYin

我也在算类似的配合物体系，请问H2O也需要按其它分子的方法计算G(H2O,g)+校正项+溶解自由能吗？

sobereva

整个过程中，相同的原子（不是指元素）必须用相同的基组，否则没法求差