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.
Hamann et al. (Tue,) studied this question.