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:1306.6067 · Physical Review B

Phase diagram of the spin- 12 J1 - J2 Heisenberg model on a honeycomb lattice

Shou-Shu Gong, D. N. Sheng, Olexei I. Motrunich, Matthew P. A. Fisher

Published Oct 29, 2013• 113 citations
arXivHTMLPDF

Abstract

We use the density matrix renormalization group (DMRG) algorithm to study the phase diagram of the spin-$\frac{1}{2}$ Heisenberg model on a honeycomb lattice with first (${J}_{1}$) and second (${J}_{2}$) neighbor antiferromagnetic interactions, where a ${Z}_{2}$ spin liquid region has been proposed. By implementing SU(2) symmetry in the DMRG code, we are able to obtain accurate results for long cylinders with a width slightly over 15 lattice spacings and a torus up to the size $N=2\ifmmode\times\else\texttimes\fi{}6\ifmmode\times\else\texttimes\fi{}6$. With increasing ${J}_{2}$, we find a N\'eel phase with a vanishing spin gap and a plaquette valence-bond (PVB) phase with a nonzero spin gap. By extrapolating the square of the staggered magnetic moment ${m}_{s}^{2}$ on finite-size cylinders to the thermodynamic limit, we find the N\'eel order vanishing at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.22$. For $0.25<{J}_{2}/{J}_{1}\ensuremath{\le}0.35$, we find a possible PVB order, which shows a fast growing PVB decay length with increasing system width. For $0.22<{J}_{2}/{J}_{1}<0.25$, both spin and dimer orders vanish in the thermodynamic limit, which is consistent with a possible spin liquid phase. We present calculations of the topological entanglement entropy, compare the DMRG results with the variational Monte Carlo, and discuss possible scenarios in the thermodynamic limit for this region.

DOI
10.1103/physrevb.88.165138

Related papers

Quantum phase diagram of the spin-1 J1−J2 Heisenberg model on the honeycomb lattice

Shou-Shu Gong, Wei Kang Zhu, D. N. Sheng

Strongly correlated systems with geometric frustrations can host the emergent phases of matter with unconventional properties. Here, we study the spin $S=1$ Heisenberg model on the honeycomb lattice with the antiferromagnetic first- (${J}_{1}$) and second-neighbor (${J}_{2}$) interactions ($0.0\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$) by means of density-matrix renormalization group. In the parameter regime ${J}_{2}/{J}_{1}\ensuremath{\lesssim}0.27$, the system sustains a N\'eel antiferromagnetic phase. At the large ${J}_{2}$ side ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.32$, a stripe antiferromagnetic phase is found. Between the two magnetic ordered phases $0.27\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.32$, we find a nonmagnetic intermediate region with a plaquette valence-bond order. Although our calculations are limited within a six unit-cell width on cylinder, we present evidence that this plaquette state could be a strong candidate for this nonmagnetic region in the thermodynamic limit. We also briefly discuss the nature of the quantum phase transitions in the system. We gain further insight into the nonmagnetic phases in the spin-1 system by comparing its phase diagram with the spin-$1/2$ system.

arXiv:1508.00515•Physical Review B•Nov 6, 2015•44 citations

Plaquette Ordered Phase and Quantum Phase Diagram in the Spin- 12 J1−J2 Square Heisenberg Model

Shou-Shu Gong, Wei Kang Zhu, D. N. Sheng +2 more

We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J(1) and J(2), which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z(2) spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J(2), we find a Néel phase and a plaquette valence-bond (PVB) phase with a finite spin gap. From the finite-size scaling of the magnetic order parameter, we estimate that the Néel order vanishes at J(2)/J(1)≃0.44. For 0.5<J(2)/J(1)<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J(2)/J(1)<0.5, spin order, dimer order, and spin gap are small on finite-size systems, which is consistent with a near-critical behavior. The critical exponents obtained from the finite-size spin and dimer correlations could be compatible with the deconfined criticality in this small region. We compare and contrast our results with earlier numerical studies.

arXiv:1311.5962•Physical Review Letters•Jul 7, 2014•275 citations

Phase diagram of a frustrated quantum antiferromagnet on the honeycomb lattice: Magnetic order versus valence-bond crystal formation

A. Fabricio Albuquerque, David R. Schwandt, Balázs Hetényi +3 more

We present a comprehensive computational study of the phase diagram of the frustrated $S=1/2$ Heisenberg antiferromagnet on the honeycomb lattice, with second-nearest $({J}_{2})$ and third-neighbor $({J}_{3})$ couplings. Using a combination of exact diagonalizations (EDs) of the original spin model, of the Hamiltonian projected into the nearest-neighbor short-range valence-bond basis, and of an effective quantum dimer model, as well as a self-consistent cluster mean-field theory, we determine the boundaries of several magnetically ordered phases in the region ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, and find a sizable magnetically disordered region in between. We characterize part of this magnetically disordered phase as a plaquette valence-bond crystal phase. At larger ${J}_{2}$, we locate a sizable region in which staggered valence-bond crystal correlations are found to be important, either due to genuine valence-bond crystal (VBC) ordering or as a consequence of magnetically ordered phases, which break lattice rotational symmetry. Furthermore, we find that a particular parameter-free Gutzwiller projected tight-binding wave function has remarkably accurate energies compared to finite-size extrapolated ED energies along the transition line from conventional N\'eel to plaquette VBC phases, a fact that points to possibly interesting critical behavior---such as a deconfined critical point---across this transition. We also comment on the relevance of this spin model to model the spin liquid region found in the half filled Hubbard model on the honeycomb lattice.

arXiv:1102.5325•Physical Review B•Jul 6, 2011•204 citations

Weak Plaquette Valence Bond Order in the S=1/2 Honeycomb J1−J2 Heisenberg Model

Zhenyue Zhu, David A. Huse, Steven R. White

Using the density matrix renormalization group, we investigate the S = 1/2 Heisenberg model on the honeycomb lattice with first (J(1)) and second (J(2)) neighbor interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For J(2)/J(1) near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for J(2) ≲ 0.26 and by a valence bond crystal for J(2) ≳ 0.36. The longest correlation length that we find in this intermediate phase is for plaquette valence bond order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak plaquette valence bond order.

arXiv:1212.6322•Physical Review Letters•Mar 18, 2013•109 citations

Competition between spin liquids and valence-bond order in the frustrated spin- 12 Heisenberg model on the honeycomb lattice

Francesco Ferrari, Samuel Bieri, Federico Becca

Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor ${J}_{1}$ and second-neighbor ${J}_{2}$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. N\'eel order is stabilized for small values of the frustrating second-neighbor coupling. Increasing the ratio ${J}_{2}/{J}_{1}$, we find strong evidence for a continuous transition to a nonmagnetic phase at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.23$. Close to the transition point, the Gutzwiller-projected uniform resonating valence-bond state gives an excellent approximation to the exact ground-state energy. For $0.23\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.36$, a gapless ${Z}_{2}$ spin liquid with Dirac nodes competes with a plaquette valence-bond solid. In contrast, the gapped spin liquid considered in previous works has significantly higher variational energy. Although the plaquette valence-bond order is expected to be present as soon as the N\'eel order melts, this ordered state becomes clearly favored only for ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.3$. Finally, for $0.36\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$, a valence-bond solid with columnar order takes over as the ground state, being also lower in energy than the magnetic state with collinear order. We perform a detailed finite-size scaling and standard data collapse analysis, and we discuss the possibility of a deconfined quantum critical point separating the N\'eel antiferromagnet from the plaquette valence-bond solid.

arXiv:1706.07810•Physical review. B./Physical review. B•Sep 1, 2017•46 citations

Spin liquid ground state of the spin- 1 2 square J 1 - J 2 Heisenberg model

Hong‐Chen Jiang, Hong Yao, Leon Balents

We perform highly accurate density matrix renormalization group (DMRG) simulations to investigate the ground-state properties of the spin-$1/2$ antiferromagnetic square lattice Heisenberg ${J}_{1}$-${J}_{2}$ model. Based on studies of numerous long cylinders with circumferences of up to 14 lattice spacings, we obtain strong evidence for a topological quantum spin liquid state in the region $0.41\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}0.62$, separating conventional N\'eel and striped antiferromagnetic states for smaller and larger ${J}_{2}/{J}_{1}$, respectively. The quantum spin liquid is characterized numerically by the absence of magnetic or valence bond solid order, and nonzero singlet and triplet energy gaps. Furthermore, we positively identify its topological nature by measuring a nonzero topological entanglement entropy $\ensuremath{\gamma}=0.70\ifmmode\pm\else\textpm\fi{}0.02$, extremely close to $\ensuremath{\gamma}=\mathrm{ln}(2)\ensuremath{\approx}0.69$ (expected for a ${Z}_{2}$ quantum spin liquid) and a nontrivial finite size dimerization effect depending upon the parity of the circumference of the cylinder. We also point out that a valence bond solid, and indeed any discrete symmetry breaking state, would be expected to show a constant correction to the entanglement entropy of opposite sign to the topological entanglement entropy.

arXiv:1112.2241•Physical Review B•Jul 18, 2012•359 citations
We use the density matrix renormalization group (DMRG) algorithm to study the phase diagram of the spin-12\frac{1}{2}21​ Heisenberg model on a honeycomb lattice with first (J1{J}_{1}J1​) and second (J2{J}_{2}J2​) neighbor antiferromagnetic interactions, where a Z2{Z}_{2}Z2​ spin liquid region has been proposed. By implementing SU(2) symmetry in the DMRG code, we are able to obtain accurate results for long cylinders with a width slightly over 15 lattice spacings and a torus up to the size N=2\ifmmode×\else\texttimes\fi6\ifmmode×\else\texttimes\fi6N=2\ifmmode\times\else\texttimes\fi{}6\ifmmode\times\else\texttimes\fi{}6N=2\ifmmode×\else\texttimes\fi6\ifmmode×\else\texttimes\fi6. With increasing J2{J}_{2}J2​, we find a N'eel phase with a vanishing spin gap and a plaquette valence-bond (PVB) phase with a nonzero spin gap. By extrapolating the square of the staggered magnetic moment ms2{m}_{s}^{2}ms2​ on finite-size cylinders to the thermodynamic limit, we find the N'eel order vanishing at J2/J1\ensuremath≃0.22{J}_{2}/{J}_{1}\ensuremath{\simeq}0.22J2​/J1​\ensuremath≃0.22. For 0.25<J2/J1\ensuremath≤0.350.25<{J}_{2}/{J}_{1}\ensuremath{\le}0.350.25<J2​/J1​\ensuremath≤0.35, we find a possible PVB order, which shows a fast growing PVB decay length with increasing system width. For 0.22<J2/J1<0.250.22<{J}_{2}/{J}_{1}<0.250.22<J2​/J1​<0.25, both spin and dimer orders vanish in the thermodynamic limit, which is consistent with a possible spin liquid phase. We present calculations of the topological entanglement entropy, compare the DMRG results with the variational Monte Carlo, and discuss possible scenarios in the thermodynamic limit for this region.
Strongly correlated systems with geometric frustrations can host the emergent phases of matter with unconventional properties. Here, we study the spin S=1S=1S=1 Heisenberg model on the honeycomb lattice with the antiferromagnetic first- (J1{J}_{1}J1​) and second-neighbor (J2{J}_{2}J2​) interactions ($0.0\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$) by means of density-matrix renormalization group. In the parameter regime ${J}_{2}/{J}_{1}\ensuremath{\lesssim}0.27$, the system sustains a N\'eel antiferromagnetic phase. At the large J2{J}_{2}J2​ side ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.32$, a stripe antiferromagnetic phase is found. Between the two magnetic ordered phases $0.27\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.32$, we find a nonmagnetic intermediate region with a plaquette valence-bond order. Although our calculations are limited within a six unit-cell width on cylinder, we present evidence that this plaquette state could be a strong candidate for this nonmagnetic region in the thermodynamic limit. We also briefly discuss the nature of the quantum phase transitions in the system. We gain further insight into the nonmagnetic phases in the spin-1 system by comparing its phase diagram with the spin-1/21/21/2 system.
We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J(1) and J(2), which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z(2) spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J(2), we find a Néel phase and a plaquette valence-bond (PVB) phase with a finite spin gap. From the finite-size scaling of the magnetic order parameter, we estimate that the Néel order vanishes at J(2)/J(1)≃0.44. For 0.5<J(2)/J(1)<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J(2)/J(1)<0.5, spin order, dimer order, and spin gap are small on finite-size systems, which is consistent with a near-critical behavior. The critical exponents obtained from the finite-size spin and dimer correlations could be compatible with the deconfined criticality in this small region. We compare and contrast our results with earlier numerical studies.
We present a comprehensive computational study of the phase diagram of the frustrated S=1/2S=1/2S=1/2 Heisenberg antiferromagnet on the honeycomb lattice, with second-nearest (J2)({J}_{2})(J2​) and third-neighbor (J3)({J}_{3})(J3​) couplings. Using a combination of exact diagonalizations (EDs) of the original spin model, of the Hamiltonian projected into the nearest-neighbor short-range valence-bond basis, and of an effective quantum dimer model, as well as a self-consistent cluster mean-field theory, we determine the boundaries of several magnetically ordered phases in the region ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, and find a sizable magnetically disordered region in between. We characterize part of this magnetically disordered phase as a plaquette valence-bond crystal phase. At larger J2{J}_{2}J2​, we locate a sizable region in which staggered valence-bond crystal correlations are found to be important, either due to genuine valence-bond crystal (VBC) ordering or as a consequence of magnetically ordered phases, which break lattice rotational symmetry. Furthermore, we find that a particular parameter-free Gutzwiller projected tight-binding wave function has remarkably accurate energies compared to finite-size extrapolated ED energies along the transition line from conventional N\'eel to plaquette VBC phases, a fact that points to possibly interesting critical behavior---such as a deconfined critical point---across this transition. We also comment on the relevance of this spin model to model the spin liquid region found in the half filled Hubbard model on the honeycomb lattice.
Using the density matrix renormalization group, we investigate the S = 1/2 Heisenberg model on the honeycomb lattice with first (J(1)) and second (J(2)) neighbor interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For J(2)/J(1) near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for J(2) ≲ 0.26 and by a valence bond crystal for J(2) ≳ 0.36. The longest correlation length that we find in this intermediate phase is for plaquette valence bond order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak plaquette valence bond order.
Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor J1{J}_{1}J1​ and second-neighbor J2{J}_{2}J2​ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. N\'eel order is stabilized for small values of the frustrating second-neighbor coupling. Increasing the ratio J2/J1{J}_{2}/{J}_{1}J2​/J1​, we find strong evidence for a continuous transition to a nonmagnetic phase at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.23$. Close to the transition point, the Gutzwiller-projected uniform resonating valence-bond state gives an excellent approximation to the exact ground-state energy. For $0.23\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.36$, a gapless Z2{Z}_{2}Z2​ spin liquid with Dirac nodes competes with a plaquette valence-bond solid. In contrast, the gapped spin liquid considered in previous works has significantly higher variational energy. Although the plaquette valence-bond order is expected to be present as soon as the N\'eel order melts, this ordered state becomes clearly favored only for ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.3$. Finally, for $0.36\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$, a valence-bond solid with columnar order takes over as the ground state, being also lower in energy than the magnetic state with collinear order. We perform a detailed finite-size scaling and standard data collapse analysis, and we discuss the possibility of a deconfined quantum critical point separating the N\'eel antiferromagnet from the plaquette valence-bond solid.
We perform highly accurate density matrix renormalization group (DMRG) simulations to investigate the ground-state properties of the spin-1/21/21/2 antiferromagnetic square lattice Heisenberg J1{J}_{1}J1​-J2{J}_{2}J2​ model. Based on studies of numerous long cylinders with circumferences of up to 14 lattice spacings, we obtain strong evidence for a topological quantum spin liquid state in the region $0.41\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}0.62$, separating conventional N\'eel and striped antiferromagnetic states for smaller and larger J2/J1{J}_{2}/{J}_{1}J2​/J1​, respectively. The quantum spin liquid is characterized numerically by the absence of magnetic or valence bond solid order, and nonzero singlet and triplet energy gaps. Furthermore, we positively identify its topological nature by measuring a nonzero topological entanglement entropy $\ensuremath{\gamma}=0.70\ifmmode\pm\else\textpm\fi{}0.02$, extremely close to $\ensuremath{\gamma}=\mathrm{ln}(2)\ensuremath{\approx}0.69$ (expected for a Z2{Z}_{2}Z2​ quantum spin liquid) and a nontrivial finite size dimerization effect depending upon the parity of the circumference of the cylinder. We also point out that a valence bond solid, and indeed any discrete symmetry breaking state, would be expected to show a constant correction to the entanglement entropy of opposite sign to the topological entanglement entropy.