Entanglement entropy from non-equilibrium Monte Carlo simulations

We study the entanglement entropy in lattice field theory using a simulation algorithm based on Jarzynski's theorem. We focus on the entropic c-function for the Ising model in two and in three dimensions: after validating our algorithm against known analytical results from conformal field theory in two dimensions, we present novel results for the three-dimensional case. We show that our algorithm, which is highly parallelized on graphics processing units, allows one to precisely determine the subleading corrections to the area law, which have been investigated in many recent works. Possible generalizations of this study to other strongly coupled theories are discussed.