Reweighting scheme for the calculation of grand-canonical expectation values in quantum Monte Carlo simulations with a fermion sign problem

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 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 grandcanonical 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 grandcanonical 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