No Arabic abstract
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.
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.
A scheme for analyzing Thomson scattering of x-rays by warm dense matter, based on the average-atom model, is developed. Emphasis is given to x-ray scattering by bound electrons. Contributions to the scattered x-ray spectrum from elastic scattering by electrons moving with the ions and from inelastic scattering by free and bound electrons are evaluated using parameters (chemical potential, average ionic charge, free electron density, bound and continuum wave functions, and occupation numbers) taken from the average-atom model. The resulting scheme provides a relatively simple diagnostic for use in connection with x-ray scattering measurements. Applications are given to dense hydrogen, beryllium, aluminum, titanium, and tin plasmas. At high momentum transfer, contributions from inelastic scattering by bound electrons are dominant features of the scattered x-ray spectrum for aluminum, titanium, and tin.
The influence of finite relaxation times on Thomson scattering from warm-dense plasmas is examined within the framework of the average-atom approximation. Presently most calculations use the collision-free Lindhard dielectric function to evaluate the free-electron contribution to the Thomson cross section. In this work, we use the Mermin dielectric function, which includes relaxation time explicitly. The relaxation time is evaluated by treating the average atom as an impurity in a uniform electron gas and depends critically on the transport cross section. The calculated relaxation rates agree well with values inferred from the Ziman formula for the static conductivity and also with rates inferred from a fit to the frequency-dependent conductivity. Transport cross sections determined by the phase-shift analysis in the average-atom potential are compared with those evaluated in the commonly used Born approximation. The Born approximation converges to the exact cross sections at high energies; however, differences that occur at low energies lead to corresponding differences in relaxation rates. The relative importance of including relaxation time when modeling X-ray Thomson scattering spectra is examined by comparing calculations of the free-electron dynamic structure function for Thomson scattering using Lindhard and Mermin dielectric functions. Applications are given to warm-dense Be plasmas, with temperatures ranging from 2 to 32~eV and densities ranging from 2 to 64~g/cc.
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 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.