Cluster expansion of the solvation free energy difference: Systematic improvements in the solvation of single ions.

The cluster expansion method has been used in the imperfect gas theory for several decades. This paper proposes a cluster expansion of the solvation free energy difference. This difference, which results from a change in the solute-solvent potential energy, can be written as the logarithm of a finite series. Similar to the Mayer function, the terms in the series are related to configurational integrals, which makes the integrand relevant only for configurations of the solvent molecules close to the solute. In addition, the terms involve interaction of solute with one, two, and so on solvent molecules. The approach could be used for hybrid quantum mechanical and molecular mechanics methods or mixed cluster-continuum approximation. A simple form of the theory was applied for prediction of pKa in methanol; the results indicated that three explicit methanol molecules and the dielectric continuum lead to a root of mean squared error (RMSE) of only 1.3 pKa units, whereas the pure continuum solvation model based on density method leads to a RMSE of 6.6 pKa units.