The importance of polarizability in the modeling of solubility: Quantifying the effect of solute polarizability on the solubility of small nonpolar solutes in popular models of water

Crustwater: Modeling Hydrophobic Solvation

We describe Crustwater, a statistical mechanical model of nonpolar solvation in water. It treats bulk water using the Cage Water model and introduces a crust, i.e., a solvation shell of coordinated partially structured waters. Crustwater is analytical and fast to compute. We compute here solvation vs temperature over the liquid range, and vs pressure and solute size. Its thermal predictions are as accurate as much more costly explicit models such as TIP4P/2005. This modeling gives new insights into the hydrophobic effect: (1) that oil-water insolubility in cold water is due to solute-water (SW) translational entropy and not water-water (WW) orientations, even while hot water is dominated by WW cage breaking, and (2) that a size transition at the Angstrom scale, not the nanometer scale, takes place as previously predicted.

Temperature, Pressure, and Concentration Derivatives of Nonpolar Gas Hydration: Impact on the Heat Capacity, Temperature of Maximum Density, and Speed of Sound of Aqueous Mixtures

The hydrophobic effect is an umbrella term encompassing a number of solvation phenomena associated with solutions of nonpolar species in water, including the following: a meager solubility opposed by entropy at room temperature; large positive hydration heat capacities; positive shifts in the temperature of maximum density of aqueous mixtures; increases in the speed of sound of dilute aqueous mixtures; and negative volumes of association between interacting solutes. Here we present a molecular simulation study of nonpolar gas hydration over the temperature range 273.15-373.15 K and a pressure range -500 to 1000 bar to investigate the interrelationships between distinct hydrophobic phenomena. We develop a new free energy correlation for the solute chemical potentials founded on the Tait equation description of the equation-of-state of liquid water. This analytical correlation is shown to provide a quantitatively accurate description of nonpolar gas hydration over the entire range of thermodynamic state points simulated, with an error of ∼0.02 k_BT or lower in the fitted chemical potentials. Our simulations and the correlation accurately reproduce many of the available experimental results for the hydration of the solutes examined here. Moreover, the correlation reproduces the characteristic entropies of hydration, temperature dependence of the hydration heat capacity, perturbations in the temperature of maximum density, and changes in the speed of sound. While negative volumes of association result from pairwise interactions in solution, beyond the limits of our simulations performed at infinite dilution, we discuss how our correlation could be supplemented with second virial coefficient information to expand to finite concentrations. In total, this work demonstrates that many distinct phenomena associated with the hydrophobic effect can be captured within a single thermodynamically consistent correlation for solute hydration free energies.

Three-Dimensional RISM Integral Equation Theory for Polarizable Solute Models

Modeling solute polarizability is a key ingredient for improving the description of solvation phenomena. In recent years, polarizable molecular mechanics force fields have emerged that circumvent the limitations of classical fixed charge force fields by the ability to adapt their electrostatic potential distribution to a polarizing environment. Solvation phenomena are characterized by the solute's excess chemical potential, which can be computed by expensive fully atomistic free energy simulations. The alternative is to employ an implicit solvent model, which poses a challenge to the formulation of the solute-solvent interaction term within a polarizable framework. Here, we adapt the three-dimensional reference interaction site model (3D RISM) integral equation theory as a solvent model, which analytically yields the chemical potential, to the polarizable AMOEBA force field using an embedding cluster (EC-RISM) strategy. The methodology is analogous to our earlier approach to the coupling of a quantum-chemical solute description with a classical 3D RISM solvent. We describe the conceptual physical and algorithmic basis as well as the performance for several benchmark cases as a proof of principle. The results consistently show reasonable agreement between AMOEBA and quantum-chemical free energies in solution in general and allow for separate assessment of energetic and solvation-related contributions. We find that, depending on the parametrization, AMOEBA reproduces the chemical potential in better agreement with reference quantum-chemical calculations than the intramolecular energies, which suggests possible routes toward systematic improvement of polarizable force fields.

Development of the trappe force field for ammonia

The transferable potentials for phase equilibria (TraPPE) force field is extended through the development of a non-polarizable five-site ammonia model. In this model, the electrostatic interactions are represented by three positive partial charges placed at the hydrogen position and a compensating partial charge placed on an M site that is located on the C3 molecular axis and displaced from the nitrogen atom toward the hydrogen atoms. The repulsive and dispersive interactions are represented by placing a single Lennard–Jones site at the position of the nitrogen atom. Starting from the five-site model by Impey and Klein (Chem. Phys. Lett. 1984, 104, 579), this work optimizes the Lennard–Jones parameters and the magnitude of the partial charges for three values of the M site displacement. This parameterization is done by fitting to the vapor–liquid coexistence curve of neat ammonia. The accuracy of the three resulting models (differing in the displacement of the M site) is assessed through computation of the binary vapor–liquid equilibria with methane, the structure and the dielectric constant of liquid ammonia. The five-site model with an intermediate displacement of 0.08 Å for the M site yields a much better value for the dielectric constant, whereas differences in the other properties are quite small.

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.