Electron binding energies of aqueous alkali and halide ions: EUV photoelectron spectroscopy of liquid solutions and combined ab initio and molecular dynamics calculations

Photoelectron spectroscopy combined with the liquid microjet technique enables the direct probing of the electronic structure of aqueous solutions. We report measured and calculated lowest vertical electron binding energies of aqueous alkali cations and halide anions. In some cases, ejection from deeper electronic levels of the solute could be observed. Electron binding energies of a given aqueous ion are found to be independent of the counterion and the salt concentration. The experimental results are complemented by ab initio calculations, at the MP2 and CCSD(T) level, of the ionization energies of these prototype ions in the aqueous phase. The solvent effect was accounted for in the electronic structure calculations in two ways. An explicit inclusion of discrete water molecules using a set of snapshots from an equilibrium classical molecular dynamics simulations and a fractional charge representation of solvent molecules give good results for halide ions. The electron binding energies of alkali cations computed with this approach tend to be overestimated. On the other hand, the polarizable continuum model, which strictly provides adiabatic binding energies, performs well for the alkali cations but fails for the halides. Photon energies in the experiment were in the EUV region (typically 100 eV) for which the technique is probing the top layers of the liquid sample. Hence, the reported energies of aqueous ions are closely connected with both structures and chemical reactivity at the liquid interface, for example, in atmospheric aerosol particles, as well as fundamental bulk solvation properties.     

Ions in water: the microscopic structure of concentrated hydroxide solutions

Neutron-diffraction data on aqueous solutions of hydroxides, at solute concentrations ranging from 1 solute per 12 water molecules to 1 solute per 3 water molecules, are analyzed by means of a Monte Carlo simulation (empirical potential structure refinement), in order to determine the hydration shell of the OH- in the presence of the smaller alkali metal ions. It is demonstrated that the symmetry argument between H+ and OH- cannot be used, at least in the liquid phase at such high concentrations, for determining the hydroxide hydration shell. Water molecules in the hydration shell of K+ orient their dipole moment at about 45 degrees from the K+-water oxygen director, instead of radially as in the case of the Li+ and Na+ hydration shells. The K+-water oxygen radial distribution function shows a shallower first minimum compared to the other cation-water oxygen functions. The influence of the solutes on the water-water radial distribution functions is shown to have an effect on the water structure equivalent to an increase in the pressure of the water, depending on both ion concentration and ionic radius. The changes of the water structure in the presence of charged solutes and the differences among the hydration shells of the different cations are used to present a qualitative explanation of the observed cation mobility.     

Cation electric field is related to hydration energy

Electronic structure calculations on ions that use dielectric continuum theory to mimic solvent around the bare ionic solutes are often prone to make large errors in the hydration energies. It is found for cations that much of the error can be accounted for by a simple linear correlation with the maximum value found anywhere on the dielectric cavity surface of the solute potential or, even better, the outgoing normal electric field, thus mirroring analogous results previously obtained for anions. This correlation allows for significantly improved estimates of cation hydration energies while still retaining cavities of physically reasonable size in determination of the bulk dielectric contributions.     

Hydration forces underlie the exclusion of salts and of neutral polar solutes from hydroxypropylcellulose

The distance dependence for the preferential exclusion of several salts and neutral solutes from hydroxypropyl cellulose (HPC) has been measured via the effect of these small molecules on the thermodynamic forces between HPC polymers in ordered arrays. The concentration of salts and neutral solutes decreases exponentially as the spacing between apposing nonpolar HPC surfaces decreases. For all solutes, the spatial decay lengths of this exclusion are remarkably similar to those observed between many macromolecules at close spacings where intermolecular forces have been ascribed to the energetics of water structuring. Exclusion magnitudes depend strongly on the nature and size of the particular salt or solute; for the three potassium salts studied, exclusion follows the anionic Hofmeister series. The change in the number of excess waters associated with HPC polymers is independent of solute concentration suggesting that the dominating interactions are between solutes and the hydrated polymer. These findings further confirm the importance of solvation interactions and reveal an unexpected unity of Hofmeister effects, preferential hydration, and hydration forces.     

Effect of ionic solutes on the hydrogen bond network dynamics of water: power spectral analysis of aqueous NaCl solutions

To understand the modifications of the hydrogen bond network of water by ionic solutes, power spectra as well as static distributions of the potential energies of tagged solvent molecules and solute ions have been computed from molecular dynamics simulations of aqueous NaCl solutions. The key power spectral features of interest are the presence of high-frequency peaks due to localized vibrational modes, the existence of a multiple time scale or 1/falpha frequency regime characteristic of networked liquids, and the frequency of crossover from 1/falpha type behavior to white noise. Hydrophilic solutes, such as the sodium cation and the chloride anion, are shown to mirror the multiple time scale behavior of the hydrogen bond network fluctuations, unlike hydrophobic solutes which display essentially white noise spectra. While the power spectra associated with tagged H2O molecules are not very sensitive to concentration in the intermediate frequency 1/falpha regime, the crossover to white noise is shifted to lower frequencies on going from pure solvent to aqueous alkali halide solutions. This suggests that new and relatively slow time scales enter the picture, possibly associated with processes such as migration of water molecules from the hydration shell to the bulk or conversion of contact ion pairs into solvent-separated ion pairs which translate into variations in equilibrium transport properties of salt solutions with concentration. For anions, cations, and solvent molecules, the trends in the alpha exponents of the multiple time scale region and the self-diffusivities are found to be strongly correlated.     

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.     

A new angle on heat capacity changes in hydrophobic solvation

The differential solubility of polar and apolar groups in water is important for the self-assembly of globular proteins, lipid membranes, nucleic acids, and other specific biological structures through hydrophobic and hydrophilic effects. The increase in water's heat capacity upon hydration of apolar compounds is one signature of the hydrophobic effect and differentiates it from the hydration of polar compounds, which cause a decrease in heat capacity. Water structuring around apolar and polar groups is an important factor in their differential solubility and heat capacity effects. Here, it is shown that joint radial/angular distribution functions of water obtained from simulations reveal quite different hydration structures around polar and apolar groups: polar and apolar groups have a deficit or excess, respectively, of "low angle hydrogen bonds". Low angle hydrogen bonds have a larger energy fluctuation than high angle bonds, and analysis of these differences provides a physical reason for the opposite changes in heat capacity and new insight into water structure around solutes and the hydrophobic effect.     

Structure-property relationships in push-pull amino/cyanovinyl end-capped oligothiophenes: quantum chemical and experimental studies

A series of push-pull chromophores built around thiophene-based pi-conjugating spacers and bearing various types of amino donors and cyanovinyl acceptors have been analyzed by means of UV-vis-NIR, IR, and Raman spectroscopic measurements in the solid state as well as in solution. The intramolecular charge transfer (ICT) of these pi-conjugated systems has also been tested by analyzing the ability of the solute molecules to undergo shifts in their fluorescence emission maxima with increasing solvent polarity. These push-pull oligomers also display an attractive electrochemical behavior since they generate stable species both upon oxidation and reduction. Oxidation mainly involves changes in the electron-rich aminooligothienyl half-part of the molecule and leads to the formation of stable cations. On the other hand, reduction to radical anions and dianions is mainly cyanovinyl-centered but also affects the pi-conjugated electron relay. Density functional theory (DFT) calculations have been carried out to help the assignment of the most relevant electronic and vibrational features and to derive useful information about the molecular structure of these NLO-phores.     

Towards an understanding of the molecular basis of hydrophobicity

Summary The modern view is stressed that the structuring of water around nonpolar solutes, a process called hydrophobic hydration, actually favors the solubility of nonpolar solutes in water, its associated positive free energy of transfer arising from the enthalpic input required to create a cavity in water to accommodate the solute. The results of a series of molecular dynamics simulations of methane in SPC/E water at different temperatures are reported. These results show the existence of a larger fraction of broken hydrogen bonds in the hydration-shell water of the nonpolar solutes with respect to the bulk water, the difference increasing with a rise in temperature. This supports Muller's modified hydration-shell hydrogen-bond model predictions, where hydration-shell water molecules have lower free energies of hydrogen-bond breaking than those in the bulk.This paper is based on a presentation given at the 14th Molecular Graphics and Modelling Society Conference, held in Cairns, Australia, August 27 September 1, 1995.     

Effect of solute size and solute-water attractive interactions on hydration water structure around hydrophobic solutes.

Using Monte Carlo simulations, we investigated the influence of solute size and solute-water attractive interactions on hydration water structure around spherical clusters of 1, 13, 57, 135, and 305 hexagonally close-packed methanes and the single hard-sphere (HS) solute analogues of these clusters. We obtain quantitative results on the density of water molecules in contact with the HS solutes as a function of solute size for HS radii between 3.25 and 16.45 A. Analysis of these results based on scaled-particle theory yields a hydration free energy/surface area coefficient equal to 139 cal/(mol A2), independent of solute size, when this coefficient is defined with respect to the van der Waals surface of the solute. The same coefficient defined with respect to the solvent-accessible surface decreases with decreasing solute size for HS radii less than approximately 10 A. We also find that solute-water attractive interactions play an important role in the hydration of the methane clusters. Water densities in the first hydration shell of the three largest clusters are greater than bulk water density and are insensitive to the cluster size. In contrast, contact water densities for the HS analogues of these clusters decrease with solute size, falling below the bulk density of water for the two largest solutes. Thus, the large HS solutes dewet, while methane clusters of the same size do not.     

The Hydrophobic Effects: Our Current Understanding

Hydrophobic interactions are involved in and believed to be the fundamental driving force of many chemical and biological phenomena in aqueous environments. This review focuses on our current understanding on hydrophobic effects. As a solute is embedded into water, the interface appears between solute and water, which mainly affects the structure of interfacial water (the topmost water layer at the solute/water interface). From our recent structural studies on water and air-water interface, hydration free energy is derived and utilized to investigate the origin of hydrophobic interactions. It is found that hydration free energy depends on the size of solute. With increasing the solute size, it is reasonably divided into initial and hydrophobic solvation processes, and various dissolved behaviors of the solutes are expected in different solvation processes, such as dispersed and accumulated distributions in solutions. Regarding the origin of hydrophobic effects, it is ascribed to the structural competition between the hydrogen bondings of interfacial and bulk water. This can be applied to understand the characteristics of hydrophobic interactions, such as the dependence of hydrophobic interactions on solute size (or concentrations), the directional natures of hydrophobic interactions, and temperature effects on hydrophobic interactions.     

Structural and dynamic properties of concentrated alkali halide solutions: a molecular dynamics simulation study

The physicochemical properties of alkali halide solutions have long been attributed to the collective interactions between ions and water molecules in the solution, yet the structure of water in these systems and its effect on the equilibrium and dynamic properties of these systems are not clearly understood. Here, we present a systematic view of water structure in concentrated alkali halide solutions using molecular dynamics simulations. The results of the simulations show that the size of univalent ions in the solution has a significant effect on the dynamics of ions and other transport properties such as the viscosity that are correlated with the structural properties of water in aqueous ionic solution. Small cations (e.g., Li+) form electrostatically stabilized hydrophilic hydration shells that are different from the hydration shells of large ions (e.g., Cs+) which behave more like neutral hydrophobic particles, encapsulated by hydrogen-bonded hydration cages. The properties of solutions with different types of ion solvation change in different ways as the ion concentration increases. Examples of this are the diffusion coefficients of the ions and the viscosities of solutions. In this paper we use molecular dynamics (MD) simulations to study the changes in the equilibrium and transport properties of LiCl, RbCl, and CsI solutions at concentrations from 0.22 to 3.97 M.     

Excess chemical potential of small solutes across water--membrane and water--hexane interfaces

The excess chemical potentials of five small, structurally related solutes, CH4, CH3F, CH2F2, CHF3, and CF4, across the water-glycerol 1-monooleate bilayer and water-hexane interfaces were calculated at 300, 310, and 340 K using the particle insertion method. The excess chemical potentials of nonpolar molecules (CH4 and CF4) decrease monotonically or nearly monotonically from water to a nonpolar phase. In contrast, for molecules that possess permanent dipole moments (CH3F, CH2F, and CHF3), the excess chemical potentials exhibit an interfacial minimum that arises from superposition of two monotonically and oppositely changing contributions: electrostatic and nonelectrostatic. The nonelectrostatic term, dominated by the reversible work of creating a cavity that accommodates the solute, decreases, whereas the electrostatic term increases across the interface from water to the membrane interior. In water, the dependence of this term on the dipole moment is accurately described by second order perturbation theory. To achieve the same accuracy at the interface, third order terms must also be included. In the interfacial region, the molecular structure of the solvent influences both the excess chemical potential and solute orientations. The excess chemical potential across the interface increases with temperature, but this effect is rather small. Our analysis indicates that a broad range of small, moderately polar molecules should be surface active at the water-membrane and water-oil interfaces. The biological and medical significance of this result, especially in relation to the mechanism of anesthetic action, is discussed.     

Solute size effects on the solvation structure and diffusion of ions in liquid methanol under normal and cold conditions

We have performed a series of molecular dynamics simulations of alkali metal (Li+, Na+, K+, Rb+, and Cs+) and halide (F-, Cl-, Br-, and I-) ions in liquid methanol at two different temperatures to investigate the effects of ion size on the hydration structure and diffusion of ions in methanol under normal and cold conditions. Simulations are also carried out for some of the larger cations such as I+, (CH3)4N+, and (C2H5)4N+ and also neutral alkali metal atoms in methanol at both temperatures. With the increase of ion size, the diffusion coefficients of both positive and negative ions are found to show anomalous behavior. For cations, it is found that the maximum of the diffusion coefficient versus ion size curve occurs at the rather large cation of (CH3)4N+ unlike in water where the maximum occurs at the relatively smaller ion of Rb+. For halide ions, the anomalous behavior, i.e., the increase of diffusion with ion size, continues up to iodide ion and no maximum is observed. These results are in good agreement with experimental observations. The diffusion coefficients of neutral atoms are found to be greater in methanol than that in water and they decrease monotonically with solute size, whereas the diffusion coefficients of the corresponding ions are found to be smaller in methanol. Accordingly, an ion experiences a smaller Stokes friction and a higher dielectric friction in methanol than in water. These contrasting effects are believed to be responsible for the shift of the maximum of ion diffusion toward a larger ion size when compared with similar anomalous size dependence in liquid water.     

Hydration of aprotic donor solvents studied by means of FTIR spectroscopy

The paper attempts to explain the mutual influence of nonpolar and electron-donor groups on solute hydration, the problem of big importance for biological aqueous systems. Aprotic organic solvents have been used as model solutes, differing in electron-donating power. Hydration of acetonitrile, acetone, 2-butanone, and triethylamine has been studied by HDO and (partially) H2O spectra. The quantitative version of difference spectra method has been applied to determine solute-affected water spectra. Analysis of the data suggests that solvent-water interaction via the donor center of the solute is averaged between water-water interactions around the solute. Such behavior can be simply explained by the model of solute rotating in a cavity of water structure, which is formed by clathratelike hydrogen-bonded water network. On the basis of the band shape of solute-affected HDO spectra and the corresponding distribution of intermolecular distances, the criterion for hydrophobic type hydration has been proposed. From that point of view, all the studied solutes could be treated as hydrophobic ones. The limiting band position and the corresponding intermolecular distance of affected water, gained with increasing electron-donating power of solutes, has been inferred from the data obtained. These observations are important for interpretation of vibrational spectra of water as well as for volumetric measurements of solutions. The simple model of hydration, proposed to better justify the results, connects the values obtained from the methods providing microscopic and macroscopic characteristics of the system studied.     

Quantized hydration energies of ions and structure of hydration shell from the experimental gas-phase data

With previous data on alkali metal and halide ions included [Rais, J.; Okada, T. Anal. Sci. 2006, 22, 533], we analyzed rather broad data on ionic hydration from the point of view of gaseous cluster energetics. We have now added alkaline earth cations, Zn(2+), H(+), OH(-), Cu(+), Ag(+), Bi(+), Pb(+), and alkylammonium cations. The present analysis revealed the octa-coordinated nature of alkaline earth cations, which is not fully pronounced for Be(2+) and Zn(2+), existence of Eigen protonium complex, which is trigonally hydrated, and particular property of the first OH-, H(2)O cluster. Whereas these findings are generally in accordance with theoretical model calculation studies, we have foreseen in addition tetrahedral hydration for halide anions and Rb(+) and Cs(+), as well as for alkylammonium ions. The obtained picture of the quantized solvation of ions is mirrored in the ionization potentials of outer electrons of pertinent atoms. This is a second independent phenomenon, and together, they invoked a common pattern formation ("Aufbau") obeying tetra- and octa-coordinated principles.     

Hydrogen Bonding in Liquid Water and in the Hydration Shell of Salts

A near‐IR spectral study on pure water and aqueous salt solutions is used to investigate stoichiometric concentrations of different types of hydrogen‐bonded water species in liquid water and in water comprising the hydration shell of salts. Analysis of the thermodynamics of hydrogen‐bond formation signifies that hydrogen‐bond making and breaking processes are dominated by enthalpy with non‐negligible heat capacity effects, as revealed by the temperature dependence of standard molar enthalpies of hydrogen‐bond formation and from analysis of the linear enthalpy–entropy compensation effects. A generalized method is proposed for the simultaneous calculation of the spectrum of water in the hydration shell and hydration number of solutes. Resolved spectra of water in the hydration shell of different salts clearly differentiate hydrogen bonding of water in the hydration shell around cations and anions. A comparison of resolved liquid water spectra and resolved hydration‐shell spectra of ions highlights that the ordering of absorption frequencies of different kinds of hydrogen‐bonded water species is also preserved in the bound state with significant changes in band position, band width, and band intensity because of the polarization of water molecules in the vicinity of ions.     

Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. 1. Approximate expression for the free energy of hydrophobic association based on a Gaussian-overlap model

A physics-based model is proposed to derive approximate analytical expressions for the cavity component of the free energy of hydrophobic association of spherical and spheroidal solutes in water. The model is based on the difference between the number and context of the water molecules in the hydration sphere of a hydrophobic dimer and of two isolated hydrophobic solutes. It is assumed that the water molecules touching the convex part of the molecular surface of the dimer and those in the hydration spheres of the monomers contribute equally to the free energy of solvation, and those touching the saddle part of the molecular surface of the dimer result in a more pronounced increase in free energy because of their more restricted mobility (entropy loss) and fewer favorable electrostatic interactions with other water molecules. The density of water in the hydration sphere around a single solute particle is approximated by the derivative of a Gaussian centered on the solute molecule with respect to its standard deviation. On the basis of this approximation, the number of water molecules in different parts of the hydration sphere of the dimer is expressed in terms of the first and the second mixed derivatives of the two Gaussians centered on the first and second solute molecules, respectively, with respect to the standard deviations of these Gaussians, and plausible analytical expressions for the cavity component of the hydrophobic-association energy of spherical and spheroidal solutes are introduced. As opposed to earlier hydration-shell models, our expressions reproduce the desolvation maxima in the potentials of mean force of pairs of nonpolar solutes in water, and their advantage over the models based on molecular-surface area is that they have continuous gradients in the coordinates of solute centers.     

"Up" versus "down" alignment and hydration structures of solutes at the air/water interface revealed by heterodyne-detected electronic sum frequency generation with classical molecular dynamics simulation

We unambiguously demonstrate the "up" versus "down" alignment of a pair of prototypical solute molecules adsorbed at the air/water interface for the first time using heterodyne-detected electronic sum frequency generation spectroscopy. This molecular alignment is also reproduced by classical molecular dynamics (MD) simulation theoretically. Furthermore, the MD simulation indicates distinctly different interface-specific hydration structures around the two solute molecules, which dictate the molecular alignment at the interface. It is concluded that the hydrophilicity difference between the terminal functional groups of the solute governs the molecular orientation and surrounding hydration structures at the interface.