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:2005.01255 · Physical review. E

Thermodynamic behavior of modified integer-spin Kitaev models on the honeycomb lattice

Owen Bradley, J. Oitmaa, Diptiman Sen, Rajiv R. P. Singh

Published Feb 5, 2021• 1 citation
arXivHTMLPDF

Abstract

We study the thermodynamic behavior of modified spin-S Kitaev models introduced by Baskaran, Sen, and Shankar [Phys. Rev. B 78, 115116 (2008)PRBMDO1098-012110.1103/PhysRevB.78.115116]. These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high-temperature expansions, and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin S=1 and S=2. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as T→0, become critical as T→0, or diverge exponentially as T→0. The Z_{2} flux variable associated with each hexagonal plaquette saturates at the value +1 as T→0 in all models except the S=1 antiferromagnet where the mean flux remains zero as T→0. We provide qualitative explanations for these results.

DOI
10.1103/physreve.103.022109

Related papers

Revised thermodynamics of the Kitaev honeycomb model

Kexin Feng, Natalia Borisovna Perkins, F. J. Burnell

Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes, and verify using large-scale Monte Carlo simulations that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. % Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically, and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.

•Aug 4, 2020•0 citations

Further insights into the thermodynamics of the Kitaev honeycomb model

Kexin Feng, Natalia Borisovna Perkins, F. J. Burnell

Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes. Using large-scale Monte Carlo simulations, we verify that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.

arXiv:2008.01822•Physical review. B./Physical review. B•Dec 1, 2020•19 citations

Majorana dynamical mean-field study of spin dynamics at finite temperatures in the honeycomb Kitaev model

Junki Yoshitake, Joji Nasu, Yasuyuki Kato +1 more

A prominent feature of quantum spin liquids is fractionalization of the spin degree of freedom. Fractionalized excitations have their own dynamics in different energy scales, and hence, affect finite-temperature ($T$) properties in a peculiar manner even in the paramagnetic state harboring the quantum spin liquid state. We here present a comprehensive theoretical study of the spin dynamics in a wide $T$ range for the Kitaev model on a honeycomb lattice, whose ground state is a quantum spin liquid. In this model, the fractionalization occurs to break up quantum spins into itinerant matter fermions and localized ${Z}_{2}$ fluxes, which results in two crossovers at very different $T$ scales. Extending the previous study for the isotropic coupling case [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016)], we calculate the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$, the NMR relaxation rate $1/{T}_{1}$, and the magnetic susceptibility $\ensuremath{\chi}$ while changing the anisotropy in the exchange coupling constants, by using the dynamical mean-field theory based on a Majorana fermion representation. We describe the details of the methodology including the continuous-time quantum Monte Carlo method for computing dynamical spin correlations and the maximum entropy method for analytic continuation. We confirm that the combined method provides accurate results in a wide $T$ range including the region where the spins are fractionalized. We find that also in the anisotropic cases the system exhibits peculiar behaviors below the high-$T$ crossover whose temperature is comparable to the average of the exchange constants: $S(\mathbf{q},\ensuremath{\omega})$ shows an inelastic response at the energy scale of the averaged exchange constant, $1/{T}_{1}$ continues to grow even though the equal-time spin correlations are saturated and almost $T$ independent, and $\ensuremath{\chi}$ deviates from the Curie-Weiss behavior. In particular, when the exchange interaction in one direction is stronger than the other two, the dynamical quantities exhibit qualitatively different $T$ dependences from the isotropic case at low $T$, reflecting the opposite parity between the flux-free ground state and the flux-excited state, and a larger energy cost for flipping a spin in the strong interaction direction. On the other hand, when the exchange anisotropy is in the opposite way, the results are qualitatively similar to those in the isotropic case. All these behaviors manifest the spin fractionalization in the paramagnetic region. Among them, the dichotomy between the static and dynamical spin correlations is unusual behavior hardly seen in conventional magnets. We discuss the relation between the dichotomy and the spatial configuration of the ${Z}_{2}$ gauge fields. Our results could stimulate further experimental and theoretical analyses of candidate materials for the Kitaev quantum spin liquids.

arXiv:1704.02707•Physical review. B./Physical review. B•Jul 25, 2017•43 citations

Spin-S Decorated Honeycomb Kitaev Model with S ≥ 1

Tomonari Mizoguchi

In the prior work [J. Phys. Soc. Jpn. 91 , 124709 (2022)], we have introduced the spin-half decorated honeycomb Kitaev model, and have elucidated its exactly-solvable nature at any temperature. In this work, we extend the model to that with an arbitrary spin quantum number S . We show that the ground state is classical spin liquid whose degeneracy is independent of S . As for the fi nite-temperature properties, the partition function cannot be obtained analytically, unlike the spin-half case. We therefore perform the Monte Carlo simulation to reveal the fi nite-temperature properties, taking S = 1 and 3/2 cases as examples. We fi nd that the temperature dependence of the energy and the specific heat exhibits typical behaviors of the classical spin liquid. In particular, the broad peak of the specific heat is accounted for by the configuration of the classical variables on the edge of the hexagons.

•Nov 5, 2024•0 citations

Z2 Spin Liquids in the Higher Spin- S Kitaev Honeycomb Model: An Exact Deconfined Z2 Gauge Structure in a Nonintegrable Model

Han Ma

The higher spin Kitaev model prominently features the extensive locally conserved quantities the same as the spin-1/2 Kitaev honeycomb model, although it is not exactly solvable. It remains an open question regarding the physical meaning of these conserved quantities in the higher spin model. In this Letter, by introducing a Majorana parton construction for a general spin-S we uncover that these conserved quantities are exactly the Z_{2} gauge fluxes in the general spin-S model, including the case of spin-1/2. Particularly, we find an even-odd effect that the Z_{2} gauge charges are fermions in the half integer spin model, but are bosons in the integer spin model. We further prove that the fermionic Z_{2} gauge charges are always deconfined; hence, the half integer spin Kitaev model would have nontrivial spin liquid ground states regardless of interaction strengths in the Hamiltonian. The bosonic Z_{2} gauge charges of the integer spin model, on the other hand, could condense, leading to a trivial product state, and this is indeed the case at the anisotropic limit of the model.

arXiv:2212.00053•Physical Review Letters•Apr 11, 2023•27 citations

Temperature evolution of spin dynamics in two- and three-dimensional Kitaev models: Influence of fluctuating Z2 flux

Junki Yoshitake, Joji Nasu, Yukitoshi Motome

The long-sought quantum spin liquid is a quantum-entangled magnetic state leading to the fractionalization of spin degrees of freedom. Quasiparticles emergent from the fractionalization affect not only the ground-state properties but also thermodynamic behavior in a peculiar manner. Here we investigate how the spin dynamics evolves from the high-temperature paramagnet to the quantum spin-liquid ground state, for the Kitaev spin model describing the fractionalization into itinerant matter fermions and localized ${Z}_{2}$ fluxes. Beyond the previous studies [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016); J. Yoshitake, J. Nasu, Y. Kato, and Y. Motome, Phys. Rev. B 96, 024438 (2017)], in which the mean-field nature of the cluster dynamical mean-field theory prevented us from studying low-temperature properties, we develop a numerical technique by applying the continuous-time quantum Monte Carlo (CTQMC) method to statistical samples generated by the quantum Monte Carlo (QMC) method in a Majorana fermion representation. This QMC+CTQMC method is fully unbiased and enables us to investigate the low-temperature spin dynamics dominated by thermally excited ${Z}_{2}$ fluxes, including the unconventional phase transition caused by ${Z}_{2}$ flux loops in three dimensions, which was unreachable by the previous methods. We apply this technique to the Kitaev model in both two and three dimensions. Our results clearly distinguish two cases: while the dynamics changes smoothly through the crossover in the two-dimensional honeycomb case, it exhibits singular behaviors at the phase transition in the three-dimensional hyperhoneycomb case. We show that the low-temperature spin dynamics is a sensitive probe for thermally fluctuating ${Z}_{2}$ fluxes that behave very differently between two and three dimensions.

arXiv:1705.07760•Physical review. B./Physical review. B•Aug 28, 2017•35 citations
We study the thermodynamic behavior of modified spin-S Kitaev models introduced by Baskaran, Sen, and Shankar [Phys. Rev. B 78, 115116 (2008)PRBMDO1098-012110.1103/PhysRevB.78.115116]. These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high-temperature expansions, and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin S=1 and S=2. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as T→0, become critical as T→0, or diverge exponentially as T→0. The Z_{2} flux variable associated with each hexagonal plaquette saturates at the value +1 as T→0 in all models except the S=1 antiferromagnet where the mean flux remains zero as T→0. We provide qualitative explanations for these results.
Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes, and verify using large-scale Monte Carlo simulations that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. % Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically, and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.
Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes. Using large-scale Monte Carlo simulations, we verify that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.
A prominent feature of quantum spin liquids is fractionalization of the spin degree of freedom. Fractionalized excitations have their own dynamics in different energy scales, and hence, affect finite-temperature (TTT) properties in a peculiar manner even in the paramagnetic state harboring the quantum spin liquid state. We here present a comprehensive theoretical study of the spin dynamics in a wide TTT range for the Kitaev model on a honeycomb lattice, whose ground state is a quantum spin liquid. In this model, the fractionalization occurs to break up quantum spins into itinerant matter fermions and localized Z2{Z}_{2}Z2​ fluxes, which results in two crossovers at very different TTT scales. Extending the previous study for the isotropic coupling case [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016)], we calculate the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$, the NMR relaxation rate 1/T11/{T}_{1}1/T1​, and the magnetic susceptibility $\ensuremath{\chi}$ while changing the anisotropy in the exchange coupling constants, by using the dynamical mean-field theory based on a Majorana fermion representation. We describe the details of the methodology including the continuous-time quantum Monte Carlo method for computing dynamical spin correlations and the maximum entropy method for analytic continuation. We confirm that the combined method provides accurate results in a wide TTT range including the region where the spins are fractionalized. We find that also in the anisotropic cases the system exhibits peculiar behaviors below the high-TTT crossover whose temperature is comparable to the average of the exchange constants: $S(\mathbf{q},\ensuremath{\omega})$ shows an inelastic response at the energy scale of the averaged exchange constant, 1/T11/{T}_{1}1/T1​ continues to grow even though the equal-time spin correlations are saturated and almost TTT independent, and $\ensuremath{\chi}$ deviates from the Curie-Weiss behavior. In particular, when the exchange interaction in one direction is stronger than the other two, the dynamical quantities exhibit qualitatively different TTT dependences from the isotropic case at low TTT, reflecting the opposite parity between the flux-free ground state and the flux-excited state, and a larger energy cost for flipping a spin in the strong interaction direction. On the other hand, when the exchange anisotropy is in the opposite way, the results are qualitatively similar to those in the isotropic case. All these behaviors manifest the spin fractionalization in the paramagnetic region. Among them, the dichotomy between the static and dynamical spin correlations is unusual behavior hardly seen in conventional magnets. We discuss the relation between the dichotomy and the spatial configuration of the Z2{Z}_{2}Z2​ gauge fields. Our results could stimulate further experimental and theoretical analyses of candidate materials for the Kitaev quantum spin liquids.
In the prior work [J. Phys. Soc. Jpn. 91 , 124709 (2022)], we have introduced the spin-half decorated honeycomb Kitaev model, and have elucidated its exactly-solvable nature at any temperature. In this work, we extend the model to that with an arbitrary spin quantum number S . We show that the ground state is classical spin liquid whose degeneracy is independent of S . As for the fi nite-temperature properties, the partition function cannot be obtained analytically, unlike the spin-half case. We therefore perform the Monte Carlo simulation to reveal the fi nite-temperature properties, taking S = 1 and 3/2 cases as examples. We fi nd that the temperature dependence of the energy and the specific heat exhibits typical behaviors of the classical spin liquid. In particular, the broad peak of the specific heat is accounted for by the configuration of the classical variables on the edge of the hexagons.
The higher spin Kitaev model prominently features the extensive locally conserved quantities the same as the spin-1/2 Kitaev honeycomb model, although it is not exactly solvable. It remains an open question regarding the physical meaning of these conserved quantities in the higher spin model. In this Letter, by introducing a Majorana parton construction for a general spin-S we uncover that these conserved quantities are exactly the Z_{2} gauge fluxes in the general spin-S model, including the case of spin-1/2. Particularly, we find an even-odd effect that the Z_{2} gauge charges are fermions in the half integer spin model, but are bosons in the integer spin model. We further prove that the fermionic Z_{2} gauge charges are always deconfined; hence, the half integer spin Kitaev model would have nontrivial spin liquid ground states regardless of interaction strengths in the Hamiltonian. The bosonic Z_{2} gauge charges of the integer spin model, on the other hand, could condense, leading to a trivial product state, and this is indeed the case at the anisotropic limit of the model.
The long-sought quantum spin liquid is a quantum-entangled magnetic state leading to the fractionalization of spin degrees of freedom. Quasiparticles emergent from the fractionalization affect not only the ground-state properties but also thermodynamic behavior in a peculiar manner. Here we investigate how the spin dynamics evolves from the high-temperature paramagnet to the quantum spin-liquid ground state, for the Kitaev spin model describing the fractionalization into itinerant matter fermions and localized Z2{Z}_{2}Z2​ fluxes. Beyond the previous studies [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016); J. Yoshitake, J. Nasu, Y. Kato, and Y. Motome, Phys. Rev. B 96, 024438 (2017)], in which the mean-field nature of the cluster dynamical mean-field theory prevented us from studying low-temperature properties, we develop a numerical technique by applying the continuous-time quantum Monte Carlo (CTQMC) method to statistical samples generated by the quantum Monte Carlo (QMC) method in a Majorana fermion representation. This QMC+CTQMC method is fully unbiased and enables us to investigate the low-temperature spin dynamics dominated by thermally excited Z2{Z}_{2}Z2​ fluxes, including the unconventional phase transition caused by Z2{Z}_{2}Z2​ flux loops in three dimensions, which was unreachable by the previous methods. We apply this technique to the Kitaev model in both two and three dimensions. Our results clearly distinguish two cases: while the dynamics changes smoothly through the crossover in the two-dimensional honeycomb case, it exhibits singular behaviors at the phase transition in the three-dimensional hyperhoneycomb case. We show that the low-temperature spin dynamics is a sensitive probe for thermally fluctuating Z2{Z}_{2}Z2​ fluxes that behave very differently between two and three dimensions.