Taylor series perspective on ab initio path integral Monte Carlo simulations with Fermi-Dirac statistics

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Dornheim, Tobias, Robles, Alexander Benedix, Hamann, Paul, Chuna, Thomas M., Svensson, Pontus, Schwalbe, Sebastian, Moldabekov, Zhandos A., Tolias, Panagiotis, Vorberger, Jan
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866912586046373888
author Dornheim, Tobias
Robles, Alexander Benedix
Hamann, Paul
Chuna, Thomas M.
Svensson, Pontus
Schwalbe, Sebastian
Moldabekov, Zhandos A.
Tolias, Panagiotis
Vorberger, Jan
author_facet Dornheim, Tobias
Robles, Alexander Benedix
Hamann, Paul
Chuna, Thomas M.
Svensson, Pontus
Schwalbe, Sebastian
Moldabekov, Zhandos A.
Tolias, Panagiotis
Vorberger, Jan
contents The fermion sign problem constitutes a fundamental computational bottleneck across a plethora of research fields in physics, quantum chemistry and related disciplines. Recently, it has been suggested to alleviate the sign problem in \emph{ab initio} path integral Molecular Dynamics and path integral Monte Carlo (PIMC) calculations based on the simulation of fictitious identical particles that are represented by a continuous quantum statistics variable $ξ$ [\textit{J.~Chem.~Phys.}~\textbf{157}, 094112 (2022)]. This idea facilitated a host of applications including the interpretation of an x-ray scattering experiment with strongly compressed beryllium at the National Ignition Facility [\textit{Nature Commun.}~\textbf{16}, 5103 (2025)]. In the present work, we express the original isothermal $ξ$-extrapolation method as a special case of a truncated Taylor series expansion around the $ξ=0$ limit of distinguishable particles. We derive new PIMC estimators that allow us to evaluate the Taylor coefficients up to arbitrary order and we carry out extensive new PIMC simulations of the warm dense electron gas to systematically analyze the sign problem from this new perspective. This gives us important insights into the applicability of the $ξ$-extrapolation method for different levels of quantum degeneracy in terms of the Taylor series radius of convergence. Moreover, the direct PIMC evaluation of the $ξ$-derivatives, in principle, removes the necessity for simulations at different values of $ξ$ and can facilitate more efficient simulations that are designed to maximize compute time in those regions of the full permutation space that contribute most to the final Taylor estimate of the fermionic expectation value of interest.
format Preprint
id arxiv_https___arxiv_org_abs_2509_11317
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Taylor series perspective on ab initio path integral Monte Carlo simulations with Fermi-Dirac statistics
Dornheim, Tobias
Robles, Alexander Benedix
Hamann, Paul
Chuna, Thomas M.
Svensson, Pontus
Schwalbe, Sebastian
Moldabekov, Zhandos A.
Tolias, Panagiotis
Vorberger, Jan
Chemical Physics
Plasma Physics
The fermion sign problem constitutes a fundamental computational bottleneck across a plethora of research fields in physics, quantum chemistry and related disciplines. Recently, it has been suggested to alleviate the sign problem in \emph{ab initio} path integral Molecular Dynamics and path integral Monte Carlo (PIMC) calculations based on the simulation of fictitious identical particles that are represented by a continuous quantum statistics variable $ξ$ [\textit{J.~Chem.~Phys.}~\textbf{157}, 094112 (2022)]. This idea facilitated a host of applications including the interpretation of an x-ray scattering experiment with strongly compressed beryllium at the National Ignition Facility [\textit{Nature Commun.}~\textbf{16}, 5103 (2025)]. In the present work, we express the original isothermal $ξ$-extrapolation method as a special case of a truncated Taylor series expansion around the $ξ=0$ limit of distinguishable particles. We derive new PIMC estimators that allow us to evaluate the Taylor coefficients up to arbitrary order and we carry out extensive new PIMC simulations of the warm dense electron gas to systematically analyze the sign problem from this new perspective. This gives us important insights into the applicability of the $ξ$-extrapolation method for different levels of quantum degeneracy in terms of the Taylor series radius of convergence. Moreover, the direct PIMC evaluation of the $ξ$-derivatives, in principle, removes the necessity for simulations at different values of $ξ$ and can facilitate more efficient simulations that are designed to maximize compute time in those regions of the full permutation space that contribute most to the final Taylor estimate of the fermionic expectation value of interest.
title Taylor series perspective on ab initio path integral Monte Carlo simulations with Fermi-Dirac statistics
topic Chemical Physics
Plasma Physics
url https://arxiv.org/abs/2509.11317