pith. sign in

arxiv: 2512.13490 · v2 · submitted 2025-12-15 · ❄️ cond-mat.str-el · cond-mat.dis-nn· cond-mat.mes-hall

Generation of chirality and orbital magnetization by Stone-Wales-type lattice defects in the Kitaev spin liquid

Pith reviewed 2026-05-16 22:10 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.dis-nncond-mat.mes-hall
keywords Kitaev spin liquidStone-Wales defectschiralityorbital magnetizationchiral spin liquidtopological gapIsing interactionsdefect density
0
0 comments X

The pith

Stone-Wales defects in the Kitaev spin liquid generate chirality that opens a topological gap of 11 times defect density and drives a finite-temperature transition to the chiral spin liquid.

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

The paper identifies a solvable Stone-Wales defect in the gapless Kitaev honeycomb spin liquid consisting of a 90-degree bond rotation that preserves bond labels and creates odd-sided plaquettes with plus or minus pi over 2 fluxes. An isolated defect supports time-reversal pairs of flux configurations carrying large net chirality, which appears in the Majorana local Chern marker, scalar spin chirality, and electronic orbital magnetization. At finite defect density, T-matrix analysis combined with numerics shows these chiralities produce a topological gap scaling as 11 times the density that protects a Chern number of plus or minus 1. The same defects generate emergent ferromagnetic Ising interactions decaying as a power law with exponent between 2 and 3, which stabilize a chiral spin liquid phase with a transition temperature proportional to defect density.

Core claim

Stone-Wales-type defects consisting of local 90-degree bond rotations in the Kitaev honeycomb spin liquid preserve exact solvability while introducing plus or minus pi over 2 fluxes on odd plaquettes; an isolated defect hosts a time-reversal pair of ground states with large net chirality that generates orbital magnetization, and at finite density these chiral defects open a topological gap of 11 n_d protecting Chern number plus or minus 1 while producing long-range ferromagnetic Ising interactions with exponent gamma between 2 and 3 that drive a finite-temperature transition into the chiral spin liquid with T_c proportional to n_d.

What carries the argument

The Stone-Wales defect, defined as a 90-degree bond rotation that preserves Kitaev bond labels for edge-sharing octahedra and thereby creates solvable odd-sided plaquettes carrying plus or minus pi over 2 fluxes, which hosts chiral ground-state configurations whose interactions are analyzed via T-matrix methods.

If this is right

  • Defect chiralities open a topological gap of 11 n_d that protects a Chern number of plus or minus 1.
  • Emergent ferromagnetic long-range Ising interactions with exponent gamma between 2 and 3 drive a finite-temperature transition into the chiral spin liquid.
  • The transition temperature T_c scales proportionally with defect density n_d and diverges as gamma approaches 2 from above.
  • Additional solvable impurity potentials can reduce gamma below 2.3 and correspondingly raise T_c.
  • The mechanism applies directly to 2D Dirac cone systems containing a finite density of fluctuating Ising magnetic impurities.

Where Pith is reading between the lines

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

  • Varying defect density in candidate Kitaev materials would provide a direct experimental knob for tuning the chiral transition temperature.
  • The orbital magnetization generated by the chiral defects offers a potential signature detectable through local magnetic probes or torque magnetometry.
  • The same defect-chirality mechanism could be tested in other exactly solvable spin models or in graphene-like Dirac systems with magnetic impurities.
  • Finite-density numerics at higher defect concentrations might reveal whether screening or clustering effects eventually cut off the power-law interactions.

Load-bearing premise

The T-matrix analysis and numerics at finite defect density correctly capture the effective long-range interactions between defect chiralities without significant higher-order corrections or screening effects that would alter the interaction exponent gamma or the gap coefficient 11.

What would settle it

A measurement showing that the spin-liquid energy gap scales linearly with defect density n_d at low densities, or that the transition temperature into the chiral phase is proportional to n_d and diverges as the interaction exponent approaches 2.

Figures

Figures reproduced from arXiv: 2512.13490 by Arnab Seth, Fay Borhani, Itamar Kimchi.

Figure 1
Figure 1. Figure 1: bottom right panel, gives the same z direction Kitaev interaction, in addition to other perturbations, thereby also preserving solvability for the Kitaev terms. These considerations require the 900 rotation to be in a “(0,0,1)” plane that is tilted relative to the “(1,1,1)” honeycomb lattice plane and imply a local out-of-plane distortion for SW defects. Such distortion already arises for Stone-Wales defec… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: shows the distribution of SSC around the SW for Lieb and PT-flux configurations. The contribution is large near the SW and decays away from the defect. Sim￾ilar to the local marker result, the SSC shows the symme￾try of a monopole, and a dipole, for the Lieb and PT-flux states, respectively. For Lieb states, SSC is contributed only with a single sign. For PT, it has both positive and negative contributions… view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: An additional advantage of this comparison is that the consideration of finite defect density automatically in￾cludes any nonlinearities coming from interactions be￾yond two-body. We will show that the power law Ising model described by 1 still remains a good approximation [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12 [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13 [PITH_FULL_IMAGE:figures/full_fig_p016_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Computation of the flux energies in the gapped Ki [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
read the original abstract

In this work we extend our study of the effect of certain crystallographic defects on the spin-1/2 Kitaev honeycomb spin liquid (arXiv:2511.19409), focusing on its gapless phase and contrasting with the gapped phase. We identify a Stone-Wales (SW) local defect consisting of a 90$^\circ$ bond rotation that preserves Kitaev bond labels for edge-sharing octahedra and thereby enables exact solvability. These SW-type defects involve odd-sided plaquettes with $\pm \pi/2$ fluxes, but can be created locally. An isolated defect hosts a time-reversal pair of ground-state flux configurations with large net chirality. Certain excitations are also chiral. The chirality manifests in Majorana local Chern marker and in scalar spin chirality, producing electronic orbital magnetization. T-matrix analysis and numerics at finite defect density $n_d$ show that defect chiralities generate a topological gap of $11 n_d$ protecting a Chern number $C=\pm 1$. Emergent ferromagnetic long range Ising interactions $r^{-\gamma}$ with $2<\gamma < 3$ between defect chiralities lead to a finite temperature $T_c$ phase transition into the chiral spin liquid. The $T_c$ is proportional to $n_d$ and diverges when $\gamma\rightarrow 2$. We also consider additional solvable impurity potentials and find that $\gamma$ can be reduced to below $2.3$ and correspondingly enhance $T_c$. Our results offer applications to 2D Dirac cone systems with a finite density of fluctuating Ising magnetic impurities and to identifying spin liquids with lattice defects.

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

3 major / 2 minor

Summary. The paper extends prior work on defects in the Kitaev honeycomb spin liquid to its gapless phase. It identifies Stone-Wales defects (90° bond rotations preserving Kitaev labels) that remain exactly solvable, host time-reversal pairs of chiral ground-state flux configurations with net chirality, and produce orbital magnetization via the Majorana local Chern marker and scalar spin chirality. At finite defect density n_d, T-matrix analysis combined with numerics is used to show that these chiralities open a topological gap of 11 n_d protecting Chern number C=±1; emergent ferromagnetic Ising interactions decaying as r^{-γ} with 2<γ<3 then drive a finite-T_c transition into the chiral spin liquid, with T_c proportional to n_d (and diverging as γ→2). Additional solvable impurity potentials are shown to reduce γ below 2.3 and enhance T_c.

Significance. If the central claims hold, the work supplies an exactly solvable route to defect-induced chirality and orbital magnetization in a gapless Kitaev spin liquid, together with concrete scaling predictions (gap linear in n_d, T_c∝n_d) that could be tested in 2D Dirac systems with magnetic impurities. The exact solvability for isolated defects and the parameter-free character of the isolated-defect chirality are genuine strengths; the finite-density T-matrix plus numerics, if fully documented, would constitute a reproducible prediction for the effective interactions.

major comments (3)
  1. [T-matrix analysis] T-matrix analysis section: the coefficient 11 in the reported topological gap Δ=11 n_d is stated as an output of the T-matrix calculation, yet no explicit T-matrix equations, matrix elements, or intermediate steps are shown; without these it is impossible to verify that the prefactor is independent of the chosen defect configurations or finite-size cutoffs.
  2. [Finite-density numerics] Finite-density numerics section: the interaction exponent range 2<γ<3 is extracted from numerics, but the manuscript provides neither error bars, finite-size extrapolation of the fitted γ, nor explicit checks for higher-order multi-defect scattering or screening of the Majorana wavefunctions that could renormalize either γ or the gap prefactor 11.
  3. [Temperature-transition discussion] Temperature-transition discussion: the claim that T_c∝n_d survives at finite density assumes the emergent Ising couplings remain ferromagnetic with γ strictly inside (2,3) after T-matrix resummation; the text does not address possible density-dependent corrections that would push γ outside this window or alter the linear gap scaling.
minor comments (2)
  1. [Abstract] The abstract states that 'certain excitations are also chiral' without identifying the excitations or directing the reader to the relevant figure or section.
  2. [Introduction] The contrast between the gapless-phase results and the gapped-phase findings of the companion preprint (arXiv:2511.19409) is mentioned only briefly; a short dedicated paragraph would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each of the major comments below and indicate the revisions we will make to strengthen the presentation.

read point-by-point responses
  1. Referee: [T-matrix analysis] T-matrix analysis section: the coefficient 11 in the reported topological gap Δ=11 n_d is stated as an output of the T-matrix calculation, yet no explicit T-matrix equations, matrix elements, or intermediate steps are shown; without these it is impossible to verify that the prefactor is independent of the chosen defect configurations or finite-size cutoffs.

    Authors: We agree that the T-matrix section would benefit from greater transparency. In the revised manuscript, we will include the explicit form of the T-matrix equations used, the matrix elements for the scattering off the Stone-Wales defects, and the intermediate steps in the calculation of the gap. The prefactor 11 is obtained from the low-energy effective theory projecting onto the Dirac fermions and is robust against variations in defect configurations within the dilute limit, as confirmed by our numerical validations. revision: yes

  2. Referee: [Finite-density numerics] Finite-density numerics section: the interaction exponent range 2<γ<3 is extracted from numerics, but the manuscript provides neither error bars, finite-size extrapolation of the fitted γ, nor explicit checks for higher-order multi-defect scattering or screening of the Majorana wavefunctions that could renormalize either γ or the gap prefactor 11.

    Authors: We acknowledge these omissions. We will add error bars to the extracted γ values from the numerical fits, include finite-size scaling analysis to extrapolate γ, and provide additional checks comparing results at different densities to assess multi-defect scattering effects. Regarding screening of Majorana wavefunctions, our approach incorporates the T-matrix resummation which accounts for multiple scatterings; we will clarify this point and discuss any potential renormalization effects. revision: yes

  3. Referee: [Temperature-transition discussion] Temperature-transition discussion: the claim that T_c∝n_d survives at finite density assumes the emergent Ising couplings remain ferromagnetic with γ strictly inside (2,3) after T-matrix resummation; the text does not address possible density-dependent corrections that would push γ outside this window or alter the linear gap scaling.

    Authors: The T-matrix analysis is performed at finite but low densities, and our results indicate that the ferromagnetic character and the exponent range are preserved. To address this concern, we will expand the discussion to include an analysis of density-dependent corrections, showing that within the regime studied, γ remains between 2 and 3 and the gap scaling stays linear. We will also note the conditions under which these assumptions hold. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper's quantitative claims—the topological gap coefficient of 11 n_d and the interaction exponent range 2<γ<3—are explicitly attributed to T-matrix analysis plus numerics performed at finite defect density. These constitute independent computational outputs on the defect model rather than quantities defined in terms of themselves or obtained by fitting a parameter and relabeling it a prediction. The self-citation to arXiv:2511.19409 supplies background on the defect construction but does not carry the load of the new gap or exponent results; no equation or step in the provided abstract reduces the reported outputs to the inputs by construction. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the exact solvability of the defected Kitaev Hamiltonian, the validity of the T-matrix approximation for dilute defects, and the emergence of effective Ising interactions from the chiral flux configurations. No new particles are postulated, but the coefficient 11 and the exponent bounds are numerical outputs whose independence from fitting is not shown in the abstract.

free parameters (2)
  • gap coefficient 11
    Numerical prefactor in the topological gap 11 n_d extracted from T-matrix analysis; its value is not derived from first principles in the abstract.
  • interaction exponent γ
    Power-law decay exponent reported in the range 2<γ<3 from numerics; bounds appear to be outputs rather than inputs but are not shown to be parameter-free.
axioms (2)
  • domain assumption The Stone-Wales defect preserves Kitaev bond labels and remains exactly solvable
    Invoked to enable exact treatment of isolated defects and flux configurations.
  • domain assumption T-matrix analysis accurately captures the effective interactions between dilute chiral defects
    Used to obtain the gap and the long-range Ising coupling without higher-order corrections.

pith-pipeline@v0.9.0 · 5617 in / 1688 out tokens · 40507 ms · 2026-05-16T22:10:53.197490+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Impurity quadrupole moments as local probes of flux sectors in the Kitaev spin liquid

    cond-mat.str-el 2026-03 unverdicted novelty 6.0

    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

66 extracted references · 66 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    and antialigned (µ z 1 =−µ z

  2. [2]

    Such a comparison is shown in Fig

    flux configura- tions. Such a comparison is shown in Fig. 8(a) for a small 48×48 system withN= 2304 sites. We consider defects separated by a vector⃗ rwhich is||x(horizontal separation ‘h’) or||y(vertical separation ‘v’), which show a distinct difference. In all cases we find that the aligned configu- ration has lower energy, consistent with a ferromagnet...

  3. [3]

    and antialigned (µz 1 =−µ z

  4. [4]

    bipartite

    chirality configurations of two SW defects in a 48×48 site system with open (OBC) (blue) or periodic (PBC) (black) boundary conditions. The energy difference is plotted after dividing by 2, asj≡∆E/2, to suggest an interpretation in terms of effective Ising spin interaction of Eq. 1. This numerically computed energy dif- ferencej(r) vanishes with increasin...

  5. [5]

    Imaginary hopping Inversion (P): ˜crA → −˜c−r,B ,˜c rB →˜c−r,A (B1) ˜ckA → −˜c−k,B ,˜c kB →˜c−k,A (B2) σx,z → −σ x,z , σ y →σ y (B3) τ x →τ x , τ y,z → −τ y,z (B4) Time-reversal (T): ˜crA →˜cr,A ,˜c rB → −˜cr,B (B5) ˜ckA →˜c−k,A ,˜c kB → −˜c−k,B (B6) σx → −σ x , σ y,z →σ y,z (B7) τ x,y →τ x,y , τ z → −τ z (B8)

  6. [6]

    The T-matrix for a defect potentialVis given by, T(E) =V(1−G 0(E)V) −1 (C1) whereG 0(E) represent the Green’s function of the unper- turbed Hamiltonian at energyE

    Real hopping Inversion (P): ¯crA →¯c−r,B ,¯c rB →¯c−r,A (B9) ¯ckA →¯c−k,B ,¯c kB →¯c−k,A (B10) σx →σ x , σ y,z → −σ y,z (B11) τ x →τ x , τ y,z → −τ y,z (B12) Time-reversal (T): ¯crA →¯cr,A ,¯c rB →¯cr,B (B13) ¯ckA →¯c−k,A ,¯c kB →¯c−k,B (B14) σx,z →σ x,z , σ y → −σ y (B15) τ x,y →τ x,y , τ z → −τ z (B16) Appendix C: Details of T-matrix computations In thi...

  7. [7]

    (C5) c1 =−9 √ 3π(t2 1 +t 2 2 −2t 1), d 1 = 18πt1 −(9 √ 3 + 6π)(t2 1 +t 2

  8. [8]

    (C6) PT-flux: ¯V SW PT = ¯c† R,A,¯c† R,B t1σ0 +it 2µz P T σz ¯cR+d+,A ¯cR−d+,B + h.c. (C7) P[ ¯V SW PT ] = 1 Nc ¯ψ† q − √ 3µz PTt2τ z −2t 1(Re , Im)e iπ/3 ·(σ x, σyτ z)−t 1(Re , Im)e −i(K−K ′)·R ·(τ x, τ y)σx ¯ψq (C8) P[ ¯T SW PT ] = 1 Nc ¯ψ† qf2 a2t2µz PTτ zσz +b 2σx +c 2σyτ z +d 2 (Re,Im)e −i(K−K ′)·R ·(τ x, τ y)σ x ¯ψq (C9) f2 = 1 (3 √ 3−4π)(t 2 2 −t 2

  9. [9]

    + 18π−12πt 1 , a 2 = 18 √ 3π , b 2 =−18πt 1 + (3π−9 √ 3)(t2 2 −t 2

  10. [10]

    (C10) c2 =−9 √ 3π(t2 2 −t 2 1 + 2t1), d 2 =−18πt 1 −(9 √ 3 + 6π)(t2 2 −t 2

  11. [11]

    (C11) 19 Uniform-flux: ¯V Uniform PT = ¯c† R,A,¯c† R,B t1σ0 +it 2µz P T σz ¯cR+d+,A ¯cR−d+,B + 2t1¯c† R,B¯cR+d−,A + h.c. (C12) P[ ¯V SW Uniform] = 1 Nc ¯ψ† q − √ 3µz Ut2τ z −2t 1(Re , Im)e iπ/3 ·(σ x, σyτ z)−t ′ 1(Re , Im)e −iπ/3 ·(σ x, σyτ z) −t1(Re , Im)e i(K−K ′)·R ·(τ x, τ y)σx −t ′ 1(Re , Im) e−i(K−K ′)·Reiπ/3 ·(τ x, τ y)σx ¯ψq (C13) P[ ¯T SW Uniform...

  12. [12]

    Kimchi, A

    I. Kimchi, A. Nahum, and T. Senthil, Valence bonds in random quantum magnets: Theory and application to YbMgGaO4, Physical Review X8, 031028 (2018)

  13. [13]

    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)

  14. [14]

    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, Physical Review Letters102, 017205 (2009)

  15. [15]

    Hermanns, I

    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)

  16. [16]

    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)

  17. [17]

    Trebst and C

    S. Trebst and C. Hickey, Kitaev materials, Physics Re- ports Kitaev Materials,950, 1 (2022)

  18. [18]

    Savary and L

    L. Savary and L. Balents, Quantum spin liquids: A re- view, Reports on Progress in Physics80, 016502 (2016)

  19. [19]

    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)

  20. [20]

    A. J. Willans, J. T. Chalker, and R. Moessner, Site dilu- tion in the Kitaev honeycomb model, Physical Review B 84, 115146 (2011). 21

  21. [21]

    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)

  22. [22]

    Nasu and Y

    J. Nasu and Y. Motome, Thermodynamic and transport properties in disordered Kitaev models, Physical Review B102, 054437 (2020)

  23. [23]

    W.-H. Kao, N. B. Perkins, and G. B. Hal´ asz, Va- cancy spectroscopy of non-Abelian Kitaev spin liquids, arXiv:2307.10376 (2023), arXiv:2307.10376 [cond-mat]

  24. [24]

    Dhochak, R

    K. Dhochak, R. Shankar, and V. Tripathi, Magnetic Im- purities in the Honeycomb Kitaev Model, Physical Re- view Letters105, 117201 (2010)

  25. [25]

    Kao and N

    W.-H. Kao and N. B. Perkins, Disorder upon disorder: Localization effects in the Kitaev spin liquid, Annals of Physics435, 168506 (2021)

  26. [26]

    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)

  27. [27]

    Singhania, J

    A. Singhania, J. Van Den Brink, and S. Nishimoto, Dis- order effects in the Kitaev-Heisenberg model, Physical Review Research5, 023009 (2023)

  28. [28]

    Zschocke and M

    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)

  29. [29]

    Knolle, R

    J. Knolle, R. Moessner, and N. B. Perkins, Bond- disordered spin liquid and the honeycomb iridate H3LiIr2O6: Abundant low-energy density of states from random majorana hopping, Physical Review Letters122, 047202 (2019)

  30. [30]

    Cassella, P

    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]

  31. [31]

    A. G. Grushin and C. Repellin, Amorphous and poly- crystalline routes toward a chiral spin liquid, Physical Review Letters130, 186702 (2023)

  32. [32]

    Petrova, P

    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)

  33. [33]

    Sanyal, K

    S. Sanyal, K. Damle, J. T. Chalker, and R. Moessner, Emergent Moments and Random Singlet Physics in a Majorana Spin Liquid, Physical Review Letters127, 127201 (2021)

  34. [34]

    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)

  35. [35]

    G. B. Hal´ asz and J. T. Chalker, Coherent hole propaga- tion in an exactly solvable gapless spin liquid, Physical Review B94, 235105 (2016)

  36. [36]

    Lahtinen, A

    V. Lahtinen, A. W. W. Ludwig, and S. Trebst, Perturbed vortex lattices and the stability of nucleated topological phases, Physical Review B89, 085121 (2014)

  37. [37]

    Udagawa, Vison-Majorana complex zero-energy res- onance in the Kitaev spin liquid, Physical Review B98, 220404 (2018)

    M. Udagawa, Vison-Majorana complex zero-energy res- onance in the Kitaev spin liquid, Physical Review B98, 220404 (2018)

  38. [38]

    Kimchi, J

    I. Kimchi, J. P. Sheckelton, T. M. McQueen, and P. A. Lee, Scaling and data collapse from local moments in frustrated disordered quantum spin systems, Nature Communications9, 4367 (2018)

  39. [39]

    Otten, A

    D. Otten, A. Roy, and F. Hassler, Dynamical structure factor in the non-Abelian phase of the Kitaev honeycomb model in the presence of quenched disorder, Physical Re- view B99, 035137 (2019)

  40. [40]

    Nasu and Y

    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)

  41. [41]

    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)

  42. [42]

    Dantas and E

    V. Dantas and E. C. Andrade, Disorder, Low-Energy Ex- citations, and Topology in the Kitaev Spin Liquid, Phys- ical Review Letters129, 037204 (2022)

  43. [43]

    Song, Y.-Z

    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)

  44. [44]

    Yatsuta and D

    I. Yatsuta and D. F. Mross, Vacancies in generic Kitaev spin liquids, arXiv:2312.00147 (2023), arXiv:2312.00147 [cond-mat]

  45. [45]

    Vojta, A

    M. Vojta, A. K. Mitchell, and F. Zschocke, Kondo Impu- rities in the Kitaev Spin Liquid: Numerical Renormaliza- tion Group Solution and Gauge-Flux-Driven Screening, Physical Review Letters117, 037202 (2016)

  46. [46]

    A. Seth, F. Borhani, and I. Kimchi, Chiral spin liquid instability of the Kitaev honeycomb model with crystal- lographic defects (2025), arXiv:2511.19409 [cond-mat]

  47. [47]

    A. J. Stone and D. J. Wales, Theoretical studies of icosa- hedral C60 and some related species, Chemical Physics Letters128, 501 (1986)

  48. [48]

    Vozmediano, M

    M. Vozmediano, M. Katsnelson, and F. Guinea, Gauge fields in graphene, Physics Reports496, 109 (2010)

  49. [49]

    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)

  50. [50]

    Borhani, A

    F. Borhani, A. Seth, and I. Kimchi, Real-space chirality from crystalline topological defects in the Kitaev spin liquid, npj Quantum Materials10, 1 (2025)

  51. [51]

    Yao and S

    H. Yao and S. A. Kivelson, Exact chiral spin liquid with non-abelian anyons, Physical Review Letters99, 247203 (2007)

  52. [52]

    d’Ornellas and J

    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)

  53. [53]

    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)

  54. [54]

    M. G. Yamada, H. Fujita, and M. Oshikawa, Designing Kitaev Spin Liquids in Metal-Organic Frameworks, Phys- ical Review Letters119, 057202 (2017)

  55. [55]

    J. Ma, D. Alf` e, A. Michaelides, and E. Wang, Stone- Wales defects in graphene and other planar sp 2-bonded materials, Physical Review B80, 033407 (2009)

  56. [56]

    E. H. Lieb, Flux Phase of the Half-Filled Band, Physical Review Letters73, 2158 (1994)

  57. [57]

    Panigrahi, P

    A. Panigrahi, P. Coleman, and A. Tsvelik, Analytic cal- culation of the vison gap in the Kitaev spin liquid, Phys- ical Review B108, 045151 (2023)

  58. [58]

    R. Xu, A. Seth, and I. Kimchi, Chirality reversal at fi- nite magnetic impurity strength and local signatures of a topological phase transition (2025), arXiv:2510.11707 22 [cond-mat]

  59. [59]

    Shindou and N

    R. Shindou and N. Nagaosa, Orbital Ferromagnetism and Anomalous Hall Effect in Antiferromagnets on the Dis- torted fcc Lattice, Physical Review Letters87, 116801 (2001)

  60. [60]

    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)

  61. [61]

    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)

  62. [62]

    Bianco and R

    R. Bianco and R. Resta, Mapping topological order in coordinate space, Physical Review B84, 241106 (2011)

  63. [63]

    d’Ornellas, R

    P. d’Ornellas, R. Barnett, and D. K. K. Lee, Quantized bulk conductivity as a local Chern marker, Physical Re- view B106, 155124 (2022)

  64. [64]

    Kogan, RKKY Interaction in Gapped or Doped Graphene, Graphene2, 8 (2013)

    E. Kogan, RKKY Interaction in Gapped or Doped Graphene, Graphene2, 8 (2013)

  65. [65]

    J. Nasu, M. Udagawa, and Y. Motome, Thermal frac- tionalization of quantum spins in a Kitaev model: Temperature-linear specific heat and coherent transport of Majorana fermions, Physical Review B92, 115122 (2015)

  66. [66]

    Rousochatzakis, S

    I. Rousochatzakis, S. Kourtis, J. Knolle, R. Moessner, and N. B. Perkins, Quantum spin liquid at finite tem- perature: Proximate dynamics and persistent typicality, Physical Review B100, 045117 (2019)