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:1311.5962 · Physical Review Letters

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, Olexei I. Motrunich, Matthew P. A. Fisher

Published Jul 7, 2014• 275 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.113.027201

Related papers

Plaquette valence bond state in the spin- 12 J1−J2 XY model on the square lattice

Yang‐Hao Chan, Hong‐Chen Jiang, Ying-Cheng Chen

We studied the ground-state phase diagram of spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model on the square lattice with first-neighbor ${J}_{1}$ and second-neighbor ${J}_{2}$ antiferromagnetic interactions using both infinite density matrix renormalization group (iDMRG) and DMRG approaches. We show that a plaquette valence bond phase is realized in an intermediate region $0.50\ensuremath{\le}{J}_{2}/{J}_{1}<0.54$ between a N\'eel magnetic ordered phase at ${J}_{2}/{J}_{1}<0.50$ and a stripy magnetic ordered phase at ${J}_{2}/{J}_{1}\ensuremath{\ge}0.54$. The plaquette valence bond phase is characterized by finite dimer orders in both the horizontal and vertical directions. Contrary to the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model, we do not find numerical evidence for a quantum spin liquid phase in the ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model.

arXiv:2212.04863•Physical review. B./Physical review. B•Jun 1, 2023•3 citations

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

Shou-Shu Gong, D. N. Sheng, Olexei I. Motrunich +1 more

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.

arXiv:1306.6067•Physical Review B•Oct 29, 2013•113 citations

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

Tensor-product state approach to spin- 12 square J1−J2 antiferromagnetic Heisenberg model: Evidence for deconfined quantum criticality

Ling Wang, Zheng‐Cheng Gu, Frank Verstraete +1 more

The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point ${J}_{2}^{{c}_{1}}=0.572(5){J}_{1}$ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at ${J}_{2}/{J}_{1}=0.6$. These results are consistent with a potential gapless $U(1)$ spin-liquid phase. However, since the $U(1)$ spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at ${J}_{2}^{{c}_{1}}$. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.

arXiv:1112.3331•Physical review. B./Physical review. B•Aug 19, 2016•102 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

Plaquette valence bond ordering in a J 1 – J 2 Heisenberg antiferromagnet on a honeycomb lattice

H Mosadeq, F Shahbazi, S A Jafari

We study an S = 1/2 Heisenberg model on the honeycomb lattice with first and second neighbor antiferromagnetic exchange (J(1)-J(2) model), employing exact diagonalization in both the S(z) = 0 basis and nearest neighbor singlet valence bond (NNVB) basis. We find that for 0.2 < J(2)/J(1) < 0.3, the NNVB basis gives a proper description of the ground state in comparison with the exact results. By analyzing the dimer-dimer as well as the plaquette-plaquette correlations and also defining appropriate structure factors, we investigate possible symmetry breaking states as the candidates for the ground state in the frustrated region. We provide a body of evidence in favor of plaquette valence bond ordering for 0.2 < J(2)/J(1) < 0.3. By further increasing the ratio J(2)/J(1), this state undergoes a transition to the staggered dimerized state.

arXiv:1007.0127•Journal of Physics Condensed Matter•May 19, 2011•76 citations
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 studied the ground-state phase diagram of spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model on the square lattice with first-neighbor J1{J}_{1}J1​ and second-neighbor J2{J}_{2}J2​ antiferromagnetic interactions using both infinite density matrix renormalization group (iDMRG) and DMRG approaches. We show that a plaquette valence bond phase is realized in an intermediate region $0.50\ensuremath{\le}{J}_{2}/{J}_{1}<0.54$ between a N\'eel magnetic ordered phase at J2/J1<0.50{J}_{2}/{J}_{1}<0.50J2​/J1​<0.50 and a stripy magnetic ordered phase at ${J}_{2}/{J}_{1}\ensuremath{\ge}0.54$. The plaquette valence bond phase is characterized by finite dimer orders in both the horizontal and vertical directions. Contrary to the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model, we do not find numerical evidence for a quantum spin liquid phase in the ${J}_{1}\text{\ensuremath{-}}{J}_{2} XY$ model.
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\times\else\texttimes\fi{}6\ifmmode\times\else\texttimes\fi{}6$. 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 ${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<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.
The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point J2c1=0.572(5)J1{J}_{2}^{{c}_{1}}=0.572(5){J}_{1}J2c1​​=0.572(5)J1​ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at J2/J1=0.6{J}_{2}/{J}_{1}=0.6J2​/J1​=0.6. These results are consistent with a potential gapless U(1)U(1)U(1) spin-liquid phase. However, since the U(1)U(1)U(1) spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at J2c1{J}_{2}^{{c}_{1}}J2c1​​. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.
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.
We study an S = 1/2 Heisenberg model on the honeycomb lattice with first and second neighbor antiferromagnetic exchange (J(1)-J(2) model), employing exact diagonalization in both the S(z) = 0 basis and nearest neighbor singlet valence bond (NNVB) basis. We find that for 0.2 < J(2)/J(1) < 0.3, the NNVB basis gives a proper description of the ground state in comparison with the exact results. By analyzing the dimer-dimer as well as the plaquette-plaquette correlations and also defining appropriate structure factors, we investigate possible symmetry breaking states as the candidates for the ground state in the frustrated region. We provide a body of evidence in favor of plaquette valence bond ordering for 0.2 < J(2)/J(1) < 0.3. By further increasing the ratio J(2)/J(1), this state undergoes a transition to the staggered dimerized state.