How hydrophobic hydration responds to solute size and attractions: Theory and simulations

We focus on the hydration of a methane and spherical single and multisite C60 and C180 solutes over a range of solute-water attractions to quantify the vicinal water structure and their hydration thermodynamics using extensive molecular dynamics simulations and theory. We show that water structure near larger solutes is more sensitive to solute-water attractions compared to that near smaller ones. To understand the sensitivity, we separate the solute-water potential of mean force into a direct solute-water interaction and an indirect or solvent contribution [omega(r)]. In the absence of omega(r), water density in the solute vicinity would increase exponentially with solute-water interactions. Instead, omega(r) becomes increasingly repulsive with strengthening of solute-water attractions thereby opposing those direct interactions. We term this phenomenon "competitive expulsion," which characterizes the repulsion of a test water molecule by the hydration shell solvent waters. We develop a physically motivated theoretical approach to predict changes in omega(r) with attractions. We call this approach the modified-EXP (M-EXP) approximation owing to the similarity of ideas and especially our final expression with that of the EXP approximation of Chandler and Andersen [J. Chem. Phys. 57, 1930 (1972)]. Solute-water radial distribution functions and chemical potentials calculated using the M-EXP approach are in good agreement with simulation data. These calculations highlight the sensitivity of hydration structure and thermodynamics of bucky ball like solutes to solute-water interactions. We find that excess chemical potentials of bucky balls with standard alkane-like carbon-water interactions parameters are negative, suggesting the need for a careful calibration of those parameters for predictions of solubility, wetting, and water-mediated interactions using molecular simulations.