We study the interplay of disorder and Heisenberg interactions in the Kitaev model on a honeycomb lattice. The effect of disorder on the transition between Kitaev spin liquid and magnetic ordered states as well as the stability of magnetic ordering is investigated. Using Lanczos exact diagonalization we discuss the consequences of two types of disorder: (i) random-coupling disorder and (ii) singular-coupling disorder. They exhibit qualitatively similar effects in the pure Kitaev-Heisenberg model without long-range interactions. The range of spin-liquid phases is reduced and the transition to magnetic ordered phases becomes more crossoverlike. Furthermore, the long-range zigzag and stripy orderings in the clean system are replaced by their three domains with different ordering direction. Especially in the crossover range the coexistence of magnetically ordered and Kitaev spin-liquid domains is possible. With increasing the disorder strength the area of domains becomes smaller and the system goes into a spin-glass state. However, the disorder effect is different in magnetically ordered phases caused by long-range interactions. The stability of such magnetic ordering is diminished by singular-coupling disorder and, accordingly, the range of the spin-liquid regime is extended. This mechanism may be relevant to materials like $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ and ${\mathrm{H}}_{3}{\mathrm{LiIr}}_{2}{\mathrm{O}}_{6}$ where the zigzag ground state is stabilized by weak long-range interactions. We also find that the flux gap closes at a critical disorder strength and vortices appears in the flux arrangement. Interestingly, the vortices tend to form kinds of commensurate ordering.
We study the interplay of disorder and Heisenberg interactions in the Kitaev model on a honeycomb lattice. The effect of disorder on the transition between Kitaev spin liquid and magnetic ordered states as well as the stability of magnetic ordering is investigated. Using Lanczos exact diagonalization we discuss the consequences of two types of disorder: (i) random-coupling disorder and (ii) singular-coupling disorder. They exhibit qualitatively similar effects in the pure Kitaev-Heisenberg model without long-range interactions. The range of spin-liquid phases is reduced and the transition to magnetic ordered phases becomes more crossoverlike. Furthermore, the long-range zigzag and stripy orderings in the clean system are replaced by their three domains with different ordering direction. Especially in the crossover range the coexistence of magnetically ordered and Kitaev spin-liquid domains is possible. With increasing the disorder strength the area of domains becomes smaller and the system goes into a spin-glass state. However, the disorder effect is different in magnetically ordered phases caused by long-range interactions. The stability of such magnetic ordering is diminished by singular-coupling disorder and, accordingly, the range of the spin-liquid regime is extended. This mechanism may be relevant to materials like \ensuremathα\ensuremath-RuCl3 and H3LiIr2O6 where the zigzag ground state is stabilized by weak long-range interactions. We also find that the flux gap closes at a critical disorder strength and vortices appears in the flux arrangement. Interestingly, the vortices tend to form kinds of commensurate ordering.
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.
We investigate the properties of the Kitaev honeycomb model with both site dilution
and exchange randomness. Embarking on this work, we review disorder in some strongly correlated electron systems, including spin-½ and spin-1 Heisenberg antiferromagnetic chains, two dimensional Heisenberg antiferromagnets, the cuprates and graphene. We outline some aspects of resonating valence bond phases, valence bond solids, spin liquids and quantum computation that are pertinent to an understanding of the Kitaev model. The properties of the Kitaev model without disorder are discussed and it is found to be a critical spin liquid, with algebraic correlations in two spin operators sigma^{alpha}_{i}sigma^{alpha}_{j}, where i and j,/em> are either end of a link of type alpha = x, y or z on the honeycomb lattice.
The Kitaev model is exactly solvable and we show that this remains so in the presence of site dilution and exchange randomness. We find that vacancies bind a flux. In the gapped phase, a vacancy forms an effective paramagnetic moment. With two or more vacancies we describe the interaction of their effective moments and show that a finite density of vacancies leads to a divergent macroscopic susceptibility at small fields.
In the gapless phase the effective moment has a susceptibility that is, to leading order at small fields, chi(h)~log(1/h). Interaction between the moments from two vacancies on opposite sublattices cuts off this divergence in susceptibility at a large but finite constant. Two vacancies on the same sublattice behave quite differently and we find the combined susceptibility is parametrically larger than that of an isolated vacancy, chi(h)sim [h(log(1/h))^{3/2}]^{-1}.
We also investigate the effects of slowly varying, quenched disorder in exchange coupling. We demonstrate that this does not qualitatively affect the susceptibility but show that the heat capacity C ~ T^{2/z}, where z is a measure of the disorder and increases from one with increasing disorder strength.
We investigate the properties of the Kitaev honeycomb model with both site dilution and exchange randomness. Embarking on this work, we review disorder in some strongly correlated electron systems, including spin-½ and spin-1 Heisenberg antiferromagnetic chains, two dimensional Heisenberg antiferromagnets, the cuprates and graphene. We outline some aspects of resonating valence bond phases, valence bond solids, spin liquids and quantum computation that are pertinent to an understanding of the Kitaev model. The properties of the Kitaev model without disorder are discussed and it is found to be a critical spin liquid, with algebraic correlations in two spin operators sigma^{alpha}_{i}sigma^{alpha}_{j}, where i and j are either end of a link of type alpha = x, y or z on the honeycomb lattice. The Kitaev model is exactly solvable and we show that this remains so in the presence of site dilution and exchange randomness. We find that vacancies bind a flux. In the gapped phase, a vacancy forms an effective paramagnetic moment. With two or more vacancies we describe the interaction of their effective moments and show that a finite density of vacancies leads to a divergent macroscopic susceptibility at small fields. In the gapless phase the effective moment has a susceptibility that is, to leading order at small fields, chi(h)~log(1/h). Interaction between the moments from two vacancies on opposite sublattices cuts off this divergence in susceptibility at a large but finite constant. Two vacancies on the same sublattice behave quite differently and we find the combined susceptibility is parametrically larger than that of an isolated vacancy, chi(h)sim [h(log(1/h))^{3/2}]^{-1}. We also investigate the effects of slowly varying, quenched disorder in exchange coupling. We demonstrate that this does not qualitatively affect the susceptibility but show that the heat capacity C ~ T^{2/z}, where z is a measure of the disorder and increases from one with increasing disorder strength.
Abstract We study the ground state of the 1D Kitaev-Heisenberg (KH) model using the density-matrix renormalization group and Lanczos exact diagonalization methods. We obtain a rich ground-state phase diagram as a function of the ratio between Heisenberg (J = cosϕ) and Kitaev (K = sinϕ) interactions. Depending on the ratio, the system exhibits four long-range ordered states: ferromagnetic-z, ferromagnetic-xy, staggered-xy, Néel-z, and two liquid states: Tomonaga-Luttinger liquid and spiral-xy. The two Kitaev points ϕ=2π ϕ=π2 and φ=23π φ=3π2 are singular. Theϕ-dependent phase diagram is similar to that for the 2D honeycomb-lattice KH model. Remarkably, all the ordered states of the honeycomb-lattice KH model can be interpreted in terms of the coupled KH chains. We also discuss the magnetic structure of the K-intercalated RuCl3, a potential Kitaev material, in the framework of the 1D KH model. Furthermore, we demonstrate that the low-lying excitations of the 1D KH Hamiltonian can be explained within the combination of the known six-vertex model and spin-wave theory.
We study the ground state of the 1D Kitaev-Heisenberg (KH) model using the density-matrix renormalization group and Lanczos exact diagonalization methods. We obtain a rich ground-state phase diagram as a function of the ratio between Heisenberg (J = cosϕ) and Kitaev (K = sinϕ) interactions. Depending on the ratio, the system exhibits four long-range ordered states: ferromagnetic-z, ferromagnetic-xy, staggered-xy, Néel-z, and two liquid states: Tomonaga-Luttinger liquid and spiral-xy. The two Kitaev points ϕ=π2 and φ=3π2 are singular. The ϕ-dependent phase diagram is similar to that for the 2D honeycomb-lattice KH model. Remarkably, all the ordered states of the honeycomb-lattice KH model can be interpreted in terms of the coupled KH chains. We also discuss the magnetic structure of the K-intercalated RuCl3, a potential Kitaev material, in the framework of the 1D KH model. Furthermore, we demonstrate that the low-lying excitations of the 1D KH Hamiltonian can be explained within the combination of the known six-vertex model and spin-wave theory.
We study quantum disordered ground states of the two-dimensional Heisenberg-Kitaev model on the triangular lattice using a Schwinger boson approach. Our aim is to identify and characterize potential gapped quantum spin liquid phases that are stabilized by anisotropic Kitaev interactions. For antiferromagnetic Heisenberg and Kitaev couplings and sufficiently small spin S, we find three different symmetric Z2 spin liquid phases, separated by two continuous quantum phase transitions. Interestingly, the gap of elementary excitations remains finite throughout the transitions. The first spin liquid phase corresponds to the well-known zero-flux state in the Heisenberg limit, which is stable with respect to small Kitaev couplings and develops ${120}^{\ensuremath{\circ}}$ order in the semiclassical limit at large S. In the opposite Kitaev limit, we find a different spin liquid ground state, which is a quantum disordered version of a magnetically ordered state with antiferromagnetic chains, in accordance with results in the classical limit. Finally, at intermediate couplings, we find a spin liquid state with unusual spin correlations. Upon spinon condensation, this state develops Bragg peaks at incommensurate momenta in close analogy to the magnetically ordered Z2 vortex crystal phase, which has been analyzed in recent theoretical works.