No Arabic abstract
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.
Average atom (AA) models allow one to efficiently compute electronic and optical properties of materials over a wide range of conditions and are often employed to interpret experimental data. However, at high pressure, predictions from AA models have been shown to disagree with results from ab initio computer simulations. Here we reconcile these deviations by developing an innovative type of AA model, AVION, that computes the electronic eigenstates with novel boundary conditions within the ion sphere. Bound and free states are derived consistently. We drop the common AA image that the free-particle spectrum starts at the potential threshold, which we found to be incompatible with ab initio calculations. We perform ab initio simulations of crystalline and liquid carbon and aluminum over a wide range of densities and show that the computed band structure is in very good agreement with predictions from AVION.
The free-free opacity in plasmas is fundamental to our understanding of energy transport in stellar interiors and for inertial confinement fusion research. However, theoretical predictions in the challenging dense plasma regime are conflicting and there is a dearth of accurate experimental data to allow for direct model validation. Here we present time-resolved transmission measurements in solid-density Al heated by an XUV free-electron laser. We use a novel functional optimization approach to extract the temperature-dependent absorption coefficient directly from an oversampled pool of single-shot measurements, and find a pronounced enhancement of the opacity as the plasma is heated to temperatures of order the Fermi energy. Plasma heating and opacity-enhancement is observed on ultrafast time scales, within the duration of the femtosecond XUV pulse. We attribute further rises in the opacity on ps timescales to melt and the formation of warm-dense matter.
Wave packet molecular dynamics (WPMD) has recently received a lot of attention as a computationally fast tool to study dynamical processes in warm dense matter beyond the Born-Oppenheimer approximation. These techniques, typically, employ many approximations to achieve computational efficiency while implementing semi-empirical scaling parameters to retain accuracy. We investigate three of the main approximations ubiquitous to WPMD: a restricted basis set, approximations to exchange, and the lack of correlation. We examine each of these approximations in atomic and molecular hydrogen in addition to a dense hydrogen plasma. We find that the biggest improvement to WPMD comes from combining a two Gaussian basis with a semi-empirical correction based on the valence-bond wave function. A single parameter scales this correction to match experimental pressures of dense hydrogen. Ultimately, we find that semi-empirical scaling parameters are necessary to correct for the main approximations in WPMD. However, reducing the scaling parameters for more ab-initio terms gives more accurate results and displays the underlying physics more readily.
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(mathbf{q},omega)$ constitutes a key quantity as it is directly measured in X-ray Thomson (XRTS) scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first textit{ab initio} results for $S(mathbf{q},omega)$ 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(mathbf{q},omega)$. In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for $S(mathbf{q},omega)$ over substantial parts of the warm dense matter regime.
In a classical plasma the momentum distribution, $n(k)$, decays exponentially, for large $k$, and the same is observed for an ideal Fermi gas. However, when quantum and correlation effects are relevant simultaneously, an algebraic decay, $n_infty(k)sim k^{-8}$ has been predicted. This is of relevance for cross sections and threshold processes in dense plasmas that depend on the number of energetic particles. Here we present extensive textit{ab initio} results for the momentum distribution of the nonideal uniform electron gas at warm dense matter conditions. Our results are based on first principle fermionic path integral Monte Carlo (CPIMC) simulations and clearly confirm the $k^{-8}$ asymptotic. This asymptotic behavior is directly linked to short-range correlations which are analyzed via the on-top pair distribution function (on-top PDF), i.e. the PDF of electrons with opposite spin. We present extensive results for the density and temperature dependence of the on-top PDF and for the momentum distribution in the entire momentum range.