Anomalies in supercooled NaCl aqueous solutions: A microscopic perspective

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials

The diffusivity of water in aqueous cesium iodide solutions is larger than that in neat liquid water and vice versa for sodium chloride solutions. Such peculiar ion-specific behavior, called anomalous diffusion, is not reproduced in typical force field based molecular dynamics (MD) simulations due to inadequate treatment of ion-water interactions. Herein, this hurdle is tackled by using machine learned atomic potentials (MLPs) trained on data from density functional theory calculations. MLP based atomistic MD simulations of aqueous salt solutions reproduce experimentally determined thermodynamic, structural, dynamical, and transport properties, including their varied trends in water diffusivities across salt concentration. This enables an examination of their intermolecular structure to unravel the microscopic underpinnings of the differences in their transport properties. While both ions in CsI solutions contribute to the faster diffusion of water molecules, the competition between the heavy retardation by Na ions and the slight acceleration by Cl ions in NaCl solutions reduces their water diffusivity.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Examining the Hofmeister Series through Activation Energies: Water Diffusion in Aqueous Alkali-Halide Solutions

The effect of ions on the properties of aqueous solutions is often categorized in terms of the Hofmeister series that ranks them from chaotropes ("structure-breakers"), which weaken the surrounding hydrogen-bond network to kosmotropes ("structure-makers"), which enhance it. Here, we investigate the Hofmeister series in ∼1 M sodium-halide solutions using molecular dynamics simulations to calculate the effect of the identity and proximity of the halide anion on both the water diffusion coefficient and its activation energy. A recently developed method for calculating the activation energy from a single-temperature simulation is used, which also permits a rigorous decomposition into contributions from different interactions and motions. The mechanisms of the salt effects on the water dynamics are explored by separately considering water molecules based on their location relative to the ions. The results show that water diffusion is accelerated moving down the halide group from F^- to I^-. The behavior of the diffusion activation energy, E_a, is more complex, indicating a significant role for entropic effects. However, water molecules in the first or second solvation shell of an ion exhibit a decrease in E_a moving down the halide series and Na⁺ exhibits a larger effect than any of the anions. The E_a for water molecules within the second solvation shell of an ion are modest, indicating a short-ranged nature of the ion influence.

Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations

One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.

Pressure-induced transformations in LiCl-H2O at 77 K

A systematic study of the properties of high-density amorphous ice (HDA) in the presence of increasing amounts of salt is missing, especially because it is challenging to avoid ice crystallization upon cooling the pressurized liquid. In order to be able to study HDA also in the presence of small amounts of salt, we have investigated the transformation behaviour of quenched aqueous LiCl solutions (mole fraction x < 0.25) upon pressurization in a piston-cylinder setup at 77 K. The sample properties were characterized by in situ dilatometry under high pressure conditions and after recovery by ex situ powder X-ray diffraction (XRD) and differential scanning calorimetry (DSC) at ambient pressure. Two regimes can be identified, with a rather sharp switch at about x = 0.12. At x < 0.12 the samples show the phenomenology also known for pure water samples. They are composed mainly of hexagonal ice (Ih) and experience pressure-induced amorphization to HDA at P > 1 GPa. The observed densification is consistent with the idea that a freeze concentrated LiCl solution of x = 0.14 (R = 6) segregates, which transforms to the glassy state upon cooling, and that the densification is only due to the Ih → HDA transition. Also the XRD patterns and DSC scans are almost unaffected by the presence of the segregated glassy LiCl solution. Upon heating at ambient pressure HDA experiences the polyamorphic transition to low-density amorphous ice (LDA) at ∼120 K, even at x ∼ 0.10. Based on the latent heat evolved in the transition we suggest that almost all water in the sample transforms to an LDA-like state, even the water in the vicinity of the ions. The glassy LiCl solution acts as a spectator that does not shift the transformation temperature significantly and experiences a glass-to-liquid transition at ∼140 K prior to the crystallization to cubic ice. By contrast, at x > 0.12 the phenomenology completely changes and is now dominated by the salt. Hexagonal ice no longer forms upon quenching the LiCl solution, but instead LDA forms. A broad pressure-induced transformation at >0.6 GPa can be attributed to the densification of LDA, the glassy LiCl solution and/or glassy hydrates.

Sudden switchover between the polyamorphic phase separation and the glass-to-liquid transition in glassy LiCl aqueous solutions

Lithium chloride aqueous solutions (LiClaq solutions) below 10 mol.% are vitrified by cooling from room temperature to 77 K at 0.3 GPa. We examine the solvent state of the glassy sample and its transformation by heating at 1 atm using low-temperature differential scanning calorimetry and Raman spectroscopy. This experimental study suggests strongly that the solvent state of the glassy LiClaq solution closely relates to the state of high-density amorphous ice. Moreover, we reconfirm that the separation into the low-density amorphous ice and the glassy highly concentrated LiClaq solution occurs in the glassy dilute LiClaq solution at ∼130 K, not the glass-to-liquid transition which is commonly observed in the glassy LiClaq solution above ∼10 mol.%. In order to interpret the sudden switchover between the glass-to-liquid transition and the phase separation at ∼10 mol.%, we propose a state diagram of LiClaq solution which connects with a polyamorphic state diagram of pure water and discuss a possibility that the electric field induces a polyamorphic transition of water.

A molecular dynamics study of the equation of state and the structure of supercooled aqueous solutions of methanol

We perform molecular dynamics computer simulations in order to study the equation of state and the structure of supercooled aqueous solutions of methanol at methanol mole fractions x(m) = 0.05 and x(m) = 0.10. We model the solvent using the TIP4P/2005 potential and the methanol using the OPLS-AA force field. We find that for x(m) = 0.05 the behavior of the equation of state, studied in the P - T and P - ρ planes, is consistent with the presence of a liquid-liquid phase transition, reminiscent of that previously found for x(m) = 0. We estimate the position of the liquid-liquid critical point to be at T = 193 K, P = 96 MPa, and ρ = 1.003 g/cm(3). When the methanol mole fraction is doubled to x(m) = 0.10 no liquid-liquid transition is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-methanol structure in the two solutions. We find that down to low temperature methanol can be incorporated into the water structure for both x(m) = 0.05 and x(m) = 0.10.

Fragile-to-strong crossover coupled to the liquid-liquid transition in hydrophobic solutions

Using discrete molecular dynamics simulations we study the relation between the thermodynamic and diffusive behaviors of a primitive model of aqueous solutions of hydrophobic solutes consisting of hard spheres in the Jagla particles solvent, close to the liquid-liquid critical point of the solvent. We find that the fragile-to-strong dynamic transition in the diffusive behavior is always coupled to the low-density-high-density liquid transition. Above the liquid-liquid critical pressure, the diffusivity crossover occurs at the Widom line, the line along which the thermodynamic response functions show maxima. Below the liquid-liquid critical pressure, the diffusivity crossover occurs when the limit of mechanical stability lines are crossed, as indicated by the hysteresis observed when going from high to low temperature and vice versa. These findings show that the strong connection between dynamics and thermodynamics found in bulk water persists in hydrophobic solutions for concentrations from low to moderate, indicating that experiments measuring the relaxation time in aqueous solutions represent a viable route for solving the open questions in the field of supercooled water.

Liquid-liquid coexistence in NaCl aqueous solutions: a simulation study of concentration effects

In this paper we investigate by means of molecular dynamics computer simulations how the hypothesized liquid-liquid critical point of water shifts in supercooled aqueous solutions of salt as a function of concentration. We study sodium chloride solutions in TIP4P water, NaCl(aq), for concentrations c = 1.36 mol/kg and c = 2.10 mol/kg. The liquid-liquid critical point is found up to the highest concentration investigated, and its position in the P-T plane shifts to higher temperatures and lower pressures upon increasing concentration. For c = 2.10 mol/kg it is also located very close to the temperature of maximum density line of the system. The results are discussed and compared with previous results for bulk TIP4P water and for c = 0.67 mol/kg NaCl(aq) and with experimental findings. We observe a progressive shrinkage of the low-density liquid region when the concentration of salt increases; this suggests an eventual disappearance of the liquid-liquid coexistence upon further increase of NaCl concentration.