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

Competing states in the SU(3) Heisenberg model on the honeycomb lattice: Plaquette valence-bond crystal versus dimerized color-ordered state

Philippe Corboz, Miklós Lajkó, Karlo Penc, Frédéric Mila, Andreas Martin Läuchli

Published May 10, 2013• 41 citations
arXivHTMLPDF

Abstract

Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al., Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory suggested a dimerized, color-ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension $D$. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small $D$ the dimer state has a lower variational energy than the plaquette state; however, for large $D$ it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are nonchiral flux states.

DOI
10.1103/physrevb.87.195113

Related papers

Plaquette order in the SU(6) Heisenberg model on the honeycomb lattice

Pierre Nataf, Miklós Lajkó, Philippe Corboz +3 more

We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611(R) (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair state simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence that the model with one particle per site and nearest-neighbor exchange actually develops plaquette order. This is further confirmed by the investigation of the model with a ring-exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring-exchange term.

arXiv:1601.00959•Physical review. B./Physical review. B•May 26, 2016•39 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

SU(3) fermions on the honeycomb lattice at 13 filling

Sangwoo S. Chung, Philippe Corboz

$\mathrm{SU}(N)$ symmetric fermions on a lattice, which can be realized in ultracold-atom-based quantum simulators, have very promising prospects for realizing exotic states of matter. Here we present the ground-state phase diagram of the repulsive SU(3) Hubbard model on a honeycomb lattice at 1/3 filling obtained from infinite projected entangled pair states tensor network calculations. In the strongly interacting limit the ground state has plaquette order. Upon decreasing the interaction strength $U/t$ we find a first-order transition at $U/t=7.2(2)$ into a dimerized, color-ordered state, which extends down to $U/t=4.5(5)$, at which the Mott transition occurs and the ground state becomes uniform. Our results may serve as a prediction and benchmark for future quantum simulators of SU(3) fermions.

arXiv:1902.10478•Physical review. B./Physical review. B•Jul 29, 2019•26 citations

Tetramerization in a SU(4) Heisenberg model on the honeycomb lattice

Miklós Lajkó, Karlo Penc

The SU(4) Heisenberg model can serve as a low-energy model of the Mott insulating state in materials where the spins and orbitals are highly symmetric or in systems of alkaline-earth atoms on an optical lattice. Recently, it has been argued that on the honeycomb lattice, the model exhibits a unique spin-orbital liquid phase with an algebraic decay of correlations [P. Corboz et al., Phys. Rev. X 2, 041013 (2012)]. Here we study the instability of the algebraic spin-orbital liquid toward the spontaneous formation of SU(4) singlet plaquettes (tetramerization). Using a variational Monte Carlo approach to evaluate the projected wave function of fermions with $\ensuremath{\pi}$-flux state, we find that the algebraic liquid is robust and that a small finite value of the next-nearest exchange is needed to induce tetramerization. We also studied the phase diagram of a model which interpolates between the nearest-neighbor Heisenberg model and a Hamiltonian for which the singlet-plaquette-product state is an exact ground state.

arXiv:1303.6297•Physical Review B•Jun 27, 2013•21 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

Symmetry breaking and competing valence bond states in the star lattice Heisenberg antiferromagnet

Pratyay Ghosh, Jan Alexander Koziol, Samuel Nyckees +2 more

We investigate the ground state phase diagram of the spin-$1/2$ antiferromagnetic Heisenberg model on the star lattice using infinite projected entangled pair states (iPEPS) and high-order series expansions. The model includes two distinct couplings: ${J}_{d}$ on the dimer bonds and ${J}_{t}$ on the trimer bonds. While it is established that the system hosts a valence bond solid (VBS) phase for ${J}_{d}\ensuremath{\ge}{J}_{t}$, the ground state phase diagram for ${J}_{d}<{J}_{t}$ has remained unsettled. Our iPEPS simulations uncover a first-order phase transition at ${J}_{d}/{J}_{t}\ensuremath{\approx}0.18$, significantly lower than previously reported estimates. Beyond this transition, we identify a close competition between two valence bond crystal (VBC) states; a columnar VBC and a $\sqrt{3}\ifmmode\times\else\texttimes\fi{}\sqrt{3}$ VBC, with the latter consistently exhibiting lower energy across all finite bond dimensions. The high-order series expansion supports this by finding that the $\sqrt{3}\ifmmode\times\else\texttimes\fi{}\sqrt{3}$ VBC state indeed becomes energetically favorable, but only at sixth order in perturbation theory, revealing the subtle nature of the competition between candidate states.

arXiv:2508.05133•Physical review. B./Physical review. B•Sep 24, 2025•3 citations
Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al., Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory suggested a dimerized, color-ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension DDD. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small DDD the dimer state has a lower variational energy than the plaquette state; however, for large DDD it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are nonchiral flux states.
We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611(R) (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair state simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence that the model with one particle per site and nearest-neighbor exchange actually develops plaquette order. This is further confirmed by the investigation of the model with a ring-exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring-exchange term.
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.
SU(N)\mathrm{SU}(N)SU(N) symmetric fermions on a lattice, which can be realized in ultracold-atom-based quantum simulators, have very promising prospects for realizing exotic states of matter. Here we present the ground-state phase diagram of the repulsive SU(3) Hubbard model on a honeycomb lattice at 1/3 filling obtained from infinite projected entangled pair states tensor network calculations. In the strongly interacting limit the ground state has plaquette order. Upon decreasing the interaction strength U/tU/tU/t we find a first-order transition at U/t=7.2(2)U/t=7.2(2)U/t=7.2(2) into a dimerized, color-ordered state, which extends down to U/t=4.5(5)U/t=4.5(5)U/t=4.5(5), at which the Mott transition occurs and the ground state becomes uniform. Our results may serve as a prediction and benchmark for future quantum simulators of SU(3) fermions.
The SU(4) Heisenberg model can serve as a low-energy model of the Mott insulating state in materials where the spins and orbitals are highly symmetric or in systems of alkaline-earth atoms on an optical lattice. Recently, it has been argued that on the honeycomb lattice, the model exhibits a unique spin-orbital liquid phase with an algebraic decay of correlations [P. Corboz et al., Phys. Rev. X 2, 041013 (2012)]. Here we study the instability of the algebraic spin-orbital liquid toward the spontaneous formation of SU(4) singlet plaquettes (tetramerization). Using a variational Monte Carlo approach to evaluate the projected wave function of fermions with $\ensuremath{\pi}$-flux state, we find that the algebraic liquid is robust and that a small finite value of the next-nearest exchange is needed to induce tetramerization. We also studied the phase diagram of a model which interpolates between the nearest-neighbor Heisenberg model and a Hamiltonian for which the singlet-plaquette-product state is an exact ground state.
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 investigate the ground state phase diagram of the spin-1/21/21/2 antiferromagnetic Heisenberg model on the star lattice using infinite projected entangled pair states (iPEPS) and high-order series expansions. The model includes two distinct couplings: Jd{J}_{d}Jd​ on the dimer bonds and Jt{J}_{t}Jt​ on the trimer bonds. While it is established that the system hosts a valence bond solid (VBS) phase for ${J}_{d}\ensuremath{\ge}{J}_{t}$, the ground state phase diagram for Jd<Jt{J}_{d}<{J}_{t}Jd​<Jt​ has remained unsettled. Our iPEPS simulations uncover a first-order phase transition at ${J}_{d}/{J}_{t}\ensuremath{\approx}0.18$, significantly lower than previously reported estimates. Beyond this transition, we identify a close competition between two valence bond crystal (VBC) states; a columnar VBC and a $\sqrt{3}\ifmmode\times\else\texttimes\fi{}\sqrt{3}$ VBC, with the latter consistently exhibiting lower energy across all finite bond dimensions. The high-order series expansion supports this by finding that the $\sqrt{3}\ifmmode\times\else\texttimes\fi{}\sqrt{3}$ VBC state indeed becomes energetically favorable, but only at sixth order in perturbation theory, revealing the subtle nature of the competition between candidate states.