  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.
