First-Principles Estimation of Electronic Temperature from X-Ray Thomson Scattering Spectrum of Isochorically Heated Warm Dense Matter

Extended Time-Dependent Density Functional Theory for Multibody Densities

Time-dependent density functional theory (TDDFT) is widely used for understanding and predicting properties and behaviors of matter. As one of the fundamental theorems in TDDFT, Van Leeuwen theorem [Phys. Rev. Lett. 82, 3863 (1999)PRLTAO0031-900710.1103/PhysRevLett.82.3863] guarantees how to construct a unique potential with the same one-body density evolution. Here we extend Van Leeuwen theorem by exploring truncation criteria in Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. Our generalized theorem demonstrates the existence of a unique nonlocal potential to accurately reconstruct the multibody density evolution in binary interacting systems. Under nonstringent conditions, truncation of the BBGKY-hierarchy equations aligns with the behavior of multibody density evolution and maintains consistency in the reduced equations. Based on this theorem, we further construct a static many-body potential and extend the well-known Casida equation for multiple-excitation energy. The extended TDDFT provides an accurate and first-principles framework capable of describing the kinetic processes of correlated systems, including strongly coupled particle transport, multiple-excitation, and ionization processes.

Ultrafast Heating-Induced Suppression of d-Band Dominance in the Electronic Excitation Spectrum of Cuprum

The combination of isochoric heating of solids by free-electron lasers (FELs) and in situ diagnostics by X-ray Thomson scattering (XRTS) allows for measurements of material properties at warm dense matter (WDM) conditions relevant for astrophysics, inertial confinement fusion, and materials science. In the case of metals, the FEL beam pumps energy directly into electrons with the lattice structure of ions being nearly unaffected. This leads to a unique transient state that gives rise to a set of interesting physical effects, which can serve as a reliable testing platform for WDM theories. In this work, we present extensive linear-response time-dependent density functional theory (TDDFT) results for the electronic dynamic structure factor of isochorically heated copper with a face-centered cubic lattice. At ambient conditions, the plasmon is heavily damped due to the presence of d-band excitations, and its position is independent of the wavenumber. In contrast, the plasmon feature starts to dominate the excitation spectrum and has a Bohm-Gross-type plasmon dispersion for temperatures T ≥ 4 eV, where the quasi-free electrons in the interstitial region are in the WDM regime. In addition, we analyze the thermal changes in the d-band excitations and outline the possibility to use future XRTS measurements of isochorically heated copper as a controlled testbed for WDM theories.

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula

We study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

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.

Accurate temperature diagnostics for matter under extreme conditions

The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.

Thermal excitation signals in the inhomogeneous warm dense electron gas

We investigate the emergence of electronic excitations from the inhomogeneous electronic structure at warm dense matter parameters based on first-principles calculations. The emerging modes are controlled by the imposed perturbation amplitude. They include satellite signals around the standard plasmon feature, transformation of plasmons to optical modes, and double-plasmon modes. These modes exhibit a pronounced dependence on the temperature. This makes them potentially invaluable for the diagnostics of plasma parameters in the warm dense matter regime. We demonstrate that these modes can be probed with present experimental techniques.

Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization

We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.

Bound electron screening effect on ion-ion potential of warm and hot dense matter

The effects of bound electron screening in warm and hot dense matter are investigated analytically and a theoretical description of screened short-range repulsion is given meanwhile. An empirical ion-ion potential including the classic charge screening and chemical bond attraction at various temperatures and densities is proposed. By solving hypernetted chain equations and comparing the obtained radial distribution function (RDF) with ab initio simulations, the proposed ion-ion potential is found to be promising over a wide range of temperatures and densities for warm dense aluminum and iron. The elastic scattering amplitude and the x-ray absorption near the edge structure of warm dense aluminum calculated from the obtained RDF are in good agreement with experiment results.

Finite-temperature density-functional-theory investigation on the nonequilibrium transient warm-dense-matter state created by laser excitation

We present a finite-temperature density-functional-theory investigation of the nonequilibrium transient electronic structure of warm dense Li, Al, Cu, and Au created by laser excitation. Photons excite electrons either from the inner shell orbitals or from the valence bands according to the photon energy, and give rise to isochoric heating of the sample. Localized states related to the 3d orbital are observed for Cu when the hole lies in the inner shell 3s orbital. The electrical conductivity for these materials at nonequilibrium states is calculated using the Kubo-Greenwood formula. The change of the electrical conductivity, compared to the equilibrium state, is different for the case of holes in inner shell orbitals or the valence band. This is attributed to the competition of two factors: the shift of the orbital energies due to reduced screening of core electrons, and the increase of chemical potential due to the excitation of electrons. The finite-temperature effect of both the electrons and the ions on the electrical conductivity is discussed in detail. This work is helpful to better understand the physics of laser excitation experiments of warm dense matter.

Effective Static Approximation: A Fast and Reliable Tool for Warm-Dense Matter Theory

We present an 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,ω), the static structure factor S(q), and the interaction energy v. The ESA combines the recent neural-net representation by T. Dornheim et al., [J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013] 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 reevaluating the results of the recent x-ray Thomson scattering experiment on aluminum by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.115001]. 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.

Ab initio dielectric response function of diamond and other relevant high pressure phases of carbon

The electronic structure and dielectric properties of the diamond, body centered cubic diamond (bc8), and hexagonal diamond (lonsdaleite) phases of carbon are computed using density functional theory and many-body perturbation theory with the emphasis on the excitonic picture of the solid phases relevant in the regimes of high-pressure physics and warm dense matter. We also discuss the capabilities of reproducing the inelastic x-ray scattering spectra in comparison with the existing models in light of recent x-ray scattering experiments on carbon and carbon bearing materials in the Megabar range.