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.

Laser-driven ramp-compression experiments on the national ignition facility.

This report details the analyses and related uncertainties in measuring longitudinal-stress-density paths in indirect laser-driven ramp equation-of-state (EOS) experiments [Smith et al., Nat. Astron. 2(6), 452-458 (2018); Smith et al., Nature 511(7509), 330-333 (2014); Fratanduono et al., Science 372(6546), 1063-1068 (2021); and Fratanduono et al., Phys. Rev. Lett. 124(1), 015701 (2020)]. Experiments were conducted at the National Ignition Facility (NIF) located at the Lawrence Livermore National Laboratory. The NIF can deliver up to 2 MJ of laser energy over 30 ns and provide the necessary laser power and control to ramp compress materials to TPa pressures (1 TPa = 10 × 106 atmospheres). These data provide low-temperature solid-state EOS data relevant to the extreme conditions found in the deep interiors of giant planets. In these experiments, multi-stepped samples with thicknesses in the range of 40-120 µm experience an initial shock compression followed by a time-dependent ramp compression to peak pressure. Interface velocity measurements from each thickness combine to place a constraint on the Lagrangian sound speed as a function of particle velocity, which in turn allows for the determination of a continuous stress-density path to high levels of compressibility. In this report, we present a detailed description of the experimental techniques and measurement uncertainties and describe how these uncertainties combine to place a final uncertainty in both stress and density. We address the effects of time-dependent deformation and the sensitivity of ramp EOS techniques to the onset of phase transformations.

High pressure phase transition and strength estimate in polycrystalline alumina during laser-driven shock compression.

Alumina (Al2O3) is an important ceramic material notable for its compressive strength and hardness. It represents one of the major oxide components of the Earth's mantle. Static compression experiments have reported evidence for phase transformations from the trigonalα-corundum phase to the orthorhombic Rh2O3(II)-type structure at ∼90 GPa, and then to the post-perovskite structure at ∼130 GPa, but these phases have yet to be directly observed under shock compression. In this work, we describe laser-driven shock compression experiments on polycrystalline alumina conducted at the Matter in Extreme Conditions endstation of the Linac Coherent Light Source. Ultrafast x-ray pulses (50 fs, 1012photons/pulse) were used to probe the atomic-level response at different times during shock propagation and subsequent pressure release. At 107 ± 8 GPa on the Hugoniot, we observe diffraction peaks that match the orthorhombic Rh2O3(II) phase with a density of 5.16 ± 0.03 g cm-3. Upon unloading, the material transforms back to theα-corundum structure. Upon release to ambient pressure, densities are lower than predicted assuming isentropic release, indicating additional lattice expansion due to plastic work heating. Using temperature values calculated from density measurements, we provide an estimate of alumina's strength on release from shock compression.

Quantitative analysis of diffraction by liquids using a pink-spectrum X-ray source.

A new approach for performing quantitative structure-factor analysis and density measurements of liquids using X-ray diffraction with a pink-spectrum X-ray source is described. The methodology corrects for the pink beam effect by performing a Taylor series expansion of the diffraction signal. The mean density, background scale factor, peak X-ray energy about which the expansion is performed, and the cutoff radius for density measurement are estimated using the derivative-free optimization scheme. The formalism is demonstrated for a simulated radial distribution function for tin. Finally, the proposed methodology is applied to experimental data on shock compressed tin recorded at the Dynamic Compression Sector at the Advanced Photon Source, with derived densities comparing favorably with other experimental results and the equations of state of tin.

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa.

Triple bonding in the nitrogen molecule (N_{2}) is among the strongest chemical bonds with a dissociation enthalpy of 9.8 eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0 g/cm^{3}. Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29 g/cm^{3}, our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

Measuring the melting curve of iron at super-Earth core conditions.

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth's interior. We used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth's inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth's mass have the longest dynamos, which provide important shielding against cosmic radiation.

Non-iterative characteristics analysis for high-pressure ramp loading.

In the canonical ramp compression experiment, a smoothly increasing load is applied to the surface of the sample, and the particle velocity history is measured at two or more different distances into the sample, at interfaces where the surface of the sample can be probed. The velocity histories are used to deduce a stress-density relation, usually using iterative Lagrangian analysis to account for the perturbing effect of the impedance mismatch at the interface. In that technique, a stress-density relation is assumed in order to correct for the perturbation and is adjusted until it becomes consistent with the deduced stress-density relation. This process is subject to the usual difficulties of nonlinear optimization, such as the existence of local minima (sensitivity to the initial guess), possible failure to converge, and relatively large computational effort. We show that, by considering the interaction of successive characteristics reaching a free surface, the stress-density relation can be deduced directly by recursion rather than iteration. This calculation is orders of magnitude faster than iterative analysis and does not require an initial guess. Direct recursion may be less suitable for very noisy data, but it was robust when applied to trial data. The stress-density relation deduced was identical to the result from iterative Lagrangian analysis.

Response to Comment on "Insulator-metal transition in dense fluid deuterium".

In their comment, Desjarlais et al claim that a small temperature drop occurs after isentropic compression of fluid deuterium through the first-order insulator-metal transition. We show that their calculations do not correspond to the experimental thermodynamic path, and that thermodynamic integrations with parameters from first-principles calculations produce results in agreement with our original estimate of the temperature drop.

Insulator-metal transition in dense fluid deuterium.

Dense fluid metallic hydrogen occupies the interiors of Jupiter, Saturn, and many extrasolar planets, where pressures reach millions of atmospheres. Planetary structure models must describe accurately the transition from the outer molecular envelopes to the interior metallic regions. We report optical measurements of dynamically compressed fluid deuterium to 600 gigapascals (GPa) that reveal an increasing refractive index, the onset of absorption of visible light near 150 GPa, and a transition to metal-like reflectivity (exceeding 30%) near 200 GPa, all at temperatures below 2000 kelvin. Our measurements and analysis address existing discrepancies between static and dynamic experiments for the insulator-metal transition in dense fluid hydrogen isotopes. They also provide new benchmarks for the theoretical calculations used to construct planetary models.

Crystal structure and equation of state of Fe-Si alloys at super-Earth core conditions.

The high-pressure behavior of Fe alloys governs the interior structure and dynamics of super-Earths, rocky extrasolar planets that could be as much as 10 times more massive than Earth. In experiments reaching up to 1300 GPa, we combine laser-driven dynamic ramp compression with in situ x-ray diffraction to study the effect of composition on the crystal structure and density of Fe-Si alloys, a potential constituent of super-Earth cores. We find that Fe-Si alloy with 7 weight % (wt %) Si adopts the hexagonal close-packed structure over the measured pressure range, whereas Fe-15wt%Si is observed in a body-centered cubic structure. This study represents the first experimental determination of the density and crystal structure of Fe-Si alloys at pressures corresponding to the center of a ~3-Earth mass terrestrial planet. Our results allow for direct determination of the effects of light elements on core radius, density, and pressures for these planets.