pith. sign in

arxiv: 2603.08710 · v2 · pith:EC4GQICKnew · submitted 2026-03-09 · ❄️ cond-mat.str-el

Thermal and quantum fluctuations in extended Kitaev-Yao-Lee spin-orbital model

Jiefu Cen , Hae-Young Kee This is my paper

Pith reviewed 2026-05-21 11:28 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords Kitaev modelYao-Lee modelspin-orbital liquidnematic phasequantum fluctuationsthermal fluctuationsdisordered phaseshoneycomb lattice
0
0 comments X

The pith

The extra spin-orbital degree of freedom produces strong thermal and quantum fluctuations that stabilize disordered phases in extended Kitaev-Yao-Lee models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper studies the extended Kitaev-Yao-Lee model that adds Kugel-Khomskii interactions to the solvable Yao-Lee spin-orbital Hamiltonian on the honeycomb lattice. The authors apply classical Monte Carlo simulations to track thermal fluctuations and a generalized spin wave theory to quantify quantum fluctuations. They show that the orbital degree of freedom amplifies both kinds of fluctuations, which in turn suppress magnetic order and support a wide nematic phase resembling a spin-orbital liquid. A reader would care because the result offers a concrete mechanism for why certain materials with spin-orbit coupling avoid conventional ordering down to low temperatures.

Core claim

The additional spin-orbital degree of freedom in the extended Kitaev-Yao-Lee model gives rise to strong thermal and quantum fluctuations. These fluctuations are responsible for the emergence and stability of disordered phases, including a broad nematic region. The authors reach this conclusion by combining classical Monte Carlo simulations that capture thermal effects with a generalized spin wave theory tailored to the spin-orbital degrees of freedom.

What carries the argument

Generalized spin wave theory adapted to spin-orbital models, used together with classical Monte Carlo simulations to measure the amplitude of fluctuations generated by the extra orbital degree of freedom.

If this is right

  • Conventional magnetic order is suppressed over wider parameter ranges than in pure spin Kitaev models.
  • The nematic phase occupies a broad region of the phase diagram because fluctuations are stronger.
  • Disordered phases become more stable against perturbations when orbital degrees of freedom are active.
  • Similar fluctuation-driven disordering can be expected in other lattices that host spin-orbital interactions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Materials with strong spin-orbit coupling may host more robust spin-orbital liquids than their spin-only counterparts.
  • The same fluctuation mechanism could be tested in three-dimensional Kitaev-like compounds.
  • Tensor-network or quantum Monte Carlo methods could provide independent checks on the spin-wave spectra inside the nematic region.

Load-bearing premise

The generalized spin wave theory remains quantitatively reliable when the nematic phase is present.

What would settle it

An exact-diagonalization calculation on a finite cluster inside the nematic phase that shows fluctuation spectra differing markedly from those obtained with the generalized spin wave theory would falsify the claim.

Figures

Figures reproduced from arXiv: 2603.08710 by Hae-Young Kee, Jiefu Cen.

Figure 1
Figure 1. Figure 1: FIG. 1. Classical phase diagrams of the spin-orbital model [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Example of a Cartesian state [ [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Structure factors of spin [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Quantum phase diagram of the spin-orbital model [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

Building upon the spin-1/2 Kitaev model on a honeycomb lattice, the Yao-Lee spin-orbital model provides exactly solvable quantum spin liquids with potentially better stability against perturbations due to the additional degree of freedom. Recently, the microscopic mechanism underlying the Yao-Lee interaction in honeycomb materials has been uncovered, leading to an extended Kitaev-Yao-Lee spin-orbital model when the celebrated Kugel-Khomskii interaction is included. Numerical studies of this model have identified various disordered phases, including a broad region of the nematic phase that is reminiscent of a spin-orbital liquid. Here, we investigate the origin and stability of this nematic phase via thermal and quantum fluctuations using classical Monte Carlo simulations and a generalized spin wave theory appropriate for the spin-orbital model. We demonstrate that the additional spin-orbital degree of freedom gives rise to strong thermal and quantum fluctuations in spin-orbital models, providing insight into the emergence of disordered phases.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript examines the extended Kitaev-Yao-Lee spin-orbital model on the honeycomb lattice. It employs classical Monte Carlo simulations to study thermal fluctuations and a generalized spin-wave theory to analyze quantum fluctuations, concluding that the additional spin-orbital degree of freedom generates strong fluctuations that stabilize disordered phases, notably a broad nematic phase reminiscent of a spin-orbital liquid.

Significance. If substantiated, the results offer useful insight into how extra orbital degrees of freedom enhance fluctuation effects and promote disordered phases in Kitaev-like spin-orbital models. The dual use of Monte Carlo for thermal effects and generalized spin-wave theory for quantum corrections is a reasonable methodological choice for this class of models.

major comments (1)
  1. Generalized spin-wave theory section: The method is used to compute fluctuation spectra and strengths inside the nematic phase, yet the manuscript provides no explicit benchmarks such as comparisons of magnon dispersions, equal-time correlations, or specific-heat features against exact diagonalization on finite clusters or quantum Monte Carlo for the same parameters. This validation is load-bearing because the nematic phase is characterized as disordered or highly degenerate, where the harmonic approximation around a classical or mean-field reference state may lose quantitative reliability.
minor comments (2)
  1. Model Hamiltonian section: The notation distinguishing spin and orbital operators could be made more explicit when the Kugel-Khomskii term is introduced, to avoid potential confusion with standard Kitaev notation.
  2. Figure captions (e.g., those showing fluctuation spectra): Inclusion of statistical error estimates from the Monte Carlo runs or convergence checks for the spin-wave expansion would improve clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback on our manuscript. The major comment raises a valid point about validation of the generalized spin-wave theory, which we address below. We will revise the manuscript to incorporate additional benchmarks and discussion as outlined.

read point-by-point responses

  1. Referee: Generalized spin-wave theory section: The method is used to compute fluctuation spectra and strengths inside the nematic phase, yet the manuscript provides no explicit benchmarks such as comparisons of magnon dispersions, equal-time correlations, or specific-heat features against exact diagonalization on finite clusters or quantum Monte Carlo for the same parameters. This validation is load-bearing because the nematic phase is characterized as disordered or highly degenerate, where the harmonic approximation around a classical or mean-field reference state may lose quantitative reliability.

    Authors: We appreciate the referee's emphasis on the need for validation of the generalized spin-wave theory (GSWT). The GSWT is employed to quantify the strength of quantum fluctuations around classical reference states and to obtain the associated spectra; in the nematic phase these calculations indicate fluctuation amplitudes large enough to destabilize long-range order, supporting our interpretation of a disordered phase. We agree that direct benchmarks against exact diagonalization (ED) or quantum Monte Carlo would strengthen the presentation. However, the four-dimensional local Hilbert space per site restricts ED to very small clusters (N ≤ 8), where finite-size effects are severe and quantitative comparison to the thermodynamic limit is limited. We have carried out ED calculations on small clusters in the exactly solvable limits of the model (pure Kitaev and Yao-Lee points) and find qualitative agreement in the overall scale of fluctuations. In the revised manuscript we will add an appendix presenting these comparisons, together with a discussion of the applicability and limitations of the harmonic approximation inside the highly degenerate nematic region. We believe this addition will address the concern. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The manuscript applies classical Monte Carlo simulations and a generalized spin-wave theory to study fluctuations in the extended Kitaev-Yao-Lee model. No equations, self-referential definitions, fitted parameters renamed as predictions, or load-bearing self-citations appear in the abstract or described text. The central demonstration of strong fluctuations from the extra spin-orbital degree of freedom is obtained from these external numerical methods rather than any reduction of a claimed result to its own inputs by construction. The derivation therefore remains self-contained against the stated computational benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the central claim rests on the applicability of the two numerical methods and on the prior identification of the nematic phase.

pith-pipeline@v0.9.0 · 5696 in / 1057 out tokens · 30018 ms · 2026-05-21T11:28:29.896166+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages

  1. [1]

    Anderson, Resonating valence bonds: A new kind of insulator?, Mater

    P. Anderson, Resonating valence bonds: A new kind of insulator?, Mater. Res. Bull.8, 153 (1973)

    work page 1973

  2. [2]

    Balents, Spin liquids in frustrated magnets, Nature 464, 199 (2010)

    L. Balents, Spin liquids in frustrated magnets, Nature 464, 199 (2010)

    work page 2010

  3. [3]

    Savary and L

    L. Savary and L. Balents, Quantum spin liquids: a re- view, Rep. Prog. Phys.80, 016502 (2016)

    work page 2016

  4. [4]

    Y. Zhou, K. Kanoda, and T.-K. Ng, Quantum spin liquid states, Rev. Mod. Phys.89, 025003 (2017)

    work page 2017

  5. [5]

    Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev

    X.-G. Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev. Mod. Phys.89, 041004 (2017)

    work page 2017

  6. [6]

    Hermanns, I

    M. Hermanns, I. Kimchi, and J. Knolle, Physics of the Kitaev model: Fractionalization, dynamic correlations, and material connections, Annu. Rev. Condens. Matter Phys.9, 17 (2018)

    work page 2018

  7. [7]

    Knolle and R

    J. Knolle and R. Moessner, A field guide to spin liquids, Annu. Rev. Condens. Matter Phys.10, 451 (2019)

    work page 2019

  8. [8]

    Takagi, T

    H. Takagi, T. Takayama, G. Jackeli, G. Khaliullin, and S. E. Nagler, Concept and realization of kitaev quantum spin liquids, Nature Reviews Physics1, 264 (2019)

    work page 2019

  9. [9]

    Broholm, R

    C. Broholm, R. J. Cava, S. A. Kivelson, D. G. Nocera, M. R. Norman, and T. Senthil, Quantum spin liquids, Science367, eaay0668 (2020)

    work page 2020

  10. [10]

    Kitaev, Anyons in an exactly solved model and be- yond, Annals of Physics321, 2 (2006)

    A. Kitaev, Anyons in an exactly solved model and be- yond, Annals of Physics321, 2 (2006)

    work page 2006

  11. [11]

    Jackeli and G

    G. Jackeli and G. Khaliullin, Mott insulators in the strong spin-orbit coupling limit: From heisenberg to a quantum compass and kitaev models, Phys. Rev. Lett. 102, 017205 (2009)

    work page 2009

  12. [12]

    J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Generic spin model for the honeycomb iridates beyond the kitaev limit, Phys. Rev. Lett.112, 077204 (2014)

    work page 2014

  13. [13]

    J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Spin-orbit physics giving rise to novel phases in correlated systems: Iridates and related materials, Annu. Rev. Condens. Mat- ter Phys.7, 195 (2016)

    work page 2016

  14. [14]

    Motome, R

    Y. Motome, R. Sano, S. Jang, Y. Sugita, and Y. Kato, Materials design of kitaev spin liquids beyond the jack- eli–khaliullin mechanism, J. Phys. Condens. Matter32, 404001 (2020)

    work page 2020

  15. [15]

    S. M. Winter, Y. Li, H. O. Jeschke, and R. Valent´ ı, Chal- lenges in design of Kitaev materials: Magnetic interac- tions from competing energy scales, Phys. Rev. B93, 214431 (2016)

    work page 2016

  16. [16]

    Liu and G

    H. Liu and G. Khaliullin, Pseudospin exchange interac- tions ind 7 cobalt compounds: Possible realization of the kitaev model, Phys. Rev. B97, 014407 (2018)

    work page 2018

  17. [17]

    H. Liu, J. c. v. Chaloupka, and G. Khaliullin, Kitaev spin liquid in 3dtransition metal compounds, Phys. Rev. Lett. 125, 047201 (2020)

    work page 2020

  18. [18]

    Trebst and C

    S. Trebst and C. Hickey, Kitaev materials, Phys. Rep. 950, 1 (2022)

    work page 2022

  19. [19]

    Liu and H.-Y

    X. Liu and H.-Y. Kee, Non-kitaev versus kitaev honey- comb cobaltates, Phys. Rev. B107, 054420 (2023)

    work page 2023

  20. [20]

    Rousochatzakis, N

    I. Rousochatzakis, N. Perkins, Q. Luo, and H.-Y. Kee, Beyond kitaev physics in spin-orbit coupled magnets, Rep. Prog. Phys.87, 026502 (2024)

    work page 2024

  21. [21]

    Matsuda, T

    Y. Matsuda, T. Shibauchi, and H.-Y. Kee, Kitaev quan- tum spin liquids, Rev. Mod. Phys.97, 045003 (2025)

    work page 2025

  22. [22]

    M¨ oller, P

    M. M¨ oller, P. A. Maksimov, S. Jiang, S. R. White, R. Valenti, and A. L. Chernyshev, The saga ofα- rucl3: Parameters, models, and phase diagrams., arXiv:2502.08698 (2025)

    work page arXiv 2025

  23. [23]

    Yao and D.-H

    H. Yao and D.-H. Lee, Fermionic magnons, non-abelian spinons, and the spin quantum hall effect from an ex- actly solvable spin-1/2 kitaev model with su(2) symme- try, Phys. Rev. Lett.107, 087205 (2011)

    work page 2011

  24. [24]

    Chulliparambil, U

    S. Chulliparambil, U. F. P. Seifert, M. Vojta, L. Janssen, and H.-H. Tu, Microscopic models for kitaev’s sixteenfold way of anyon theories, Phys. Rev. B102, 201111 (2020)

    work page 2020

  25. [25]

    U. F. P. Seifert, X.-Y. Dong, S. Chulliparambil, M. Vo- jta, H.-H. Tu, and L. Janssen, Fractionalized fermionic quantum criticality in spin-orbital mott insulators, Phys. Rev. Lett.125, 257202 (2020)

    work page 2020

  26. [26]

    Chulliparambil, L

    S. Chulliparambil, L. Janssen, M. Vojta, H.-H. Tu, and U. F. P. Seifert, Flux crystals, majorana metals, and flat bands in exactly solvable spin-orbital liquids, Phys. Rev. B103, 075144 (2021)

    work page 2021

  27. [27]

    Akram, E

    M. Akram, E. M. Nica, Y.-M. Lu, and O. Erten, Vi- son crystals, chiral, and crystalline phases in the Yao-Lee model, Phys. Rev. B108, 224427 (2023)

    work page 2023

  28. [28]

    E. Nica, M. Akram, A. Vijayvargia, R. Moessner, and O. Erten, Kitaev spin-orbital bilayers and their moir´ e superlattices, npj Quantum Mater. 8, 9 (2023)

    work page 2023

  29. [29]

    Poliakov, W.-H

    V. Poliakov, W.-H. Kao, and N. B. Perkins, Topological transitions in the Yao-Lee spin-orbital model and effects of site disorder, Phys. Rev. B110, 054418 (2024)

    work page 2024

  30. [30]

    Wu, J.-Y

    Z. Wu, J.-Y. Zhang, and H. Yao, Exact deconfined gauge structures in the higher-spin yao-lee model: A quantum spin-orbital liquid with spin fractionalization and non- abelian anyons, Phys. Rev. Lett.133, 236504 (2024)

    work page 2024

  31. [31]

    Churchill, E

    D. Churchill, E. Z. Zhang, and H.-Y. Kee, Microscopic roadmap to a kitaev-yao-lee spin-orbital liquid, npj Quantum Mater.10, 26 (2025)

    work page 2025

  32. [32]

    K. I. Kugel’ and D. I. Khomskii, Crystal structure and magnetic properties of substances with orbital degen- eracy, Zhurnal Eksperimental’noi i Teoreticheskoi Fiziki 64, 1429 (1973)

    work page 1973

  33. [33]

    P. P. Stavropoulos, D. Pereira, and H.-Y. Kee, Micro- scopic mechanism for a higher-spin kitaev model, Phys. Rev. Lett.123, 037203 (2019)

    work page 2019

  34. [34]

    K. I. Kugel and D. I. Khomskii, The jahn-teller effect and magnetism: transition metal compounds, Soviet Physics Uspekhi25, 231 (1982)

    work page 1982

  35. [35]

    van den Brink and D

    J. van den Brink and D. Khomskii, Orbital ordering of complex orbitals in doped mott insulators, Phys. Rev. B 63, 140416 (2001). 8

    work page 2001

  36. [36]

    M. V. Mostovoy and D. I. Khomskii, Orbital ordering in frustrated jahn-teller systems with 90°exchange, Phys. Rev. Lett.89, 227203 (2002)

    work page 2002

  37. [37]

    D. I. Khomskii, K. I. Kugel, A. O. Sboychakov, and S. V. Streltsov, Role of local geometry in the spin and or- bital structure of transition metal compounds, Sov. Phys. JETP122, 484 (2016)

    work page 2016

  38. [38]

    Hukushima and K

    K. Hukushima and K. Nemoto, Exchange monte carlo method and application to spin glass simulations, J. Phys. Soc. Jpn.65, 1604 (1996)

    work page 1996

  39. [39]

    H. G. Katzgraber, S. Trebst, D. A. Huse, and M. Troyer, Feedback-optimized parallel tempering monte carlo, J. Stat. Mech.2006, P03018 (2006)

    work page 2006

  40. [40]

    Hamze, N

    F. Hamze, N. Dickson, and K. Karimi, Robust parameter selection for parallel tempering, Int. J. Mod. Phys. C21, 603 (2010)

    work page 2010

  41. [41]

    Baskaran, D

    G. Baskaran, D. Sen, and R. Shankar, Spin-skitaev model: Classical ground states, order from disorder, and exact correlation functions, Phys. Rev. B78, 115116 (2008)

    work page 2008

  42. [42]

    Chandra, K

    S. Chandra, K. Ramola, and D. Dhar, Classical heisen- berg spins on a hexagonal lattice with kitaev couplings, Phys. Rev. E82, 031113 (2010)

    work page 2010

  43. [43]

    Rousochatzakis, N

    I. Rousochatzakis, N. B. Perkins, Q. Luo, and H.-Y. Kee, Beyond kitaev physics in strong spin-orbit coupled mag- nets, Rep. Prog. Phys.87, 026502 (2024)

    work page 2024

  44. [44]

    Holstein and H

    T. Holstein and H. Primakoff, Field dependence of the intrinsic domain magnetization of a ferromagnet, Phys. Rev.58, 1098 (1940)

    work page 1940

  45. [45]

    Y. Q. Li, M. Ma, D. N. Shi, and F. C. Zhang, Su(4) theory for spin systems with orbital degeneracy, Phys. Rev. Lett.81, 3527 (1998)

    work page 1998

  46. [46]

    R. A. Muniz, Y. Kato, and C. D. Batista, Gen- eralized spin-wave theory: Application to the bilin- ear–biquadratic model, Prog. Theor. Exp. Phys.2014, 083I01 (2014)

    work page 2014

  47. [47]

    Joshi, M

    A. Joshi, M. Ma, F. Mila, D. N. Shi, and F. C. Zhang, Elementary excitations in magnetically ordered systems with orbital degeneracy, Phys. Rev. B60, 6584 (1999)

    work page 1999

  48. [48]

    Onufrieva, Low-temperature properties of spin sys- tems with tensor order parameters, Zh

    F. Onufrieva, Low-temperature properties of spin sys- tems with tensor order parameters, Zh. Eksp. Teor. Fiz 89, 2270 (1985)

    work page 1985

  49. [49]

    Colpa, Diagonalization of the quadratic boson hamil- tonian, Physica A93, 327 (1978)

    J. Colpa, Diagonalization of the quadratic boson hamil- tonian, Physica A93, 327 (1978)

    work page 1978

  50. [50]

    J. c. v. Chaloupka and G. Khaliullin, Hidden symmetries of the extended kitaev-heisenberg model: Implications for the honeycomb-lattice iridatesA 2iro3, Phys. Rev. B 92, 024413 (2015)

    work page 2015

  51. [51]

    Rastelli, L

    E. Rastelli, L. Reatto, and A. Tassi, Quantum fluctua- tions in helimagnets, J. Phys. C: Solid State Phys.18, 353 (1985)

    work page 1985

  52. [52]

    A. V. Chubukov, S. Sachdev, and T. Senthil, Large-s ex- pansion for quantum antiferromagnets on a triangular lattice, J. Phys. Condens. Matter6, 8891 (1994)

    work page 1994

  53. [53]

    A. L. Chernyshev and M. E. Zhitomirsky, Spin waves in a triangular lattice antiferromagnet: Decays, spec- trum renormalization, and singularities, Phys. Rev. B79, 144416 (2009)

    work page 2009

  54. [54]

    D. P. Arovas and A. Auerbach, Functional integral the- ories of low-dimensional quantum heisenberg models, Phys. Rev. B38, 316 (1988)

    work page 1988

  55. [55]

    Auerbach,Interacting electrons and quantum mag- netism(Springer Science & Business Media, 2012)

    A. Auerbach,Interacting electrons and quantum mag- netism(Springer Science & Business Media, 2012)

    work page 2012

  56. [56]

    Zhang, E

    S.-S. Zhang, E. A. Ghioldi, L. O. Manuel, A. E. Trumper, and C. D. Batista, Schwinger boson theory of ordered magnets, Phys. Rev. B105, 224404 (2022)

    work page 2022