Previous investigations have suggested that the simplest spin-orbital model on the simplest frustrated lattice can host a nematic quantum spin-orbital liquid state. Namely, the orbital degeneracy of the SU(4) Kugel-Khomskii model tends to enhance quantum fluctuations and stabilize a quantum spin-orbital liquid exhibiting stripy features on the triangular lattice, as revealed by the state-of-the-art method of the density matrix renormalization group boosted by Gutzwiller projected wave functions. In this work, using the variational quantum Monte Carlo method, we study several spin-orbital liquid states, including a uniform $\ensuremath{\pi}$-flux state, three stripy states, and a plaquette state, on the $L\ifmmode\times\else\texttimes\fi{}L$ torus up to
L=24. It turns out that one of these stripy states, called the ``stripe-II'' state, is energetically favored. This ground state breaks the
C6 symmetry of the lattice, resulting in a reduced
C2 symmetry and doubled unit cells while preserving the SU(4) spin-orbital rotation symmetry. Such a nematic quantum spin-orbital liquid state can be characterized by a parton Fermi surface (FS) consisting of open orbits in the Brillouin zone, in contrast to the circular FS of the uniform $\ensuremath{\pi}$-flux state.
Previous investigations have suggested that the simplest spin-orbital model on the simplest frustrated lattice can host a nematic quantum spin-orbital liquid state. Namely, the orbital degeneracy of the SU(4) Kugel-Khomskii (KK) model tends to enhance quantum fluctuations and stabilize a quantum spin-orbital liquid exhibiting stripy features on the triangular lattice, as revealed by the state-of-the-art method of the density matrix renormalization group boosted by Gutzwiller projected wave functions. In this work, using the variational quantum Monte Carlo method, we have studied several spin-orbital liquid states, including a uniform
π flux state, three stripy states, and a plaquette state, on the
L×L torus up to
L=24. It turns out that one of these stripy states, called the "stripe-II" state, is energetically favored. This ground state breaks the
C6 symmetry of the lattice, resulting in a reduced
C2 symmetry and doubled unit cells, while preserving the SU(4) spin-orbital rotation symmetry. Such a nematic quantum spin-orbital liquid state can be characterized by a parton Fermi surface (FS) consisting of open orbits in the Brillouin zone, in contrast to the circular FS of the uniform
π-flux state.
By combining the density matrix renormalization group (DMRG) method with Gutzwiller projected wave functions, we study the SU(4) symmetric spin-orbital model on the honeycomb lattice. We find that the ground states can be well described by a Gutzwiller projected $\ensuremath{\pi}$-flux state with Dirac-type gapless excitations at one quarter filling. Although these Dirac points are gapped by emergent gauge fluxes on finite cylinders, they govern the critical behavior in the thermodynamic limit. By inserting a $\ensuremath{\theta}=\ensuremath{\pi}$ spin flux to twist the boundary condition, we can shift the gapless sector to the ground state, which provides compelling evidence for the presence of a gapless Dirac spin-orbital liquid.
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.
Dirac spin liquids represent a class of highly-entangled quantum phases in two dimensional Mott insulators, featuring exotic properties such as critical correlation functions and absence of well-defined low energy quasi-particles. Existing numerical works suggest that the spin-orbital
SU(4) symmetric Kugel-Khomskii model of Mott insulators on the honeycomb lattice realizes a Dirac spin-orbital liquid, described at low energy by
(2+1)d quantum electrodynamics (QED
3) with
Nf=8 Dirac fermions. We generalize methods previously developed for
SU(2) spin systems to analyze the symmetry properties and stability of the Dirac spin-orbital liquid. We conclude that the standard Dirac state in the
SU(4) honeycomb system, based on a simple parton mean-field ansatz, is intrinsically unstable at low energy due to the existence of a monopole perturbation that is allowed by physical symmetries and relevant under renormalization group flow. We propose two plausible alternative scenarios compatible with existing numerics. In the first scenario, we propose an alternative
U(1) Dirac spin-orbital liquid, which is similar to the standard one except for its monopole symmetry quantum numbers. This alternative
U(1) state represents a stable gapless phase. In the second scenario, we start from the standard
U(1) Dirac liquid and Higgs the
U(1) gauge symmetry down to
Z4. The resulting
Z4 Dirac spin-orbital liquid is stable. We also discuss the continuous quantum phase transitions from the
Z4 Dirac liquids to conventional symmetry-breaking orders, described by the QED
3 theory with
Nf=8 supplemented with a critical charge-
4 Higgs field. We discuss possible ways to distinguish these scenarios in numerics. We also extend previous calculations of the quantum anomalies of QED
3 and match with generalized lattice Lieb-Schultz-Mattis constraints.