Communication: Mean-field theory of water-water correlations in electrolyte solutions

Dissolving salt is not equivalent to applying a pressure on water

Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

By advanced machine learning techniques, first-principles simulations find that dissolving salt in water does not change water structure drastically. It is contrary to the notion of “pressure effect” which has been widely applied over past 25 years.

Hyaluronan orders water molecules in its nanoscale extended hydration shells

Hyaluronan (HA) is an anionic, highly hydrated bio-polyelectrolyte found in the extracellular environment, like the synovial fluid between joints. We explore the extended hydration shell structure of HA in water using femtosecond elastic second-harmonic scattering (fs-ESHS). HA enhances orientational water-water correlations. Angle-resolved fs-ESHS measurements and nonlinear optical modeling show that HA behaves like a flexible chain surrounded by extended shells of orientationally correlated water. We describe several ways to determine the concentration-dependent size and shape of a polyelectrolyte in water, using the amount of water oriented by the polyelectrolyte charges as a contrast agent. The spatial extent of the hydration shell is determined via temperature-dependent measurements and can reach up to 475 nm, corresponding to a length of 1600 water molecules. A strong isotope effect, stemming from nuclear quantum effects, is observed when light water (H2O) is replaced by heavy water (D2O), amounting to a factor of 4.3 in the scattered SH intensity.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Polyelectrolytes induce water-water correlations that result in dramatic viscosity changes and nuclear quantum effects

Ions interact with water via short-ranged ion-dipole interactions. Recently, an additional unexpected long-ranged interaction was found: The total electric field of ions influences water-water correlations over tens of hydration shells, leading to the Jones Ray effect, a 0.3% surface tension depression. Here, we report such long-range interactions contributing substantially to both molecular and macroscopic properties. Femtosecond elastic second harmonic scattering (fs-ESHS) shows that long-range electrostatic interactions are remarkably strong in aqueous polyelectrolyte solutions, leading to an increase in water-water correlations. This increase plays a role in the reduced viscosity, which changes more than two orders of magnitude with polyelectrolyte concentration. Using D₂O instead of H₂O shifts both the fs-ESHS and the viscosity curve by a factor of ~10 and reduces the maximum viscosity value by 20 to 300%, depending on the polyelectrolyte. These phenomena cannot be explained using a mean-field approximation of the solvent and point to nuclear quantum effects.

The Jones-Ray Effect Is Not Caused by Surface-Active Impurities

Pure aqueous electrolyte solutions display a minimum in surface tension at concentrations of 2 ± 1 mM. This effect has been a source of controversy since it was first reported by Jones and Ray in the 1930s. The Jones-Ray effect has frequently been dismissed as an artifact linked to the presence of surface-active impurities. Herein we systematically consider the effect of surface-active impurities by purposely adding nanomolar concentrations of surfactants to dilute electrolyte solutions. Trace amounts of surfactant are indeed found to decrease the surface tension and influence the surface chemistry. However, surfactants can be removed by repeated aspiration and stirring cycles, which eventually deplete the surfactant from solution, creating a pristine surface. Upon following this cleaning procedure, a reduction in the surface tension by millimolar concentrations of salt is still observed. Consequently, we demonstrate that the Jones-Ray effect is not caused by surface-active impurities.

What Does Second-Harmonic Scattering Measure in Diluted Electrolytes?

We derive a theoretical expression of the second harmonic scattering signal in diluted electrolytes compared with bulk water. We show that the enhancement of the signal with respect to pure water observed recently for electrolytes at very low dilution in the micromolar range is a mere manifestation of the Debye screening that makes the infinite-range dipole-dipole solvent correlations in 1/ r³ disappear as soon as the ionic concentration becomes finite. In q space, this translates into a correlation function having a well known singular behavior around q = 0, which drives the observed ionic effects. We find that the signal is independent of the ion-induced long-range behavior of the function ⟨cos ϕ( r)⟩ that has been recently discussed. We find also that the enhancement depends on the experimental geometry and occurs only for in-plane polarization detection, as observed experimentally. On the contrary, the measured isotope effect between light and heavy water cannot be fully explained.

The spatial range of protein hydration

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

Salt-induced Long-to-Short Range Orientational Transition in Water

We report the long-range orientational organization of water using polarization-resolved second harmonic scattering operated in a right-angle configuration. A transition is observed between the neat water orientational organization involving an azimuthal molecular orientation distribution towards a radial molecular orientation distribution when salt is added. These two orientational phases are quantitatively described using a molecular model of the second harmonic scattering response. It is observed that the long-range correlation present in the neat water phase abruptly disappears and is replaced by a shorter range correlation centered around the ions as the salt concentration is increased.

Water at surfaces with tunable surface chemistries

Aqueous interfaces are ubiquitous in natural environments, spanning atmospheric, geological, oceanographic, and biological systems, as well as in technical applications, such as fuel cells and membrane filtration. Where liquid water terminates at a surface, an interfacial region is formed, which exhibits distinct properties from the bulk aqueous phase. The unique properties of water are governed by the hydrogen-bonded network. The chemical and physical properties of the surface dictate the boundary conditions of the bulk hydrogen-bonded network and thus the interfacial properties of the water and any molecules in that region. Understanding the properties of interfacial water requires systematically characterizing the structure and dynamics of interfacial water as a function of the surface chemistry. In this review, we focus on the use of experimental surface-specific spectroscopic methods to understand the properties of interfacial water as a function of surface chemistry. Investigations of the air-water interface, as well as efforts in tuning the properties of the air-water interface by adding solutes or surfactants, are briefly discussed. Buried aqueous interfaces can be accessed with careful selection of spectroscopic technique and sample configuration, further expanding the range of chemical environments that can be probed, including solid inorganic materials, polymers, and water immiscible liquids. Solid substrates can be finely tuned by functionalization with self-assembled monolayers, polymers, or biomolecules. These variables provide a platform for systematically tuning the chemical nature of the interface and examining the resulting water structure. Finally, time-resolved methods to probe the dynamics of interfacial water are briefly summarized before discussing the current status and future directions in studying the structure and dynamics of interfacial water.

Orientational ordering of water in extended hydration shells of cations is ion-specific and is correlated directly with viscosity and hydration free energy

Specific ion effects in aqueous solutions are investigated at the molecular, nanoscopic and macroscopic levels.