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

Quantum Monte Carlo for Gauge Fields and Matter without the Fermion Determinant

Quantum Monte Carlo for gauge fields and matter without the fermion determinant

Meron-Cluster Algorithms for Quantum Link Models

Cluster-Algorithm-Amenable Models of Gauge Fields and Matter

Solving the Fermion Sign Problem in Quantum Monte Carlo Simulations by Majorana Representation

A variational Monte Carlo algorithm for lattice gauge theories with continuous gauge groups: a study of (2+1)-dimensional compact QED with dynamical fermions at finite density

Some bijective correspondences involving domino tableaux

Quantum simulations of quantum electrodynamics in Coulomb gauge

Monte Carlo Study of Lattice Compact Quantum Electrodynamics with Fermionic Matter: the Parent State of Quantum Phases

Fermion-qudit quantum processors for simulating lattice gauge theories with matter

Pfaffian quantum Monte Carlo: solution to Majorana sign ambiguity and applications

Quantum-assisted Monte Carlo algorithms for fermions

Resource-Efficient Quantum Simulation of Lattice Gauge Theories in Arbitrary Dimensions: Solving for Gauss' Law and Fermion Elimination

Computation of Dynamical Correlation Functions for Many Fermion Systems with Auxiliary-Field Quantum Monte Carlo

Sign-Problem-Free Fermionic Quantum Monte Carlo: Developments and Applications

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

General quantum algorithms for Hamiltonian simulation with applications to a non-Abelian lattice gauge theory

Quantum computation of thermal averages for a non-Abelian lattice gauge theory via quantum Metropolis sampling

Post-measurement Quantum Monte Carlo

Preempting Fermion Sign Problem: Unveiling Quantum Criticality through Nonequilibrium Dynamics

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