Finite Size Scaling of Mutual Information: A Scalable Simulation

We develop a quantum Monte Carlo procedure to compute the Renyi mutual information of an interacting quantum many-body system at non-zero temperature. Performing simulations on a spin-1/2 XXZ model, we observe that for a subregion of fixed size embedded in a system of size L, the mutual information converges at large L to a limiting function which displays non-monotonic temperature behavior corresponding to the onset of correlations. For a region of size L/2 embedded in a system of size L, the mutual information divided by L converges to a limiting function of temperature, with apparently nontrivial corrections near critical points.Comment: 4 pages, 4 figure

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 non-interacting Boson modes. And, third, even the finite temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that 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.Comment: 12 pages, 7 figures, 2 tables. Numerical values in Table I correcte

Entanglement scaling in two-dimensional gapless systems

We numerically determine subleading scaling terms in the ground-state entanglement entropy of several two-dimensional (2D) gapless systems, including a Heisenberg model with N\'eel order, a free Dirac fermion in the {\pi}-flux phase, and the nearest-neighbor resonating-valence-bond wavefunction. For these models, we show that the entanglement entropy between cylindrical regions of length x and L - x, extending around a torus of length L, depends upon the dimensionless ratio x/L. This can be well-approximated on finite-size lattices by a function ln(sin({\pi}x/L)), akin to the familiar chord-length dependence in one dimension. We provide evidence, however, that the precise form of this bulk-dependent contribution is a more general function in the 2D thermodynamic limit.Comment: 5 pages, 5 figure