Abstract:The \emph{ab initio} path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion sign problem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the sign problem by carrying out PIMC simulations of fictitious systems which contain an interpolating continuous variable $\xi\in[-1,1]$ in their partition function, with the physical Fermi- and Bose-statistics corresponding to the endpoint limits $\xi=-1$ and $\xi=1$. It has been proposed that thermodynamic information about the fermionic limit might be obtained by path integral calculations within the bosonic sector $\xi>0$ combined with a quadratic $\xi$ extrapolation throughout the fermionic sector $\xi<0$, essentially bypassing the sign problem. In this work, we show how the inclusion of the artificial parameter $\xi$ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We empirically demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the archetypal uniform electron gas (UEG), where we find excellent agreement ($\sim0.5\%$) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding eleven orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG).

Application of Deconvolution in Path Integral Simulations

Application of Deconvolution in Path Integral Simulations

Path integral Lindblad master equation through transfer tensor method & the generalized quantum master equation

A path integral methodology for obtaining thermodynamic properties of nonadiabatic systems using Gaussian mixture distributions

Path-integral methodology and simulations of quantum thermal transport: Full counting statistics approach

Hybrid Quantum Algorithm for Simulating Real-Time Thermal Correlation Functions

Padé and Padé-Laplace Methods for masses and matrix elements

Certain aspects of prestack deconvolution

Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Numerical implementation of some reweighted path integral methods

On the Practical Applications of Information Field Dynamics

Path Integral Methods and Applications

Fermionic physics from ab initio path integral Monte Carlo simulations of fictitious identical particles

Ab Initio Path Integral Monte Carlo Simulations of the Uniform Electron Gas on Large Length Scales

Gaussian deconvolution on $\mathbb R^d$ with application to self-repellent Brownian motion

Construction of stochastic hybrid path integrals using "quantum-mechanical'' operators

Nonperturbative Time-Convolutionless Quantum Master Equation From The Path Integral Approach

Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water

Introduction to Configuration Path Integral Monte Carlo

Classical-path integral adaptive resolution in molecular simulation: towards a smooth quantum-classical coupling

Grassmann time-evolving matrix product operators: An efficient numerical approach for fermionic path integral simulations