Topological spin liquids in two spatial dimensions are stable phases in the presence of a small magnetic field, but may give way to field-induced phenomena at intermediate field strengths. Sandwiched between the low-field spin liquid physics and the
high-field spin-polarized phase, the exploration of magnetic phenomena in this intermediate regime however often remains elusive to controlled analytical approaches. Here we numerically study such intermediate-field magnetic phenomena for two representative Kitaev models (on the square-octagon and decorated honeycomb lattice) that exhibit either Abelian or non-Abelian topological order in the low-field limit. Using a combination of exact diagonalization and density matrix renormalization group techniques, as well as linear spin-wave theory, we establish the generic features of Kitaev spin liquids in an external magnetic field. While ferromagnetic models typically exhibit a direct transition to the polarized state at a relatively low field strength, antiferromagnetic couplings not only substantially stabilizes the topological spin liquid phase, but generically lead to the emergence of a distinct field-induced intermediate regime, separated by a crossover from the high-field polarized regime. Our results suggest that, for most lattice geometries, this regime generically exhibits significant spin canting, antiferromagnetic spin-spin correlations, and an extended proximate spin liquid regime at finite temperatures. Notably, we identify a symmetry obstruction in the original honeycomb Kitaev model that prevents, at least for certain field directions, the formation of such canted magnetism without breaking symmetries -- consistent with the recent numerical observation of an extended gapless spin liquid in this case.
The appearance of nontrivial phases in Kitaev materials exposed to an external magnetic field has recently been a subject of intensive studies. Here, we elucidate the relation between the field-induced ground states of the classical and quantum spin
models proposed for such materials, by using the infinite density matrix renormalization group (iDMRG) and the linear spin wave theory (LSWT). We consider the $K Gamma Gamma$ model, where $Gamma$ and $Gamma$ are off-diagonal spin exchanges on top of the dominant Kitaev interaction $K$. Focusing on the magnetic field along the $[111]$ direction, we explain the origin of the nematic paramagnet, which breaks the lattice-rotational symmetry and exists in an extended window of magnetic field, in the quantum model. This phenomenon can be understood as the effect of quantum order-by-disorder in the frustrated ferromagnet with a continuous manifold of degenerate ground states discovered in the corresponding classical model. We compute the dynamical spin structure factors using a matrix operator based time evolution and compare them with the predictions from LSWT. We, thus, provide predictions for future inelastic neutron scattering experiments on Kitaev materials in an external magnetic field along the $[111]$ direction. In particular, the nematic paramagnet exhibits a characteristic pseudo-Goldstone mode which results from the lifting of a continuous degeneracy via quantum fluctuations.
Motivated by recent experiments on $alpha$-RuCl$_3$, we investigate a possible quantum spin liquid ground state of the honeycomb-lattice spin model with bond-dependent interactions. We consider the $K-Gamma$ model, where $K$ and $Gamma$ represent the
Kitaev and symmetric-anisotropic interactions between spin-1/2 moments on the honeycomb lattice. Using the infinite density matrix renormalization group (iDMRG), we provide compelling evidence for the existence of quantum spin liquid phases in an extended region of the phase diagram. In particular, we use transfer matrix spectra to show the evolution of two-particle excitations with well-defined two-dimensional dispersion, which is a strong signature of quantum spin liquid. These results are compared with predictions from Majorana mean-field theory and used to infer the quasiparticle excitation spectra. Further, we compute the dynamical structure factor using finite size cluster computations and show that the results resemble the scattering continuum seen in neutron scattering experiments on $alpha$-RuCl$_3$. We discuss these results in light of recent and future experiments.
