Entanglement, noise, and the cumulant expansion

We put forward a variation of a recently proposed method to overcome the signal-to-noise problem found in Monte Carlo calculations of the entanglement entropy of interacting fermions. The present method takes advantage of the approximate lognormal distributions that characterize the signal-to-noise properties of other approaches. We test our technique by implementing it in combination with the hybrid Monte Carlo algorithm and calculating the second R\'enyi entanglement entropy of the 1D Hubbard model.