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.

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.

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.

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.

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.

Superionic Phases of the 1:1 Water-Ammonia Mixture.

We report four structures for the 1:1 water-ammonia mixture showing superionic behavior at high temperature with the space groups P4/nmm, Ima2, Pma2, and Pm, which have been identified from evolutionary random structure search calculations at 0 K. Analyzing the respective pair distribution functions and diffusive properties the superionic phase is found to be stable in a temperature range between 1000 and 6000 K for pressures up to 800 GPa. We propose a high-pressure phase diagram of the water-ammonia mixture for the first time and compare the self-diffusion coefficients in the mixture to the ones found in water and ammonia. Finally, possible implications on the interior structure of the giant planets Uranus and Neptune 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.

Ab initio simulations for the ion-ion structure factor of warm dense aluminum.

We perform ab initio simulations based on finite-temperature density functional theory in order to determine the static and dynamic ion-ion structure factor in aluminum. We calculate the dynamic structure factor via the intermediate scattering function and extract the dispersion relation for the collective excitations. The results are compared with available experimental x-ray scattering data. Very good agreement is obtained for the liquid metal domain. In addition we perform simulations for warm dense aluminum in order to obtain the ion dynamics in this strongly correlated quantum regime. We determine the sound velocity for both liquid and warm dense aluminum which can be checked experimentally using narrow-bandwidth free electron laser radiation.

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.

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.

Equation of state and phase diagram of ammonia at high pressures from ab initio simulations.

We present an equation of state as well as a phase diagram of ammonia at high pressures and high temperatures derived from ab initio molecular dynamics simulations. The predicted phases of ammonia are characterized by analyzing diffusion coefficients and structural properties. Both the phase diagram and the subsequently computed Hugoniot curves are compared to experimental results. Furthermore, we discuss two methods that allow us to take into account nuclear quantum effects, which are of considerable importance in molecular fluids. Our data cover pressures up to 330 GPa and a temperature range from 500 K to 10,000 K. This regime is of great interest for interior models of the giant planets Uranus and Neptune, which contain, besides water and methane, significant amounts of ammonia.

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.

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.

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.

Dynamical screening and ionic conductivity in water from ab initio simulations.

We present a method to calculate ionic conductivities of complex fluids from ab initio simulations. This is achieved by combining density functional theory molecular dynamics simulations with polarization theory. Conductivities are then obtained via a Green-Kubo formula using time-dependent effective charges of electronically screened ions. The method is applied to two different phases of warm dense water. We observe large fluctuations in the effective charges; protons can transport effective charges greater than +e for ultrashort time scales. Furthermore, we compare our results with a simpler model of ionic conductivity in water that is based on diffusion coefficients. Our approach can be directly applied to study ionic conductivities of electronically insulating materials of arbitrary composition, e.g., complex molecular mixtures under such extreme conditions that occur deep inside giant planets.

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.

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.

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.

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.

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.