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

Melting of B1‐Phase MgO From Simultaneous True Radiative Shock Temperature and Sound Speed Measurements to 250 GPa on Samples Preheated to 2300 K

To refine the melting curve, equation of state, and physical properties of MgO we performed plate impact experiments spanning 170–250 GPa on MgO single crystals, preheated to 2300 K. A controlled thermal gradient in 20 mm long samples enabled radiative temperature (3%–4%) and rarefaction overtake observations (yielding sound speed 2%) close to the hot Mo driver with a free surface below 2000 K that minimized evaporation. Ta flyers were launched by two‐stage light‐gas gun up to 7.6 km/s and sample radiance was recorded with a 6‐channel (500–850 nm) pyrometer. Shock front reflectivity was measured at 198 and 243 GPa using 50/50 sapphire beam‐splitters. Most experiments show monotonic increases of shock temperature with pressure, from (168 GPa, 7100 K) to (243 GPa, 9400 K), in good agreement with predictions of our MgO B1 phase equation of state. Measured sound speeds are parallel to but 10% higher than model predictions for bulk sound speed of solid B1 MgO, confirming ductile behavior of preheated MgO. Two experiments, at 238 and 246 GPa, showed anomalously low shock temperature and sound speed, suggesting melting. Using reported MgO melting data up to 120 GPa and our data at 232–246 GPa, we constructed a maximum‐likelihood Simon‐Glatzel fit. At Earth's core‐mantle boundary pressure (135 GPa), our best‐fit interpolated MgO melting point is K. Our proposed melting line falls within the envelope of theoretical predictions but does not completely agree with any particular model curve. Our results reduce the uncertainty on MgO melting temperature at Earth's core‐mantle boundary by a factor of 17 and provide an anchor for extension to multicomponent systems.

B1-B2 transition in shock-compressed MgO

Magnesium oxide (MgO) is a major component of the Earth's mantle and is expected to play a similar role in the mantles of large rocky exoplanets. At extreme pressures, MgO transitions from the NaCl B1 crystal structure to a CsCl B2 structure, which may have implications for exoplanetary deep mantle dynamics. In this study, we constrain the phase diagram of MgO with laser-compression along the shock Hugoniot, with simultaneous measurements of crystal structure, density, pressure, and temperature. We identify the B1 to B2 phase transition between 397 and 425 gigapascal (around 9700 kelvin), in agreement with recent theory that accounts for phonon anharmonicity. From 425 to 493 gigapascal, we observe a mixed-phase region of B1 and B2 coexistence. The transformation follows the Watanabe-Tokonami-Morimoto mechanism. Our data are consistent with B2-liquid coexistence above 500 gigapascal and complete melting at 634 gigapascal. This study bridges the gap between previous theoretical and experimental studies, providing insights into the timescale of this phase transition.

Davemaoite as the mantle mineral with the highest melting temperature

Knowledge of high-pressure melting curves of silicate minerals is critical for modeling the thermal-chemical evolution of rocky planets. However, the melting temperature of davemaoite, the third most abundant mineral in Earth's lower mantle, is still controversial. Here, we investigate the melting curves of two minerals, MgSiO3 bridgmanite and CaSiO3 davemaoite, under their stability field in the mantle by performing first-principles molecular dynamics simulations based on the density functional theory. The melting curve of bridgmanite is in excellent agreement with previous studies, confirming a general consensus on its melting temperature. However, we predict a much higher melting curve of davemaoite than almost all previous estimates. Melting temperature of davemaoite at the pressure of core-mantle boundary (~136 gigapascals) is about 7700(150) K, which is approximately 2000 K higher than that of bridgmanite. The ultrarefractory nature of davemaoite is critical to reconsider many models in the deep planetary interior, for instance, solidification of early magma ocean and geodynamical behavior of mantle rocks.

Extended X-ray absorption fine structure of dynamically-compressed copper up to 1 terapascal

Large laser facilities have recently enabled material characterization at the pressures of Earth and Super-Earth cores. However, the temperature of the compressed materials has been largely unknown, or solely relied on models and simulations, due to lack of diagnostics under these challenging conditions. Here, we report on temperature, density, pressure, and local structure of copper determined from extended x-ray absorption fine structure and velocimetry up to 1 Terapascal. These results nearly double the highest pressure at which extended x-ray absorption fine structure has been reported in any material. In this work, the copper temperature is unexpectedly found to be much higher than predicted when adjacent to diamond layer(s), demonstrating the important influence of the sample environment on the thermal state of materials; this effect may introduce additional temperature uncertainties in some previous experiments using diamond and provides new guidance for future experimental design.

Measuring the principal Hugoniot of inertial-confinement-fusion-relevant TMPTA plastic foams

Wetted-foam layers are of significant interest for inertial-confinement-fusion capsules, due to the control they provide over the convergence ratio of the implosion and the opportunity this affords to minimize hydrodynamic instability growth. However, the equation of state for fusion-relevant foams are not well characterized, and many simulations rely on modeling such foams as a homogeneous medium with the foam average density. To address this issue, an experiment was performed using the VULCAN Nd:glass laser at the Central Laser Facility. The aim was to measure the principal Hugoniot of TMPTA plastic foams at 260mg/cm^{3}, corresponding to the density of liquid DT-wetted-foam layers, and their "hydrodynamic equivalent" capsules. A VISAR was used to obtain the shock velocity of both the foam and an α-quartz reference layer, while streaked optical pyrometry provided the temperature of the shocked material. The measurements confirm that, for the 20-120 GPa pressure range accessed, this material can indeed be well described using the equation of state of the homogeneous medium at the foam density.

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.

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.

Metallic Aluminum Suboxides with Ultrahigh Electrical Conductivity at High Pressure

Aluminum, as the most abundant metallic elemental content in the Earth's crust, usually exists in the form of alumina (Al₂O₃). However, the oxidation state of aluminum and the crystal structures of aluminum oxides in the pressure range of planetary interiors are not well established. Here, we predicted two aluminum suboxides (Al₂O, AlO) and two superoxides (Al₄O₇, AlO₃) with uncommon stoichiometries at high pressures using first-principle calculations and crystal structure prediction methods. We find that the P4/nmm Al₂O becomes stable above ~765 GPa and may survive in the deep mantles or cores of giant planets such as Neptune. Interestingly, the Al₂O and AlO are metallic and have electride features, in which some electrons are localized in the interstitials between atoms. We find that Al₂O has an electrical conductivity one order of magnitude higher than that of iron under the same pressure-temperature conditions, which may influence the total conductivity of giant planets. Our findings enrich the high-pressure phase diagram of aluminum oxides and improve our understanding of the interior structure of giant planets.

Liquid structure under extreme conditions: high-pressure x-ray diffraction studies

Under extreme conditions of high pressure and temperature, liquids can undergo substantial structural transformations as their atoms rearrange to minimise energy within a more confined volume. Understanding the structural response of liquids under extreme conditions is important across a variety of disciplines, from fundamental physics and exotic chemistry to materials and planetary science.In situexperiments and atomistic simulations can provide crucial insight into the nature of liquid-liquid phase transitions and the complex phase diagrams and melting relations of high-pressure materials. Structural changes in natural magmas at the high-pressures experienced in deep planetary interiors can have a profound impact on their physical properties, knowledge of which is important to inform geochemical models of magmatic processes. Generating the extreme conditions required to melt samples at high-pressure, whilst simultaneously measuring their liquid structure, is a considerable challenge. The measurement, analysis, and interpretation of structural data is further complicated by the inherent disordered nature of liquids at the atomic-scale. However, recent advances in high-pressure technology mean that liquid diffraction measurements are becoming more routinely feasible at synchrotron facilities around the world. This topical review examines methods for high pressure synchrotron x-ray diffraction of liquids and the wide variety of systems which have been studied by them, from simple liquid metals and their remarkable complex behaviour at high-pressure, to molecular-polymeric liquid-liquid transitions in pnicogen and chalcogen liquids, and density-driven structural transformations in water and silicate melts.

The High Energy Density Scientific Instrument at the European XFEL

The European XFEL delivers up to 27000 intense (>1012 photons) pulses per second, of ultrashort (≤50 fs) and transversely coherent X-ray radiation, at a maximum repetition rate of 4.5 MHz. Its unique X-ray beam parameters enable groundbreaking experiments in matter at extreme conditions at the High Energy Density (HED) scientific instrument. The performance of the HED instrument during its first two years of operation, its scientific remit, as well as ongoing installations towards full operation are presented. Scientific goals of HED include the investigation of extreme states of matter created by intense laser pulses, diamond anvil cells, or pulsed magnets, and ultrafast X-ray methods that allow their diagnosis using self-amplified spontaneous emission between 5 and 25 keV, coupled with X-ray monochromators and optional seeded beam operation. The HED instrument provides two target chambers, X-ray spectrometers for emission and scattering, X-ray detectors, and a timing tool to correct for residual timing jitter between laser and X-ray pulses.

Measuring the structure and equation of state of polyethylene terephthalate at megabar pressures

We present structure and equation of state (EOS) measurements of biaxially orientated polyethylene terephthalate (PET, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({\hbox {C}}_{10} {\hbox {H}}_8 {\hbox {O}}_4)_n$$\end{document}(C10H8O4)n, also called mylar) shock-compressed to (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$155 \pm 20$$\end{document}155±20) GPa and (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6000 \pm 1000$$\end{document}6000±1000) K using in situ X-ray diffraction, Doppler velocimetry, and optical pyrometry. Comparing to density functional theory molecular dynamics (DFT-MD) simulations, we find a highly correlated liquid at conditions differing from predictions by some equations of state tables, which underlines the influence of complex chemical interactions in this regime. EOS calculations from ab initio DFT-MD simulations and shock Hugoniot measurements of density, pressure and temperature confirm the discrepancy to these tables and present an experimentally benchmarked correction to the description of PET as an exemplary material to represent the mixture of light elements at planetary interior conditions.

Temperature and Density on the Forsterite Liquid-Vapor Phase Boundary

The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid-vapor phase boundary. The entire shock-and-release thermodynamic path must be calculated accurately in order to predict the post-impact structures of planetary bodies. Forsterite (Mg2SiO4) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid-vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under-predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.

Melting and density of MgSiO3 determined by shock compression of bridgmanite to 1254GPa

The essential data for interior and thermal evolution models of the Earth and super-Earths are the density and melting of mantle silicate under extreme conditions. Here, we report an unprecedently high melting temperature of MgSiO3 at 500 GPa by direct shockwave loading of pre-synthesized dense MgSiO3 (bridgmanite) using the Z Pulsed Power Facility. We also present the first high-precision density data of crystalline MgSiO3 to 422 GPa and 7200 K and of silicate melt to 1254 GPa. The experimental density measurements support our density functional theory based molecular dynamics calculations, providing benchmarks for theoretical calculations under extreme conditions. The excellent agreement between experiment and theory provides a reliable reference density profile for super-Earth mantles. Furthermore, the observed upper bound of melting temperature, 9430 K at 500 GPa, provides a critical constraint on the accretion energy required to melt the mantle and the prospect of driving a dynamo in massive rocky planets.

Electrical conductivity of warm dense silica from double-shock experiments

Understanding materials behaviour under extreme thermodynamic conditions is fundamental in many branches of science, including High-Energy-Density physics, fusion research, material and planetary science. Silica (SiO₂) is of primary importance as a key component of rocky planets' mantles. Dynamic compression is the most promising approach to explore molten silicates under extreme conditions. Although most experimental studies are restricted to the Hugoniot curve, a wider range of conditions must be reached to distill temperature and pressure effects. Here we present direct measurements of equation of state and two-colour reflectivity of double-shocked α-quartz on a large ensemble of thermodynamic conditions, which were until now unexplored. Combining experimental reflectivity data with numerical simulations we determine the electrical conductivity. The latter is almost constant with pressure while highly dependent on temperature, which is consistent with simulations results. Based on our findings, we conclude that dynamo processes are likely in Super-Earths' mantles.

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.

Metallization of Shock-Compressed Liquid Ammonia

Ammonia is predicted to be one of the major components in the depths of the ice giant planets Uranus and Neptune. Their dynamics, evolution, and interior structure are insufficiently understood and models rely imperatively on data for equation of state and transport properties. Despite its great significance, the experimentally accessed region of the ammonia phase diagram today is still very limited in pressure and temperature. Here we push the probed regime to unprecedented conditions, up to ∼350  GPa and ∼40 000  K. Along the Hugoniot, the temperature measured as a function of pressure shows a subtle change in slope at ∼7000  K and ∼90  GPa, in agreement with ab initio simulations we have performed. This feature coincides with the gradual transition from a molecular liquid to a plasma state. Additionally, we performed reflectivity measurements, providing the first experimental evidence of electronic conduction in high-pressure ammonia. Shock reflectance continuously rises with pressure above 50 GPa and reaches saturation values above 120 GPa. Corresponding electrical conductivity values are up to 1 order of magnitude higher than in water in the 100 GPa regime, with possible significant contributions of the predicted ammonia-rich layers to the generation of magnetic dynamos in ice giant interiors.

Neptune and Uranus: ice or rock giants?

Existing observations of Uranus and Neptune's fundamental physical properties can be fitted with a wide range of interior models. A key parameter in these models is the bulk rock:ice ratio and models broadly fall into ice-dominated (ice giant) and rock-dominated (rock giant) categories. Here we consider how observations of Neptune's atmospheric temperature and composition (H2, He, D/H, CO, CH4, H2O and CS) can provide further constraints. The tropospheric CO profile in particular is highly diagnostic of interior ice content, but is also controversial, with deep values ranging from zero to 0.5 parts per million. Most existing CO profiles imply extreme O/H enrichments of >250 times solar composition, thus favouring an ice giant. However, such high O/H enrichment is not consistent with D/H observations for a fully mixed and equilibrated Neptune. CO and D/H measurements can be reconciled if there is incomplete interior mixing (ice giant) or if tropospheric CO has a solely external source and only exists in the upper troposphere (rock giant). An interior with more rock than ice is also more compatible with likely outer solar system ice sources. We primarily consider Neptune, but similar arguments apply to Uranus, which has comparable C/H and D/H enrichment, but no observed tropospheric CO. While both ice- and rock-dominated models are viable, we suggest a rock giant provides a more consistent match to available atmospheric observations. This article is part of a discussion meeting issue 'Future exploration of ice giant systems'.

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Time-resolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 PPa (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. The techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petapascal range, and can be used to deduce the energy flow in a wide variety of implosion systems to petapascal energy densities.

Green chemical mechanical polishing of sapphire wafers using a novel slurry

Toxic and corrosive solutions are widely used in the preparation of abrasives and chemical mechanical polishing (CMP) of sapphire wafers, resulting in potential environmental pollution. Developing a novel green CMP technique to achieve light-emitting diode sapphire wafers is a significant challenge. In this study, a novel green CMP slurry, consisting of silica, sorbitol, aminomethyl propanol, and deionized water was developed for sapphire wafers. After CMP, the sapphire wafers were cleaned with deionized water and dried with compressed air, which is a green process. After CMP, the surface roughness Ra of the sapphire wafer surface with an area of 5 × 5 μm2 was 0.098 nm, which is the lowest surface roughness reported to date for sapphire wafers. Tetrahydroxy-coordinated Al(OH)4- ions were produced in the alkaline CMP slurry, and chelation occurred between sorbitol and these ions. The proposed green CMP has potential applications in the semiconductor and microelectronics industries.

Evidence of shock-compressed stishovite above 300 GPa

SiO₂ is one of the most fundamental constituents in planetary bodies, being an essential building block of major mineral phases in the crust and mantle of terrestrial planets (1-10 M_E). Silica at depths greater than 300 km may be present in the form of the rutile-type, high pressure polymorph stishovite (P4₂/mnm) and its thermodynamic stability is of great interest for understanding the seismic and dynamic structure of planetary interiors. Previous studies on stishovite via static and dynamic (shock) compression techniques are contradictory and the observed differences in the lattice-level response is still not clearly understood. Here, laser-induced shock compression experiments at the LCLS- and SACLA XFEL light-sources elucidate the high-pressure behavior of stishovite on the lattice-level under in situ conditions on the Hugoniot to pressures above 300 GPa. We find stishovite is still (meta-)stable at these conditions, and does not undergo any phase transitions. This contradicts static experiments showing structural transformations to the CaCl₂, α-PbO₂ and pyrite-type structures. However, rate-limited kinetic hindrance may explain our observations. These results are important to our understanding into the validity of EOS data from nanosecond experiments for geophysical applications.

Gated detection of supercontinuum pulses enables optical probing of solid and molten silicates at extreme pressure-temperature conditions

Optical studies of materials at high pressure-temperature (P-T) conditions provide insights into their physical properties that may be inaccessible to direct determination at extreme conditions. Incandescent light sources, however, are insufficiently bright to optically probe samples with radiative temperatures above ∼1000 K. Here we report on a system to perform optical absorption experiments in a laser-heated diamond anvil cell at T up to at least 4000 K. This setup is based on a pulsed supercontinuum (broadband) light probe and a gated CCD detector. Precise and tight synchronization of the detector gates (3 ns) to the bright probe pulses (1 ns) diminishes the recorded thermal background and preserves an excellent probe signal at high temperature. We demonstrate the efficiency of this spectroscopic setup by measuring the optical absorbance of solid and molten (Mg,Fe)SiO₃, an important constituent of planetary mantles, at P ∼30 GPa and T ∼1200 K to 4150 K. Optical absorbance of the hot solid (Mg,Fe)SiO₃ is moderately sensitive to temperature but increases abruptly upon melting and acquires a strong temperature dependence. Our results enable quantitative estimates of the opacity of planetary mantles with implications to their thermal and electrical conductivities, all of which have never been constrained at representative P-T conditions, and call for an optical detection of melting in silicate-bearing systems to resolve the extant ambiguity in their high-pressure melting curves.

Laser-driven shock compression of "synthetic planetary mixtures" of water, ethanol, and ammonia

Water, methane, and ammonia are commonly considered to be the key components of the interiors of Uranus and Neptune. Modelling the planets' internal structure, evolution, and dynamo heavily relies on the properties of the complex mixtures with uncertain exact composition in their deep interiors. Therefore, characterising icy mixtures with varying composition at planetary conditions of several hundred gigapascal and a few thousand Kelvin is crucial to improve our understanding of the ice giants. In this work, pure water, a water-ethanol mixture, and a water-ethanol-ammonia "synthetic planetary mixture" (SPM) have been compressed through laser-driven decaying shocks along their principal Hugoniot curves up to 270, 280, and 260 GPa, respectively. Measured temperatures spanned from 4000 to 25000 K, just above the coldest predicted adiabatic Uranus and Neptune profiles (3000-4000 K) but more similar to those predicted by more recent models including a thermal boundary layer (7000-14000 K). The experiments were performed at the GEKKO XII and LULI2000 laser facilities using standard optical diagnostics (Doppler velocimetry and optical pyrometry) to measure the thermodynamic state and the shock-front reflectivity at two different wavelengths. The results show that water and the mixtures undergo a similar compression path under single shock loading in agreement with Density Functional Theory Molecular Dynamics (DFT-MD) calculations using the Linear Mixing Approximation (LMA). On the contrary, their shock-front reflectivities behave differently by what concerns both the onset pressures and the saturation values, with possible impact on planetary dynamos.

Evidence for Fe-Si-O liquid immiscibility at deep Earth pressures

Seismic observations suggest that the uppermost region of Earth's liquid outer core is buoyant, with slower velocities than the bulk outer core. One possible mechanism for the formation of a stably stratified layer is immiscibility in molten iron alloy systems, which has yet to be demonstrated at core pressures. We find immiscibility between liquid Fe-Si and Fe-Si-O persisting to at least 140 GPa through a combination of laser-heated diamond-anvil cell experiments and first-principles molecular dynamics simulations. High-pressure immiscibility in the Fe-Si-O system may explain a stratified layer atop the outer core, complicate differentiation and evolution of the deep Earth, and affect the structure and intensity of Earth's magnetic field. Our results support silicon and oxygen as coexisting light elements in the core and suggest that [Formula: see text] does not crystallize out of molten Fe-Si-O at the core-mantle boundary.

An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA

An experimental platform using X-ray free-electron laser (XFEL) pulses with high-intensity optical laser pulses is open for early users' experiments at the SACLA XFEL facility after completion of the commissioning. The combination of the hard XFEL and the high-intensity laser provides capabilities to open new frontiers of laser-based high-energy-density science. During the commissioning phase, characterization of the XFEL and the laser at the platform has been carried out for the combinative utilization as well as the development of instruments and basic diagnostics for user experiments. An overview of the commissioning and the current capabilities of the experimental platform is presented.

Insights into the Hydrothermal Metastability of Stishovite and Coesite

Hydrothermal experiments aiming at the crystal growth of stishovite near ambient pressure and temperature were performed in conventional autoclave systems using 1 M (molar) NaOH, 0.8 M Na₂CO₃, and pure water as a mineralizing agent. It was found that the hydrothermal metastability of stishovite and coesite is very different from the thermal metastability in all mineralizing agents and that because of this fact crystals could not be grown. While stishovite and coesite are thermally metastable up to 500 and >1000 °C, respectively, their hydrothermal metastability is below 150 and 200 °C, respectively. The thermally induced conversion of stishovite and coesite leads to amorphous products, whereas the hydrothermally induced conversion leads to crystalline quartz. Both stishovite and coesite are minerals occurring in nature where they can be exposed to hydrothermal conditions. The low hydrothermal stability of these phases may be an important factor to explain the rarity of these minerals in nature.

Nanosecond Melting and Recrystallization in Shock-Compressed Silicon

In situ, time-resolved, x-ray diffraction and simultaneous continuum measurements were used to examine structural changes in Si shock compressed to 54 GPa. Shock melting was unambiguously established above ∼31-33  GPa, through the vanishing of all sharp crystalline diffraction peaks and the emergence of a single broad diffraction ring. Reshock from the melt boundary results in rapid (nanosecond) recrystallization to the hexagonal-close-packed Si phase and further supports melting. Our results also provide new constraints on the high-temperature, high-pressure Si phase diagram.

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.

Continuous Sound Velocity Measurements along the Shock Hugoniot Curve of Quartz

We report continuous measurements of the sound velocity along the principal Hugoniot curve of α quartz between 0.25 and 1.45 TPa, as determined from lateral release waves intersecting the shock front as a function of time in decaying-shock experiments. The measured sound velocities are lower than predicted by prior models, based on the properties of stishovite at densities below ∼7  g/cm^{3}, but agree with density functional theory molecular dynamics calculations and an empirical wide-regime equation of state presented here. The Grüneisen parameter calculated from the sound velocity decreases from γ∼1.3 at 0.25 TPa to 0.66 at 1.45 TPa. In combination with evidence for increased (configurational) specific heat and decreased bulk modulus, the values of γ suggest a high thermal expansion coefficient at ∼0.25-0.65  TPa, where SiO_{2} is thought to be a bonded liquid. From our measurements, dissociation of the molecular bonds persists to ∼0.65-1.0  TPa, consistent with estimates by other methods. At higher densities, the sound velocity is close to predictions from previous models, and the Grüneisen parameter approaches the ideal gas value.

Local structure of molten AuGa2 under pressure: Evidence for coordination change and planetary implications

In situ x-ray diffraction measurements and inverse Monte Carlo simulations of pair distribution functions were used to characterize the local structure of molten AuGa₂ up to 16 GPa and 940 K. Our results document systematic changes in liquid structure due to a combination of bond compression and coordination increase. Empirical potential structure refinement shows the first-neighbor coordination of Ga around Au and of Au around Ga to increase from about 8 to 10 and 4 to 5, respectively between 0 and 16 GPa, and the inferred changes in liquid structure can explain the observed melting-point depression of AuGa₂ up to 5 GPa. As intermetallic AuGa₂ is an analogue for metallic SiO₂ at much higher pressures, our results imply that structural changes documented for non-metallic silicate melts below 100 GPa are followed by additional coordination changes in the metallic state at pressures in the 0.2–1 TPa range achieved inside large planets.

In situ X-Ray Diffraction of Shock-Compressed Fused Silica

Because of its widespread applications in materials science and geophysics, SiO_{2} has been extensively examined under shock compression. Both quartz and fused silica transform through a so-called "mixed-phase region" to a dense, low compressibility high-pressure phase. For decades, the nature of this phase has been a subject of debate. Proposed structures include crystalline stishovite, another high-pressure crystalline phase, or a dense amorphous phase. Here we use plate-impact experiments and pulsed synchrotron x-ray diffraction to examine the structure of fused silica shock compressed to 63 GPa. In contrast to recent laser-driven compression experiments, we find that fused silica adopts a dense amorphous structure at 34 GPa and below. When compressed above 34 GPa, fused silica transforms to untextured polycrystalline stishovite. Our results can explain previously ambiguous features of the shock-compression behavior of fused silica and are consistent with recent molecular dynamics simulations. Stishovite grain sizes are estimated to be ∼5-30  nm for compression over a few hundred nanosecond time scale.

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.

Characterizations of dynamic material properties on compact pulsed power generator CQ-4

Over last two decades, the techniques of magnetically driven quasi-isentropic compression and launching high velocity flyer plates based on pulsed high current generators have being extensively used to do dynamic material experiments under extreme conditions, such as high pressure, high temperature and high strain rate. A compact pulsed power generator CQ-4 was developed to do quasi-isentropic compression experiments of materials at Institute of Fluid Physics of CAEP, which can deliver maximum peak current of about 4 MA to short-circuit loads and produce approximate 100 GPa pressure on the metallic samples. On CQ-4, several types of dynamic material experiments have being conducted for equation of states, phase transitions, constitutive relationships, micro-structure evolutions of matter under quasi-isentropic compression and shock loadings. Meanwhile the dynamic behaviors of solid plastic bonded explosives and their components have also being researched for better understanding the interaction of explosive components under stress waves and the hot spot originations and evolutions mechanism of PBX explosives under dynamic loadings. Several typical applications in dynamic material properties were shown in this paper to exhibit the capabilities of CQ-4.

SiO_{2} Glass Density to Lower-Mantle Pressures

The convection or settling of matter in the deep Earth's interior is mostly constrained by density variations between the different reservoirs. Knowledge of the density contrast between solid and molten silicates is thus of prime importance to understand and model the dynamic behavior of the past and present Earth. SiO_{2} is the main constituent of Earth's mantle and is the reference model system for the behavior of silicate melts at high pressure. Here, we apply our recently developed x-ray absorption technique to the density of SiO_{2} glass up to 110 GPa, doubling the pressure range for such measurements. Our density data validate recent molecular dynamics simulations and are in good agreement with previous experimental studies conducted at lower pressure. Silica glass rapidly densifies up to 40 GPa, but the density trend then flattens to become asymptotic to the density of SiO_{2} minerals above 60 GPa. The density data present two discontinuities at ∼17 and ∼60  GPa that can be related to a silicon coordination increase from 4 to a mixed 5/6 coordination and from 5/6 to sixfold, respectively. SiO_{2} glass becomes denser than MgSiO_{3} glass at ∼40  GPa, and its density becomes identical to that of MgSiO_{3} glass above 80 GPa. Our results on SiO_{2} glass may suggest that a variation of SiO_{2} content in a basaltic or pyrolitic melt with pressure has at most a minor effect on the final melt density, and iron partitioning between the melts and residual solids is the predominant factor that controls melt buoyancy in the lowermost mantle.

Conductivity and dissociation in liquid metallic hydrogen and implications for planetary interiors

Liquid metallic hydrogen (LMH) is the most abundant form of condensed matter in our solar planetary structure. The electronic and thermal transport properties of this metallic fluid are of fundamental interest to understanding hydrogen's mechanism of conduction, atomic or pairing structure, as well as the key input for the magnetic dynamo action and thermal models of gas giants. Here, we report spectrally resolved measurements of the optical reflectance of LMH in the pressure region of 1.4-1.7 Mbar. We analyze the data, as well as previously reported measurements, using the free-electron model. Fitting the energy dependence of the reflectance data yields a dissociation fraction of 65 ± 15%, supporting theoretical models that LMH is an atomic metallic liquid. We determine the optical conductivity of LMH and find metallic hydrogen's static electrical conductivity to be 11,000-15,000 S/cm, substantially higher than the only earlier reported experimental values. The higher electrical conductivity implies that the Jovian and Saturnian dynamo are likely to operate out to shallower depths than previously assumed, while the inferred thermal conductivity should provide a crucial experimental constraint to heat transport models.

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.

Miniature diamond anvils for X-ray Raman scattering spectroscopy experiments at high pressure

X-ray Raman scattering (XRS) spectroscopy is an inelastic scattering method that uses hard X-rays of the order of 10 keV to measure energy-loss spectra at absorption edges of light elements (Si, Mg, O etc.), with an energy resolution below 1 eV. The high-energy X-rays employed with this technique can penetrate thick or dense sample containers such as the diamond anvils employed in high-pressure cells. Here, we describe the use of custom-made conical miniature diamond anvils of less than 500 µm thickness which allow pressure generation of up to 70 GPa. This set-up overcomes the limitations of the XRS technique in very high-pressure measurements (>10 GPa) by drastically improving the signal-to-noise ratio. The conical shape of the base of the diamonds gives a 70° opening angle, enabling measurements in both low- and high-angle scattering geometry. This reduction of the diamond thickness to one-third of the classical diamond anvils considerably lowers the attenuation of the incoming and the scattered beams and thus enhances the signal-to-noise ratio significantly. A further improvement of the signal-to-background ratio is obtained by a recess of ∼20 µm that is milled in the culet of the miniature anvils. This recess increases the sample scattering volume by a factor of three at a pressure of 60 GPa. Examples of X-ray Raman spectra collected at the O K-edge and Si L-edge in SiO₂ glass at high pressures up to 47 GPa demonstrate the significant improvement and potential for spectroscopic studies of low-Z elements at high pressure.

Absolute calibration of the OMEGA streaked optical pyrometer for temperature measurements of compressed materials

Experiments in high-energy-density physics often use optical pyrometry to determine temperatures of dynamically compressed materials. In combination with simultaneous shock-velocity and optical-reflectivity measurements using velocity interferometry, these experiments provide accurate equation-of-state data at extreme pressures (P > 1 Mbar) and temperatures (T > 0.5 eV). This paper reports on the absolute calibration of the streaked optical pyrometer (SOP) at the Omega Laser Facility. The wavelength-dependent system response was determined by measuring the optical emission from a National Institute of Standards and Technology-traceable tungsten-filament lamp through various narrowband (40-nm-wide) filters. The integrated signal over the SOP's ∼250-nm operating range is then related to that of a blackbody radiator using the calibrated response. We present a simple closed-form equation for the brightness temperature as a function of streak-camera signal derived from this calibration. Error estimates indicate that brightness temperature can be inferred to a precision of <5%.

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.

Shock compression response of forsterite above 250 GPa

Forsterite (Mg2SiO4) is one of the major planetary materials, and its behavior under extreme conditions is important to understand the interior structure of large planets, such as super-Earths, and large-scale planetary impact events. Previous shock compression measurements of forsterite indicate that it may melt below 200 GPa, but these measurements did not go beyond 200 GPa. We report the shock response of forsterite above ~250 GPa, obtained using the laser shock wave technique. We simultaneously measured the Hugoniot and temperature of shocked forsterite and interpreted the results to suggest the following: (i) incongruent crystallization of MgO at 271 to 285 GPa, (ii) phase transition of MgO at 285 to 344 GPa, and (iii) remelting above ~470 to 500 GPa. These exothermic and endothermic reactions are seen to occur under extreme conditions of pressure and temperature. They indicate complex structural and chemical changes in the system MgO-SiO2 at extreme pressures and temperatures and will affect the way we understand the interior processes of large rocky planets as well as material transformation by impacts in the formation of planetary systems.

Terapascal static pressure generation with ultrahigh yield strength nanodiamond

Studies of materials' properties at high and ultrahigh pressures lead to discoveries of unique physical and chemical phenomena and a deeper understanding of matter. In high-pressure research, an achievable static pressure limit is imposed by the strength of available strong materials and design of high-pressure devices. Using a high-pressure and high-temperature technique, we synthesized optically transparent microballs of bulk nanocrystalline diamond, which were found to have an exceptional yield strength (~460 GPa at a confining pressure of ~70 GPa) due to the unique microstructure of bulk nanocrystalline diamond. We used the nanodiamond balls in a double-stage diamond anvil cell high-pressure device that allowed us to generate static pressures beyond 1 TPa, as demonstrated by synchrotron x-ray diffraction. Outstanding mechanical properties (strain-dependent elasticity, very high hardness, and unprecedented yield strength) make the nanodiamond balls a unique device for ultrahigh static pressure generation. Structurally isotropic, homogeneous, and made of a low-Z material, they are promising in the field of x-ray optical applications.

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.

What We Talk About When We Talk About Light

UNESCO (the United Nations Educational, Scientific, and Cultural Organization) has declared 2015 the "International Year of Light and Light-Based Technologies". In celebration of this proclamation, this Outlook provides a general history of light and its applications, from the earliest moments of the Big Bang through its present impact on all forms of life on the planet. Special emphasis is placed on fundamental advances in the generation and use of artificial light, as well as the harvesting and use of light from the Sun and other natural sources. During the past century, the role of light in the fields of physics, chemistry, and biology has expanded to include emerging fields such as environmental engineering, agriculture, materials science, and biomedicine. In this regard, future research challenges and new potential applications in these areas, in the context of "the central science", are presented and discussed.