
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39313509
52505
10.1038/s41467-024-52505-0
Article
Phase transition kinetics of superionic H2O ice phases revealed by Megahertz X-ray free-electron laser-heating experiments
http://orcid.org/0000-0002-7666-401X
Husband R. J. rachel.husband@desy.de

1
http://orcid.org/0000-0001-5039-1183
Liermann H. P. 1
http://orcid.org/0000-0002-2630-8092
McHardy J. D. 2
http://orcid.org/0000-0002-3730-8661
McWilliams R. S. 2
http://orcid.org/0000-0002-6422-8819
Goncharov A. F. 3
http://orcid.org/0000-0001-9270-2330
Prakapenka V. B. 4
http://orcid.org/0000-0003-4363-7434
Edmund E. 3
Chariton S. 4
Konôpková Z. 5
Strohm C. 1
http://orcid.org/0000-0001-5046-1612
Sanchez-Valle C. 6
http://orcid.org/0000-0001-6879-0422
Frost M. 7
Andriambariarijaona L. 8
http://orcid.org/0000-0002-2902-2102
Appel K. 5
http://orcid.org/0000-0003-1480-511X
Baehtz C. 9
http://orcid.org/0000-0002-5215-0153
Ball O. B. 2
Briggs R. 10
http://orcid.org/0000-0001-5671-5214
Buchen J. 1115
http://orcid.org/0000-0002-2808-2963
Cerantola V. 516
http://orcid.org/0000-0002-9143-4810
Choi J. 12
http://orcid.org/0000-0002-5692-4400
Coleman A. L. 10
http://orcid.org/0000-0003-4658-5764
Cynn H. 10
Dwivedi A. 5
http://orcid.org/0000-0003-2304-667X
Graafsma H. 1
http://orcid.org/0000-0002-8498-3811
Hwang H. 117
Koemets E. 1118
http://orcid.org/0000-0002-2258-2123
Laurus T. 1
http://orcid.org/0000-0002-2043-0804
Lee Y. 12
Li X. 119
http://orcid.org/0000-0003-1784-6515
Marquardt H. 11
http://orcid.org/0000-0002-3424-5693
Mondal A. 6
http://orcid.org/0000-0003-0868-4745
Nakatsutsumi M. 5
Ninet S. 8
http://orcid.org/0000-0002-7328-1099
Pace E. 2
http://orcid.org/0000-0002-9638-3303
Pepin C. 13
Prescher C. 14
Stern S. 120
Sztuk-Dambietz J. 5
http://orcid.org/0000-0002-3575-4449
Zastrau U. 5
McMahon M. I. M.I.McMahon@ed.ac.uk

2
1 https://ror.org/01js2sh04 grid.7683.a 0000 0004 0492 0453 Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
2 https://ror.org/01nrxwf90 grid.4305.2 0000 0004 1936 7988 SUPA, School of Physics and Astronomy, and Centre for Science at Extreme Conditions, The University of Edinburgh, Edinburgh, UK
3 grid.531655.1 Carnegie Science, Earth and Planets Laboratory, Washington, DC USA
4 https://ror.org/024mw5h28 grid.170205.1 0000 0004 1936 7822 The University of Chicago, Center for Advanced Radiation Sources, Chicago, IL USA
5 https://ror.org/01wp2jz98 grid.434729.f 0000 0004 0590 2900 European XFEL, Schenefeld, Germany
6 https://ror.org/00pd74e08 grid.5949.1 0000 0001 2172 9288 Universität Münster, Institut für Mineralogie, Corrensstraße 24, Münster, Germany
7 https://ror.org/05gzmn429 grid.445003.6 0000 0001 0725 7771 SLAC National Accelerator Laboratory, California, USA
8 https://ror.org/02en5vm52 grid.462844.8 0000 0001 2308 1657 Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Sorbonne Université, Paris, France
9 https://ror.org/01zy2cs03 grid.40602.30 0000 0001 2158 0612 Institute of Radiation Physics, Helmholtz-Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, Dresden, Germany
10 https://ror.org/041nk4h53 grid.250008.f 0000 0001 2160 9702 Lawrence Livermore National Laboratory, Livermore, CA USA
11 https://ror.org/052gg0110 grid.4991.5 0000 0004 1936 8948 Department of Earth Sciences, University of Oxford, Oxford, UK
12 https://ror.org/01wjejq96 grid.15444.30 0000 0004 0470 5454 Department of Earth System Sciences, Yonsei University, Seoul, Korea
13 https://ror.org/03xjwb503 grid.460789.4 0000 0004 4910 6535 CEA, DAM, DIF, 91297 Arpajon, France; Université Paris-Saclay, CEA, Laboratoire Matière en Conditions Extrêmes, Bruyères-le-Châtel, France
14 https://ror.org/0245cg223 grid.5963.9 0000 0004 0491 7203 Institute of Earth and Environmental Sciences, University of Freiburg, Freiburg, Germany
15 https://ror.org/0234wmv40 grid.7384.8 0000 0004 0467 6972 Present Address: Bayerisches Geoinstitut, Universität Bayreuth, Universitätsstraße 30, Bayreuth, Germany
16 grid.7563.7 0000 0001 2174 1754 Present Address: Department of Earth and Environmental Sciences (DISAT), University of Milano-Bicocca, Milan, Italy
17 https://ror.org/024kbgz78 grid.61221.36 0000 0001 1033 9831 Present Address: School of Earth Sciences and Environmental Engineering, Gwangju Institute of Science and Technology, Gwangju, Republic of Korea
18 https://ror.org/05etxs293 grid.18785.33 0000 0004 1764 0696 Present Address: Diamond Light Source Ltd, Harwell Science and Innovation Campus, Didcot, Oxfordshire UK
19 grid.64924.3d 0000 0004 1760 5735 Present Address: Synergetic Extreme Condition High-Pressure Science Center, State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun, China
20 grid.523979.3 Present Address: X-Spectrum GmbH, Luruper Hauptstraße 1, Hamburg, Germany
23 9 2024
23 9 2024
2024
15 82567 6 2024
11 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
H2O transforms to two forms of superionic (SI) ice at high pressures and temperatures, which contain highly mobile protons within a solid oxygen sublattice. Yet the stability field of both phases remains debated. Here, we present the results of an ultrafast X-ray heating study utilizing MHz pulse trains produced by the European X-ray Free Electron Laser to create high temperature states of H2O, which were probed using X-ray diffraction during dynamic cooling. We confirm an isostructural transition during heating in the 26-69 GPa range, consistent with the formation of SI-bcc. In contrast to prior work, SI-fcc was observed exclusively above ~50 GPa, despite evidence of melting at lower pressures. The absence of SI-fcc in lower pressure runs is attributed to short heating timescales and the pressure-temperature path induced by the pump-probe heating scheme in which H2O was heated above its melting temperature before the observation of quenched crystalline states, based on the earlier theoretical prediction that SI-bcc nucleates more readily from the fluid than SI-fcc. Our results may have implications for the stability of SI phases in ice-rich planets, for example during dynamic freezing, where the preferential crystallization of SI-bcc may result in distinct physical properties across mantle ice layers.

The authors perform heating experiments using femtosecond X-ray free electron laser pulses to explore the phase stability of superionic H2O. The absence of a face-centered cubic phase below 50 GPa, where superionic ice forms from the melt, is attributed to the short heating time and may help understanding the stability of superionic phases in ice-rich planets.

Subject terms

Planetary science
Phase transitions and critical phenomena
https://doi.org/10.13039/501100000766 Atomic Weapons Establishment (AWE) AWE CASE studentship P030463429 McHardy J. D. https://doi.org/10.13039/100000001 National Science Foundation (NSF) EAR-2049127 EAR - 1634415 EAR-2049127 EAR - 1634415 Goncharov A. F. Prakapenka V. B. Edmund E. Chariton S. https://doi.org/10.13039/100000015 U.S. Department of Energy (DOE) FWP100182 DE-AC52-07NA27344 Frost M. Cynn H. https://doi.org/10.13039/100010661 EC | Horizon 2020 Framework Programme (EU Framework Programme for Research and Innovation H2020) 864877 Marquardt H. https://doi.org/10.13039/100012946 Association Sorbonne Université (Sorbonne University association) Emergence HP-XFEL Ninet S. https://doi.org/10.13039/501100000266 RCUK | Engineering and Physical Sciences Research Council (EPSRC) EP/R02927X/1 McMahon M. I. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The high pressure (P) phase diagram of H2O is remarkably complex, containing numerous phases that differ in their crystal structures1–3, chemical bonding4–6, and degree of orientational7–9 and dynamical disorder6,10. In recent years, particular interest has been given to the high temperature (T), superionic (SI) forms of solid H2O (ice)11–15, which are characterized by the fast diffusion of protons within a crystalline oxygen sublattice. In particular, studies of the stability fields, structural behavior, and protonic conductivity of SI ice are motivated by the theorized presence of vast quantities of these phases in the interiors of Uranus and Neptune16,17, where the unusual properties of these SI ice phases have been proposed to contribute to the non-dipolar and non-axisymmetric magnetic fields of these planets18. Since its initial prediction19, SI ice has been the focus of numerous theoretical studies which proposed that it adopts a body-centered cubic (SI-bcc) oxygen sublattice at low pressures20–22, transforming to a face-centered (SI-fcc) lattice above ~100 GPa23–25. The SI-bcc oxygen lattice is isostructural with that of cubic ice VII, which is stable at lower temperatures, but the boundary between the two phases is first order in nature26. Experimental studies, however, are comparably scarce due to the challenges associated with creating and probing SI states using both static and dynamic compression methods. The stability field of SI ice cannot be reached by single shock compression of the ambient liquid due to the steep rise in temperature along the principal Hugoniot at low pressures, and the first experimental verification of SI ice was achieved using shock compression of pre-compressed samples11 and a reverberating shock compression scheme14,27,28 based on its high (total) conductivity (>100 S/cm)27,28 and low optical conductivity («100 S/cm)11 above 100 GPa. The crystal structure in the SI region (at 160 GPa and 3200 K) was first identified as fcc using nanosecond X-ray diffraction (XRD)14, with the fcc phase remaining stable up to pressures >400 GPa; however, the subsequent observation of a bcc crystal structure at 3300 K and 5500 K at 205 GPa29 suggests that the boundary between the bcc and fcc phases is strongly dependent on the time-dependent P-T path, and may lie at higher pressures than previously anticipated.

Static compression experiments using diamond anvil cells (DACs) and synchrotron XRD have mapped out the H2O phase diagram at lower pressures. However, experiments conducted above ~1300 K, which are performed either by direct heating of H2O using a CO2 laser (wavelength λ ~ 10.6 µm) or by indirect heating of an embedded laser absorber which interacts with near IR radiation (λ ~ 1 µm), can encounter additional complexities arising from temperature fluctuations (e.g. induced by changes in the optical properties of the sample during a phase transition13) or the presence of temperature gradients within the sample chamber30. Indirect laser heating experiments using laser couplers are particularly susceptible to temperature gradients, resulting in a discernible signal originating from the colder regions of the sample in XRD patterns. Although CO2 laser heating provides more spatially uniform heating as the laser is directly absorbed by the sample, CO2 lasers typically exhibit large power fluctuations at microsecond timescales due to the use of pulse width modulation to control the average laser power31. In addition, direct CO2 laser heating of H2O is not possible above ~60 GPa due to the reduced absorption of λ ~ 10 µm radiation above this pressure12,32. Consequently, despite extensive experimental efforts, the location of the melting curve and the SI-bcc/fcc boundary remain controversial, with recent studies reporting melting temperatures which exhibit discrepancies of ~1000 K at 50 GPa13,15,32–34.

Ionic conductivity measurements performed using resistively-heated DACs initially identified the onset of superionic conduction at T ~ 740 K at 56 GPa35; however, the maximum reported conductivity (0.12 S/cm) was ~2 orders of magnitude lower than values determined in shock compression experiments at higher pressures27,28. A subsequent XRD study reported the observation of an isostructural phase transition on compression above 15.6(2) GPa at 905 K15, which was identified as the SI-bcc to ice VII transition based on the large entropy increase at the transition and agreement with results from ab initio molecular dynamics simulations. The same work also saw evidence of SI-bcc during laser heating in the 27–44 GPa pressure range, with the appearance of SI-bcc coinciding with the first observation of fluid diffraction15. In more comprehensive studies, the SI-bcc to SI-fcc transition was observed in two independent laser heating studies12,13, but with significant discrepancies in the reported phases diagrams. In particular, Prakapenka et al. 13 observed SI-fcc during heating at pressures as low as ~30 GPa, whereas Weck et al. 12 did not observe evidence of SI-fcc below 50 GPa. Instead, Weck et al. observed fluid/SI-bcc coexistence on heating in the 27–47 GPa range at temperatures below the SI-bcc to SI-fcc transition temperatures reported by Prakapenka et al. The SI-fcc phase boundaries reported by Prakapenka et al. consistently lie at higher temperatures than those reported by Weck et al., with reported SI-bcc/SI-fcc transition temperatures differing by ~600 K at 60 GPa.

In this work, the structural properties and phase stability region of SI ice were investigated using trains of MHz pulses produced by the European X-ray Free Electron Laser (European XFEL), which were utilized to create and probe high temperature states of statically-compressed H2O. This experiment used a femtosecond serial pump-probe X-ray heating scheme, where each pulse deposits energy into the sample and simultaneously captures a snapshot of the sample as it cools from the hot state created by the previous pulse36,37. The timescales associated with this dynamic heating approach are vastly different from those experienced in conventional laser heating DAC experiments; the near-instantaneous energy deposition in the sample during single pulse irradiation produces ultrafast (sub-ns) heating followed by rapid cooling and relaxation due to heat dissipation38. Due to the long X-ray absorption length of H2O (~4000 µm at 40 GPa and 18 keV), the necessarily thin (<30 µm) samples were heated using various inert couplers which were directly heated by the XFEL beam, many of which are not suitable for IR laser heating due to their high reflectivity in the IR spectrum. This ability of indirect XFEL heating to induce bulk heating of H2O was recently demonstrated at PAL XFEL38; however, previous conclusions relied on indirect methods of melt detection due to the low repetition rate of the PAL source (30 Hz), which was insufficient to probe the hot state in situ. In our study, the structural evolution of H2O during the heating process was monitored using pulse-resolved, MHz XRD data collected using an Adaptive Gain Integrating Pixel Detector (AGIPD)39, which enabled the tracking of structural changes, including phase transformations and chemical reactivity, with sub-microsecond temporal resolution. The development of a data analysis methodology in which diffraction spots from individual crystallites were identified in XRD images allowed us to harness the extensive volume of XRD data collected from the hot sample, which significantly improved the diffracted signal of SI ice in comparison to time-integrating methods. H2O was observed to transform from ice VII to SI-bcc on heating in the 26–69 GPa pressure range, whereas SI-fcc was only observed at pressures ≥50 GPa. The absence of SI-fcc in the lower pressure runs, despite evidence of melting and SI-bcc/fluid coexistence, is in good agreement with previous CO2 laser heating studies12,15 but contrasts with results from IR laser heating experiments in which SI-fcc was observed at pressures as low as 30 GPa13. The discrepancies in the reported SI-fcc phase stability fields are discussed in terms of the timescales and P-T paths associated with different experimental approaches.

Results and discussion

X-ray heating

High-pressure H2O samples were indirectly heated using a variety of mid- and high-Z couplers (see “Methods” section, Supplementary Table 1, and Supplementary Fig. 1) which were irradiated with 300 XFEL pulses at a 2.2 MHz repetition rate. A number of different coupler geometries were employed, the majority of which were based on a doughnut-type crucible and two of which consisted of a dispersed nanopowder. One example of a doughnut-type design is shown in Fig. 1a-c, in which multiple couplers were embedded in a rhenium (Re) gasket and insulated from the diamonds by a thin self-insulating layer of H2O. Samples were subjected to multiple heating runs at different X-ray fluences by varying the beamline transport transmission to the target, starting from a low X-ray transmission (≤1%) and increasing in steps of 0.5–10%, where a run comprises a single train of X-ray pulses striking the sample. Data from doughnut-type crucibles were collected with the XFEL beam aligned to the center of the coupler hole so that heating of the coupler was primarily performed by the tails of the focused beam (outer region of the beam spot), whereas multiple spots were targeted on nanopowder couplers. Data were collected using three different X-ray focal spot sizes of <8 µm, 13 µm, and 26 µm (FWHM), which are specified in Supplementary Table 1. For some samples, the temporal dependence of the coupler temperature was constrained by streaked optical pyrometry (SOP) measurements (Fig. 1d–f).Fig. 1 Sample configuration and diagnostics used for X-ray heating of H2O at ~67 GPa.

a Photomicrograph showing the three doughnut-type couplers (Au, Cu, and Ag) in DAC 8, which were used to indirectly heat the H2O sample during XFEL irradiation. The couplers were imbedded in the Re gasket and insulated from the diamond anvils by a thin layer of H2O, as illustrated in (b) and (c). Data were collected with the XFEL beam (<8 µm FWHM) aligned to the center of the coupler hole so that the coupler was heated by the tails of the beam. d–f Data from run 533, which was collected from the Cu coupler using 70% X-ray transmission. d Pulse energy as a function of time for the 300 pulses in the train. e SOP spectrogram after fluorescence removal showing the thermal emission from the coupler during the run, where time 0 corresponds to the arrival of the first XFEL pulse. f Temporal evolution of the total SOP intensity and temperature determined from a Planck fit to the thermal emission in a 9.07 µs time window, where the horizontal error bars indicate the bin width and the vertical error bars correspond to one half of the standard deviation confidence from the Planck fit. The temporal resolution of the SOP is not sufficient to resolve temperature oscillations during the heating/cooling process (Fig. 2a); instead, it is sensitive to the hottest part of the run where thermal emittance is the brightest. Source data for panels (d)–(f) are provided as a Source Data file.

Despite its weak scattering power, the use of intense XFEL pulses enabled the collection of high-quality XRD images from H2O using a single X-ray pulse even at the lowest X-ray transmission (0.3%). Due to the short duration of individual pulses ( < 50 fs), diffraction occurs before the thermal expansion of the lattice, providing an essentially instantaneous snapshot of the sample before it is heated by the XFEL beam. XRD from the first pulse in the train therefore probes the room temperature state, and was used for pressure determination and to confirm the integrity of the sample after the previous heating run. On longer timescales, energy transferred to the lattice via electron-phonon thermalization processes results in sub-ns heating followed by more gradual cooling via thermal conduction38,40, where cooling rates depend on material properties (i.e. thermal conductivity) and geometry of the sample/coupler assembly. Serial pulse irradiation produces incremental heating with a saw-tooth like temporal temperature profile until a limiting state of thermal balance is achieved40, with each XRD image providing a snapshot of the sample as it cools from the hot state created by the previous pulse. Finite Element Analysis (FEA) simulations (Fig. 2) suggest that indirectly heated materials experience reduced thermal oscillations due to the damping effect of heat diffusion; however, oscillations are still significant in the vicinity of the coupler. Overall, this dynamic heating approach reduces the total heating time to ~130 µs per run, which is ~10,000 times shorter than typical XRD exposure times at synchrotron sources12,13,15.Fig. 2 Results of a FEA model used to simulate X-ray heating.

The model simulates X-ray irradiation of the Cu coupler in DAC 8 (run 533). a Temperature evolution at the coupler edge, in the crucible center, and 0.5 µm from the edge of the coupler. Although the coupler temperature undergoes continuous oscillations, these are strongly damped in the center of the crucible and the temperature stabilizes after ~20 pulses. However, temperature oscillations in the H2O are still significant in the vicinity of the Cu coupler (~500 K at 0.5 µm from the coupler edge). b Temperature distribution in the sample at the time of the 82nd XFEL pulse, showing conditions detected with XRD. Black lines show possible boundaries between ice VII and SI-ice, taken to be at 1500 K12 or 2000 K13. The XFEL beam is incident from below. Temperatures in (a) are reported at the blue, green, and orange points. Source data for panel (a) are provided as a Source Data file.

Observation of body-centered cubic superionic ice

Evidence of SI-bcc was observed in data collected from 9 DACs between 26 and 69 GPa (Supplementary Table 1). The evolution of XRD patterns collected during irradiation was similar for all runs (Supplementary Figs. 2 and 3), illustrated here by a dataset collected from an Ag doughnut-type coupler at 69.3 GPa (Fig. 3). In the Ag run, heating resulted in a splitting of the ice VII and Ag reflections which increased during the first ~20 pulses, with both sets of reflections shifting to lower diffraction angles due to thermal expansion. At the same time, the growth of a new reflection at ~19.2° from pulse 3 onwards is consistent with the expected d-spacing of the (110) SI-bcc reflection at this pressure12,13. No other SI-bcc reflections could be observed due to the limited 2θ detector coverage, which prevented a definitive identification of the crystal structure. However, the possibility that this peak originated from SI-fcc was ruled out because of the absence of the (200) fcc reflection, which would fall within the angular range of our detector, and the fact that indexing this new peak as the (111) fcc reflection would imply an unphysical density almost identical to cold ice VII at the start of the run (Supplementary Table 2). The observation of SI-bcc alongside hot and cold ice VII is indicative of a significant spatial temperature gradient in the probed sample volume, which is typical in such X-ray heating experiments41. Although the SI-bcc (110) reflection overlaps with the weak signal from the (101) Re reflection, it is clearly identifiable in the 2D diffraction images due to the different texture of their Debye Scherrer rings. The possibility that the spots originated from the recrystallization of Re at high temperatures was ruled out because the texture of the (100) Re ring remained unchanged during the run, and because the spots identified as SI-bcc were absent in the first pulse pattern of the subsequent run.Fig. 3 X-ray diffraction collected during XFEL irradiation of H2O at 69.3 GPa.

Data were collected using an Ag doughnut coupler (DAC 8) and an X-ray beam diameter of <8 µm FWHM. a, b Integrated XRD patterns and (c–h) unwrapped (2θ-φ) XRD images collected during irradiation with 300 XFEL pulses at 29% transmission (run 505). In (a), darker shades correspond to higher intensity. In (c–h), lighter shades correspond to higher intensity. Heating results in a splitting of the ice (110), Ag (111), and Ag (200) reflections, which is indicative of two distinct temperatures in the probed volume, predominantly originating from the lack of thermal insulation between the H2O and the diamond anvils. Panel (a) illustrates how heating occurs during the first ~20 pulses, after which the diffracted signal from H2O remains essentially constant. In (d), the diameters of the orange circles are unrelated to individual spot sizes. Source data are provided as a Source Data file.

Following the disappearance of the hot ice VII reflection around pulse 20, no significant changes were observed for the remainder of the run (Fig. 3a and Supplementary Fig. 4), suggesting that the system reached a balance between X-ray heating and heat loss via thermal conduction40. Based on this observation, images from pulses 51-300 were summed and integrated to produce a 1D diffraction profile from the hot sample (Fig. 4a) in which the (110) SI-bcc reflection is clearly visible alongside the broad, hot (110) ice VII reflection. Due to the spotty nature of the H2O signal, it was possible to analyze the unwrapped (2θ-φ) 2D diffraction images using a spot-finding algorithm (see Methods) to determine the 2θ position of H2O diffraction spots in each image, which was used to produce a histogram showing the total number of ice spots in images collected from pulses 51-300 (Fig. 4a and Supplementary Movie 1). The resultant histogram is in excellent agreement with the summed pattern after subtracting the Re signal, validating this approach. This analysis procedure was repeated for all runs from this sample in which SI ice was observed, and the datasets were combined to produce a single integrated profile and histogram (Fig. 4b). Although these runs were collected with different X-ray transmissions (25.5–37%) which are expected to produce different temperatures, the SI-bcc reflection is clearly distinct from that of ice VII, confirming that the observed peak separation is not an artifact relating to the large temperature gradient within the heated area. The narrow width of the SI reflection is most likely due to the narrow temperature stability field of this phase at 69 GPa, which is above the upper pressure limit of the SI-bcc stability field reported in previous work12,13.Fig. 4 Integrated patterns from summed images and histograms showing the results of the spot finding approach.

Data are shown for (a) run 505, and for (b) all runs collected from the Ag coupler in DAC 8 where SI was observed. The construction of the histogram for run 505 is illustrated in Supplementary Movie 1. When the Re contribution (estimated from the first pulse pattern) is subtracted from integrated patterns, they are in good agreement with the histogram. Source data are provided as a Source Data file.

Observation of face-centered cubic superionic ice

The results from the Ag coupler in DAC 8 (Figs. 3 and 4) are compared to runs collected from Au and Cu couplers in the same DAC (Supplementary Fig. 5), which were analyzed using the summed and spot-finding methods described above (Fig. 5). The SI-bcc reflection is clearly visible in both the summed XRD pattern and histogram produced from the Au runs, and the lattice parameters determined using each approach are in excellent agreement with each other (a = 2.856 Å and 2.854 Å from the summed pattern and histogram, respectively), and with those determined from the Ag runs (a = 2.844 Å and 2.843 Å) where the pressure was 5.2 GPa higher. Although the summed pattern from the Cu runs is difficult to interpret due to the broad background originating from the detector behavior at high signal levels, the SI-bcc reflection is clearly visible in the histogram alongside a second, well separated, lower-angle reflection. Indexing this peak as the (111) SI-fcc reflection determines a density of 2.618 g/cm3, which is almost identical to the SI-bcc density of 2.619 g/cm3, in agreement with previous work which found that both SI phases can be described by the same equation of state (EoS)13. In addition to the observation of the (111) SI-fcc reflection in the 90 and 100% transmissions runs (Fig. 6), the (200) SI-fcc reflection is also present in the 80% and 100% runs (Supplementary Fig. 6), further confirming the nature of this phase. On subsequent re-examination of the Ag data, a small number of isolated SI-fcc diffraction spots were located in several XRD images from 2 different runs (Supplementary Fig. 7), confirming that the fcc reflections in the Cu runs originated from H2O rather than a chemical reaction between the sample and the coupler. However, no evidence of SI-fcc was located in any images from the Au runs from this DAC.Fig. 5 Histograms showing the number of H2O diffraction spots in runs collected from different samples.

For each sample, histograms were compiled from XRD data collected from pulses 51-300 in runs where SI ice was observed. The width of each bin is 0.01 degrees. Histograms are compared with integrated XRD patterns produced from the summed Adaptive Gain Integrated Pixel Detector images of the hot (pulses 51–300) and cold (pulse 1) sample. The y-axis refers to the histogram plot, and the integrated patterns from the hot and cold sample were normalized to a single arbitrary scaling factor for the figure. The coupler material is specified in each case. The high-angle signal in the DAC 3 and 5 runs are from the hot coupler. With the exception of DAC 1 and DAC 2, in which the coupler was in the form of a dispersed nanopowder, all couplers were of a doughnut-type design. Source data are provided as a Source Data file.

Fig. 6 XRD data collected during X-ray heating of H2O at 67.3 GPa.

Data were collected during irradiation of the embedded Cu doughnut coupler in DAC 8 using an X-ray beam size of <8 µm FWHM. Integrated XRD data are shown from runs (a) 532 and (b) 536, and single pulse, unwrapped (2θ-φ) XRD images from runs (c) 532 and (d) 536 are shown for a reduced 2θ range. In (a–b), darker shades correspond to higher intensity, whereas in (c–d), lighter shades correspond to higher intensity. The data shown in (a) and (c) were collected using 60% transmission, and the (110) SI-bcc reflection is observed from pulse ~5 onwards (see insert), The SI-bcc diffracted signal is predominantly from several intense spots, such as the one highlighted in (c). The data in (b) and (d) were collected at 100% transmission, and both SI phases (bcc and fcc) are observed. Source data are provided as a Source Data file.

The difficulty in forming SI-fcc in these runs is surprising when SI-bcc formed so readily. Although SI-bcc is stable at lower temperatures than SI-fcc (Fig. 7), these runs were performed at pressures of 64-69 GPa, which are close to (or above) the previously-reported upper limit of the SI-bcc pressure stability field. In particular, Weck et al. 12 did not observed SI-bcc above 57 GPa, and although Prakapenka et al. 13 reported a single SI-bcc data point at 67 GPa and ~2000 K, there is no evidence of SI-bcc in their 2500 K diffraction pattern at 67 GPa, despite clear evidence of SI-fcc and ice VII. Although the limited azimuthal coverage of the AGIPD may account for the absence of SI-fcc in certain runs, it is unlikely to be the only contributing factor given the consistently strong signal observed from SI-bcc. Instead, the difficulty in forming SI-fcc points to an intrinsic feature of the pump-probe X-ray heating scheme and the specific coupler materials and geometry used in these runs, combined with a large energy barrier to form SI-fcc. As XRD provides a snapshot of the sample during cooling, the sample temperature is heavily dependent on heat distribution to the surrounding media, and the high thermal conductivities of Au, Ag, and Cu (318, 429, and 401 W/mK at ambient conditions, respectively42) limit heating of H2O by fast heat transfer away from the sample. Although the diffracted signal from hot Ag (Supplementary Fig. 8) suggests that a portion of the coupler remained close to its melting point after 443 ns, the coupler reflections in the Cu and Au runs remain essentially unchanged from the cold pattern (Supplementary Fig. 9) - despite SOP determining an average coupler temperature (Supplementary Table 3) which is well above the previously-reported lower temperature limit of SI-fcc12,13 – suggesting that the coupler cooled significantly between subsequent XFEL pulses. These observations are supported by FEA simulations (Fig. 2a) which suggest that the Cu becomes colder than the sample during the cooling process. Enhanced cooling in the real sample compared to the FEA is attributed to the significant reduction in the H2O insulation layer thickness at high pressure, potentially leading to direct contact between the coupler and the diamond anvil.Fig. 7 Average coupler temperatures in experiments where SI ice was observed.

Data points correspond to the average SOP temperature during the run, and the error bars show the standard deviation. Source data are provided as a Source Data file. Thermal emission was not detectable for all samples, and data points are shown for runs in which thermal emission was observed. Due to the low emissivity of H2O, SOP measures thermal emission from the coupler surface, which does not necessarily correspond to the temperature in the H2O sample. Our data points are compared to the phase diagram of H2O reported in previous studies. The text labels indicate the phase stability regions reported in ref. 13. For simplicity, we do not make any distinction between ice VII and VII’ in the main text. SI-bcc and SI-fcc data points from previous static12,15,50 and shock compression14,29 studies are included for comparison, as well as ice VII data points from shock compression work14,29 and the H2O melting curves from ref. 15,32. Ice VII data points from static compression studies are not included to avoid overcomplicating the figure. The data points from ref. 15. indicate the SI-bcc/ice VII phase line determined on isothermal compression.

Optimal parameters for indirect X-ray heating

The summed XRD profiles and diffraction spot histograms are shown in Fig. 5 for DACs 1 - 8. The ice VII and SI-bcc peaks are very clearly separated at all pressures, particularly in the histogram plots, due to the lower photon energy (18 keV) compared to previous DAC experiments (>30 keV). Bulk heating is particularly pronounced in runs collected from the Rh doughnut-type coupler in DAC 9 (Supplementary Fig. 10), evidenced by the strong SI-bcc signal and almost complete disappearance of ice VII signal, most likely due to a combination of optimal geometrical parameters and coupler properties. In particular, the lower thermal conductivity of Rh (150 W/mK at ambient conditions) compared to Au, Ag, and Cu is expected to reduce heat dissipation in the unheated coupler, resulting in a more uniform temperature distribution in the sample. However, although Rh was the most promising coupler in terms of heating efficiency, the appearance of a number of unidentified diffraction spots after the first observation of SI-bcc in some runs suggests that Rh reacts with SI H2O. Data from DAC 9 were therefore not analyzed using the spot finding algorithm, but are presented to illustrate that uniform heating of low-Z samples using a high-Z coupler is possible by careful choice of X-ray parameters, sample geometry, and material parameters. Although chemical reactivity was a general problem for Rh couplers, one run from DAC 7 was included in the analysis because reaction products were not observed until the end of the run, after all traces of SI-bcc had disappeared.

Absence of face-centered cubic superionic ice at low pressures

Examination of individual XRD images from the lower pressure (<60 GPa) runs revealed a single (111) SI-fcc crystallite (Supplementary Fig. 11) at 49.9 GPa (DAC 6). Although this does not provide conclusive evidence for the existence of SI-fcc, it is consistent with the expected stability field of this phase reported by both ref. 13,12 and so was included in further analysis. No evidence of SI-fcc was observed in data collected from DACs 1-5 (26–38 GPa); however, the dramatic reduction in the H2O diffracted signal in the highest fluence runs collected from DAC 4 (~38 GPa) indicates that the bulk of the sample had melted. This is supported by SOP measurements, which determined the average coupler temperature to be above the ice melting temperature at this pressure. Melting is particularly striking in the 90% Ag run due to the strong signal from ice VII at the start of the run (Fig. 8 and Supplementary Movie 2), which almost completely disappears by pulse 30. Individual diffraction spots from both bcc phases reappear throughout the run, suggestive of solid-fluid coexistence. Although no evidence of diffuse scattering from the fluid was observed due to the weak scattering power of H2O and the limited azimuthal coverage of the detector, the possibility that the disappearance of the H2O solid signal originated from chemical reactivity or sample loss was ruled out by the strong ice VII signal and absence of chemical reaction products in the first pulse pattern of the subsequent run (Supplementary Fig. 12).Fig. 8 Unwrapped (2θ-φ) XRD images showing evidence of melting of H2O.

Data were collected using 90% X-ray transmission from the embedded Ag coupler in DAC 4 (36.7 GPa) using an X-ray beam size of <8 µm FWHM (run 917:17). Lighter shades correspond to higher intensity. Source data are provided as a Source Data file. The significant reduction in the intensity of the ice VII/SI-bcc diffracted signal during the run is attributed to melting of H2O. No evidence of the SI fcc (111) reflection was observed, which would be expected to be present at approximately 18.4°. Unwrapped images from the full run are shown in Supplementary Movie 2.

The absence of SI-fcc in the 38–40 GPa runs is surprising due to the strong evidence for melting and fluid/SI-bcc coexistence, and the fact that previous work13 observed SI-fcc at pressures as low as 29 GPa. Instead, our results are in agreement with previous CO2 laser heating studies12,15 which did not observe SI-fcc during heating at pressures up to 45 GPa, but instead observed SI-bcc to coexist with the fluid. In order to resolve the discrepancy between results obtained using different heating techniques (X-ray and CO2 vs. IR laser heating), it is necessary to consider the differences and similarities between each experimental approach. The results do not appear to be correlated with the heating method (i.e. direct or indirect heating), as both Prakapenka et al. 13 and this work performed indirect heating using embedded absorbers and obtained different results. Heating durations of ≥5 s in the CO2 and IR heating experiments are 4 orders of magnitude longer than XFEL heating (~130 µs), which seemingly suggests that timescale is not a factor. However, as the output power of CO2 lasers is typically controlled by pulse width modulation which produces fluctuations in laser intensity31, the time that the sample lies within SI-fcc stability field could be reduced by the resultant thermal fluctuations which could potentially cause the sample to oscillate across the SI-fcc phase boundary. Although SI-fcc may be thermodynamically stable below the melting line, previous theoretical studies calculated a lower solid/fluid interfacial free energy for the SI-bcc lattice in comparison to SI-fcc43, suggesting that the fcc lattice may have a longer nucleation time than SI-bcc when crystallizing from the fluid. We therefore propose that the absence of SI-fcc in the 38–40 GPa XFEL heating runs is due to short heating timescales and the P-T path associated with XFEL pump-probe experiments, combined with kinetic hindrance associated with the formation of SI-fcc from the melt. This hypothesis is consistent with the results of IR laser heating experiments using the same approach as Prakapenka et al. (Supplementary Fig. 13), which were performed at the GSECARs beamline at the Advanced Photon Source using Au flakes as a laser absorber. After initially heating H2O above its melting temperature, molten ice was first observed to coexist with SI-bcc when the laser power was reduced, whereas the observation of SI-fcc on further decrease of the laser power coincided with the disappearance of the diffuse scattering signal. The observation of liquid/SI-bcc coexistence is attributed to localized heating and melting of the H2O in the vicinity of the Au flakes, which drives the movement of the Au particles and results in rapid melting and recrystallization of H2O. Consequently, the sample does not remain in the SI-fcc stability field long enough for SI-fcc to form. Similarly, thermal fluctuations in the vicinity of the coupler in the XFEL heating experiments (Fig. 2a) due to the pulsed nature of the XFEL beam can account for fast recrystallization of SI-bcc from the fluid in the 90% transmission run collected from DAC 4 (Fig. 8), suggesting an approximate lower bound of 443 ns for the nucleation time. The observation of SI-fcc at higher pressures in the XFEL experiment is attributed to the higher melting temperature of H2O at this pressure and the fact that the average coupler temperature lies in the stability field of SI-fcc, which avoids transforming a portion of the sample into the fluid phase. This is consistent with the observation of large SI-fcc crystallites at the same azimuthal position for multiple frames in the 100% transmission run collected from DAC 8 (Fig. 6), rather than small, recrystallizing spots which are typically formed during fast quenching. The possibility of timescale- and path-dependent transitions in high-temperature H2O should therefore be taken into account when comparing results obtained using different experimental techniques; in particular, considering the conflicting results of recent shock compression experiments (Millot et al. observed SI-fcc near 160 GPa and 3200 K14, whereas Gleason et al.29 saw SI-bcc as pressures as high as 207(10) GPa at 5500(500) K).

Density of superionic ice

The histograms in Fig. 5 were used to determine the densities of both SI phases, which are compared with the results of previous studies in Fig. 9. Our SI-bcc points are in good agreement with the results previous studies12,13; our highest-pressure data points lie closer to the Prakapenka et al. 13 EoS which includes a thermal pressure correction, despite the fact that such a correction was not included in our analysis, whereas our lower-pressure points lie closer to their uncorrected curve. The SI-fcc densities determined in this study are in excellent agreement with those of SI-bcc at the same pressure, in agreement with Prakapenka et al. 13, but lie at lower densities than those reported by Weck et al. 12 at the SI-bcc/SI-fcc transition temperature. Weck et al.12 reported the density of SI-fcc to be strongly temperature dependent, and our SI-fcc data are in better agreement with the density they report at 57 GPa and 1927 K, which is ~500 K above their reported SI-bcc/SI-fcc transition temperature. Although they report a corresponding SI-bcc density for the 57 GPa run which is significantly larger (~2.67 g/cm3)33 than our measured SI-bcc density at the same pressure (~2.55 g/cm3), our data are in good agreement with the SI-bcc density determined from the 1927 K XRD pattern in Fig. S4 of their supplementary material (~2.55 g/cm3), suggesting that the reported SI-bcc density actually corresponds to that of expanded ice VII. The possibility that SI-fcc was not observed in our XFEL experiments until ~500 K above the transition temperature is intriguing considering the difficulty in forming SI-fcc, and is consistent with a large energy barrier associated its formation. However, due to the lack of accurate temperature measurements in our XFEL experiments, it is not possible for us to resolve discrepancies between previous work. Overall, the large density difference between ice VII and SI-bcc reported in this work contrasts with the 108 GPa SI-bcc data point from the shock compression study of Gleason et al. 29, which suggests that the two phases have similar densities. However, as their SI-bcc data point lies outside the SI stability field reported by ref. 13, it is possible that that this sample was actually in the ice VII stability field.Fig. 9 Density as a function of pressure for different phases of H2O.

The data points from this work show the SI-bcc and SI-fcc densities calculated from the histograms in Fig. 5. Source data are provided as a Source Data file. The data point from DAC 9 (34.2 GPa) is included to illustrate the agreement with other data, despite weak evidence of chemical reactivity in this sample. Our results are compared with SI ice data from previous DAC experiments12,13,15, which are identified as using laser heating (LH-DAC) or resistive heating (RH-DAC) techniques. Ice VII data from ambient temperature12,13,51,52 and high temperature35 DAC experiments, shock compression experiments14,29, as well as data from the high temperature fluid53, are shown for comparison. In addition to the data points taken from the main paper of Weck et al.12, which indicate the SI-bcc and SI-fcc densities at the transition temperature, their SI-fcc data point at 57 GPa and 1927 K is also included to indicate the volume of thermally-expanded SI-fcc at higher temperature.

In conclusion, an indirect XFEL heating technique was used to investigate the structural properties of SI ice, demonstrating that stepwise XFEL heating is a viable method to study low-Z materials such as H2O at high pressures and temperatures. Evidence of SI-bcc was identified in a large number of runs performed at 27-69 GPa, illustrating that this approach is able to detect isostructural transitions, and SI-fcc was observed in runs ≥50 GPa. The large volume of XRD data collected using a pulse-resolved, MHz detector, combined with the highly textured nature of the SI-bcc diffraction signal, motivated the development of a new data analysis approach in which the identification of diffraction from individual crystallites in XRD images was used to enhance the signal from high-temperature SI phases and determine their density as a function of pressure. No evidence of SI-fcc was observed in lower pressure runs, despite clear evidence of melting and the simultaneous observation of SI-bcc and fluid at ~40 GPa. Based on previous theoretical studies which calculated a lower solid/fluid interfacial free energy for the SI-bcc lattice in comparison to SI-fcc43, and considering prior observations of SI-fcc forming at these P-T conditions upon heating13, we attribute the absence of SI-fcc in these experiments to the short heating timescale combined with the P-T path of the pump-probe approach in which SI ice is formed on cooling from the melt. The results may have important implications for the stability of SI phases in solar and extrasolar ice-rich planets during dynamical freezing such as in internal convection processes, where the preferential crystallization of SI-bcc from the fluid may result in different physical properties of the solid ice (e.g. electrical and thermal conductivity) across the same internal layer that in turn, may affect interior dynamics and magnetic fields.

Methods

Sample preparation

H2O samples were loaded into a total of 9 symmetric DACs (DACs 1–9) with standard design diamonds mounted on the upstream side and Boehler Almax diamonds facing downstream. Microscope images of the samples before the XFEL experiment are shown in Fig. 1a and Supplementary Fig. 1, which were collected after the sample was compressed to the final pressure for data collection. A total of 5 different coupler materials were employed (Au, Ag, Cu, Rh, and Fe3O4) in different coupler geometries (nanopowder, doughnut, and embedded doughnut) which were designed with the aim of minimizing the temperature gradient in the indirectly heated H2O samples (for preparation details see Supplementary Note 1). With the exception of DAC 7, all samples were loaded without a thermal insulation layer, relying instead on a self-insulating H2O layer close to the diamonds to act as a thermal barrier. This avoided the need for commonly-used insulation layers such as Al2O3, LiF etc., which can complicate the interpretation of diffraction patterns and/or result in unwanted chemical reaction products. After loading, samples were pre-compressed to the desired pressure for the XFEL experiment

Experimental details: sample screening at PETRA III

After compression to the desired pressure for the European XFEL experiment, samples were screened at the P02.2 beamline at the PETRA III synchrotron source to confirm the integrity of the sample loading using XRD. The screening was performed using 25.6 keV X-rays focused to 3(v) × 8(h) µm2 (FWHM) using a series of compound refractive lenses (CRLs), and XRD data were collected using an XRD1621 area detector (PerkinElmer). After the XFEL experiment, a 2D diffraction map collected from DAC 8 found no evidence of chemical contamination. This was performed at the P02.2 beamline using 42.7 keV X-rays focused to 2(v) x 2(h) µm2 (FWHM) using Kirkpatrick–Baez mirrors, and data were collected using an XRD1621 area detector.

Experimental details: European XFEL parameters

X-ray heating experiments were performed at the High Energy Density (HED) Instrument44 at the European XFEL using the dedicated DAC platform in Interaction Chamber 236. The experiment was performed using an X-ray photon energy of 18.047 keV and a 2.2 MHz intra-train repetition rate, which corresponds to a 443 ns spacing between consecutive pulses, where the photon energy was calibrated using the HIgh REsolution hard X-ray single-shot (HIREX) spectrometer45. X-rays were focused using a series of CRLs using three different configurations to achieve focal spot sizes of <8 µm, 13 µm, and 26 µm (FWHM), where the beam size was estimated from edge scans using a polished W rod. Accurate determination of the beam diameter for the smallest focal spot size was not possible based on the edge scan data due to erosion of the W rod.

Experimental details: X-ray diffraction at European XFEL

Pulse-resolved XRD data were collected using a 500k AGIPD detector positioned outside of the vacuum chamber at a sample-to-detector distance (SDD) of ~430 mm, which provided an angular coverage of ~7.5 to 27 degrees (q = 4.27-1.20 Å−1). The detector geometry (SDD, tilt, and rotation) was calibrated using a Cr2O3 diffraction standard (NIST SRM 674b) using the Dioptas software46. Radial integration of diffraction images was performed using Dioptas and azimuthally unwrapped (2θ-φ) diffraction images were produced using FIT2D47. Prior to integration, the intensity of individual diffraction images was scaled using the pulse-resolved intensity and position monitor (IPM)36 positioned upstream of the sample to account for fluctuations in the pulse energy across the pulse train, and all images within the run were normalized to the mean IPM value within the run to allow for comparison of runs collected at different levels of X-ray attenuation. The X-ray fluence incident on the sample during the run was controlled using a series of solid attenuators, which enabled the X-ray transmission to be varied from 0.3 to 100%. The pulse energy on target was determined using the IPM, which was calibrated using an X-ray gas monitor (XGM)36 at the start of the experiment.

Alignment of the sample to the X-ray beam was performed based on visual observation using an optical imaging system. This alignment method resulted in an offset in the sample position along the X-ray beam direction with respect to the diffraction calibrant due to refraction from the downstream diamond anvil, which was corrected by changing the SDD used in the detector calibration based on the measured diamond thickness and refractive index. The X-ray position on the optical camera was determined by exposing an area on the gasket for several seconds at low fluence, which produced a small dark region resulting from X-ray damage. This approach was found to be sufficient to align the XFEL beam with small holes (6–22 µm diameter) in doughnut couplers, where the XFEL beam was positioned at the center of the hole to maximize the diffracted signal from the H2O sample. However, the positional jitter in the XFEL beam position (approximately equal to the X-ray focal spot size) meant that the beam position varied from run to run, which could potentially result in different degrees of X-ray heating in subsequent runs due to the large difference in absorption lengths of H2O and the mid/high-Z couplers.

Experimental details: streaked optical pyrometry

The temporal evolution of the temperature profile during the heating run was determined from streaked optical pyrometry (SOP) measurements performed using a streak camera (C13410-01A, Hamamatsu) coupled to an optical spectrometer (IsoPlane 160, Princeton Instruments), which collected thermal emission from a 50 µm diameter region on the upstream side of the sample. Full details of the SOP system are given in ref. 49. In all cases, data were collected using a 200 µs streak window to cover the entire duration of the 132.9 µs long X-ray pulse train. SOP data collected using a low X-ray fluence were dominated by a fluorescence signal originating from the diamond anvils and/or H2O media, which was most prominent in the short wavelength (~500 nm) range. The presence of thermal emission at higher X-ray fluence was identified by a discontinuous increase in the total SOP intensity as a function of incident pulse energy when all runs from the same coupler were compared. For runs in which thermal emission was observed, the fluorescence signal was removed by subtracting the fluorescence-only spectrogram collected at the highest X-ray transmission. After the fluorescence correction was applied, the temperature was determined by fitting a Planck function to the 600–775 nm wavelength region of the thermal emission spectrum from each 9.1 µs temporal window assuming gray-body emission. The reliability of each temperature measurement was evaluated based on the criteria outlined in ref. 49. Examples of SOP data treatment, including the fluorescence correction, are shown in Supplementary Fig. 14. A summary of SOP temperatures from all runs are given in Supplementary Table 3, and correspond to the average over the entire run. Due to the low emissivity of H2O in the temperature region investigated in this work, SOP provided a measurement of the coupler surface temperature at the hot spot, rather than a direct measurement of the sample temperature, which is expected to be lower than that of the coupler (Fig. 2a). As the total radiated energy is proportional to T4, the strongest contribution to the thermal radiation spectrum is from the hottest region on the coupler surface. Thermal emission was not detectable for all samples, which we attribute to variations in the volume of heated coupler material due to a combination of factors such as X-ray focal spot size and coupler geometry.

Experimental details: IR laser heating at GSECARS

IR laser heating experiments were performed at the GSECARS undulator beamline (sector 13, APS, ANL). The experiment was performed using 37.07 keV X-rays focused to 3 ×4 µm2 (FWHM), and XRD data were collected using a MAR-165 CCD (charge-coupled device). Doubled-sided, coaxial heating of samples was performed using a near-IR (1064 nm) laser with a 10 µm diameter flat-top focal spot49. The temperature was determined using spectroradiography measurements performed using a Princeton grating spectrometer (300 mm focal length) combined with PIXIS and PI-MAX3 CCD array detectors, and thermal emission was collected from both sides of the heated sample. For this experiment, the H2O sample was thermally insulated from the diamond by a layer of SiO2, and small flakes of Au were used as the laser absorber.

Data analysis: X-ray diffraction

Out of a total of 400 runs collected from H2O samples, evidence of SI-bcc was identified in 68 runs (Table SI) which were collected using a range of different coupler materials, geometries, and X-ray transmissions. In all cases, the (110) SI-bcc reflection was clearly distinct from the broad, hot ice VII reflection, consistent with the expected behavior for an isostructural transition. Two different approaches were used to evaluate the XRD patterns collected from the hot sample. In the first approach, images collected from pulse 51–300 from each run contained SI were summed to produce a single integrated XRD profile for each sample. For DACs in which multiple couplers were used (DACs 4 and 8), individual diffraction profiles were produced for each coupler type. For each sample, the cold pattern was produced by summing the first pulse XRD images from the runs used for data analysis, which determined the average pressure across all runs. In order to avoid systematic uncertainties introduced by the choice of EoS for different coupler types, the sample pressure was estimated based on the position of the (110) ice VII reflection in the cold pattern using the EoS of Prakapenka et al. 13. The result of the summed approach is illustrated for DAC 1 in Supplementary Fig. 15.

For construction of the histograms shown in Fig. 5, a spot finding algorithm was used to identify the 2θ-φ positions of individual diffraction spots in the unwrapped (2θ-φ) diffraction images using a custom python code. First, a threshold of ībk+nσibk was applied, where ībk and σibk are the mean and standard deviation, respectively, of the background intensity of a region on the detector in which no diffraction lines were observed. A value of n = 6 was used in all cases except for the Cu coupler in DAC 8, where n = 5 was used due to an increase in detector noise at the highest X-ray fluence due to issues related to gain switching. Individual reflections were identified as connected regions in the threshold image, and those with a pixel size of <2 were assumed to due to detector noise and discarded. The number of local maxima was computed for each of the identified regions. If multiple local maxima were present in a single identified region and separated by at least 3 pixels, the region was split into multiple reflections using a watershed algorithm. Finally, the angular coordinates (2θ-φ) of each reflection was determined from its center of mass, and the resultant list of 2θ positions was used to produce the histogram. The value of n was chosen by visual inspection of the number of spots located in the image (Supplementary Fig. 16). However, although the choice of n determined the number of spots which were identified, histograms produced using different values of n were found to be in good agreement with one another (Supplementary Fig. 17).

Finite Element Analysis

Finite Element Analysis (FEA) of sample heating, accounting for sample geometry, the materials used and their thermochemical parameters at relevant high pressure conditions, were performed using previously described techniques38,40,41,48. FEA was used primarily in planning experiments, to establish beam properties and sample designs to avoid coupler melting and damage, achieve sufficiently high temperature in the water sample, optimize heat transfer to the water, and minimize temperature gradients in the crucible. While major predictions of the models were confirmed by pyrometry measurements during experiments, including absolute temperatures reached and the achievement of a quasi-steady state at longer timescales, models predict a gradual rise in temperature at early times which contrasted with an observed decrease with time from elevated initial temperature. This effect could be related to beam misalignment with holes, coupler melting and movement as well as changes in sample optical properties at high temperature41. Due to the use of longer sweep windows with lower time resolution in SOP, the individual heating and cooling events seen in the models are not resolved48. For the model in Fig. 2, representative high-pressure parameters of thermal conductivity k, heat capacity at constant pressure CP, density ρ, and absorptivity α were used for the H2O sample (k = 20 W m−1 K−1, CP = 3050 J kg−1 K−1, ρ = 2694 kg m−3, α = 201 m−1) and the Cu coupler (k = 450 W m−1 K−1, CP = 472 J kg−1 K−1, ρ = 9500 kg m−3, α = 4.41 × 104 m−1), along with standard parameters for the gasket and diamonds40.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Movie 1

Supplementary Movie 2

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52505-0.

Acknowledgements

We acknowledge European XFEL in Schenefeld, Germany, for provision of X-ray free-electron laser beam time at Scientific Instrument HED (High Energy Density Science) and would like to thank the staff for their assistance. These data were collected as part of a DAC community proposal (#2590, by McMahon and Husband: 10.22003/XFEL.EU-DATA-002590-00). The authors are indebted to the HIBEF user consortium for the provision of instrumentation and staff that enabled this experiment. We acknowledge DESY (Hamburg, Germany), a member of the Helmholtz Association HGF, for the provision of experimental facilities. Parts of this research were carried out at PETRA III (beamline P02.2). V.B.P. and S.C. acknowledge the support of GeoSoilEnviroCARS and National Science Foundation – Earth Sciences (EAR - 1634415). We acknowledge E. Shevchenko (CNM, ANL), who synthesized the nano-Fe3O4 samples. J.D.M. acknowledges support from AWE CASE studentship P030463429. Support is acknowledged from the U.K. Engineering and Physical Sciences Research Council (EPSRC) Grant Nos. EP/R02927X/1 (E.J.P. and M.I.M.) and EP/P024513/1 (R.S.M). M.F. Acknowledges DOE FES funding FWP100182. Y.L. is grateful for support from the Leader Researcher programme (NRF-2018R1A3B1052042) of the Korean Ministry of Science and ICT (MSIT). This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344 (H.C.). This research was supported through the European Union’s Horizon 2020 research and innovation Programme (ERC grant 864877, H.M., and 101002868, R.S.M.) as well as UKRI STFC grant ST/V000527/1 (H.M.). A.F.G. and E.E. are grateful for the support of Carnegie Science and NSF EAR-2049127. A.F.G. is grateful for the support of NSF CHE- 2302437. S.N. and L.A. acknowledge financial support from Sorbonne University under grant Emergence HP-XFEL. We acknowledge support from the Deutsche Forschungsgemeinschaft (DFG) Research Unit FOR 2440 grants SA2585/5-1 (R.J.H, A.M., C.S.V., and H.P.L) and AP262/2-2 (K.A.).

Author contributions

R.J.H., M.I.M., H.P.L., J.D.M., R.S.M., A.F.G., V.B.P., Z.K., C.S., C.S.V., M.F., K.A., O.B.B., R.B., A.L.C., H.C., Y.L., H.M., S.N., E.P., C.Pe., C.Pr, and U.Z. were involved in the conception of the experiment and writing of the proposal. R.J.H., M.I.M., H.P.L., J.D.M., R.S.M., A.F.G., V.B.P., E.E., S.C., Z.K., C.S., C.S.V., M.F., L.A., K.A., O.B.B., R.B., J.B., V.C., J.C., A.L.C., H.C., H.H., E.K., Y.L., X.L., H.M., A.M., M.N., S.N., E.P., C.Pe., C.Pr, and U.Z. participated in discussions of the experimental approach and data analysis. R.J.H., J.D.M., V.P., S.C., E.E., and A.F.G. prepared the samples. R.J.H. performed the data analysis and wrote the manuscript with input from all authors. R.S.M. performed the FEA calculations. R.J.H, H.P.L., J.D.M, C.S.V., L.A., O.B.B. A.M., S.N., C.Pr., and M.I.M. performed the experiment at European XFEL, and R.S.M, A.F.G., V.B.P., E.E., S.C., R.B., J.B., J.C., A.L.C., H.C., H.H., E.K., Y.L., X.L., H.M., E.P., and C.Pe. provided remote data analysis support. Z.K., C.S., K.A., V.C., C.B., A.D., and M.N. provided support at the HED instrument (European XFEL). H.G., T.L., S.S., and J.S.-D. provided support for the AGIPD detector. R.J.H., J.D.M., and H.H. performed screening experiments P02.2 beamline (PETRA III), and H.P.L. provided support at P02.2.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

The DOI for the original European XFEL data is: 10.22003/XFEL.EU-DATA-002590-00, and will be publicly available after the embargo period of 3 years. Source data are provided with this paper.

Code availability

The computer code used to generate the results reported in this study is available from the corresponding author upon request.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Mishima O Calvert LD Whalley E ‘Melting ice’ I at 77 K and 10 kbar: a new method of making amorphous solids Nature 1984 310 393 395 10.1038/310393a0
Mishima, O., Calvert, L. D. & Whalley, E. ‘Melting ice’ I at 77 K and 10 kbar: a new method of making amorphous solids. Nature 310, 393–395 (1984).
2. Mishima O Calvert LD Whalley E An apparently first-order transition between two amorphous phases of ice induced by pressure Nature 1985 314 76 78 10.1038/314076a0
Mishima, O., Calvert, L. D. & Whalley, E. An apparently first-order transition between two amorphous phases of ice induced by pressure. Nature 314, 76–78 (1985).
3. Salzmann CG Radaelli PG Slater B Finney JL The polymorphism of ice: five unresolved questions Phys. Chem. Chem. Phys. 2011 13 18468 18480 10.1039/c1cp21712g 21946782
Salzmann, C. G., Radaelli, P. G., Slater, B. & Finney, J. L. The polymorphism of ice: five unresolved questions. Phys. Chem. Chem. Phys. 13, 18468–18480 (2011).21946782
4. Polian A Grimsditch M New high-pressure phase of H2O: Ice X Phys. Rev. Lett. 1984 52 1312 1314 10.1103/PhysRevLett.52.1312
Polian, A. & Grimsditch, M. New high-pressure phase of H2O: Ice X. Phys. Rev. Lett. 52, 1312–1314 (1984).
5. Li J Ross DK Evidence for two kinds of hydrogen bond in ice Nature 1993 365 327 329 10.1038/365327a0
Li, J. & Ross, D. K. Evidence for two kinds of hydrogen bond in ice. Nature 365, 327–329 (1993).
6. Goncharov AF Dynamic ionization of water under extreme conditions Phys. Rev. Lett. 2005 94 125508 10.1103/PhysRevLett.94.125508 15903935
Goncharov, A. F. et al. Dynamic ionization of water under extreme conditions. Phys. Rev. Lett. 94, 125508 (2005).15903935
7. Kuhs WF Finney JL Vettier C Bliss DV Structure and hydrogen ordering in ices VI, VII, and VIII by neutron powder diffraction J. Chem. Phys. 1984 81 3612 3623 10.1063/1.448109
Kuhs, W. F., Finney, J. L., Vettier, C. & Bliss, D. V. Structure and hydrogen ordering in ices VI, VII, and VIII by neutron powder diffraction. J. Chem. Phys. 81, 3612–3623 (1984).
8. Jackson SM Nield VM Whitworth RW Oguro M Wilson CC Single-crystal neutron diffraction studies of the structure of ice XI J. Phys. Chem. B 1997 101 6142 6145 10.1021/jp9632551
Jackson, S. M., Nield, V. M., Whitworth, R. W., Oguro, M. & Wilson, C. C. Single-crystal neutron diffraction studies of the structure of ice XI. J. Phys. Chem. B 101, 6142–6145 (1997).
9. Benoit M Romero AH Marx D Reassigning hydrogen-bond centering in dense ice Phys. Rev. Lett. 2002 89 145501 10.1103/PhysRevLett.89.145501 12366053
Benoit, M., Romero, A. H. & Marx, D. Reassigning hydrogen-bond centering in dense ice. Phys. Rev. Lett. 89, 145501 (2002).12366053
10. Shi W Single-crystal elasticity of high-pressure ice up to 98 GPa by Brillouin Scattering Geophys. Res. Lett. 2021 48 e2021GL092514 10.1029/2021GL092514
Shi, W. et al. Single-crystal elasticity of high-pressure ice up to 98 GPa by Brillouin Scattering. Geophys. Res. Lett. 48, e2021GL092514 (2021).
11. Millot M Experimental evidence for superionic water ice using shock compression Nat. Phys. 2018 14 297 302 10.1038/s41567-017-0017-4
Millot, M. et al. Experimental evidence for superionic water ice using shock compression. Nat. Phys. 14, 297–302 (2018).
12. Weck G Evidence and stability field of fcc superionic water ice using static compression Phys. Rev. Lett. 2022 128 165701 10.1103/PhysRevLett.128.165701 35522490
Weck, G. et al. Evidence and stability field of fcc superionic water ice using static compression. Phys. Rev. Lett. 128, 165701 (2022).35522490
13. Prakapenka VB Holtgrewe N Lobanov SS Goncharov AF Structure and properties of two superionic ice phases Nat. Phys. 2021 17 1233 1238 10.1038/s41567-021-01351-8
Prakapenka, V. B., Holtgrewe, N., Lobanov, S. S. & Goncharov, A. F. Structure and properties of two superionic ice phases. Nat. Phys. 17, 1233–1238 (2021).
14. Millot M Nanosecond X-ray diffraction of shock-compressed superionic water ice Nature 2019 569 251 255 10.1038/s41586-019-1114-6 31068720
Millot, M. et al. Nanosecond X-ray diffraction of shock-compressed superionic water ice. Nature 569, 251–255 (2019).31068720
15. Queyroux J-A Melting curve and isostructural solid transition in superionic ice Phys. Rev. Lett. 2020 125 195501 10.1103/PhysRevLett.125.195501 33216588
Queyroux, J.-A. et al. Melting curve and isostructural solid transition in superionic ice. Phys. Rev. Lett. 125, 195501 (2020).33216588
16. Hubbard WB Interiors of the giant planets Science 1981 214 145 149 10.1126/science.214.4517.145 17733985
Hubbard, W. B. Interiors of the giant planets. Science 214, 145–149 (1981).17733985
17. Hubbard WB MacFarlane JJ Structure and evolution of Uranus and Neptune J. Geophys. Res.: Solid Earth 1980 85 225 234 10.1029/JB085iB01p00225
Hubbard, W. B. & MacFarlane, J. J. Structure and evolution of Uranus and Neptune. J. Geophys. Res.: Solid Earth 85, 225–234 (1980).
18. Redmer R Mattsson TR Nettelmann N French M The phase diagram of water and the magnetic fields of Uranus and Neptune Icarus 2011 211 798 803 10.1016/j.icarus.2010.08.008
Redmer, R., Mattsson, T. R., Nettelmann, N. & French, M. The phase diagram of water and the magnetic fields of Uranus and Neptune. Icarus 211, 798–803 (2011).
19. Demontis P LeSar R Klein ML New high-pressure phases of ice Phys. Rev. Lett. 1988 60 2284 2287 10.1103/PhysRevLett.60.2284 10038311
Demontis, P., LeSar, R. & Klein, M. L. New high-pressure phases of ice. Phys. Rev. Lett. 60, 2284–2287 (1988).10038311
20. Cavazzoni C Superionic and metallic states of water and ammonia at giant planet conditions Science 1999 283 44 46 10.1126/science.283.5398.44 9872734
Cavazzoni, C. et al. Superionic and metallic states of water and ammonia at giant planet conditions. Science 283, 44–46 (1999).9872734
21. Schwegler E Sharma M Gygi F Galli G Melting of ice under pressure Proc. Natl Acad. Sci. USA 2008 105 14779 14783 10.1073/pnas.0808137105 18809909
Schwegler, E., Sharma, M., Gygi, F. & Galli, G. Melting of ice under pressure. Proc. Natl Acad. Sci. USA 105, 14779–14783 (2008).18809909
22. Goldman N Fried LE Kuo I-FW Mundy CJ Bonding in the superionic phase of water Phys. Rev. Lett. 2005 94 217801 10.1103/PhysRevLett.94.217801 16090349
Goldman, N., Fried, L. E., Kuo, I.-F. W. & Mundy, C. J. Bonding in the superionic phase of water. Phys. Rev. Lett. 94, 217801 (2005).16090349
23. Wilson HF Wong ML Militzer B Superionic to superionic phase change in water: consequences for the interiors of uranus and neptune Phys. Rev. Lett. 2013 110 151102 10.1103/PhysRevLett.110.151102 25167242
Wilson, H. F., Wong, M. L. & Militzer, B. Superionic to superionic phase change in water: consequences for the interiors of uranus and neptune. Phys. Rev. Lett. 110, 151102 (2013).25167242
24. French M Desjarlais MP Redmer R Ab initio calculation of thermodynamic potentials and entropies for superionic water Phys. Rev. E 2016 93 022140 10.1103/PhysRevE.93.022140 26986321
French, M., Desjarlais, M. P. & Redmer, R. Ab initio calculation of thermodynamic potentials and entropies for superionic water. Phys. Rev. E 93, 022140 (2016).26986321
25. Sun J Clark BK Torquato S Car R The phase diagram of high-pressure superionic ice Nat. Commun. 2015 6 8156 10.1038/ncomms9156 26315260
Sun, J., Clark, B. K., Torquato, S. & Car, R. The phase diagram of high-pressure superionic ice. Nat. Commun. 6, 8156 (2015).26315260
26. Reinhardt A Thermodynamics of high-pressure ice phases explored with atomistic simulations Nat. Commun. 2022 13 4707 10.1038/s41467-022-32374-1 35948550
Reinhardt, A. et al. Thermodynamics of high-pressure ice phases explored with atomistic simulations. Nat. Commun. 13, 4707 (2022).35948550
27. Yakushev VV Postnov VI Fortov VE Yakysheva TI Electrical conductivity of water during quasi-isentropic compression to 130 GPa J. Exp. Theor. Phys. 2000 90 617 622 10.1134/1.559145
Yakushev, V. V., Postnov, V. I., Fortov, V. E. & Yakysheva, T. I. Electrical conductivity of water during quasi-isentropic compression to 130 GPa. J. Exp. Theor. Phys. 90, 617–622 (2000).
28. Chau R Mitchell AC Minich RW Nellis WJ Electrical conductivity of water compressed dynamically to pressures of 70–180 GPa (0.7–1.8 Mbar) J. Chem. Phys. 2001 114 1361 1365 10.1063/1.1332079
Chau, R., Mitchell, A. C., Minich, R. W. & Nellis, W. J. Electrical conductivity of water compressed dynamically to pressures of 70–180 GPa (0.7–1.8 Mbar). J. Chem. Phys. 114, 1361–1365 (2001).
29. Gleason AE Dynamic compression of water to conditions in ice giant interiors Sci. Rep. 2022 12 715 10.1038/s41598-021-04687-6 35027608
Gleason, A. E. et al. Dynamic compression of water to conditions in ice giant interiors. Sci. Rep. 12, 715 (2022).35027608
30. Anzellini S Boccato S A practical review of the laser-heated diamond anvil cell for university laboratories and synchrotron applications Crystals 2020 10 459 10.3390/cryst10060459
Anzellini, S. & Boccato, S. A practical review of the laser-heated diamond anvil cell for university laboratories and synchrotron applications. Crystals 10, 459 (2020).
31. Childs, C., O’Donnell, W., Ellison, P. B., Shelton, D. P. & Salamat, A. Optical and electronic solutions for power stabilization of CO2 lasers. Rev. Sci. Instrum. 91, 103003 (2020).
32. Kimura T Murakami M Revisiting the melting curve of H2O by Brillouin spectroscopy to 54 GPa J. Chem. Phys. 2023 158 134504 10.1063/5.0137943 37031134
Kimura, T. & Murakami, M. Revisiting the melting curve of H2O by Brillouin spectroscopy to 54 GPa. J. Chem. Phys. 158, 134504 (2023).37031134
33. Goncharov AF Prakapenka VB Comment on “evidence and stability field of fcc superionic water ice using static compression” Phys. Rev. Lett. 2023 131 049601 10.1103/PhysRevLett.131.049601 37566869
Goncharov, A. F. & Prakapenka, V. B. Comment on “evidence and stability field of fcc superionic water ice using static compression”. Phys. Rev. Lett. 131, 049601 (2023).37566869
34. Weck G Weck et al. Reply Phys. Rev. Lett. 2023 131 049602 10.1103/PhysRevLett.131.049602 37566873
Weck, G. et al. Weck et al. Reply. Phys. Rev. Lett. 131, 049602 (2023).37566873
35. Sugimura E Experimental evidence of superionic conduction in H2O ice J. Chem. Phys. 2012 137 194505 10.1063/1.4766816 23181324
Sugimura, E. et al. Experimental evidence of superionic conduction in H2O ice. J. Chem. Phys. 137, 194505 (2012).23181324
36. Liermann HP Novel experimental setup for megahertz X-ray diffraction in a diamond anvil cell at the High Energy Density (HED) instrument of the European X-ray Free-Electron Laser (EuXFEL) J. Synchrotron Rad. 2021 28 688 706 10.1107/S1600577521002551
Liermann, H. P. et al. Novel experimental setup for megahertz X-ray diffraction in a diamond anvil cell at the High Energy Density (HED) instrument of the European X-ray Free-Electron Laser (EuXFEL). J. Synchrotron Rad. 28, 688–706 (2021).
37. Hwang H X-ray free electron laser-induced synthesis of ε-iron nitride at high pressures J. Phys. Chem. Lett. 2021 12 3246 3252 10.1021/acs.jpclett.1c00150 33764078
Hwang, H. et al. X-ray free electron laser-induced synthesis of ε-iron nitride at high pressures. J. Phys. Chem. Lett. 12, 3246–3252 (2021).33764078
38. Husband RJ X-ray free electron laser heating of water and gold at high static pressure Commun. Mater. 2021 2 1 9 10.1038/s43246-021-00158-7
Husband, R. J. et al. X-ray free electron laser heating of water and gold at high static pressure. Commun. Mater. 2, 1–9 (2021).
39. Allahgholi A Megapixels @ megahertz – the AGIPD high-speed cameras for the European XFEL Nucl. Instrum. Methods Phys. Res. Sect. A: Accel. Spectrom. Detect. Assoc. Equip. 2019 942 162324 10.1016/j.nima.2019.06.065
Allahgholi, A. et al. Megapixels @ megahertz – the AGIPD high-speed cameras for the European XFEL. Nucl. Instrum. Methods Phys. Res. Sect. A: Accel. Spectrom. Detect. Assoc. Equip. 942, 162324 (2019).
40. Meza-Galvez J Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation J. Appl. Phys. 2020 127 195902 10.1063/1.5141360
Meza-Galvez, J. et al. Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation. J. Appl. Phys. 127, 195902 (2020).
41. Frost M Diamond precipitation dynamics from hydrocarbons at icy planet interior conditions Nat. Astron. 2024 8 174 181 10.1038/s41550-023-02147-x
Frost, M. et al. Diamond precipitation dynamics from hydrocarbons at icy planet interior conditions. Nat. Astron. 8, 174–181 (2024).
42. Ho CY Powell RW Liley PE Thermal conductivity of the elements J. Phys. Chem. Ref. Data 1972 1 279 421 10.1063/1.3253100
Ho, C. Y., Powell, R. W. & Liley, P. E. Thermal conductivity of the elements. J. Phys. Chem. Ref. Data 1, 279–421 (1972).
43. Cheng B Bethkenhagen M Pickard CJ Hamel S Phase behaviours of superionic water at planetary conditions Nat. Phys. 2021 17 1228 1232 10.1038/s41567-021-01334-9
Cheng, B., Bethkenhagen, M., Pickard, C. J. & Hamel, S. Phase behaviours of superionic water at planetary conditions. Nat. Phys. 17, 1228–1232 (2021).
44. Zastrau U The high energy density scientific instrument at the european XFEL J. Synchrotron Rad. 2021 28 1393 1416 10.1107/S1600577521007335
Zastrau, U. et al. The high energy density scientific instrument at the european XFEL. J. Synchrotron Rad. 28, 1393–1416 (2021).
45. Kujala N Hard x-ray single-shot spectrometer at the European X-ray Free-Electron Laser Rev. Sci. Instrum. 2020 91 103101 10.1063/5.0019935 33138553
Kujala, N. et al. Hard x-ray single-shot spectrometer at the European X-ray Free-Electron Laser. Rev. Sci. Instrum. 91, 103101 (2020).33138553
46. Prescher C Prakapenka VB DIOPTAS: a program for reduction of two-dimensional X-ray diffraction data and data exploration High. Press. Res. 2015 35 223 230 10.1080/08957959.2015.1059835
Prescher, C. & Prakapenka, V. B. DIOPTAS: a program for reduction of two-dimensional X-ray diffraction data and data exploration. High. Press. Res. 35, 223–230 (2015).
47. Hammersley AP FIT2D: a multi-purpose data reduction, analysis and visualization program J. Appl Cryst. 2016 49 646 652 10.1107/S1600576716000455
Hammersley, A. P. FIT2D: a multi-purpose data reduction, analysis and visualization program. J. Appl Cryst. 49, 646–652 (2016).
48. Ball OB Dynamic optical spectroscopy and pyrometry of static targets under optical and x-ray laser heating at the European XFEL J. Appl. Phys. 2023 134 055901 10.1063/5.0142196
Ball, O. B. et al. Dynamic optical spectroscopy and pyrometry of static targets under optical and x-ray laser heating at the European XFEL. J. Appl. Phys. 134, 055901 (2023).
49. Prakapenka VB Advanced flat top laser heating system for high pressure research at GSECARS: application to the melting behavior of germanium High. Press. Res. 2008 28 225 235 10.1080/08957950802050718
Prakapenka, V. B. et al. Advanced flat top laser heating system for high pressure research at GSECARS: application to the melting behavior of germanium. High. Press. Res. 28, 225–235 (2008).
50. Kim M Evidence for superionic H2O and diffusive He–H2O at high temperature and high pressure J. Phys. Condens. Matter 2022 34 394001 10.1088/1361-648X/ac8134
Kim, M. et al. Evidence for superionic H2O and diffusive He–H2O at high temperature and high pressure. J. Phys. Condens. Matter 34, 394001 (2022).
51. Frank MR Fei Y Hu J Constraining the equation of state of fluid H2O to 80 GPa using the melting curve, bulk modulus, and thermal expansivity of Ice VII1 1Associate editor: D. Sverjensky Geochim. Cosmochim. Acta 2004 68 2781 2790 10.1016/j.gca.2003.12.007
Frank, M. R., Fei, Y. & Hu, J. Constraining the equation of state of fluid H2O to 80 GPa using the melting curve, bulk modulus, and thermal expansivity of Ice VII1 1Associate editor: D. Sverjensky. Geochim. Cosmochim. Acta 68, 2781–2790 (2004).
52. Hemley RJ Static compression of H2O-ice to 128 GPa (1.28 Mbar) Nature 1987 330 737 740 10.1038/330737a0
Hemley, R. J. et al. Static compression of H2O-ice to 128 GPa (1.28 Mbar). Nature 330, 737–740 (1987).
53. Goncharov AF Dissociative melting of ice VII at high pressure J. Chem. Phys. 2009 130 124514 10.1063/1.3100771 19334858
Goncharov, A. F. et al. Dissociative melting of ice VII at high pressure. J. Chem. Phys. 130, 124514 (2009).19334858
