Effective electronic forces and potentials from ab initio path integral Monte Carlo simulations

Modeling of warm dense hydrogen via explicit real-time electron dynamics: Dynamic structure factors

We present two methods for computing the dynamic structure factor for warm dense hydrogen without invoking either the Born-Oppenheimer approximation or the Chihara decomposition, by employing a wave-packet description that resolves the electron dynamics during ion evolution. First, a semiclassical method is discussed, which is corrected based on known quantum constraints, and second, a direct computation of the density response function within the molecular dynamics. The wave-packet models are compared to PIMC and DFT-MD for the static and low-frequency behavior. For the high-frequency behavior the models recover the expected behavior in the limits of small and large momentum transfers and show the characteristic flattening of the plasmon dispersion for intermediate momentum transfers due to interactions, in agreement with commonly used models for x-ray Thomson scattering. By modeling the electrons and ions on an equal footing, both the ion and free electron part of the spectrum can now be treated within a single framework where we simultaneously resolve the ion-acoustic and plasmon mode, with a self-consistent description of collisions and screening.

Ab Initio Path Integral Monte Carlo Simulations of the Uniform Electron Gas on Large Length Scales

The accurate description of non-ideal quantum many-body systems is of prime importance for a host of applications within physics, quantum chemistry, materials science, and related disciplines. At finite temperatures, the gold standard is given by ab initio path integral Monte Carlo (PIMC) simulations, which do not require any empirical input but exhibit an exponential increase in the required computation time for Fermionic systems with an increase in system size N. Very recently, computing Fermionic properties without this bottleneck based on PIMC simulations of fictitious identical particles has been suggested. In our work, we use this technique to perform very large (N ≤ 1000) PIMC simulations of the warm dense electron gas and demonstrate that it is capable of providing a highly accurate description of the investigated properties, i.e., the static structure factor, the static density response function, and the local field correction, over the entire range of length scales.

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.

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'.

Quantum version of the integral equation theory-based dielectric scheme for strongly coupled electron liquids

A novel dielectric scheme is proposed for strongly coupled electron liquids, which handles quantum mechanical effects beyond the random phase approximation level and treats electronic correlations within the integral equation theory of classical liquids. The self-consistent scheme features a complicated dynamic local field correction functional and its formulation is guided by ab initio path integral Monte Carlo simulations. Remarkably, our scheme is capable of providing unprecedently accurate results for the static structure factor with the exception of the Wigner crystallization vicinity, despite the absence of adjustable or empirical parameters.

Assessing the accuracy of hybrid exchange-correlation functionals for the density response of warm dense electrons

We assess the accuracy of common hybrid exchange-correlation (XC) functionals (PBE0, PBE0-1/3, HSE06, HSE03, and B3LYP) within the Kohn-Sham density functional theory for the harmonically perturbed electron gas at parameters relevant for the challenging conditions of the warm dense matter. Generated by laser-induced compression and heating in the laboratory, the warm dense matter is a state of matter that also occurs in white dwarfs and planetary interiors. We consider both weak and strong degrees of density inhomogeneity induced by the external field at various wavenumbers. We perform an error analysis by comparing with the exact quantum Monte Carlo results. In the case of a weak perturbation, we report the static linear density response function and the static XC kernel at a metallic density for both the degenerate ground-state limit and for partial degeneracy at the electronic Fermi temperature. Overall, we observe an improvement in the density response when the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are used, compared with the previously reported results for the PBE, PBEsol, local-density approximation, and AM05 functionals; B3LYP, on the other hand, does not perform well for the considered system. Additionally, the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are more accurate for the density response properties than SCAN in the regime of partial degeneracy.

Non-empirical Mixing Coefficient for Hybrid XC Functionals from Analysis of the XC Kernel

We present an analysis of the static exchange-correlation (XC) kernel computed from hybrid functionals with a single mixing coefficient such as PBE0 and PBE0-1/3. We break down the hybrid XC kernels into the exchange and correlation parts using the Hartree-Fock functional, the exchange-only PBE, and the correlation-only PBE. This decomposition is combined with exact data for the static XC kernel of the uniform electron gas and an Airy gas model within a subsystem functional approach. This gives us a tool for the non-empirical choice of the mixing coefficient under ambient and extreme conditions. Our analysis provides physical insights into the effect of the variation of the mixing coefficient in hybrid functionals, which is of immense practical value. The presented approach is general and can be used for other types of functionals like screened hybrids.

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.

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.