Liquid-liquid phase transition in hydrogen by coupled electron-ion Monte Carlo simulations

Bound-State Breaking and the Importance of Thermal Exchange-Correlation Effects in Warm Dense Hydrogen

Hydrogen at extreme temperatures and pressures is of key relevance for cutting-edge technological applications, with inertial confinement fusion research being a prime example. In addition, it is ubiquitous throughout our universe and naturally occurs in a variety of astrophysical objects. In the present work, we present exact ab initio path integral Monte Carlo (PIMC) results for the electronic density of warm dense hydrogen along a line of constant degeneracy across a broad range of densities. Using the well-known concept of reduced density gradients, we develop a new framework to identify the breaking of bound states due to pressure ionization in bulk hydrogen. Moreover, we use our PIMC results as a reference to rigorously assess the accuracy of a variety of exchange-correlation (XC) functionals in density functional theory calculations for different density regions. Here, a key finding is the importance of thermal XC effects for the accurate description of density gradients in high-energy-density systems. Our exact PIMC test set is freely available online and can be used to guide the development of new methodologies for the simulation of warm dense matter and beyond.

Deep Variational Free Energy Approach to Dense Hydrogen

We developed a deep generative model-based variational free energy approach to the equations of state of dense hydrogen. We employ a normalizing flow network to model the proton Boltzmann distribution and a fermionic neural network to model the electron wave function at given proton positions. By jointly optimizing the two neural networks we reached a comparable variational free energy to the previous coupled electron-ion Monte Carlo calculation. The predicted equation of state of dense hydrogen under planetary conditions is denser than the findings of ab initio molecular dynamics calculation and empirical chemical model. Moreover, direct access to the entropy and free energy of dense hydrogen opens new opportunities in planetary modeling and high-pressure physics research.

Electronic density response of warm dense hydrogen on the nanoscale

The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive ab initio path integral Monte Carlo results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us unique insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Exciton Nature of Plasma Phase Transition in Warm Dense Fluid Hydrogen: ROKS Simulation

The transition of warm dense fluid hydrogen from an insulator to a conducting state at pressures of about 20-400 GPa and temperatures of 500-5000 K has been the subject of active scientific research over the past few decades. However, various experimental and theoretical methods do not provide consistent results. In this work, we have applied the restricted open-shell Kohn-Sham (ROKS) method for first principles molecular dynamics of dense hydrogen after thermal excitation to the first singlet excited state. The Wannier localization method has allowed us to analyze the exciton dynamics in this system. The model shows that a key mechanism of the transition is associated with the dissociation of electron-hole pairs, which allows explaining several stages of the transition of fluid H₂ from molecular state to plasma. This mechanism is able to give a quantitative description of several experimental results as well as to resolve the discrepancies between experimental studies.

Stable Solid Molecular Hydrogen above 900 K from a Machine-Learned Potential Trained with Diffusion Quantum Monte Carlo

We survey the phase diagram of high-pressure molecular hydrogen with path integral molecular dynamics using a machine-learned interatomic potential trained with quantum Monte Carlo forces and energies. Besides the HCP and C2/c-24 phases, we find two new stable phases both with molecular centers in the Fmmm-4 structure, separated by a molecular orientation transition with temperature. The high temperature isotropic Fmmm-4 phase has a reentrant melting line with a maximum at higher temperature (1450 K at 150 GPa) than previously estimated and crosses the liquid-liquid transition line around 1200 K and 200 GPa.

Ab Initio Static Exchange-Correlation Kernel across Jacob's Ladder without Functional Derivatives

The electronic exchange─correlation (XC) kernel constitutes a fundamental input for the estimation of a gamut of properties such as the dielectric characteristics, the thermal and electrical conductivity, or the response to an external perturbation. In this work, we present a formally exact methodology for the computation of the system specific static XC kernel exclusively within the framework of density functional theory (DFT) and without employing functional derivatives─no external input apart from the usual XC-functional is required. We compare our new results with exact quantum Monte Carlo (QMC) data for the archetypical uniform electron gas model under both ambient and warm dense matter conditions. This gives us unprecedented insights into the performance of different XC functionals, and it has important implications for the development of new functionals that are designed for the application at extreme temperatures. In addition, we obtain new DFT results for the XC kernel of warm dense hydrogen as it occurs in fusion applications and astrophysical objects. The observed excellent agreement to the QMC reference data demonstrates that presented framework is capable to capture nontrivial effects such as XC-induced isotropy breaking in the density response of hydrogen at large wave numbers.

Ab initio path integral Monte Carlo simulations of hydrogen snapshots at warm dense matter conditions

We combine ab initio path integral Monte Carlo (PIMC) simulations with fixed ion configurations from density functional theory molecular dynamics (DFT-MD) simulations to solve the electronic problem for hydrogen under warm dense matter conditions [Böhme et al., Phys. Rev. Lett. 129, 066402 (2022)0031-900710.1103/PhysRevLett.129.066402]. The problem of path collapse due to the Coulomb attraction is avoided by utilizing the pair approximation, which is compared against the simpler Kelbg pair potential. We find very favorable convergence behavior towards the former. Since we do not impose any nodal restrictions, our PIMC simulations are afflicted with the notorious fermion sign problem, which we analyze in detail. While computationally demanding, our results constitute an exact benchmark for other methods and approximations within DFT. Our setup gives us the unique capability to study important properties of warm dense hydrogen such as the electronic static density response and exchange-correlation kernel without any model assumptions, which will be very valuable for a variety of applications such as the interpretation of experiments and the development of new XC functionals.

Anomalous properties in the potential energy landscape of a monatomic liquid across the liquid-gas and liquid-liquid phase transitions

As a liquid approaches the gas state, the properties of the potential energy landscape (PEL) sampled by the system become anomalous. Specifically, (i) the mechanically stable local minima of the PEL [inherent structures (IS)] can exhibit cavitation above the so-called Sastry volume, vS, before the liquid enters the gas phase. In addition, (ii) the pressure of the liquid at the sampled IS [i.e., the PEL equation of state, PIS(v)] develops a spinodal-like minimum at vS. We perform molecular dynamics simulations of a monatomic water-like liquid and verify that points (i) and (ii) hold at high temperatures. However, at low temperatures, cavitation in the liquid and the corresponding IS occurs simultaneously and a Sastry volume cannot be defined. Remarkably, at intermediate/high temperatures, the IS of the liquid can exhibit crystallization, i.e., the liquid regularly visits the regions of the PEL that belong to the crystal state. The model liquid considered also exhibits a liquid-liquid phase transition (LLPT) between a low-density and a high-density liquid (LDL and HDL). By studying the behavior of PIS(v) during the LLPT, we identify a Sastry volume for both LDL and HDL. The HDL Sastry volume marks the onset above which IS are heterogeneous (composed of LDL and HDL particles), analogous to points (i) and (ii) above. However, the relationship between the LDL Sastry volume and the onset of heterogeneous IS is less evident. We conclude by presenting a thermodynamic argument that can explain the behavior of the PEL equation of state PIS(v) across both the liquid-gas phase transition and LLPT.

Thermodynamic modeling of fluid polyamorphism in hydrogen at extreme conditions

Fluid polyamorphism, the existence of multiple amorphous fluid states in a single-component system, has been observed or predicted in a variety of substances. A remarkable example of this phenomenon is the fluid–fluid phase transition (FFPT) in high-pressure hydrogen between insulating and conducting high-density fluids. This transition is induced by the reversible dimerization/dissociation of the molecular and atomistic states of hydrogen. In this work, we present the first attempt to thermodynamically model the FFPT in hydrogen at extreme conditions. Our predictions for the phase coexistence and the reaction equilibrium of the two alternative forms of fluid hydrogen are based on experimental data and supported by the results of simulations. Remarkably, we find that the law of corresponding states can be utilized to construct a unified equation of state combining the available computational results for different models of hydrogen and the experimental data.

Static Electronic Density Response of Warm Dense Hydrogen: Ab Initio Path Integral Monte Carlo Simulations

The properties of hydrogen under extreme conditions are important for many applications, including inertial confinement fusion and astrophysical models. A key quantity is given by the electronic density response to an external perturbation, which is probed in x-ray Thomson scattering experiments-the state of the art diagnostics from which system parameters like the free electron density n_{e}, the electronic temperature T_{e}, and the charge state Z can be inferred. In this work, we present highly accurate path integral Monte Carlo results for the static electronic density response of hydrogen. We obtain the static exchange-correlation (XC) kernel K_{XC}, which is of central relevance for many applications, such as time-dependent density functional theory. This gives us a first unbiased look into the electronic density response of hydrogen in the warm-dense matter regime, thereby opening up a gamut of avenues for future research.

Modeling fluid polyamorphism through a maximum-valence approach

We suggest a simple model to describe polyamorphism in single-component fluids using a maximum-valence approach. The model contains three types of interactions: (i) Atoms attract each other by van der Waals forces that generate a liquid-gas transition at low pressures, (ii) atoms may form covalent bonds that induce association, and (iii) atoms with maximal valence attract or repel each other stronger than other atoms, thus generating liquid-liquid separation. As an example, we qualitatively compare this model with the behavior of liquid sulfur and show that condition (iii) generates a liquid-liquid phase transition in addition to the liquid-gas phase transition.

Nucleation of Helium in Liquid Lithium at 843 K and High Pressures

Fusion energy stands out as a promising alternative for a future decarbonised energy system. In order to be sustainable, future fusion nuclear reactors will have to produce their own tritium. In the so-called breeding blanket of a reactor, the neutron bombardment of lithium will produce the desired tritium, but also helium, which can trigger nucleation mechanisms owing to the very low solubility of helium in liquid metals. An understanding of the underlying microscopic processes is important for improving the efficiency, sustainability and reliability of the fusion energy conversion process. The spontaneous creation of helium droplets or bubbles in the liquid metal used as breeding material in some designs may be a serious issue for the performance of the breeding blankets. This phenomenon has yet to be fully studied and understood. This work aims to provide some insight on the behaviour of lithium and helium mixtures at experimentally corresponding operating conditions (843 K and pressures between 108 and 1010 Pa). We report a microscopic study of the thermodynamic, structural and dynamical properties of lithium–helium mixtures, as a first step to the simulation of the environment in a nuclear fusion power plant. We introduce a new microscopic model devised to describe the formation of helium droplets in the thermodynamic range considered. Our model predicts the formation of helium droplets at pressures around 109 Pa, with radii between 1 and 2 Å. The diffusion coefficient of lithium (2 Å2/ps) is in excellent agreement with reference experimental data, whereas the diffusion coefficient of helium is in the range of 1 Å2/ps and tends to decrease as pressure increases.

The 2021 room-temperature superconductivity roadmap

Designing materials with advanced functionalities is the main focus of contemporary solid-state physics and chemistry. Research efforts worldwide are funneled into a few high-end goals, one of the oldest, and most fascinating of which is the search for an ambient temperature superconductor (A-SC). The reason is clear: superconductivity at ambient conditions implies being able to handle, measure and access a single, coherent, macroscopic quantum mechanical state without the limitations associated with cryogenics and pressurization. This would not only open exciting avenues for fundamental research, but also pave the road for a wide range of technological applications, affecting strategic areas such as energy conservation and climate change. In this roadmap we have collected contributions from many of the main actors working on superconductivity, and asked them to share their personal viewpoint on the field. The hope is that this article will serve not only as an instantaneous picture of the status of research, but also as a true roadmap defining the main long-term theoretical and experimental challenges that lie ahead. Interestingly, although the current research in superconductor design is dominated by conventional (phonon-mediated) superconductors, there seems to be a widespread consensus that achieving A-SC may require different pairing mechanisms.In memoriam, to Neil Ashcroft, who inspired us all.

Isotope Quantum Effects in the Metallization Transition in Liquid Hydrogen

Quantum effects in condensed matter normally only occur at low temperatures. Here we show a large quantum effect in high-pressure liquid hydrogen at thousands of Kelvins. We show that the metallization transition in hydrogen is subject to a very large isotope effect, occurring hundreds of degrees lower than the equivalent transition in deuterium. We examined this using path integral molecular dynamics simulations which identify a liquid-liquid transition involving atomization, metallization, and changes in viscosity, specific heat, and compressibility. The difference between H_{2} and D_{2} is a quantum mechanical effect that can be associated with the larger zero-point energy in H_{2} weakening the covalent bond. Our results mean that experimental results on deuterium must be corrected before they are relevant to understanding hydrogen at planetary conditions.

First-principles equation of state of liquid hydrogen and dissociative transition

The properties of dense hot hydrogen, in particular the phase transition between the molecular insulating and atomic conductive states, are important in the fields of astrophysics and high-pressure physics. Previous ab initio calculations suggested the metallization in liquid hydrogen, accompanied by dissociation, is a first-order phase transition and ends at a critical point in temperature range between 1500 and 2000 K and pressure close to 100 GPa. Using density functional theoretical molecular dynamics simulations, we report a first-principles equation of state of hydrogen that covers dissociation transition conditions at densities ranging from 0.20 to 1.00 g/cc and temperatures of 600-9000 K. Our results clearly indicate that a drop in pressure and a sharp structural change still occur as the system transforms from a diatomic to monoatomic phase at temperatures above 2000 K, and support the first-order phase transition in liquid hydrogen would end in the temperature about 4500 K.

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

Energy Gap Closure of Crystalline Molecular Hydrogen with Pressure

We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by coupled electron ion Monte Carlo methods. Nuclear zero point effects cause a large reduction in the gap (∼2  eV). Depending on the structure, the fundamental indirect gap closes between 380 and 530 GPa for ideal crystals and 330-380 GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to ∼450-500  GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

A Spectroscopic Study of the Insulator-Metal Transition in Liquid Hydrogen and Deuterium

The insulator-to-metal transition in dense fluid hydrogen is an essential phenomenon in the study of gas giant planetary interiors and the physical and chemical behavior of highly compressed condensed matter. Using direct fast laser spectroscopy techniques to probe hydrogen and deuterium precompressed in a diamond anvil cell and laser heated on microsecond timescales, an onset of metal-like reflectance is observed in the visible spectral range at P >150 GPa and T ≥ 3000 K. The reflectance increases rapidly with decreasing photon energy indicating free-electron metallic behavior with a plasma edge in the visible spectral range at high temperatures. The reflectance spectra also suggest much longer electronic collision time (≥1 fs) than previously inferred, implying that metallic hydrogen at the conditions studied is not in the regime of saturated conductivity (Mott-Ioffe-Regel limit). The results confirm the existence of a semiconducting intermediate fluid hydrogen state en route to metallization.

First-order liquid-liquid phase transition in compressed hydrogen and critical point

One of the fundamental problems relating to the properties of hydrogen is that of insulator-metal transition. Recent theoretical and experimental studies show that the metallization in liquid hydrogen could be a first-order phase transition and involve molecular to atomic transition. However, the location of the critical point is still an unresolved question. Earlier studies reported the critical point at a temperature of 1500-2000 K, but recent experimental observations on diamond-anvil cells show that the discontinuous transition still persists at temperatures well above 2000 K. We have carried out a detailed study on the liquid-liquid phase transition in dense hydrogen by uisng ab initio molecular dynamics simulations and found new evidence for the abrupt metallization between weakly dissociated and strongly dissociated fluid phases at temperatures as high as 3000 and 4000 K. Also, the predicted phase boundary is in excellent agreement with the recent experiments. Our results suggest that this first-order transition in liquid hydrogen likely ends in a critical point around 4000 K, which is significantly higher than the previous theoretical predictions.

Optical properties of high-pressure fluid hydrogen across molecular dissociation

Optical properties of compressed fluid hydrogen in the region where dissociation and metallization is observed are computed by ab initio methods and compared with recent experimental results. We confirm that at T > 3,000 K, both processes are continuous, while at T < 1,500 K, the first-order phase transition is accompanied by a discontinuity of the dc conductivity and the thermal conductivity, while both the reflectivity and absorption coefficient vary rapidly but continuously. Our results support the recent analysis of National Ignition Facility (NIF) experiments [Celliers PM, et al. (2018) Science 361:677-682], which assigned the inception of metallization to pressures where the reflectivity is ∼0.3. Our results also support the conclusion that the temperature plateau seen in laser-heated diamond-anvil cell (DAC) experiments at temperatures higher than 1,500 K corresponds to the onset of optical absorption, not to the phase transition.

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.

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.

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.

Theoretical and experimental investigation of the equation of state of boron plasmas

We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1×10^{4}-5.2×10^{8} K) and densities (0.25-49 g/cm^{3}) and experimental shock Hugoniot data at unprecedented high pressures (5608±118 GPa). The calculations are performed with first-principles methods combining path-integral Monte Carlo (PIMC) at high temperatures and density-functional-theory molecular-dynamics (DFT-MD) methods at lower temperatures. PIMC and DFT-MD cross-validate each other by providing coherent EOS (difference <1.5 Hartree/boron in energy and <5% in pressure) at 5.1×10^{5} K. The Hugoniot measurement is conducted at the National Ignition Facility using a planar shock platform. The pressure-density relation found in our shock experiment is on top of the shock Hugoniot profile predicted with our first-principles EOS and a semiempirical EOS table (LEOS 50). We investigate the self-diffusivity and the effect of thermal and pressure-driven ionization on the EOS and shock compression behavior in high-pressure and -temperature conditions. We also study the sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on applying pressure multipliers to LEOS 50 and by utilizing a new EOS model based on our ab initio simulations via one-dimensional radiation-hydrodynamic calculations. The results are valuable for future theoretical and experimental studies and engineering design in high-energy density research.

Metallization and molecular dissociation of dense fluid nitrogen

Diatomic nitrogen is an archetypal molecular system known for its exceptional stability and complex behavior at high pressures and temperatures, including rich solid polymorphism, formation of energetic states, and an insulator-to-metal transformation coupled to a change in chemical bonding. However, the thermobaric conditions of the fluid molecular–polymer phase boundary and associated metallization have not been experimentally established. Here, by applying dynamic laser heating of compressed nitrogen and using fast optical spectroscopy to study electronic properties, we observe a transformation from insulating (molecular) to conducting dense fluid nitrogen at temperatures that decrease with pressure and establish that metallization, and presumably fluid polymerization, occurs above 125 GPa at 2500 K. Our observations create a better understanding of the interplay between molecular dissociation, melting, and metallization revealing features that are common in simple molecular systems.

Nitrogen is a model system still presenting unknown behaviors at the pressures and temperatures typical of deep planets’ interiors. Here the authors explore, by pulsed laser heating in a diamond anvil cell and optical measurements, the metallization and non-molecular states of nitrogen in a previously unexplored domain above 1 Mbar and at 2000-7000K.

Metal-Insulator Transition of Solid Hydrogen by the Antisymmetric Shadow Wave Function

We revisit the pressure-induced molecular-atomic metal-insulator transition of solid hydrogen by means of variational quantum Monte Carlo simulations based on the antisymmetric shadow wave function. For the purpose of facilitating the study of the electronic structure of large-scale fermionic systems, the shadow wave function formalism is extended by a series of technical advancements as implemented in our HswfQMC code. Among others, these improvements include a revised optimization method for the employed shadow wave function and an enhanced treatment of periodic systems with long-range interactions. It is found that the superior accuracy of the antisymmetric shadow wave function results in a significantly increased transition pressure with respect to previous theoretical estimates.

Anomalous phase behavior of first-order fluid-liquid phase transition in phosphorus

Although the existence of liquid-liquid phase transition has become more and more convincing, whether it will terminate at a critical point and what is the order parameter are still open. To explore these questions, we revisit the fluid-liquid phase transition (FLPT) in phosphorus (P) and study its phase behavior by performing extensive first-principles molecular dynamics simulations. The FLPT observed in experiments is well reproduced, and a fluid-liquid critical point (FLCP) at T = 3000 ∼ 3500 K, P = 1.5-2.0 Kbar is found. With decreasing temperature from the FLCP along the transition line, the density difference (Δρ) between two coexisting phases first increases from zero and then anomalously decreases; however, the entropy difference (ΔS) continuously increases from zero. These features suggest that an order parameter containing contributions from both the density and the entropy is needed to describe the FLPT in P, and at least at low temperatures, the entropy, instead of the density, governs the FLPT.

Dynamic conductivity and partial ionization in dense fluid hydrogen

A theoretical description for optical conduction experiments in dense fluid hydrogen is presented. Different quantum statistical approaches are used to describe the mechanism of electronic transport in hydrogen's high-temperature dense phase. We show that at the onset of the metallic transition, optical conduction could be described by a strong rise in atomic polarizability, due to increased ionization, whereas in the highly degenerate limit, the Ziman weak scattering model better accounts for the observed saturation of reflectance. The inclusion of effects of partial ionization in the highly degenerate region provides great agreement with experimental results. Hydrogen's fluid metallic state is revealed to be a partially ionized free-electron plasma. Our results provide some of the first theoretical transport models that are experimentally benchmarked, as well as an important guide for future studies.

Dynamic electron-ion collisions and nuclear quantum effects in quantum simulation of warm dense matter

The structural, thermodynamic and transport properties of warm dense matter (WDM) are crucial to the fields of astrophysics and planet science, as well as inertial confinement fusion. WDM refers to the states of matter in a regime of temperature and density between cold condensed matter and hot ideal plasmas, where the density is from near-solid up to ten times solid density, and the temperature between 0.1 and 100 eV. In the WDM regime, matter exhibits moderately or strongly coupled, partially degenerate properties. Therefore, the methods used to deal with condensed matter and isolated atoms need to be properly validated for WDM. It is therefore a big challenge to understand WDM within a unified theoretical description with reliable accuracy. Here, we review the progress in the theoretical study of WDM with state-of-the-art simulations, i.e. quantum Langevin molecular dynamics and first principles path integral molecular dynamics. The related applications for WDM are also included.

Phase Diagram of Hydrogen and a Hydrogen-Helium Mixture at Planetary Conditions by Quantum Monte Carlo Simulations

Understanding planetary interiors is directly linked to our ability of simulating exotic quantum mechanical systems such as hydrogen (H) and hydrogen-helium (H-He) mixtures at high pressures and temperatures. Equation of state (EOS) tables based on density functional theory are commonly used by planetary scientists, although this method allows only for a qualitative description of the phase diagram. Here we report quantum Monte Carlo (QMC) molecular dynamics simulations of pure H and H-He mixture. We calculate the first QMC EOS at 6000 K for a H-He mixture of a protosolar composition, and show the crucial influence of He on the H metallization pressure. Our results can be used to calibrate other EOS calculations and are very timely given the accurate determination of Jupiter's gravitational field from the NASA Juno mission and the effort to determine its structure.

Ab initio Exchange-Correlation Free Energy of the Uniform Electron Gas at Warm Dense Matter Conditions

In a recent Letter [T. Dornheim et al., Phys. Rev. Lett. 117, 156403 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.156403], we presented the first quantum Monte Carlo (QMC) results for the warm dense electron gas in the thermodynamic limit. However, a complete parametrization of the exchange-correlation free energy with respect to density, temperature, and spin polarization remained out of reach due to the absence of (i) accurate QMC results below θ=k_{B}T/E_{F}=0.5 and (ii) QMC results for spin polarizations different from the paramagnetic case. Here we overcome both remaining limitations. By closing the gap to the ground state and by performing extensive QMC simulations for different spin polarizations, we are able to obtain the first completely ab initio exchange-correlation free energy functional; the accuracy achieved is an unprecedented ∼0.3%. This also allows us to quantify the accuracy and systematic errors of various previous approximate functionals.

Melting and High P-T Transitions of Hydrogen up to 300 GPa

High P-T Raman spectra of hydrogen in the vibron and lattice mode regions were measured up to 300 GPa and 900 K using externally heated diamond anvil cell techniques. A new melting line determined from the disappearance of lattice mode excitations was measured directly for the first time above 140 GPa. The results differ from theoretical predictions and extrapolations from lower pressure melting relations. In addition, discontinuities in Raman frequencies are observed as a function of pressure and temperature indicative of phase transition at these conditions. The appearance of a new Raman feature near 2700  cm^{-1} at ∼300  GPa and 370 K indicates the transformation to a new crystalline phase. Theoretical calculations of the spectrum suggest the new phase is the proposed Cmca-4 metallic phase. The transition pressure is close to that of a recently reported transition observed on dynamic compression.

Accelerating ab initio Molecular Dynamics and Probing the Weak Dispersive Forces in Dense Liquid Hydrogen

We propose an ab initio molecular dynamics method, capable of dramatically reducing the autocorrelation time required for the simulation of classical and quantum particles at finite temperatures. The method is based on an efficient implementation of a first order Langevin dynamics modified by means of a suitable, position dependent acceleration matrix S. Here, we apply this technique to both Lennard-Jones models, to demonstrate the accuracy and speeding-up of the sampling, and within a quantum Monte Carlo based wave function approach, for determining the phase diagram of high-pressure hydrogen with simulations much longer than the autocorrelation time. With the proposed method, we are able to equilibrate in a few hundred steps even close to the liquid-liquid phase transition (LLT). Within our approach, we find that the LLT transition is consistent with recent density functionals predicting a much larger transition pressure when the long range dispersive forces are taken into account.