Real-time imaging of acoustic waves in bulk materials with X-ray microscopy.

The dynamics of lattice vibrations govern many material processes, such as acoustic wave propagation, displacive phase transitions, and ballistic thermal transport. The maximum velocity of these processes and their effects is determined by the speed of sound, which therefore defines the temporal resolution (picoseconds) needed to resolve these phenomena on their characteristic length scales (nanometers). Here, we present an X-ray microscope capable of imaging acoustic waves with subpicosecond resolution within mm-sized crystals. We directly visualize the generation, propagation, branching, and energy dissipation of longitudinal and transverse acoustic waves in diamond, demonstrating how mechanical energy thermalizes from picosecond to microsecond timescales. Bulk characterization techniques capable of resolving this level of structural detail have previously been available on millisecond time scales-orders of magnitude too slow to capture these fundamental phenomena in solid-state physics and geoscience. As such, the reported results provide broad insights into the interaction of acoustic waves with the structure of materials, and the availability of ultrafast time-resolved dark-field X-ray microscopy opens a vista of new opportunities for 3D imaging of materials dynamics on their intrinsic submicrosecond time scales.

Nanosecond X-ray diffraction of shock-compressed superionic water ice.

Since Bridgman's discovery of five solid water (H2O) ice phases1 in 1912, studies on the extraordinary polymorphism of H2O have documented more than seventeen crystalline and several amorphous ice structures2,3, as well as rich metastability and kinetic effects4,5. This unique behaviour is due in part to the geometrical frustration of the weak intermolecular hydrogen bonds and the sizeable quantum motion of the light hydrogen ions (protons). Particularly intriguing is the prediction that H2O becomes superionic6-12-with liquid-like protons diffusing through the solid lattice of oxygen-when subjected to extreme pressures exceeding 100 gigapascals and high temperatures above 2,000 kelvin. Numerical simulations suggest that the characteristic diffusion of the protons through the empty sites of the oxygen solid lattice (1) gives rise to a surprisingly high ionic conductivity above 100 Siemens per centimetre, that is, almost as high as typical metallic (electronic) conductivity, (2) greatly increases the ice melting temperature7-13 to several thousand kelvin, and (3) favours new ice structures with a close-packed oxygen lattice13-15. Because confining such hot and dense H2O in the laboratory is extremely challenging, experimental data are scarce. Recent optical measurements along the Hugoniot curve (locus of shock states) of water ice VII showed evidence of superionic conduction and thermodynamic signatures for melting16, but did not confirm the microscopic structure of superionic ice. Here we use laser-driven shockwaves to simultaneously compress and heat liquid water samples to 100-400 gigapascals and 2,000-3,000 kelvin. In situ X-ray diffraction measurements show that under these conditions, water solidifies within a few nanoseconds into nanometre-sized ice grains that exhibit unambiguous evidence for the crystalline oxygen lattice of superionic water ice. The X-ray diffraction data also allow us to document the compressibility of ice at these extreme conditions and a temperature- and pressure-induced phase transformation from a body-centred-cubic ice phase (probably ice X) to a novel face-centred-cubic, superionic ice phase, which we name ice XVIII2,17.

A MHz X-ray diffraction set-up for dynamic compression experiments in the diamond anvil cell.

An experimental platform for dynamic diamond anvil cell (dDAC) research has been developed at the High Energy Density (HED) Instrument at the European X-ray Free Electron Laser (European XFEL). Advantage was taken of the high repetition rate of the European XFEL (up to 4.5 MHz) to collect pulse-resolved MHz X-ray diffraction data from samples as they are dynamically compressed at intermediate strain rates (≤103 s-1), where up to 352 diffraction images can be collected from a single pulse train. The set-up employs piezo-driven dDACs capable of compressing samples in ≥340 µs, compatible with the maximum length of the pulse train (550 µs). Results from rapid compression experiments on a wide range of sample systems with different X-ray scattering powers are presented. A maximum compression rate of 87 TPa s-1 was observed during the fast compression of Au, while a strain rate of ∼1100 s-1 was achieved during the rapid compression of N2 at 23 TPa s-1.

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.

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.

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon.

Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

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.

Simultaneous bright- and dark-field X-ray microscopy at X-ray free electron lasers.

The structures, strain fields, and defect distributions in solid materials underlie the mechanical and physical properties across numerous applications. Many modern microstructural microscopy tools characterize crystal grains, domains and defects required to map lattice distortions or deformation, but are limited to studies of the (near) surface. Generally speaking, such tools cannot probe the structural dynamics in a way that is representative of bulk behavior. Synchrotron X-ray diffraction based imaging has long mapped the deeply embedded structural elements, and with enhanced resolution, dark field X-ray microscopy (DFXM) can now map those features with the requisite nm-resolution. However, these techniques still suffer from the required integration times due to limitations from the source and optics. This work extends DFXM to X-ray free electron lasers, showing how the [Formula: see text] photons per pulse available at these sources offer structural characterization down to 100 fs resolution (orders of magnitude faster than current synchrotron images). We introduce the XFEL DFXM setup with simultaneous bright field microscopy to probe density changes within the same volume. This work presents a comprehensive guide to the multi-modal ultrafast high-resolution X-ray microscope that we constructed and tested at two XFELs, and shows initial data demonstrating two timing strategies to study associated reversible or irreversible lattice dynamics.

Structural complexity in ramp-compressed sodium to 480 GPa.

The properties of all materials at one atmosphere of pressure are controlled by the configurations of their valence electrons. At extreme pressures, neighboring atoms approach so close that core-electron orbitals overlap, and theory predicts the emergence of unusual quantum behavior. We ramp-compress monovalent elemental sodium, a prototypical metal at ambient conditions, to nearly 500 GPa (5 million atmospheres). The 7-fold increase of density brings the interatomic distance to 1.74 Å well within the initial 2.03 Å of the Na+ ionic diameter, and squeezes the valence electrons into the interstitial voids suggesting the formation of an electride phase. The laser-driven compression results in pressure-driven melting and recrystallization in a billionth of a second. In situ x-ray diffraction reveals a series of unexpected phase transitions upon recrystallization, and optical reflectivity measurements show a precipitous decrease throughout the liquid and solid phases, where the liquid is predicted to have electronic localization. These data reveal the presence of a rich, temperature-driven polymorphism where core electron overlap is thought to stabilize the formation of peculiar electride states.

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.

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.

Interferometric measurements of refractive index and dispersion at high pressure.

We describe a high precision interferometer system to measure the pressure dependence of the refractive index and its dispersion in the diamond anvil cell (DAC). The reflective Fabry-Perot fringe patterns created by both a white light and a monochromatic beam are recorded to determine both the sample thickness and its index at the laser wavelength and to characterize the dispersion in the visible range. Advances in sample preparation, optical setup, and data analysis enable us to achieve [Formula: see text] random uncertainty, demonstrated with an air sample, a factor of five improvement over the best previous DAC measurement. New data on [Formula: see text] liquid water and ice VI up to 2.21 GPa at room temperature illustrate how higher precision measurements of the index and its optical dispersion open up new opportunities to reveal subtle changes in the electronic structure of water at high pressure.

In situ visualization of long-range defect interactions at the edge of melting.

Connecting a bulk material's microscopic defects to its macroscopic properties is an age-old problem in materials science. Long-range interactions between dislocations (line defects) are known to play a key role in how materials deform or melt, but we lack the tools to connect these dynamics to the macroscopic properties. We introduce time-resolved dark-field x-ray microscopy to directly visualize how dislocations move and interact over hundreds of micrometers deep inside bulk aluminum. With real-time movies, we reveal the thermally activated motion and interactions of dislocations that comprise a boundary and show how weakened binding forces destabilize the structure at 99% of the melting temperature. Connecting dynamics of the microstructure to its stability, we provide important opportunities to guide and validate multiscale models that are yet untested.

X-ray Free Electron Laser-Induced Synthesis of ε-Iron Nitride at High Pressures.

The ultrafast synthesis of ε-Fe3N1+x in a diamond-anvil cell (DAC) from Fe and N2 under pressure was observed using serial exposures of an X-ray free electron laser (XFEL). When the sample at 5 GPa was irradiated by a pulse train separated by 443 ns, the estimated sample temperature at the delay time was above 1400 K, confirmed by in situ transformation of α- to γ-iron. Ultimately, the Fe and N2 reacted uniformly throughout the beam path to form Fe3N1.33, as deduced from its established equation of state (EOS). We thus demonstrate that the activation energy provided by intense X-ray exposures in an XFEL can be coupled with the source time structure to enable exploration of the time-dependence of reactions under high-pressure conditions.

Heat-release equation of state and thermal conductivity of warm dense carbon by proton differential heating.

Warm dense carbon is generated at 0.3-2.0 g/cc and 1-7 eV by proton heating. The release equation of state (EOS) after heating and thermal conductivity of warm dense carbon are studied experimentally in this regime using a Au/C dual-layer target to initiate a temperature gradient and two picosecond time-resolved diagnostics to probe the surface expansion and heat flow. Comparison between the data and simulations using various EOSs and thermal conductivity models is quantified with a statistical χ^{2} analysis. Out of seven EOS tables and five thermal conductivity models, only L9061 with the Lee-More model provides a probability above 50% to match all data.

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.

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.

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.