The exact solution of Kitaev's spin-
1/2 honeycomb spin-liquid model has sparked an intense search for Mott insulators hosting bond-dependent Kitaev interactions, of which
Na2IrO3 and $\ensuremath{\alpha}\text{\ensuremath{-}}\mathrm{Ru}{\mathrm{Cl}}_{3}$ are prime examples. Subsequently, it has been proposed that also spin-1 and spin-
3/2 analogs of Kitaev interactions may occur in materials with strong spin-orbit coupling. As a minimal model to describe these Kitaev materials, we study the Heisenberg-Kitaev Hamiltonian in a consistent
1/S expansion, with
S being the spin size. We present a comprehensive study of this model in the presence of an external magnetic field applied along two different directions, [001] and [111], for which an intricate classical phase diagram has been reported. In both settings, we employ spin-wave theory in a number of ordered phases to compute phase boundaries at the next-to-leading order in
1/S and show that quantum corrections substantially modify the classical phase diagram. More broadly, our work presents a consistent route to investigate the leading quantum corrections in spin models that break spin-rotational symmetry.
We study the extended Kitaev-Heisenberg (EKH) quantum spin model by adding bond-dependent off-diagonal Heisenberg term into the original KH model, which was recently proposed to describe the honeycomb Iridates. A rigorous mathematical mapping of spin operators reveals the intrinsic symmetry of the model Hamiltonian. By employing an unbiased numerical entanglement renormalization method based on tensor network ansatz, we obtain the global phase diagram containing eight distinct quantum phases. By using the dual mapping of spin operators, each of the individual magnetic phase in the global phase diagram can be clearly understood. At last, we show that a valence solid state emerges as the ground state in the quadro-critical region where multiple magnetic phases compete most intensively.
We derive and study a spin one-half Hamiltonian on a honeycomb lattice describing the exchange interactions between Ir4+ ions in a family of layered iridates A2IrO3 (A=Li,Na). Depending on the microscopic parameters, the Hamiltonian interpolates between the Heisenberg and exactly solvable Kitaev models. Exact diagonalization and a complementary spin-wave analysis reveal the presence of an extended spin-liquid phase near the Kitaev limit and a conventional Néel state close to the Heisenberg limit. The two phases are separated by an unusual stripy antiferromagnetic state, which is the exact ground state of the model at the midpoint between two limits.
We explore the phase diagram of the Kitaev-Heisenberg model with nearest neighbor interactions on the honeycomb lattice using the exact diagonalization of finite systems combined with the cluster mean field approximation, and supplemented by the insights from analytic approaches: the linear spin-wave and second-order perturbation theories. This study confirms that by varying the balance between the Heisenberg and Kitaev term, frustrated exchange interactions stabilize in this model either one of four phases with magnetic long range order: N\'eel phase, ferromagnetic phase, and two other phases with coexisting antiferromagnetic and ferromagnetic bonds, zigzag and stripy phase, or one of two distinct spin-liquid phases. Out of these latter disordered phases, the one with ferromagnetic Kitaev interactions has a substantially broader range of stability as the neighboring competing ordered phases, ferromagnetic and stripy, have very weak quantum fluctuations. Focusing on the quantum spin-liquid phases, we study spatial spin correlations and dynamic spin structure factor of the model by the exact diagonalization technique, and discuss the evolution of gapped low-energy spin response across the quantum phase transitions between the disordered spin liquid and phases with long range magnetic order.
The quest for Kitaev quantum spin liquids has led to great interest in honeycomb quantum magnets with strong spin-orbit coupling. It has been recently proposed that even Mott insulators with
3d transition metal ions, having nominally weak spin-orbit coupling, can realize such exotic physics. Motivated by this, we study the rhombohedral honeycomb cobaltates CoTiO
3, BaCo
2(PO
4)
2, and BaCo
2(AsO
4)
2, using
ab initio density functional theory, which takes into account realistic crystal field distortions and chemical information, in conjunction with exact diagonalization numerics. We show that these Co
2+ magnets host
jeff=1/2 local moments with highly anisotropic
g-factors, and we extract their full spin Hamiltonians including longer-range and anisotropic exchange couplings. For CoTiO
3, we find a nearest-neighbor easy-plane ferromagnetic
XXZ model with additional bond-dependent anisotropies and interlayer exchange, which supports three-dimensional (3D) Dirac nodal line magnons. In contrast, for BaCo
2(PO
4)
2 and BaCo
2(AsO
4)
2, we find a strongly suppressed interlayer coupling, and significant frustration from additional third-neighbor antiferromagnetic exchange mediated by P/As. Such bond-anisotropic
J1-
J3 spin models can support collinear zig-zag or coplanar spiral ground states; we discuss their dynamical spin correlations which reveal a gapped Goldstone mode, and argue that the effective parameters of these pseudospin-
1/2 models may be strongly renormalized by coupling to a low energy spin-exciton. Our results call for re-examining proposals for realizing Kitaev spin liquids in the honeycomb cobaltates.
In the quest for materials hosting Kitaev spin liquids, much of the efforts have been focused on the fourth- and fifth-row transition metal compounds, which are spin-orbit coupling--assisted Mott insulators. Here we study the structural and magnetic properties of
3d transition metal oxides,
Na2Co2TeO6 and
Na3Co2SbO6. The partial occupancy of sodium in former compound is addressed using a cluster expansion, and a honeycomb lattice of sodiums is found to be energetically favored. Starting from the ab initio band structures, a many-body second-order perturbation theory leads to a pseudospin-
21 Hamiltonian with estimated magnetic interactions. We show that the experimentally observed zigzag magnetic state is stabilized only when the first-neighbor Kitaev coupling dominates over the Heisenberg term, both of which are highly suppressed due to presence of
eg orbitals. A third-neighbor Heisenberg interaction is found dominant in both these compounds. We also present a phase diagram for
Na2Co2TeO6 by varying the electron-electron and spin-orbit interactions. The computed spin excitation spectra are found to capture essential features of recent experimental magnon spectrum, lending support to our results.