Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2005.07745 · Physical review. B./Physical review. B

Spin-wave study of entanglement and Rényi entropy for coplanar and collinear magnetic orders in two-dimensional quantum Heisenberg antiferromagnets

Dag-Vidar Bauer, John Ove Fjærestad

Published May 13, 2020• 4 citations
arXivHTMLPDF

Abstract

We use modified linear spin-wave theory (MLSWT) to study ground-state entanglement for a length-$L$ line subsystem in $L\ifmmode\times\else\texttimes\fi{}L$ square- and triangular-lattice quantum Heisenberg antiferromagnets with coplanar spiral magnetic order with ordering vector $\mathbit{Q}=(q,q)$ and ${N}_{G}=3$ Goldstone modes, except if $q=\ensuremath{\pi}$ (collinear order, ${N}_{G}=2$). Generalizing earlier MLSWT results for $q=\ensuremath{\pi}$ to commensurate spiral order with $s\ensuremath{\ge}3$ sublattices ($q=2\ensuremath{\pi}r/s$ with $r$ and $s$ coprime), we find analytically for large $L$ a universal and $n$-independent subleading term $({N}_{G}/2)lnL$ in the R\'enyi entropy ${S}_{n}$, associated with ${L}^{1/2}$ scaling of ${\ensuremath{\lambda}}_{0}$ and ${\ensuremath{\lambda}}_{\ifmmode\pm\else\textpm\fi{}q}$, with ${\ensuremath{\lambda}}_{0}\ensuremath{\ne}{\ensuremath{\lambda}}_{\ifmmode\pm\else\textpm\fi{}q}$ for spiral order; here ${{\ensuremath{\lambda}}_{{k}_{y}}}$ are the $L$ mode occupation numbers of the entanglement Hamiltonian. The term $(3/2)lnL$ in ${S}_{n}$ agrees with a nonlinear sigma model study of $s=3$ spiral order ($q=2\ensuremath{\pi}/3$). These and other properties of ${S}_{n}$ and ${\ensuremath{\lambda}}_{{k}_{y}}$ are explored numerically for an anisotropic nearest-neighbor triangular-lattice model for which $q$ varies in the spiral phase.

DOI
10.1103/physrevb.101.195124

Related papers

Spin-wave approach for entanglement entropies of the J1−J2 Heisenberg antiferromagnet on the square lattice

Nicolas Laflorencie, David J. Luitz, Fabien Alet

Using a modified spin-wave theory which artificially restores zero sublattice magnetization on finite lattices, we investigate the entanglement properties of the N\'eel ordered ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg antiferromagnet on the square lattice. Different kinds of subsystem geometries are studied, either corner-free (line, strip) or with sharp corners (square). Contributions from the ${n}_{G}=2$ Nambu-Goldstone modes give additive logarithmic corrections with a prefactor ${n}_{G}/2$ independent of the R\'enyi index. On the other hand, $\ensuremath{\pi}/2$ corners lead to additional (negative) logarithmic corrections with a prefactor ${l}_{q}^{c}$ which does depend on both ${n}_{G}$ and the R\'enyi index $q$, in good agreement with scalar field theory predictions. By varying the second neighbor coupling ${J}_{2}$ we also explore universality across the N\'eel ordered side of the phase diagram of the ${J}_{1}\ensuremath{-}{J}_{2}$ antiferromagnet, from the frustrated side $0<{J}_{2}/{J}_{1}<1/2$ where the area law term is maximal, to the strongly ferromagnetic regime $\ensuremath{-}{J}_{2}/{J}_{1}\ensuremath{\gg}1$ with a purely logarithmic growth ${S}_{q}=\frac{{n}_{G}}{2}lnN$, thus recovering the mean-field limit for a subsystem of $N$ sites. Finally, a universal subleading constant term ${\ensuremath{\gamma}}_{q}^{\mathrm{ord}}$ is extracted in the case of strip subsystems, and a direct relation is found (in the large-$S$ limit) with the same constant extracted from free lattice systems. The singular limit of vanishing aspect ratios is also explored, where we identify for ${\ensuremath{\gamma}}_{q}^{\text{ord}}$ a regular part and a singular component, explaining the discrepancy of the linear scaling term for fixed width vs fixed aspect ratio subsystems.

arXiv:1506.03703•Physical Review B•Sep 11, 2015•32 citations

Numerical study of incommensurability of the spiral state on spin- 12 spatially anisotropic triangular antiferromagnets using entanglement renormalization

Kenji Harada

The ground state of an $S=1/2$ antiferromagnetic Heisenberg model on a spatially anisotropic triangular lattice, which is an effective model of Mott insulators on a triangular layer of organic charge transfer salts or Cs${}_{2}$CuCl${}_{4}$, is numerically studied. We apply a numerical variational method by using a tensor network with entanglement renormalization, which improves the capability of describing a quantum state. Magnetic ground states are identified for $0.7\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}1$ in the thermodynamic limit, where ${J}_{1}$ and ${J}_{2}$ denote the innerchain and interchain coupling constants, respectively. Except for the isotropic case (${J}_{1}={J}_{2}$), the magnetic structure is spiral with an incommensurate wave vector that is different from the classical one. The quantum fluctuation weakens the effective coupling between chains, but the magnetic order remains in the thermodynamic limit. In addition, the incommensurate wave number is in good agreement with that of the series expansion method.

arXiv:1208.4306•Physical Review B•Nov 16, 2012•28 citations

Shannon-R\\'enyi entropy for Nambu-Goldstone modes in two dimensions

Grégoire Misguich, Vincent Pasquier, Masaki Oshikawa

We study the scaling of the (basis dependent) Shannon entropy for two-dimensional quantum antiferromagnets with N\'eel long-range order. We use a massless free-field description of the gapless spin wave modes and phase space arguments to treat the fact that the finite-size ground state is rotationally symmetric, while there are degenerate physical ground states which break the symmetry. Our results show that the Shannon entropy (and its R\'enyi generalizations) possesses some universal logarithmic term proportional to the number $N_\text{NG}$ of Nambu-Goldstone modes. In the case of a torus, we show that $S_{n>1} \simeq \mathcal{O}(N)+ \frac{N_\text{NG}}{4}\frac{n}{n-1} \log{N}$ and $S_1 \simeq \mathcal{O}(N) - \frac{N_\text{NG}}{4} \log{N}$, where $N$ is the total number of sites and $n$ the R\'enyi index. The result for $n>1$ is in reasonable agreement with the quantum Monte Carlo results of Luitz et al. [Phys. Rev. Lett. 112, 057203 (2014)], and qualitatively similar to those obtained previously for the entanglement entropy. The Shannon entropy of a line subsystem (embedded in the two-dimensional system) is also considered. Finally, we present some density-matrix renormalization group (DMRG) calculations for a spin$\frac{1}{2}$ XY model on the square lattice in a cylinder geometry. These numerical data confirm our findings for logarithmic terms in the $n=\infty$ R\'enyi entropy (also called $-\log{p_{\rm max}}$). They also reveal some universal dependence on the cylinder aspect ratio, in good agreement with the fact that, in that case, $p_{\rm max}$ is related to a non-compact free-boson partition function in dimension 1+1.

•Jan 1, 2016•0 citations

Entanglement entropy of the two-dimensional Heisenberg antiferromagnet

Hao Song, Nicolas Laflorencie, Stephan Rachel +1 more

We compute the von Neumann and generalized R\'enyi entanglement entropies in the ground-state of the spin-$1/2$ antiferromagnetic Heisenberg model on the square lattice using the modified spin-wave theory for finite lattices. The addition of a staggered magnetic field to regularize the zero modes associated with symmetry breaking is shown to be essential for obtaining well-behaved values for the entanglement entropy. The von Neumann and R\'enyi entropies obey an area law with additive logarithmic corrections, and are in good quantitative agreement with numerical results from valence bond quantum Monte Carlo and density matrix renormalization group calculations. We also compute the spin fluctuations and observe a multiplicative logarithmic correction to the area law in excellent agreement with quantum Monte Carlo calculations.

arXiv:1103.1636•Physical Review B•Jun 21, 2011•96 citations

A numerical study for two-dimensional spin 1/2 antiferromagnets: a generalization of Entanglement Perturbation Theory to two-dimensional systems

Lihua Wang, Tomonori Shirakawa, Hiroshi Watanabe +1 more

Two-dimensional spin 1/2 antiferromagnetic Heisenberg models are numerically studied using Entanglement Perturbation Theory, where the ground state wave function is described by a product of local matrices defined at every site, and each matrix is optimized variationally to minimize the energy. We first apply this method for the spin 1/2 antiferromagnetic Heisenberg model on the square lattices to reproduce correctly the known ground state energy and the spin structure factors. Then, we study the spin 1/2 antiferromagnetic Heisenberg model on the triangular lattice with spatially anisotropic nearest neighbor couplings, J and J ', to establish the ground state phase diagram as a function of J'/ J.

Journal of Physics Conference Series•Jan 1, 2011•0 citations

Entanglement entropy for the valence bond solid phases of two-dimensional dimerized Heisenberg antiferromagnets

Leonardo S. G. Leite, R. L. Doretto

We calculate the bipartite von Neumann and second R\'enyi 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}^{\ensuremath{'}}\phantom{\rule{0.16em}{0ex}}=\phantom{\rule{0.16em}{0ex}}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 they seem to not diverge as the dimerization decreases and the system approaches the N\'eel-VBS quantum phase transition. Finally, the entanglement spectra associated with the VBS ground states are presented.

arXiv:1902.04355•Physical review. B./Physical review. B•Jul 10, 2019•4 citations
We use modified linear spin-wave theory (MLSWT) to study ground-state entanglement for a length-LLL line subsystem in L\ifmmode×\else\texttimes\fiLL\ifmmode\times\else\texttimes\fi{}LL\ifmmode×\else\texttimes\fiL square- and triangular-lattice quantum Heisenberg antiferromagnets with coplanar spiral magnetic order with ordering vector \mathbitQ=(q,q)\mathbit{Q}=(q,q)\mathbitQ=(q,q) and NG=3{N}_{G}=3NG​=3 Goldstone modes, except if q=\ensuremathπq=\ensuremath{\pi}q=\ensuremathπ (collinear order, NG=2{N}_{G}=2NG​=2). Generalizing earlier MLSWT results for q=\ensuremathπq=\ensuremath{\pi}q=\ensuremathπ to commensurate spiral order with s\ensuremath≥3s\ensuremath{\ge}3s\ensuremath≥3 sublattices (q=2\ensuremathπr/sq=2\ensuremath{\pi}r/sq=2\ensuremathπr/s with rrr and sss coprime), we find analytically for large LLL a universal and nnn-independent subleading term (NG/2)lnL({N}_{G}/2)lnL(NG​/2)lnL in the R'enyi entropy Sn{S}_{n}Sn​, associated with L1/2{L}^{1/2}L1/2 scaling of \ensuremathλ0{\ensuremath{\lambda}}_{0}\ensuremathλ0​ and \ensuremathλ\ifmmode±\else\textpm\fiq{\ensuremath{\lambda}}_{\ifmmode\pm\else\textpm\fi{}q}\ensuremathλ\ifmmode±\else\textpm\fiq​, with \ensuremathλ0\ensuremath≠\ensuremathλ\ifmmode±\else\textpm\fiq{\ensuremath{\lambda}}_{0}\ensuremath{\ne}{\ensuremath{\lambda}}_{\ifmmode\pm\else\textpm\fi{}q}\ensuremathλ0​\ensuremath=\ensuremathλ\ifmmode±\else\textpm\fiq​ for spiral order; here \ensuremathλky{{\ensuremath{\lambda}}_{{k}_{y}}}\ensuremathλky​​ are the LLL mode occupation numbers of the entanglement Hamiltonian. The term (3/2)lnL(3/2)lnL(3/2)lnL in Sn{S}_{n}Sn​ agrees with a nonlinear sigma model study of s=3s=3s=3 spiral order (q=2\ensuremathπ/3q=2\ensuremath{\pi}/3q=2\ensuremathπ/3). These and other properties of Sn{S}_{n}Sn​ and \ensuremathλky{\ensuremath{\lambda}}_{{k}_{y}}\ensuremathλky​​ are explored numerically for an anisotropic nearest-neighbor triangular-lattice model for which qqq varies in the spiral phase.
Using a modified spin-wave theory which artificially restores zero sublattice magnetization on finite lattices, we investigate the entanglement properties of the N\'eel ordered ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg antiferromagnet on the square lattice. Different kinds of subsystem geometries are studied, either corner-free (line, strip) or with sharp corners (square). Contributions from the nG=2{n}_{G}=2nG​=2 Nambu-Goldstone modes give additive logarithmic corrections with a prefactor nG/2{n}_{G}/2nG​/2 independent of the R\'enyi index. On the other hand, $\ensuremath{\pi}/2$ corners lead to additional (negative) logarithmic corrections with a prefactor lqc{l}_{q}^{c}lqc​ which does depend on both nG{n}_{G}nG​ and the R\'enyi index qqq, in good agreement with scalar field theory predictions. By varying the second neighbor coupling J2{J}_{2}J2​ we also explore universality across the N\'eel ordered side of the phase diagram of the ${J}_{1}\ensuremath{-}{J}_{2}$ antiferromagnet, from the frustrated side 0<J2/J1<1/20<{J}_{2}/{J}_{1}<1/20<J2​/J1​<1/2 where the area law term is maximal, to the strongly ferromagnetic regime $\ensuremath{-}{J}_{2}/{J}_{1}\ensuremath{\gg}1$ with a purely logarithmic growth Sq=nG2lnN{S}_{q}=\frac{{n}_{G}}{2}lnNSq​=2nG​​lnN, thus recovering the mean-field limit for a subsystem of NNN sites. Finally, a universal subleading constant term ${\ensuremath{\gamma}}_{q}^{\mathrm{ord}}$ is extracted in the case of strip subsystems, and a direct relation is found (in the large-SSS limit) with the same constant extracted from free lattice systems. The singular limit of vanishing aspect ratios is also explored, where we identify for ${\ensuremath{\gamma}}_{q}^{\text{ord}}$ a regular part and a singular component, explaining the discrepancy of the linear scaling term for fixed width vs fixed aspect ratio subsystems.
The ground state of an S=1/2S=1/2S=1/2 antiferromagnetic Heisenberg model on a spatially anisotropic triangular lattice, which is an effective model of Mott insulators on a triangular layer of organic charge transfer salts or Cs2{}_{2}2​CuCl4{}_{4}4​, is numerically studied. We apply a numerical variational method by using a tensor network with entanglement renormalization, which improves the capability of describing a quantum state. Magnetic ground states are identified for $0.7\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}1$ in the thermodynamic limit, where J1{J}_{1}J1​ and J2{J}_{2}J2​ denote the innerchain and interchain coupling constants, respectively. Except for the isotropic case (J1=J2{J}_{1}={J}_{2}J1​=J2​), the magnetic structure is spiral with an incommensurate wave vector that is different from the classical one. The quantum fluctuation weakens the effective coupling between chains, but the magnetic order remains in the thermodynamic limit. In addition, the incommensurate wave number is in good agreement with that of the series expansion method.
We study the scaling of the (basis dependent) Shannon entropy for two-dimensional quantum antiferromagnets with N\'eel long-range order. We use a massless free-field description of the gapless spin wave modes and phase space arguments to treat the fact that the finite-size ground state is rotationally symmetric, while there are degenerate physical ground states which break the symmetry. Our results show that the Shannon entropy (and its R\'enyi generalizations) possesses some universal logarithmic term proportional to the number NNGN_\text{NG}NNG​ of Nambu-Goldstone modes. In the case of a torus, we show that Sn>1≃O(N)+NNG4nn−1log⁡NS_{n>1} \simeq \mathcal{O}(N)+ \frac{N_\text{NG}}{4}\frac{n}{n-1} \log{N}Sn>1​≃O(N)+4NNG​​n−1n​logN and S1≃O(N)−NNG4log⁡NS_1 \simeq \mathcal{O}(N) - \frac{N_\text{NG}}{4} \log{N}S1​≃O(N)−4NNG​​logN, where NNN is the total number of sites and nnn the R\'enyi index. The result for n>1n>1n>1 is in reasonable agreement with the quantum Monte Carlo results of Luitz et al. [Phys. Rev. Lett. 112, 057203 (2014)], and qualitatively similar to those obtained previously for the entanglement entropy. The Shannon entropy of a line subsystem (embedded in the two-dimensional system) is also considered. Finally, we present some density-matrix renormalization group (DMRG) calculations for a spin12\frac{1}{2}21​ XY model on the square lattice in a cylinder geometry. These numerical data confirm our findings for logarithmic terms in the n=∞n=\inftyn=∞ R\'enyi entropy (also called −log⁡pmax-\log{p_{\rm max}}−logpmax​). They also reveal some universal dependence on the cylinder aspect ratio, in good agreement with the fact that, in that case, pmaxp_{\rm max}pmax​ is related to a non-compact free-boson partition function in dimension 1+1.
We compute the von Neumann and generalized R\'enyi entanglement entropies in the ground-state of the spin-1/21/21/2 antiferromagnetic Heisenberg model on the square lattice using the modified spin-wave theory for finite lattices. The addition of a staggered magnetic field to regularize the zero modes associated with symmetry breaking is shown to be essential for obtaining well-behaved values for the entanglement entropy. The von Neumann and R\'enyi entropies obey an area law with additive logarithmic corrections, and are in good quantitative agreement with numerical results from valence bond quantum Monte Carlo and density matrix renormalization group calculations. We also compute the spin fluctuations and observe a multiplicative logarithmic correction to the area law in excellent agreement with quantum Monte Carlo calculations.
Two-dimensional spin 1/2 antiferromagnetic Heisenberg models are numerically studied using Entanglement Perturbation Theory, where the ground state wave function is described by a product of local matrices defined at every site, and each matrix is optimized variationally to minimize the energy. We first apply this method for the spin 1/2 antiferromagnetic Heisenberg model on the square lattices to reproduce correctly the known ground state energy and the spin structure factors. Then, we study the spin 1/2 antiferromagnetic Heisenberg model on the triangular lattice with spatially anisotropic nearest neighbor couplings, J and J ', to establish the ground state phase diagram as a function of J'/ J.
We calculate the bipartite von Neumann and second R\'enyi 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}^{\ensuremath{'}}\phantom{\rule{0.16em}{0ex}}=\phantom{\rule{0.16em}{0ex}}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 they seem to not diverge as the dimerization decreases and the system approaches the N\'eel-VBS quantum phase transition. Finally, the entanglement spectra associated with the VBS ground states are presented.