Anomalies in the entanglement properties of the square-lattice Heisenberg model

We compute the bipartite entanglement properties of the spin-half square-lattice Heisenberg model by a variety of numerical techniques that include valence-bond quantum Monte Carlo (QMC), stochastic series expansion QMC, high-temperature series expansions, and zero-temperature coupling constant expansions around the Ising limit. We find that the area law is always satisfied, but in addition to the entanglement entropy per unit boundary length, there are other terms that depend logarithmically on the subregion size, arising from broken symmetry in the bulk and from the existence of corners at the boundary. We find that the numerical results are anomalous in several ways. First, the bulk term arising from broken symmetry deviates from an exact calculation that can be done for a mean-field Neel state. Second, the corner logs do not agree with the known results for noninteracting Boson modes. And, third, even the finite-temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that the T -> 0 and L -> infinity limits do not commute. These calculations show that entanglement entropy demonstrates a very rich behavior in d > 1, which deserves further attention.