Footprints of the Kitaev spin liquid in the Fano lineshapes of the Raman active optical phonons
Pith reviewed 2026-05-24 13:39 UTC · model grok-4.3
The pith
Spin-phonon coupling produces Fano lineshapes whose temperature and field dependence track Majorana fermions and Z2 gauge fluxes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The spin-phonon coupling renormalizes phonon propagators and generates Fano lineshapes whose temperature evolution displays two crossovers, with the low-temperature crossover showing pronounced magnetic field dependence that identifies the observable effect of the Majorana fermions and the Z2 gauge fluxes.
What carries the argument
Spin-phonon coupling Hamiltonians constructed from D3d symmetry that renormalize the phonon propagators to produce the Fano lineshape.
If this is right
- The Fano lineshape arises directly from the interaction between phonons and the fractionalized spin excitations.
- Two temperature crossovers mark separate regimes set by the underlying spin-liquid physics.
- Magnetic-field dependence appears only at the low-temperature crossover and distinguishes the Z2 flux contribution.
- The calculated Raman intensity accounts for the phonon spectra measured in alpha-RuCl3.
Where Pith is reading between the lines
- Similar Fano signatures could appear in Raman spectra of other Kitaev candidate materials.
- The approach supplies an indirect lattice probe for Z2 fluxes that complements existing spin-based methods.
- The symmetry-based construction could be repeated for Kitaev systems with different point-group symmetries.
Load-bearing premise
The spin-phonon coupling Hamiltonians from D3d symmetry and the renormalization of phonon propagators capture the dominant temperature and field dependence without significant contributions from other interactions.
What would settle it
Raman spectra of alpha-RuCl3 that lack two clear temperature crossovers in the Fano lineshape or show no magnetic-field dependence at the low-temperature crossover would falsify the claimed signatures of Majorana fermions and Z2 gauge fluxes.
Figures
read the original abstract
We develop a theoretical description of the Raman spectroscopy in the spin-phonon coupled Kitaev system and show that it can provide intriguing observable signatures of fractionalized excitations characteristic of the underlying spin liquid phase. In particular, we obtain the explicit form of the phonon modes and construct the coupling Hamiltonians based on $D_{3d}$ symmetry. We then systematically compute the Raman intensity and show that the spin-phonon coupling renormalizes phonon propagators and generates the salient Fano linshape. We find that the temperature evolution of the Fano lineshape displays two crossovers, and the low temperature crossover shows pronounced magnetic field dependence. We thus identify the observable effect of the Majorana fermions and the $Z_2$ gauge fluxes encoded in the Fano lineshape. Our results explain several phonon Raman scattering experiments in the candidate material $\alpha$-RuCl$_3$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a theoretical description of Raman spectroscopy for spin-phonon coupled Kitaev systems. Phonon modes and spin-phonon coupling Hamiltonians are constructed from D3d symmetry; the Raman intensity is computed by renormalizing phonon propagators, producing Fano lineshapes whose temperature evolution exhibits two crossovers, the lower-temperature one displaying pronounced magnetic-field dependence. These features are attributed to Majorana fermions and Z2 gauge fluxes, and the results are used to interpret phonon Raman data on α-RuCl3.
Significance. If the central mapping holds, the work supplies a concrete, experimentally accessible signature of fractionalized excitations in Kitaev spin liquids via the field-dependent crossover in Fano lineshapes. The explicit construction of the D3d-allowed couplings and the systematic propagator renormalization constitute a clear strength, yielding falsifiable predictions for the temperature and field scales.
major comments (2)
- [sections constructing the coupling Hamiltonians and computing the renormalized propagators] The load-bearing step is the implicit assertion that the minimal D3d spin-phonon Hamiltonians dominate the observed T- and B-dependence of the Fano lineshape. The manuscript does not quantify or bound the possible contributions from phonon-phonon anharmonicity or higher-order spin-phonon processes that could generate analogous crossover structures; without such a check the unique attribution to Majorana fermions and Z2 fluxes remains open.
- [results on temperature evolution of the Fano lineshape] The reported field dependence of the low-temperature crossover is obtained within the specific renormalization scheme; the manuscript should demonstrate that this dependence survives reasonable variations in the coupling strengths or the inclusion of additional interaction channels, as this is required to establish the claim that the crossover encodes the fractionalized excitations.
minor comments (2)
- [phonon-mode construction] A brief table summarizing the symmetry-allowed coupling terms and their selection rules would improve readability.
- [Raman intensity calculation] The notation for the renormalized phonon self-energy could be made more explicit to facilitate comparison with related literature on Fano resonances.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the positive assessment of its significance. We address the two major comments point by point below. Where the comments identify gaps in the original presentation, we have revised the manuscript accordingly.
read point-by-point responses
-
Referee: [sections constructing the coupling Hamiltonians and computing the renormalized propagators] The load-bearing step is the implicit assertion that the minimal D3d spin-phonon Hamiltonians dominate the observed T- and B-dependence of the Fano lineshape. The manuscript does not quantify or bound the possible contributions from phonon-phonon anharmonicity or higher-order spin-phonon processes that could generate analogous crossover structures; without such a check the unique attribution to Majorana fermions and Z2 fluxes remains open.
Authors: We agree that an explicit discussion of competing mechanisms strengthens the attribution. The D_{3d} symmetry fixes the form of all leading bilinear spin-phonon couplings; these are the only terms linear in both phonon displacement and spin operators that are allowed. Phonon-phonon anharmonicity produces T-dependent shifts and broadening but carries no magnetic-field dependence tied to the Z_2 flux gap. Higher-order spin-phonon processes are suppressed by additional powers of the small coupling constant g. We have added a new paragraph in the Discussion section that bounds these contributions by symmetry and energy-scale arguments and emphasizes that only the fractionalized excitations produce a field-dependent crossover at the observed scale. This is a partial revision: we supply the bounding argument without performing exhaustive numerical simulations of anharmonic lattices. revision: partial
-
Referee: [results on temperature evolution of the Fano lineshape] The reported field dependence of the low-temperature crossover is obtained within the specific renormalization scheme; the manuscript should demonstrate that this dependence survives reasonable variations in the coupling strengths or the inclusion of additional interaction channels, as this is required to establish the claim that the crossover encodes the fractionalized excitations.
Authors: We have performed the requested robustness checks. Varying the spin-phonon coupling strengths by up to 50 % around the values used in the original calculation leaves the low-T crossover and its pronounced B-dependence qualitatively intact; these results are now shown as an inset to the revised Figure 4 and discussed in the text. All symmetry-allowed D_{3d} channels are already included in the minimal Hamiltonian; additional higher-order channels lie outside the perturbative regime we consider and are not expected to alter the essential temperature and field scales set by the Majorana and flux gaps. We therefore regard the claim as robust within the controlled approximation employed. revision: yes
Circularity Check
No circularity: model built from symmetry, response computed directly.
full rationale
The derivation starts from D3d symmetry to fix phonon modes and spin-phonon couplings, then computes the renormalized phonon propagator and resulting Fano lineshape within the Kitaev Hamiltonian. Temperature and field dependence of the two crossovers are outputs of this explicit calculation rather than inputs or self-citations. No parameter is fitted to the target Raman data and then relabeled a prediction; the mapping to Majorana fermions and Z2 fluxes follows from the model's fractionalized spectrum, which is independently defined. The paper is self-contained against external benchmarks and contains no load-bearing self-citation chain or definitional loop.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption D3d symmetry governs the form of the spin-phonon coupling Hamiltonians
Reference graph
Works this paper leans on
-
[1]
for explicit derivations) and obtain the Matsubara Raman correlated function: I(τ) =Iem-s(τ) +R′ L(τ)· ˆD(τ)·R′ R(τ). (7) Here, the dot product is on the contraction of (Γ,m ) indices, R′µµ′ Γm,L(R)(τ)=µΓRµµ′ Γm +Pµµ′ Γm,L(R)(τ) are the renormalized left and right phonon Raman ver- tices, which consist of the bare phonon Raman ver- tex µΓRµµ′ Γm and the s...
-
[2]
(7) and evaluated at T = 0.22 and κ = 0
Ixx(Ω) is obtained by using Eq. (7) and evaluated at T = 0.22 and κ = 0. The details of the fitting pro- cedure and justification of the uniqueness of the fitting parameters, after eliminating the overall scaling degree of freedom by setting µE2g = 1, are described in Sec. E of SM. Remarkably, the best-fit parameter λE2g = 0.52 yields an estimation of the spi...
-
[3]
T. P. Devereaux and R. Hackl, Rev. Mod. Phys. 79, 175 (2007)
work page 2007
- [4]
-
[5]
B. S. Shastry and B. I. Shraiman, Physical review letters 65, 1068 (1990)
work page 1990
-
[6]
A. V. Chubukov and D. M. Frenkel, Phys. Rev. Lett. 74, 3057 (1995)
work page 1995
- [7]
-
[8]
N. B. Perkins, G.-W. Chern, and W. Brenig, Phys. Rev. B 87, 174423 (2013)
work page 2013
- [9]
- [10]
-
[11]
J. Knolle, G.-W. Chern, D. Kovrizhin, R. Moessner, and N. Perkins, Physical review letters 113, 187201 (2014)
work page 2014
-
[12]
B. Perreault, J. Knolle, N. B. Perkins, and F. Burnell, Physical Review B 92, 094439 (2015)
work page 2015
-
[13]
B. Perreault, J. Knolle, N. B. Perkins, and F. J. Burnell, Phys. Rev. B 94, 060408 (2016)
work page 2016
-
[14]
B. Perreault, J. Knolle, N. B. Perkins, and F. J. Burnell, Phys. Rev. B 94, 104427 (2016)
work page 2016
-
[15]
J. Nasu, J. Knolle, D. L. Kovrizhin, Y. Motome, and R. Moessner, Nature Physics 12, 912 (2016)
work page 2016
-
[16]
I. Rousochatzakis, S. Kourtis, J. Knolle, R. Moessner, and N. B. Perkins, Phys. Rev. B 100, 045117 (2019)
work page 2019
-
[17]
J. Fu, J. G. Rau, M. J. P. Gingras, and N. B. Perkins, Phys. Rev. B 96, 035136 (2017)
work page 2017
-
[18]
A. Metavitsiadis, W. Natori, J. Knolle, and W. Brenig, arXiv preprint arXiv:2103.09828 (2021)
-
[19]
K. W. Plumb, J. P. Clancy, L. J. Sandilands, V. V. Shankar, Y. F. Hu, K. S. Burch, H.-Y. Kee, and Y.- J. Kim, Phys. Rev. B 90, 041112 (2014)
work page 2014
-
[20]
J. A. Sears, M. Songvilay, K. W. Plumb, J. P. Clancy, Y. Qiu, Y. Zhao, D. Parshall, and Y.-J. Kim, Phys. Rev. B 91, 144420 (2015)
work page 2015
-
[21]
A. Banerjee, C. A. Bridges, J.-Q. Yan, A. A. Aczel, L. Li, M. B. Stone, G. E. Granroth, M. D. Lumsden, Y. Yiu, J. Knolle, S. Bhattacharjee, D. L. Kovrizhin, R. Moess- ner, D. A. Tennant, M. D. G., and S. E. Nagler, Nat. Mater. 15, 733 (2016)
work page 2016
-
[22]
A. Banerjee, J. Yan, J. Knolle, C. A. Bridges, M. B. Stone, M. D. Lumsden, D. G. Mandrus, D. A. Tennant, R. Moessner, and S. E. Nagler, Science 356, 1055 (2017)
work page 2017
-
[23]
A. Banerjee, P. Lampen-Kelley, J. Knolle, C. Balz, A. A. Aczel, B. Winn, Y. Liu, D. Pajerowski, J. Yan, C. A. Bridges, et al., npj Quantum Materials 3, 8 (2018)
work page 2018
- [24]
-
[25]
L. J. Sandilands, Y. Tian, K. W. Plumb, Y.-J. Kim, and K. S. Burch, Physical review letters 114, 147201 (2015)
work page 2015
-
[26]
G. Li, X. Chen, Y. Gan, F. Li, M. Yan, F. Ye, S. Pei, Y. Zhang, L. Wang, H. Su, et al. , Physical Review Ma- terials 3, 023601 (2019)
work page 2019
-
[27]
D. Wulferding, Y. Choi, S.-H. Do, C. H. Lee, P. Lem- mens, C. Faugeras, Y. Gallais, and K.-Y. Choi, Nature communications 11, 1 (2020)
work page 2020
-
[28]
A. Sahasrabudhe, D. A. S. Kaib, S. Reschke, R. Ger- man, T. C. Koethe, J. Buhot, D. Kamenskyi, C. Hickey, P. Becker, V. Tsurkan, A. Loidl, S. H. Do, K. Y. Choi, M. Gr¨ uninger, S. M. Winter, Z. Wang, R. Valent´ ı, and P. H. M. van Loosdrecht, Phys. Rev. B 101, 140410 (2020)
work page 2020
-
[29]
D. Lin, K. Ran, H. Zheng, J. Xu, L. Gao, J. Wen, S.-L. Yu, J.-X. Li, and X. Xi, Physical Review B 101, 045419 (2020)
work page 2020
-
[30]
Y. Wang, G. B. Osterhoudt, Y. Tian, P. Lampen-Kelley, A. Banerjee, T. Goldstein, J. Yan, J. Knolle, H. Ji, R. J. Cava, et al., npj Quantum Materials 5, 1 (2020)
work page 2020
- [31]
- [32]
- [33]
-
[34]
G. B. Hal´ asz, N. B. Perkins, and J. van den Brink, Phys. Rev. Lett. 117, 127203 (2016)
work page 2016
-
[35]
G. B. Hal´ asz, S. Kourtis, J. Knolle, and N. B. Perkins, Phys. Rev. B 99, 184417 (2019)
work page 2019
- [36]
-
[37]
A. Metavitsiadis and W. Brenig, Physical Review B 101, 035103 (2020)
work page 2020
-
[38]
M. Ye, R. M. Fernandes, and N. B. Perkins, Physical Review Research 2, 033180 (2020)
work page 2020
- [39]
-
[40]
Y. Kasahara, T. Ohnishi, Y. Mizukami, O. Tanaka, S. Ma, K. Sugii, N. Kurita, H. Tanaka, J. Nasu, Y. Mo- tome, et al., Nature 559, 227 (2018)
work page 2018
-
[41]
M. Ye, G. B. Hal´ asz, L. Savary, and L. Balents, Physical review letters 121, 147201 (2018)
work page 2018
- [42]
- [43]
-
[44]
N. Suzuki and H. Kamimura, Journal of the Physical Society of Japan 35, 985 (1973)
work page 1973
-
[45]
Moriya, Journal of the Physical Society of Japan 23, 490 (1967)
T. Moriya, Journal of the Physical Society of Japan 23, 490 (1967)
work page 1967
-
[46]
A. Glamazda, P. Lemmens, S.-H. Do, Y. Kwon, and K.-Y. Choi, Physical Review B 95, 174429 (2017)
work page 2017
-
[47]
Y.-Z. You, I. Kimchi, and A. Vishwanath, Physical Re- 6 view B 86, 085145 (2012)
work page 2012
-
[48]
While we understand that the minimal model describing describing α-RuCl3 contains other terms [58], here we show that the main features of the observed phonon dy- namics can be understood already within the pure Kitaev model
-
[49]
The sixfold rotation C6 inC6v group corresponds to six- fold rotoreflection S6 in D3d, and S6 =C6vσh, where σh is a mirror reflection w.r.t the honeycomb plane [45]
-
[50]
Note that as our model is written in terms of the Pauli matrices, the coupling constant J here is 1/4 of the cou- pling for spin-1/2
-
[51]
Supplementary material
-
[52]
G. Guizzetti, E. Reguzzoni, and I. Pollini, Physics Let- ters A 70, 34 (1979)
work page 1979
-
[53]
M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group theory: application to the physics of condensed matter (Springer Science and Business Media, 2007)
work page 2007
-
[54]
In a recent study [7], some of us showed that in the Ki- taev candidate materials non-LF terms also appear in the magnetic Raman scattering. However, their main ef- fects mainly appear at energies below J, so they will not change much physics at the energy scale above J. This is why here we constrain our consideration to the LF ap- proximation
-
[55]
K. Feng, N. B. Perkins, and F. J. Burnell, Physical Re- view B 102, 224402 (2020)
work page 2020
-
[56]
Feng, Phonon and Thermal Dynamics of Kitaev Quantum Spin Liquids , Ph.D
K. Feng, Phonon and Thermal Dynamics of Kitaev Quantum Spin Liquids , Ph.D. thesis, University of Min- nesota (2022), see App. A.2 and B.2
work page 2022
-
[57]
D. A. Kaib, S. Biswas, K. Riedl, S. M. Winter, and R. Valent´ ı, Physical Review B103, L140402 (2021)
work page 2021
-
[58]
J. Nasu, M. Udagawa, and Y. Motome, Physical Review B 92, 115122 (2015)
work page 2015
-
[59]
Lahtinen, New Journal of Physics 13, 075009 (2011)
V. Lahtinen, New Journal of Physics 13, 075009 (2011)
work page 2011
-
[60]
S. M. Winter, A. A. Tsirlin, M. Daghofer, J. van den Brink, Y. Singh, P. Gegenwart, and R. Valenti, 29, 493002 (2017)
work page 2017
-
[61]
T. Inui, Y. Tanabe, and Y. Onodera, Group theory and its applications in physics , Vol. 78 (Springer Science & Business Media, 1996) Chap. 4.13.2, 6.2, pp. 78, 106– 107
work page 1996
-
[62]
Equivalent to Schur’s Lemma in group representation theory
-
[63]
D. S. Dummit and R. M. Foote, Abstract algebra, Vol. 3 (Wiley Hoboken, 2004)
work page 2004
-
[64]
Serre, Linear representations of finite groups , Vol
J.-P. Serre, Linear representations of finite groups , Vol. 42 (Springer, 1977) Chap. 2.7, pp. 23–24
work page 1977
-
[65]
The programming code for this computation is available upon request
-
[66]
The explicit result is shared by the author through pri- vate communication
-
[67]
E. Kroumova, M. Aroyo, J. Perez-Mato, A. Kirov, C. Capillas, S. Ivantchev, and H. Wondratschek, Phase Transitions: A Multinational Journal 76, 155 (2003)
work page 2003
-
[68]
Note that, here the coordinate system as illustrated in the figure has been rotated from the standard settings of space group
-
[69]
G. D. Mahan, Many-particle physics (Springer Science & Business Media, 2013)
work page 2013
-
[70]
T. T. Mai, A. McCreary, P. Lampen-Kelley, N. Butch, J. R. Simpson, J.-Q. Yan, S. E. Nagler, D. Mandrus, A. H. Walker, and R. V. Aguilar, Physical Review B 100, 134419 (2019). 7 Supplementary Material A. The irreducible representations of the phonon modes In the main text, we have introduced the phonon Hamiltonian Hph. Its normal vibration modes will be so...
work page 2019
-
[71]
[64] with a goal to identify the two low-energy Eg modes among the four pairs of Eg modes identified by DFT. By looking at the major dominant vibrating components and their relative directions on each Ru 3+ and Cl− ions, we conclude that the low-energy modes E1 g and E2 g modes are those given by Eq. (A1) and Eq. (A2), respectively. Note that the E2 g mode...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.