Search

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).

We perform path-integral molecular dynamics (PIMD) simulations of H 2 O and D 2 O 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 H 2 O and D 2 O 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 C P (T), and static dielectric constant ε(T). While these properties are in excellent [α P (T) and κ T (T)] or semi-quantitative agreement [C P (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 C P (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 C P (T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H 2 O and D 2 O 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 C P (T), D(T), and structural properties being the only (expected) exceptions.

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.

Polyamorphic substances have the ability to exist in more than one liquid and/or glass states. Examples include water, silicon, and hydrogen. In many of these substances, nuclear quantum effects may become important in the proximity of the liquid-liquid and glass-glass transformation. Here, we study the nuclear quantum effects on a monatomic liquid that exhibits water-like anomalous properties and a liquid-liquid phase transition (LLPT) ending at a liquid-liquid critical point (LLCP). By performing path integral Monte Carlo simulations with different values of the Planck's constant h, we are able to explore how the location of the LLCP/LLPT in the P-T plane shifts as the system evolves from classical, h = 0, to quantum, h > 0. We find that, as the quantum nature of the liquid (as quantified by h) increases, and the atoms in the liquid become more delocalized, the LLCP pressure increases, the LLCP temperature decreases, and the LLCP volume remains constant. In addition, the crystallization temperature decreases with increasing h. For large values of h, the LLCP is not accessible due to rapid crystallization. The structure of the liquids studied at different values of h are also investigated.

We report a detailed experimental study of (i) pressure-induced transformations in glycerol-water mixtures at T = 77 K and P = 0-1.8 GPa, and (ii) heating-induced transformations of glycerol-water mixtures recovered at 1 atm and T = 77 K. Our samples are prepared by cooling the solutions at ambient pressure at various cooling rates (100 K s(-1)-10 K h(-1)) and for the whole range of glycerol mole fractions, χ(g). Depending on concentration and cooling rates, cooling leads to samples containing amorphous ice (χg ≥ 0.20), ice (χ(g) ≤ 0.32), and/or "distorted ice" (0 < χ(g) ≤ 0.38). Upon compression, we find that (a) fully vitrified samples at χ(g) ≥ 0.20 do not show glass polymorphism, in agreement with previous works; (b) samples containing ice show pressure-induced amorphization (PIA) leading to the formation of high-density amorphous ice (HDA). PIA of ice domains within the glycerol-water mixtures is shown to be possible only up to χ(g) ≈ 0.32 (T = 77 K). This is rather surprising since it has been known that at χ(g) < 0.38, cooling leads to phase-separated samples with ice and maximally freeze-concentrated solution of χ(g) ≈ 0.38. Accordingly, in the range 0.32 < χ(g) < 0.38, we suggest that the water domains freeze into an interfacial ice, i.e., a highly-distorted form of layered ice, which is unable to transform to HDA upon compression. Upon heating samples recovered at 1 atm, we observe a rich phase behavior. Differential scanning calorimetry indicates that only at χ(g) ≤ 0.15, the water domains within the sample exhibit polyamorphism, i.e., the HDA-to-LDA (low-density amorphous ice) transformation. At 0.15 < χ(g) ≤ 0.38, samples contain ice, interfacial ice, and/or HDA domains. All samples (χ(g) ≤ 0.38) show: the crystallization of amorphous ice domains, followed by the glass transition of the vitrified glycerol-water domains and, finally, the melting of ice at high temperatures. Our work exemplifies the complex "phase" behavior of glassy binary mixtures due to phase-separation (ice formation) and polyamorphism, and the relevance of sample preparation, concentration as well as cooling rates. The presence of the distorted ice (called "interphase" by us) also explains the debated "drift anomaly" upon melting. These results are compatible with the high-pressure study by Suzuki and Mishima indicating disappearance of polyamorphism at P ≈ 0.03-0.05 GPa at χ(g) ≈ 0.12-0.15 [J. Chem. Phys., 2014, 141, 094505].

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.

We measure the solvation free energy, Δμ*, for hard spheres and Lennard-Jones particles in a number of artificial liquids made from modified water models. These liquids have reduced hydrogen bond strengths or altered bond angles. By measuring Δμ* for a number of state points at P = 1 bar and different temperatures, we obtain solvation entropies and enthalpies, which are related to the temperature dependence of the solubilities. By resolving the solvation entropy into the sum of the direct solute-solvent interaction and a term depending on the solvent reorganisation enthalpy we show that, although the hydrophobic effect in water at 300 K arises mainly from the small molecular size, its temperature dependence is anomalously low because the reorganisation enthalpy of liquid water is unusually small. We attribute this to the strong tetrahedral network which results from both the molecular geometry and the hydrogen bond strength.