Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water

Toward mitigating the impact of non-bulk defects on describing water structure in salt aqueous solutions: Characterizing solution density with a network-based structural indicator

Structural indicators, also known as structural descriptors, including order parameters, have been proposed to quantify the structural properties of water to account for its anomalous behaviors. However, these indicators, mainly designed for bulk water, are not naturally transferrable to the vicinity of ions due to disruptions in the immediate neighboring space and a resulting loss of feature completeness. To address these non-bulk defects, we introduced a structural indicator that draws on the concept of clique number from graph theory and the criterion in agglomerative clustering, denoted as the average cluster number. This structural indicator aims to discern intrinsic structural characteristics within the water molecules regardless of the ions occupying the neighboring space, without requiring additional corrections. From molecular dynamics simulation results for neat water and salt aqueous solutions utilizing the TIP4P/2005 water model and the Madrid-2019 force field, we characterized the variations in densities with temperature using this network-based indicator, thereby demonstrating its practical utility. The findings suggest that at lower temperatures, the addition of ions disrupts the intrinsic structure of water molecules, with this effect diminishing as the temperature rises. Cations with larger charge density tend to induce stronger disruptions. This study highlights the importance of mitigating the impact of non-bulk defects before applying the indicators to analyze water's intrinsic structural properties in solutions. By doing so, the relationship between changes in water structure and solution behaviors can be more accurately assessed.

Supercritical density fluctuations and structural heterogeneity in supercooled water-glycerol microdroplets

Recent experiments and theoretical studies strongly indicate that water exhibits a liquid-liquid phase transition (LLPT) in the supercooled domain. An open question is how the LLPT of water can affect the properties of aqueous solutions. Here, we study the structural and thermodynamic properties of supercooled glycerol-water microdroplets at dilute conditions (χg = 3.2% glycerol mole fraction). The combination of rapid evaporative cooling with femtosecond X-ray scattering allows us to outrun crystallization and gain access to the deeply supercooled regime down to T = 229.3 K. We find that the density fluctuations of the glycerol-water solution or, equivalently, its isothermal compressibility, κT, increases upon cooling. This is confirmed by molecular dynamics simulations, which indicate that the presence of glycerol shifts the temperature of maximum κT from T = 230 K in pure water down to T = 223 K in the solution. Our findings elucidate the interplay between the complex behavior of water, including its LLPT, and the properties of aqueous solutions at low temperatures, which can have practical consequences in cryogenic biological applications and cryopreservation techniques.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

A continuum of amorphous ices between low-density and high-density amorphous ice

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.

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.