Prediction of a hexagonal SiO2 phase affecting stabilities of MgSiO3 and CaSiO3 at multimegabar pressures

An ultra-short-period super-Earth with an extremely high density and an outer companion

We present the discovery and characterization of a new multi-planetary system around the Sun-like star K2-360 (EPIC 201595106). K2-360 was first identified in K2 photometry as the host of an ultra-short-period (USP) planet candidate with a period of 0.88 d. We obtained follow-up transit photometry, confirming the star as the host of the signal. High precision radial velocity measurements from HARPS and HARPS-N confirm the transiting USP planet and reveal the existence of an outer (non-transiting) planet with an orbital period of ∼ 10 d. We measure a mass of 7.67 ± 0.75   M ⊕ and a radius of 1.57 ± 0.08   R ⊕ for the transiting planet, yielding a high mean density of 11 ± 2  g  cm - 3 , making it the densest well-characterized USP super-Earth discovered to date. We measure a minimum mass of 15.2 ± 1.8   M ⊕ for the outer planet, and explore a migration formation pathway via the von Zeipel-Kozai-Lidov (ZKL) mechanism and tidal dissipation.

High Pressure-Temperature Study of MgF2, CaF2, and BaF2 by Raman Spectroscopy: Phase Transitions and Vibrational Properties of AF2 Difluorides

The phase transition of AF2 difluorides strongly depends on pressure, temperature, and cationic radius. Here, we have investigated the phase transition of three difluorides, including MgF2, CaF2, and BaF2, at simultaneously high pressures and temperatures using Raman spectroscopy and X-ray diffraction in externally heated diamond anvil cells up to 55 GPa at 300-700 K. Rutile-type difluoride MgF2 with a small cationic radius undergoes a transition to the CaCl2-type phase at 9.9(1) GPa and 300 K, to the HP-PdF2-type phase at 21.0(2) GPa, and to the cotunnite-type phase at 44.2(2) GPa. The phase transition pressure to the HP-PdF2 and cotunnite structure at 300 K for our single crystal was found to be higher than that in previous studies using polycrystalline samples. Elevating the temperature increases the transition pressure from rutile- to the CaCl2-type phase but has a negative influence on the transition pressure when MgF2 transforms from the HP-PdF2- to cotunnite-type phase. Meanwhile, the transition pressure from the CaCl2- to HP-PdF2-type phase for MgF2 was identified to be independent of the temperature. Raman peaks suspected to belong to the α-PbO2-type phase were observed at 14.6-21.0(1) GPa and 400-700 K. At 300 K, difluorides CaF2 and BaF2 in the fluorite structure with larger cationic radii transform to the cotunnite-type phase at 9.6(3) and 3.0(3) GPa at 300 K, respectively, and BaF2 further undergoes a transition to the Ni2In-type phase at 15.5(4) GPa. For both CaF2 and BaF2, elevating the temperature leads to a lower transition pressure from fluorite- to the cotunnite-type phase but has little influence on the transition to the Ni2In structure. Raman data provide valuable insights for mode Grüneisen parameters. We note that the mode Grüneisen parameters for both difluorides and dioxides vary linearly with the cation radius. Further calculations for the mode Grüneisen parameters at high pressures for MgF2, CaF2, and BaF2 yield a deeper understanding of the thermodynamic properties of the difluorides.

A new high-pressure structure of SiO2directly converted fromα-quartz under nonhydrostatic compression

High-pressure behavior of SiO₂is one of the prototypical subjects in several research areas including condensed matter physics, inorganic chemistry, mineralogy, materials science, and crystallography. Therefore, numerous studies have been performed on the structure evolution of SiO₂under pressure. Here, we show a new structure directly converted fromα-quartz under uniaxial compression. Ourab initiocalculations elucidate a simple transition pathway fromα-quartz to the Fe₂P-type phase, and an intermediate state with the Li₂ZrF₆-type structure appears in this structure conversion. Some interesting properties are found on this intermediate state. (1) The Li₂ZrF₆-type phase is metastable probably due to a volumetric unbalance between the Li and Zr sites but becomes more energetically stable thanα-quartz over ∼12 GPa. (2) It is vibrationally stable at 0 GPa, suggesting that this phase can be recovered down to ambient condition once synthesized. (3) The crystal structures of Li₂ZrF₆-type SiO₂and phase D, one of dense magnesium hydrous silicates, are found identical, suggesting the stabilization of their solid solution under high-P,Tcondition.

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.

From Two-, One-, to Zero-Dimensional Vacancies: A Densification Pattern for a Typcial Transition-Metal Dichalcogenide of TiSe2

In the areas of condensed matter physics, geoscience, material science, and inorganic chemistry, how the crystal structures evolve under an external field such as high-pressure is a fundamental question. By taking TiSe₂ as the case, we investigate the phase transformations of the layered transition-metal dichalcogenides (TMDs) under high-pressure. The ambient 6-fold P-3m1 TiSe₂ undergoes a transformation into the monoclinic 8-fold coordinated C2/m phase at 15 GPa and then into the hexagonal 9-fold Fe₂P-type structure at 34 GPa. The above phase transitions can be unitedly described as the evolution of the vacancies: from a layered structure with two-dimensional (2D) vacancies to the structure with one-dimensional (1D) and zero-dimensional (0D) vacancies. The proposed densification model of TiSe₂ reveals the processes how the symmetry breaking phase of spatial chemical bonding restores the symmetry under the isotropic external pressure.

Mixed Coordination Silica at Megabar Pressure

Silica (SiO_{2}), as a raw material of silicon, glass, ceramics, abrasive, and refractory substances, etc., is of significant importance in industrial applications and fundamental research such as electronics and planetary science. Here, using a crystal structure searching method and first-principles calculations, we predicted that a ground state crystalline phase of silica with R3[over ¯] symmetry is stable at around 645-890 GPa, which contains six-, eight-, and nine-coordinated silicon atoms and results in an average coordination number of eight. This mixed-coordination silica fills in the density, electronic band gap, and coordination number gaps between the previously known sixfold pyrite-type and ninefold Fe_{2}P-type phases, and may appear in the core or mantle of super-Earth exoplanets, or even the solar giant planets such as the Neptune. In addition, we also found that some silicon superoxides, Cmcm SiO_{3} and Ccce SiO_{6}, are stable in this pressure range and may appear in an oxygen-rich environment. Our finding enriches the high-pressure phase diagram of silicon oxides and improves understanding of the interior structure of giant planets in our solar system.

Exploring the structures and properties of nickel silicides at the pressures of the Earth's core

Given the highly possible existence of nickel and silicon in the Earth's core, the study of the reaction between Ni and Si and the resulting structures at the pressure corresponding to that of the Earth's core is highly required. Therefore, we have investigated the crystal structures of Ni-Si compounds at pressures of 0-350 GPa by adopting a crystal structure search algorithm in conjunction with first-principles calculations. We uncover two high Ni-content Ni5Si and Ni6Si compounds with 12-coordination Si bonded with Ni, with both showing strong chemical stability in the Earth's core. Bonding analysis reveals that the Ni atoms in these Ni-Si compounds present oxidant features and act as electron acceptors. This distinctive anomaly is the natural result of the energy shifts of the Ni 3d and Si 3p bands, resulting in charge transfer from Si to Ni. By examining the key properties (e.g., density and sound velocities) of the Ni5Si and Ni6Si compounds, the obtained density lies within the range of the Earth's inner core, and the estimated sound velocities are found to be consistent with seismic data. These results indicate that these two compounds could be considered as possible core constituents. Our findings provide valuable insights into the enigmatic Earth's core as well as geophysical and geochemical processes.

The Critical Point and the Supercritical State of Alkali Feldspars: Implications for the Behavior of the Crust During Impacts

The position of the vapor-liquid dome and of the critical point determine the evolution of the outermost parts of the protolunar disk during cooling and condensation after the Giant Impact. The parts of the disk in supercritical or liquid state evolve as a single thermodynamic phase; when the thermal trajectory of the disk reaches the liquid-vapor dome, gas and melt separate leading to heterogeneous convection and phase separation due to friction. Different layers of the proto-Earth behaved differently during the Giant Impact depending on their constituent materials and initial thermodynamic conditions. Here we use first-principles molecular dynamics to determine the position of the critical point for NaAlSi3O8 and KAlSi3O8 feldspars, major minerals of the Earth and Moon crusts. The variations of the pressure calculated at various volumes along isotherms yield the position of the critical points: 0.5-0.8 g cm-3 and 5500-6000 K range for the Na-feldspar, 0.5-0.9 g cm-3 and 5000-5500 K range for the K-feldspar. The simulations suggest that the vaporization is incongruent, with a degassing of O2 starting at 4000 K and gas component made mostly of free Na and K cations, O2, SiO and SiO2 species for densities below 1.5 g cm-3. The Hugoniot equations of state imply that low-velocity impactors (<8.3 km s-1) would at most melt a cold feldspathic crust, whereas large impacts in molten crust would see temperatures raise up to 30000 K.

Observation of 9-Fold Coordinated Amorphous TiO2 at High Pressure

Knowledge of the structure in amorphous dioxides is important in many fields of science and engineering. Here we report new experimental results of high-pressure polyamorphism in amorphous TiO₂ (a-TiO₂). Our data show that the Ti coordination number (CN) increases from 7.2 ± 0.3 at ∼16 GPa to 8.8 ± 0.3 at ∼70 GPa and finally reaches a plateau at 8.9 ± 0.3 at ≲86 GPa. The evolution of the structural changes under pressure is rationalized by the ratio (γ) of the ionic radius of Ti to that of O. It appears that the CN ≈ 9 plateau correlates with the two 9-fold coordinated polymorphs (cotunnite, Fe₂P) with different γ values. This CN-γ relationship is compared with those of a-SiO₂ and a-GeO₂, displaying remarkably consistent behavior between CN and γ. The unified CN-γ relationship may be generally used to predict the compression behavior of amorphous AO₂ compounds under extreme conditions.

Melting of CaSiO3 Perovskite at High Pressure

Ab initio molecular dynamics simulations predict that CaSiO3 perovskite melts at 5600 K at 136 GPa, and 6400 K at 300 GPa, significantly higher than MgSiO3 perovskite. The entropy of melting (1.8 kB per atom) is much larger than that of many silicates at ambient pressure and of simple liquids and varies little with pressure. The volume of melting decreases rapidly with increasing pressure, to 3 % at 136 GPa, producing a melting slope that diminishes rapidly with pressure. We determine the melting temperature via the ZW method, combining the Z method, for which we clarify the theoretical basis, with a waiting time analysis. The ZW method results are internally confirmed by integrating the Clausius-Clapeyron equation, which also yields our results for the entropy and volume of melting. We find the eutectic composition on the MgSiO3-CaSiO3 join to be x Ca = 0.26 at 136 GPa and that metasilicate melt is denser than coexisting silicates.

Electrical conductivity and magnetic dynamos in magma oceans of Super-Earths

Super-Earths are extremely common among the numerous exoplanets that have been discovered. The high pressures and temperatures in their interiors are likely to lead to long-lived magma oceans. If their electrical conductivity is sufficiently high, the mantles of Super-Earth would generate their own magnetic fields. With ab initio simulations, we show that upon melting, the behavior of typical mantle silicates changes from semi-conducting to semi-metallic. The electrical conductivity increases and the optical properties are substantially modified. Melting could thus be detected with high-precision reflectivity measurements during the short time scales of shock experiments. We estimate the electrical conductivity of mantle silicates to be of the order of 100 Ω⁻¹ cm⁻¹, which implies that a magnetic dynamo process would develop in the magma oceans of Super-Earths if their convective velocities have typical values of 1 mm/s or higher. We predict exoplanets with rotation periods longer than 2 days to have multipolar magnetic fields.

With the discovery of large rocky exoplanets called Super-Earths, questions have arisen regarding the properties of their interiors and their ability to produce a magnetic field. Here, the authors show that under high pressure, molten silicates are semi-metallic and that magma oceans would host a dynamo process.

Structural, magnetic and electronic properties of CrO2 at multimegabar pressures

As the only half-metallic ferromagnetic material in 3d transition metal dioxides, CrO2 has attracted great scientific interest from materials science to physical chemistry. Here, an investigation into the structural, magnetic and electronic properties of CrO2 under high pressure has been conducted by first-principles calculations based on density functional theory. Static calculations have predicted that CrO2 undergoes structural transitions with the sequence of rutile-type → CaCl2-type → pyrite-type → Pnma → (Fe2P-type→) I4/mmm at high pressures. In addition, a transition from the ferromagnetic state to the non-magnetic state with the magnetic collapse of Cr is observed in CrO2 at the pyrite-Pnma transition. This transition also delocalizes the 3d electrons of Cr and leads to a metallic character of CrO2. The equation of state, elasticity and band gap for each energetically favorable phase of CrO2 are determined. Our results not only bridge the gap about the high-pressure behavior of CrO2 in previous studies but also extend our understanding of its properties up to multimegabar conditions. According to previous data and present results, we further discuss and summarize the high-pressure behavior of various AO2 compounds. This can contribute to investigating properties of other AO2 compounds or exploring novel materials at high pressures.

Correlated High-Pressure Phase Sequence of VO2 under Strong Compression

Understanding how the structures of a crystal behave under compression is a fundamental issue both for condensed matter physics and for geoscience. Traditional description of a crystal as the stacking of a unit cell with special symmetry has gained much success on the analysis of physical properties. Unfortunately, it is hard to reveal the relationship between the compressed phases. Taking the family of metal dioxides (MO₂) as an example, the structural evolution, subject to fixed chemical formula and highly confined space, often appears as a set of random and uncorrelated events. Here we provide an alternative way to treat the crystal as the stacking of the coordination polyhedron and then discover a unified structure transition pattern, in our case VO₂. X-ray diffraction (XRD) experiments and first-principles calculations show that the coordination increase happens only at one apex of the V-centered octahedron in an orderly fashion, leaving the base plane and the other apex topologically intact. The polyhedron evolves toward increasing their sharing, indicating a general rule for the chemical bonds of MO₂ to give away the ionicity in exchange for covalency under pressure.

The melting points of MgO up to 4 TPa predicted based on ab initio thermodynamic integration molecular dynamics

The melting curve of MgO is extended up to 4 TPa, corresponding to the Jovian core pressure, based on the one-step thermodynamic integration method implemented on ab initio molecular dynamics. The calculated melting temperatures are 3100 and 16 000 K at 0 and 500 GPa, respectively, which are consistent with previous experimental results, and 20 600 K at 3900 GPa, which is inconsistent with a recent experimental extrapolation, which implies the molten Jovian core. A quite small Clapeyron slope ([Formula: see text]) of [Formula: see text] is found at 3900 GPa due to comparable densities of the liquid and B2 phases under extreme compression. The Mg-O coordination number in the liquid phase is saturated at around 7.5 above 1 TPa and remains smaller than that in the B2 phase (8) even at 4 TPa, suggesting no density crossover between liquid and crystal and thus no further denser crystalline phases. Dynamical properties (atomic diffusivity and viscosity) are also investigated along the melting curve to understand these behaviors in greater detail.

Beyond sixfold coordinated Si in SiO2 glass at ultrahigh pressures

We investigated the structure of SiO2 glass up to 172 GPa using high-energy X-ray diffraction. The combination of a multichannel collimator with diamond anvil cells enabled the measurement of structural changes in silica glass with total X-ray diffraction to previously unachievable pressures. We show that SiO2 first undergoes a change in Si-O coordination number from fourfold to sixfold between 15 and 50 GPa, in agreement with previous investigations. Above 50 GPa, the estimated coordination number continuously increases from 6 to 6.8 at 172 GPa. Si-O bond length shows first an increase due to the fourfold to sixfold coordination change and then a smaller linear decrease up to 172 GPa. We reconcile the changes in relation to the oxygen-packing fraction, showing that oxygen packing decreases at ultrahigh pressures to accommodate the higher than sixfold Si-O coordination. These results give experimental insight into the structural changes of silicate glasses as analogue materials for silicate melts at ultrahigh pressures.

The phase diagrams of KCaF3 and NaMgF3 by ab initio simulations

ABF3 compounds have been found to make valuable low-pressure analogues for high-pressure silicate phases that are present in the Earth's deep interior and that may also occur in the interiors of exoplanets. The phase diagrams of two of these materials, KCaF3 and NaMgF3, have been investigated in detail by static ab initio computer simulations based on density functional theory. Six ABF3 polymorphs were considered, as follows: the orthorhombic perovskite structure (GdFeO3-type; space group Pbnm); the orthorhombic CaIrO3 structure (Cmcm; commonly referred to as the "post-perovskite" structure); the orthorhombic Sb2S3 and La2S3 structures (both Pmcn); the hexagonal structure previously suggested in computer simulations of NaMgF3 (P63/mmc); the monoclinic structure found to be intermediate between the perovskite and CaIrO3 structures in CaRhO3 (P21/m). Volumetric and axial equations of state of all phases considered are presented. For KCaF3, as expected, the perovskite phase is shown to be the most thermodynamically stable at atmospheric pressure. With increasing pressure, the relative stability of the KCaF3 phases then follows the sequence: perovskite → La2S3 structure → Sb2S3 structure → P63/mmc structure; the CaIrO3 structure is never the most stable form. Above about 2.6 GPa, however, none of the KCaF3 polymorphs are stable with respect to dissociation into KF and CaF2. The possibility that high-pressure KCaF3 polymorphs might exist metastably at 300 K, or might be stabilised by chemical substitution so as to occur within the standard operating range of a multi-anvil press, is briefly discussed. For NaMgF3, the transitions to the high-pressure phases occur at pressures outside the normal range of a multi-anvil press. Two different sequences of transitions had previously been suggested from computer simulations. With increasing pressure, we find that the relative stability of the NaMgF3 phases follows the sequence: perovskite → CaIrO3 structure → Sb2S3 structure → P63/mmc structure. However, only the perovskite and CaIrO3 structures are stable with respect to dissociation into NaF and MgF2.

Exploring the free energy surface using ab initio molecular dynamics

Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti.

High-pressure structural changes in liquid silica

The structural properties of liquid silica at high pressure and moderate temperature conditions, also referred to as the warm dense matter regime, were investigated using time-resolved K-edge x-ray absorption spectroscopy and ab initio calculations. We used a nanosecond laser beam to compress uniformly a solid SiO_{2} target and a picosecond laser beam to generate a broadband x-ray source. We obtained x-ray absorption spectra at the Si K edge over a large pressure-temperature domain to probe the liquid phase up to 3.6 times the normal solid density. Using ab initio simulations, we are able to interpret the changes in the x-ray absorption near-edge structure with increasing densities as an increase in the coordination number of silicon by oxygen atoms from 4 to 9. This indicates that, up to significant temperatures, the liquid structure becomes akin to what is found in the solid SiO_{2} phases.

Melting curve of SiO2 at multimegabar pressures: implications for gas giants and super-Earths

Ultrahigh-pressure phase boundary between solid and liquid SiO₂ is still quite unclear. Here we present predictions of silica melting curve for the multimegabar pressure regime, as obtained from first principles molecular dynamics simulations. We calculate the melting temperatures from three high pressure phases of silica (pyrite-, cotunnite-, and Fe₂P-type SiO₂) at different pressures using the Z method. The computed melting curve is found to rise abruptly around 330 GPa, an increase not previously reported by any melting simulations. This is in close agreement with recent experiments reporting the α-PbO₂–pyrite transition around this pressure. The predicted phase diagram indicates that silica could be one of the dominant components of the rocky cores of gas giants, as it remains solid at the core of our Solar System’s gas giants. These results are also relevant to model the interior structure and evolution of massive super-Earths.

Experimental and theoretical thermal equations of state of MgSiO3 post-perovskite at multi-megabar pressures

The MgSiO₃ post-perovskite phase is the most abundant silicate phase in a super-Earth’s mantle, although it only exists within the Earth’s lowermost mantle. In this study, we established the thermal equation of state (EoS) of the MgSiO₃ post-perovskite phase, which were determined by using both laser-heated diamond anvil cell and density-functional theoretical techniques, within a multi-megabar pressure range, corresponding to the conditions of a super-Earth’s mantle. The Keane and AP2 EoS models were adopted for the first time to extract meaningful physical properties. The experimentally determined Grüneisen parameter, which is one of the thermal EoS parameters, and its volume dependence were found to be consistent with their theoretically obtained values. This reduced the previously reported discrepancy observed between experiment and theory. Both the experimental and theoretical EoS were also found to be in very good agreement for volumes at pressures and temperatures of up to 300 GPa and 5000 K, respectively. Our newly developed EoS should be applicable to a super-Earth’s mantle, as well as the Earth’s core-mantle boundary region.

Prediction of novel stable compounds in the Mg-Si-O system under exoplanet pressures

The Mg-Si-O system is the major Earth and rocky planet-forming system. Here, through quantum variable-composition evolutionary structure explorations, we have discovered several unexpected stable binary and ternary compounds in the Mg-Si-O system. Besides the well-known SiO₂ phases, we have found two extraordinary silicon oxides, SiO₃ and SiO, which become stable at pressures above 0.51 TPa and 1.89 TPa, respectively. In the Mg-O system, we have found one new compound, MgO₃, which becomes stable at 0.89 TPa. We find that not only the (MgO)_x·(SiO₂)_y compounds, but also two (MgO₃)_x·(SiO₃)_y compounds, MgSi₃O₁₂ and MgSiO₆, have stability fields above 2.41 TPa and 2.95 TPa, respectively. The highly oxidized MgSi₃O₁₂ can form in deep mantles of mega-Earths with masses above 20 M_⊕ (M_⊕:Earth’s mass). Furthermore, the dissociation pathways of pPv-MgSiO₃ are also clarified, and found to be different at low and high temperatures. The low-temperature pathway is MgSiO₃ ⇒ Mg₂SiO₄ + MgSi₂O₅ ⇒ SiO₂ + Mg₂SiO₄ ⇒ MgO + SiO₂, while the high-temperature pathway is MgSiO₃ ⇒ Mg₂SiO₄ + MgSi₂O₅ ⇒ MgO + MgSi₂O₅ ⇒ MgO + SiO₂. Present results are relevant for models of the internal structure of giant exoplanets, and for understanding the high-pressure behavior of materials.

Laser-shock compression of magnesium oxide in the warm-dense-matter regime

Magnesium oxide has been experimentally and computationally investigated in the warm-dense solid and liquid ranges from 200 GPa to 1 TPa along the principal Hugoniot. The linear approximation between shock velocity and particle velocity is validated up to a shock velocity of 15 km/s from the experimental data, this suggesting that the MgO B1 structure is stable up to the corresponding shock pressure of ∼350 GPa. Moreover, our Hugoniot data, combined with ab initio simulations, show two crossovers between MgO Hugoniot and the extrapolation of the linear approximation line, occurring at a shock pressures of approximately 350 and 650 GPa, with shock temperatures of 8000 and 14,000 K, respectively. These crossover regions are consistent with the solid-solid (B1-B2) and the solid-liquid (B2-melt) phase boundaries predicted by the ab initio calculations.

Prediction of 10-fold coordinated TiO2 and SiO2 structures at multimegabar pressures

We predict by first-principles methods a phase transition in TiO2 at 6.5 Mbar from the Fe2P-type polymorph to a ten-coordinated structure with space group I4/mmm. This is the first report, to our knowledge, of the pressure-induced phase transition to the I4/mmm structure among all dioxide compounds. The I4/mmm structure was found to be up to 3.3% denser across all pressures investigated. Significant differences were found in the electronic properties of the two structures, and the metallization of TiO2 was calculated to occur concomitantly with the phase transition to I4/mmm. The implications of our findings were extended to SiO2, and an analogous Fe2P-type to I4/mmm transition was found to occur at 10 TPa. This is consistent with the lower-pressure phase transitions of TiO2, which are well-established models for the phase transitions in other AX2 compounds, including SiO2. As in TiO2, the transition to I4/mmm corresponds to the metallization of SiO2. This transformation is in the pressure range reached in the interiors of recently discovered extrasolar planets and calls for a reformulation of the equations of state used to model them.

Planetary science. Shock compression of stishovite and melting of silica at planetary interior conditions

Deep inside planets, extreme density, pressure, and temperature strongly modify the properties of the constituent materials. In particular, how much heat solids can sustain before melting under pressure is key to determining a planet's internal structure and evolution. We report laser-driven shock experiments on fused silica, α-quartz, and stishovite yielding equation-of-state and electronic conductivity data at unprecedented conditions and showing that the melting temperature of SiO2 rises to 8300 K at a pressure of 500 gigapascals, comparable to the core-mantle boundary conditions for a 5-Earth mass super-Earth. We show that mantle silicates and core metal have comparable melting temperatures above 500 to 700 gigapascals, which could favor long-lived magma oceans for large terrestrial planets with implications for planetary magnetic-field generation in silicate magma layers deep inside such planets.

Metallization of warm dense SiO(2) studied by XANES spectroscopy

We investigate the evolution of the electronic structure of fused silica in a dense plasma regime using time-resolved x-ray absorption spectroscopy. We use a nanosecond (ns) laser beam to generate a strong uniform shock wave in the sample and a picosecond (ps) pulse to produce a broadband x-ray source near the Si K edge. By varying the delay between the two laser beams and the intensity of the ns beam, we explore a large thermodynamical domain with densities varying from 1 to 5  g/cm^{3} and temperatures up to 5 eV. In contrast to normal conditions where silica is a well-known insulator with a wide band gap of 8.9 eV, we find that shocked silica exhibits a pseudogap as a semimetal throughout this thermodynamical domain. This is in quantitative agreement with density functional theory predictions performed using the generalized gradient approximation.

An adaptive genetic algorithm for crystal structure prediction

We present a genetic algorithm (GA) for structural search that combines the speed of structure exploration by classical potentials with the accuracy of density functional theory (DFT) calculations in an adaptive and iterative way. This strategy increases the efficiency of the DFT-based GA by several orders of magnitude. This gain allows a considerable increase in the size and complexity of systems that can be studied by first principles. The performance of the method is illustrated by successful structure identifications of complex binary and ternary intermetallic compounds with 36 and 54 atoms per cell, respectively. The discovery of a multi-TPa Mg-silicate phase with unit cell containing up to 56 atoms is also reported. Such a phase is likely to be an essential component of terrestrial exoplanetary mantles.

Theoretical and experimental evidence for a new post-cotunnite phase of titanium dioxide with significant optical absorption

We report the discovery of a post-cotunnite phase of TiO2 by both density-functional ab initio calculations and high-pressure experiments. A pressure-induced phase transition to a hexagonal Fe2P-type structure (space group P62m) was predicted to occur at 161 GPa and 0 K and successfully observed by in situ synchrotron x-ray diffraction measurements at 210 GPa and 4000 K with a significant increase in opacity. This change in opacity is attributed to a reduction of band gap from 3.0 to 1.9 eV across the phase change. The Fe2P-type structure is proved to be the densest phase in major metal dioxides.