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

Real-space spectral simulation of quantum spin models: Application to the Kitaev-Heisenberg model

Francisco M. O. Brito, Aires Ferreira

Published Oct 4, 2021• 0 citations
OpenAlex

Abstract

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 Chebyshev iterative method that gives access to the thermodynamic properties and critical behavior of frustrated quantum spin models with good accuracy. The computational complexity scales linearly with the Hilbert space dimension and the number of Chebyshev iterations used to approximate the eigenstates. Using this approach, we calculate the spin correlations of the Kitaev-Heisenberg model, a paradigmatic model of honeycomb iridates that exhibits a rich phase diagram including a quantum spin liquid phase. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum (TPQ) states. All methods accurately predict a transition to a stripy (spin-liquid) phase for the critical value of the Kitaev coupling $ J_K \approx -1.3 J_H $ ($J_K \approx -8.0 J_H$) for honeycomb layers with ferromagnetic Heisenberg interactions ($J_{H}>0$). Our findings suggest that a hybrid Chebyshev-TPQ approach could open the door to previously unattainable studies of quantum spin models in two dimensions.

Related papers

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

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

Dorota Gotfryd, Juraj Rusnačko, Krzysztof Wohlfeld +3 more

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.

arXiv:1608.05333•Physical review. B./Physical review. B•Jan 23, 2017•75 citations

Spin- S Kitaev-Heisenberg model on the honeycomb lattice: A high-order treatment via the many-body coupled cluster method

M. Georgiou, Ioannis Rousochatzakis, DAMIAN J. J. FARNELL +2 more

We study the spin- S Kitaev-Heisenberg model on the honeycomb lattice for S = 1 / 2 , 1, and 3 / 2 , by using the coupled cluster method (CCM) of microscopic quantum many-body theory. This system is one of the earliest extensions of the Kitaev model and is believed to contain two extended spin liquid phases for any value of the spin quantum number S . We show that the CCM delivers accurate estimates for the phase boundaries of these spin liquid phases, as well as other transition points in the phase diagram. Moreover, we find evidence of two unexpected narrow phases for S = 1 / 2 , one sandwiched between the zigzag and ferromagnetic phases and the other between the Néel and the stripy phases. The results establish the CCM as a versatile numerical technique that can capture the strong quantum-mechanical fluctuations that are inherently present in generalized Kitaev models with competing bond-dependent anisotropies. Published by the American Physical Society 2024

arXiv:2405.14378•Physical Review Research•Aug 13, 2024•14 citations

Kitaev-like honeycomb magnets: Global phase behavior and emergent effective models

Juraj Rusnačko, Dorota Gotfryd, Jiří Chaloupka

Compounds of transition metal ions with strong spin-orbit coupling recently attracted attention due to the possibility to host frustrated bond-dependent anisotropic magnetic interactions. In general, such interactions lead to complex phase diagrams that may include exotic phases, e.g., the Kitaev spin liquid. Here we report on our comprehensive analysis of the global phase diagram of the extended Kitaev-Heisenberg model relevant to honeycomb lattice compounds ${\mathrm{Na}}_{2}{\mathrm{IrO}}_{3}$ and $\ensuremath{\alpha}\ensuremath{-}{\mathrm{RuCl}}_{3}$. We have utilized recently developed method based on spin coherent states that enabled us to resolve arbitrary spin patterns in the cluster ground states obtained by exact diagonalization. Global trends in the phase diagram are understood in combination with the analytical mappings of the Hamiltonian that uncover peculiar links to known models---Heisenberg, Ising, Kitaev, or compass models on the honeycomb lattice---or reveal entire manifolds of exact fluctuation-free ground states. Finally, our study can serve as a methodological example that can be applied to other spin models with complex bond-dependent non-Heisenberg interactions.

arXiv:1902.07451•Physical review. B./Physical review. B•Feb 19, 2019•55 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

A Study of Highly Frustrated Spin Systems with mixed PEPS in Infinite Honeycomb Lattice

Huan He, Zhen Wang, Chuan‐Feng Li +2 more

Highly frustrated spin systems represent a central and challenging problem in condensed mater physics. To this problem, we introduce an algorithm based on mixed projected entangled pair states (m-PEPS), which is a novel type of tensor network. We use the famous Kitaev model on an infinite honeycomb lattice, which can be solved exactly, as a benchmark. With very limited parameters and finite scaling, our calculation results are in good agreement with the exact results, indicating the efficiency of our algorithm. After presenting the benchmark, we investigate the Kitaev-Heisenberg model, which was proposed to describe the effective magnetic momentum interaction in iridate Na$_2$IrO$_3$ which may be used to realize the spin liquid phase. However, our calculations suggest that the gapless spin liquid phase is not robust at the thermodynamic limit, and thus this phenomenon is very difficult to observe.

arXiv:1303.2431•Mar 11, 2013•0 citations
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 Chebyshev iterative method that gives access to the thermodynamic properties and critical behavior of frustrated quantum spin models with good accuracy. The computational complexity scales linearly with the Hilbert space dimension and the number of Chebyshev iterations used to approximate the eigenstates. Using this approach, we calculate the spin correlations of the Kitaev-Heisenberg model, a paradigmatic model of honeycomb iridates that exhibits a rich phase diagram including a quantum spin liquid phase. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum (TPQ) states. All methods accurately predict a transition to a stripy (spin-liquid) phase for the critical value of the Kitaev coupling JK≈−1.3JHJ_K \approx -1.3 J_HJK​≈−1.3JH​ (JK≈−8.0JHJ_K \approx -8.0 J_HJK​≈−8.0JH​) for honeycomb layers with ferromagnetic Heisenberg interactions (JH>0J_{H}>0JH​>0). Our findings suggest that a hybrid Chebyshev-TPQ approach could open the door to previously unattainable studies of quantum spin models in two dimensions.
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.
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 spin- S Kitaev-Heisenberg model on the honeycomb lattice for S = 1 / 2 , 1, and 3 / 2 , by using the coupled cluster method (CCM) of microscopic quantum many-body theory. This system is one of the earliest extensions of the Kitaev model and is believed to contain two extended spin liquid phases for any value of the spin quantum number S . We show that the CCM delivers accurate estimates for the phase boundaries of these spin liquid phases, as well as other transition points in the phase diagram. Moreover, we find evidence of two unexpected narrow phases for S = 1 / 2 , one sandwiched between the zigzag and ferromagnetic phases and the other between the Néel and the stripy phases. The results establish the CCM as a versatile numerical technique that can capture the strong quantum-mechanical fluctuations that are inherently present in generalized Kitaev models with competing bond-dependent anisotropies. Published by the American Physical Society 2024
Compounds of transition metal ions with strong spin-orbit coupling recently attracted attention due to the possibility to host frustrated bond-dependent anisotropic magnetic interactions. In general, such interactions lead to complex phase diagrams that may include exotic phases, e.g., the Kitaev spin liquid. Here we report on our comprehensive analysis of the global phase diagram of the extended Kitaev-Heisenberg model relevant to honeycomb lattice compounds Na2IrO3{\mathrm{Na}}_{2}{\mathrm{IrO}}_{3}Na2​IrO3​ and $\ensuremath{\alpha}\ensuremath{-}{\mathrm{RuCl}}_{3}$. We have utilized recently developed method based on spin coherent states that enabled us to resolve arbitrary spin patterns in the cluster ground states obtained by exact diagonalization. Global trends in the phase diagram are understood in combination with the analytical mappings of the Hamiltonian that uncover peculiar links to known models---Heisenberg, Ising, Kitaev, or compass models on the honeycomb lattice---or reveal entire manifolds of exact fluctuation-free ground states. Finally, our study can serve as a methodological example that can be applied to other spin models with complex bond-dependent non-Heisenberg interactions.
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.
Highly frustrated spin systems represent a central and challenging problem in condensed mater physics. To this problem, we introduce an algorithm based on mixed projected entangled pair states (m-PEPS), which is a novel type of tensor network. We use the famous Kitaev model on an infinite honeycomb lattice, which can be solved exactly, as a benchmark. With very limited parameters and finite scaling, our calculation results are in good agreement with the exact results, indicating the efficiency of our algorithm. After presenting the benchmark, we investigate the Kitaev-Heisenberg model, which was proposed to describe the effective magnetic momentum interaction in iridate Na2_22​IrO3_33​ which may be used to realize the spin liquid phase. However, our calculations suggest that the gapless spin liquid phase is not robust at the thermodynamic limit, and thus this phenomenon is very difficult to observe.