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:2011.11577 · Physical Review Research

Thermodynamics of three-dimensional Kitaev quantum spin liquids via tensor networks

Saeed S. Jahromi, Hadi Yarloo, Román Orús

Published Sep 2, 2021• 13 citations
arXivHTMLPDF

Abstract

We study the 3D Kitaev and Kitaev-Heisenberg models, respectively, on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based on graph projected entangled-pair states (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high temperature leads to fractionalization of spins to Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations, and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the nontrivial Majorana band topology and the presence of Weyl nodes, which manifests itself via nonvanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.

DOI
10.1103/physrevresearch.3.033205

Related papers

Phase diagram and thermal Hall conductivity of the spin-liquid Kekulé-Kitaev model

Fatemeh Mirmojarabian, Mehdi Kargarian, A. Langari

In this work we study the phase diagram of the Kekul\'e-Kitaev model. The model is defined on a honeycomb lattice with bond-dependent anisotropic exchange interactions, making it exactly solvable in terms of Majorana representation of spins, in close analogy to the Kitaev model. However, the energy spectrum of Majorana fermions has a multiband structure characterized by Chern numbers 0, $\ifmmode\pm\else\textpm\fi{}1$, and $\ifmmode\pm\else\textpm\fi{}2$. We obtained the phase diagram of the model in the plane of exchange couplings and in the presence of a magnetic field and found chiral topological and trivial spin-liquid ground states. In the absence of magnetic field most of the phase diagram is a trivial gapped phase continuously connected to an Abelian phase, while in the presence of the magnetic field a topological phase arises. Furthermore, motivated by recent thermal measurements on the spin-liquid candidate $\ensuremath{\alpha}\ensuremath{-}{\mathrm{RuCl}}_{3}$, we calculated the thermal Hall conductivity at different regimes of parameters and temperatures and found the latter is quantized over a wide range of temperatures.

arXiv:1909.13336•Physical review. B./Physical review. B•Mar 11, 2020•4 citations

Thermodynamic classification of three-dimensional Kitaev spin liquids

Tim Eschmann, Petr A. Mishchenko, Kevin M. O'Brien +5 more

In the field of frustrated magnetism, Kitaev models provide a unique framework to study the phenomena of spin fractionalization and emergent lattice gauge theories in two and three spatial dimensions. Their ground states are quantum spin liquids, which can typically be described in terms of a Majorana band structure and an ordering of the underlying ${\mathbb{Z}}_{2}$ gauge structure. Here we provide a comprehensive classification of the ``gauge physics'' of a family of elementary three-dimensional Kitaev models, discussing how their thermodynamics and ground state order depends on the underlying lattice geometry. Using large-scale, sign-free quantum Monte Carlo simulations we show that the ground-state gauge order can generally be understood in terms of the length of elementary plaquettes---a result which extends the applicability of Lieb's theorem to lattice geometries beyond its original scope. At finite temperatures, the proliferation of (gapped) vison excitations destroys the gauge order at a critical temperature scale, which we show to correlate with the size of vison gap for the family of three-dimensional Kitaev models. We also discuss two notable exceptions where the lattice structure gives rise to ``gauge frustration'' or intertwines the gauge ordering with time-reversal symmetry breaking. In a more general context, the thermodynamic gauge transitions in such 3D Kitaev models are one of the most natural settings for phase transitions beyond the standard Landau-Ginzburg-Wilson paradigm.

arXiv:2006.07386•Physical review. B./Physical review. B•Aug 18, 2020•36 citations

Three-dimensional quantum spin liquids in models of harmonic-honeycomb iridates and phase diagram in an infinite- D approximation

Itamar Kimchi, James G. Analytis, Ashvin Vishwanath

Motivated by the recent synthesis of two insulating ${\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}$ polymorphs, where ${\mathrm{Ir}}^{4+}$ ${S}_{\mathrm{eff}}=1/2$ moments form 3D (``harmonic'') honeycomb structures with threefold coordination, we study magnetic Hamiltonians on the resulting $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}$ hyperhoneycomb lattice and $\ensuremath{\gamma}\text{\ensuremath{-}}{\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}$ stripyhoneycomb lattice. Experimentally measured magnetic susceptibilities suggest that Kitaev interactions, predicted for the ideal ${90}^{\ensuremath{\circ}}$ Ir-O-Ir bonds, are sizable in these materials. We first consider pure Kitaev interactions, which lead to an exactly soluble 3D quantum spin liquid (QSL) with emergent Majorana fermions and ${\mathrm{Z}}_{2}$ flux loops. Unlike 2D QSLs, the 3D QSL is stable to finite temperature, with ${T}_{c}\ensuremath{\approx}|K|/100$. On including Heisenberg couplings, exact solubility is lost. However, by noting that the shortest closed loop $\ensuremath{\ell}$ is relatively large in these structures, we construct an $\ensuremath{\ell}\ensuremath{\rightarrow}\ensuremath{\infty}$ approximation by defining the model on the Bethe lattice. The phase diagram of the Kitaev-Heisenberg model on this lattice is obtained directly in the thermodynamic limit, using tensor network states and the infinite-system time-evolving-block-decimation (iTEBD) algorithm. Both magnetically ordered and gapped QSL phases are found, the latter being identified by an entanglement fingerprint.

arXiv:1309.1171•Physical Review B•Nov 19, 2014•137 citations

Majorana-fermion mean field theories of Kitaev quantum spin liquids

Shahnam Ghanbari Saheli, Jennifer Shu-Jen Lin, Huanzhi Hu +1 more

We determine the phase diagrams of anisotropic Kitaev-Heisenberg models on the honeycomb lattice using parton mean-field theories based on different Majorana fermion representations of the S=1/2 spin operators. First, we use a two-dimensional Jordan-Wigner transformation (JWT) involving a semi-infinite snake string operator. To ensure that the fermionized Hamiltonian remains local, we consider the limit of extreme Ising exchange anisotropy in the Heisenberg sector. Second, we use the conventional Kitaev representation in terms of four Majorana fermions subject to local constraints, which we enforce through Lagrange multipliers. For both representations, we self-consistently decouple the interaction terms in the bond and magnetization channels and determine the phase diagrams as a function of the anisotropy of the Kitaev couplings and the relative strength of the Ising exchange. While both mean-field theories produce identical phase boundaries for the topological phase transition between the gapless and gapped Kitaev quantum spin liquids, the JWT fails to correctly describe the magnetic instability and finite-temperature behavior. Our results show that the magnetic phase transition is first order at low temperatures but becomes continuous above a certain temperature. At this energy scale, we also observe a finite-temperature crossover on the quantum-spin-liquid side, from a fractionalized paramagnet at low temperatures, in which gapped flux excitations are frozen out, to a conventional paramagnet at high temperatures. Published by the American Physical Society 2024

arXiv:2310.10230•Physical review. B./Physical review. B•Jan 9, 2024•7 citations

Low-energy Majorana states in spin liquid transitions in a three-dimensional Kitaev model

Joji Nasu, Masafumi Udagawa, Yukitoshi Motome

A three-dimensional Kitaev model on a hyperhoneycomb lattice is investigated numerically at finite temperature. The Kitaev model is one of the solvable quantum spin models, where the ground state is given by gapped and gapless spin liquids, depending on the anisotropy of the interactions. This model can be rewritten as a free Majorana fermion system coupled with Z 2 variables. The density of states of Majorana fermions shows an excitation gap in the gapped region, while it is semimetallic in the gapless region reflecting the Dirac node. Performing the Monte Carlo simulation, we calculate the temperature dependence of the Majorana spectra. We find that the semimetallic dip is filled as temperature increases in the gapless region, but surprisingly, the spectrum develops an excitation gap in the region near the gapless-gapped boundary. Such changes of the low-energy spectrum appear sharply at the transition temperature from the spin liquid to the paramagnetic state. The results indicate that thermal fluctuations of the Z 2 fields significantly influence the low-energy state of Majorana fermions, especially in the spin liquid formation.

arXiv:1409.4865•Journal of Physics Conference Series•Mar 18, 2015•8 citations

Finite-temperature phase transition to a Kitaev spin liquid phase on a hyperoctagon lattice: A large-scale quantum Monte Carlo study

Petr A. Mishchenko, Yasuyuki Kato, Yukitoshi Motome

The quantum spin liquid is an enigmatic quantum state in insulating magnets, in which conventional long-range order is suppressed by strong quantum fluctuations. Recently, an unconventional phase transition was reported between the low-temperature quantum spin liquid and the high-temperature paramagnet in the Kitaev model on a three-dimensional hyperhoneycomb lattice. Here, we show that a similar ``liquid-gas'' transition takes place in another three-dimensional lattice, the hyperoctagon lattice. We investigate the critical phenomena by adopting the Green-function based Monte Carlo technique with the kernel polynomial method, which enables systematic analysis of up to 2048 sites. The critical temperature is lower than that in the hyperhoneycomb case, reflecting the smaller flux gap. We also discuss the transition on the basis of an effective model in the anisotropic limit.

arXiv:1706.05057•Physical review. B./Physical review. B•Sep 15, 2017•22 citations
We study the 3D Kitaev and Kitaev-Heisenberg models, respectively, on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based on graph projected entangled-pair states (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high temperature leads to fractionalization of spins to Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations, and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the nontrivial Majorana band topology and the presence of Weyl nodes, which manifests itself via nonvanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.
In this work we study the phase diagram of the Kekul\'e-Kitaev model. The model is defined on a honeycomb lattice with bond-dependent anisotropic exchange interactions, making it exactly solvable in terms of Majorana representation of spins, in close analogy to the Kitaev model. However, the energy spectrum of Majorana fermions has a multiband structure characterized by Chern numbers 0, $\ifmmode\pm\else\textpm\fi{}1$, and $\ifmmode\pm\else\textpm\fi{}2$. We obtained the phase diagram of the model in the plane of exchange couplings and in the presence of a magnetic field and found chiral topological and trivial spin-liquid ground states. In the absence of magnetic field most of the phase diagram is a trivial gapped phase continuously connected to an Abelian phase, while in the presence of the magnetic field a topological phase arises. Furthermore, motivated by recent thermal measurements on the spin-liquid candidate $\ensuremath{\alpha}\ensuremath{-}{\mathrm{RuCl}}_{3}$, we calculated the thermal Hall conductivity at different regimes of parameters and temperatures and found the latter is quantized over a wide range of temperatures.
In the field of frustrated magnetism, Kitaev models provide a unique framework to study the phenomena of spin fractionalization and emergent lattice gauge theories in two and three spatial dimensions. Their ground states are quantum spin liquids, which can typically be described in terms of a Majorana band structure and an ordering of the underlying Z2{\mathbb{Z}}_{2}Z2​ gauge structure. Here we provide a comprehensive classification of the ``gauge physics'' of a family of elementary three-dimensional Kitaev models, discussing how their thermodynamics and ground state order depends on the underlying lattice geometry. Using large-scale, sign-free quantum Monte Carlo simulations we show that the ground-state gauge order can generally be understood in terms of the length of elementary plaquettes---a result which extends the applicability of Lieb's theorem to lattice geometries beyond its original scope. At finite temperatures, the proliferation of (gapped) vison excitations destroys the gauge order at a critical temperature scale, which we show to correlate with the size of vison gap for the family of three-dimensional Kitaev models. We also discuss two notable exceptions where the lattice structure gives rise to ``gauge frustration'' or intertwines the gauge ordering with time-reversal symmetry breaking. In a more general context, the thermodynamic gauge transitions in such 3D Kitaev models are one of the most natural settings for phase transitions beyond the standard Landau-Ginzburg-Wilson paradigm.
Motivated by the recent synthesis of two insulating Li2IrO3{\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}Li2​IrO3​ polymorphs, where Ir4+{\mathrm{Ir}}^{4+}Ir4+ Seff=1/2{S}_{\mathrm{eff}}=1/2Seff​=1/2 moments form 3D (``harmonic'') honeycomb structures with threefold coordination, we study magnetic Hamiltonians on the resulting $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}$ hyperhoneycomb lattice and $\ensuremath{\gamma}\text{\ensuremath{-}}{\mathrm{Li}}_{2}{\mathrm{IrO}}_{3}$ stripyhoneycomb lattice. Experimentally measured magnetic susceptibilities suggest that Kitaev interactions, predicted for the ideal ${90}^{\ensuremath{\circ}}$ Ir-O-Ir bonds, are sizable in these materials. We first consider pure Kitaev interactions, which lead to an exactly soluble 3D quantum spin liquid (QSL) with emergent Majorana fermions and Z2{\mathrm{Z}}_{2}Z2​ flux loops. Unlike 2D QSLs, the 3D QSL is stable to finite temperature, with ${T}_{c}\ensuremath{\approx}|K|/100$. On including Heisenberg couplings, exact solubility is lost. However, by noting that the shortest closed loop $\ensuremath{\ell}$ is relatively large in these structures, we construct an $\ensuremath{\ell}\ensuremath{\rightarrow}\ensuremath{\infty}$ approximation by defining the model on the Bethe lattice. The phase diagram of the Kitaev-Heisenberg model on this lattice is obtained directly in the thermodynamic limit, using tensor network states and the infinite-system time-evolving-block-decimation (iTEBD) algorithm. Both magnetically ordered and gapped QSL phases are found, the latter being identified by an entanglement fingerprint.
We determine the phase diagrams of anisotropic Kitaev-Heisenberg models on the honeycomb lattice using parton mean-field theories based on different Majorana fermion representations of the S=1/2 spin operators. First, we use a two-dimensional Jordan-Wigner transformation (JWT) involving a semi-infinite snake string operator. To ensure that the fermionized Hamiltonian remains local, we consider the limit of extreme Ising exchange anisotropy in the Heisenberg sector. Second, we use the conventional Kitaev representation in terms of four Majorana fermions subject to local constraints, which we enforce through Lagrange multipliers. For both representations, we self-consistently decouple the interaction terms in the bond and magnetization channels and determine the phase diagrams as a function of the anisotropy of the Kitaev couplings and the relative strength of the Ising exchange. While both mean-field theories produce identical phase boundaries for the topological phase transition between the gapless and gapped Kitaev quantum spin liquids, the JWT fails to correctly describe the magnetic instability and finite-temperature behavior. Our results show that the magnetic phase transition is first order at low temperatures but becomes continuous above a certain temperature. At this energy scale, we also observe a finite-temperature crossover on the quantum-spin-liquid side, from a fractionalized paramagnet at low temperatures, in which gapped flux excitations are frozen out, to a conventional paramagnet at high temperatures. Published by the American Physical Society 2024
A three-dimensional Kitaev model on a hyperhoneycomb lattice is investigated numerically at finite temperature. The Kitaev model is one of the solvable quantum spin models, where the ground state is given by gapped and gapless spin liquids, depending on the anisotropy of the interactions. This model can be rewritten as a free Majorana fermion system coupled with Z 2 variables. The density of states of Majorana fermions shows an excitation gap in the gapped region, while it is semimetallic in the gapless region reflecting the Dirac node. Performing the Monte Carlo simulation, we calculate the temperature dependence of the Majorana spectra. We find that the semimetallic dip is filled as temperature increases in the gapless region, but surprisingly, the spectrum develops an excitation gap in the region near the gapless-gapped boundary. Such changes of the low-energy spectrum appear sharply at the transition temperature from the spin liquid to the paramagnetic state. The results indicate that thermal fluctuations of the Z 2 fields significantly influence the low-energy state of Majorana fermions, especially in the spin liquid formation.
The quantum spin liquid is an enigmatic quantum state in insulating magnets, in which conventional long-range order is suppressed by strong quantum fluctuations. Recently, an unconventional phase transition was reported between the low-temperature quantum spin liquid and the high-temperature paramagnet in the Kitaev model on a three-dimensional hyperhoneycomb lattice. Here, we show that a similar ``liquid-gas'' transition takes place in another three-dimensional lattice, the hyperoctagon lattice. We investigate the critical phenomena by adopting the Green-function based Monte Carlo technique with the kernel polynomial method, which enables systematic analysis of up to 2048 sites. The critical temperature is lower than that in the hyperhoneycomb case, reflecting the smaller flux gap. We also discuss the transition on the basis of an effective model in the anisotropic limit.