Pressure-induced transformations in glassy water: A computer simulation study using the TIP4P/2005 model

Link to Densify: Topological Transitions and Origin of Hysteresis During the Compression and Decompression of Amorphous Ices

In this Letter, we study the phase transition between amorphous ices and the nature of the hysteresis cycle separating them. We discover that a topological transition takes place as the system transforms from low-density amorphous ice (LDA) at low pressures to high-density amorphous ice (HDA) at high pressures. Specifically, we uncover that the hydrogen bond network (HBN) displays qualitatively different topologies in the LDA and HDA phases: the former characterized by disentangled loop motifs, with the latter displaying topologically complex long-lived Hopf-linked and knotted configurations. At the phase transition, the transient opening of the HBN topological motifs yields mechanical fragility on the macroscale. Our results provide a detailed microscopic description of the topological nature of the phase transition and the hysteresis cycle between amorphous ices. We argue that the topological transition discovered in this work may not only improve our understanding of amorphous ices, but also represent a generic mechanism for the densification of network-forming materials.

Nuclear quantum effects on glassy water under pressure: Vitrification and pressure-induced transformations

We perform classical molecular dynamics (MD) and path-integral MD (PIMD) simulations of H2O and D2O using the q-TIP4P/F model over a wide range of temperatures and pressures to study the nuclear quantum effects (NQEs) on (i) the vitrification of liquid water upon isobaric cooling at different pressures and (ii) pressure-induced transformations at constant temperature between low-density amorphous and high-density amorphous ice (LDA and HDA) and hexagonal ice Ih and HDA. Upon isobaric cooling, classical and quantum H2O and D2O vitrify into a continuum of intermediate amorphous ices (IA), with densities in-between those of LDA and HDA (depending on pressure). Importantly, the density of the IA varies considerably if NQEs are included (similar conclusions hold for ice Ih at all pressures studied). While the structure of the IA is not very sensitive to NQE, the geometry of the hydrogen-bond (HB) is. NQE leads to longer and less linear HB in LDA, HDA, and ice Ih than found in the classical case. Interestingly, the delocalization of the H/D atoms is non-negligible and identical in LDA, HDA, and ice Ih at all pressures studied. Our isothermal compression/decompression MD/PIMD simulations show that classical and quantum H2O and D2O all exhibit LDA-HDA and ice Ih-HDA transformations, consistent with experiments. The inclusion of NQE leads to a softer HB-network, which lowers slightly the LDA/ice Ih-to-HDA transformation pressures. Interestingly, the HB in HDA is longer and less linear than in LDA, which is counterintuitive given that HDA is ≈25% denser than LDA. Overall, our results show that, while classical computer simulations provide the correct qualitative phenomenology of ice and glassy water, NQEs are necessary for a quantitative description.

Density isobar of water and melting temperature of ice: Assessing common density functionals

We investigate the density isobar of water and the melting temperature of ice using six different density functionals. Machine-learning potentials are employed to ensure computational affordability. Our findings reveal significant discrepancies between various base functionals. Notably, even the choice of damping can result in substantial differences. Overall, the outcomes obtained through density functional theory are not entirely satisfactory across most utilized functionals. All functionals exhibit significant deviations either in the melting temperature or equilibrium volume, with most of them even predicting an incorrect volume difference between ice and water. Our heuristic analysis indicates that a hybrid functional with 25% exact exchange and van der Waals damping averaged between zero and Becke-Johnson dampings yields the closest agreement with experimental data. This study underscores the necessity for further enhancements in the treatment of van der Waals interactions and, more broadly, density functional theory to enable accurate quantitative predictions for molecular liquids.

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.

The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice

We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (I_h) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H₂O and D₂O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice I_h. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice I_h increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice I_h. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient α_P(T) and bulk modulus B(T) of both LDA and ice I_h. Remarkably, the ρ(T), α_P(T), and B(T) of LDA are practically identical to those of ice I_h. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice I_h. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.

Kinetics and Mechanisms of Pressure-Induced Ice Amorphization and Polyamorphic Transitions in a Machine-Learned Coarse-Grained Water Model

Water glasses have attracted considerable attention due to their potential connection to a liquid-liquid transition in supercooled water. Here we use molecular simulations to investigate the formation and phase behavior of water glasses using the machine-learned bond-order parameter (ML-BOP) water model. We produce glasses through hyperquenching of water, pressure-induced amorphization (PIA) of ice, and pressure-induced polyamorphic transformations. We find that PIA of polycrystalline ice occurs at a lower pressure than that of monocrystalline ice and through a different mechanism. The temperature dependence of the amorphization pressure of polycrystalline ice for ML-BOP agrees with that in experiments. We also find that ML-BOP accurately reproduces the density, coordination number, and structural features of low-density (LDA), high-density (HDA), and very high-density (VHDA) amorphous water glasses. ML-BOP accurately reproduces the experimental radial distribution function of LDA but overpredicts the minimum between the first two shells in high-density glasses. We examine the kinetics and mechanism of the transformation between low-density and high-density glasses and find that the sharp nature of these transitions in ML-BOP is similar to that in experiments and all-atom water models with a liquid-liquid transition. Transitions between ML-BOP glasses occur through a spinodal-like mechanism, similar to ice crystallization from LDA. Both glass-to-glass and glass-to-ice transformations have Avrami-Kolmogorov kinetics with exponent n = 1.5 ± 0.2 in experiments and simulations. Importantly, ML-BOP reproduces the competition between crystallization and HDA→LDA transition above the glass transition temperature T_g, and separation of their time scales below T_g, observed also in experiments. These findings demonstrate the ability of ML-BOP to accurately reproduce water properties across various regimes, making it a promising model for addressing the competition between polyamorphic transitions and crystallization in water and solutions.

Free energy calculations and unbiased molecular dynamics targeting the liquid-liquid transition in water no man's land

The existence of a first-order phase transition between a low-density liquid (LDL) and a high-density liquid (HDL) form of supercooled water has been a central and highly debated issue of physics and chemistry for the last three decades. We present a computational study that allows us to determine the free-energy landscapes of supercooled water over a wide range of pressure and temperature conditions using the TIP4P/2005 force field. Our approach combines topology-based structural transformation coordinates, state-of-the-art free-energy calculation methods, and extensive unbiased molecular dynamics. All our diverse simulations cannot detect any barrier within the investigated timescales and system size, for a discontinuous transition between the LDL and HDL forms throughout the so-called "no man's land," until the onset of the solid, non-diffusive amorphous forms.

Stability and Metastability of Liquid Water in a Machine-Learned Coarse-Grained Model with Short-Range Interactions

Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice-liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice-liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice-water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.

Molecular dynamics analysis of elastic properties and new phase formation during amorphous ices transformations

Unlike conventional first-order phase transitions, the kinetics of amorphous-amorphous transitions has been much less studied. The ultrasonic experiments on the transformations between low-density and high-density amorphous ice induced by pressure or heating provided the pressure and temperature dependencies of elastic moduli. In this article, we make an attempt to build a microscopic picture of these experimentally studied transformations using the molecular dynamics method with the TIP4P/Ice water model. We study carefully the dependence of the results of elastic constants calculations on the deformation rates. The system size effects are considered as well. The comparison with the experimental data enriches our understanding of the transitions observed. Our modeling gives new information about the formation mechanisms of new phase clusters during the transition between low-density and high-density amorphous ices. We analyse the applicability of the term "nucleation" for these processes.

Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid-liquid phase transition and liquid-liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., Tg LDL(P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., Tg HDL(P) increases upon compression. (iii) NQE shift both Tg LDL(P) and Tg HDL(P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA-HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA-HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA-HDA transformation.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

The role of high-density and low-density amorphous ice on biomolecules at cryogenic temperatures: a case study with polyalanine

Experimental techniques, such as cryo-electron microscopy, require biological samples to be recovered at cryogenic temperatures (T ≈ 100 K) with water being in an amorphous ice state. However, (bulk) water can exist in two amorphous ices at P < 1 GPa, low-density amorphous (LDA) ice at low pressures and high-density amorphous ice (HDA) at high pressures; HDA is ≈20-25% denser than LDA. While fast/plunge cooling at 1 bar brings the sample into LDA, high-pressure cooling (HPC), at sufficiently high pressure, produces HDA. HDA can also be produced by isothermal compression of LDA at cryogenic temperatures. Here, we perform classical molecular dynamics simulations to study the effects of LDA, HDA, and the LDA-HDA transformation on the structure and hydration of a small peptide, polyalanine. We follow thermodynamic paths corresponding to (i) fast/plunge cooling at 1 bar, (ii) HPC at P = 400 MPa, and (iii) compression/decompression cycles at T = 80 K. While process (i) produced LDA in the system, path (iii) produces HDA. Interestingly, the amorphous ice produced in process (ii) is an intermediate amorphous ice (IA) with properties that fall in-between those of LDA and HDA. Remarkably, the structural changes in polyalanine are negligible at all conditions studied (0-2000 MPa, 80-300 K) even when water changes among the low and high-density liquid states as well as the amorphous solids LDA, IA, and HDA. The similarities and differences in the hydration of polyalanine vitrified in LDA, IA, and HDA are described. Since the studied thermodynamic paths are suitable for the cryopreservation of biomolecules, we also study the structure and hydration of polyalanine along isobaric and isochoric heating paths, which can be followed experimentally for the recovery of cryopreserved samples. Upon heating, the structure of polyalanine remains practically unchanged. We conclude with a brief discussion of the practical advantages of (a) using HDA and IA as a cryoprotectant environment (as opposed to LDA), and (b) the use of isochoric heating as a recovery process (as opposed to isobaric heating).

Experimental evidence for glass polymorphism in vitrified water droplets

The question of whether a first-order liquid-to-liquid transition is at the origin of water’s anomalous properties has been controversial since the pioneering experiments by Mishima et al. in 1985 and molecular simulations by Poole et al. in 1992. Since then, experiments aimed at shedding light on this question have been performed using amorphous ices made from crystalline ice, fueling criticism about their crystal-like nature. In the present study, we avoid crystalline ice at any time of the experiment yet still observe a first-order glass-to-glass transition in vitrified liquid droplets. This makes the strong case for glass polymorphism and the direct thermodynamic connection to the liquid-to-liquid transition at higher temperatures, dismissing the criticism voiced for three decades.

The nature of amorphous ices has been debated for more than 35 years. In essence, the question is whether they are related to ice polymorphs or to liquids. The fact that amorphous ices are traditionally prepared from crystalline ice via pressure-induced amorphization has made a clear distinction tricky. In this work, we vitrify liquid droplets through cooling at ≥10⁶ K ⋅ s⁻¹ and pressurize the glassy deposit. We observe a first order–like densification upon pressurization and recover a high-density glass. The two glasses resemble low- and high-density amorphous ice in terms of both structure and thermal properties. Vitrified water shows all features that have been reported for amorphous ices made from crystalline ice. The only difference is that the hyperquenched and pressurized deposit shows slightly different crystallization kinetics to ice I upon heating at ambient pressure. This implies a thermodynamically continuous connection of amorphous ices with liquids, not crystals.

Quantum polyamorphism in compressed distinguishable helium-4

We demonstrate that two amorphous solid states can exist in ⁴He consisting of distinguishable Boltzmann atoms under compressed conditions. The isothermal compression of normal or supercritical fluid ⁴He was conducted at 3-25 K using the isobaric-isothermal path integral centroid molecular dynamics simulation. The compression of fluid first produced the low-dispersion amorphous (LDA) state possessing modest extension of atomic necklaces. Further isothermal compression up to the order of 10 kbar to 1 Mbar or an isobaric cooling of LDA induced the transition to the high-dispersion amorphous (HDA) state. The HDA was characterized by long quantum wavelengths of atoms extended over several Angstroms and the promotion of atomic residual diffusion. They were related to the quantum tunneling of atoms bestriding the potential saddle points in this glass. The change in pressure or temperature induced the LDA-HDA transition reversibly with hysteresis, while it resembled the coil-globule transition of classical polymers. The HDA had lower kinetic and higher Gibbs free energies than the LDA at close temperature. The HDA was absent at T ≥ 13 K, while the LDA-HDA transition pressure significantly decreased with lowering temperature. The LDA and HDA correspond to the trapped and tunneling regimes proposed by Markland et al. [J. Chem. Phys. 136, 074511 (2012)], respectively. The same reentrant behavior as they found was observed for the expansion factor of the quantum wavelength as well as for atomic diffusivity.

Glass polyamorphism in gallium: Two amorphous solid states and their transformation on the potential energy landscape

Using the potential energy landscape (PEL) formalism and molecular dynamics simulations, we investigate a phase transformation between two amorphous solid states of gallium, namely, a low-density amorphous solid (LDA) and a high-density amorphous solid (HDA), and compare with its equilibrium counterpart, the liquid-liquid phase transition (LLPT). It is found that on the PEL, the signatures of the out-of-equilibrium LDA-HDA transition are reminiscent of those of the equilibrium LLPT in terms of pressure, inherent structure pressure, inherent structure energy, and shape function, indicating that the LDA-HDA transformation is a first-order-like transition. However, differences are also found between the out-of-equilibrium phase transition and the equilibrium one, for example, the path from LDA to HDA on the PEL cannot be accessed by the path from LDL to HDL. Our results also suggest that the signatures of the out-of-equilibrium transition in gallium are rather general features of systems with an accessible LLPT-not only systems with pairwise interactions but also those with many-body interactions. This finding is of crucial importance for obtaining a deeper understanding of the nature of transitions in the polyamorphic family.

Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water

We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.

Medium-range order in amorphous ices revealed by persistent homology

Despite the amorphous nature of glassy water, x-ray or neutron scattering experiments reveal sharp peaks in the structure factor, indicating the existence of medium-range order (MRO) in the system. However the real space origin of the peaks has yet to be disclosed. Herein, we use a combined approach of molecular dynamics simulations and persistent homology (PH) to investigate two types of glassy water, low-density amorphous (LDA) and high-density amorphous (HDA) ices. We present prominent MRO ring structures in each type of the ices, distinguished by their size and shape as well as the number of their components: MRO rings in HDA are observed smaller, less planar and more membered, compared to those in LDA. The PH-extracted MRO rings successfully reproduce the quantitative features, including the position and width, of the first sharp diffraction peaks in the structure factor, hence suitably serving as the origin of experimental MRO signatures in the amorphous ices. Our study supports that PH is an effective tool to identify hidden MRO in amorphous configurations.

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.

Glass polymorphism in TIP4P/2005 water: A description based on the potential energy landscape formalism

The potential energy landscape (PEL) formalism is a statistical mechanical approach to describe supercooled liquids and glasses. Here, we use the PEL formalism to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) using computer simulations of the TIP4P/2005 molecular model of water. We find that the properties of the PEL sampled by the system during the LDA-HDA transformation exhibit anomalous behavior. In particular, at conditions where the change in density during the LDA-HDA transformation is approximately discontinuous, reminiscent of a first-order phase transition, we find that (i) the inherent structure (IS) energy, e_IS(V), is a concave function of the volume and (ii) the IS pressure, P_IS(V), exhibits a van der Waals-like loop. In addition, the curvature of the PEL at the IS is anomalous, a nonmonotonic function of V. In agreement with previous studies, our work suggests that conditions (i) and (ii) are necessary (but not sufficient) signatures of the PEL for the LDA-HDA transformation to be reminiscent of a first-order phase transition. We also find that one can identify two different regions of the PEL, one associated with LDA and another with HDA. Our computer simulations are performed using a wide range of compression/decompression and cooling rates. In particular, our slowest cooling rate (0.01 K/ns) is within the experimental rates employed in hyperquenching experiments to produce LDA. Interestingly, the LDA-HDA transformation pressure that we obtain at T = 80 K and at different rates extrapolates remarkably well to the corresponding experimental pressure.

Temperature Effects on Water-Mediated Interactions at the Nanoscale

We perform molecular dynamics simulations to study the effects of temperature on the water-mediated interactions between nanoscale apolar solutes. Specifically, we calculate the potential of mean force (PMF) between two graphene plates immersed in water at 240 ≤ T ≤ 400 K and P = 0.1 MPa. These are thermodynamic conditions relevant to cooling- and heating-induced protein denaturation. It is found that both cooling and heating tend to suppress the attraction, and ultimate collapse, of the graphene plates. However, the underlying role played by water upon heating and cooling is different. Isobaric heating reduces the strength and range of the interactions between the plates. Instead, isobaric cooling stabilizes the plates separations that can accommodate an integer number of water layers between the graphene plates. In particular, the energy barriers separating these plate separations increase linearly with 1/ T. We also explore the sensitivity of the plates PMF to the water model employed. In the case of the TIP4P/2005 model, water confined between the plates crystallizes into a defective bilayer ice at low temperatures, whereas in the case of the SPC/E model, water remains in the liquid state at same thermodynamic conditions. The effects of varying water-graphene interactions on the plates PMF are also studied.

Heating- and pressure-induced transformations in amorphous and hexagonal ice: A computer simulation study using the TIP4P/2005 model

We characterize the phase behavior of glassy water by performing extensive out-of-equilibrium molecular dynamics simulations using the TIP4P/2005 water model. Specifically, we study (i) the pressure-induced transformations between low-density (LDA) and high-density amorphous ice (HDA), (ii) the pressure-induced amorphization (PIA) of hexagonal ice (I_h), (iii) the heating-induced LDA-to-HDA transformation at high pressures, (iv) the heating-induced HDA-to-LDA transformation at low and negative pressures, (v) the glass transition temperatures of LDA and HDA as a function of pressure, and (vi) the limit of stability of LDA upon isobaric heating and isothermal decompression (at negative pressures). These transformations are studied systematically, over a wide range of temperatures and pressures, allowing us to construct a P-T phase diagram for glassy TIP4P/2005 water. Our results are in qualitative agreement with experimental observations and with the P-T phase diagram obtained for glassy ST2 water that exhibits a liquid-liquid phase transition and critical point. We also discuss the mechanism for PIA of ice I_h and show that this is a two-step process where first, the hydrogen-bond network (HBN) is distorted and then the HBN abruptly collapses. Remarkably, the collapse of the HB in ice I_h occurs when the average molecular orientations order, a measure of the tetrahedrality of the HBN, is of the same order as in LDA, suggesting a common mechanism for the LDA-to-HDA and I_h-to-HDA transformations.

On the time required to freeze water

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.

Electro-nucleation of water nano-droplets in No Man's Land to fault-free ice Ic

Elucidating water-to-ice freezing, especially in "No Man's Land" (150 K < T < 235 K), is fundamentally important (e.g., predicting upper-troposphere cirrus-cloud formation) - and elusive. An oft-neglected aspect of tropospheric ice-crystallite formation lies in inevitably-present electric fields' role. Exploring nucleation in No Man's Land is technically demanding, owing to rapid nucleation rates, to mention nothing of difficulties of applying relevant electric fields thereto. Here, we tackle these intriguing open questions, via non-equilibrium molecular-dynamics simulations of sub-microsecond formation of rhombus-shaped ice I_c nano-crystallites from aggressively-quenched supercooled water nano-droplets in the gas phase, in external static electric fields. We explore droplets' nano-confined geometries and the entropic-ordering agent of external electric fields as a means of realising cubic-ice formation, especially with very few stacking faults and defects.

Large-Scale Structure and Hyperuniformity of Amorphous Ices

We investigate the large-scale structure of amorphous ices and transitions between their different forms by quantifying their large-scale density fluctuations. Specifically, we simulate the isothermal compression of low-density amorphous ice (LDA) and hexagonal ice to produce high-density amorphous ice (HDA). Both HDA and LDA are nearly hyperuniform; i.e., they are characterized by an anomalous suppression of large-scale density fluctuations. By contrast, in correspondence with the nonequilibrium phase transitions to HDA, the presence of structural heterogeneities strongly suppresses the hyperuniformity and the system becomes hyposurficial (devoid of "surface-area fluctuations"). Our investigation challenges the largely accepted "frozen-liquid" picture, which views glasses as structurally arrested liquids. Beyond implications for water, our findings enrich our understanding of pressure-induced structural transformations in glasses.

Glass polymorphism in glycerol-water mixtures: I. A computer simulation study

We perform out-of-equilibrium molecular dynamics (MD) simulations of water-glycerol mixtures in the glass state. Specifically, we study the transformations between low-density (LDA) and high-density amorphous (HDA) forms of these mixtures induced by compression/decompression at constant temperature. Our MD simulations reproduce qualitatively the density changes observed in experiments. Specifically, the LDA-HDA transformation becomes (i) smoother and (ii) the hysteresis in a compression/decompression cycle decreases as T and/or glycerol content increase. This is surprising given the fast compression/decompression rates (relative to experiments) accessible in MD simulations. We study mixtures with glycerol molar concentration χ(g) = 0-13% and find that, for the present mixture models and rates, the LDA-HDA transformation is detectable up to χ(g) ≈ 5%. As the concentration increases, the density of the starting glass (i.e., LDA at approximately χ(g) ≤ 5%) rapidly increases while, instead, the density of HDA remains practically constant. Accordingly, the LDA state and hence glass polymorphism become inaccessible for glassy mixtures with approximately χ(g) > 5%. We present an analysis of the molecular-level changes underlying the LDA-HDA transformation. As observed in pure glassy water, during the LDA-to-HDA transformation, water molecules within the mixture approach each other, moving from the second to the first hydration shell and filling the first interstitial shell of water molecules. Interestingly, similar changes also occur around glycerol OH groups. It follows that glycerol OH groups contribute to the density increase during the LDA-HDA transformation. An analysis of the hydrogen bond (HB)-network of the mixtures shows that the LDA-HDA transformation is accompanied by minor changes in the number of HBs of water and glycerol. Instead, large changes in glycerol and water coordination numbers occur. We also perform a detailed analysis of the effects that the glycerol force field (FF) has on our results. By comparing MD simulations using two different glycerol models, we find that glycerol conformations indeed depend on the FF employed. Yet, the thermodynamic and microscopic mechanisms accompanying the LDA-HDA transformation and hence, our main results, do not. This work is accompanied by an experimental report where we study the glass polymorphism in glycerol-water mixtures prepared by isobaric cooling at 1 bar.