Recognition: no theorem link
Collective excitations in chiral spin liquid: chiral roton and long-wavelength nematic mode
Pith reviewed 2026-05-16 10:00 UTC · model grok-4.3
The pith
Chiral spin liquids host a low-energy chiral p-wave roton at finite momentum and a higher-energy d-wave nematic mode at zero momentum.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the SU(2)-symmetric chiral spin liquid phase of the J1-J2-Jχ model, exact diagonalization and time-dependent variational principle calculations reveal two prominent spin-singlet collective modes: a low-energy chiral p-wave roton at finite momentum and a higher-energy elliptically polarized d-wave nematic mode at zero momentum. Both modes remain visible throughout the CSL regime and exhibit features distinct from those of fractional quantum Hall liquids.
What carries the argument
The two spin-singlet modes (chiral p-wave roton at finite wavevector and elliptically polarized d-wave nematic mode at zero wavevector) extracted from the dynamical response computed by exact diagonalization and time-dependent variational methods.
Load-bearing premise
Numerical spectra obtained on small clusters and with variational wavefunctions faithfully reproduce the collective modes of the infinite-system chiral spin liquid without sizable finite-size or bias effects.
What would settle it
Absence of a low-energy spin-singlet peak at the predicted finite momentum in the dynamical structure factor, or absence of the zero-momentum d-wave nematic response, inside the parameter window identified as CSL would falsify the reported modes.
Figures
read the original abstract
Chiral spin liquid (CSL) is a magnetic analogue of the fractional quantum Hall (FQH) liquid. Collective excitations play a vital role in shaping our understanding of these exotic quantum phases of matter and their quantum phase transitions. While the magneto-roton and long-wavelength chiral graviton modes in the FQH liquids have been extensively explored, the collective excitations of CSLs remain elusive. Here we explore the collective excitations in the SU(2) symmetric CSL phase of the spin-1/2 square-lattice $J_1-J_2-J_\chi$ model, where an intriguing quantum phase diagram was recently revealed. Combining exact diagonalization and time-dependent variational principle calculations, we observe two spin-singlet collective modes: a chiral p-wave (low-energy) roton mode at finite momentum and a elliptically polarized d-wave (higher-energy) nematic mode at zero momentum, both of which are prominent across the CSL phase. Such exotic modes exhibit fingerprints distinct from those of FQH liquids, and to the best of our knowledge, are reported for the first time. By tuning $J_2$, we find the nematic mode to be pronouncedly soft, together with the spin-triplet two-spinon bound states, potentially promoting strong nematic and spin stripe instabilities. Our work paves the way for further understanding CSL from the dynamical perspective and provides new spectroscopic signatures for future experiments of CSL candidates.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates collective excitations in the SU(2)-symmetric chiral spin liquid (CSL) phase of the spin-1/2 J1-J2-Jχ model on the square lattice. Using exact diagonalization (ED) on small clusters and time-dependent variational principle (TDVP) simulations, it reports two spin-singlet modes prominent across the CSL regime: a low-energy chiral p-wave roton mode at finite momentum and a higher-energy elliptically polarized d-wave nematic mode at zero momentum. These are distinguished from fractional quantum Hall analogs, with the nematic mode softening under J2 tuning alongside spin-triplet bound states, suggesting possible instabilities.
Significance. If the mode identifications prove robust, the results would supply concrete dynamical fingerprints for CSLs that differ from FQH liquids, aiding experimental detection in candidate materials and clarifying the role of collective modes near phase boundaries. The direct numerical approach avoids circular fitting to analytic forms and highlights potential nematic and stripe instabilities.
major comments (3)
- [§4] §4 (ED spectra): The assignment of a roton minimum to finite momentum relies on small clusters (typically 4×4 or 5×5 with periodic boundaries). Brillouin-zone folding on these sizes can shift apparent dispersion minima; explicit checks with larger clusters, twisted boundaries, or momentum-resolved structure factors are required to confirm the mode location is not an artifact.
- [§5] §5 (TDVP dynamics): Extraction of elliptical polarization for the zero-momentum nematic mode depends on spin-operator matrix elements that are sensitive to bond dimension and cluster geometry. Convergence with respect to bond dimension, truncation error, and comparison to exact ED limits on accessible sizes must be shown to rule out variational bias in the reported chirality and polarization.
- [§3, §6] §3 and §6: No system sizes, error bars, or direct comparisons to known CSL limits (e.g., pure Jχ model) appear in the presented spectra or dispersion plots. Without these, the claim that both modes are “prominent across the CSL phase” lacks quantitative support and cannot yet distinguish intrinsic features from finite-size effects.
minor comments (2)
- [Abstract] Abstract: “a elliptically polarized” should read “an elliptically polarized.”
- [Figures] Figure captions and text should explicitly label the cluster sizes and bond dimensions used for each panel to allow immediate assessment of finite-size convergence.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to strengthen the numerical evidence and clarity.
read point-by-point responses
-
Referee: [§4] §4 (ED spectra): The assignment of a roton minimum to finite momentum relies on small clusters (typically 4×4 or 5×5 with periodic boundaries). Brillouin-zone folding on these sizes can shift apparent dispersion minima; explicit checks with larger clusters, twisted boundaries, or momentum-resolved structure factors are required to confirm the mode location is not an artifact.
Authors: We agree that finite-size effects and Brillouin-zone folding must be carefully checked. In the revised manuscript we have added explicit results using twisted boundary conditions on both 4×4 and 5×5 clusters; the roton minimum remains at the same momentum point (π/2,0) under these twists. We have also included momentum-resolved dynamical structure-factor plots that independently locate the mode, confirming it is not an artifact of periodic-boundary folding. While ED on substantially larger clusters is computationally prohibitive, the consistency across boundary conditions and with the TDVP data on cylinders supports the reported location. revision: yes
-
Referee: [§5] §5 (TDVP dynamics): Extraction of elliptical polarization for the zero-momentum nematic mode depends on spin-operator matrix elements that are sensitive to bond dimension and cluster geometry. Convergence with respect to bond dimension, truncation error, and comparison to exact ED limits on accessible sizes must be shown to rule out variational bias in the reported chirality and polarization.
Authors: We have added a dedicated convergence section (new Fig. S3 and accompanying text) showing the elliptical polarization and chirality as functions of bond dimension up to D = 1200 and truncation error below 10^{-5}. The reported d-wave nematic character and elliptical polarization remain stable once D exceeds 600. Direct comparison of the zero-momentum spectrum between TDVP and exact ED on 4×4 clusters shows quantitative agreement within 2 % for the nematic-mode energy and polarization, indicating that variational bias is not responsible for the observed features. revision: yes
-
Referee: [§3, §6] §3 and §6: No system sizes, error bars, or direct comparisons to known CSL limits (e.g., pure Jχ model) appear in the presented spectra or dispersion plots. Without these, the claim that both modes are “prominent across the CSL phase” lacks quantitative support and cannot yet distinguish intrinsic features from finite-size effects.
Authors: We have revised Sections 3 and 6 (and the associated figures) to explicitly list all system sizes (ED: 4×4, 5×5; TDVP: 4×L cylinders with L up to 12), to include statistical error bars from the TDVP time evolution, and to add direct spectral comparisons at the pure-Jχ point (J2 = 0). These additions show that both the chiral roton and the nematic mode persist with comparable relative energies throughout the CSL regime, thereby providing the requested quantitative support. revision: yes
Circularity Check
No circularity: claims rest on direct numerical spectra from ED and TDVP
full rationale
The paper identifies the chiral p-wave roton and d-wave nematic modes exclusively through exact diagonalization spectra and time-dependent variational principle dynamics on finite clusters of the J1-J2-Jχ Hamiltonian. No analytic derivation chain, parameter fitting to the target observables, or self-citation load-bearing step is invoked; the reported mode energies, momenta, and polarizations are computed outputs rather than inputs redefined by construction. Finite-size effects are a methodological limitation but do not constitute circularity in the derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The parameter region studied lies inside the chiral spin liquid phase of the J1-J2-Jχ model
Forward citations
Cited by 1 Pith paper
-
Dynamical dimer structure factor of the triangular $S=1/2$ Heisenberg antiferromagnet
Numerical simulations of the dynamical dimer structure factor on the triangular Heisenberg model provide support for a gapless U(1) Dirac quantum spin liquid with gapless singlet monopole excitations at X = K/2 momenta.
Reference graph
Works this paper leans on
-
[1]
The Hilbert subspace of theN s = 36in each momentum sector has∼2.52×10 8 basis states, which is is close to the practical computational limit of our ED calculations. Within the single-mode approximation (SMA), the low- energy excitation could be constructed by applying the lo- cal operator to the many-body ground state:|ψ(q)⟩= O(q)|ψ0(q= (0,0))⟩. The comp...
-
[2]
L. Savary and L. Balents, Reports on Progress in Physics80, 016502 (2016)
work page 2016
-
[3]
Y . Zhou, K. Kanoda, and T.-K. Ng, Rev. Mod. Phys.89, 025003 (2017)
work page 2017
-
[4]
C. Broholm, R. J. Cava, S. A. Kivelson, D. G. Nocera, M. R. Norman, and T. Senthil, Science367, eaay0668 (2020), https://www.science.org/doi/pdf/10.1126/science.aay0668
- [5]
-
[6]
X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B39, 11413 (1989)
work page 1989
-
[7]
H. L. Stormer, D. C. Tsui, and A. C. Gossard, Rev. Mod. Phys. 71, S298 (1999)
work page 1999
-
[8]
K. Yang, L. K. Warman, and S. M. Girvin, Phys. Rev. Lett.70, 2641 (1993)
work page 1993
-
[9]
F. D. M. Haldane and D. P. Arovas, Phys. Rev. B52, 4223 (1995)
work page 1995
-
[10]
A. E. B. Nielsen, G. Sierra, and J. I. Cirac, Nat. Commun.4, 1 (2013)
work page 2013
-
[11]
Y .-C. He, D. N. Sheng, and Y . Chen, Phys. Rev. Lett.112, 137202 (2014)
work page 2014
-
[12]
S.-S. Gong, W. Zhu, and D. N. Sheng, Scientific Reports4, 6317 (2014)
work page 2014
- [13]
- [14]
-
[15]
W.-J. Hu, W. Zhu, Y . Zhang, S. Gong, F. Becca, and D. N. Sheng, Physical Review B91, 10.1103/physrevb.91.041124 (2015)
-
[16]
R. B. Laughlin, Phys. Rev. Lett.50, 1395 (1983)
work page 1983
-
[17]
W. Zhu, S. shu Gong, and D. N. Sheng, Proceedings of the National Academy of Sciences116, 5437 (2019), https://www.pnas.org/doi/pdf/10.1073/pnas.1807840116
-
[18]
S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Phys. Rev. Lett.54, 581 (1985)
work page 1985
-
[19]
S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Phys. Rev. B33, 2481 (1986)
work page 1986
- [20]
-
[21]
F. D. M. Haldane, Physical Review Letters107, 10.1103/phys- revlett.107.116801 (2011)
-
[22]
B. Yang, Z.-X. Hu, Z. Papi ´c, and F. D. M. Haldane, Phys. Rev. Lett.108, 256807 (2012)
work page 2012
-
[23]
S. Golkar, D. X. Nguyen, and D. T. Son, Journal of High Energy Physics2016, 10.1007/jhep01(2016)021 (2016)
-
[24]
Z. Liu, A. Gromov, and Z. Papi ´c, Phys. Rev. B98, 155140 (2018)
work page 2018
-
[25]
S.-F. Liou, F. D. M. Haldane, K. Yang, and E. H. Rezayi, Phys. Rev. Lett.123, 146801 (2019)
work page 2019
- [26]
-
[27]
Y . Liu, T. Zhao, and T. Xiang, Phys. Rev. B110, 195137 (2024)
work page 2024
-
[28]
I. V . Kukushkin, J. H. Smet, V . W. Scarola, V . Uman- sky, and K. von Klitzing, Science324, 1044 (2009), https://www.science.org/doi/pdf/10.1126/science.1171472
- [29]
-
[30]
C. Repellin, T. Neupert, Z. Papi ´c, and N. Regnault, Phys. Rev. B90, 045114 (2014)
work page 2014
- [31]
-
[32]
B. M. Kousa, N. Morales-Dur ´an, T. M. R. Wolf, E. Khalaf, and A. H. MacDonald, Theory of magnetoroton bands in moir´e ma- terials (2025), arXiv:2502.17574 [cond-mat.mes-hall]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[33]
M. Long, H. Lu, H.-Q. Wu, and Z. Y . Meng, Phys. Rev. B113, L041108 (2026)
work page 2026
- [34]
- [35]
-
[36]
H. B. Xavier, Z. Bacciconi, T. Chanda, D. T. Son, and M. Dal- monte, Phys. Rev. Lett.135, 196501 (2025)
work page 2025
- [37]
-
[38]
R. Samajdar, M. S. Scheurer, S. Chatterjee, H. Guo, C. Xu, and S. Sachdev, Nature Physics15, 1290 (2019)
work page 2019
- [39]
-
[40]
J. Yang, Z. Liu, and L. Wang, Phys. Rev. B110, 224404 (2024)
work page 2024
- [41]
-
[42]
F. Pichler, W. Kadow, C. Kuhlenkamp, and M. Knap, Phys. Rev. B110, 045116 (2024)
work page 2024
-
[43]
J. Willsher and J. Knolle, Dynamics and stability of u(1) spin liquids beyond mean-field theory: Triangular-latticej 1-j2 heisenberg model (2025), arXiv:2503.13831 [cond-mat.str-el]
-
[44]
Additional results in the Supplemental Information
- [45]
-
[46]
For cross-sector overlap estimates, the6×6cluster is not in- cluded due to the larger Hilbert-space dimension
-
[47]
F. D. M. Haldane, Phys. Rev. Lett.107, 116801 (2011)
work page 2011
- [48]
-
[49]
A. W. Sandvik, S. Capponi, D. Poilblanc, and E. Dagotto, Phys. Rev. B57, 8478 (1998)
work page 1998
-
[50]
T. P. Devereaux and R. Hackl, Rev. Mod. Phys.79, 175 (2007)
work page 2007
- [51]
-
[52]
L. J. P. Ament, M. van Veenendaal, T. P. Devereaux, J. P. Hill, and J. van den Brink, Rev. Mod. Phys.83, 705 (2011)
work page 2011
- [53]
-
[54]
L. Nie, A. V . Maharaj, E. Fradkin, and S. A. Kivelson, Phys. Rev. B96, 085142 (2017)
work page 2017
-
[55]
R. M. Fernandes, P. P. Orth, and J. Schmalian, Annual Review of Condensed Matter Physics10, 133 (2019)
work page 2019
- [56]
-
[57]
N. Regnault, J. Maciejko, S. A. Kivelson, and S. L. Sondhi, Phys. Rev. B96, 035150 (2017)
work page 2017
- [58]
-
[59]
S. Pu, A. C. Balram, J. Taylor, E. Fradkin, and Z. Papi ´c, Phys. Rev. Lett.132, 236503 (2024)
work page 2024
- [60]
-
[61]
H. Lu, H.-Q. Wu, B.-B. Chen, and Z. Y . Meng, Newton 10.1016/j.newton.2025.100300 (2025)
-
[62]
A. W. Sandvik, AIP Conference Proceedings1297, 135 (2010),
work page 2010
-
[63]
S. R. White, Phys. Rev. Lett.69, 2863 (1992)
work page 1992
-
[64]
J. Haegeman, J. I. Cirac, T. J. Osborne, I. Piˇzorn, H. Verschelde, and F. Verstraete, Phys. Rev. Lett.107, 070601 (2011)
work page 2011
- [65]
-
[66]
J. Hauschild, J. Unfried, S. Anand, B. Andrews, M. Bintz, U. Borla, S. Divic, M. Drescher, J. Geiger, M. Hefel, K. H ´emery, W. Kadow, J. Kemp, N. Kirchner, V . S. Liu, G. M ¨oller, D. Parker, M. Rader, A. Romen, S. Scalet, L. Schoonderwoerd, M. Schulz, T. Soejima, P. Thoma, Y . Wu, P. Zechmann, L. Zweng, R. S. K. Mong, M. P. Zaletel, and F. Pollmann, Sci...
work page 2024
- [67]
-
[68]
X.-L. Qi, H. Katsura, and A. W. W. Ludwig, Phys. Rev. Lett. 108, 196402 (2012). 8 The supplementary information contains: I. the ground-state entanglement spectrum; II. additional ED and TDVP results of the collective modes; III. additional information of TDVP simulations. I. GROUND-STA TE ENTANGLEMENT SPECTRUM Here, we show the entanglement spectrum from...
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.