Christina M.
Tonauer
,
Marion
Bauer‡
and
Thomas
Loerting
*
Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, Innsbruck A-6020, Austria. E-mail: thomas.loerting@uibk.ac.at
First published on 6th December 2021
For slowly compressed hexagonal ice pressure-induced amorphisation to high-density amorphous ice (HDA) takes place below and at 130 K, but polymorphic transformation to ice IX takes place at 140–170 K. Stable ice II only forms above 170 K. Ice IX impurities trigger ice IX growth even at 120 K. HDA and ice IX are equally long-lived, where both can be regarded as metastable phases.
We here take a stand on the controversy emerging from these findings, where they challenge the idea that HDA is a glassy form thermodynamically related to the liquid. Their results contradict the results obtained in the pioneering work by Mishima et al.1,10 that were later confirmed by many groups, including Klotz et al.,11 Strässle et al.,12 Gromnitskaya et al.13 and our own group.14 While failing to reproduce the results in these literature studies, the results of Tulk et al. are quite similar to results in our own group obtained at much higher compression temperature.15,16 Specifically, we have shown in M. Bauer et al. that compression rate as well as compression temperature are of crucial importance to determine the outcome of compression experiments of ice Ih.15,16 In that work we have shown that ice III/IX forms as the major product quite easily at 170 K, but ice II predominates at higher temperature. Faster compression favours the metastable form ice III/IX, and slower compression favours the thermodynamically stable form ice II. In order to obtain thermodynamically stable, pure ice II upon compression of ice Ih, the compression needs to be slower than 2 MPa min−1 at 170 K, but can be as fast as 100 MPa min−1 at 220 K (see Fig. 4 in ref. 16). In other words, metastable ice III/IX is favoured at lower temperatures and higher compression rates over stable ice II. Exactly the same principle applies for the formation of HDA vs. ice III/IX, where in this case HDA is metastable with respect to ice III/IX, which itself is metastable with respect to stable ice II. Throughout this manuscript, we define the term “metastable” in the sense that relaxation times are shorter than transformation times by orders of magnitude. This implies that equilibration is achieved prior to transformation. Based on the dielectric relaxation times measured by Lemke et al.17 this criterion is fulfilled for HDA, where dielectric relaxation times are approximately ten thousand times shorter than the polyamorphic transition to LDA. The metastable HDA only forms at lower temperature whereas ice III/IX (or even ice II) forms at higher temperature. The compression rate dependence and the phases emerging found by Tulk et al.9 are exactly the same as that reported by us in ref. 15 and 16, but at 100 K instead of 170 K.
R. Bauer et al.18 repeated the experiment by Tulk et al.9 at a different high-pressure neutron diffraction setup (J-PARC MLF), where they find that HDA indeed does form at 100 K in slow compression experiments, not ice III/IX. In spite of the contradicting results, R. Bauer et al. come up with an interpretation similar to the one of Tulk et al., i.e., formation of HDA is the result of an incomplete transition between two crystalline forms of ice, that is, ice Ih and ice VI/XV.18 In order to clarify what the experimental findings actually are and how to understand the nature of HDA we conducted slow compression experiments following their protocol, but with direct measurement of pressure and temperature in our piston-cylinder setup and using powder X-ray diffraction for sample characterisation. This extends our earlier work done at 170–230 K15,16 to lower temperature. Furthermore, we have purposefully contaminated ice Ih samples with ice III/IX to study the impact of these impurities on the outcome of the slow compression experiment.
On the methodological side, the slow compression experiments were done following the procedures established over the last 20 years in our group.19,20 400 μL of ultrapure liquid H2O were pipetted into a precooled container made of ∼0.7 g indium foil. After crystallisation of water to hexagonal ice, the sample container was inserted into the 8 mm bore of our custom-made high-pressure piston cylinder setup that is pressurised using a commercial “universal material testing machine” (Zwick, model BZ100/TL3S). The high-pressure cell was cooled by liquid nitrogen flowing through a copper loop around the cylinder. Heating elements as well as a Pt-100 temperature sensor were inserted in additional bores in the cell, 1 cm from the sample position (see Fig. 3.1 in ref. 21 or Fig. 10 in ref. 22 for a schematic drawing of our cell/setup). This apparatus enables the simultaneous detection and control of volume change, temperature and pressure in real time. The samples were pre-compressed (except for the sample compressed at 4 MPa min−1 shown in ESI,† Fig. S1a and S2a) to 0.4 GPa at 100 MPa min−1 and 77 K in order to remove air trapped within the ice. Then the sample was decompressed to 0.02 GPa at 20 MPa min−1 and heated to the desired temperature. Ice IX impurities were introduced on purpose into some ice Ih samples as follows: hexagonal ice was isobarically heated to 170 K at 0.02 GPa and then compressed at 100 MPa min−1. At ∼0.32 GPa, in the process of transformation to IX, the sample was quenched by pouring liquid nitrogen over the cell. Depending on the progress of ice IX formation until quenching, the amount of ice IX impurity in a mixed sample could be up to ∼50%. After reaching the end pressure at a certain temperature, samples were quench-recovered and analysed using X-ray diffraction (Cu-Kα) at temperatures <80 K and subambient pressure.
Fig. 1a shows volume curves upon compression of ice Ih at 100 K at ∼4 MPa min−1, which corresponds to the linear equivalent of the stepwise compression profile used by Tulk et al.9 These curves exhibit linear densification up to 0.94 GPa when a steplike densification of about 20% occurs. X-Ray diffractograms (see Fig. 1b) reveal that ice Ih is stable up to 0.8 GPa, whereas Tulk et al. have already observed ice IX at 0.3 GPa. At 1.0 GPa, first signs of the HDA halo peak (maximum at d = 2.9 Å) become apparent (emerging halo peak marked by vertical arrow), and at 1.4 GPa only the halo peak is observed, with no signs of any Bragg peaks. That is, pressure-induced amorphisation is complete and HDA easily survives for hours at p > 1 GPa at 100 K, whereas ice IX does not appear at all.
In order to mimic the experiment by Tulk et al. even more closely, we have also used a stepwise compression protocol at 100 K. For this purpose, pressure was increased using a step-profile, using 50 MPa min−1 for each step, and maintaining the plateaus for 1 hour (see incrementally coloured ΔV(p) line in ESI,† Fig. S1a). No matter whether we use the linear ramp or the stepwise protocol (see ESI,† Fig. S1b), HDA is produced, not crystalline ices. X-Ray diffractograms of the resulting samples feature the broad halo peak typical of HDA, but no sharp Bragg peaks typical of crystalline ice (shown as insets in ESI,† Fig. S1a). Even at half the compression rate used by Tulk et al. (2 MPa min−1) we do not find any signs of crystalline ice IX at 100 K. The respective (a) ΔV(p) and (b) p(t) lines (red) are compared with the ones of an experiment conducted at 4 MPa min−1 (black) in ESI,† Fig. S2. Again, the X-ray diffractograms (inset in ESI,† Fig. S2a) exhibit the formation of only HDA (with a minor amount of ice Ih). Even though this experiment takes about 10 hours, HDA survives. There are no signs of any crystalline high-pressure ice phases after such a long experiment. This demonstrates the longevity of HDA and implies that at 100 K the activation barrier is high enough and the prefactor small enough so that none of the attempts to cross the barrier for HDA crystallisation is successful. There might be an even slower compression experiment that allows formation of metastable ice III/IX (or even stable ice II) at 100 K from ice Ih. However, such an experiment would certainly exceed the limit of available time resources for a PhD thesis. We surmise that the ice Ih sample would need to be kept for years at 100 K to actually reach conversion to an ice that is thermodynamically more stable than HDA. Put in other words, HDA is an ice form of longevity at 100 K, maybe even an infinitely stable phase.
The situation changes at 140 K (see Fig. 1c and d): instead of amorphisation at 1.0 GPa we observe crystallisation to ice IX at 0.65 GPa (Fig. 1c). The steplike densification occurs at much lower pressure at 140 K. In the middle of the densification, at 0.7 GPa (red ΔV(p) line in Fig. 1c), we observe Bragg peaks both from ice Ih and ice IX, indicating that the transition is in progress, but incomplete at 0.7 GPa (see red middle diffractogram in Fig. 1d). At 1.2 GPa, far above the offset of densification (orange ΔV(p) line in Fig. 1c), the diffractograms no longer show any signs of a broad halo peak, but only ice IX Bragg peaks (see orange bottom diffractogram in Fig. 1d). That is, at this elevated temperature we observe the transition Ih → IX at 4 MPa min−1. Yet, the onset pressure of the transformation (0.65 GPa) is significantly higher compared to the one reported by Tulk et al. (0.3 GPa). Finally, at 170 K (see Fig. 1e and f) we observe a transition starting at ∼0.3 GPa. For a compression rate of 5 MPa min−1 (purple line in Fig. 1e) stable ice II and some metastable ice IX form (see upper diffractogram in Fig. 1f). Only for a rate of 2 MPa min−1 (pink line in Fig. 1e) at 170 K thermodynamics ultimately win, and the stable ice II is formed (see lower diffractogram in Fig. 1f). Again, our results at 170 K are quite similar to the ones by Tulk et al. at 100 K, i.e., there is a huge offset of 70 K.
That is, sample temperature and its measurement could be one source of the discrepancy. However, a stable temperature gradient of ∼70 K between sample and sensor, unnoticed by Tulk et al. seems unrealistic. In the interest of full transparency, we here show our temperature sensor data for all the experiments considered in the present study in ESI,† Fig. S3. In general the temperature remains constant to better than ±0.3 K for many hours. At some points in the experiment the flow through the cooling loops suddenly changed (e.g., through dirt or ice in the lines). In such situations the temperature fluctuations increase, and a bit larger deviations from the set temperature appear for a few minutes, after which the set temperature is reached again (e.g., yellow line for set temperature 130 K in ESI,† Fig. S3b). Since the temperature sensor is very close to the sample and remains constant for many hours we regard the sensor temperature to exactly reflect the sample temperature, so that we trust our temperatures to better than ±1 K.
In order to rule out, that an isotope effect (H2O ↔ D2O) is at the origin of the difference between the results obtained by Tulk et al. and ourselves, we conducted a control experiment compressing ∼400 μL D2O ice Ih at 100 K at a rate of 4 MPa min−1 (see ESI,† Fig. S4). As expected, we observed the same transformation as for H2O, i.e., pressure-induced amorphisation. The main difference is a slightly higher onset pressure for the transformation. Such small isotope effects are typical for phase-transformations, in which the oxygen-lattice is changing, such as pressure-amorphisation.
An actual source for the deviating results of Tulk et al. might be the presence of nano-sized ice IX seeds in their ice Ih samples, too small to be detected by neutron diffraction. They mention precautions to avoid contamination with other ice phases upon closing their cell.9 Yet, slight fluctuations of gas pressure during the increase of load on the Paris–Edinburgh press could induce abrupt changes of the compression rate during the experiment, causing uncontrolled shock-wave heating.23 Such a sudden heating event, taking place within a millisecond, could remain unnoticed. Shock-wave heating is avoided in our piston-cylinder experiments by the use of indium as a low-temperature lubricant and controlled build-up of pressure from the first moment on. This concept goes back to the work of Mishima.1 In the absence of indium a loud bang and a massive drop in pressure indicate the release of the shock-wave. Compression experiments of ice Ih inflicted with shock-wave heating were conducted by Kohl et al.24 They have shown that such shock-waves may raise the sample temperature by up to 100 K, causing formation of ice XII, provided some initial HDA has formed through PIA at >0.9 GPa. Shock-wave heating of ice Ih at lower pressure might produce other ice phases, e.g., ice IX at <0.2 GPa.
For this reason, we studied the influence of ice IX contamination by slowly compressing mixtures of Ih/IX at 100–140 K (see Fig. 2). Depending on temperature, the volume curves in Fig. 2a and the X-ray diffractograms in Fig. 2b show three distinct outcomes: below 120 K, the ΔV(p) curves reveal only a single Ih → HDA densification step at ∼1 GPa (purple and blue curves in Fig. 2a). The presence of ice IX contamination barely affects pressure-amorphisation; compression curves are almost identical to the pure ice Ih case (see Fig. 1a). There is only a very small change of onset pressure when comparing compression of pure ice Ih and an ice Ih/IX mixture (see arrow directed upwards in Fig. 2c at 100 K). At 120–130 K both transitions, Ih → IX and Ih → HDA, occur well separated in terms of pressure. The compression curves (light-blue, green and yellow in Fig. 2a) show two kinks, one near 0.4 GPa, the other near 0.9 GPa. The first kink indicates the onset of ice IX formation, the second kink indicates formation of HDA from the remainder of ice Ih. That is, there is a competition between amorphisation and polymorphic transformation under these conditions. The horizontal dotted lines in Fig. 2a indicate the volume change incurred through ice IX formation from ice Ih (ΔVIh→IX) as well as the volume change incurred through HDA formation (ΔVIh→HDA). Fig. 2d shows that the HDA-fraction is 100% below 120 K and drops to about 50% at 130 K. The fraction of ice IX increases with temperature, at the expense of HDA. This is also evident in the X-ray diffractograms in Fig. 2b (coloured traces in the middle). It is clearly seen that the HDA halo marked by an arrow becomes less and less pronounced with higher compression temperature. Above 130 K there is only ice IX, but no HDA. At 135 K and 140 K (orange and red curves in Fig. 2a) there is a single densification step again, where the transition is Ih → IX now. When comparing the onset pressure of ice IX growth at 140 K between compression of pure ice Ih (Fig. 1e) and compression of ice Ih/IX (Fig. 2a, red curve) a marked difference is evident. Ice IX impurities act as seeds for ice IX formation and lower the onset pressure from 0.65 GPa to 0.40 GPa at 140 K (downward-arrow in Fig. 2c). This has implications for our understanding of the issues in the experiment by Tulk et al. on PIA of ice Ih: a short-term temperature gradient between sensor and sample of about 30–40 K and some nano-scaled ice IX contamination (undetected by neutron diffraction) as a result of shock-wave heating could explain the incompatibility with the experiments in other groups.
![]() | ||
Fig. 2 (a) Volume curves ΔV(p) obtained upon compression of mixtures of ices Ih/IX at temperatures 100–140 K and a constant compression rate of 4 MPa min−1. (b) Powder X-ray diffractograms of quench-recovered samples after compression to 1.2 GPa at 4 GPa min−1. The top diffractogram represents the starting material for this series of experiments, that is, mixtures of Ih/IX. Diffraction angles at top x-axis correspond to Cu-Kα radiation. (c) Onset pressures for the transitions Ih → IX and Ih → HDA in samples of Ih/IX as a function of temperature ponset(T). The reproducibility is better than 0.01 GPa. Onset pressures of the transitions of pure ice Ih at 100 K and 140 K (extracted from Fig. 1a and c) are shown (blue and red circles) for comparison. The onset pressures are defined as the intersection point between two linear functions fitting the linear section of the ΔV(p) curves in (a) just before and during the transition. (d) Evaluation of the densification step of the volume curves ΔV(p) in (a). At low and high temperatures, either Ih → HDA or Ih → IX is operative. The total change of volume ΔVtotal in the intermediate temperature region is composed of the volume changes of both transition modes, ΔVIh→IX and ΔVIh→HDA. |
![]() | ||
Fig. 3 Internal energy sketches for the scenarios proposed by (a) us and (b) Tulk et al.9 |
This view is inconsistent with all experimental findings on PIA of ice Ih, except for the work from Tulk et al., which we suggest to be afflicted with issues of temperature gradients and ice IX contamination based on our findings. This sets the experimental record straight and confirms earlier work on pressure-induced amorphisation of ice Ih. Furthermore, we point out that Tulk et al. themselves observe a metastable ice, namely ice IX, rather than the stable ice II at 0.3 GPa. This means that they do not reach the thermodynamic limit and the equilibrium form of ice, but instead an ice that is kinetically stabilised. In their case the barrier of polymorphic conversion from ice IX to ice II is high enough to prevent the most stable ice form (see transition state between ice IX and ice II in Fig. 3, marked by black-dotted line). Rather than observing stable ice II they observe that metastable ice IX persists. Our earlier work, M. Bauer et al.,15 shows that a temperature of 230 K or above is necessary to overcome this barrier for the polymorphic transition and to reach ice II in the slow compression limit. Just like ice IX at 0.3 GPa also HDA is kinetically stable at >1 GPa, where 140 K should not be exceeded to avoid HDA crystallisation on the time scale of hours. Finally, Tulk et al. have not accounted for the huge enhancement of HDA lifetime after its annealing above 1 GPa.27,28 Annealing removes nano-sized crystalline domains serving as growth sites for crystalline ice and thus enhances stability.29 What Tulk et al. produce by PIA at 100 K is referred to as unannealed HDA (uHDA),27 while only the much more stable expanded HDA (eHDA) was found to be directly related to ultraviscous water.25 This leads to the conclusion that the activation barriers for HDA crystallisation are significantly higher than claimed by Tulk et al.9 as well as R. Bauer et al.18 In fact, they are high enough to keep HDA stable above its glass-to-liquid transition temperature, so that high-density liquid water (HDL) may equilibrate instead of crystallisation to ice IV and ice XII. This has been demonstrated in the recent work of our group by Stern & Loerting30 as well as Handle & Loerting.31
While the recent experimental findings presented by R. Bauer et al.18 are consistent with the results presented in Fig. 1, there is a significant difference about the interpretation of the nature of HDA. They claim HDA to be a frustrated state between ice Ih and ice VI/XV, contrasting our claim above of HDA being a metastable phase related to the high-pressure liquid state. Similarly Shephard et al. claim HDA to be a derailed state along the ice Ih to ice IV conversion.32 The resistance of HDA against crystallisation for many hours at 120 K speaks against both ideas of HDA being a transient state in a polymorphic conversion. In slow heating experiments HDA even resists crystallisation up to 161 K.30 Also the conjecture by R. Bauer et al.18 that amorphisation represents the onset from a single H-bond network to an interpenetrating H-bond network defies experimental evidence. In reality, HDA always crystallises first to the single H-bond network ice IV and/or ice XII,33,34 but not directly to ice VI/XV with interpenetrating H-bond networks. In fact, the transition to interpenetrating networks does not require an amorphous intermediate, but directly takes place in the polymorphic ice IV → ice VI or XII → ice VI transition.
Taken together this experimental evidence allows regarding HDA/HDL as a “metastable phase” in the thermodynamic sense rather than a “transient state”, in spite of its glassy nature. Unlike “simple amorphous states”, metastable HDA/HDL can be reached through many different thermodynamic paths, where pressure-induced amorphisation studied in the present work is only one of them. Other ways include cooling of the pressurised liquid,35 temperature-induced amorphisation of ice,36 decompression-induced melting of high-pressure ice,37 or compression of hyperquenched water,38 to name just a few. This is only possible if HDA/HDL is describable as a deep basin (as sketched in Fig. 3a), where many paths lead to the same, local minimum on the surface.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cp03922a |
‡ Current address: Sandoz GmbH, Biochemiestr. 10, A-6250 Kundl, Austria. |
This journal is © the Owner Societies 2022 |