Superionic and metallic states of water and ammonia at giant planet conditions

On the growth rate of ices: Effect of pressure and ice phase

Ice, the solid phase of water, can adopt a wide range of structures, making it of interest from both fundamental and applied perspectives. In this study, we used molecular dynamics simulations to compute the growth rates of four ice phases: Ih, III, V, and VI. To enable comparisons at the same temperature, different pressures were applied to each phase. Our analysis of pressure effects on the growth rate of ice Ih revealed only a minor influence, allowing us to attribute variations in growth rates primarily to structural differences among ice phases. We observed that ices Ih and VI exhibit similar growth rates, whereas ices III and V grow significantly faster. Rapidly growing ice phases exhibit a high growth efficiency, meaning that less molecular motion is needed to form the solid. We hypothesize that proton ordering may influence ice growth, as partially ordered ice phases (III and V) exhibit faster growth rates than fully disordered phases (Ih and VI). Alternative explanations for ice growth rate trends, such as unit cell complexity or melting entropy, are ruled out. Finally, we assessed the predictive capability of the Wilson-Frenkel theory. While the theory does not inherently account for structural complexity-resulting in similar growth rates across all phases-we found that introducing a phase-specific characteristic length enables it to accurately reproduce our simulation results.

Observation of plastic ice VII by quasi-elastic neutron scattering

Water is the third most abundant molecule in the universe and a key component in the interiors of icy moons, giant planets and Uranus- and Neptune-like exoplanets1-3. Owing to its distinct molecular structure and flexible hydrogen bonds that readily adapt to a wide range of pressures and temperatures, water forms numerous crystalline and amorphous phases4-6. Most relevant for the high pressures and temperatures of planetary interiors is ice VII (ref. 4), and simulations have identified along its melting curve the existence of a so-called plastic phase7-12 in which individual molecules occupy fixed positions as in a solid yet are able to rotate as in a liquid. Such plastic ice has not yet been directly observed in experiments. Here we present quasi-elastic neutron scattering measurements, conducted at temperatures between 450 and 600 K and pressures up to 6 GPa, that reveal the existence of a body-centred cubic structure, as found in ice VII, with water molecules showing picosecond rotational dynamics typical for liquid water. Comparison with molecular dynamics simulations indicates that this plastic ice VII does not conform to a free rotor phase but rather shows rapid orientational jumps, as observed in jump-rotor plastic crystals13,14. We anticipate that our observation of plastic ice VII will affect our understanding of the geodynamics of icy planets and the differentiation processes of large icy moons.

Water uptake of solids and its impact on ion transport

The interaction modes of water with (polar) solids are manifold, comprising surface adsorption and incorporation into the bulk, both in molecular and in dissociated form. This Review discusses these processes and the respective pronounced effects on the ionic transport properties. The concentration as well as the mobility of ionic carriers can vary by orders of magnitude depending on the water content on or within a solid. Selected materials examples, which are relevant for electrochemical devices (for example, low- and intermediate-temperature fuel cells) or which are of fundamental interest (such as molecular water acting as dopant in a lithium halide), are treated in more detail. Interrelations between hydration and electronic defects are also briefly touched upon.

Superconducting and superionic behaviors of electride Na6C under moderate pressure

Electrides exhibiting diverse electride states, superconductivity, and superionic behavior have attracted considerable attention for elucidating intricate chemical bonding and particle interactions. However, due to the scarcity of electrides exhibiting three-dimensional electride states, the associated properties abovementioned remain elusive. Here, we propose an electride C2/m-Na6C that exhibits peculiar three-dimensional electride states under 30 GPa and maintains dynamic stability at ambient conditions. Our electron-phonon interaction calculations reveal that the T c of 0.051 K at 30 GPa arises from the scattering Na/C sp-hybridized electrons by Na/C-derived low-frequency phonons, rather than three-dimensional electride states. Further molecular dynamics simulations indicate that this structure exhibits superionic states at 1,200 K, where the unusually heavy sodium atoms exhibit diffusive behavior. This anomalous behavior can be attributed to the duality of formation process of three-dimensional electride states and their multi-center bonding effect. Our study provides theoretical guidance for further investigation into the diverse physical characteristics of electrides.

Stability of Proton Superoxide and its Superionic Transition Under High Pressure

Under extreme conditions, condensed matters are subject to undergo a phase transition and there have been many attempts to find another form of hydroxide stabilized over H2O. Here, using Density Functional Theory (DFT)-based crystal structure prediction including zero-point energy, it is that proton superoxide (HO2), the lightest superoxide, can be stabilized energetically at high pressure and temperature conditions. HO2 is metallic at high pressure, which originates from the 𝜋* orbitals overlap between adjacent superoxide anions (O2 -). By lowering pressure, it undergoes a metal-to-insulator transition similar to LiO2. Ab initio molecular dynamics (AIMD) calculations reveal that HO2 becomes superionic with high electrical conductivity. The possibility of creating hydrogen-mixed superoxide at lower pressure using a (Lix,H1-x)O2 hypothetical structure is also proposed. This discovery bridges gaps in superoxide and superionicity, guiding the design of various H-O compounds under high pressure.

X-Ray Signature of the Superionic Transition in Warm Dense fcc Water Ice

The fcc superionic phase of ice is a key component of the warm dense water phase diagram. While a few x-ray diffraction studies, under dynamic and static compressions, have reported the stability of the fcc structure, the transition to the superionic state has not been investigated in detail. Here, a remarkable thermal volume expansion is disclosed, which is interpreted as being directly related to the superionic transition. This could be achieved by implementing a heating capsule geometry within the laser-heated diamond anvil cell. Fcc ice is recovered metastable at ambient temperature, allowing us to observe that superionicity in the fcc phase emerges at a slightly lower temperature than for the bcc-fcc structural transition. The crossover in volume thermal expansion at the superionic transition agrees with recent ab initio calculations; however, its magnitude is larger than predicted.

A dynamical system approach to relaxation in glass-forming liquids

The "classical" thermodynamic and statistical mechanical theories of Gibbs and Boltzmann are both predicated on axiomatic assumptions whose applicability is hard to ascertain. Theoretical objections and an increasing number of observed deviations from these theories have led to sustained efforts to develop an improved mathematical and physical foundation for them, and the search for appropriate extensions that are generally applicable to condensed materials at low temperatures (T) and high material densities where the assumptions of these theories start to become particularly questionable. These theoretical efforts have largely focused on minimal models of condensed material systems, such as the Fermi-Ulam-Pasta-Tsingou model, and other simplified models of condensed materials that are amenable to numerical and analytic treatments and that can serve to illuminate essential features of relaxation processes in condensed materials under conditions approaching integrable dynamics where clear departures from classical thermodynamics and dynamics can be generally expected. These studies indicate an apparently general multi-step relaxation process, corresponding to an initial "fast" relaxation process (termed the fast β-relaxation in the context of cooled liquids), followed by a longer "equipartition time", namely, the α-relaxation time τα in the context of cooled liquids. This relaxation timescale can be enormously longer than the fast β-relaxation time τβ so that τα is the primary parameter governing the rate at which the material comes into equilibrium, and thus is a natural focus of theoretical attention. Since the dynamics of these simplified dynamical systems, originally intended as simplified models of real crystalline materials exhibiting anharmonic interactions, greatly resemble the observed relaxation dynamics of both heated crystals and cooled liquids, we adapt this dynamical system approach to the practical matter of estimating relaxation times in both cooled liquids and crystals at elevated temperatures, which we identify as weakly non-integrable dynamical systems.

New metallic ice phase under high pressure

Crystal materials can exhibit novel properties under high pressure, which are completely different from properties under ambient conditions. Water ice has an exceptionally rich phase diagram with at least 20 known crystalline ice phases from experiments, where the high-pressure ice X and ice XVIII behave as an ionic state and a superionic state, respectively. Thus, the ice structures stabilized under high pressure are very likely to possess other novel properties. Herein, we constructed a sequence of hypothetical high-pressure ices whose structures were generated according to the topological frameworks of common metal oxides. Based on density functional theory calculations, the pressure-induced phase transition sequence is in order that the known Ag2O-Pn3̄m structure (ice X) transformed into a previously undiscovered TiO2_brookite-Pbca structure at a pressure of 300 GPa, followed by a transition to a new NaO2-Pa3 structure at a pressure of 2120 GPa. Hitherto unreported NaO2-Pa3 ice with a cubic structure is in the ionic state, where the oxygen atoms in NaO2-Pa3 have a face-centered cubic (fcc) sublattice, and the coordination number of H atoms increases to 3. These two structures are dynamically stable according to phonon spectrum analysis and remain stable at temperature of 100 K based on ab initio molecular dynamics (AIMD) simulations. More importantly, the NaO2-Pa3 ice exhibits novel metallic properties with a closing band gap above a pressure of 2600 GPa, owing to the electron orbital coupling of oxygen atoms in close proximity induced by pressure.

Phase transition kinetics of superionic H2O ice phases revealed by Megahertz X-ray free-electron laser-heating experiments

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.

Orientational Disorder Drives Site Disorder in Plastic Ammonia Hemihydrate

In the 2-10 GPa pressure range, ammonia hemihydrate H_{2}O:(NH_{3})_{2} (AHH) is a molecular solid in which intermolecular interactions are ruled by distinct types of hydrogen bonds. Upon heating, the low-temperature ordered P2_{1}/c crystal (AHH-II) transits to a bcc phase (AHH-pbcc) where each site is randomly occupied by water or ammonia. In addition to the site disorder, experiments suggest that AHH-pbcc is a plastic solid, but the physical origin and mechanisms at play for the rotational and site disordering remain unknown. Using large-scale (∼10^{5}  atoms) and long-time (>10  ns) simulations, we show that, as temperature rises above the transition line, orientational disorder sets in, breaking the strongest hydrogen bonds that provide the largest contribution to the cohesion of the ordered AHH-II phase and enabling the molecules to migrate from a crystal site to a neighboring one. This generates a plastic molecular alloy with site disorder while the solid state is overall maintained until melting at a higher temperature. The case of high (P,T) plastic ammonia hemihydrate can be extended to other water-ammonia alloys where a similar interplay between distinct hydrogen bonds occurs.

Superionic iron hydride shapes ultralow-velocity zones at Earth's core-mantle boundary

Seismological studies have exposed numerous ultralow velocity zones (ULVZs) exhibiting extraordinary physical attributes at Earth's core-mantle boundary, yet their compositions and origins remain controversial. Water-iron reaction can generate unique phases under lowermost-mantle conditions and likely plays a crucial role in forming ULVZs. Through first-principles molecular dynamic simulations with machine learning techniques, we determine that iron hydride, the product of water-iron reaction, is stable as a superionic phase at the core-mantle boundary. This superionic iron hydride has much slower velocities and a higher density than the ambient mantle under lowermost-mantle conditions. Accumulation of iron hydride, created through either a chemical reaction between subducted water and iron or solidification of core material entrained in the lower mantle by convection, can explain the seismic observations of ULVZs particularly those associated with subduction. This work suggests that water may have a substantial role in creating seismic heterogeneities at the core-mantle boundary.

Proton Collective Quantum Tunneling Induces Anomalous Thermal Conductivity of Ice under Pressure

Proton tunneling is believed to be nonlocal in ice, but its range has been shown to be limited to only a few molecules. Here, we measured the thermal conductivity of ice under pressure up to 50 GPa and found it increases with pressure until 20 GPa but decreases at higher pressures. We attribute this nonmonotonic thermal conductivity to the collective tunneling of protons at high pressures, supported by large-scale quantum molecular dynamics simulations. The collective tunneling loops span several picoseconds in time and are as large as nanometers in space, which match the phonon periods and wavelengths, leading to strong phonon scattering at high pressures. Our results show direct evidence of global quantum motion existing in high-pressure ice and provide a new perspective to understanding the coupling between phonon propagation and atomic tunneling.

Mixture of hydrogen and methane under planetary interior conditions

We employ first-principles molecular dynamics simulations to provide equation-of-state data, pair distribution functions (PDFs), diffusion coefficients, and band gaps of a mixture of hydrogen and methane under planetary interior conditions as relevant for Uranus, Neptune, and similar icy exoplanets. We test the linear mixing approximation, which is fulfilled within a few percent for the chosen P-T conditions. Evaluation of the PDFs reveals that methane molecules dissociate into carbon clusters and free hydrogen atoms at temperatures greater than 3000 K. At high temperatures, the clusters are found to be short-lived. Furthermore, we calculate the electrical conductivity from which we derive the non-metal-to-metal transition region of the mixture. We also calculate the electrical conductivity along the P-T profile of Uranus [N. Nettelmann et al., Planet. Space Sci., 2013, 77, 143-151] and observe the transition of the mixture from a molecular to an atomic fluid as a function of the radius of the planet. The density and temperature ranges chosen in our study can be achieved using dynamic shock compression experiments and seek to aid such future experiments. Our work also provides a relevant data set for a better understanding of the interior, evolution, luminosity, and magnetic field of the ice giants in our solar system and beyond.

Raman and IR spectra of water under graphene nanoconfinement at ambient and extreme pressure-temperature conditions: a first-principles study

The nanoconfinement of water can result in dramatic differences in its physical and chemical properties compared to bulk water. However, a detailed molecular-level understanding of these properties is still lacking. Vibrational spectroscopy, such as Raman and infrared, is a popular experimental tool for studying the structure and dynamics of water, and is often complemented by atomistic simulations to interpret experimental spectra, but there have been few theoretical spectroscopy studies of nanoconfined water using first-principles methods at ambient conditions, let alone under extreme pressure-temperature conditions. Here, we compute the Raman and IR spectra of water nanoconfined by graphene at ambient and extreme pressure-temperature conditions using ab initio simulations. Our results revealed alterations in the Raman stretching and low-frequency bands due to the graphene confinement. We also found spectroscopic evidence indicating that nanoconfinement considerably changes the tetrahedral hydrogen bond network, which is typically found in bulk water. Furthermore, we observed an unusual bending band in the Raman spectrum at ∼10 GPa and 1000 K, which is attributed to the unique molecular structure of confined ionic water. Additionally, we found that at ∼20 GPa and 1000 K, confined water transformed into a superionic fluid, making it challenging to identify the IR stretching band. Finally, we computed the ionic conductivity of confined water in the ionic and superionic phases. Our results highlight the efficacy of Raman and IR spectroscopy in studying the structure and dynamics of nanoconfined water in a large pressure-temperature range. Our predicted Raman and IR spectra can serve as a valuable guide for future experiments.

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.

Electrolyte Permeability in Plastic Ice VII

Deep brines in water-rich planets form when electrolytes diffuse from the rocky interior through layers of thick dense ice such as ice VII and the hypothesized plastic ice VII. We perform classical molecular dynamics simulations of Li⁺, Na⁺, and K⁺ alkali ions and F^- and Cl^- halide ions in plastic ice VII at conditions similar to water-rich super-Earths, icy moons, and ocean worlds. We find that plastic ice VII is permeable to electrolytes on geological timescales. Diffusion occurs via jumps between adjacent voids in the bcc crystal structure and is governed by molecular rotations. An exception to this mechanism is Na⁺ which, at variance with other ions, can substitute water molecules on lattice positions. The bulk modulus of pristine plastic ice VII is dependent on the pace of molecular rotations: when the rotations are slow, the bulk modulus is 1 order of magnitude lower compared to the bulk modulus at conditions of fast rotations, hence providing direct evidence of the role of molecular rotations in determining elastic properties. Electrolytes affect the bulk modulus only at high-concentration conditions and slow molecular rotations. Our results show that plastic ice VII may facilitate the development of brines in water-rich planets and ocean worlds, with a clear significance for their potential to support exobiology and for the chemical evolution of their aqueous reservoirs.

Colloidal superionic conductors

Nanoparticles with highly asymmetric sizes and charges that self-assemble into crystals via electrostatics may exhibit behaviors reminiscent of those of metals or superionic materials. Here, we use coarse-grained molecular simulations with underdamped Langevin dynamics to explore how a binary charged colloidal crystal reacts to an external electric field. As the field strength increases, we find transitions from insulator (ionic state), to superionic (conductive state), to laning, to complete melting (liquid state). In the superionic state, the resistivity decreases with increasing temperature, which is contrary to metals, yet the increment decreases as the electric field becomes stronger. Additionally, we verify that the dissipation of the system and the fluctuation of charge currents obey recently developed thermodynamic uncertainty relation. Our results describe charge transport mechanisms in colloidal superionic conductors.

Molecular rotations trigger a glass-to-plastic fcc heterogeneous crystallization in high-pressure water

We report a molecular dynamics study of the heterogeneous crystallization of high-pressure glassy water using (plastic) ice VII as a substrate. We focus on the thermodynamic conditions P ∈ [6-8] GPa and T ∈ [100-500] K, at which (plastic) ice VII and glassy water are supposed to coexist in several (exo)planets and icy moons. We find that (plastic) ice VII undergoes a martensitic phase transition to a (plastic) fcc crystal. Depending on the molecular rotational lifetime τ, we identify three rotational regimes: for τ > 20 ps, crystallization does not occur; for τ ∼ 15 ps, we observe a very sluggish crystallization and the formation of a considerable amount of icosahedral environments trapped in a highly defective crystal or in the residual glassy matrix; and for τ < 10 ps, crystallization takes place smoothly, resulting in an almost defect-free plastic fcc solid. The presence of icosahedral environments at intermediate τ is of particular interest as it shows that such a geometry, otherwise ephemeral at lower pressures, is, indeed, present in water. We justify the presence of icosahedral structures based on geometrical arguments. Our results represent the first study of heterogeneous crystallization occurring at thermodynamic conditions of relevance for planetary science and unveil the role of molecular rotations in achieving it. Our findings (i) show that the stability of plastic ice VII, widely reported in the literature, should be reconsidered in favor of plastic fcc, (ii) provide a rationale for the role of molecular rotations in achieving heterogeneous crystallization, and (iii) represent the first evidence of long-living icosahedral structures in water. Therefore, our work pushes forward our understanding of the properties of water.

Observation of a Plastic Crystal in Water-Ammonia Mixtures under High Pressure and Temperature

Solid mixtures of ammonia and water, the so-called ammonia hydrates, are thought to be major components of solar and extra-solar icy planets. We present here a thorough characterization of the recently reported high pressure (P)-temperature (T) phase VII of ammonia monohydrate (AMH) using Raman spectroscopy, X-ray diffraction, and quasi-elastic neutron scattering (QENS) experiments in the ranges 4-10 GPa, 450-600 K. Our results show that AMH-VII exhibits common structural features with the disordered ionico-molecular alloy (DIMA) phase, stable above 7.5 GPa at 300 K: both present a substitutional disorder of water and ammonia over the sites of a body-centered cubic lattice and are partially ionic. The two phases however markedly differ in their hydrogen dynamics, and QENS measurements show that AMH-VII is characterized by free molecular rotations around the lattice positions which are quenched in the DIMA phase. AMH-VII is thus a peculiar crystalline solid in that it combines three types of disorder: substitutional, compositional, and rotational.

Cooperative diffusion in body-centered cubic iron in Earth and super-Earths' inner core conditions

The physical chemistry of iron at the inner-core conditions is key to understanding the evolution and habitability of Earth and super-Earth planets. Based on full first-principles simulations, we report cooperative diffusion along the longitudinally fast⟨111⟩directions of body-centered cubic (bcc) iron in temperature ranges of up to 2000-4000 K below melting and pressures of ∼300-4000 GPa. The diffusion is due to the low energy barrier in the corresponding direction and is accompanied by mechanical and dynamical stability, as well as strong elastic anisotropy of bcc iron. These findings provide a possible explanation for seismological signatures of the Earth's inner core, particularly the positive correlation between P wave velocity and attenuation. The diffusion can also change the detailed mechanism of core convection by increasing the diffusivity and electrical conductivity and lowering the viscosity. The results need to be considered in future geophysical and planetary models and should motivate future studies of materials under extreme conditions.

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.

Imaging velocity interferometer system for any reflector (VISAR) diagnostics for high energy density sciences

Two variants of optical imaging velocimetry, specifically the one-dimensional streaked line-imaging and the two-dimensional time-resolved area-imaging versions of the Velocity Interferometer System for Any Reflector (VISAR), have become important diagnostics in high energy density sciences, including inertial confinement fusion and dynamic compression of condensed matter. Here, we give a brief review of the historical development of these techniques, then describe the current implementations at major high energy density (HED) facilities worldwide, including the OMEGA Laser Facility and the National Ignition Facility. We illustrate the versatility and power of these techniques by reviewing diverse applications of imaging VISARs for gas-gun and laser-driven dynamic compression experiments for materials science, shock physics, condensed matter physics, chemical physics, plasma physics, planetary science and astronomy, as well as a broad range of HED experiments and laser-driven inertial confinement fusion research.

In operando measurements of high explosives

In operando studies of high explosives involve dynamic extreme conditions produced as a shock wave travels through the explosive to produce a detonation. Here, we describe a method to safely produce detonations and dynamic extreme conditions in high explosives and in inert solids and liquids on a tabletop in a high-throughput format. This method uses a shock compression microscope, a microscope with a pulsed laser that can launch a hypervelocity flyer plate along with a velocimeter, an optical pyrometer, and a nanosecond camera that together can measure pressures, densities, and temperatures with high time and space resolution (2 ns and 2 µm). We discuss how a detonation builds up in liquid nitromethane and show that we can produce and study detonations in sample volumes close to the theoretical minimum. We then discuss how a detonation builds up from a shock in a plastic-bonded explosive (PBX) based on HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane), where the initial steps are hotspot formation and deflagration growth in the shocked microstructure. A method is demonstrated where we can measure thermal emission from high-temperature reactions in every HMX crystal in the PBX, with the intent of determining which configurations produce the critical hot spots that grow and ignite the entire PBX.

Liquid-like atoms in dense-packed solid glasses

Revealing the microscopic structural and dynamic pictures of glasses is a long-standing challenge for scientists^1,2. Extensive studies on the structure and relaxation dynamics of glasses have constructed the current classical picture^3-5: glasses consist of some 'soft zones' of loosely bound atoms embedded in a tightly bound atomic matrix. Recent experiments have found an additional fast process in the relaxation spectra^6-9, but the underlying physics of this process remains unclear. Here, combining extensive dynamic experiments and computer simulations, we reveal that this fast relaxation is associated with string-like diffusion of liquid-like atoms, which are inherited from the high-temperature liquids. Even at room temperature, some atoms in dense-packed metallic glasses can diffuse just as easily as they would in liquid states, with an experimentally determined viscosity as low as 10⁷ Pa·s. This finding extends our current microscopic picture of glass solids and might help establish the dynamics-property relationship of glasses⁴.

Electrofreezing of Liquid Ammonia

Here we prove that, in addition to temperature and pressure, another important thermodynamic variable permits the exploration of the phase diagram of ammonia: the electric field. By means of (path integral) ab initio molecular dynamics simulations, we predict that, upon applying intense electric fields on ammonia, the electrofreezing phenomenon occurs, leading the liquid toward a novel ferroelectric solid phase. This study proves that electric fields can generally be exploited as the access key to otherwise-unreachable regions in phase diagrams, unveiling the existence of new condensed-phase structures. Furthermore, the reported findings have manifold practical implications, from the safe storage and transportation of ammonia to the understanding of the solid structures this compound forms in planetary contexts.

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

Evidence for superionic H2O and diffusive He-H2O at high temperature and high pressure

We present the evidence of superionic phase formed in H₂O and, for the first time, diffusive H₂O-He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H₂O and diffusive He-H₂O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H₂O phase forms gradually in pure H₂O over the temperature range of 1350-1400 K at 23 GPa, but the diffusive He-H₂O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He-H₂O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He-H₂O.

The pressure induced phase diagram of double-layer ice under confinement: a first-principles study

Here, we present double-layer ice confined within various carbon nanotubes (CNTs) using state-of-the-art pressure induced (-5 GPa to 5 GPa) dispersion corrected density functional theory (DFT) calculations. We find that the double-layer ice exhibits remarkably rich and diverse phase behaviors as a function of pressure with varying CNT diameters. The lattice cohesive energies for various pure double layer ice phases follow the order of hexagonal > pentagonal > square tube > hexagonal-close-packed (HCP) > square > buckled-rhombic (b-RH). The confinement width was found to play a crucial role in the square and square tube phases in the intermediate pressure range of about 0-1 GPa. Unlike the phase transition in pure bilayer ice structures, the relative enthalpies demonstrate that the pentagonal phase, rather than the hexagonal structure, is the most stable ice polymorph at ambient pressure as well as in a deep negative pressure region, whereas the b-RH phase dominates under high pressure. The relatively short O⋯O distance of b-RH demonstrates the presence of a strong hydrogen bonding network, which is responsible for stabilizing the system.

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.

Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials

The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely, the fluorite-structure materials PbF₂, CaF₂, and UO₂ (type II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several properties, at increasing size and in a wide temperature range. This work shows that, for the considered type-II superionic conductors, the diffusivity dramatically depends on the system size and that the superionic regime is shifted to larger temperatures in smaller cells. Furthermore, only simulations of several hundred atoms are able to capture the experimentally observed, characteristic change in the activation energy of the diffusion process, occurring at the order-disorder transition to the superionic regime. Finite-size effects on ion diffusion are instead much weaker in α-AgI. The thermal conductivity is found generally smaller for smaller cells, where the temperature-independent (Allen-Feldman) regime is also reached at significantly lower temperatures. The finite-size effects on the thermal motion of the non-mobile ions composing the solid matrix follow the simple law that holds for solids.

Superionic states formation in group III oxides irradiated with ultrafast lasers

After ultrafast laser irradiation, a target enters a poorly explored regime where physics of a solid state overlaps with plasma physics and chemistry, creating an unusual synergy—a warm dense matter state (WDM). We study theoretically the WDM kinetics and chemistry in a number of group III-metal oxides with highly excited electronic system. We employ density functional theory to investigate a possibility of nonthermal transition of the materials into a superionic state under these conditions. Atomic and electronic properties of the materials are analyzed during the transitions to acquire insights into physical mechanisms guiding such transformations.

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.

Ultrahigh-Pressure Magnesium Hydrosilicates as Reservoirs of Water in Early Earth

The origin of water on the Earth is a long-standing mystery, requiring a comprehensive search for hydrous compounds, stable at conditions of the deep Earth and made of Earth-abundant elements. Previous studies usually focused on the current range of pressure-temperature conditions in the Earth's mantle and ignored a possible difference in the past, such as the stage of the core-mantle separation. Here, using ab initio evolutionary structure prediction, we find that only two magnesium hydrosilicate phases are stable at megabar pressures, α-Mg_{2}SiO_{5}H_{2} and β-Mg_{2}SiO_{5}H_{2}, stable at 262-338 GPa and >338  GPa, respectively (all these pressures now lie within the Earth's iron core). Both are superionic conductors with quasi-one-dimensional proton diffusion at relevant conditions. In the first 30 million years of Earth's history, before the Earth's core was formed, these must have existed in the Earth, hosting much of Earth's water. As dense iron alloys segregated to form the Earth's core, Mg_{2}SiO_{5}H_{2} phases decomposed and released water. Thus, now-extinct Mg_{2}SiO_{5}H_{2} phases have likely contributed in a major way to the evolution of our planet.

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.

Faraday law, oxidation numbers, and ionic conductivity: The role of topology

Faraday's experiment measures-within a modern view-the charge adiabatically transported over a macroscopic distance by a given nuclear species in insulating liquids: the reason why it is an integer is deeply rooted in topology. Whole numbers enter chemistry in a different form: atomic oxidation states. They are not directly measurable and are determined instead from an agreed set of rules. Insulating liquids are a remarkable exception; Faraday's experiment indeed measures the oxidation numbers of each dissociated component in the liquid phase, whose topological values are unambiguous. Ionic conductivity in insulating liquids is expressed in terms of the autocorrelation function of the fluctuating charge current at a given temperature in a zero electric field; topology plays a major role in this important observable as well. The existing literature deals with the above issues by adopting the independent-electron framework; here, I provide the many-body generalization of all the above findings, which, furthermore, allows for compact and very transparent notations and formulas.

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.

Pressure-Induced Superionicity of H- in Hypervalent Sodium Silicon Hydrides

Superionic states simultaneously exhibit properties of a fluid and a solid. Proton (H⁺) superionicity in ice, H₃O, He-H₂O, and He-NH₃ compounds is well-studied. However, hydride (H^-) superionicity in H-rich compounds is rare, being associated with instability and strongly reducing conditions. Silicon, sodium, and hydrogen are abundant elements in many astrophysical bodies. Here, we use first-principles calculations to show that, at high pressure, Na, Si, and H can form several hypervalent compounds. A previously unreported superionic state of Na₂SiH₆ results from unconstrained H^- in the hypervalent [SiH₆]^2- unit. Na₂SiH₆ is dynamically stable at low pressure (3 GPa), becoming superionic at 5 GPa, and re-entering solid/fluid states at about 25 GPa. Our observation of H^- transport opens up a new field of H^- conductors. It also has implications for the formation of conducting layers at depth in exotic carbon exoplanets, potentially enhancing the habitability of such planets.

Spectroscopic evidence for a gold-coloured metallic water solution

Insulating materials can in principle be made metallic by applying pressure. In the case of pure water, this is estimated¹ to require a pressure of 48 megabar, which is beyond current experimental capabilities and may only exist in the interior of large planets or stars^2-4. Indeed, recent estimates and experiments indicate that water at pressures accessible in the laboratory will at best be superionic with high protonic conductivity⁵, but not metallic with conductive electrons¹. Here we show that a metallic water solution can be prepared by massive doping with electrons upon reacting water with alkali metals. Although analogous metallic solutions of liquid ammonia with high concentrations of solvated electrons have long been known and characterized^6-9, the explosive interaction between alkali metals and water^10,11 has so far only permitted the preparation of aqueous solutions with low, submetallic electron concentrations^12-14. We found that the explosive behaviour of the water-alkali metal reaction can be suppressed by adsorbing water vapour at a low pressure of about 10^-4 millibar onto liquid sodium-potassium alloy drops ejected into a vacuum chamber. This set-up leads to the formation of a transient gold-coloured layer of a metallic water solution covering the metal alloy drops. The metallic character of this layer, doped with around 5 × 10²¹ electrons per cubic centimetre, is confirmed using optical reflection and synchrotron X-ray photoelectron spectroscopies.

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.

First principles study of dense and metallic nitric sulfur hydrides

Studies of molecular mixtures containing hydrogen sulfide (H₂S) could open up new routes towards hydrogen-rich high-temperature superconductors under pressure. H₂S and ammonia (NH₃) form hydrogen-bonded molecular mixtures at ambient conditions, but their phase behavior and propensity towards mixing under pressure is not well understood. Here, we show stable phases in the H₂S-NH₃ system under extreme pressure conditions to 4 Mbar from first-principles crystal structure prediction methods. We identify four stable compositions, two of which, (H₂S) (NH₃) and (H₂S) (NH₃)₄, are stable in a sequence of structures to the Mbar regime. A re-entrant stabilization of (H₂S) (NH₃)₄ above 300 GPa is driven by a marked reversal of sulfur-hydrogen chemistry. Several stable phases exhibit metallic character. Electron-phonon coupling calculations predict superconducting temperatures up to 50 K, in the Cmma phase of (H₂S) (NH₃) at 150 GPa. The present findings shed light on how sulfur hydride bonding and superconductivity are affected in molecular mixtures. They also suggest a reservoir for hydrogen sulfide in the upper mantle regions of icy planets in a potentially metallic mixture, which could have implications for their magnetic field formation.

Electronic-structure methods for materials design

The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.

Crystal Structure of Symmetric Ice X in H2O-H2 and H2O-He under Pressure

Ice VII and ice X are the two most dominant phases, stable over a large pressure range between 2 and 150 GPa and made of fundamentally different chemical bonding. Yet, the two ice phases share a similar bcc-based crystal structure and lattice constants, resulting in a challenge to discern the crystal structure of ice VII and ice X. Here, we present well-resolved X-ray diffraction data of H₂O in quasi-hydrostatic H₂ and He pressure media, clearly resolving the two ice phases to 130 GPa and the dissociative nature of ice VII to X transition occurring at 20-50 GPa in H₂O-H₂ and 60-70 GPa in H₂O-He. The present diffraction data permits, for the first time, the accurate determination of the bulk moduli B₀ of 225 (or 228) GPa for ice X and 6.2 (or 4.5) GPa for ice VII, in H₂O-H₂ (or H₂O-He), which can provide new constraints for Giant planetary models.

Water-Gas Shift Reaction Produces Formate at Extreme Pressures and Temperatures in Deep Earth Fluids

The water-gas shift reaction is one of the most important reactions in industrial hydrogen production and plays a key role in Fischer-Tropsch-type synthesis, which is widely believed to generate hydrocarbons in the deep carbon cycle but is little known at extreme pressure-temperature conditions found in the Earth's upper mantle. Here, we performed extensive ab initio molecular dynamics simulations and free energy calculations to study the water-gas shift reaction. We found the direct formation of formic acid from CO and supercritical water at 10-13 GPa and 1400 K without any catalyst. Contrary to the common assumption that formic acid or formate is an intermediate product, we found that HCOOH is thermodynamically more stable than the products of the water-gas shift reaction above 3 GPa and at 1000-1400 K. Our study suggests that the water-gas shift reaction may not happen in the Earth's upper mantle, and formic acid or formate may be an important carbon carrier in reducing environments, participating in many geochemical processes in deep Earth.

Unusual Chemistry of the C-H-N-O System under Pressure and Implications for Giant Planets

C-H-N-O system is central for organic chemistry and biochemistry and plays a major role in planetary science (dominating the composition of "ice giants" Uranus and Neptune). The inexhaustible chemical diversity of this system at normal conditions explains its role as the basis of all known life, but the chemistry of this system at high pressures and temperatures of planetary interiors is poorly known. Using ab initio evolutionary algorithm USPEX, we performed an extensive study of the phase diagram of the C-H-N-O system at pressures of 50, 200, and 400 GPa and temperatures up to 3000 K. Seven novel thermodynamically stable phases were predicted, including quaternary polymeric crystal C₂H₂N₂O₂ and several new N-O and H-N-O compounds. We describe the main patterns of changes in the chemistry of the C-H-N-O system under pressure and confirm that diamond should be formed at conditions of the middle-ice layers of Uranus and Neptune. We also provide the detailed CH₄-NH₃-H₂O phase diagrams at high pressures, which are important for further improvement of the models of ice giants, and point out that current models are clearly deficient. In particular, in the existing models, Uranus and Neptune are assumed to have identical composition, nearly identical pressure-temperature profiles, and a single convecting middle layer ("mantle") made of a mixture of H₂O/CH₄/NH₃ in the ratio of 56.5:32.5:11. Here, we provide new insights, shedding light into the difference of heat flows from Uranus and Neptune, which require them to have different compositions, pressure-temperature conditions, and a more complex internal structure.

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.