Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements

Abstract

We propose a method to mitigate heavy-tailed distributions in fermion Quantum Monte Carlo simulations originating from zeros of the fermion determinant. We show that by merely changing the synchronization between local updates and computation of observables, one can reduce the prefactor of the heavy-tailed distribution, thus substantially suppressing statistical fluctuations of observables. The method is especially suitable for local observables similar to e.g. double occupancy, where the resulting speedup can reach two orders of magnitude. For observables, containing spatial correlators, the speedup is more moderate, but still ranges between five and ten. Our results are independent on the nature of the auxiliary field, discrete or continuous, and pave the way to improve measurement strategies for Hybrid Monte Carlo simulations.