Entanglement entropy and topological order in resonating valence-bond quantum spin liquids

On the triangular and kagome lattices, short-ranged resonating valence-bond wave functions can be sampled without the sign problem using a recently developed Pfaffian Monte Carlo scheme. In this Rapid Communication, we study the Renyi entanglement entropy in these wave functions using a replica-trick method. Using various spatial bipartitions, including the Levin-Wen construction, our finite-size scaled Renyi entropy gives a topological contribution consistent with γ=ln(2), as expected for a gapped Z2 quantum spin liquid. We prove that the mutual statistics is consistent with the toric code anyon model and rule out any other quasiparticle statistics such as the double semion model.