Mode-coupling study on the dynamics of hydrophobic hydration II: Aqueous solutions of benzene and rare gases

The dynamic properties of both the solute and solvent of the aqueous solution of benzene, xenon and neon are calculated by the mode-coupling theory for molecular liquids based on the interaction-site model. The B-coefficients of the reorientational relaxation and the translational diffusion of the solvent are evaluated from their dependence on the concentration of the solute, and the reorientational relaxation time of water within the hydration shell is estimated based on the two-state model. The reorientational relaxation times of water in the bulk and within the hydration shell, that of solute, and the translational diffusion coefficients of solute and solvent, are calculated at 0-30 degrees C. The temperature dependence of these dynamic properties is in qualitative agreement with that of NMR experiment reported by Nakahara et al. (M. Nakahara, C. Wakai, Y. Yoshimoto and N. Matubayasi, J. Phys. Chem., 1996, 100, 1345-1349, ref. 36), although the agreement of the absolute values is not so good. The B-coefficients of the reorientational relaxation times for benzene, xenon and neon solution are correlated with the hydration number and the partial molar volume of the solute. The proportionality with the latter is better than that with the former. These results support the mechanism that the retardation of the mobility of water is caused by the cavity formation of the solute, as previously suggested by us (T. Yamaguchi, T. Matsuoka and S. Koda, J. Chem. Phys., 2004, 120, 7590-7601, ref. 34), rather than the conventional one that the rigid hydration structure formed around the hydrophobic solute reduces the mobility of water.     

The hydration of hydrophobic substances

A model of the hydration of hydrophobic substances in water is suggested. The models of fluctuation formation of empty cavities in water as a stage of hydration extensively used in the literature were shown to be at variance with experiment. The fundamental role played by the interphase boundary surface was emphasized. On this surface, the successive addition of water molecules with the formation of capsules around hydrophobic molecules occurred. The physical meaning of the Ostwald equation was revealed. This equation characterized the distribution of hydrophobic volatile substances between the gas and aqueous phases. The method of optical probes (hydrophobic aromatic molecules) was used to reveal the synergistic character of autocorrelation of dispersion interactions between water and hydrophobic substance molecules. This synergism was at variance with the Lennard-Jones potential. The synergism (superadditivity) of dispersion attraction forces, which strengthened their directional character, caused the self-organization and enhanced stability of hydration capsules with encapsulated hydrophobic molecules. Computer models were used to show that the spatially directional character of dispersion interactions necessary for the self-organization of hydrated aggregates could be simulated by the molecular mechanics method on the basis of orientational correlation of water molecules and hydrophobic substances in the starting system.     

Group contributions and hydrophobic hydration

The group contribution concept proposed by Savage and Wood is used for characterization of hydrophobic hydration. A parameter is proposed to measure hydrophobic hydration that has additive properties and can be applied to both nonelectrolytes and electrolytes. Theoretical and practical arguments are given for its use. The given results are in agreement with the published experimental data and theoretical results concerning hydrophobic hydration and interaction.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

Raman spectroscopic study of hydrophobic hydration of organic molecules in aqueous solution

Raman spectra of several organic molecules which form hydrate clathrates have been measured in aqueous solution, and the hydration structure around those molecules has been investigated from the analysis of Raman linewidth data.     

Structural contribution to the effect of hydrophobic hydration of noble gases

Within the concepts of structurally-thermodynamic characteristics of solvation and pseudo-chemical potential, the sample collection of the most authentic experimental data on solubility of gaseous He, Ne, Ar, Kr, Xe, and Rn in H₂O and D₂O is analysed at ≈0.1 MPa and T = 278–318 K. The conclusion is drawn that at deuteration of water molecules and also with increasing molar mass of noble gas, the relative contribution of effect of its hydrophobic hydration decreases. However in case of pass from lightweight noble gases (He, Ne, Ar) to heavy ones (Kr, Xe, Rn), structural transformations in their aqueous solutions become more expressed as a whole due to strengthening interaction between dissolved substance and solvent.     

Evaluation of the contribution to hydration of nonelectrolytes from the hydrophobic effect

A new method was suggested for estimating the hydrophobic effect of contributions to the Gibbs energies and enthalpies of hydration of hydrocarbons, inorganic gases and rare gases. In accordance with this method the hydrophobic effect contribution to the Gibbs energy was evaluated from the difference between the hydration Gibbs energy of a solute and the non hydrophobic contribution. To estimate the latter value, the known dependence connecting the Gibbs energies of solvation of a solute in a number of aprotic solvents to the Hildebrand solubility parameter for these solvents was used. The non hydrophobic contribution to the Gibbs energy of hydration was calculated for various solutes from such dependences extended to water as solvent. The Hildebrand solubility parameter for water used in the calculation was corrected for the effect of association through hydrogen bonding. This correction was made by subtraction of the water self-association enthalpy from the enthalpy of vaporization of water. The evaluated Gibbs energies of the hydrophobic effect are positive for saturated hydrocarbons, inorganic gases and rare gases and linearly depend on the solute molecular refraction. The hydrophobic contribution to the hydration enthalpies of the solutes was calculated in the same manner as was made to calculate the hydrophobic contribution to Gibbs energies of hydration. Enthalpies of the hydrophobic effect for the solutes under study are negative.     

The effect of pressure on the hydration structure around hydrophobic solute: a molecular dynamics simulation study

Molecular dynamics simulations are performed to study the effects of pressure on the hydrophobic interactions between neopentane molecules immersed in water. Simulations are carried out for five different pressure values ranging from 1 atm to 8000 atm. From potential of mean force calculations, we find that with enhancement of pressure, there is decrease in the well depth of contact minimum (CM) and the relative stability of solvent separated minimum over CM increases. Lower clustering of neopentane at high pressure is also observed in association constant and cluster-structure analysis. Selected site-site radial distribution functions suggest efficient packing of water molecules around neopentane molecules at elevated pressure. The orientational profile calculations of water molecules show that the orientation of water molecules in the vicinity of solute molecule is anisotropic and this distribution becomes flatter as we move away from the solute. Increasing pressure slightly changes the water distribution. Our hydrogen bond properties and dynamics calculations reveal pressure-induced formation of more and more number of water molecules with five and four hydrogen bond at the expense of breaking of two and three hydrogen bonded water molecules. We also find lowering of water-water continuous hydrogen bond lifetime on application of pressure. Implication of these results for relative dispersion of hydrophobic molecules at high pressure are discussed.     

Long-time correlations and hydrophobe-modified hydrogen-bonding dynamics in hydrophobic hydration

The physical mechanisms behind hydrophobic hydration have been debated for over 65 years. Spectroscopic techniques have the ability to probe the dynamics of water in increasing detail, but many fundamental issues remain controversial. We have performed systematic first-principles ab initio Car-Parrinello molecular dynamics simulations over a broad temperature range and provide a detailed microscopic view on the dynamics of hydration water around a hydrophobic molecule, tetramethylurea. Our simulations provide a unifying view and resolve some of the controversies concerning femtosecond-infrared, THz-GHz dielectric relaxation, and nuclear magnetic resonance experiments and classical molecular dynamics simulations. Our computational results are in good quantitative agreement with experiments, and we provide a physical picture of the long-debated "iceberg" model; we show that the slow, long-time component is present within the hydration shell and that molecular jumps and over-coordination play important roles. We show that the structure and dynamics of hydration water around an organic molecule are non-uniform.     

Global thermodynamics of hydrophobic cavitation, dewetting, and hydration

Pure water experimental and simulation results are combined to predict the thermodynamics of cavity formation, spanning atomic to macroscopic length scales, over the entire ambient liquid temperature range. The resulting cavity equation of state is used to quantify dewetting excess contributions to cavity formation thermodynamics and construct a thermodynamic perturbation theory of hydrophobic hydration. Predictions are compared with large cavity simulations and experimental rare-gas hydration thermodynamics data (for He, Ne, Ar, Kr, Xe, and Rn). Key findings include the strong temperature dependence of the critical length scale for hydrophobic dewetting and the evaluation of fundamental solute-solvent interaction contributions to rare-gas hydration chemical potentials.     

A study of the hydration of ribonuclease A using isothermal calorimetry

This work is part of a systematic study undertaken to find the excess thermodynamic functions of binary protein–water systems. Isothermal calorimetry and water sorption measurements were applied to characterize the hydration dependencies of the excess thermodynamic functions. The advantages of our methodology are (i) we are able to simultaneously determine the excess partial quantities of water and proteins; (ii) these thermodynamic quantities can be determined in the entire range of water content. Here, in particular, the excess partial enthalpies of water and bovine pancreatic ribonuclease A (RNase A) have been determined. The excess partial enthalpies for RNase A are compared with the published data for several unrelated globular proteins (lysozyme, chymotrypsinogen A, serum albumin, lactoglobulin). These biomacromolecules represent a series of proteins in which the hydrophobicity of proteins is gradually changed in a wide range. It was found that the excess partial quantities for the studied proteins are determined by the hydration of the hydrophilic and hydrophobic protein groups. The more hydrophilic a protein, the more significant a hydrophilic hydration contribution is and vice versa. RNase A is the most hydrophilic protein under the study. This protein has the most significant hydrophilic hydration contribution. Lactoglobulin is the most hydrophobic protein under the study. This protein has the most significant hydrophobic hydration contribution.     

Cooperativity of hydrophillic and hydrophobic effects in the aqueous hydration of polyfunctional solutes

Liquid state computer simulation is reported on a dilute aqueous solution of formamide at 25°C. Analysis of the results reveals essential characteristics of hydrophilic and hydrophobic hydration in the first shell of the hydration complex. A novel cooperativity between the hydrophilic hydration of the carbonyl group and the hydrophobic hydration of the methyne group is observed.     

Two-dimensional hydration shells of alkali metal ions at a hydrophobic surface

We study the hydration shell formation of alkali metal ions at a graphite surface. Two-dimensional shell structures are found in the initial stage of hydration, in contrast to the three-dimensional structures in bulk water and clusters. Comparison of vibrational spectra with experiments identifies the shell structures and the thermally induced transition from the first to the second shell. We also found intriguing competition between hydration and ion-surface interaction, leading to different solvation dynamics between K and Na. Implications of these results in ionic processes at interfaces are elaborated.     

A Monte Carlo study of liquid benzene

A Monte Carlo study for liquid benzene, where the pair potential is described as a sum of Lennard-Jones (12, 6) interactions acting between six equivalent centres of each molecule, has been carried out. The potential parameters have been chosen in such a way that the internal energy and virial pressure were in agreement with the experimental values. The liquid structure has been examined in detail by computing several distribution functions: the radial, the angular correlations and the radial-angular distribution between molecular planes.     

Microcalorimetric study on micellization of nonionic surfactants with a benzene ring or adamantane in their hydrophobic chains

The micellization of a novel family of nonionic surfactants poly(oxyethylene) glycol alkyl ethers has been studied by microcalorimetry. One of the surfactants has adamantane, and the other nonionic surfactants have a benzene ring in their hydrophobic chains, which moves from the terminal of the hydrophobic chain toward the headgroup. Moreover, the alkyl chain of the nonionic surfactants is straight or branched. Both the critical micelle concentration (cmc) and the thermodynamic parameters associated with the micelle formation have been obtained. The cmc decreases and the enthalpy of micelle formation (deltaH(mic)) becomes less positive gradually as the length of hydrophobic chain increases, whereas the values of cmc and deltaH(mic) tend to increase for the surfactants with a longer ethylene oxide chain. However, the deltaH(mic) value of the surfactant with seven carbon atoms in a hydrophobic chain is more positive than that of the surfactant with six carbon atoms in a hydrophobic chain. Comparing with the nonionic surfactant with a methylene hydrophobic chain, the surfactants with benzene rings and adamantane groups have larger cmc values and the cmc values increase with the size of the groups. Furthermore, moving the phenyl group from the terminal of the hydrophobic chain to the neighbor of the hydrophilic headgroup leads to the decreased cmc. Both the variation of hydrophobic interaction from the movement of phenyl group and pi-pi interaction among adjacent phenyl groups affect deltaH(mic) values.