Chiral spin liquid instability of the Kitaev honeycomb model with crystallographic defects
Pith reviewed 2026-05-17 04:59 UTC · model grok-4.3
The pith
Finite density of lattice defects drives a phase transition to a chiral quantum spin liquid in the Kitaev model at temperature set by defect density.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Introducing a finite defect density n_d ≈ 10^{-4}–10^{-2} produces a true phase transition with a sizeable T_c ≈ 2 n_d in units of the Kitaev exchange. The resulting non-Abelian chiral quantum spin liquid exhibits scalar spin chirality and electron orbital magnetization which peak near lattice defects. This disorder-driven instability relies on an emergent long range ferromagnetic interaction r^{-γ} (γ ≈ 2.7) between defect chiralities, mediated by the nearly-gapless fermions.
What carries the argument
Emergent long-range ferromagnetic interaction r^{-γ} (γ≈2.7) between defect chiralities, mediated by the nearly-gapless fermions of the Kitaev spin liquid.
If this is right
- Scalar spin chirality becomes finite and peaks near each defect site.
- Electron orbital magnetization appears and is largest in the vicinity of the defects.
- The chiral phase provides a concrete route to topology generation inside Dirac cones that carry fluctuating mass terms.
- The transition temperature remains sizeable even when the defect density is as low as one part in ten thousand.
Where Pith is reading between the lines
- Similar defect-induced ordering may appear in other gapless Dirac spin liquids once odd-plaquette defects are introduced.
- The same mechanism could be tested by deliberately engineering controlled densities of Stone-Wales defects in candidate Kitaev materials.
- The r^{-2.7} decay suggests that the effective interaction is marginally long-ranged, which may place the transition in a distinct universality class from short-range Ising models.
Load-bearing premise
The long-range ferromagnetic interaction between defect chiralities remains strong enough at low defect densities to produce a genuine thermodynamic phase transition rather than a crossover.
What would settle it
A Monte Carlo simulation or finite-size scaling analysis that finds no diverging susceptibility or order parameter at the predicted T_c ≈ 2 n_d when the long-range interaction is artificially removed or cut off.
Figures
read the original abstract
We study the spin-1/2 Kitaev honeycomb gapless spin liquid in the presence of Stone-Wales-type local lattice defects with odd-sided plaquettes. While the clean Kitaev model has no finite-temperature phase transitions, we find that introducing a finite defect density $n_d\approx 10^{-4}$--$10^{-2}$ produces a true phase transition with a sizeable $T_c \approx 2 n_d$ in units of the Kitaev exchange. The resulting non-Abelian chiral quantum spin liquid exhibits scalar spin chirality and electron orbital magnetization which peak near lattice defects. This disorder-driven instability relies on an emergent long range ferromagnetic interaction $r^{-\gamma}$ ($\gamma \approx 2.7$) between defect chiralities, mediated by the nearly-gapless fermions, with implications for topology generation in Dirac cones with fluctuating mass terms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines the gapless Kitaev honeycomb spin liquid subject to a dilute density of Stone-Wales crystallographic defects that introduce odd-sided plaquettes. It reports that these defects generate an emergent long-range ferromagnetic interaction between local chiralities, mediated by the nearly gapless Majorana fermions and decaying as r^{-γ} with γ≈2.7, which in turn drives a finite-temperature transition to a non-Abelian chiral quantum spin liquid at T_c≈2 n_d (in units of the Kitaev exchange) for defect densities n_d≈10^{-4}–10^{-2}. The resulting phase is characterized by scalar spin chirality and orbital magnetization that are peaked near the defects.
Significance. If the central numerical claim holds, the work identifies a concrete, disorder-driven route to chiral topological order in an otherwise gapless Kitaev liquid, with direct implications for materials realizations that inevitably contain lattice defects. The fermion-mediated power-law interaction provides a mechanism that can be tested in other Dirac or Majorana systems with fluctuating mass terms.
major comments (2)
- [Numerical evidence for the phase transition and interaction extraction] The headline result—that a true thermodynamic transition occurs at the quoted densities with T_c≈2 n_d—rests on the assertion that the derived r^{-γ} (γ≈2.7) interaction is sufficiently long-ranged and strong to overcome 2D thermal fluctuations. No finite-size scaling of the order parameter, Binder-cumulant crossings, or explicit extrapolation to the thermodynamic limit is described; at mean defect separations of 10–100 lattice spacings this leaves open the possibility that the observed signatures are crossovers driven by rare defect pairs rather than a bulk transition.
- [Derivation of the defect–defect interaction J(r)] The exponent γ≈2.7 and the linear relation T_c≈2 n_d appear to be extracted from the same set of simulations used to identify the transition itself. This introduces a circularity that must be broken by an independent, parameter-free derivation of the interaction prefactor or by a separate scaling analysis that does not presuppose the transition temperature.
minor comments (2)
- [Introduction and model definition] The abstract and introduction state that the Majorana spectrum remains “nearly gapless,” but the manuscript should quantify the gap opened by the defects (e.g., as a function of n_d) to justify the long-range character of the mediated interaction.
- [Figures showing scalar spin chirality and orbital magnetization] Figure captions and axis labels for the chirality and magnetization profiles should explicitly state the system size and boundary conditions used, as these quantities are reported to peak near defects.
Simulated Author's Rebuttal
We thank the referee for the careful reading and valuable comments on our work. We address the two major concerns point by point below, providing clarifications on our numerical procedures and indicating the revisions we will implement to strengthen the evidence for a thermodynamic transition.
read point-by-point responses
-
Referee: No finite-size scaling of the order parameter, Binder-cumulant crossings, or explicit extrapolation to the thermodynamic limit is described; at mean defect separations of 10–100 lattice spacings this leaves open the possibility that the observed signatures are crossovers driven by rare defect pairs rather than a bulk transition.
Authors: We agree that explicit finite-size scaling and Binder cumulant analysis are important to confirm a true thermodynamic transition rather than a crossover. In the revised manuscript we will add Binder cumulant crossings for the scalar spin chirality as a function of temperature for multiple system sizes at fixed defect densities, together with data collapse of the order parameter. These additions will demonstrate that the transition temperature remains stable with increasing system size and is not dominated by rare defect pairs. revision: yes
-
Referee: The exponent γ≈2.7 and the linear relation T_c≈2 n_d appear to be extracted from the same set of simulations used to identify the transition itself. This introduces a circularity that must be broken by an independent, parameter-free derivation of the interaction prefactor or by a separate scaling analysis.
Authors: The exponent γ≈2.7 is obtained from an independent zero-temperature calculation of the effective interaction J(r) between two isolated defects. This uses the static Majorana fermion susceptibility of the clean Kitaev model via second-order perturbation theory in the defect-induced mass terms and does not involve the finite-temperature many-defect simulations. The relation T_c≈2 n_d follows from a subsequent mean-field estimate using this J(r). We will expand the methods section to present this derivation explicitly and separate it from the finite-T results. revision: yes
Circularity Check
Numerical fit of γ≈2.7 and T_c≈2 n_d from defect simulations introduces moderate circularity in transition claim
specific steps
-
fitted input called prediction
[Abstract]
"introducing a finite defect density n_d≈10^{-4}--10^{-2} produces a true phase transition with a sizeable T_c ≈ 2 n_d in units of the Kitaev exchange. [...] emergent long range ferromagnetic interaction r^{-γ} (γ ≈ 2.7) between defect chiralities, mediated by the nearly-gapless fermions"
The quoted numerical values for the interaction exponent γ≈2.7 and the transition temperature scaling T_c≈2 n_d are extracted from the same finite-density defect simulations that are used to assert the existence of the thermodynamic transition; the 'prediction' of a true phase transition is therefore statistically tied to the input data used to measure those quantities rather than independently derived.
full rationale
The paper derives the long-range interaction and reports a phase transition directly from numerical data on defected lattices at the quoted densities. While the underlying Majorana-mediated mechanism has independent content, the specific exponent and linear T_c(n_d) relation are obtained from the same simulations used to identify the transition, without separate analytic derivation or external validation. This creates partial circularity but does not reduce the entire claim to a tautology. No self-citation load-bearing or self-definitional steps were identified in the provided text.
Axiom & Free-Parameter Ledger
free parameters (2)
- defect density range n_d
- interaction exponent γ
axioms (1)
- domain assumption Kitaev honeycomb Hamiltonian supports a gapless Majorana spin liquid in the clean limit
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
emergent long range ferromagnetic interaction r^{-γ} (γ≈2.7) between defect chiralities, mediated by the nearly-gapless fermions
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
T_c ≈ 2 n_d in units of the Kitaev exchange
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 2 Pith papers
-
Generation of chirality and orbital magnetization by Stone-Wales-type lattice defects in the Kitaev spin liquid
Stone-Wales defects in the gapless Kitaev spin liquid generate net chirality that opens a topological gap of 11 n_d and drives a finite-temperature transition to a chiral spin liquid via long-range Ising interactions.
-
Impurity quadrupole moments as local probes of flux sectors in the Kitaev spin liquid
Impurity quadrupole moments exhibit discontinuous jumps at flux sector transitions in the Kitaev spin liquid, serving as a local probe of flux configurations.
Reference graph
Works this paper leans on
-
[1]
A. J. Willans, J. T. Chalker, and R. Moessner, Disorder in a quantum spin liquid: Flux binding and local moment formation, Physical Review Letters104, 237203 (2010)
work page 2010
-
[2]
A. J. Willans, J. T. Chalker, and R. Moessner, Site dilu- tion in the Kitaev honeycomb model, Physical Review B 84, 115146 (2011)
work page 2011
-
[3]
G. J. Sreejith, S. Bhattacharjee, and R. Moessner, Va- cancies in Kitaev quantum spin liquids on the three- dimensional hyperhoneycomb lattice, Physical Review B 93, 064433 (2016)
work page 2016
-
[4]
J. Nasu and Y. Motome, Thermodynamic and transport properties in disordered Kitaev models, Physical Review B102, 054437 (2020)
work page 2020
- [5]
-
[6]
K. Dhochak, R. Shankar, and V. Tripathi, Magnetic Im- purities in the Honeycomb Kitaev Model, Physical Re- view Letters105, 117201 (2010)
work page 2010
- [7]
-
[8]
W.-H. Kao, J. Knolle, G. B. Hal´ asz, R. Moessner, and N. B. Perkins, Vacancy-induced low-energy density of states in the kitaev spin liquid, Physical Review X11, 011034 (2021)
work page 2021
-
[9]
A. Singhania, J. Van Den Brink, and S. Nishimoto, Dis- order effects in the Kitaev-Heisenberg model, Physical Review Research5, 023009 (2023)
work page 2023
-
[10]
F. Zschocke and M. Vojta, Physical states and finite-size effects in Kitaev’s honeycomb model: Bond disorder, spin excitations, and NMR line shape, Physical Review B92, 014403 (2015)
work page 2015
- [11]
-
[12]
G. Cassella, P. D’Ornellas, T. Hodson, W. M. H. Na- tori, and J. Knolle, An exact chiral amorphous spin liq- uid, Nature Commun.14, 6663 (2023), arXiv:2208.08246 [cond-mat.str-el]
-
[13]
A. G. Grushin and C. Repellin, Amorphous and poly- crystalline routes toward a chiral spin liquid, Physical Review Letters130, 186702 (2023)
work page 2023
-
[14]
O. Petrova, P. Mellado, and O. Tchernyshyov, Unpaired Majorana modes on dislocations and string defects in Ki- taev’s honeycomb model, Physical Review B90, 134404 (2014)
work page 2014
- [15]
-
[16]
G. B. Hal´ asz, J. T. Chalker, and R. Moessner, Doping a topological quantum spin liquid: Slow holes in the Kitaev honeycomb model, Physical Review B90, 035145 (2014)
work page 2014
-
[17]
G. B. Hal´ asz and J. T. Chalker, Coherent hole propaga- tion in an exactly solvable gapless spin liquid, Physical Review B94, 235105 (2016)
work page 2016
-
[18]
V. Lahtinen, A. W. W. Ludwig, and S. Trebst, Perturbed vortex lattices and the stability of nucleated topological phases, Physical Review B89, 085121 (2014)
work page 2014
-
[19]
M. Udagawa, Vison-Majorana complex zero-energy res- onance in the Kitaev spin liquid, Physical Review B98, 220404 (2018)
work page 2018
- [20]
- [21]
- [22]
-
[23]
J. Nasu and Y. Motome, Spin dynamics in the Kitaev model with disorder: Quantum Monte Carlo study of dynamical spin structure factor, magnetic susceptibil- ity, and NMR relaxation rate, Physical Review B104, 035116 (2021)
work page 2021
-
[24]
L. R. D. Freitas and R. G. Pereira, Gapless excitations in non-Abelian Kitaev spin liquids with line defects, Physi- cal Review B105, L041104 (2022)
work page 2022
-
[25]
V. Dantas and E. C. Andrade, Disorder, Low-Energy Ex- citations, and Topology in the Kitaev Spin Liquid, Phys- ical Review Letters129, 037204 (2022)
work page 2022
-
[26]
X.-Y. Song, Y.-Z. You, and L. Balents, Low-energy spin dynamics of the honeycomb spin liquid beyond the kitaev limit, Physical Review Letters117, 037209 (2016)
work page 2016
-
[27]
I. Yatsuta and D. F. Mross, Vacancies in generic Kitaev spin liquids, arXiv:2312.00147 (2023), arXiv:2312.00147 [cond-mat]
- [28]
-
[29]
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
-
[30]
G. Jackeli and G. Khaliullin, Mott insulators in the strong spin-orbit coupling limit: From heisenberg to a quantum compass and kitaev models, Physical Review Letters102, 017205 (2009)
work page 2009
-
[31]
M. Hermanns, I. Kimchi, and J. Knolle, Physics of the Kitaev Model: Fractionalization, Dynamic Correlations, and Material Connections, Annual Review of Condensed Matter Physics9, 17 (2018)
work page 2018
- [32]
-
[33]
S. Trebst and C. Hickey, Kitaev materials, Physics Re- ports Kitaev Materials,950, 1 (2022)
work page 2022
-
[34]
L. Savary and L. Balents, Quantum spin liquids: A re- view, Reports on Progress in Physics80, 016502 (2016)
work page 2016
- [35]
-
[36]
P. d’Ornellas and J. Knolle, Kitaev-Heisenberg model on the star lattice: From chiral Majorana fermions to chiral triplons, Physical Review B109, 094421 (2024)
work page 2024
-
[37]
V. Peri, S. Ok, S. S. Tsirkin, T. Neupert, G. Baskaran, M. Greiter, R. Moessner, and R. Thomale, Non-Abelian chiral spin liquid on a simple non-Archimedean lattice, Physical Review B101, 041114 (2020)
work page 2020
-
[38]
F. Borhani, A. Seth, and I. Kimchi, Real-space chirality from crystalline topological defects in the Kitaev spin liquid, npj Quantum Materials10, 1 (2025)
work page 2025
-
[39]
A. J. Stone and D. J. Wales, Theoretical studies of icosa- hedral C60 and some related species, Chemical Physics Letters128, 501 (1986)
work page 1986
-
[40]
M. Vozmediano, M. Katsnelson, and F. Guinea, Gauge fields in graphene, Physics Reports496, 109 (2010)
work page 2010
-
[41]
P. Kot, J. Parnell, S. Habibian, C. Straßer, P. M. Os- trovsky, and C. R. Ast, Band dispersion of graphene with structural defects, Physical Review B101, 235116 (2020)
work page 2020
-
[42]
A. Seth, F. Borhani, and I. Kimchi, Generation of chiral- ity and orbital magnetization by Stone-Wales-type lattice defects in the Kitaev spin liquid (2025)
work page 2025
-
[43]
E. H. Lieb, Flux Phase of the Half-Filled Band, Physical Review Letters73, 2158 (1994)
work page 1994
-
[44]
This limit is not preserved by the low energy projection, requiring us to keep track of the sense of rotationν
-
[45]
F. D. M. Haldane, Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the ”Parity Anomaly”, Physical Review Letters61, 2015 (1988)
work page 2015
- [46]
- [47]
-
[48]
T. A. Loring and M. B. Hastings, Disordered topological insulators via C*-algebras, Europhysics Letters92, 67004 (2011)
work page 2011
-
[49]
M. B. Hastings and T. A. Loring, Almost commuting matrices, localized Wannier functions, and the quantum Hall effect, Journal of Mathematical Physics51, 015214 (2010)
work page 2010
-
[50]
R. Shindou and N. Nagaosa, Orbital Ferromagnetism and Anomalous Hall Effect in Antiferromagnets on the Dis- torted fcc Lattice, Physical Review Letters87, 116801 (2001)
work page 2001
-
[51]
O. I. Motrunich, Orbital magnetic field effects in spin liquid with spinon Fermi sea: Possible application to k- (ET)2Cu2(CN)3, Physical Review B73, 155115 (2006)
work page 2006
-
[52]
L. N. Bulaevskii, C. D. Batista, M. V. Mostovoy, and D. I. Khomskii, Electronic orbital currents and polarization in Mott insulators, Physical Review B78, 024402 (2008)
work page 2008
-
[53]
We also note that if the system is doped with itiner- ant electrons then SSC also produces Hall conductivity at zero frequency, through multiple mechanisms, includ- ing the topological Hall effect well known in skyrmion lattices [67] as well as anomalous Hall effect induced by spin chirality skew scattering [68]
-
[54]
Kogan, RKKY interaction in graphene, Physical Re- view B84, 115119 (2011)
E. Kogan, RKKY interaction in graphene, Physical Re- view B84, 115119 (2011)
work page 2011
-
[55]
Smoothly cutting off thiskintegral byn 1/2 d can produce an approximater −γ form with 2< γ <3
The observation thatγis suppressed rather than en- hanced can be interpreted by considering the graphene r−3 RKKY as arising from the energy or equivalentlyk integral of ther −2 spin LDOS for energies near the Dirac point [69]. Smoothly cutting off thiskintegral byn 1/2 d can produce an approximater −γ form with 2< γ <3
-
[56]
M. E. Fisher, S.-k. Ma, and B. G. Nickel, Critical Expo- nents for Long-Range Interactions, Physical Review Let- ters29, 917 (1972)
work page 1972
-
[57]
M. Aizenman and R. Fern´ andez, Critical exponents for long-range interactions, Letters in Mathematical Physics 16, 39 (1988)
work page 1988
- [58]
- [59]
-
[60]
E. Koller, V. Leeb, N. B. Perkins, and J. Knolle, Raman Circular Dichroism and Quantum Geometry of Chiral Quantum Spin Liquids (2025), arXiv:2503.14091 [cond- mat]
work page internal anchor Pith review arXiv 2025
- [61]
-
[62]
D. Vasyukov, Y. Anahory, L. Embon, D. Halbertal, J. Cuppens, L. Neeman, A. Finkler, Y. Segev, Y. Mya- soedov, M. L. Rappaport, M. E. Huber, and E. Zeldov, A scanning superconducting quantum interference device with single electron spin sensitivity, Nature Nanotechnol- ogy8, 639 (2013)
work page 2013
-
[63]
A. Finkler, Y. Segev, Y. Myasoedov, M. L. Rappaport, L. Ne’eman, D. Vasyukov, E. Zeldov, M. E. Huber, J. Martin, and A. Yacoby, Self-Aligned Nanoscale SQUID on a Tip, Nano Letters10, 1046 (2010)
work page 2010
-
[64]
N. J. McLaughlin, S. Li, J. A. Brock, S. Zhang, H. Lu, M. Huang, Y. Xiao, J. Zhou, Y. Tserkovnyak, E. E. Fullerton, H. Wang, and C. R. Du, Local Control of a Sin- gle Nitrogen-Vacancy Center by Nanoscale Engineered Magnetic Domain Wall Motion, ACS Nano17, 25689 (2023)
work page 2023
-
[65]
C. Du, T. van der Sar, T. X. Zhou, P. Upadhyaya, F. Casola, H. Zhang, M. C. Onbasli, C. A. Ross, R. L. Walsworth, Y. Tserkovnyak, and A. Yacoby, Control and local measurement of the spin chemical potential in a magnetic insulator, Science357, 195 (2017)
work page 2017
-
[66]
S. Grover, M. Bocarsly, A. Uri, P. Stepanov, G. Di Bat- tista, I. Roy, J. Xiao, A. Y. Meltzer, Y. Myasoedov, K. Pareek, K. Watanabe, T. Taniguchi, B. Yan, A. Stern, E. Berg, D. K. Efetov, and E. Zeldov, Chern mosaic and Berry-curvature magnetism in magic-angle graphene, Na- ture Physics18, 885 (2022)
work page 2022
-
[67]
T. Kurumaji, T. Nakajima, M. Hirschberger, A. Kikkawa, Y. Yamasaki, H. Sagayama, H. Nakao, Y. Taguchi, T.-h. Arima, and Y. Tokura, Skyrmion lattice with a giant topological Hall effect in a frustrated triangular-lattice magnet, Science365, 914 (2019), https://www.science.org/doi/pdf/10.1126/science.aau0968
-
[68]
H. Ishizuka and N. Nagaosa, Spin chirality induced skew scattering and anomalous Hall effect in chi- ral magnets, Science Advances4, eaap9962 (2018), https://www.science.org/doi/pdf/10.1126/sciadv.aap9962
- [69]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.