Entanglement spectroscopy using quantum Monte Carlo

We present a numerical scheme to reconstruct a subset of the entanglement spectrum of quantum many body systems using quantum Monte Carlo. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector and typically performs better when the eigenvalues are not too different.