Abstract

Quasichemical theory (QCT) provides a framework that can be used to partition the influence of the solvent surrounding an ion into near and distant contributions. Within QCT, the solvation properties of the ion are expressed as a sum of configurational integrals comprising only the ion and a small number of solvent molecules. QCT adopts a particularly simple form if it is assumed that the clusters undergo only small thermal fluctuations around a welldefined energy minimum and are affected exclusively in a meanfield sense by the surrounding bulk solvent. The fluctuations can then be integrated out via a simple vibrational analysis, leading to a closedform expression for the solvation free energy of the ion. This constitutes the primitive form of quasichemical theory (pQCT), which is an approximate mathematical formulation aimed at reproducing the results from the full manybody configurational averages of statistical mechanics. While the results from pQCT from previous applications are reasonable, the accuracy of the approach has not been fully characterized and its range of validity remains unclear. Here, a direct test of pQCT for a set of ion models is carried out by comparing with the results of free energy simulations with explicit solvent. The influence of the distant surrounding bulk on the cluster comprising the ion and the nearest solvent molecule is treated both with a continuum dielectric approximation and with free energy perturbation molecular dynamics simulations with explicit solvent. The analysis shows that pQCT can provide an accurate framework in the case of a small cation such as Li+. However, the approximation encounters increasing difficulties when applied to larger cations such as Na+, and particularly for K+. This suggests that results from pQCT should be interpreted with caution when comparing ions of different sizes.