Entanglement spectra of interacting fermions in quantum Monte Carlo simulations

In a recent article, 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 develop this approach further 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.