Superionic-Superionic Phase Transitions in Body-Centered Cubic H_{2}O Ice

Dynamic-to-static switch of hydrogen bonds induces a metal-insulator transition in an organic-inorganic superlattice

Hydrogen bonds profoundly influence the fundamental chemical, physical and biological properties of molecules and materials. Owing to their relatively weaker interactions compared to other chemical bonds, hydrogen bonds alone are generally insufficient to induce substantial changes in electrical properties, thus imposing severe constraints on their applications in related devices. Here we report a metal-insulator transition controlled by hydrogen bonds for an organic-inorganic (1,3-diaminopropane)0.5SnSe2 superlattice that exhibits a colossal on-off ratio of 107 in electrical resistivity. The key to inducing the transition is a change in the amino group's hydrogen-bonding structure from dynamic to static. In the dynamic state, thermally activated free rotation continuously breaks and forms transient hydrogen bonds with adjacent Se anions. In the static state, the amino group forms three fixed-angle positions, each separated by 120°. Our findings contribute to the understanding of electrical phenomena in organic-inorganic hybrid materials and may be used for the design of future molecule-based electronic materials.

Ab initio determination of melting and sound velocity of neon up to the deep interior of the Earth

The thermophysical properties and elemental abundances of the noble gases in terrestrial materials can provide unique insights into the Earth's evolution and mantle dynamics. Here, we perform extensive ab initio molecular dynamics simulations to determine the melting temperature and sound velocity of neon up to 370 GPa and 7500 K to constrain its physical state and storage capacity, together with to reveal its implications for the deep interior of the Earth. It is found that solid neon can exist stably under the lower mantle and inner core conditions, and the abnormal melting of neon is not observed under the entire temperature (T) and pressure (P) region inside the Earth owing to its peculiar electronic structure, which is substantially distinct from other heavier noble gases. An inspection of the reduction for sound velocity along the Earth's geotherm evidences that neon can be used as a light element to account for the low-velocity anomaly and density deficit in the deep Earth. A comparison of the pair distribution functions and mean square displacements of MgSiO3-Ne and Fe-Ne alloys further reveals that MgSiO3 has a larger neon storage capacity than the liquid iron under the deep Earth condition, indicating that the lower mantle may be a natural deep noble gas storage reservoir. Our results provide valuable information for studying the fundamental behavior and phase transition of neon in a higher T-P regime, and further enhance our understanding for the interior structure and evolution processes inside the Earth.

Synthesis of large-sized spherical Co-C alloys with soft magnetic properties though a high-pressure solid-state metathesis reaction

In this work, we report a novel high-pressure solid-state metathesis (HSM) reaction to produce spherical bulk (diameters 2-4 mm) Co-C alloys (Co3C and Co1-xCx). At 2-5 GPa and 1300 °C, C atoms preferentially occupy the interstitial sites of the face-centered cubic (fcc) Co lattice, leading to the formation of metastable Pnma Co3C. The Co3C decomposes above 1400 °C at 2-5 GPa, C atoms infiltrate the interstitial sites of the fcc Co lattice, saturating the C content in Co, forming an fcc Co1-xCx solid solution while the C atoms in excess are found to precipitate in the form of graphite. The Vickers hardness of the Co-C alloys is approximately 6.1 GPa, representing a 19.6% increase compared to hexagonal close-packed (hcp) Co. First-principles calculations indicate that the presence of C atoms in the Pnma Co3C structure leads to a relative decrease in the magnetic moments of the two distinct Co atom occupancies. The Co-C alloys exhibited a soft magnetic behavior with saturation magnetization up to 93.71 emu g-1 and coercivity of 74.8 Oe; coercivity increased as the synthesis pressure rises.

Double superionicity in icy compounds at planetary interior conditions

The elements hydrogen, carbon, nitrogen and oxygen are assumed to comprise the bulk of the interiors of the ice giant planets Uranus, Neptune, and sub-Neptune exoplanets. The details of their interior structures have remained largely unknown because it is not understood how the compounds H2O, NH3 and CH4 behave and react once they have been accreted and exposed to high pressures and temperatures. Here we study thirteen H-C-N-O compounds with ab initio computer simulations and demonstrate that they assume a superionic state at elevated temperatures, in which the hydrogen ions diffuse through a stable sublattice that is provided by the larger nuclei. At yet higher temperatures, four of the thirteen compounds undergo a second transition to a novel doubly superionic state, in which the smallest of the heavy nuclei diffuse simultaneously with hydrogen ions through the remaining sublattice. Since this transition and the melting transition at yet higher temperatures are both of first order, this may introduce additional layers in the mantle of ice giant planets and alter their convective patterns.

Magnesium oxide-water compounds at megabar pressure and implications on planetary interiors

Magnesium Oxide (MgO) and water (H₂O) are abundant in the interior of planets. Their properties, and in particular their interaction, significantly affect the planet interior structure and thermal evolution. Here, using crystal structure predictions and ab initio molecular dynamics simulations, we find that MgO and H₂O can react again at ultrahigh pressure, although Mg(OH)₂ decomposes at low pressure. The reemergent MgO-H₂O compounds are: Mg₂O₃H₂ above 400 GPa, MgO₃H₄ above 600 GPa, and MgO₄H₆ in the pressure range of 270-600 GPa. Importantly, MgO₄H₆ contains 57.3 wt % of water, which is a much higher water content than any reported hydrous mineral. Our results suggest that a substantial amount of water can be stored in MgO rock in the deep interiors of Earth to Neptune mass planets. Based on molecular dynamics simulations we show that these three compounds exhibit superionic behavior at the pressure-temperature conditions as in the interiors of Uranus and Neptune. Moreover, the water-rich compound MgO₄H₆ could be stable inside the early Earth and therefore may serve as a possible early Earth water reservoir. Our findings, in the poorly explored megabar pressure regime, provide constraints for interior and evolution models of wet planets in our solar system and beyond.

Melting conditions and entropies of superionic water ice: Free-energy calculations based on hybrid solid/liquid reference systems

Superionic (SI) water ices-high-temperature, high-pressure phases of water in which oxygen ions occupy a regular crystal lattice whereas the protons flow in a liquid-like manner-have attracted a growing amount of attention over the past few years, in particular due to their possible role in the magnetic anomalies of the ice giants Neptune and Uranus. In this paper, we consider the calculation of the free energies of such phases, exploring hybrid reference systems consisting of a combination of an Einstein solid for the oxygen ions occupying a crystal lattice and a Uhlenbeck-Ford potential for the protonic fluid that avoids irregularities associated with possible particle overlaps. Applying this approach to a recent neural-network potential-energy landscape for SI water ice, we compute Gibbs free energies as a function of temperature for the SI fcc and liquid phases to determine the melting temperature T_m at 340 GPa. The results are consistent with previous estimates and indicate that the entropy difference between both phases is comparatively small, in particular due to the large amplitude of vibration of the oxygen ions in the fcc phase at the melting temperature.

Ultralow Melting Temperature of High-Pressure Face-Centered Cubic Superionic Ice

Superionic ice with oxygen in a face-centered cubic (fcc) sublattice is ascribed to the origin of magnetic fields of Uranus and Neptune, since the melting temperature (T_m) of fcc-superionic ice is believed to be higher than the isentropes of ice giants. However, precisely measuring the fcc-superionic phase experimentally remains a difficult task. The majority of the systematic investigations of its T_m were performed using perfect oxygen fcc-sublattice computations, which could result in superheating and overestimation of T_m. On the basis of the ab initio molecular dynamics method and the model with H₂O vacancy, we avoid superheating and obtain a much lower T_m than previous reports, indicating that fcc-superionic ice cannot exist in the interiors of Uranus and Neptune. Further simulations with the two-phase method justify the conclusion. The results suggest that superheating should be seriously treated when simulating the phase diagram of other hydrogen-related superionic states, which are widely used to understand the properties of ice giants, Earth, and Venus.

Thermodynamics of high-pressure ice phases explored with atomistic simulations

Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our large-scale molecular-dynamics simulations with a machine-learning potential indicate that, amongst these bcc ice phases, ices VII, VII′ and X are the same thermodynamic phase under different conditions, whereas superionic ice VII″ has a first-order phase boundary with ice VII′. Moreover, at about 300 GPa, the transformation between ice X and the Pbcm phase has a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure ices and reveals peculiar solid–solid transition mechanisms not known in other systems.

Many experimentally known high-pressure ice phase are structurally very similar. Here authors elucidate the phase behaviour of the high-pressure insulating ices and reveal solid-solid transition mechanisms not known in other systems.

Temperature- and pressure-dependence of the hydrogen bond network in plastic ice VII

We model, via classical molecular dynamics simulations, the plastic phase of ice VII across a wide range of the phase diagram of interest for planetary investigations. Although structural and dynamical properties of plastic ice VII are mostly independent on the thermodynamic conditions, the hydrogen bond network (HBN) acquires a diverse spectrum of topologies distinctly different from that of liquid water and of ice VII simulated at the same pressure. We observe that the HBN topology of plastic ice carries some degree of similarity with the crystal phase, stronger at thermodynamic conditions proximal to ice VII, and gradually lessening upon approaching the liquid state. Our results enrich our understanding of the properties of water at high pressure and high temperature, and may help in rationalizing the geology of

Stability of high-temperature salty ice suggests electrolyte permeability in water-rich exoplanet icy mantles

Electrolytes play an important role in the internal structure and dynamics of water-rich satellites and potentially water-rich exoplanets. However, in planets, the presence of a large high-pressure ice mantle is thought to hinder the exchange and transport of electrolytes between various liquid and solid deep layers. Here we show, using first-principles simulations, that up to 2.5 wt% NaCl can be dissolved in dense water ice at interior conditions of water-rich super-Earths and mini-Neptunes. The salt impurities enhance the diffusion of H atoms, extending the stability field of recently discovered superionic ice, and push towards higher pressures the transition to the stiffer ice X phase. Scaling laws for thermo-compositional convection show that salts entering the high pressure ice layer can be readily transported across. These findings suggest that the high-pressure ice mantle of water-rich exoplanets is permeable to the convective transport of electrolytes between the inner rocky core and the outer liquid layer.

Hot cubic ice is shown to retain dissolved salt in its lattice, suggesting the mantle of water-rich exoplanets is more permeable to electrolytes than assumed, which has implications on its properties and on the element cycles inside such planets.

Evidence and Stability Field of fcc Superionic Water Ice Using Static Compression

Structural transformation of hot dense water ice is investigated by combining synchrotron x-ray diffraction and a laser-heating diamond anvil cell above 25 GPa. A transition from the body-centered-cubic (bcc) to face-centered-cubic (fcc) oxygen atoms sublattices is observed from 57 GPa and 1500 K to 166 GPa and 2500 K. That is the structural signature of the transition to fcc superionic (fcc SI) ice. The sign of the density discontinuity at the transition is obtained and a phase diagram is disclosed, showing an extended fcc SI stability field. Present data also constrain the stability field of the bcc superionic (bcc SI) ice up to 100 GPa at least. The current understanding of warm dense water ice based on ab initio simulations is discussed in the light of present data.

Superionic iron alloys and their seismic velocities in Earth's inner core

Earth's inner core (IC) is less dense than pure iron, indicating the existence of light elements within it¹. Silicon, sulfur, carbon, oxygen and hydrogen have been suggested to be the candidates^2,3, and the properties of iron-light-element alloys have been studied to constrain the IC composition^4-19. Light elements have a substantial influence on the seismic velocities^4-13, the melting temperatures^14-17 and the thermal conductivities^18,19 of iron alloys. However, the state of the light elements in the IC is rarely considered. Here, using ab initio molecular dynamics simulations, we find that hydrogen, oxygen and carbon in hexagonal close-packed iron transform to a superionic state under the IC conditions, showing high diffusion coefficients like a liquid. This suggests that the IC can be in a superionic state rather than a normal solid state. The liquid-like light elements lead to a substantial reduction in the seismic velocities, which approach the seismological observations of the IC^20,21. The substantial decrease in shear-wave velocity provides an explanation for the soft IC²¹. In addition, the light-element convection has a potential influence on the IC seismological structure and magnetic field.

Superionic Silica-Water and Silica-Hydrogen Compounds in the Deep Interiors of Uranus and Neptune

Silica, water, and hydrogen are known to be the major components of celestial bodies, and have significant influence on the formation and evolution of giant planets, such as Uranus and Neptune. Thus, it is of fundamental importance to investigate their states and possible reactions under the planetary conditions. Here, using advanced crystal structure searches and first-principles calculations in the Si-O-H system, we find that a silica-water compound (SiO_{2})_{2}(H_{2}O) and a silica-hydrogen compound SiO_{2}H_{2} can exist under high pressures above 450 and 650 GPa, respectively. Further simulations reveal that, at high pressure and high temperature conditions corresponding to the interiors of Uranus and Neptune, these compounds exhibit superionic behavior, in which protons diffuse freely like liquid while the silicon and oxygen framework is fixed as solid. Therefore, these superionic silica-water and silica-hydrogen compounds could be regarded as important components of the deep mantle or core of giants, which also provides an alternative origin for their anomalous magnetic fields. These unexpected physical and chemical properties of the most common natural materials at high pressure offer key clues to understand some abstruse issues including demixing and erosion of the core in giant planets, and shed light on building reliable models for solar giants and exoplanets.

Dynamic compression of water to conditions in ice giant interiors

Recent discoveries of water-rich Neptune-like exoplanets require a more detailed understanding of the phase diagram of H₂O at pressure–temperature conditions relevant to their planetary interiors. The unusual non-dipolar magnetic fields of ice giant planets, produced by convecting liquid ionic water, are influenced by exotic high-pressure states of H₂O—yet the structure of ice in this state is challenging to determine experimentally. Here we present X-ray diffraction evidence of a body-centered cubic (BCC) structured H₂O ice at 200 GPa and ~ 5000 K, deemed ice XIX, using the X-ray Free Electron Laser of the Linac Coherent Light Source to probe the structure of the oxygen sub-lattice during dynamic compression. Although several cubic or orthorhombic structures have been predicted to be the stable structure at these conditions, we show this BCC ice phase is stable to multi-Mbar pressures and temperatures near the melt boundary. This suggests variable and increased electrical conductivity to greater depths in ice giant planets that may promote the generation of multipolar magnetic fields.

Dielectric properties of ice VII under the influence of time-alternating external electric fields

The high-pressure solid phase of water known as ice VII has recently attracted a lot of attention when its presence was detected in large exoplanets, their icy satellites, and even in Earth's mantle. Moreover, a transition of ice VII to the superionic phase can be triggered by external electric fields. Here, we investigate the dielectric responses of ice VII to applied oscillating electric fields of various frequencies employing non-equilibrium ab initio molecular dynamics. We focus on the dynamical properties of a dipole-ordered ice VII structure, for which we explored external-field-induced electronic polarisation and the vibrational spectral density of states (VDOS). These analyses are important for the understanding of collective motions in the ice-VII lattice and the electronic properties of this exotic water phase.

Superionicity, disorder, and bandgap closure in dense hydrogen chloride

Hydrogen bond networks play a crucial role in biomolecules and molecular materials such as ices. How these networks react to pressure directs their properties at extreme conditions. We have studied one of the simplest hydrogen bond formers, hydrogen chloride, from crystallization to metallization, covering a pressure range of more than 2.5 million atmospheres. Following hydrogen bond symmetrization, we identify a previously unknown phase by the appearance of new Raman modes and changes to x-ray diffraction patterns that contradict previous predictions. On further compression, a broad Raman band supersedes the well-defined excitations of phase V, despite retaining a crystalline chlorine substructure. We propose that this mode has its origin in proton (H⁺) mobility and disorder. Above 100 GPa, the optical bandgap closes linearly with extrapolated metallization at 240(10) GPa. Our findings suggest that proton dynamics can drive changes in these networks even at very high densities.

Strain-Driven Switchable Thermal Conductivity in Ferroelastic PdSe2

Materials with either high or low lattice thermal conductivity are remarkable for thermal management with applying in high-power electronic, optoelectronic, and thermoelectric devices. The realization of thermal switch between high and low thermal conductivities can greatly promote the ability of thermal energy control. Here, based on first-principles calculations, we propose that ferroelastic PdSe₂ can achieve continuous switchable thermal conductivity through strain-driven structural phase transition. Thermal switch we explored mainly stems from soft mechanical properties and strong anharmonicity of the structure after ferroelastic phase transition. We demonstrate that the maximum ratio of thermal switch can reach an order of magnitude, indicating PdSe₂ as a promising candidate in thermal devices.

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. The fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII^{''} and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

Dynamic Nature of High-Pressure Ice VII

Starting from Shannon's definition of dynamic entropy, we propose a theory to describe the rare-event-determined dynamic states in condensed matter and their transitions and apply it to high-pressure ice VII. A dynamic intensive quantity named dynamic field, rather than the conventional thermodynamic intensive quantities such as temperature and pressure, is taken as the controlling variable. The dynamic entropy versus dynamic field curve demonstrates two dynamic states in the stability region of ice VII and dynamic ice VII. Their microscopic differences were assigned to the dynamic patterns of proton transfer. This study puts a similar dynamical theory used in earlier studies of glass models on a simpler and more fundamental basis, which could be applied to describe the dynamic states of more realistic condensed matter systems.

Fluid-like elastic response of superionic NH3 in Uranus and Neptune

Nondipolar magnetic fields exhibited at Uranus and Neptune may be derived from a unique geometry of their icy mantle with a thin convective layer on top of a stratified nonconvective layer. The presence of superionic H₂O and NH₃ has been thought as an explanation to stabilize such nonconvective regions. However, a lack of experimental data on the physical properties of those superionic phases has prevented the clarification of this matter. Here, our Brillouin measurements for NH₃ show a two-stage reduction in longitudinal wave velocity (V _p) by ∼9% and ∼20% relative to the molecular solid in the temperature range of 1,500 K and 2,000 K above 47 GPa. While the first V _p reduction observed at the boundary to the superionic α phase was most likely due to the onset of the hydrogen diffusion, the further one was likely attributed to the transition to another superionic phase, denoted γ phase, exhibiting the higher diffusivity. The reduction rate of V _p in the superionic γ phase, comparable to that of the liquid, implies that this phase elastically behaves almost like a liquid. Our measurements show that superionic NH₃ becomes convective and cannot contribute to the internal stratification.

Melting Curve and Isostructural Solid Transition in Superionic Ice

The phase diagram and melting curve of water ice is investigated up to 45 GPa and 1600 K by synchrotron x-ray diffraction in the resistively and laser heated diamond anvil cell. Our melting data evidence a triple point at 14.6 GPa, 850 K. The latter is shown to be related to a first-order solid transition from the dynamically disordered form of ice VII, denoted ice VII^{'}, toward a high-temperature phase with the same bcc oxygen lattice but larger volume and higher entropy. Our experiments are compared to ab initio molecular dynamics simulations, enabling us to identify the high-temperature bcc phase with the predicted superionic ice VII^{''} phase [J.-A. Hernandez and R. Caracas, Phys. Rev. Lett. 117, 135503 (2016).PRLTAO0031-900710.1103/PhysRevLett.117.135503].

Structure and phase behavior of high-density ice from molecular-dynamics simulations with the ReaxFF potential

We report a molecular dynamics simulation study of dense ice modeled by the reactive force field (ReaxFF) potential, focusing on the possibility of phase changes between crystalline and plastic phases as observed in earlier simulation studies with rigid water models. It is demonstrated that the present model system exhibits phase transitions, or crossovers, among ice VII and two plastic ices with face-centered cubic (fcc) and body-centered cubic (bcc) lattice structures. The phase diagram derived from the ReaxFF potential is different from those of the rigid water models in that the bcc plastic phase lies on the high-pressure side of ice VII and does the fcc plastic phase on the low-pressure side of ice VII. The phase boundary between the fcc and bcc plastic phases on the pressure, temperature plane extends to the high-temperature region from the triple point of ice VII, fcc plastic, and bcc plastic phases. Proton hopping, i.e., delocalization of a proton, along between two neighboring oxygen atoms in dense ice is observed for the ReaxFF potential but only at pressures and temperatures both much higher than those at which ice VII-plastic ice transitions are observed.

Computational Prediction of Novel Ice Phases: A Perspective

Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.

Sublattice melting in binary superionic colloidal crystals

In superionic compounds one component premelts, providing high ionic conductivity to solid-state electrolytes. Here we find sublattice melting in colloidal crystals of oppositely charged particles that are highly asymmetric in size and charge in salt solutions. The small particles in ionic compounds melt when the temperature increases, forming a superionic phase. These delocalized small particles in a crystal of large oppositely charged particles, in contrast to superionic phases in atomic systems, form crystals with nonelectroneutral stoichiometric ratios. This generates structures with multiple domains of ionic crystals in percolated superionic phases with adjustable stoichiometries.

Anomalous hydrogen dynamics of the ice VII-VIII transition revealed by high-pressure neutron diffraction

Above 2 GPa the phase diagram of water simplifies considerably and exhibits only two solid phases up to 60 GPa, ice VII and ice VIII. The two phases are related to each other by hydrogen ordering, with the oxygen sublattice being essentially the same. Here we present neutron diffraction data to 15 GPa which reveal that the rate of hydrogen ordering at the ice VII-VIII transition decreases strongly with pressure to reach timescales of minutes at 10 GPa. Surprisingly, the ordering process becomes more rapid again upon further compression. We show that such an unusual change in transition rate can be explained by a slowing down of the rotational dynamics of water molecules with a simultaneous increase of translational motion of hydrogen under pressure, as previously suspected. The observed cross-over in the hydrogen dynamics in ice is likely the origin of various hitherto unexplained anomalies of ice VII in the 10-15 GPa range reported by Raman spectroscopy, X-ray diffraction, and proton conductivity.

Stability of H3O at extreme conditions and implications for the magnetic fields of Uranus and Neptune

The anomalous nondipolar and nonaxisymmetric magnetic fields of Uranus and Neptune have long challenged conventional views of planetary dynamos. A thin-shell dynamo conjecture captures the observed phenomena but leaves unexplained the fundamental material basis and underlying mechanism. Here we report extensive quantum-mechanical calculations of polymorphism in the hydrogen-oxygen system at the pressures and temperatures of the deep interiors of these ice giant planets (to >600 GPa and 7,000 K). The results reveal the surprising stability of solid and fluid trihydrogen oxide (H₃O) at these extreme conditions. Fluid H₃O is metallic and calculated to be stable near the cores of Uranus and Neptune. As a convecting fluid, the material could give rise to the magnetic field consistent with the thin-shell dynamo model proposed for these planets. H₃O could also be a major component in both solid and superionic forms in other (e.g., nonconvecting) layers. The results thus provide a materials basis for understanding the enigmatic magnetic-field anomalies and other aspects of the interiors of Uranus and Neptune. These findings have direct implications for the internal structure, composition, and dynamos of related exoplanets.

Nanosecond X-ray diffraction of shock-compressed superionic water ice

Since Bridgman's discovery of five solid water (H₂O) ice phases¹ in 1912, studies on the extraordinary polymorphism of H₂O have documented more than seventeen crystalline and several amorphous ice structures^2,3, as well as rich metastability and kinetic effects^4,5. This unique behaviour is due in part to the geometrical frustration of the weak intermolecular hydrogen bonds and the sizeable quantum motion of the light hydrogen ions (protons). Particularly intriguing is the prediction that H₂O becomes superionic^6-12-with liquid-like protons diffusing through the solid lattice of oxygen-when subjected to extreme pressures exceeding 100 gigapascals and high temperatures above 2,000 kelvin. Numerical simulations suggest that the characteristic diffusion of the protons through the empty sites of the oxygen solid lattice (1) gives rise to a surprisingly high ionic conductivity above 100 Siemens per centimetre, that is, almost as high as typical metallic (electronic) conductivity, (2) greatly increases the ice melting temperature^7-13 to several thousand kelvin, and (3) favours new ice structures with a close-packed oxygen lattice^13-15. Because confining such hot and dense H₂O in the laboratory is extremely challenging, experimental data are scarce. Recent optical measurements along the Hugoniot curve (locus of shock states) of water ice VII showed evidence of superionic conduction and thermodynamic signatures for melting¹⁶, but did not confirm the microscopic structure of superionic ice. Here we use laser-driven shockwaves to simultaneously compress and heat liquid water samples to 100-400 gigapascals and 2,000-3,000 kelvin. In situ X-ray diffraction measurements show that under these conditions, water solidifies within a few nanoseconds into nanometre-sized ice grains that exhibit unambiguous evidence for the crystalline oxygen lattice of superionic water ice. The X-ray diffraction data also allow us to document the compressibility of ice at these extreme conditions and a temperature- and pressure-induced phase transformation from a body-centred-cubic ice phase (probably ice X) to a novel face-centred-cubic, superionic ice phase, which we name ice XVIII^2,17.