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 molecular dynamics simulations for Green's function in a system of identical bosons

Quadratic Scaling Bosonic Path Integral Molecular Dynamics

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

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 molecular dynamics for anyons, bosons and fermions

Green's Functions for Spin Boson Systems: Beyond Conventional Perturbation Theories

Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials

Approach to Solving Spin-Boson Dynamics Via Non-Markovian Quantum Trajectories

GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics

Ab initio path integral Monte Carlo simulations of warm dense two-component systems without fixed nodes: Structural properties

Ab initio simulations of the thermodynamic properties and phase transition of Fermi systems based on fictitious identical particles and physics-informed neural networks

Performance of wave function and Green's functions based methods for non equilibrium many-body dynamics

Many-Body Green's Function Theory for Electronic Excitations in Complex Chemical Systems

Pauli-Potential and Green Function Monte-Carlo Method for Many-Fermion Systems

Path integral molecular dynamics approximations of quantum canonical observables

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

An Advective-Spectral-Mixed Method for Time-Dependent Many-Body Wigner Simulations.

Non-equilibrium Green's function approach to inhomogeneous quantum many-body systems using the Generalized Kadanoff Baym Ansatz