Laser-shock compression of diamond and evidence of a negative-slope melting curve

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.

Double-Shock Compression Pathways from Diamond to BC8 Carbon

Carbon is one of the most important elements for both industrial applications and fundamental research, including life, physics, chemistry, materials, and even planetary science. Although theoretical predictions on the transition from diamond to the BC8 (Ia3[over ¯]) carbon were made more than thirty years ago, after tremendous experimental efforts, direct evidence for the existence of BC8 carbon is still lacking. In this study, a machine learning potential was developed for high-pressure carbon fitted from first-principles calculations, which exhibited great capabilities in modeling the melting and Hugoniot line. Using the molecular dynamics based on this machine learning potential, we designed a thermodynamic pathway that is achievable for the double shock compression experiment to obtain the elusive BC8 carbon. Diamond was compressed up to 584 GPa after the first shock at 20.5  km/s. Subsequently, in the second shock compression at 24.8 or 25.0  km/s, diamond was compressed to a supercooled liquid and then solidified to BC8 in around 1 ns. Furthermore, the critical nucleus size and nucleation rate of BC8 were calculated, which are crucial for nano-second x-ray diffraction measurements to observe BC8 carbon during shock compressions. The key to obtaining BC8 carbon lies in the formation of liquid at a sufficient supercooling. Our work provides a feasible pathway by which the long-sought BC8 phase of carbon can be reached in experiments.

Diagnosis of ultrafast surface dynamics of thin foil targets irradiated by intense laser pulses

The temporal modulation of an electron bunch train accelerated from a foil target irradiated by an intense laser pulse is studied by measuring the coherent transition radiation (CTR) from the rear surface of a target. We experimentally obtained CTR spectra from a 1 µm thick foil target irradiated at a maximum intensity of 6.5 × 10¹⁹ W/cm². Spectral redshifts of the emitted radiation corresponding to increases in laser intensity were observed. These measurements were compared with the theoretical calculation of CTR spectra considering ultrafast surface dynamics, such as plasma surface oscillation and relativistically induced transparency. Plasma surface oscillations induce a spectral redshift, while relativistic transparency causes a spectral blueshift. Both effects are required to find reasonable agreement with the experiment over the entire range of laser intensities.

Probing extreme states of matter using ultra-intense x-ray radiation

Extreme states of matter, that is, matter at extremes of density (pressure) and temperature, can be created in the laboratory either statically or dynamically. In the former, the pressure-temperature state can be maintained for relatively long periods of time, but the sample volume is necessarily extremely small. When the extreme states are generated dynamically, the sample volumes can be larger, but the pressure-temperature conditions are maintained for only short periods of time (ps toμs). In either case, structural information can be obtained from the extreme states by the use of x-ray scattering techniques, but the x-ray beam must be extremely intense in order to obtain sufficient signal from the extremely-small or short-lived sample. In this article I describe the use of x-ray diffraction at synchrotrons and XFELs to investigate how crystal structures evolve as a function of density and temperature. After a brief historical introduction, I describe the developments made at the Synchrotron Radiation Source in the 1990s which enabled the almost routine determination of crystal structure at high pressures, while also revealing that the structural behaviour of materials was much more complex than previously believed. I will then describe how these techniques are used at the current generation of synchrotron and XFEL sources, and then discuss how they might develop further in the future at the next generation of x-ray lightsources.

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ∼0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

Direct observation of imploded core heating via fast electrons with super-penetration scheme

Fast ignition (FI) is a promising approach for high-energy-gain inertial confinement fusion in the laboratory. To achieve ignition, the energy of a short-pulse laser is required to be delivered efficiently to the pre-compressed fuel core via a high-energy electron beam. Therefore, understanding the transport and energy deposition of this electron beam inside the pre-compressed core is the key for FI. Here we report on the direct observation of the electron beam transport and deposition in a compressed core through the stimulated Cu Kα emission in the super-penetration scheme. Simulations reproducing the experimental measurements indicate that, at the time of peak compression, about 1% of the short-pulse energy is coupled to a relatively low-density core with a radius of 70 μm. Analysis with the support of 2D particle-in-cell simulations uncovers the key factors improving this coupling efficiency. Our findings are of critical importance for optimizing FI experiments in a super-penetration scheme.

Guiding of relativistic electron beams in dense matter by laser-driven magnetostatic fields

Intense lasers interacting with dense targets accelerate relativistic electron beams, which transport part of the laser energy into the target depth. However, the overall laser-to-target energy coupling efficiency is impaired by the large divergence of the electron beam, intrinsic to the laser-plasma interaction. Here we demonstrate that an efficient guiding of MeV electrons with about 30 MA current in solid matter is obtained by imposing a laser-driven longitudinal magnetostatic field of 600 T. In the magnetized conditions the transported energy density and the peak background electron temperature at the 60-μm-thick target's rear surface rise by about a factor of five, as unfolded from benchmarked simulations. Such an improvement of energy-density flux through dense matter paves the ground for advances in laser-driven intense sources of energetic particles and radiation, driving matter to extreme temperatures, reaching states relevant for planetary or stellar science as yet inaccessible at the laboratory scale and achieving high-gain laser-driven thermonuclear fusion.

Free energy model for solid high-pressure phases of carbon

Analytic free energy models for three solid high-pressure phases--diamond, body centered cubic phase with eight atoms in the unit cell (BC8), and simple cubic (SC)--are developed using density functional theory. We explicitly include anharmonic effects by performing molecular dynamics simulations and investigate their density and temperature dependence in detail. Anharmonicity in the nuclear motion shifts the phase transitions significantly compared to the harmonic approximation. Furthermore, we apply a thermodynamically constrained correction that brings the equation of state in accordance with diamond anvil cell experiments. The performance of our thermodynamic functions is validated against Hugoniot experiments.

X-Ray Diffraction of Solid Tin to 1.2 TPa

We report direct in situ measurements of the crystal structure of tin between 0.12 and 1.2 TPa, the highest stress at which a crystal structure has ever been observed. Using angle-dispersive powder x-ray diffraction, we find that dynamically compressed Sn transforms to the body-centered-cubic (bcc) structure previously identified by ambient-temperature quasistatic-compression studies and by zero-kelvin density-functional theory predictions between 0.06 and 0.16 TPa. However, we observe no evidence for the hexagonal close-packed (hcp) phase found by those studies to be stable above 0.16 TPa. Instead, our results are consistent with bcc up to 1.2 TPa. We conjecture that at high temperature bcc is stabilized relative to hcp due to differences in vibrational free energy.

Experimental evidence of new tetragonal polymorphs of silicon formed through ultrafast laser-induced confined microexplosion

Ordinary materials can transform into novel phases at extraordinary high pressure and temperature. The recently developed method of ultrashort laser-induced confined microexplosions initiates a non-equilibrium disordered plasma state. Ultra-high quenching rates overcome kinetic barriers to the formation of new metastable phases, which are preserved in the surrounding pristine crystal for subsequent exploitation. Here we demonstrate that confined microexplosions in silicon produce several metastable end phases. Comparison with an ab initio random structure search reveals six energetically competitive potential phases, four tetragonal and two monoclinic structures. We show the presence of bt8 and st12, which have been predicted theoretically previously, but have not been observed in nature or in laboratory experiments. In addition, the presence of the as yet unidentified silicon phase, Si-VIII and two of our other predicted tetragonal phases are highly likely within laser-affected zones. These findings may pave the way for new materials with novel and exotic properties.

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.

Ramp compression of diamond to five terapascals

The recent discovery of more than a thousand planets outside our Solar System, together with the significant push to achieve inertially confined fusion in the laboratory, has prompted a renewed interest in how dense matter behaves at millions to billions of atmospheres of pressure. The theoretical description of such electron-degenerate matter has matured since the early quantum statistical model of Thomas and Fermi, and now suggests that new complexities can emerge at pressures where core electrons (not only valence electrons) influence the structure and bonding of matter. Recent developments in shock-free dynamic (ramp) compression now allow laboratory access to this dense matter regime. Here we describe ramp-compression measurements for diamond, achieving 3.7-fold compression at a peak pressure of 5 terapascals (equivalent to 50 million atmospheres). These equation-of-state data can now be compared to first-principles density functional calculations and theories long used to describe matter present in the interiors of giant planets, in stars, and in inertial-confinement fusion experiments. Our data also provide new constraints on mass-radius relationships for carbon-rich planets.

Ramp compression of diamond to five terapascals

The recent discovery of more than a thousand planets outside our Solar System, together with the significant push to achieve inertially confined fusion in the laboratory, has prompted a renewed interest in how dense matter behaves at millions to billions of atmospheres of pressure. The theoretical description of such electron-degenerate matter has matured since the early quantum statistical model of Thomas and Fermi, and now suggests that new complexities can emerge at pressures where core electrons (not only valence electrons) influence the structure and bonding of matter. Recent developments in shock-free dynamic (ramp) compression now allow laboratory access to this dense matter regime. Here we describe ramp-compression measurements for diamond, achieving 3.7-fold compression at a peak pressure of 5 terapascals (equivalent to 50 million atmospheres). These equation-of-state data can now be compared to first-principles density functional calculations and theories long used to describe matter present in the interiors of giant planets, in stars, and in inertial-confinement fusion experiments. Our data also provide new constraints on mass-radius relationships for carbon-rich planets.

Anomalous sound propagation and slow kinetics in dynamically compressed amorphous carbon

We have performed molecular-dynamics simulations of dynamic compression waves propagating through amorphous carbon using the Tersoff potential and find that a variety of dynamic compression features appear for two different initial densities. These features include steady elastic shocks, steady chemically reactive shocks, unsteady elastic waves, and unsteady chemically reactive waves. We show how these features can be distinguished by analyzing time-dependent propagation speeds, time-dependent sound speeds, and comparison to multiscale shock technique (MSST) simulations. Understanding such features is a key challenge in quasi-isentropic experiments involving phase transformations. In addition to direct simulations of dynamic compression, we employ the MSST and find agreement with the direct method for this system for the shocks observed. We show how the MSST can be extended to include explicit material viscosity and demonstrate on an amorphous Lennard-Jones system.

Diamond at 800 GPa

A new compression technique, which enables the study of solids into the TPa regime, is described and used to ramp (or quasi-isentropically) compress diamond to a peak pressure of 1400 GPa. Diamond stress versus density data are reported to 800 GPa and suggest that the diamond phase is stable and has significant material strength up to at least this stress level. Data presented here are the highest ramp compression pressures by more than a factor of 5 and the highest-pressure solid equation-of-state data ever reported.

Ab initio simulations of the K-edge shift along the aluminum Hugoniot

We develop a first-principles approach to calculate the near-edge absorption spectrum of dense plasmas based on density functional electronic structure calculations and molecular dynamics simulations. We apply the method to the calculation of the K-edge shift along the aluminum shock compressed Hugoniot. We obtain a good agreement with measurements performed at moderate compression and find that the variation of the XANES spectra could be used as a signature for melting along the Hugoniot. We also show that the calculation of the K-edge shift along the Hugoniot formally requires a fully self-consistent calculation beyond the frozen-core approximation and provides an opportunity to test the accuracy of first principle simulation methods in the high-pressure high-temperature regime.