Release dynamics of nanodiamonds created by laser-driven shock-compression of polyethylene terephthalate.

Laser-driven dynamic compression experiments of plastic materials have found surprisingly fast formation of nanodiamonds (ND) via X-ray probing. This mechanism is relevant for planetary models, but could also open efficient synthesis routes for tailored NDs. We investigate the release mechanics of compressed NDs by molecular dynamics simulation of the isotropic expansion of finite size diamond from different P-T states. Analysing the structural integrity along different release paths via molecular dynamic simulations, we found substantial disintegration rates upon shock release, increasing with the on-Hugnoiot shock temperature. We also find that recrystallization can occur after the expansion and hence during the release, depending on subsequent cooling mechanisms. Our study suggests higher ND recovery rates from off-Hugoniot states, e.g., via double-shocks, due to faster cooling. Laser-driven shock compression experiments of polyethylene terephthalate (PET) samples with in situ X-ray probing at the simulated conditions found diamond signal that persists up to 11 ns after breakout. In the diffraction pattern, we observed peak shifts, which we attribute to thermal expansion of the NDs and thus a total release of pressure, which indicates the stability of the released NDs.

Mixture of hydrogen and methane under planetary interior conditions.

We employ first-principles molecular dynamics simulations to provide equation-of-state data, pair distribution functions (PDFs), diffusion coefficients, and band gaps of a mixture of hydrogen and methane under planetary interior conditions as relevant for Uranus, Neptune, and similar icy exoplanets. We test the linear mixing approximation, which is fulfilled within a few percent for the chosen P-T conditions. Evaluation of the PDFs reveals that methane molecules dissociate into carbon clusters and free hydrogen atoms at temperatures greater than 3000 K. At high temperatures, the clusters are found to be short-lived. Furthermore, we calculate the electrical conductivity from which we derive the non-metal-to-metal transition region of the mixture. We also calculate the electrical conductivity along the P-T profile of Uranus [N. Nettelmann et al., Planet. Space Sci., 2013, 77, 143-151] and observe the transition of the mixture from a molecular to an atomic fluid as a function of the radius of the planet. The density and temperature ranges chosen in our study can be achieved using dynamic shock compression experiments and seek to aid such future experiments. Our work also provides a relevant data set for a better understanding of the interior, evolution, luminosity, and magnetic field of the ice giants in our solar system and beyond.

Evidence for phonon hardening in laser-excited gold using x-ray diffraction at a hard x-ray free electron laser.

Studies of laser-heated materials on femtosecond timescales have shown that the interatomic potential can be perturbed at sufficiently high laser intensities. For gold, it has been postulated to undergo a strong stiffening leading to an increase of the phonon energies, known as phonon hardening. Despite efforts to investigate this behavior, only measurements at low absorbed energy density have been performed, for which the interpretation of the experimental data remains ambiguous. By using in situ single-shot x-ray diffraction at a hard x-ray free-electron laser, the evolution of diffraction line intensities of laser-excited Au to a higher energy density provides evidence for phonon hardening.

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula.

We study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

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.

Diamond formation kinetics in shock-compressed C─H─O samples recorded by small-angle x-ray scattering and x-ray diffraction.

Extreme conditions inside ice giants such as Uranus and Neptune can result in peculiar chemistry and structural transitions, e.g., the precipitation of diamonds or superionic water, as so far experimentally observed only for pure C─H and H2O systems, respectively. Here, we investigate a stoichiometric mixture of C and H2O by shock-compressing polyethylene terephthalate (PET) plastics and performing in situ x-ray probing. We observe diamond formation at pressures between 72 ± 7 and 125 ± 13 GPa at temperatures ranging from ~3500 to ~6000 K. Combining x-ray diffraction and small-angle x-ray scattering, we access the kinetics of this exotic reaction. The observed demixing of C and H2O suggests that diamond precipitation inside the ice giants is enhanced by oxygen, which can lead to isolated water and thus the formation of superionic structures relevant to the planets' magnetic fields. Moreover, our measurements indicate a way of producing nanodiamonds by simple laser-driven shock compression of cheap PET plastics.

Electronic transport coefficients from density functional theory across the plasma plane.

We investigate the thermopower and Lorenz number of hydrogen with Kohn-Sham density functional theory (DFT) across the plasma plane toward the near-classical limit, i.e., weakly degenerate and weakly coupled states. Our results are in concordance with certain limiting values for the Lorentz plasma, a model system which only considers electron-ion scattering. Thereby, we clearly show that the widely used method of calculating transport properties via the Kubo-Greenwood (KG) formalism does not capture electron-electron scattering processes. Our discussion also addresses the inadequateness of assuming a Drude-like frequency behavior for the conductivity of nondegenerate plasmas by revisiting the relaxation time approximation within kinetic theory.

Ionization and transport in partially ionized multicomponent plasmas: Application to atmospheres of hot Jupiters.

We study ionization and transport processes in partially ionized multicomponent plasmas. The plasma composition is calculated via a system of coupled mass-action laws. The electronic transport properties are determined by the electron-ion and electron-neutral transport cross sections. The influence of electron-electron scattering is considered via a correction factor to the electron-ion contribution. Based on these data, the electrical and thermal conductivities as well as the Lorenz number are calculated. For the thermal conductivity, we consider also the contributions of the translational motion of neutral particles and of the dissociation, ionization, and recombination reactions. We apply our approach to a partially ionized plasma composed of hydrogen, helium, and a small fraction of metals (Li, Na, Ca, Fe, K, Rb, and Cs) as typical for atmospheres of hot Jupiters. We present results for the plasma composition and the transport properties as a function of density and temperature and then along typical P-T profiles for the outer part of the hot Jupiter HD 209458b. The electrical conductivity profile allows revising the Ohmic heating power related to the fierce winds in the planet's atmosphere. We show that the higher temperatures suggested by recent interior models could boost the conductivity and thus the Ohmic heating power to values large enough to explain the observed inflation of HD 209458b.

Gibbs-ensemble Monte Carlo simulation of H2-H2O mixtures.

The miscibility gap in hydrogen-water mixtures is investigated by conducting Gibbs-ensemble Monte Carlo simulations with analytical two-body interaction potentials between the molecular species. We calculate several demixing curves at pressures below 150 kbar and temperatures of 1000 K ≤T≤ 2000 K. Despite the approximations introduced by the two-body interaction potentials, our results predict a large miscibility gap in hydrogen-water mixtures at similar conditions as found in experiments. Our findings are in contrast to those from ab initio simulations and provide a renewed indication that hydrogen-water immiscibility regions may have a significant impact on the structure and evolution of ice giant planets like Uranus and Neptune.

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.

Gibbs-ensemble Monte Carlo simulation of H_{2}-He mixtures.

We explore the performance of the Gibbs-ensemble Monte Carlo simulation technique by calculating the miscibility gap of H_{2}-He mixtures with analytical exponential-six potentials. We calculate several demixing curves for pressures up to 500 kbar and for temperatures up to 1800K and predict a H_{2}-He miscibility diagram for the solar He abundance for temperatures up to 1500K and determine the demixing region. Our results are in good agreement with ab initio simulations in the nondissociated region of the phase diagram. However, the particle number necessary to converge the Gibbs-ensemble Monte Carlo method is yet too large to offer a feasible combination with ab initio electronic structure calculation techniques, which would be necessary at conditions where dissociation or ionization occurs.

Evidence of shock-compressed stishovite above 300 GPa.

SiO2 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 ME). Silica at depths greater than 300 km may be present in the form of the rutile-type, high pressure polymorph stishovite (P42/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 CaCl2, α-PbO2 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.

Paramagnetic-to-Diamagnetic Transition in Dense Liquid Iron and Its Influence on Electronic Transport Properties.

The electrical σ and thermal conductivity λ of liquid iron are calculated with spin-polarized density-functional-theory-based simulations over a significant pressure and temperature range using the Kubo-Greenwood formalism. We show that a paramagnetic state is stable in the liquid up to high temperatures at ambient pressure and that the discrepancy between experimental results and spin-degenerate simulations for σ and λ of more than 30% are reduced to within 10% with lower values resulting from the spin-polarized simulations. Along the 3700 K isotherm, we explore the persistence of magnetic fluctuations toward high densities, and beyond 20-50 GPa the liquid becomes diamagnetic, which suggests the existence of a continuous paramagnetic-to-diamagnetic transition. This transition exerts a significant influence on the physical properties of liquid iron, especially on σ and λ, and is potentially of high relevance for dynamo processes in Mercury and Mars.

Ionization potential depression and Pauli blocking in degenerate plasmas at extreme densities.

New facilities explore warm dense matter (WDM) at conditions with extreme densities (exceeding ten times condensed matter densities) so that electrons are degenerate even at temperatures of 10-100 eV. Whereas in the nondegenerate region correlation effects such as Debye screening are relevant for the ionization potential depression (IPD), new effects have to be considered in degenerate plasmas. In addition to the Fock shift of the self-energies, the bound-state Pauli blocking becomes important with increasing density. Standard approaches to IPD such as Stewart-Pyatt and widely used opacity tables (e.g., OPAL) do not contain Pauli blocking effects for bound states. The consideration of degeneracy effects leads to a reduction of the ionization potential and to a higher degree of ionization. As an example, we present calculations for the ionization degree of carbon plasmas at T = 100 eV and extreme densities up to 40 g/cm^{3}, which are relevant to experiments that are currently scheduled at the National Ignition Facility.

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

Celliers et al (Reports, 17 August 2018, p. 677), in an attempt to reconcile differences in inferred metallization pressures, provide an alternative temperature analysis of the Knudson et al experiments (Reports, 26 June 2015, p. 1455). We show that this reanalysis implies an anomalously low specific heat for the metallic fluid that is clearly inconsistent with first-principles calculations.

Ab Initio Calculation of the Miscibility Diagram for Hydrogen-Helium Mixtures.

We use finite-temperature density functional theory coupled to classical molecular dynamics simulation to calculate the miscibility gap of hydrogen-helium mixtures. The van der Waals density functional (vdW-DF) theory is used, which leads to lower demixing temperatures compared to computations using the Perdew-Burke-Ernzerhof functional. Our calculations suggest that current Jupiter models are most likely too hot to allow demixing in the interior. A Jupiter isentrope based on our vdW-DF data is presented. Our demixing phase diagram still predicts phase separation in Saturn, but in a significantly reduced fraction of its volume.

Density-functional calculations of transport properties in the nondegenerate limit and the role of electron-electron scattering.

We compute electrical and thermal conductivities of hydrogen plasmas in the nondegenerate regime using Kohn-Sham density functional theory (DFT) and an application of the Kubo-Greenwood response formula, and demonstrate that for thermal conductivity, the mean-field treatment of the electron-electron (e-e) interaction therein is insufficient to reproduce the weak-coupling limit obtained by plasma kinetic theories. An explicit e-e scattering correction to the DFT is posited by appealing to Matthiessen's Rule and the results of our computations of conductivities with the quantum Lenard-Balescu (QLB) equation. Further motivation of our correction is provided by an argument arising from the Zubarev quantum kinetic theory approach. Significant emphasis is placed on our efforts to produce properly converged results for plasma transport using Kohn-Sham DFT, so that an accurate assessment of the importance and efficacy of our e-e scattering corrections to the thermal conductivity can be made.

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.

Ab initio calculation of thermodynamic potentials and entropies for superionic water.

We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. Differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.

Bonding and structure in dense multi-component molecular mixtures.

We have performed finite-temperature density functional theory molecular dynamics simulations on dense methane, ammonia, and water mixtures (CH4:NH3:H2O) for various compositions and temperatures (2000 K ≤ T ≤ 10,000 K) that span a set of possible conditions in the interiors of ice-giant exoplanets. The equation-of-state, pair distribution functions, and bond autocorrelation functions (BACF) were used to probe the structure and dynamics of these complex fluids. In particular, an improvement to the choice of the cutoff in the BACF was developed that allowed analysis refinements for density and temperature effects. We note the relative changes in the nature of these systems engendered by variations in the concentration ratios. A basic tenet emerges from all these comparisons that varying the relative amounts of the three heavy components (C,N,O) can effect considerable changes in the nature of the fluid and may in turn have ramifications for the structure and composition of various planetary layers.