Topological and Entanglement Properties of Resonating Valence Bond wavefunctions

We examine in details the connections between topological and entanglement properties of short-range resonating valence bond (RVB) wave functions using projected entangled pair states (PEPS) on kagome and square lattices on (quasi)infinite cylinders with generalized boundary conditions (and perimeters with up to 20 lattice spacings). By making use of disconnected topological sectors in the space of dimer lattice coverings, we explicitly derive (orthogonal) “minimally entangled” PEPS RVB states. For the kagome lattice, using the quantum Heisenberg antiferromagnet as a reference model, we obtain the finite-size scaling with increasing cylinder perimeter of the vanishing energy separations between these states. In particular, we extract two separate (vanishing) energy scales corresponding (i) to insert a vison line between the two ends of the cylinder and (ii) to pull out and freeze a spin at either end. We also investigate the relations between bulk and boundary properties and show that, for a bipartition of the cylinder, the boundary Hamiltonian defined on the edge can be written as a product of a highly nonlocal projector, which fundamentally depends upon boundary conditions, with an emergent (local) SU(2)-invariant one-dimensional (superfluid) t-J Hamiltonian, which arises due to the symmetry properties of the auxiliary spins at the edge. This multiplicative structure, a consequence of the disconnected topological sectors in the space of dimer lattice coverings, is characteristic of the topological nature of the states. For minimally entangled RVB states, it is shown that the entanglement spectrum, which reflects the properties of the (gapless or gapped) edge modes, is a subset of the spectrum of the local Hamiltonian, e.g., half of it for the kagome RVB state, providing a simple argument on the origin of the topological entanglement entropy S_0=−ln2 of the Z_2 spin liquid. We propose to use these features to probe topological phases in microscopic Hamiltonians, and some results are compared to existing density matrix renormalization group data