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

Implicit Purification for Temperature-Dependent Density Matrices

Positivity Preserving Density Matrix Minimization at Finite Temperatures via Square Root

Finite-Temperature Evaluation of the Fermi Density Operator

Matrix product purifications for canonical ensembles and quantum number distributions

Accelerated density matrix expansions for Born-Oppenheimer molecular dynamics

Saturation of Thermal Complexity of Purification

Self-consistent Density Functional Calculations without Diagonalization

Grassmann Time-Evolving Matrix Product Operators for Equilibrium Quantum Impurity Problems

Efficient temperature-dependent Green's functions methods for realistic systems: compact grids for orthogonal polynomial transforms

Real-time Impurity Solver Using Grassmann Time-Evolving Matrix Product Operators

Preserving the Hermiticity of the One-Body Density Matrix for a Non-Interacting Fermi Gas

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

Integrating by parts at finite density

Frozen Gaussian Series Representation of the Imaginary Time Propagator Theory and Numerical Tests

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

Efficient Simulation of 1D Long-Range Interacting Systems at Any Temperature

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

Infinite Grassmann time-evolving matrix product operators for non-equilibrium quantum impurity problems