Efficient Test Molecule Sampling in Molecular Simulation

Order and correlation contributions to the entropy of hydrophobic solvation

The entropy of hydrophobic solvation has been explained as the result of ordered solvation structures, of hydrogen bonds, of the small size of the water molecule, of dispersion forces, and of solvent density fluctuations. We report a new approach to the calculation of the entropy of hydrophobic solvation, along with tests of and comparisons to several other methods. The methods are assessed in the light of the available thermodynamic and spectroscopic information on the effects of temperature on hydrophobic solvation. Five model hydrophobes in SPC/E water give benchmark solvation entropies via Widom's test-particle insertion method, and other methods and models are tested against these particle-insertion results. Entropies associated with distributions of tetrahedral order, of electric field, and of solvent dipole orientations are examined. We find these contributions are small compared to the benchmark particle-insertion entropy. Competitive with or better than other theories in accuracy, but with no free parameters, is the new estimate of the entropy contributed by correlations between dipole moments. Dipole correlations account for most of the hydrophobic solvation entropy for all models studied and capture the distinctive temperature dependence seen in thermodynamic and spectroscopic experiments. Entropies based on pair and many-body correlations in number density approach the correct magnitudes but fail to describe temperature and size dependences, respectively. Hydrogen-bond definitions and free energies that best reproduce entropies from simulations are reported, but it is difficult to choose one hydrogen bond model that fits a variety of experiments. The use of information theory, scaled-particle theory, and related methods is discussed briefly. Our results provide a test of the Frank-Evans hypothesis that the negative solvation entropy is due to structured water near the solute, complement the spectroscopic detection of that solvation structure by identifying the structural feature responsible for the entropy change, and point to a possible explanation for the observed dependence on length scale. Our key results are that the hydrophobic effect, i.e. the signature, temperature-dependent, solvation entropy of nonpolar molecules in water, is largely due to a dispersion force arising from correlations between rotating permanent dipole moments, that the strength of this force depends on the Kirkwood g-factor, and that the strength of this force may be obtained exactly without simulation.

How Does Solute-Polarization Affect the Hydrophobic Hydration of Methane?

We simulate the solubility and solvation free energy of methane dissolved in water at infinite dilution. Molecular dynamics simulations of TIP4P-Ew model water are carried out at ambient pressure conditions over a large temperature interval, ranging from 250 K to 370 K. Solvation free energies are determined using the Widom particle insertion method. The fitted temperature dependent data is used to calculate solvation enthalpies, entropies, as well as the heat capacity of solvation. In particular we study the effect of polarizability of methane on those thermodynamic parameters. Solute polarization leads to a lowering of the solvation free energy at 298 K to 8.3 kJ mol−1, almost exactly matching the experimental value. A close inspection of the enthalpic and entropic contributions, however, reveals that this coincidence is a consequence of a compensation of enthalpic and entropic contributions, each of them deviating even larger from their respective experimental values. Surprisingly, the solute polarizability is apparently affecting the solvation entropy more strongly than the solvation enthalpy, leading to an about 5 J K−1 mol−1 smaller (less negative) solvation entropy compared to the non-polarizable model. The solute-water radial distribution functions of the polarizable particle reveals significant modifications, favoring small distances, as well as structural changes, very similar to those caused by a temperature increase. This is suggesting that the reduced negative solvation entropy of a polarizable methane particle is related to a more disordered, "high-temperature"-like hydration shell.

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: A detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.

Temperature dependence of the hydrophobic hydration and interaction of simple solutes: An examination of five popular water models

We examine five different popular rigid water models (SPC, SPCE, TIP3P, TIP4P, and TIP5P) using molecular dynamics simulations in order to investigate the hydrophobic hydration and interaction of apolar Lennard-Jones solutes as a function of temperature in the range between 275 and 375 K along the 0.1 MPa isobar. For all investigated models and state points we calculate the excess chemical potential for the noble gases and methane employing the Widom particle insertion technique. All water models exhibit too small hydration entropies, but show a clear hierarchy. TIP3P shows poorest agreement with experiment, whereas TIP5P is closest to the experimental data at lower temperatures and SPCE is closest at higher temperatures. As a first approximation, this behavior can be rationalized as a temperature shift with respect to the solvation behavior found in real water. A rescaling procedure inspired by the information theory model of Hummer et al. [Chem. Phys. 258, 349 (2000)] suggests that the different solubility curves for the different models and real water can be largely explained on the basis of the different density curves at constant pressure. In addition, the models that give a good representation of the water structure at ambient conditions (TIP5P, SPCE, and TIP4P) show considerably better agreement with the experimental data than the ones which exhibit less structured O-O correlation functions (SPC and TIP3P). In the second part of the paper we calculate the hydrophobic interaction between xenon particles directly from a series of 60 ns simulation runs. We find that the temperature dependence of the association is to a large extent related to the strength of the solvation entropy. Nevertheless, differences between the models seem to require a more detailed molecular picture. The TIP5P model shows by far the strongest temperature dependence. The suggested density rescaling is also applied to the chemical potential in the xenon-xenon contact-pair configuration, indicating the presence of a temperature where the hydrophobic interaction turns into purely repulsive. The predicted association for xenon in real water suggests the presence of a strong variation with temperature, comparable to the behavior found for TIP5P water. Comparing different water models and experimental data we conclude that a proper description of density effects is an important requirement for a water model to account correctly for the correct description of the hydrophobic effects. A water model exhibiting a density maximum at the correct temperature is desirable.