3D Spin-orbital liquids

Spin-orbital liquids provide an exactly solvable route to three-dimensional Z2 quantum spin liquids beyond the original Kitaev setting. Built from higher-dimensional Clifford-algebra representations, spin-orbital Hamiltonians can be realized on both three- and four-coordinated lattices, giving rise to phases with 3 and 2 itinerant Majorana flavors. We demonstrate that these models host a rich set of gapless Majorana metals, characterized, in particular, by topological Fermi surfaces, nodal lines, and Weyl semimetal phases. We analyze the stability of these structures under physically motivated perturbations and identify generic splitting patterns and topological transitions driven by symmetry breaking and flavor mixing. This yields a unified organizing framework for three-dimensional Majorana metals in fractionalized spin liquids.

💡 Research Summary

This paper introduces a broad class of exactly solvable three‑dimensional (3D) quantum spin liquids (QSLs) built from spin‑orbital liquids (SOLs) that exploit higher‑dimensional representations of the Clifford algebra. By replacing the Pauli matrices of the original Kitaev honeycomb model with a set of five Γ‑matrices (the q = 1 case), the authors enlarge the on‑site Hilbert space to four dimensions, naturally interpreted as a combined spin‑orbital degree of freedom. Each Γ‑matrix is expressed as a tensor product of ordinary spin (σ) and orbital (τ) Pauli matrices, allowing a systematic construction of Hamiltonians on any lattice with coordination number γₘ.

The key structural insight is that the number of itinerant Majorana flavors ν is tied to the lattice coordination: ν = 6 − γₘ. Consequently, three‑coordinated lattices host ν = 3 identical itinerant Majoranas, while four‑coordinated lattices host ν = 2. The remaining γₘ Majoranas become static Z₂ gauge fields u_{ij}=±1 defined on bonds. After fixing a gauge sector (typically the Lieb flux sector, determined by mirror symmetries or numerical minimization), the problem reduces to free Majorana fermions hopping in a static Z₂ background.

For three‑coordinated lattices the spin part of the interaction is isotropic (σ·σ)⊗τ_iτ_j, giving an emergent SO(3) symmetry that guarantees the three flavors share identical band structures. For four‑coordinated lattices the interaction reduces to (σ_xσ_x+σ_yσ_y)⊗τ_iτ_j, leaving only an SO(2) symmetry and two distinct Majorana bands. The authors explicitly write the Hamiltonian in terms of six Majoranas per site (b₁…b₅, c) and show how the bond operators \hat u_{ij}=i b_i b_j commute with the Hamiltonian, thus acting as static Z₂ gauge fields.

The paper then analyses the gauge flux sector using Wilson loop operators \hat W_p = (−i)^N ∏{⟨jj+1⟩∈p} \hat u{jj+1}. For lattices possessing mirror planes (e.g., hyperhexagon, layered honeycomb) Lieb’s theorem predicts that loops of length 6 or 10 carry zero flux while length‑8 loops carry π‑flux. Numerical checks confirm that this “Lieb flux sector” indeed yields the lowest Majorana ground‑state energy for all lattices considered, including the four‑coordinated ones where no analytical theorem applies.

With the gauge fixed, the authors compute the Majorana band structures and classify the resulting gapless phases. For ν = 3 they find threefold Fermi surfaces, intertwined nodal‑line networks, and hybrid structures where nodal lines intersect. For ν = 2 they uncover Weyl semimetal phases with pairs of Weyl points carrying opposite Berry monopole charge; breaking time‑reversal symmetry gaps the Weyl points and produces Chern insulating states. The topological invariants (Chern numbers, Weyl charges) are directly linked to the underlying lattice symmetries (time‑reversal, inversion, rotation).

A major part of the work is devoted to the stability of these gapless phases under physically motivated perturbations. The authors consider (i) spin‑orbital anisotropies that reduce the emergent SO(3) symmetry to SO(2), (ii) external magnetic fields that break time‑reversal, and (iii) bond‑direction dependent exchange variations that mix the Majorana flavors. They show that symmetry reduction typically splits a three‑flavor Fermi surface into smaller pockets or converts it into nodal lines, while flavor mixing can shift Weyl points in momentum space or cause pair annihilation, leading to topological phase transitions. The transitions are tracked by monitoring changes in Berry curvature and the evolution of Wilson loop eigenvalues, providing a clear diagnostic for experimental probes.

Thermal stability is also discussed. Because the Z₂ gauge field is static in the exactly solvable limit, vison excitations (π‑fluxes) are gapped, and the low‑temperature phase remains a deconfined Z₂ spin liquid. However, strong perturbations or pressure could proliferate visons, potentially driving the system into a confined phase or into a different flux sector with altered Majorana band topology. The authors suggest that such transitions could be detected via NMR line broadening, Raman scattering signatures, or non‑linear optical responses that are sensitive to the presence of gapless Majorana excitations.

In summary, the paper provides a unified theoretical framework for 3D spin‑orbital liquids with multiple Majorana flavors, systematically classifies the resulting gapless Majorana metals (Fermi‑surface, nodal‑line, Weyl) across a variety of three‑ and four‑coordinated lattices, and elucidates how symmetry breaking and flavor mixing drive topological transitions. This work significantly expands the landscape of exactly solvable 3D QSL models and offers concrete guidance for future material realizations and experimental detection of exotic Majorana‑based quantum phases.