REVIEW 2 major objections 4 minor 6 cited by
In two-dimensional Heisenberg spin glasses, free spin waves around classical glass states are weakly localized, so hydrodynamics requires spin-wave interactions.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
Leading-order spin waves of the 2D Heisenberg spin glass are weakly localized in the glass phase but delocalized at low energy in magnetically ordered phases, so hydrodynamics needs interactions.
T0 review reviewed 2026-07-13 challenge →
load-bearing objection Solid free-spin-wave localization dichotomy (SG weak localization vs FM/AFM mobility edge) with coherent RG/symmetry reading; hydrodynamics restoration is openly conjectural. the 2 major comments →
Semiclassical picture of the Heisenberg spin glass in two dimensions: from weak localization to hydrodynamics
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
Despite belonging to the same Altland–Zirnbauer class D, the localization of spin waves is controlled by the renormalization-group relevance of classical-order correlations in the Bogoliubov–de Gennes matrix: residual SO(2) order produces a mobility edge with delocalized low-energy modes, while spin-glass order renders the correlations irrelevant and the entire spectrum weakly localized.
What carries the argument
The quadratic spin-wave Hamiltonian obtained by Holstein–Primakoff expansion about classical minima, whose correlated matrix elements encode the broken-symmetry pattern of the classical state and thereby select the universality class of the localization problem.
Load-bearing premise
That the classical minima found by the mixed Monte Carlo algorithm are true global minima whose local spin orientations correctly encode the only correlations that matter for the spin-wave localization class.
What would settle it
A direct computation of the participation entropy or level-spacing ratio for the true S=1/2 Heisenberg spin glass that shows either a mobility edge inside the glass phase or a finite-size scaling of the localization length incompatible with the class-D beta function.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the 2D Heisenberg spin glass via a leading-order Holstein–Primakoff expansion around classical minima obtained by mixed Monte Carlo. It constructs the quadratic Bogoliubov–de Gennes spin-wave Hamiltonian (class D), carefully treats zero modes, and diagnoses localization with the r-parameter and a Bogoliubov-adapted fractal dimension. The central claim is that classical-order correlations are RG-irrelevant in the SG phase, so the entire spectrum is weakly localized (one-parameter scaling with α_D≈2 and ξ∼exp(ω₀²/ω²)), whereas residual SO(2) correlations in the FM/AFM phases are relevant and produce a low-energy mobility edge. Free spin waves therefore cannot realize Halperin–Saslow hydrodynamics in the SG phase; interactions at higher 1/S are conjectured to restore ergodicity because the localization length is unbounded as ω→0.
Significance. If the free-spin-wave dichotomy holds, the work cleanly separates the non-interacting localization problem from the interacting hydrodynamics problem for 2D Heisenberg spin glasses and supplies a concrete RG explanation for why the same AZ class can yield opposite outcomes. Strengths include careful zero-mode handling (App. B), a physically motivated definition of the participation measure for Bogoliubov modes (App. C), quantitative data collapses for ω_c and ν, and consistency checks of one-parameter scaling (α_D≈2, relation among α_D, α_ϕ, γ). The hydrodynamics-restoration argument is explicitly labeled conjectural and is not required for the free-wave claim. The results are of clear interest to the disordered quantum magnetism and localization communities.
major comments (2)
- Sec. II.A and the subsequent localization classification rest on the mixed microcanonical–canonical Monte Carlo minima (plus Hessian check) being representative global minima whose local rotation matrices encode the only relevant correlations in the BdG matrix. The manuscript itself notes that the algorithm is trapped by an exponential number of metastable states and defers a full analysis to a follow-up [21]. Because the free-wave claim (SG weak localization vs FM/AFM mobility edge) is load-bearing on those correlations being RG-irrelevant/relevant, the paper needs either stronger evidence that the sampled states are typical of the true ground-state ensemble (e.g., energy histograms, overlap distributions across independent runs, or comparison with known classical SG benchmarks) or an explicit statement of the limitation and a robustness check on a controlled subset of lower-energy conf
- Sec. IV.C–D and App. E: the one-parameter scaling analysis extracts α_D≈2, α_ϕ∈[1.2,1.5], γ≈2.6 and c_1≈1.3 from finite-size data up to L=60. While the internal consistency among these exponents is reassuring, the accessible sizes remain modest for a 2D class-D problem with only logarithmic localization. A clearer quantification of residual finite-size corrections (or an additional larger-L check for at least one energy window) would strengthen the claim that the correlations are truly irrelevant and that the system sits in the standard class-D flow rather than a long crossover.
minor comments (4)
- Fig. 2 and App. D: the mobility-edge points ω_c(p) are central to the phase diagram; a short table of the fitted (ω_c, ν) pairs for the p values shown would improve reproducibility.
- Eqs. (19)–(20) and App. C: the action-based definition of R_α(i)² is well motivated, but a one-sentence comparison with the conventional |X|²+|Y|² weight (already used for numerics) would help readers less familiar with bosonic BdG localization.
- Sec. V is clearly labeled conjectural; a brief remark that the decay-rate estimates (Eqs. 31–32) are order-of-magnitude only would further prevent over-reading of the hydrodynamics-restoration scenario.
- Typographical: “N´ eel” and “N` eel” appear inconsistently; “spin-w ave” hyphenation in headings; a few missing spaces after commas in the abstract and introduction.
Circularity Check
No significant circularity: free-spin-wave localization is computed directly from the constructed Bogoliubov matrices and compared to external AZ class-D / one-parameter-scaling benchmarks; self-citation to prior work supplies only phase-diagram context.
full rationale
The derivation chain is: (i) locate classical minima by mixed Monte-Carlo + Hessian positivity, (ii) expand the Heisenberg Hamiltonian to quadratic order in Holstein-Primakoff bosons, obtaining a correlated Bogoliubov-de Gennes matrix whose entries encode the classical rotation matrices, (iii) diagonalize and measure r-parameter and fractal dimension of the resulting spectrum, (iv) extract beta-functions and check consistency with the known one-parameter scaling of AZ class D (or with residual SO(2) relevance in the ordered phases). None of these steps is definitional: the localization diagnostics are independent numerical observables, the scaling exponents (nu, alpha_D, gamma, omega_0) are fitted and then cross-validated against each other and against external field-theory predictions, and the hydrodynamics-restoration argument of Sec. V is explicitly labeled a conjecture. The only self-citation ([7]) is used for the classical/quantum phase boundaries and for the accuracy of the 1/S order-parameter correction; it does not define or force the localization classification. Hence the central free-wave dichotomy stands on independent computation rather than on a circular reduction.
Axiom & Free-Parameter Ledger
free parameters (5)
- mobility edge ωc(p) and critical exponent ν =
p=0.1: ωc≈3.18, ν≈2.1; Dc≈0.75
- scaling exponents αD, αϕ, γ and prefactor c1 =
αD≈2 (1.9(2)); αϕ∈[1.2,1.5]; γ=2.6(3); c1=1.3(1)
- localization-length scale ω0 =
ω0=3.7(1)
- classical phase boundaries p1, p2 =
p1≃0.2, p2≃0.8
- Monte Carlo stopping threshold and replica protocol for classical minima
axioms (6)
- domain assumption Leading-order Holstein–Primakoff 1/S expansion around classical minima yields a faithful quadratic spin-wave Hamiltonian for localization of low-energy excitations.
- domain assumption One-parameter scaling and AZ class D β-function govern the 2D bosonic BdG problem when classical-order correlations are irrelevant.
- domain assumption Positive-semidefinite Hessian classical configurations found by mixed Monte Carlo are the appropriate expansion points (global or representative minima).
- ad hoc to paper Zero modes of the dynamical matrix may be treated as c-numbers set to zero for thermodynamic-limit observables without changing localization diagnostics.
- ad hoc to paper Interactions at O(1/√S) restore ergodicity/hydrodynamics because ξ(ω) is unbounded as ω→0, precluding a BAA-style MBL phase in d=2.
- standard math Standard linear algebra / symplectic Bogoliubov diagonalization and spectral statistics (Poisson vs GUE r-values).
Cite this review
Pith. "Pith review of Semiclassical picture of the Heisenberg spin glass in two dimensions: from weak localization to hydrodynamics." pith.science (2026). https://pith.science/paper/GH6XYPRM
@misc{pith2026260322077,
author = {Pith},
title = {Pith review of: Semiclassical picture of the Heisenberg spin glass in two dimensions: from weak localization to hydrodynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/GH6XYPRM}},
note = {Machine review of arXiv:2603.22077}
}
read the original abstract
The two-dimensional Heisenberg spin-glass model is investigated by means of a semiclassical expansion around classical states. At leading order, we obtain an effective quadratic spin-wave Hamiltonian and study the localization properties of its spectrum and eigenfunctions. We find that the nature of the spin-wave excitations, whether they are hydrodynamic or localized modes, depends crucially on the relevance/irrelevance -- in the renormalization group sense -- of the correlations induced by the underlying classical order in the spin-wave Hamiltonian matrix elements: low-energy excitations around magnetically ordered states are delocalized, whereas those around spin-glass ordered states are localized, albeit weakly. We interpret this phenomenology by relating the spontaneous breaking of spin-rotation symmetry in the original Heisenberg model to the symmetry and universality class of the resulting quadratic spin-wave Hamiltonian. We conjecture that the hydrodynamic picture can be recovered through the inclusion of interactions among the spin-wave excitations at higher order in the semiclassical expansion, favoring the onset of ergodic behavior.
Forward citations
Cited by 6 Pith papers
-
Vanishing spin stiffness in weakly disordered two-dimensional Heisenberg ferromagnets
Weak bond frustration in 2D Heisenberg ferromagnets produces logarithmically correlated spin-stiffness fluctuations whose RG flow drives stiffness to zero, giving soft magnons with scale-dependent z>2.
-
A New Robust Constraint on the Self-interaction Cross-section of Dark Matter with Double Radio Relic Clusters
Shock-to-shock distances in double radio relic clusters serve as a merger chronometer to set a 68% upper limit of sigma/m < 0.22 cm²/g on dark matter self-interaction, marginalizing over mass, angle, speed, phase, imp...
-
A New Robust Constraint on the Self-interaction Cross-section of Dark Matter with Double Radio Relic Clusters
Variance of graph-energy centrality jumps at known weak ETH-breaking transitions in RPM and QSM and flags a glassy crossover in the triangular lattice gas.
-
Graph-theory measures capture weak ergodicity breaking on large quantum systems
Graph-energy centrality applied to Fock-space graphs captures weak ergodicity-breaking transitions in quantum many-body systems and scales to hundreds of sites or the thermodynamic limit.
-
Quantum spin-glass criticality in disordered frustrated dimer magnets
In a disordered triangular bilayer Heisenberg magnet, theory predicts a quantum spin glass whose near-critical order is sparse and nearly collinear — 'doubly weak' — with strongly suppressed amplitude (Higgs) mode weight.
-
Graph-theory measures capture weak ergodicity breaking on large quantum systems
Graph-energy centrality detects weak ergodicity-breaking transitions in large quantum many-body systems via changes in its distribution and applies to kinetically constrained models showing glassy dynamics.
Reference graph
Works this paper leans on
-
[1]
high disorder trap
and write the beta function as βt = dt dlnL 2 ,(21) whereLis the largest scale in the problem (i.e.the sys- tem size). In the present case, which falls in the AZ symmetry class D, the conductance is the thermal one, since the particle number is not conserved,i.e.there is no electric charge. A key consequence of the one-parameter scaling hypothesis is that...
2022
-
[2]
B. I. Halperin and W. M. Saslow, Phys. Rev. B16, 2154 (1977)
1977
-
[3]
B. I. Halperin and P. C. Hohenberg, Phys. Rev.188, 898 (1969)
1969
-
[4]
P. W. Anderson, Phys. Rev.109, 1492 (1958)
1958
-
[5]
D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Annals of physics321, 1126 (2006)
2006
-
[6]
Evers and A
F. Evers and A. D. Mirlin, Rev. Mod. Phys.80, 1355 (2008)
2008
-
[7]
Sierant, M
P. Sierant, M. Lewenstein, A. Scardicchio, L. Vidmar, and J. Zakrzewski, Reports on Progress in Physics88, 026502 (2025)
2025
-
[8]
L. L. Viteritti, R. Rende, G. Bracci-Testasecca, J. Niedda, R. Moessner, G. Carleo, and A. Scardicchio, arXiv:2507.05073 (2024)
arXiv 2024
-
[9]
A. J. Bray and M. A. Moore, Journal of Physics C: Solid State Physics14, 2629 (1981)
1981
-
[10]
P. W. Anderson, Phys. Rev.86, 694 (1952)
1952
-
[11]
Oguchi, Phys
T. Oguchi, Phys. Rev.117, 117 (1960)
1960
-
[12]
G. E. Castilla and S. Chakravarty, Phys. Rev. B43, 13687 (1991)
1991
-
[13]
Altland and M
A. Altland and M. R. Zirnbauer, Phys. Rev. B55, 1142 (1997)
1997
-
[14]
Gurarie and J
V. Gurarie and J. T. Chalker, Phys. Rev. B68, 134207 (2003)
2003
-
[15]
A. C. Potter, R. Vasseur, and S. A. Parameswaran, Phys. Rev. X5, 031033 (2015)
2015
-
[16]
De Roeck and F
W. De Roeck and F. Huveneers, Phys. Rev. B95, 155129 (2017)
2017
-
[17]
Thiery, F
T. Thiery, F. Huveneers, M. M¨ uller, and W. De Roeck, Phys. Rev. Lett.121, 140601 (2018)
2018
-
[18]
Morningstar, L
A. Morningstar, L. Colmenarez, V. Khemani, D. J. Luitz, and D. A. Huse, Phys. Rev. B105, 174205 (2022)
2022
-
[19]
J. F. Fernandez, Journal of Physics C: Solid State Physics 10, L441 (1977)
1977
-
[20]
Baity-Jesi, V
M. Baity-Jesi, V. Mart´ ın-Mayor, G. Parisi, and S. Perez- Gaviro, Phys. Rev. Lett.115, 267205 (2015)
2015
-
[21]
Baity-Jesi and G
M. Baity-Jesi and G. Parisi, Phys. Rev. B91, 134203 (2015)
2015
-
[22]
Coraggio, J
A. Coraggio, J. Niedda, G. Bracci-Testasecca, and A. Scardicchio,in preparation(2026)
2026
-
[23]
Binder and A
K. Binder and A. P. Young, Rev. Mod. Phys.58, 801 (1986)
1986
-
[24]
In fact, 1 N 2 X i,j X α,β Sα i Sα j Sβ i Sβ j =S 4 X α,β 1 3 δαβ 1 3 δαβ = S4 3
-
[25]
One can obtain an AFM minimum of the classical en- ergy from a FM one, by performing a staggered rotation, i.e.flipping all the spins on a sublattice, and changing sing to all the couplings [21], meaning that a minimum configuration atpcorresponds to another one at 1−p
-
[26]
Blaizot and G
J.-P. Blaizot and G. Ripka,Quantum Theory of Finite Systems(MIT Press, Cambridge, Massachusetts, 1986)
1986
-
[27]
Nielsen and S
H. Nielsen and S. Chadha, Nuclear Physics B105, 445 (1976)
1976
-
[28]
Watanabe and T
H. Watanabe and T. c. v. Brauner, Phys. Rev. D84, 125013 (2011)
2011
-
[29]
Watanabe and H
H. Watanabe and H. Murayama, Phys. Rev. Lett.108, 251602 (2012)
2012
-
[30]
L. D. Landau and E. M. Lifshitz,Statistical Physics, Part 1, Course of Theoretical Physics, Vol. 5(Pergamon Press, 1980)
1980
-
[31]
N. N. Bogoljubov, Journal of Physics (USSR)11, 23 (1947)
1947
-
[32]
Seibold, L
G. Seibold, L. Benfatto, C. Castellani, and J. Lorenzana, Phys. Rev. Lett.108, 207004 (2012)
2012
-
[33]
Seibold, L
G. Seibold, L. Benfatto, C. Castellani, and J. Lorenzana, Phys. Rev. B92, 064512 (2015)
2015
-
[34]
Q.-R. Xu, V. P. Flynn, A. Alase, E. Cobanera, L. Viola, and G. Ortiz, Phys. Rev. B102, 125127 (2020)
2020
-
[35]
M. L. Mehta,Random Matrices, 3rd ed. (Elsevier, 2004)
2004
-
[36]
Colpa, Physica A: Statistical Mechanics and its Appli- cations93, 327 (1978)
J. Colpa, Physica A: Statistical Mechanics and its Appli- cations93, 327 (1978)
1978
-
[37]
As such, each spin in a box of sizeλ gives to the effective field anO(1) contribution
The effective spin does not correspond to the coarse graining of a (disordered) spin configuration, but of the spins asdynamicalvariables, canonically conjugated to the rotation field. As such, each spin in a box of sizeλ gives to the effective field anO(1) contribution
-
[38]
Oganesyan and D
V. Oganesyan and D. A. Huse, Phys. Rev. B75, 155111 (2007)
2007
-
[39]
Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux, Phys. Rev. Lett.110, 084101 (2013)
2013
-
[40]
In fact, the same difficulties in studying eigenstate lo- calization in real space also are found in the many-body case, even if in that case they are more severe
From this point of view, Bogoljubov modes resemble many-body wavefunctions living in the Fock space, and this reflects the absence of particle number conservation. In fact, the same difficulties in studying eigenstate lo- calization in real space also are found in the many-body case, even if in that case they are more severe
-
[41]
Jiricek, M
S. Jiricek, M. Hopjan, V. Kravtsov, B. Altshuler, and L. Vidmar, Proceedings of the National Academy of Sci- ences of the United States of America123, e2518027123 (2026)
2026
-
[42]
C. Vanoni, B. L. Altshuler, V. E. Kravtsov, and A. Scardicchio, Proc. Natl. Acad. Sci. (USA)121, 10.1073/pnas.2401955121 (2024)
-
[43]
This is beyond the scope of this work and is left for following research
Notice: this is just a proxy ofz, whose correct value can be extracted from the correlation functions as a function of time. This is beyond the scope of this work and is left for following research
-
[44]
Abrahams, P
E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. Lett.42, 673 (1979)
1979
-
[45]
B. L. Altshuler, V. E. Kravtsov, A. Scardicchio, P. Sier- ant, and C. Vanoni, Proceedings of the National Academy of Sciences122, e2423763122 (2025)
2025
-
[46]
Niedda, G
J. Niedda, G. Bracci Testasecca, G. Magnifico, F. Bal- ducci, C. Vanoni, and A. Scardicchio, Phys. Rev. B112, 144201 (2025)
2025
-
[47]
Balducci, G
F. Balducci, G. Bracci-Testasecca, J. Niedda, A. Scardic- chio, and C. Vanoni, Phys. Rev. B111, 214206 (2025)
2025
-
[48]
To cut a long story short, we consider the situation set- tled by the theorems of [53] based on the work [54], but see also [49, 55] for the connections with the perturbation theory of [4]
-
[49]
S. John, H. Sompolinsky, and M. J. Stephen, Phys. Rev. B27, 5592 (1983)
1983
-
[50]
V. Ros, M. M¨ uller, and A. Scardicchio, Nuclear Physics B891, 420 (2015)
2015
-
[51]
De Roeck, F
W. De Roeck, F. Huveneers, M. M¨ uller, and M. Schiulaz, Physical Review B93, 014203 (2016)
2016
-
[52]
Chandran, A
A. Chandran, A. Pal, C. Laumann, and A. Scardicchio, 18 Physical Review B94, 144203 (2016)
2016
-
[53]
Eisenmenger and A
W. Eisenmenger and A. A. Kaplianski˘ ı,Nonequilibrium phonons in nonmetallic crystals., Modern problems in condensed matter sciences. v.16. (North-Holland, Ams- terdam, 1986)
1986
-
[54]
W. De Roeck, L. Giacomin, F. Huveneers, and O. Pros- niak, arXiv:2408.04338 (2024)
Pith/arXiv arXiv 2024
-
[55]
J. Z. Imbrie, Journal of Statistical Physics163, 998 (2016)
2016
-
[56]
J. Z. Imbrie, V. Ros, and A. Scardicchio, Annalen der Physik529, 1600278 (2017)
2017
This paper was first reviewed by grok-4.5 on July 13, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.