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

Theoretical evidence of H-He demixing under Jupiter and Saturn conditions

The immiscibility of hydrogen-helium mixture under the temperature and pressure conditions of planetary interiors is crucial for understanding the structures of gas giant planets (e.g., Jupiter and Saturn). While the experimental probe at such extreme conditions is challenging, theoretical simulation is heavily relied in an effort to unravel the mixing behavior of hydrogen and helium. Here we develop a method via a machine learning accelerated molecular dynamics simulation to quantify the physical separation of hydrogen and helium under the conditions of planetary interiors. The immiscibility line achieved with the developed method yields substantially higher demixing temperatures at pressure above 1.5 Mbar than earlier theoretical data, but matches better to the experimental estimate. Our results suggest a possibility that H-He demixing takes place in a large fraction of the interior radii of Jupiter and Saturn, i.e., 27.5% in Jupiter and 48.3% in Saturn. This indication of an H-He immiscible layer hints at the formation of helium rain and offers a potential explanation for the decrease of helium in the atmospheres of Jupiter and Saturn.

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [Nat. Commun. 13, 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Magnetic damping of jet flows in quasi-two-dimensional Rayleigh-Bénard convection

The mechanism responsible for the damping of the large-scale, azimuthally directed jets observed at Jupiter's surface is not well known, but electromagnetic forces are suspected to play a role as the planet's electrical conductivity increases radially with depth. To isolate the jet damping process, we carry out a suite of direct numerical simulations of quasi-two-dimensional, horizontally periodic Rayleigh-Bénard convection with stress-free boundary conditions in the presence of an external, vertical magnetic field. Jets, punctuated by intermittent convective bursts, develop at Rayleigh numbers (Ra, ratio of buoyancy to diffusion) beyond 10^{5} when the magnetic field is relatively weak. Five primary flow regimes are found by varying 10^{3}≤Ra≤10^{10} and the Chandrasekhar number (Ch, ratio of Lorentz to viscosity) 0≤Ch≤10^{6}: (i) steady convection rolls, (ii) steady magneto-columns, (iii) unsteady to turbulent magneto-plumes, (iv) horizontally drifting magneto-plumes, and (v) jets with intermittent turbulent convective bursts. We parse the parameter space using transitions derived from the interaction parameter (N, ratio of Lorentz to inertia). The transition to the regime dominated by jets has the most immediate applications to the magnetic damping of Jovian jet flows, where the separation between jets and a magnetically constrained system occurs at a jet-based interaction parameter value of N_{J}≈1. We approximate the value of the Jovian interaction parameter as a function of depth, and find that the jets may brake at ≈6000 km below the surface, which is deeper than recent estimates from NASA's Juno mission. This suggests that mechanisms in addition to electromagnetic forces are likely required to fully truncate the jets.

Ab Initio Study of Structure and Transport Properties of Warm Dense Nitric Oxide

The structure, equation of state and transport properties of warm dense nitric oxide (NO) were investigated in wide density and temperature ranges by ab initio molecular dynamics simulations. Both the Perdew–Burke–Ernzerhof (PBE) and the strongly constrained and appropriately normed functional with revised Vydrov–van Voorhis nonlocal correlation (SCAN−rVV10) functionals were used in the simulations, and the pressures predicted by the SCAN−rVV10 functional were found to be systematically lower than those predicted using PBE and experimental data along the shock Hugoniot curve. Along the Hugoniot curve, as density increased, we found that the system transformed towards a mixture of atomic nitrogen and oxygen liquids with molecular NO that remained present up to the highest densities explored. The electrical conductivity along Hugoniot indicated that nonmetal to metal transition had taken place. We also calculated the electrical and thermal conductivities of nitric oxide in the warm dense matter regime, and used them to compute the Lorentz number. In addition, we also report the electronic density of states.

Thermal excitation signals in the inhomogeneous warm dense electron gas

We investigate the emergence of electronic excitations from the inhomogeneous electronic structure at warm dense matter parameters based on first-principles calculations. The emerging modes are controlled by the imposed perturbation amplitude. They include satellite signals around the standard plasmon feature, transformation of plasmons to optical modes, and double-plasmon modes. These modes exhibit a pronounced dependence on the temperature. This makes them potentially invaluable for the diagnostics of plasma parameters in the warm dense matter regime. We demonstrate that these modes can be probed with present experimental techniques.

The relevance of electronic perturbations in the warm dense electron gas

Warm dense matter (WDM) has emerged as one of the frontiers of both experimental physics and theoretical physics and is a challenging traditional concept of plasma, atomic, and condensed-matter physics. While it has become common practice to model correlated electrons in WDM within the framework of Kohn-Sham density functional theory, quantitative benchmarks of exchange-correlation (XC) functionals under WDM conditions are yet incomplete. Here, we present the first assessment of common XC functionals against exact path-integral Monte Carlo calculations of the harmonically perturbed thermal electron gas. This system is directly related to the numerical modeling of x-ray scattering experiments on warm dense samples. Our assessment yields the parameter space where common XC functionals are applicable. More importantly, we pinpoint where the tested XC functionals fail when perturbations on the electronic structure are imposed. We indicate the lack of XC functionals that take into account the needs of WDM physics in terms of perturbed electronic structures.

Direct H-He chemical association in superionic FeO2H2He at Deep-Earth conditions

Hydrogen and helium are known to play crucial roles in geological and astrophysical environments; however, they are inert toward each other across wide pressure-temperature (P-T) conditions. Given their prominent presence and influence on the formation and evolution of celestial bodies, it is of fundamental interest to explore the nature of interactions between hydrogen and helium. Using an advanced crystal structure search method, we have identified a quaternary compound FeO₂H₂He stabilized in a wide range of P-T conditions. Ab initio molecular dynamics simulations further reveal a novel superionic state of FeO₂H₂He hosting liquid-like diffusive hydrogen in the FeO₂He sublattice, creating a conducive environment for H-He chemical association, at P-T conditions corresponding to the Earth's lowest mantle regions. To our surprise, this chemically facilitated coalescence of otherwise immiscible molecular species highlights a promising avenue for exploring this long-sought but hitherto unattainable state of matter. This finding raises strong prospects for exotic H-He mixtures inside Earth and possibly also in other astronomical bodies.

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.

Evidence of hydrogen-helium immiscibility at Jupiter-interior conditions

The phase behaviour of warm dense hydrogen-helium (H-He) mixtures affects our understanding of the evolution of Jupiter and Saturn and their interior structures^1,2. For example, precipitation of He from a H-He atmosphere at about 1-10 megabar and a few thousand kelvin has been invoked to explain both the excess luminosity of Saturn^1,3, and the depletion of He and neon (Ne) in Jupiter's atmosphere as observed by the Galileo probe^4,5. But despite its importance, H-He phase behaviour under relevant planetary conditions remains poorly constrained because it is challenging to determine computationally and because the extremes of temperature and pressure are difficult to reach experimentally. Here we report that appropriate temperatures and pressures can be reached through laser-driven shock compression of H₂-He samples that have been pre-compressed in diamond-anvil cells. This allows us to probe the properties of H-He mixtures under Jovian interior conditions, revealing a region of immiscibility along the Hugoniot. A clear discontinuous change in sample reflectivity indicates that this region ends above 150 gigapascals at 10,200 kelvin and that a more subtle reflectivity change occurs above 93 gigapascals at 4,700 kelvin. Considering pressure-temperature profiles for Jupiter, these experimental immiscibility constraints for a near-protosolar mixture suggest that H-He phase separation affects a large fraction-we estimate about 15 per cent of the radius-of Jupiter's interior. This finding provides microphysical support for Jupiter models that invoke a layered interior to explain Juno and Galileo spacecraft observations^1,4,6-8.

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.

Multishock to Quasi-Isentropic Compression of Dense Gaseous Deuterium-Helium Mixtures up to 120 GPa: Probing the Sound Velocities Relevant to Planetary Interiors

Shock reverberation compression experiments on dense gaseous deuterium-helium mixtures are carried out to provide thermodynamic parameters relevant to the conditions in planetary interiors. The multishock pressures are determined up to 120 GPa and reshock temperatures to 7400 K. Furthermore, the unique compression path from shock-adiabatic to quasi-isentropic compressions enables a direct estimation of the high-pressure sound velocities in the unexplored range of 50-120 GPa. The equation of state and sound velocity provide particular dual perspectives to validate the theoretical models. Our experimental data are found to agree with several equation of state models widely used in astrophysics within the probed pressure range. The current data improve the experimental constraints on sound velocities in the Jovian insulating-to-metallic transition layer.

Effective Static Approximation: A Fast and Reliable Tool for Warm-Dense Matter Theory

We present an effective static approximation (ESA) to the local field correction (LFC) of the electron gas that enables highly accurate calculations of electronic properties like the dynamic structure factor S(q,ω), the static structure factor S(q), and the interaction energy v. The ESA combines the recent neural-net representation by T. Dornheim et al., [J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013] of the temperature-dependent LFC in the exact static limit with a consistent large wave-number limit obtained from quantum Monte Carlo data of the on-top pair distribution function g(0). It is suited for a straightforward integration into existing codes. We demonstrate the importance of the LFC for practical applications by reevaluating the results of the recent x-ray Thomson scattering experiment on aluminum by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.115001]. We find that an accurate incorporation of electronic correlations in terms of the ESA leads to a different prediction of the inelastic scattering spectrum than obtained from state-of-the-art models like the Mermin approach or linear-response time-dependent density functional theory. Furthermore, the ESA scheme is particularly relevant for the development of advanced exchange-correlation functionals in density functional theory.

Coexistence of plastic and partially diffusive phases in a helium-methane compound

Helium and methane are major components of giant icy planets and are abundant in the universe. However, helium is the most inert element in the periodic table and methane is one of the most hydrophobic molecules, thus whether they can react with each other is of fundamental importance. Here, our crystal structure searches and first-principles calculations predict that a He3CH4 compound is stable over a wide range of pressures from 55 to 155 GPa and a HeCH4 compound becomes stable around 105 GPa. As nice examples of pure van der Waals crystals, the insertion of helium atoms changes the original packing of pure methane molecules and also largely hinders the polymerization of methane at higher pressures. After analyzing the diffusive properties during the melting of He3CH4 at high pressure and high temperature, in addition to a plastic methane phase, we have discovered an unusual phase which exhibits coexistence of diffusive helium and plastic methane. In addition, the range of the diffusive behavior within the helium-methane phase diagram is found to be much narrower compared to that of previously predicted helium-water compounds. This may be due to the weaker van der Waals interactions between methane molecules compared to those in helium-water compounds, and that the helium-methane compound melts more easily.

Ion energy-loss characteristics and friction in a free-electron gas at warm dense matter and nonideal dense plasma conditions

We investigate the energy-loss characteristics of an ion in warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio quantum Monte Carlo- (QMC) based machine learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of γ/ω_{pi}=0.01, where ω_{pi} is the ionic plasma frequency. This analysis is performed by comparing QMC-based data to results from the random-phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjölander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss of ions in a correlated partially degenerate electrons gas. Moreover, by comparing QMC-based data with the results obtained using STLS, we find that the ion energy-loss properties are not sensitive to the inaccuracy of the static local field correction (LFC) at large wave numbers, k/k_{F}>2 (with k_{F} being the Fermi wave number), but that a correct description of the static LFC at k/k_{F}≲1.5 is important.

Ab initio study of the miscibility for solid hydrogen-helium mixtures at high pressure

Understanding the behavior of H₂-He binary mixtures at high pressure is of great importance. Two more recent experiments [J. Lim and C. S. Yoo, Phys. Rev. Lett. 120, 165301 (2018) and R. Turnbull et al., ibid. 121, 195702 (2018)] are in conflict, regarding the miscibility between H₂ and He in solids at high pressure. On the basis of first-principles calculations combined with the structure prediction method, we investigate the miscibility for solid H₂-He mixtures at pressures from 0 GPa to 200 GPa. It is found that there is no sign of miscibility and chemical reactivity in H₂-He mixtures with any H:He ratio. Moreover, instead of H₂-He mixtures, the calculated Raman modes of the N-H mixtures can better explain the characteristic peaks observed experimentally, which were claimed to be the H-He vibrational modes. These calculation results are more in line with the experimental findings by Turnbull et al. [Phys. Rev. Lett. 121, 195702 (2018)].

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.

Ab initio Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter

The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor S(q,ω) constitutes a key quantity as it is directly measured in x-ray Thomson scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first ab initio results for S(q,ω) by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction G(q,ω). In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for S(q,ω) over substantial parts of the warm dense matter regime.

Liquid Structure of Shock-Compressed Hydrocarbons at Megabar Pressures

We present results for the ionic structure in hydrocarbons (polystyrene, polyethylene) that were shock compressed to pressures of up to 190 GPa, inducing rapid melting of the samples. The structure of the resulting liquid is then probed using in situ diffraction by an x-ray free electron laser beam, demonstrating the capability to obtain reliable diffraction data in a single shot, even for low-Z samples without long range order. The data agree well with ab initio simulations, validating the ability of such approaches to model mixed samples in states where complex interparticle bonds remain, and showing that simpler models are not necessarily valid. While the results clearly exclude the possibility of complete carbon-hydrogen demixing at the conditions probed, they also, in contrast to previous predictions, indicate that diffraction is not always a sufficient diagnostic for this phenomenon.

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.

Quantum molecular dynamics study on the proton exchange, ionic structures, and transport properties of warm dense hydrogen-deuterium mixtures

Comprehensive knowledge of physical properties such as equation of state (EOS), proton exchange, dynamic structures, diffusion coefficients, and viscosities of hydrogen-deuterium mixtures with densities from 0.1 to 5 g/cm^{3} and temperatures from 1 to 50 kK has been presented via quantum molecular dynamics (QMD) simulations. The existing multi-shock experimental EOS provides an important benchmark to evaluate exchange-correlation functionals. The comparison of simulations with experiments indicates that a nonlocal van der Waals density functional (vdW-DF1) produces excellent results. Fraction analysis of molecules using a weighted integral over pair distribution functions was performed. A dissociation diagram together with a boundary where the proton exchange (H_{2}+D_{2}⇌2HD) occurs was generated, which shows evidence that the HD molecules form as the H_{2} and D_{2} molecules are almost 50% dissociated. The mechanism of proton exchange can be interpreted as a process of dissociation followed by recombination. The ionic structures at extreme conditions were analyzed by the effective coordination number model. High-order cluster, circle, and chain structures can be founded in the strongly coupled warm dense regime. The present QMD diffusion coefficient and viscosity can be used to benchmark two analytical one-component plasma (OCP) models: the Coulomb and Yukawa OCP models.