Hydration free energy of a Model Lennard-Jones solute particle: microscopic Monte Carlo simulation studies, and interpretation based on mesoscopic models

In this study, the hydration of a model Lennard-Jones solute particle and the analytical approximations of the free energy of hydration as functions of solute microscopic parameters are analyzed. The control parameters of the solute particle are the charge, the Lennard-Jones diameter, and also the potential well depth. The obtained multivariate free energy functions of hydration were parametrized based on Metropolis Monte Carlo simulations in the extended NpT ensemble, and interpreted based on mesoscopic solvation models proposed by Gallicchio and Levy J. Comput. Chem. 25, 479 (2004)], and Wagoner and Baker Proc. Natl. Acad. Sci. U.S.A. 103, 8331 (2006)]. Regarding the charge and the solute diameter, the dependence of the free energy on these parameters is in qualitative agreement with former studies. The role of the third parameter, the potential well depth not previously considered, appeared to be significant for sufficiently precise bivariate solvation free energy fits. The free energy fits for cations and neutral solute particles were merged, resulting in a compact manifold of the free energy of solvation. The free energy of hydration for anions forms two separate manifolds, which most likely results from an abrupt change of the coordination number when changing the size of the anion particle.