Abstract The pursuit of the quantum spin liquid (QSL) predicted by the Kitaev honeycomb model remains a central yet challenging goal in quantum magnetism. Realistic candidate materials inevitably host non-Kitaev interactions that often lead to antiferromagnetic (AFM) order at low temperature. While magnetic fields can suppress such AFM order and may stabilize a QSL phase, conclusive evidence for a field-induced Kitaev QSL has remained elusive despite a decade of research. Here, we establish an effective K - J -Γ-
Γ′ Γ ′ model with a dominant AFM Kitaev interaction for Na 2 Co 2 TeO 6 , which quantitatively explains its key experimental measurements. Using high-precision tensor-network calculations, we reveal a QSL phase under intermediate [111] magnetic fields, which is possibly gapless, and find it can be adiabatically connected to the intensively studied intermediate-field QSL in the pure AFM Kitaev model. This correspondence confirms that the Kitaev model as the origin of the intermediate-field QSL in cobalt-based magnet Na 2 Co 2 TeO 6 , thus offering a concrete platform for exploring Kitaev-derived QSL in realistic materials.
The possible existence of a quantum spin liquid (QSL) phase, an exotic state of matter with long-range quantum entanglement and fractionalized excitations, in
α-RuCl
3 has sparked widespread interests in exploring QSLs in various Kitaev models under magnetic fields. Recently, a
K-
J-
Γ-
Γ′ model has been proposed to accurately describe the compound by fitting the measured thermodynamic data, where
J is Heisenberg interaction, and
Γ,
Γ′ are off-diagonal exchanges on top of the dominant Kitaev coupling
K. Based on this effective model, an intermediate QSL phase in presence of an out-of-plane fields along the
[111] direction, between the low-field zigzag order and high-field polarized phase, has been predicted. By combining density matrix renormalization group (DMRG), exponential tensor renormalization group (XTRG), and variational Monte Carlo (VMC) calculations, we address the nature of this QSL phase in the honeycomb
K-
J-
Γ-
Γ′ model under the
[111]-direction field. Our DMRG calculations find the algebraic-like decay of spin correlation function, the finite spin scalar chiral order, and lattice nematic order. Together with the XTRG results of power-law specific heat at low temperature, our findings naturally suggest a gapless nematic chiral spin liquid. On the other hand, our VMC study finds a gapped nematic chiral spin liquid with the variational energy very close to that obtained in DMRG. As VMC finds a very small gap that is beyond the current resolution of both ground-state DMRG and finite-temperature XTRG calculations on finite-width cylinders, we resort the full clarification for the nature of the intermediate QSL to future studies. Lastly, we discuss the implications of our results to the recent experiment on
α-RuCl
3 and the QSL-like phase in a similar
K-
Γ-
Γ′ model.
Motivated by the three-dimensional structure of Kitaev materials we explore multilayer Kitaev models. The magnetic properties of a multilayer of an arbitrary number of Kitaev honeycomb layers stacked on top of each other coupled through a Heisenberg interaction, J, is analyzed through Abrikosov fermion mean-field theory. The system sustains quantum spin liquid (QSL) solutions which have different character depending on parity of the number of layers. While in even layered Kitaev models a gapped QSL emerges, odd-layered models host gapless QSLs. The projective symmetry group analysis of these solutions unravel a layer-to-layer inversion symmetry rather than an expected reflection. Although these QSLs retain features of the single layer Kitaev spin liquid (KSL), they should be regarded hybrid QSLs consisting on several KSLs. The good agreement at large-J between the energy of the Gutzwiller projected mean-field QSL and the exact energy indicates that such QSL Ansatz is adiabatically connected to the exact ground state. We also find that the Kitaev gapped chiral quantum spin liquid induced by external magnetic fields is stabilized by an antiferromagnetic interlayer coupling. Our results are relevant to the physics of alpha-RuCl3 and H3LiIr2O6 which are examples of magnetically coupled multilayer Kitaev models.