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 $\ensuremath{\pi}$-flux phase, and the nearest-neighbor resonating-valence-bond wave function. For these models, we show that the entanglement entropy between cylindrical regions of length
x and $L\ensuremath{-}x$, extending around a torus of length
L, depends on the dimensionless ratio
x/L. This can be well approximated on finite-size lattices by a function $\mathrm{ln}(\mathrm{sin}(\ensuremath{\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.
We systematically investigate the finite size scaling behavior of the Rényi entanglement entropy (EE) of several representative 2d quantum many-body systems between a subregion and its complement, with smooth boundaries as well as boundaries with corners. In order to reveal the subleading correction, we investigate the quantity “subtracted EE” S^s(l) = S(2l) - 2S(l) S s ( l ) = S ( 2 l ) − 2 S ( l ) for each model, which is designed to cancel out the leading perimeter law. We find that (1) for a spin-1/2 model on a 2d square lattice whose ground state is the Neel order, the coefficient of the logarithmic correction to the perimeter law is consistent with the prediction based on the Goldstone modes; (2) for the (2+1)d ( 2 + 1 ) d O(3) Wilson-Fisher quantum critical point (QCP), realized with the bilayer antiferromagnetic Heisenberg model, a logarithmic subleading correction exists when there is sharp corner of the subregion, but for subregion with a smooth boundary our data suggests the absence of the logarithmic correction to the best of our efforts; (3) for the (2+1)d ( 2 + 1 ) d SU(2) J-Q _2 2 and J-Q _3 3 model for the deconfined quantum critical point (DQCP), we find a logarithmic correction for the EE even with smooth boundary.
The multi-scale entanglement renormalization ansatz (MERA) is a tensor network that can efficiently parameterize critical ground states on a 1D lattice, and also suggestively implement some aspects of the holographic correspondence of string theory on a lattice. Extending our recent work [S. Singh, Physical Review D 97, 026012 (2018); S. Singh, N. A. McMahon, and G. K. Brennen, Phys. Rev. D 97, 026013 (2018)], we show how the MERA representation of a 1D critical ground state---which has long range entanglement---can be viewed as a strange correlator: the overlap of a 2D state with short range entanglement and a 2D product state. Strange correlators were recently introduced to map 2D symmetry protected or topologically ordered quantum states to critical systems in one lower dimension. The 2D quantum state dual to the input 1D critical state is obtained by lifting the MERA, a procedure which introduces bulk quantum degrees of freedom by inserting intertwiner tensors on each bond of the MERA tensor network. We show how this dual 2D bulk state exhibits several features of holography, for example, appearance of horizon-like holographic screens and bulk gauging of global on-site symmetries at the boundary. We also derive a quantum corrected Ryu-Takayanagi formula relating boundary entanglement entropy to bulk geodesic lengths---as measured by bulk entropy---and numerically test it for ground states of a set of unitary minimal model CFTs, as realized by 1D anyonic Heisenberg models.