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).

Path integral Monte Carlo ground state approach: formalism, implementation, and applications

Path Integral Approach for Multi-Dimensional Polarons in A Symmetric Quantum Dot

Properties of fermionic systems with the path-integral ground state method

Exact ground state Monte Carlo method for Bosons without importance sampling

Introduction to Configuration Path Integral Monte Carlo

Numerical and Monte Carlo Bethe Ansatz Method: 1d Heisenberg Model

Monte Carlo Based Techniques for Quantum Magnets with Long-Range Interactions

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

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

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

Path Integral Monte Carlo Simulations for Fermion Systems: Pairing in the Electron-Hole Plasma

Path-integral approaches to strongly-coupled quantum many-body systems

Applications of the Monte Carlo Method in Statistical Physics

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

Determinant and Derivative-Free Quantum Monte Carlo Within the Stochastic Representation of Wavefunctions

Determinant- and Derivative-Free Quantum Monte Carlo Within the Stochastic Representation of Wavefunctions

Path-integral simulation of solids

Wasserstein Quantum Monte Carlo: A Novel Approach for Solving the Quantum Many-Body Schrödinger Equation