Structural Properties of High and Low Density Water in a Supercooled Aqueous Solution of Salt

Glassy dynamics of water in TIP4P/Ice aqueous solutions of trehalose in comparison with the bulk phase

We perform molecular dynamics simulations of TIP4P/Ice water in solution with trehalose for 3.65 and 18.57 wt. % concentrations and of bulk TIP4P/Ice water at ambient pressure, to characterize the structure and dynamics of water in a sugar aqueous solution in the supercooled region. We find here that TIP4P/Ice water in solution with trehalose molecules follows the Mode Coupling Theory and undergoes a fragile to strong transition up to the highest concentration investigated, similar to the bulk. Moreover, we perform a Mode Coupling Theory test, showing that the Time Temperature Superposition principle holds for both bulk TIP4P/Ice water and for TIP4P/Ice water in the solutions and we calculate the exponents of the theory. The direct comparison of the dynamical results for bulk water and water in the solutions shows upon cooling along the isobar a fastening of water dynamics for lower temperatures, T < 240 K. We found that the counter-intuitive behavior for the low temperature solutions can be explained with the diffusion anomaly of water leading us to the conclusion that the fastening observed below T = 240 K in water dynamics is only fictitious, due to the fact that the density of water molecules in the solutions is higher than the density of the bulk at the same temperature and pressure. This result should be taken into account in experimental investigations which are often carried out at constant pressure.

Supercooled solutions of sodium perchlorate in TIP4P/2005 water: The effect of martian solutes on thermodynamics and structure

We study the thermodynamic behavior of sodium perchlorate solutions in supercooled water through molecular dynamics numerical simulations. These solutions are of special interest because of the recent experimental results that led to hypothesize the presence of liquid water in perchlorate solutions beneath the Martian soil. We model water using the TIP4P/2005 potential. The results we obtain for solutions with concentrations 1.63 and 15.4 wt% are in agreement with those of a system undergoing a liquid-liquid phase transition where the liquid-liquid critical point shifts to slightly higher temperatures and lower pressures. The structure of the system is also analyzed, and we come to the conclusion that, even at the highest concentration considered, water retains its anomalous behavior.

Insights into the Emerging Networks of Voids in Simulated Supercooled Water

The structural evolution of supercooled liquid water as we approach the glass transition temperature continues to be an active area of research. Here, we use molecular dynamics simulations of TIP4P/ice water to study the changes in the connected regions of empty space within the liquid, which we investigate using the Voronoi-voids network. We observe two important features: supercooling enhances the fraction of nonspherical voids and different sizes of voids tend to cluster forming a percolating network. By examining order parameters such as the local structure index (LSI), tetrahedrality and topological defects, we show that water molecules near large void clusters tend to be slightly more tetrahedral than those near small voids, with a lower population of under- and overcoordinated defects. We show further that the distribution of closed rings of water molecules around small and large void clusters maintain a balance between 6 and 7 membered rings. Our results highlight the changes of the dual voids and water network as a structural hallmark of supercooling and provide insights into the molecular origins of cooperative effects underlying density fluctuations on the subnanometer and nanometer length scale. In addition, the percolation of the voids and the hydrogen bond network around the voids may serve as useful order parameters to investigate density fluctuations in supercooled water.

Brine rejection and hydrate formation upon freezing of NaCl aqueous solutions

Studying the freezing of saltwater on a molecular level is of fundamental importance for improving freeze desalination techniques.

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.

Unravelling the contribution of local structures to the anomalies of water: The synergistic action of several factors

We investigate the microscopic origin of water's anomalies by inspecting the hydrogen bond network (HBN) and the spatial organization of low-density-liquid (LDL) like and high-density-liquid (HDL) like environments. Specifically, we simulate-via classical molecular dynamics simulations-the isobaric cooling of a sample composed of 512 water molecules from ambient to deeply undercooled conditions at three pressures, namely, 1 bar, 400 bars, and 1000 bars. In correspondence with the Widom line (WL), (i) the HDL-like dominating cluster undergoes fragmentation caused by the percolation of LDL-like aggregates following a spinodal-like kinetics; (ii) such fragmentation always occurs at a "critical" concentration of ∼20%-30% in LDL; (iii) the HBN within LDL-like environments is characterized by an equal number of pentagonal and hexagonal rings that create a state of maximal frustration between a configuration that promotes crystallization (hexagonal ring) and a configuration that hinders it (pentagonal ring); (iv) the spatial organization of HDL-like environments shows a marked variation. Moreover, the inspection of the global symmetry shows that the intermediate-range order decreases in correspondence with the WL and such a decrease becomes more pronounced upon increasing the pressure, hence supporting the hypothesis of a liquid-liquid critical point. Our results reveal and rationalize the complex microscopic origin of water's anomalies as the cooperative effect of several factors acting synergistically. Beyond implications for water, our findings may be extended to other materials displaying anomalous behaviours.

Calorimetric study of water's two glass transitions in the presence of LiCl

A DSC study of dilute glassy LiCl aqueous solutions in the water-dominated regime provides direct evidence of a glass-to-liquid transition in expanded high density amorphous (eHDA)-type solutions. Similarly, low density amorphous ice (LDA) exhibits a glass transition prior to crystallization to ice Ic. Both glass transition temperatures are independent of the salt concentration, whereas the magnitude of the heat capacity increase differs. By contrast to pure water, the glass transition endpoint for LDA can be accessed in LiCl aqueous solutions above 0.01 mole fraction. Furthermore, we also reveal the endpoint for HDA's glass transition, solving the question on the width of both glass transitions. This suggests that both equilibrated HDL and LDL can be accessed in dilute LiCl solutions, supporting the liquid-liquid transition scenario to understand water's anomalies.

Effect of solute nature on the polyamorphic transition in glassy polyol aqueous solutions

I examined the polyamorphic behavior of glassy dilute aqueous solutions of polyols (ethylene glycol, glycerol, meso-erythritol, xylitol, and D-sorbitol) under pressure at low temperatures. Although the volume change of the glassy aqueous solution varied continuously against pressure, the rate of the volume change appeared to vary discontinuously at the onset pressure of the gradual polyamorphic transition. It is thought that low-density liquid-like solvent water and high-density liquid-like solvent water coexist during the transition. Moreover, the existence of a solute induces the shift of polyamorphic transition to the lower-pressure side. The effect of a solute on the polyamorphic transition becomes larger in the order ethylene glycol, glycerol, meso-erythritol, xylitol, and D-sorbitol. Therefore, the solute can become a variable controlling the polyamorphic state of liquid water. This experimental result suggests that the metastable-equilibrium phase boundary between the low-density and the high-density amorphs for pure water is likely to be located at 0.22-0.23 GPa at about 150 K, which is slightly larger than the previously estimated pressure. Moreover, the solute-nature dependence on the polyamorphic transition seems to connect to that on the homogeneous nucleation temperature of polyol aqueous solution at ambient pressure. The region in which a low-density liquid appears coincides with the region in which the nucleus of ice Ih appears, suggesting that the formation of a low-density liquid is a precursory phenomenon of the nucleation of ice Ih.

Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route

The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.

Water: A Tale of Two Liquids

Water is the most abundant liquid on earth and also the substance with the largest number of anomalies in its properties. It is a prerequisite for life and as such a most important subject of current research in chemical physics and physical chemistry. In spite of its simplicity as a liquid, it has an enormously rich phase diagram where different types of ices, amorphous phases, and anomalies disclose a path that points to unique thermodynamics of its supercooled liquid state that still hides many unraveled secrets. In this review we describe the behavior of water in the regime from ambient conditions to the deeply supercooled region. The review describes simulations and experiments on this anomalous liquid. Several scenarios have been proposed to explain the anomalous properties that become strongly enhanced in the supercooled region. Among those, the second critical-point scenario has been investigated extensively, and at present most experimental evidence point to this scenario. Starting from very low temperatures, a coexistence line between a high-density amorphous phase and a low-density amorphous phase would continue in a coexistence line between a high-density and a low-density liquid phase terminating in a liquid–liquid critical point, LLCP. On approaching this LLCP from the one-phase region, a crossover in thermodynamics and dynamics can be found. This is discussed based on a picture of a temperature-dependent balance between a high-density liquid and a low-density liquid favored by, respectively, entropy and enthalpy, leading to a consistent picture of the thermodynamics of bulk water. Ice nucleation is also discussed, since this is what severely impedes experimental investigation of the vicinity of the proposed LLCP. Experimental investigation of stretched water, i.e., water at negative pressure, gives access to a different regime of the complex water diagram. Different ways to inhibit crystallization through confinement and aqueous solutions are discussed through results from experiments and simulations using the most sophisticated and advanced techniques. These findings represent tiles of a global picture that still needs to be completed. Some of the possible experimental lines of research that are essential to complete this picture are explored.

Free energetics of carbon nanotube association in aqueous inorganic NaI salt solutions: Temperature effects using all-atom molecular dynamics simulations

In this study, we examine the temperature dependence of free energetics of nanotube association using graphical processing unit-enabled all-atom molecular dynamics simulations (FEN ZI) with two (10,10) single-walled carbon nanotubes in 3 m NaI aqueous salt solution. Results suggest that the free energy, enthalpy and entropy changes for the association process are all reduced at the high temperature, in agreement with previous investigations using other hydrophobes. Via the decomposition of free energy into individual components, we found that solvent contribution (including water, anion, and cation contributions) is correlated with the spatial distribution of the corresponding species and is influenced distinctly by the temperature. We studied the spatial distribution and the structure of the solvent in different regions: intertube, intratube and the bulk solvent. By calculating the fluctuation of coarse-grained tube-solvent surfaces, we found that tube-water interfacial fluctuation exhibits the strongest temperature dependence. By taking ions to be a solvent-like medium in the absence of water, tube-anion interfacial fluctuation shows similar but weaker dependence on temperature, while tube-cation interfacial fluctuation shows no dependence in general. These characteristics are discussed via the malleability of their corresponding solvation shells relative to the nanotube surface. Hydrogen bonding profiles and tetrahedrality of water arrangement are also computed to compare the structure of solvent in the solvent bulk and intertube region. The hydrophobic confinement induces a relatively lower concentration environment in the intertube region, therefore causing different intertube solvent structures which depend on the tube separation. This study is relevant in the continuing discourse on hydrophobic interactions (as they impact generally a broad class of phenomena in biology, biochemistry, and materials science and soft condensed matter research), and interpretations of hydrophobicity in terms of alternative but parallel signatures such as interfacial fluctuations, dewetting transitions, and enhanced fluctuation probabilities at interfaces.

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.

A different view of structure-making and structure-breaking in alkali halide aqueous solutions through x-ray absorption spectroscopy

X-ray absorption spectroscopy measured in transmission mode was used to study the effect of alkali and halide ions on the hydrogen-bonding (H-bonding) network of water. Cl(-) and Br(-) are shown to have insignificant effect on the structure of water while I(-) locally weakens the H-bonding, as indicated by a sharp increase of the main-edge feature in the x-ray absorption spectra. All alkali cations act as structure-breakers in water, weakening the H-bonding network. The spectral changes are similar to spectra of high density ices where the 2nd shell has collapsed due to a break-down of the tetrahedral structures, although here, around the ions, the breakdown of the local tetrahedrality is rather due to non-directional H-bonding to the larger anions. In addition, results from temperature-dependent x-ray Raman scattering measurements of NaCl solution confirm the H-bond breaking effect of Na(+) and the effect on the liquid as similar to an increase in temperature.

The associating lattice gas in the presence of interacting solutes

We have investigated the phase diagram of a statistical model for hydrogen-bonding solutions for polar solutes. The structured solvent is represented by an associating lattice gas, which presents anomalous density and liquid-liquid coexistence. Polar solute particles and solvent particles interact attractively, while the solvent-solvent interaction is made directional through bonding arms, which mimic hydrogen bonds. The model behavior is obtained via Monte Carlo simulations in the grand-canonical ensemble, for different sets of parameters. For small solute chemical potential and weak attraction between solute and solvent particles, addition of solute yields a shift in the transition lines of the pure solvent. This is the scenario explored by different authors, in the pursuit of stabilizing the water liquid-liquid coexistence line. However, as we show, in the case of larger solute chemical potentials, or of stronger solute-solvent attractions, new phases may arise.

Relation between the two-body entropy and the relaxation time in supercooled water

The two-body excess entropy of supercooled water is calculated from the radial distribution functions obtained from computer simulation of the TIP4P model for different densities upon supercooling. This quantity is considered in connection with the relaxation time of the self intermediate scattering function. The relaxation time shows a mode coupling theory (MCT) behavior in the region of mild supercooling and a strong behavior in the deep supercooled region. We find here that the two-body entropy is connected to the relaxation time and shows a logarithmic behavior with an apparent asymptotic divergence at the mode coupling crossover temperature. There is also evidence of a change in behavior of the two-body entropy upon crossing from the fragile (hopping-free) state to the strong (hopping-dominated) state of supercooled water, and the relation that connects the two-body entropy and the relxation time in the MCT region no longer holds.

Volumetric determination of apolipoprotein stoichiometry of circulating HDL subspecies

Although HDL is inversely correlated with coronary heart disease, elevated HDL-cholesterol is not always protective. Additionally, HDL has biological functions that transcend any antiatherogenic role: shotgun proteomics show that HDL particles contain 84 proteins (latest count), many correlating with antioxidant and anti-inflammatory properties of HDL. ApoA-I has been suggested to serve as a platform for the assembly of these protein components on HDL with specific functions - the HDL proteome. However, the stoichiometry of apoA-I in HDL subspecies is poorly understood. Here we use a combination of immunoaffinity chromatography data and volumetric analysis to evaluate the size and stoichiometry of LpA-I and LpA-I,A-II particles. We conclude that there are three major LpA-I subspecies: two major particles, HDL[4] in the HDL3 size range (d = 85.0 ± 1.2 Å) and HDL[7] in the HDL2 size range (d = 108.5 ± 3.8 Å) with apoA-I stoichiometries of 3 and 4, respectively, and a small minor particle, HDL[1] (d = 73.8 ± 2.1Å) with an apoA-I stoichiometry of 2. Additionally, we conclude that the molar ratio of apolipoprotein to surface lipid is significantly higher in circulating HDL subspecies than in reconstituted spherical HDL particles, presumably reflecting a lack of phospholipid transfer protein in reconstitution protocols.

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.

Effect of hydrophobic environments on the hypothesized liquid-liquid critical point of water

The complex behavior of liquid water, along with its anomalies and their crucial role in the existence of life, continue to attract the attention of researchers. The anomalous behavior of water is more pronounced at subfreezing temperatures and numerous theoretical and experimental studies are directed towards developing a coherent thermodynamic and dynamic framework for understanding supercooled water. The existence of a liquid-liquid critical point in the deep supercooled region has been related to the anomalous behavior of water. However, the experimental study of supercooled water at very low temperatures is hampered by the homogeneous nucleation of the crystal. Recently, water confined in nanoscopic structures or in solutions has attracted interest because nucleation can be delayed. These systems have a tremendous relevance also for current biological advances; e.g., supercooled water is often confined in cell membranes and acts as a solvent for biological molecules. In particular, considerable attention has been recently devoted to understanding hydrophobic interactions or the behavior of water in the presence of apolar interfaces due to their fundamental role in self-assembly of micelles, membrane formation and protein folding. This article reviews and compares two very recent computational works aimed at elucidating the changes in the thermodynamic behavior in the supercooled region and the liquid-liquid critical point phenomenon for water in contact with hydrophobic environments. The results are also compared to previous reports for water in hydrophobic environments.

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.

Melting of the precipitated ice IV in LiCl aqueous solution and polyamorphism of water

Melting of the precipitated ice IV in supercooled LiCl-H(2)O solution was studied in the range of 0-0.6 MPa and 160-270 K. Emulsified solution was used to detect this metastable transition. Ice IV was precipitated from the aqueous solution of 2.0 mol % LiCl (or 4.8 mol % LiCl) in each emulsion particle at low-temperature and high-pressure conditions, and the emulsion was decompressed at different temperatures. The melting of ice IV was detected from the temperature change of the emulsified sample during the decompression. There was an apparently sudden change in the slope of the ice IV melting curve (liquidus) in the pressure-temperature diagram. At the high-pressure and high-temperature side of the change, the solute-induced freezing point depression was observed. At the low-pressure and low-temperature side, ice IV transformed into ice Ih on the decompression, and the transition was almost unrelated to the concentration of LiCl. These experimental results were roughly explained by the presumed existence of two kinds of liquid water (low-density liquid water and high-density liquid water), or polyamorphism in water, and by the simple assumption that LiCl dissolved maily in high-density liquid water.