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

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.

Purely single-bonded spiral nitrogen chains stabilized by trivalent lanthanum ions

Inspired by the single-bonded nitrogen chains stabilized by tetravalent cerium, pentavalent tantalum, and hexavalent tungsten atoms, we explored the possibility of single-bonded nitrogen polymorphs stabilized by trivalent lanthanum ions. To achieve this, we utilized the crystal structure search method on the phase diagram of binary La-N compounds. We identified three novel thermodynamically stable phases, the C2/c LaN3, P-1 LaN4, and P-1 LaN8. Among them, the C2/c phase with infinite helical poly-N6 chains becomes thermodynamically stable above 50 GPa. Each nitrogen atom in the poly-N6 chain acquires one extra electron, and the spiral chain is purely single-bonded. The C2/c phase has an indirect band gap of ∼1.6 eV at 60 GPa. Notably, the band gap exhibits non-monotonic behavior, decreases first and then increases with increasing pressure. This abnormal behavior is attributed to the significant bonding of two La-N bonds at around 35 GPa. Phonon spectrum calculations and AIMD simulations have confirmed that the C2/c phase can be quenched to ambient conditions with slight distortion, and it exhibits excellent detonation properties. Additionally, we also discovered armchair-like nitrogen chains in LaN4 and the armchair and zigzag-like mixed nitrogen chains in LaN8. These results provide valuable insights into the electronic and bonding properties of nitrides under high pressure and may have important implications for the design and development of novel functional materials.

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. We compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

MAGUS: machine learning and graph theory assisted universal structure searcher

Crystal structure predictions based on first-principles calculations have gained great success in materials science and solid state physics. However, the remaining challenges still limit their applications in systems with a large number of atoms, especially the complexity of conformational space and the cost of local optimizations for big systems. Here, we introduce a crystal structure prediction method, MAGUS, based on the evolutionary algorithm, which addresses the above challenges with machine learning and graph theory. Techniques used in the program are summarized in detail and benchmark tests are provided. With intensive tests, we demonstrate that on-the-fly machine-learning potentials can be used to significantly reduce the number of expensive first-principles calculations, and the crystal decomposition based on graph theory can efficiently decrease the required configurations in order to find the target structures. We also summarized the representative applications of this method on several research topics, including unexpected compounds in the interior of planets and their exotic states at high pressure and high temperature (superionic, plastic, partially diffusive state, etc.); new functional materials (superhard, high-energy-density, superconducting, photoelectric materials), etc. These successful applications demonstrated that MAGUS code can help to accelerate the discovery of interesting materials and phenomena, as well as the significant value of crystal structure predictions in general.

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.

Miscibility of rock and ice in the interiors of water worlds

Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO₃), a major end-member silicate phase, and water (H₂O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO₃. Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

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.