  In a recent article T. Grover [Phys. Rev. Lett. 111, 130402 (2013)]
introduced a simple method to compute Renyi entanglement entropies in the realm
of the auxiliary field quantum Monte Carlo algorithm. Here, we further develop
this approach and provide a stabilization scheme to compute higher order Renyi
entropies and an extension to access the entanglement spectrum. The method is
tested on systems of correlated topological insulators.
