  Entanglement entropy encodes important features of strongly interacting
quantum many-body systems and gauge theories, but its analytical study is still
limited to systems with high level of symmetry. This motivates the search for
efficient techniques to investigate this quantity numerically, through Monte
Carlo calculations on the lattice. In this contribution, we discuss the
computation of the entropic c-function by means of an algorithm based on
Jarzynski's equality, which is an exact theorem from non equilibrium
statistical mechanics. After presenting benchmark results for the Ising model
in two dimensions, where our algorithm successfully reproduces the analytical
predictions from conformal field theory, we discuss its generalization to the
three-dimensional Ising model, for which we were able to extract universal
terms beyond the area law. Finally we point out some future generalizations of
this calculation.
