No Arabic abstract
We present an emph{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,omega)$, the static structure factor $S(q)$, and the interaction energy $v$. The ESA combines the recent neural-net representation [textit{J. Chem. Phys.} textbf{151}, 194104 (2019)] 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 re-evaluating the results of the recent {X-ray Thomson scattering experiment on aluminum} by Sperling textit{et al.}~[textit{Phys. Rev. Lett.} textbf{115}, 115001 (2015)]. 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.
The energy band structures caused by self-energy shifting that results in bound energy levels broadening and merging in warm dense aluminum and beryllium are observed. An energy band theory for warm dense matter (WDM) is proposed and a new code based on the energy band theory is developed to improve the traditional density functional method. Massive data of the equation of state and transport coefficients for WDM in medium and low Z have been simulated. The transition from fully degenerate to partially degenerate (related to WDM) and finally to non-degenerate state is investigated using the Lorenz number varying with the degeneracy parameter, and the lower and upper parameter boundaries for WDM are achieved. It is shown that the pressure ionization results in the Wiedemann-Franz law no longer available for WDM.
In a recent Letter [T.~Dornheim emph{et al.}, Phys.~Rev.~Lett.~textbf{125}, 085001 (2020)], we have presented the first emph{ab initio} results for the nonlinear density response of electrons in the warm dense matter regime. In the present work, we extend these efforts by carrying out extensive new path integral Monte Carlo (PIMC) simulations of a emph{ferromagnetic} electron gas that is subject to an external harmonic perturbation. This allows us to unambiguously quantify the impact of spin-effects on the nonlinear density response of the warm dense electron gas. In addition to their utility for the description of warm dense matter in an external magnetic field, our results further advance our current understanding of the uniform electron gas as a fundamental model system, which is important in its own right.
Understanding many processes, e.g. fusion experiments, planetary interiors and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter (WDM) and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly-free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive emph{ab initio} atomistic computation, typically based on quantum mechanical Kohn-Sham Density Functional Theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the WDM regime. A recently-developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (<50 eV). We have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.
We present a model for electron-ion transport in Warm Dense Matter that incorporates Coulomb coupling effects into the quantum Boltzmann equation of Uehling and Uhlenbeck through the use of a statistical potential of mean force. Although this model has been derived rigorously in the classical limit [S.D. Baalrud and J. Daligault, Physics of Plasmas 26, 8, 082106 (2019)], its quantum generalization is complicated by the uncertainty principle. Here we apply an existing model for the potential of mean force based on the quantum Ornstein-Zernike equation coupled with an average-atom model [C. E. Starrett, High Energy Density Phys. 25, 8 (2017)]. This potential contains correlations due to both Coulomb coupling and exchange, and the collision kernel of the kinetic theory enforces Pauli blocking while allowing for electron diffraction and large-angle collisions. By solving the Uehling-Uhlenbeck equation for electron-ion relaxation rates, we predict the momentum and temperature relaxation time and electrical conductivity of solid density aluminum plasma based on electron-ion collisions. We present results for density and temperature conditions that span the transition from classical weakly-coupled plasma to degenerate moderately-coupled plasma. Our findings agree well with recent quantum molecular dynamics simulations.
Warm dense matter (WDM) -- an exotic state of highly compressed matter -- has attracted high interest in recent years in astrophysics and for dense laboratory systems. At the same time, this state is extremely difficult to treat theoretically. This is due to the simultaneous appearance of quantum degeneracy, Coulomb correlations and thermal effects, as well as the overlap of plasma and condensed phases. Recent breakthroughs are due to the successful application of density functional theory (DFT) methods which, however, often lack the necessary accuracy and predictive capability for WDM applications. The situation has changed with the availability of the first textit{ab initio} data for the exchange-correlation free energy of the warm dense uniform electron gas (UEG) that were obtained by quantum Monte Carlo (QMC) simulations, for recent reviews, see Dornheim textit{et al.}, Phys. Plasmas textbf{24}, 056303 (2017) and Phys. Rep. textbf{744}, 1-86 (2018). In the present article we review recent further progress in QMC simulations of the warm dense UEG: namely, textit{ab initio} results for the static local field correction $G(q)$ and for the dynamic structure factor $S(q,omega)$. These data are of key relevance for the comparison with x-ray scattering experiments at free electron laser facilities and for the improvement of theoretical models. In the second part of this paper we discuss simulations of WDM out of equilibrium. The theoretical approaches include Born-Oppenheimer molecular dynamics, quantum kinetic theory, time-dependent DFT and hydrodynamics. Here we analyze strengths and limitations of these methods and argue that progress in WDM simulations will require a suitable combination of all methods. A particular role might be played by quantum hydrodynamics, and we concentrate on problems, recent progress, and possible improvements of this method.