Atomic-scale analysis of the solvation thermodynamics of hydrophobic hydration

Molecular dynamics simulations are used to model the transfer thermodynamics of krypton from the gas phase into water. Extra long, nanosecond simulations are required to reduce the statistical uncertainty of the calculated "solvation" enthalpy to an acceptable level. Thermodynamic integration is used to calculate the "solvation" free energy, which together with the enthalpy is used to calculate the "solvation" entropy. A comparison series of simulations are conducted using a single Lennard-Jones sphere model of water to identify the contribution of hydrogen bonding to the thermodynamic quantities. In contrast to the classical "iceberg" model of hydrophobic hydration, the favorable enthalpy change for the transfer process at room temperature is found to be due primarily to the strong van der Waals interaction between the solute and solvent. Although some stabilization of hydrogen bonding does occur in the solvation shell, this is overshadowed by a destabilization due to packing constraints. Similarly, whereas some of the unfavorable change in entropy is attributed to the reduced rotational motion of the solvation shell waters, the major component is due to a decrease in the number of positional arrangements associated with the translational motions.