The response of the uniform electron gas (UEG) to an external perturbation is of paramount importance for many applications. Recently, highly accurate results for the static density response function and the corresponding local field correction have been provided both for warm dense matter [J. Chem. Phys. 151, 194104 (2019)] and strongly coupled electron liquid [Phys. Rev. B 101, 045129 (2020)] conditions based on exact ab initio path integral Monte Carlo (PIMC) simulations. In the present work, we further complete our current description of the UEG by exploring the high energy density regime, which is relevant for, e.g., astrophysical applications and inertial confinement fusion experiments. To this end, we present extensive new PIMC results for the static density response in the range of 0.05 ≤ r s ≤ 0.5 and 0.85 ≤ θ ≤ 8. These data are subsequently used to benchmark the accuracy of the widely used random phase approximation and the dielectric theory by Singwi, Tosi, Land, and Sjölander (STLS). Moreover, we compare our results to configuration PIMC data where they are available and find perfect agreement with a relative accuracy of 0.001 − 0.01%. All PIMC data are available online.The Uniform Electron Gas in the High Density Regime 2
IntroductionThe uniform electron gas (UEG), also known as jellium or quantum one-component plasma, is defined as a system of Coulomb interacting electrons in a neutralizing homogeneous background and constitutes one of the most import model systems in physics and chemistry [1,2,3]. Having been introduced as a simple model description for conduction electrons in metals [4], the UEG exhibits a variety of interesting effects, such as Wigner crystallization [5,6,7] and a possible incipient excitonic mode which has been predicted to appear at low density [8,9,10,11].Moreover, the availability of highly accurate quantum Monte Carlo (QMC) data [12,13,14,15,16] at metallic densities and zero-temperature has sparked remarkable developments in many fields, most notably the hitherto unequaled success of density functional theory (DFT) regarding the description of real materials [17,18].Recently, the interest in matter under extreme conditions [19] has led to the need for analogous quantum Monte Carlo data at finite temperature. Of particular importance is the so-called warm dense matter (WDM) regime, which is defined by two characteristic parameters that are both of the order of one: i) the density parameter r s = r/a B (with r and a B being the average inter-particle distance and first Bohr radius) and ii) the reduced temperature θ = k B T /E F (with k B and E F being the Boltzmann constant and Fermi energy [1]), cf. Fig. 1 below. These conditions occur in astrophysical objects like giant planet interiors [20,21,22] and are expected to manifest on the pathway towards inertial confinement fusion [23]. Furthermore, WDM is nowadays routinely realized in large research facilities around the globe like the Linac Coherent Light Source (LCLS) [24], the National Ignition Facility (NIF) [25], and th...