Explicit treatment of many-body Fermi statistics in path integral Monte Carlo (PIMC) results in exponentially scaling computational cost due to the near cancellation of contributions to observables from even and odd permutations. Through direct analysis of exchange statistics we find that individual exchange probabilities in homogeneous systems are, except for finite size effects, independent of the configuration of other permutations present. For two representative systems, 3-He and the homogeneous electron gas, we show that this allows the entire antisymmetrized density matrix to be generated from a simple model depending on only a few parameters obtainable directly from a standard PIMC simulation. The result is a polynomial scaling algorithm and up to a 10 order of magnitude increase in efficiency in measuring fermionic observables for the systems considered.
Explicit treatment of many-body Fermi statistics in path integral Monte Carlo (PIMC) results in exponentially scaling computational cost due to the near cancellation of contributions to observables from even and odd permutations. Through direct analysis of exchange statistics we find that individual exchange probabilities in homogeneous systems are, except for finite size effects, independent of the configuration of other permutations present. For two representative systems, 3-He and the homogeneous electron gas, we show that this allows the entire antisymmetrized density matrix to be generated from a simple model depending on only a few parameters obtainable directly from a standard PIMC simulation. The result is a polynomial scaling algorithm and up to a 10 order of magnitude increase in efficiency in measuring fermionic observables for the systems considered.
Being motivated by the surge of fermionic quantum Monte Carlo simulations at finite temperature, we present a detailed analysis of the permutation-cycle properties of path integral Monte Carlo (PIMC) simulations of degenerate electrons. Particular emphasis is put onto the uniform electron gas in the warm dense matter regime. We carry out PIMC simulations of up to N = 100 electrons and investigate exchange-cycle frequencies, which are found not to follow any simple exponential law even in the case of ideal fermions due to the finite size of the simulation box. Moreover, we introduce a permutation-cycle correlation function, which allows us to analyze the joint probability to simultaneously find cycles of different lengths within a single configuration. Again, we find that finite-size effects predominate the observed behavior. Finally, we briefly consider an inhomogeneous system, namely, electrons in a 2D harmonic trap. We expect our results to be of interest for the further development of fermionic PIMC methods, in particular, to alleviate the notorious fermion sign problem.
The 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 simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the sign problem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We 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 uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 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).
Correlated fermions are of high interest in condensed matter (Fermi liquids, Wigner molecules), cold atomic gases and dense plasmas. Here we propose a novel approach to path integral Monte Carlo (PIMC) simulations of strongly degenerate non-ideal fermions at finite temperature by combining a fourth-order factorization of the density matrix with antisymmetric propagators, i.e., determinants, between all imaginary time slices. To efficiently run through the modified configuration space, we introduce a modification of the widely used continuous space worm algorithm, which allows for an efficient sampling at arbitrary system parameters. We demonstrate how the application of determinants achieves an effective blocking of permutations with opposite signs, leading to a significant relieve of the fermion sign problem. To benchmark the capability of our method regarding the simulation of degenerate fermions, we consider multiple electrons in a quantum dot and compare our results with other ab initio techniques, where they are available. The present permutation blocking PIMC approach allows us to obtain accurate results even for N = 20 electrons at low temperature and arbitrary coupling, where no other ab initio results have been reported, so far.
It is widely known that there is no sign problem in path integral Monte Carlo (PIMC) simulations of fermions in one dimension. As far as the author is aware, there is no direct proof of this in the literature. This work shows that the sign of the N-fermion antisymmetric free propagator is given by the product of all possible pairs of particle separations, or relative displacements. For a nonvanishing closed-loop product of such propagators, as required by PIMC, all relative displacements from adjacent propagators are paired into perfect squares, and therefore the loop product must be positive, but only in one dimension. By comparison, permutation sampling, which does not evaluate the determinant of the antisymmetric propagator exactly, remains plagued by a low-level sign problem, even in one dimension.
ABSTRACT Ab initio path integral Monte Carlo (PIMC) simulations constitute the gold standard for the estimation of a broad range of equilibrium properties of a host of interacting quantum many‐body systems spanning a broad range of conditions from ultracold atoms to warm dense quantum plasmas. A key practical limitation is given by the notorious fermion sign problem, which manifests as an exponential computational bottleneck with respect to system size and inverse temperature. In practice, the sign problem is particularly severe in the grand‐canonical ensemble, where the bosonic and fermionic configuration spaces differ not only with respect to the symmetry of the thermal density matrix but, crucially, also with respect to the particle number distribution for a given chemical potential μ [T. Dornheim, J. Phys. A 54 , 335001 (2021)]. Here, we present a simple reweighting scheme that basically allows one to retain access to grand‐canonical expectation values at the cost of fermionic PIMC simulations in the canonical ensemble for the largest significant particle number in the fermionic sector. As a practical example, we consider the warm dense electron gas, which has attracted considerable recent attention due to its relevance, e.g., for the modeling of compact astrophysical objects and inertial fusion energy applications.
Recent advances in warm dense matter physics, e.g. laser compressed matter, lead to an increasing interest in the description of correlated, degenerate electrons at finite temperatures. Path integral Monte Carlo (PIMC) methods cannot correctly describe weakly to moderately coupled and strongly degenerate Fermi systems due to the so-called fermion sign problem. The Configuration Path Integral Monte Carlo (CPIMC) approach greatly reduces the sign problem and allows for the exact computation of thermodynamic properties in this regime. In addition, the first successful implementation of the diagrammatic worm algorithm for a general Hamiltonian in Fock space with arbitrary pair interactions gives direct access to the Matsubara Green function.
This thesis demonstrates the capabilities of the CPIMC approach for a model system of Coulomb interacting fermions in a two-dimensional harmonic trap. The correctness of the CPIMC implementation is verified by rigorous comparisons with an exact diagonalization method. Benchmark results are presented, which reveal large errors of the Hartree-Fock approximation in open shell configurations even for weak coupling and a significant deviation of multi-level blocking PIMC data in the complete basis set limit.
The application of the CPIMC method to the warm dense homogeneous electron gas (HEG) quantifies the accuracy of recently published restricted PIMC (RPIMC) results, which have been the basis for the construction of exchange-correlation free energy functionals to be used in finite-temperature density functional theory calculations of warm dense matter. It is shown that the errors of the RPIMC data exceed 10 % at intermediate densities. Additionally, highly accurate data for the exchange-correlation energy at high densities, which are inaccessible by the RPIMC method, are provided. These results are needed to significantly increase the quality of future exchange-correlation functionals to be used in finite-temperature applications.