Search

A goal in computational chemistry is computing hydration free energies of nonpolar and charged solutes accurately, but with much greater computational speeds than in today's explicit-water simulations. Here, we take one step in that direction: a simple model of solvating waters that is analytical and thus essentially instantaneous to compute. Each water molecule is a 2-dimensional dipolar hydrogen-bonding disk that interacts around small circular solutes with different nonpolar and charge interactions. The model gives good qualitative agreement with experiments. As a function of the solute radius, it gives the solvation free energy, enthalpy and entropy as a function of temperature for the inert gas series Ne, Ar, Kr, and Xe. For anions and cations, it captures relatively well the trends versus ion radius. This approach should be readily generalizable to three dimensions.

Water is an unusual liquid in its solvation properties. Here, we model the process of transferring a nonpolar solute into water. Our goal was to capture the physical balance between water's hydrogen bonding and van der Waals interactions in a model that is simple enough to be nearly analytical and not heavily computational. We develop a 2-dimensional Mercedes-Benz-like model of water with which we compute the free energy, enthalpy, entropy, and the heat capacity of transfer as a function of temperature, pressure, and solute size. As validation, we find that this model gives the same trends as Monte Carlo simulations of the underlying 2D model and gives qualitative agreement with experiments. The advantages of this model are that it gives simple insights and that computational time is negligible. It may provide a useful starting point for developing more efficient and more realistic 3D models of aqueous solvation.

We study water that is confined within small geometric spaces. We use the Mercedes-Benz (MB) model of water, in NVT and muVT Monte Carlo computer simulations. For MB water molecules between two planes separated by a distance d, we explore the structures, hydrogen bond networks, and thermodynamics as a function of d, temperature T, and water chemical potential mu. We find that squeezing the planes close enough together leads to a vaporization of waters out of the cavity. This vaporization transition has a corresponding peak in the heat capacity of the water. We also find that, in small pores, hydrogen bonding is not isotropic but, rather, it preferentially forms chains along the axis of the cavity. This may be relevant for fast proton transport in pores. Our simulations show oscillations in the forces between the inert plates, due to water structure, even for plate separations of 5-10 water diameters, consistent with experiments by Israelachvili et al. [Nature 1983, 306, 249]. Finally, we find that confinement affects water's heat capacity, consistent with recent experiments of Tombari et al. on Vycor nanopores [J. Chem. Phys. 2005, 122, 104712].