Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1608.05333 · Physical review. B./Physical review. B

Phase diagram and spin correlations of the Kitaev-Heisenberg model: Importance of quantum effects

Dorota Gotfryd, Juraj Rusnačko, Krzysztof Wohlfeld, George Jackeli, Jiří Chaloupka, Andrzej M. Oleś

Published Jan 23, 2017• 75 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.95.024426

Related papers

Disorder effects in the Kitaev-Heisenberg model

Ayushi Singhania, Jeroen van den Brink, Satoshi Nishimoto

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.

arXiv:2210.17198•Physical Review Research•Apr 5, 2023•10 citations

Quantum phase diagram of the bilayer Kitaev-Heisenberg model

Elahe Samimi, Antonia Duft, Patrick Adelhardt +2 more

We study the ground-state phase diagram of the spin-$1/2$ Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.

arXiv:2412.17495•Physical review. B./Physical review. B•Apr 16, 2026•0 citations

Real-space spectral simulation of quantum spin models: Application to generalized Kitaev models

Francisco M. O. Brito, Aires Ferreira

The proliferation of quantum fluctuations and long-range entanglement presents an outstanding challenge for the numerical simulation of interacting spin systems with exotic ground states. Here, we present a toolset of Chebyshev polynomial-based iterative methods that provides a unified framework to study the thermodynamical properties, critical behavior and dynamics of frustrated quantum spin models with controlled accuracy. Similar to previous applications of the Chebyshev spectral methods to condensed matter systems, the algorithmic complexity scales linearly with the Hilbert space dimension and the Chebyshev truncation order. Using this approach, we study two paradigmatic quantum spin models on the honeycomb lattice: the Kitaev-Heisenberg (K-H) and the Kitaev-Ising (K-I) models. We start by applying the Chebyshev toolset to compute nearest-neighbor spin correlations, specific heat and entropy of the K-H model on a 24-spin cluster. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum states. The transitions between a variety of magnetic phases, namely ferromagnetic, Néel, zigzag and stripy antiferromagnetic and quantum spin liquid phases are obtained accurately and efficiently. We also determine the temperature dependence of the spin correlations, over more than three decades in temperature, by means of a finite temperature Chebyshev polynomial method introduced here. Finally, we report novel dynamical signatures of the quantum phase transitions in the K-I model. Our findings suggest that the efficiency, versatility and low-temperature stability of the Chebyshev framework developed here could pave the way for previously unattainable studies of quantum spin models in two dimensions.

arXiv:2110.01494•SciPost Physics Core•Feb 21, 2024•1 citation

Ordered states in the Kitaev-Heisenberg model: From 1D chains to 2D honeycomb

Cliò Efthimia Agrapidis, Jeroen van den Brink, Satoshi Nishimoto

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 $${\boldsymbol{\phi }}{\boldsymbol{=}}\frac{{\boldsymbol{\pi }}}{{\bf{2}}}$$ ϕ=π2 and $${\boldsymbol{\varphi }}{\boldsymbol{=}}\frac{{\bf{3}}{\boldsymbol{\pi }}}{{\bf{2}}}$$ φ=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.

arXiv:1712.07418•Scientific Reports•Jan 23, 2018•38 citations

Ordered states in the Kitaev-Heisenberg model: From 1D chains to 2D honeycomb

Cliò Efthimia Agrapidis, Jeroen van den Brink, Satoshi Nishimoto

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.

•Jan 1, 2018•0 citations

Quantum phases in the frustrated Heisenberg model on the bilayer honeycomb lattice

Hao Zhang, Marcelo Arlego, Carlos Alberto Lamas

We use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling ${J}_{2}$ and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin $S$, the frustration, and the interlayer couplings.

arXiv:1312.6324•Physical Review B•Jan 6, 2014•31 citations
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 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 H3LiIr2O6{\mathrm{H}}_{3}{\mathrm{LiIr}}_{2}{\mathrm{O}}_{6}H3​LiIr2​O6​ 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 ground-state phase diagram of the spin-1/21/21/2 Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.
The proliferation of quantum fluctuations and long-range entanglement presents an outstanding challenge for the numerical simulation of interacting spin systems with exotic ground states. Here, we present a toolset of Chebyshev polynomial-based iterative methods that provides a unified framework to study the thermodynamical properties, critical behavior and dynamics of frustrated quantum spin models with controlled accuracy. Similar to previous applications of the Chebyshev spectral methods to condensed matter systems, the algorithmic complexity scales linearly with the Hilbert space dimension and the Chebyshev truncation order. Using this approach, we study two paradigmatic quantum spin models on the honeycomb lattice: the Kitaev-Heisenberg (K-H) and the Kitaev-Ising (K-I) models. We start by applying the Chebyshev toolset to compute nearest-neighbor spin correlations, specific heat and entropy of the K-H model on a 24-spin cluster. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum states. The transitions between a variety of magnetic phases, namely ferromagnetic, Néel, zigzag and stripy antiferromagnetic and quantum spin liquid phases are obtained accurately and efficiently. We also determine the temperature dependence of the spin correlations, over more than three decades in temperature, by means of a finite temperature Chebyshev polynomial method introduced here. Finally, we report novel dynamical signatures of the quantum phase transitions in the K-I model. Our findings suggest that the efficiency, versatility and low-temperature stability of the Chebyshev framework developed here could pave the way for previously unattainable studies of quantum spin models in two dimensions.
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{\boldsymbol{\phi }}{\boldsymbol{=}}\frac{{\boldsymbol{\pi }}}{{\bf{2}}}ϕ=2π​ ϕ=π2 and φ=3π2{\boldsymbol{\varphi }}{\boldsymbol{=}}\frac{{\bf{3}}{\boldsymbol{\pi }}}{{\bf{2}}}φ=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 use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling J2{J}_{2}J2​ and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin SSS, the frustration, and the interlayer couplings.