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 for Exact Quantum Statistics of Multi-Electronic-State Systems

Efficient “Middle” Thermostat Scheme for the Quantum/Classical Canonical Ensemble via Molecular Dynamics

Quasiclassical Trajectory Simulation as a Protocol to Build Locally Accurate Machine Learning Potentials

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Different Flavors of Exact-Factorization-Based Mixed Quantum-Classical Methods for Multistate Dynamics

Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics

A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations

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

Quantum Circuit Matrix Product State Ansatz for Large-Scale Simulations of Molecules

Calculating potential energy surfaces with quantum computers by measuring only the density along adiabatic transitions

Accurate Calculation of Equilibrium Reduced Density Matrix for the System-Bath Model: A Multilayer Multiconfiguration Time-Dependent Hartree Approach and Its Comparison to a Multi-Electronic-state Path Integral Molecular Dynamics Approach

Potential energy landscape formalism for quantum molecular liquids

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

Path integral molecular dynamics approximations of quantum canonical observables

Elective-momentum Ultrasize Quantum Monte Carlo Method

Ab initio path integral Monte Carlo simulations of the uniform electron gas on large length scales

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

A Quantum-Mechanics Molecular-Mechanics scheme for extended systems

Dimension-free Ergodicity of Path Integral Molecular Dynamics

Heat transfer statistics in mixed quantum-classical systems