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

arXiv:2603.08710 · cond-mat.str-el · Submitted 2026-03-09 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: We're looking at the paper: "Thermal and quantum fluctuations in extended Kitaev-Yao-Lee spin-orbital model".

Mira: Here's the gist of the paper. Thermal and quantum fluctuations in the extended Kitaev-Yao-Lee spin-orbital model investigates how thermal and quantum fluctuations give rise to disordered phases in this extended model, providing insight into their emergence Model Background The paper builds upon the spin-1/2 Kitaev model on a honeycomb lattice, using the Yao-Lee spin-orbital model which provides exactly solvable quantum spin liquids with potentially better stability against perturbations due to the additional degree of freedom The microscopic mechanism underlying the YaoLee 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 This extended model is derived when the Kugel-Khomskii interaction is included, leading to a broad region of a disordered phase named the nematic phase (NP) reminiscent of a spin-orbital liquid in numerical studies Classical Phase Diagram Analysis The classical Monte Carlo simulation (CMC) studies reveal that the NP exhibits macroscopic degeneracies arising from the spin-orbit degree of freedom and size-independent low energy costs, causing it to emerge at low temperatures due to thermal effects The zero-temperature phase diagram shows that when a = 0 and b = 0, the YL phase is found, while when a > 0 and b = 0, the NP emerges with an energy per bond given by epsilonNP = -a + 2/3 S squared At finite temperature, the Kitaev phase and the YL point can all extend into the fully ordered a > 0 and b > 0 region, as shown by Fig.

Kai: 1(d) Quantum Fluctuations via Generalized SWT The study employs a generalized spin wave theory (SWT) appropriate for the spin-orbital model to study quantum effects in the phase diagram This generalized SWT includes all cubic and quartic terms up to order Λ−1, which are important for general spin-orbital models since the cubic terms are present even at order Λ0 The analysis shows that the NP is stabilized by significant quantum fluctuations in the generalized SWT Role of Higher-Order Terms The inclusion of higher-order interaction terms, particularly cubic terms, significantly enhances quantum fluctuations and produces a phase boundary at large b values The cubic corrections lead to the gap closing at a much larger b value, suggesting the importance of the cubic terms for the extended nematic phase The quartic terms simply renormalize the non-interacting magnon energies and do not change the gap closing at b = 0 The generalized SWT results reveal strong quantum effects in the disordered phases of both the Kitaev and the NP regions The study suggests that materials with strong spin-orbital interactions, such as 3d7 or 3d9 compounds hosting one electron or hole in a two-orbital manifold and exhibiting weak or quenched spin–orbit coupling at the transition-metal sites, are promising candidates for realizing stable spin–orbital liquids The use of alternative theoretical frameworks, such as the SU(4) Schwinger boson approach to investigate nematic and other disordered phases in spin-orbital systems is an excellent project for future study Appendix A: Transition between Fully Ordered Phases The ground state energy of the antiferromagnetic (AFM) spin with antiferromagnetic orbital (AFO) (AFM×AFO) phase per bond is given by epsilon AFM×AFO = -S squared / (3a − 2/3 + b S squared × 1 + b S 2) Comparing this with the energy of the SS×AFO phase, epsilon SS×AFO = -S squared / (a + 2/3 + b S squared × 1 + b S 2), we find that the transition of the two ordered states at a = 2 is due to the competition between the KugelKhomskii term a(Si · Sj)(Ti · Tj) and the Yao-Lee term 2(Sγ i · Sγ j)(Ti · Tj), which is independent of b Appendix B: Classical Monte Carlo Simulation with Parallel Tempering Parallel tempering allows replicas to be simulated at higher temperatures, where the system can more easily escape metastable states and subsequently relax at lower temperatures on much shorter time scales than simulated annealing The key parameters affecting equilibration efficiency are the simulation temperatures of the replicas For the CMC simulations in this work, we initialize the feedback-optimized algorithm with 200 temperature steps, exponentially increasing from T = 0.001K˜ to T = 0.1K Appendix C: Generalized Spin Wave Theory When expanded to the order Λ−1, the transformation reads O 1 n = √Λb† n / (√b† n 1 - 1/2ΛX4 l=2 b† l bl + O(Λ−2)!, which generates both cubic and quartic terms The quartic terms are treated with the mean-field Hartree-Fock method, while the cubic terms are treated with second order perturbation to the non-interacting states The cubic corrections lead to the gap closing at a much larger b, suggesting the importance of the cubic terms for the extended nematic phase Appendix D: Static Structure Factor The NP has a maximum ST(q = M) but smaller S(q = M) and T(q = 0) that depend on the lattice size, as shown by Fig.

Mira: 3(b) with a = 0.5 and b = 0.

Paper summary: Kai: To put it plainly, they are looking at this extended model and showing how thermal noise and quantum uncertainty cause disorder to show up where you might expect perfect order.

Mira: The paper claims that these fluctuations favor the nematic phase over other disordered phases at low temperatures, which is a key finding for understanding the competition between different states in this system.

Lev: They use classical Monte Carlo simulations and a generalized spin wave theory appropriate for the spin-orbital model to explore this phase diagram at both zero and finite temperatures.

Kai: The specific results they highlight are that the nematic phase is stabilized by significant quantum fluctuations when using the generalized SWT, which includes those cubic and quartic terms up to order lambda minus one.

Mira: Furthermore, they find that including those higher-order interaction terms, particularly the cubic ones, produces a phase boundary at large 'b' values because those corrections lead to the gap closing at a much larger 'b' value.

Lev: The paper also shows that the zero-temperature phase diagram has two fully ordered phases denoted by AFM times AFO and SS times AFO when 'a' and 'b' are both greater than two, and that the transition between these two occurs at 'a equals two', regardless of the finite value of 'b' in the classical model <ref:2603.08710#pg1>.

Kai: And they also note something about the static structure factor: it shows a maximum ST at q equals M but a smaller S at q equals M and T at q equals zero that depends on the lattice size when 'a' is zero point five and 'b' is zero.

Mira: So, what this all means for us is that the existence and stability of this nematic phase are very sensitive to these fluctuations, especially the quantum ones described by those higher-order terms.

Lev: And they admit a limitation: their generalized SWT only goes up to order lambda minus one, so the results might not capture everything if there are even higher-order effects at play in real materials.

Conclusion: Kai: So, the main idea is that these thermal and quantum effects aren't just minor corrections; they are what give rise to this broad region of disorder, which they link to spin-orbital liquid behavior.

Mira: The authors argue that understanding this stabilization requires looking at the cubic terms in the generalized spin wave theory because those terms are responsible for the gap closing at larger 'b' values and thus define where we see this nematic phase.

Lev: For those of us building hardware, it’s important to know that these quantum effects are what determine how stable any spin-orbital liquid state might be on real devices, and the paper suggests materials like 3d7 or 3d9 compounds could be good candidates <ref:2603.08710#pg3>.

Kai: And they're suggesting that future work should look at alternative approaches, like the SU(four) Schwinger boson approach to investigate nematic phases in spin-orbital systems more deeply <ref:2603.08710#pg3>.

Mira: It really boils down to the fact that the competition between different interactions—like Kugel-Khomskii versus Yao-Lee—is what sets up these complex phase diagrams where disorder can exist.

Lev: So, while this paper gives us a good map of what's possible in theory, it points toward needing more rigorous methods to predict exactly how strong those fluctuations will be in a real system.

Department of Physics, University of Toronto · Canadian Institute for Advanced Research

cond-mat.str-el

Submitted: 2026-03-09

Updated: 2026-05-20

Comments: 10 pages, 5 figures

Journal ref: Phys. Rev. B 113, 174430 (2026)

DOI: 10.1103/dx3p-s4f3

License: http://creativecommons.org/licenses/by-nc-nd/4.0/

Importance score: 89/100

The gist: Thermal and quantum fluctuations in the extended Kitaev-Yao-Lee spin-orbital model investigates how thermal and quantum fluctuations give rise to disordered phases in this extended model, providing

Key concepts

Kitaev Model: The base model used is the Kitaev spin-1/2 model on a honeycomb lattice. This simple model describes quantum magnetism and is known for having exactly solvable quantum spin liquids. It serves as the foundation upon which the more complex extended version of this system is built.
Yao-Lee Spin-Orbital Model: This specific model extends the Kitaev framework by adding a spin-orbital degree of freedom. This addition allows for exactly solvable quantum spin liquids and provides better stability against external perturbations compared to simpler models, making it suitable for studying complex magnetic materials.
Nematic Phase (NP): The nematic phase is a disordered region found in the phase diagram of this extended model. It arises because of macroscopic degeneracies caused by the spin-orbit degree of freedom and low energy costs. This phase is reminiscent of a spin-orbital liquid observed in numerical studies.
Generalized Spin Wave Theory (SWT): This advanced theoretical tool is used to analyze quantum effects in the model. It goes beyond simpler theories by including cubic and quartic terms up to a certain order, which are essential for accurately describing the quantum fluctuations that stabilize the disordered phases.

Terminology

Quantum fluctuations — 24× in this paper · explained in 3 episodes
Kitaev Model — 12× in this paper · explained in 2 episodes
Degeneracy — 8× in this paper · explained in 4 episodes
Degeneracies — 7× in this paper · explained in 2 episodes
Phase Transition — 6× in this paper · explained in 3 episodes
Generator — 5× in this paper · explained in 2 episodes
Mott Insulator — 4× in this paper · explained in 2 episodes
Anyons — 2× in this paper · explained in 2 episodes
Frustration — 2× in this paper · explained in 3 episodes
Majorana Fermions — 2× in this paper · explained in 2 episodes
Phase Transitions — 2× in this paper · explained in 2 episodes
Accuracy — 1× in this paper · explained in 1 episode
Anisotropy — 1× in this paper · explained in 4 episodes
Core — 1× in this paper · explained in 2 episodes
Exchange Interaction — 1× in this paper · explained in 2 episodes
Monte Carlo Sampling — 1× in this paper · explained in 2 episodes

Summary

Thermal and quantum fluctuations in the extended Kitaev-Yao-Lee spin-orbital model investigates how thermal and quantum fluctuations give rise to disordered phases in this extended model, providing insight into their emergence

Model Background

The paper builds upon the spin-1/2 Kitaev model on a honeycomb lattice, using the Yao-Lee spin-orbital model which provides exactly solvable quantum spin liquids with potentially better stability against perturbations due to the additional degree of freedom The microscopic mechanism underlying the YaoLee 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 This extended model is derived when the Kugel-Khomskii interaction is included, leading to a broad region of a disordered phase named the nematic phase (NP) reminiscent of a spin-orbital liquid in numerical studies

Classical Phase Diagram Analysis

The classical Monte Carlo simulation (CMC) studies reveal that the NP exhibits macroscopic degeneracies arising from the spin-orbit degree of freedom and size-independent low energy costs, causing it to emerge at low temperatures due to thermal effects The zero-temperature phase diagram shows that when a = 0 and b = 0, the YL phase is found, while when a > 0 and b = 0, the NP emerges with an energy per bond given by epsilonNP = -a + 2/3 S squared At finite temperature, the Kitaev phase and the YL point can all extend into the fully ordered a > 0 and b > 0 region, as shown by Fig. 1(d)

Quantum Fluctuations via Generalized SWT

The study employs a generalized spin wave theory (SWT) appropriate for the spin-orbital model to study quantum effects in the phase diagram This generalized SWT includes all cubic and quartic terms up to order Λ−1, which are important for general spin-orbital models since the cubic terms are present even at order Λ0 The analysis shows that the NP is stabilized by significant quantum fluctuations in the generalized SWT

Role of Higher-Order Terms

The inclusion of higher-order interaction terms, particularly cubic terms, significantly enhances quantum fluctuations and produces a phase boundary at large b values The cubic corrections lead to the gap closing at a much larger b value, suggesting the importance of the cubic terms for the extended nematic phase The quartic terms simply renormalize the non-interacting magnon energies and do not change the gap closing at b = 0

Conclusion and Outlook

The generalized SWT results reveal strong quantum effects in the disordered phases of both the Kitaev and the NP regions The study suggests that materials with strong spin-orbital interactions, such as 3d7 or 3d9 compounds hosting one electron or hole in a two-orbital manifold and exhibiting weak or quenched spin–orbit coupling at the transition-metal sites, are promising candidates for realizing stable spin–orbital liquids The use of alternative theoretical frameworks, such as the SU(4) Schwinger boson approach to investigate nematic and other disordered phases in spin-orbital systems is an excellent project for future study

Appendix A: Transition between Fully Ordered Phases

The ground state energy of the antiferromagnetic (AFM) spin with antiferromagnetic orbital (AFO) (AFM×AFO) phase per bond is given by epsilon AFM×AFO = -S squared / (3a − 2/3 + b S squared × 1 + b S 2) Comparing this with the energy of the SS×AFO phase, epsilon SS×AFO = -S squared / (a + 2/3 + b S squared × 1 + b S 2), we find that the transition of the two ordered states at a = 2 is due to the competition between the KugelKhomskii term a(Si · Sj)(Ti · Tj) and the Yao-Lee term 2(Sγ i · Sγ j)(Ti · Tj), which is independent of b

Appendix B: Classical Monte Carlo Simulation with Parallel Tempering

Parallel tempering allows replicas to be simulated at higher temperatures, where the system can more easily escape metastable states and subsequently relax at lower temperatures on much shorter time scales than simulated annealing The key parameters affecting equilibration efficiency are the simulation temperatures of the replicas For the CMC simulations in this work, we initialize the feedback-optimized algorithm with 200 temperature steps, exponentially increasing from T = 0.001K˜ to T = 0.1K

Appendix C: Generalized Spin Wave Theory

When expanded to the order Λ−1, the transformation reads O 1 n = √Λb† n / (√b† n 1 - 1/2ΛX4 l=2 b† l bl + O(Λ−2)!, which generates both cubic and quartic terms The quartic terms are treated with the mean-field Hartree-Fock method, while the cubic terms are treated with second order perturbation to the non-interacting states The cubic corrections lead to the gap closing at a much larger b, suggesting the importance of the cubic terms for the extended nematic phase

Appendix D: Static Structure Factor

The NP has a maximum ST(q = M) but smaller S(q = M) and T(q = 0) that depend on the lattice size, as shown by Fig. 3(b) with a = 0.5 and b = 0.

Improvements for AI systems

Quantum phase diagram prediction for spin-orbital liquids: The improved system can predict which disordered phases, such as nematic phase (NP) or Kitaev phase, are stable at specific temperature and interaction parameter values by combining classical Monte Carlo simulations with generalized spin wave theory.
Enhanced understanding of quantum fluctuations in disordered states: The system can quantify how the additional spin-orbital degree of freedom gives rise to strong thermal and quantum fluctuations in spin-orbital models, providing a deeper microscopic explanation for the emergence of these phases.
Accurate determination of phase boundaries using higher-order corrections: The system can calculate phase boundaries by incorporating cubic and quartic terms up to order Lambda minus one in generalized spin wave theory, which is necessary because the cubic terms are present even at order Lambda zero.
Distinguishing thermal selection from quantum stabilization: The improved system can differentiate between phases stabilized by thermal effects and those stabilized by quantum fluctuations, specifically showing that the NP exhibits macroscopic degeneracies arising from the spin-orbit degree of freedom and size-independent low energy costs, so it emerges at low temperatures due to thermal effects versus its stabilization by the significant quantum fluctuations in the generalized SWT.
Material candidate screening: The system can suggest promising material candidates for realizing stable spin-orbital liquids by identifying those with strong spinorbital interactions, such as 3d7 or 3d9 compounds hosting one electron or hole in a two-orbital manifold and exhibiting weak or quenched spin–orbit coupling at the transition-metal sites.

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.