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

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

Miklós Lajkó, Karlo Penc

Published Jun 27, 2013• 21 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.87.224428

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

SU(4) Heisenberg model on the honeycomb lattice with exchange-frustrated perturbations: Implications for twistronics and Mott insulators

W. M. H. Natori, Rajah P. Nutakki, Rodrigo G. Pereira +1 more

The SU(4)-symmetric spin-orbital model on the honeycomb lattice was recently studied in connection to correlated insulators such as the ${e}_{g}$ Mott insulator ${\mathrm{Ba}}_{3}{\mathrm{CuSb}}_{2}{\mathrm{O}}_{9}$ and the insulating phase of magic-angle twisted bilayer graphene at quarter filling. Here we provide a unified discussion of these systems by investigating an extended model that includes the effects of Hund's coupling and anisotropic, orbital-dependent exchange interactions. Using a combination of mean-field theory, linear flavor-wave theory, and variational Monte Carlo, we show that this model harbors a quantum spin-orbital liquid over a wide parameter regime around the SU(4)-symmetric point. For large Hund's coupling, a ferromagnetic antiferro-orbital ordered state appears, while a valence-bond crystal combined with a vortex orbital state is stabilized by dominant orbital-dependent exchange interactions.

arXiv:1908.09224•Physical review. B./Physical review. B•Nov 21, 2019•43 citations

Gapped Topological Spin-Orbital Liquid on the Honeycomb Lattice

Masahiko G. Yamada

We perform large-scale density matrix renormalization group simulations of the SU(4) Heisenberg model on the honeycomb lattice to address the long-standing question of its ground state in an unbiased and quantitatively controlled manner. We find reliable numerical evidence that the ground state is a gapped spin-orbital liquid, presumably with a Z_{4} topological order, characterized by a finite topological entanglement entropy close to ln(4), the absence of both SU(4) and lattice symmetry breaking, and a variationally optimized ground-state energy well below the previously proposed π-flux variational state. By exploiting full SU(4) symmetry and keeping up to 12 800 SU(4) multiplets, corresponding to more than 1×10^{6} U(1) states, we achieve unprecedented accuracy for two-dimensional SU(4) quantum magnets. Finite-size scaling of energies and entanglement entropies supports a robust gapped phase in the two-dimensional limit, while a gapless critical state on narrow cylinders is identified as a proximate remnant of a Dirac spin-orbital liquid. Our results find the SU(4) honeycomb Heisenberg model a realization of a gapped topological spin-orbital liquid and provide convincing numerical evidence for topological order in a highly symmetric two-dimensional quantum magnet.

arXiv:2601.06549•Physical Review Letters•Sep 10, 2026•0 citations

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 +2 more

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.

arXiv:1302.1108•Physical Review B•May 10, 2013•41 citations

SU(4)-symmetric spin-orbital liquids on the hyperhoneycomb lattice

W. M. H. Natori, Eric C. Andrade, Rodrigo G. Pereira

We study the effective spin-orbital model that describes the magnetism of $4{d}^{1}$ or $5{d}^{1}$ Mott insulators in ideal tricoordinated lattices. In the limit of vanishing Hund's coupling, the model has an emergent SU(4) symmetry, which is made explicit by means of a Klein transformation on pseudospin degrees of freedom. Taking the hyperhoneycomb lattice as an example, we employ parton constructions with fermionic representations of the pseudospin operators to investigate possible quantum spin-orbital liquid states. We then use variational Monte Carlo (VMC) methods to compute the energies of the projected wave functions. Our numerical results show that the lowest-energy quantum liquid corresponds to a zero-flux state with a Fermi surface of four-color fermionic partons. In spite of the Fermi surface, we demonstrate that this state is stable against tetramerization. A combination of linear flavor wave theory and VMC applied to the complete microscopic model also shows that this liquid state is stable against the formation of collinear long-range order.

arXiv:1802.00044•Physical review. B./Physical review. B•Nov 13, 2018•33 citations

SU(4)-symmetric quantum spin-orbital liquids on various lattices

Masahiko G. Yamada, Masaki Oshikawa, George Jackeli

An emergent SU(4) symmetry discovered in the microscopic model for d1 honeycomb materials [M. G. Yamada, M. Oshikawa, and G. Jackeli, ] has enabled us to tailor exotic SU(4) models in real materials. In the honeycomb structure, the emergent SU(4) Heisenberg model would potentially have a quantum spin-orbital liquid ground state due to the , and we can also expect similar spin-orbital liquids in three-dimensional versions of the honeycomb lattice. In such quantum spin-orbital liquids, both the spin and orbital degrees of freedom become fractionalized and entangled together due to the strong frustrated interactions between them. Similarly to spinons in pure quantum spin liquids, quantum spin-orbital liquids can host not only spinon excitations, but also fermionic excitations at low temperature.

arXiv:2103.17076•Physical review. B./Physical review. B•Dec 30, 2021•17 citations
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π\ensuremath{\pi}\ensuremathπ-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.
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.
The SU(4)-symmetric spin-orbital model on the honeycomb lattice was recently studied in connection to correlated insulators such as the eg{e}_{g}eg​ Mott insulator Ba3CuSb2O9{\mathrm{Ba}}_{3}{\mathrm{CuSb}}_{2}{\mathrm{O}}_{9}Ba3​CuSb2​O9​ and the insulating phase of magic-angle twisted bilayer graphene at quarter filling. Here we provide a unified discussion of these systems by investigating an extended model that includes the effects of Hund's coupling and anisotropic, orbital-dependent exchange interactions. Using a combination of mean-field theory, linear flavor-wave theory, and variational Monte Carlo, we show that this model harbors a quantum spin-orbital liquid over a wide parameter regime around the SU(4)-symmetric point. For large Hund's coupling, a ferromagnetic antiferro-orbital ordered state appears, while a valence-bond crystal combined with a vortex orbital state is stabilized by dominant orbital-dependent exchange interactions.
We perform large-scale density matrix renormalization group simulations of the SU(4) Heisenberg model on the honeycomb lattice to address the long-standing question of its ground state in an unbiased and quantitatively controlled manner. We find reliable numerical evidence that the ground state is a gapped spin-orbital liquid, presumably with a Z_{4} topological order, characterized by a finite topological entanglement entropy close to ln(4), the absence of both SU(4) and lattice symmetry breaking, and a variationally optimized ground-state energy well below the previously proposed π-flux variational state. By exploiting full SU(4) symmetry and keeping up to 12 800 SU(4) multiplets, corresponding to more than 1×10^{6} U(1) states, we achieve unprecedented accuracy for two-dimensional SU(4) quantum magnets. Finite-size scaling of energies and entanglement entropies supports a robust gapped phase in the two-dimensional limit, while a gapless critical state on narrow cylinders is identified as a proximate remnant of a Dirac spin-orbital liquid. Our results find the SU(4) honeycomb Heisenberg model a realization of a gapped topological spin-orbital liquid and provide convincing numerical evidence for topological order in a highly symmetric two-dimensional quantum magnet.
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 study the effective spin-orbital model that describes the magnetism of 4d14{d}^{1}4d1 or 5d15{d}^{1}5d1 Mott insulators in ideal tricoordinated lattices. In the limit of vanishing Hund's coupling, the model has an emergent SU(4) symmetry, which is made explicit by means of a Klein transformation on pseudospin degrees of freedom. Taking the hyperhoneycomb lattice as an example, we employ parton constructions with fermionic representations of the pseudospin operators to investigate possible quantum spin-orbital liquid states. We then use variational Monte Carlo (VMC) methods to compute the energies of the projected wave functions. Our numerical results show that the lowest-energy quantum liquid corresponds to a zero-flux state with a Fermi surface of four-color fermionic partons. In spite of the Fermi surface, we demonstrate that this state is stable against tetramerization. A combination of linear flavor wave theory and VMC applied to the complete microscopic model also shows that this liquid state is stable against the formation of collinear long-range order.
An emergent SU(4) symmetry discovered in the microscopic model for d1 honeycomb materials [M. G. Yamada, M. Oshikawa, and G. Jackeli, ] has enabled us to tailor exotic SU(4) models in real materials. In the honeycomb structure, the emergent SU(4) Heisenberg model would potentially have a quantum spin-orbital liquid ground state due to the , and we can also expect similar spin-orbital liquids in three-dimensional versions of the honeycomb lattice. In such quantum spin-orbital liquids, both the spin and orbital degrees of freedom become fractionalized and entangled together due to the strong frustrated interactions between them. Similarly to spinons in pure quantum spin liquids, quantum spin-orbital liquids can host not only spinon excitations, but also fermionic excitations at low temperature.