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

Introduction to Configuration Path Integral Monte Carlo

Configuration Path Integral Monte Carlo Approach to the Static Density Response of the Warm Dense Electron Gas

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

\emph{Ab initio} Quantum Monte Carlo simulation of the warm dense electron gas

Monte Carlo Collision method for low temperature plasma simulation

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

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

Ab-initio study of interacting fermions at finite temperature with neural canonical transformation

Coupling Fluid Plasma and Kinetic Neutral Models using Correlated Monte Carlo Methods

Particle-in-cell Simulation Method for Macroscopic Degenerate Plasmas.

Precision Thermodynamics of the Fermi polaron at strong coupling

Path-integral simulation of solids

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

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

Stability of a Nondegenerate Two--Component Weakly Coupled Plasma

Strong coupling quantum impurity solver on the real and imaginary axis

A brief introduction to the diffusion Monte Carlo method and the fixed-node approximation

Dynamic exchange-correlation effects in the strongly coupled electron liquid

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