We introduce for SU(2) quantum spin systems the Valence Bond Entanglement Entropy as a counting of valence bond spin singlets shared by two subsystems. For a large class of antiferromagnetic systems, it can be calculated in all dimensions with Quantum Monte Carlo simulations in the valence bond basis. We show numerically that this quantity displays all features of the von Neumann entanglement entropy for several one-dimensional systems. For two-dimensional Heisenberg models, we find a strict area law for a Valence Bond Solid state and multiplicative logarithmic corrections for the Neel phase.
We calculate the bipartite von Neumann and second Renyi entanglement entropies of the ground states of spin-1/2 dimerized Heisenberg antiferromagnets on a square lattice. Two distinct dimerization patterns are considered: columnar and staggered. In both cases, we concentrate on the valence bond solid (VBS) phase and describe such a phase with the bond-operator representation. Within this formalism, the original spin Hamiltonian is mapped into an effective interacting boson model for the triplet excitations. We study the effective Hamiltonian at the harmonic approximation and determine the spectrum of the elementary triplet excitations. We then follow an analytical procedure, which is based on a modified spin-wave theory for finite systems and was originally employed to calculate the entanglement entropies of magnetic ordered phases, and calculate the entanglement entropies of the VBS ground states. In particular, we consider one-dimensional (line) subsystems within the square lattice, a choice that allows us to consider line subsystems with sizes up to $L = 1000$. We combine such a procedure with the results of the bond-operator formalism at the harmonic level and show that, for both dimerized Heisenberg models, the entanglement entropies of the corresponding VBS ground states obey an area law as expected for gapped phases. For both columnar-dimer and staggered-dimer models, we also show that the entanglement entropies increase but do not diverge as the dimerization decreases and the system approaches the Neel--VBS quantum phase transition. Finally, the entanglement spectra associated with the VBS ground states are presented.
By relating the ground state of Temperley-Lieb hamiltonians to partition functions of 2D statistical mechanics systems on a half plane, and using a boundary Coulomb gas formalism, we obtain in closed form the valence bond entanglement entropy as well as the valence bond probability distribution in these ground states. We find in particular that for the XXX spin chain, the number N_c of valence bonds connecting a subsystem of size L to the outside goes, in the thermodynamic limit, as <N_c> = (4/pi^2) ln L, disproving a recent conjecture that this should be related with the von Neumann entropy, and thus equal to 1/(3 ln 2) ln L. Our results generalize to the Q-state Potts model.
An isotropic anti-ferromagnetic quantum state on a square lattice is characterized by symmetry arguments only. By construction, this quantum state is the result of an underlying valence bond structure without breaking any symmetry in the lattice or spin spaces. A detailed analysis of the correlations of the quantum state is given (using a mapping to a 2D classical statistical model and methods in field theory like mapping to the non-linear sigma model or bosonization techniques) as well as the results of numerical treatments (regarding exact diagonalization and variational methods). Finally, the physical relevance of the model is motivated. A comparison of the model to known anti-ferromagnetic Mott-Hubbard insulators is given by means of the two-point equal-time correlation function obtained i) numerically from the suggested state and ii) experimentally from neutron scattering on cuprates in the anti-ferromagnetic insulator phase.
We calculate exactly the von Neumann and topological entropies of the toric code as a function of system size and temperature. We do so for systems with infinite energy scale separation between magnetic and electric excitations, so that the magnetic closed loop structure is fully preserved while the electric loop structure is tampered with by thermally excited electric charges. We find that the entanglement entropy is a singular function of temperature and system size, and that the limit of zero temperature and the limit of infinite system size do not commute. From the entanglement entropy we obtain the topological entropy, which is shown to drop to half its zero-temperature value for any infinitesimal temperature in the thermodynamic limit, and remains constant as the temperature is further increased. Such discontinuous behavior is replaced by a smooth decreasing function in finite-size systems. If the separation of energy scales in the system is large but finite, we argue that our results hold at small enough temperature and finite system size, and a second drop in the topological entropy should occur as the temperature is raised so as to disrupt the magnetic loop structure by allowing the appearance of free magnetic charges. We interpret our results as an indication that the underlying magnetic and electric closed loop structures contribute equally to the topological entropy (and therefore to the topological order) in the system. Since each loop structure emph{per se} is a classical object, we interpret the quantum topological order in our system as arising from the ability of the two structures to be superimposed and appear simultaneously.
The trimer resonating valence bond (tRVB) state consisting of an equal-weight superposition of trimer coverings on a square lattice is proposed. A model Hamiltonian of the Rokhsar-Kivelson type for which the tRVB becomes the exact ground state is written. The state is shown to have $9^g$ topological degeneracy on genus g surface and support $Z_3$ vortex excitations. Correlation functions show exponential behavior with a very short correlation length consistent with the gapped spectrum. The classical problem of the degeneracy of trimer configurations is investigated by the transfer matrix method.