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.
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.
The
SU(4) Heisenberg model on the honeycomb lattice is expected to host a quantum spin-orbital liquid at low temperature with an astonishing candidate material,
α-ZrCl
3. We employed the canonical thermal pure quantum state method to investigate the finite-temperature phase of this model. Exploiting the full symmetry of
SU(4), the calculation up to a 24-site cluster, which is equivalent to 48 sites in the spin-1/2 language, is possible. This state-of-the-art computation with large-scale parallelization enables us to capture the thermodynamic properties of the
SU(4) Heisenberg model on the honeycomb lattice. In particular, the specific heat shows a characteristic peak-and-shoulder structure, which should be related to the nature of the low-temperature quantum spin-orbital liquid phase. We also discuss what can be concluded from the assumption that the ground state is gapped and symmetric in view of the generalized Lieb-Schultz-Mattis theorem.
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.
While the enhancement of spin-space symmetry from the usual SU(2) to SU(N) is promising for finding nontrivial quantum spin liquids, its realization in magnetic materials remains challenging. Here, we propose a new mechanism by which SU(4) symmetry emerges in the strong spin-orbit coupling limit. In d^{1} transition metal compounds with edge-sharing anion octahedra, the spin-orbit coupling gives rise to strongly bond-dependent and apparently SU(4)-breaking hopping between the J_{eff}=3/2 quartets. However, in the honeycomb structure, a gauge transformation maps the system to an SU(4)-symmetric Hubbard model. In the strong repulsion limit at quarter filling, as realized in α-ZrCl_{3}, the low-energy effective model is the SU(4) Heisenberg model on the honeycomb lattice, which cannot have a trivial gapped ground state and is expected to host a gapless spin-orbital liquid. By generalizing this model to other three-dimensional lattices, we also propose crystalline spin-orbital liquids protected by this emergent SU(4) symmetry and space group symmetries.
We study a
Z3 Kitaev model on the honeycomb lattice with nearest neighbor interactions. Based on matrix product state simulations and symmetry considerations, we find evidence that, with ferromagnetic isotropic couplings, the model realizes a chiral spin liquid, characterized by a possible
U(1)12 chiral topological order. This is supported by simulations on both cylinder and strip geometries. On infinitely long cylinders with various widths, scaling analysis of entanglement entropy and maximal correlation length suggests that the model has a gapped two-dimensional bulk. The topological entanglement entropy is extracted and found to be in agreement with the
U(1)12 topological order. On infinitely long strips with moderate widths, we find the model is critical with a central charge consistent with the chiral edge theory of the
U(1)12 topological phase. We conclude by discussing several open questions.