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

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Dornheim, Tobias, Chuna, Thomas M., Bellenbaum, Hannah M., Moldabekov, Zhandos, Tolias, Panagiotis, Vorberger, Jan
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866916669696245760
author Dornheim, Tobias
Chuna, Thomas M.
Bellenbaum, Hannah M.
Moldabekov, Zhandos
Tolias, Panagiotis
Vorberger, Jan
author_facet Dornheim, Tobias
Chuna, Thomas M.
Bellenbaum, Hannah M.
Moldabekov, Zhandos
Tolias, Panagiotis
Vorberger, Jan
contents 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 $χ(\mathbf{q})$ and the static quadratic density response $χ^{(2)}(\mathbf{q},0)$ fluctuate for small $q$. In addition, we perform an analytic continuation to compute the dynamic structure factor $S(\mathbf{q},ω)$ from PIMC results of the imaginary-time density--density correlation function $F(\mathbf{q},τ)$ 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.
format Preprint
id arxiv_https___arxiv_org_abs_2504_00737
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Application of a spherically averaged pair potential in \emph{ab initio} path integral Monte Carlo simulations of the warm dense electron gas
Dornheim, Tobias
Chuna, Thomas M.
Bellenbaum, Hannah M.
Moldabekov, Zhandos
Tolias, Panagiotis
Vorberger, Jan
Plasma Physics
Chemical Physics
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 $χ(\mathbf{q})$ and the static quadratic density response $χ^{(2)}(\mathbf{q},0)$ fluctuate for small $q$. In addition, we perform an analytic continuation to compute the dynamic structure factor $S(\mathbf{q},ω)$ from PIMC results of the imaginary-time density--density correlation function $F(\mathbf{q},τ)$ 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.
title Application of a spherically averaged pair potential in \emph{ab initio} path integral Monte Carlo simulations of the warm dense electron gas
topic Plasma Physics
Chemical Physics
url https://arxiv.org/abs/2504.00737