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

Exact ground state Monte Carlo method for Bosons without importance sampling

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

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

Path integral ground state with a fourth-order propagator: Application to condensed helium

Squeezed and fragmented states of strongly interacting bosons in a double well. Part II: Quantum Monte Carlo simulations

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

On importance sampling and independent Metropolis-Hastings with an unbounded weight function

Efficient Quasiparticle Determination beyond the Diagonal Approximation via Random Compression

Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality

Monte Carlo study of two-dimensional Bose-Hubbard model

Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations

Path integral molecular dynamics simulations for Green's function in a system of identical bosons

Population size bias in diffusion Monte Carlo

Path Integral Monte Carlo Study of a Doubly-Dipolar Bose Gas

Quadratic Scaling Bosonic Path Integral Molecular Dynamics

Dynamical computation of the density of states and Bayes factors using nonequilibrium importance sampling

Path integral Monte Carlo in a discrete variable representation with Gibbs sampling: dipolar planar rotor chain

An efficient and exact noncommutative quantum Gibbs sampler

Quantum Monte Carlo and perturbative study of two-dimensional Bose-Fermi mixtures

Sampling strategies for expectation values within the Herman--Kluk approximation

Unbiasing time-dependent Variational Monte Carlo by projected quantum evolution