Application of a spherically averaged pair potential in \emph{ab initio} path integral Monte Carlo simulations of the warm dense electron gas

Abstract

Spherically averaged periodic pair potentials offer the enticing promise to provide accurate results at a drastically reduced computational cost compared to the traditional Ewald sum. In this work, we employ the pair potential by Yakub and Ronchi [\textit{J.~Chem.~Phys.}~\textbf{119}, 11556 (2003)] in \emph{ab initio} path integral Monte Carlo (PIMC) simulations of the warm dense uniform electron gas. Overall, we find very accurate results with respect to Ewald reference data for integrated properties such as the kinetic and potential energy, whereas wavenumber resolved properties such as the static structure factor $S(\mathbf{q})$, the static linear density response $\chi(\mathbf{q})$ and the static quadratic density response $\chi^{(2)}(\mathbf{q},0)$ fluctuate for small $q$. In addition, we perform an analytic continuation to compute the dynamic structure factor $S(\mathbf{q},\omega)$ from PIMC results of the imaginary-time density--density correlation function $F(\mathbf{q},\tau)$ for both pair potentials. Our results have important implications for future PIMC calculations, which can be sped up significantly using the YR potential for the estimation of equation-of-state properties or $q$-resolved observables in the non-collective regime, whereas a full Ewald treatment is mandatory to accurately resolve physical effects manifesting for smaller $q$, including the evaluation of compressibility sum rules, the interpretation of x-ray scattering experiments at small scattering angles, and the estimation of optical and transport properties.