{"id":"b6fed588-0c4c-4f9f-82a9-4777d0516886","arxiv_id":"1908.02271","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a pyrochlore Heisenberg antiferromagnet with weak bond disorder, the spin-glass Goldstone magnons are gauge-field excitations with speed c ~ aΔ, independent of the exchange energy J.","lead":"This paper predicts that the low-energy magnetic vibrations in a weakly disordered pyrochlore antiferromagnet travel at a speed set only by the disorder strength, not by the main magnetic exchange. The result unifies two separate descriptions of these materials: emergent gauge fields and spin-glass hydrodynamics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (15) rests on an unproven Δ/J estimate for the smooth magnetization; a direct ratio test is needed to confirm the c∼aΔ claim.","rationale":"The reader's weakest_assumption identifies exactly the same step: the estimate |m_J| ∼ (Δ/J)|m_Δ| in Section IV is asserted without a controlled derivation and is load-bearing for Eq. (15). My stress-test agrees with that identification rather than finding a separate flaw. The numerical evidence in Fig. 2(c) is genuinely supportive: the scaling collapse of (J/Δ)⟨m²_{J+}⟩^{1/2} versus ω/Δ is consistent with the Δ/J amplitude and would likely fail if the amplitude were O(1) or √(Δ/J). The paper also provides independent support through the Hessian and dynamical density scalings N(λ) ∝ (λ/Δ)^{3/2} and D(ω) ∝ (ω/Δ)^3, which match the linear-dispersion prediction. For these reasons I do not think the concern is severe enough to change the reader's ACCEPT verdict; it is a caveat, not a demonstrated contradiction. The proposed direct ratio measurement would settle the remaining ambiguity by removing normalization choices and isolating the coefficient of Eq. (15).","tokens_in":9388,"tokens_out":15313,"duration_ms":173644,"concrete_test":"On the same disorder realizations, for L = 7 and Δ/J = 2^-6, 2^-8, 2^-10, compute the lowest non-trivial dynamical eigenvector and decompose m_r into m_Δ (zero tetrahedron average) and m_J (smooth part) using the Hessian eigenvector projection described in Section V. Measure the ratio R = |m_J|/|m_Δ| for a normalized eigenvector and verify R ∝ Δ/J at fixed ω/Δ; equivalently, extract the coefficient c_m in ṁ_J = c_m ∇²θ_Δ from the linearized equations of motion and check that c_m = a²Δ²/J within a factor of two. If R is instead O(1) or ∝ √(Δ/J), the estimate fails and the speed reverts to a√(JΔ).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim c ∼ aΔ is controlled by Eq. (15), whose coefficient a²Δ²/J is obtained in Section IV from the estimate |m_r,J| ∼ (Δ/J)|m_r,Δ|, stated without derivation in the paragraph following Eq. (13). This estimate fixes the amplitude of the smooth magnetization that survives coarse-graining; if the true ratio were O(1) or O(√(Δ/J)) instead of O(Δ/J), the coefficient in Eq. (15) would become a²Δ or a²Δ^{3/2}/J, and the mode speed would revert to a√(JΔ) or an intermediate form. The numerical scaling collapse in Fig. 2(c) supports the estimate, but it plots (J/Δ)⟨m²_{J+}⟩^{1/2} at fixed ω/Δ and does not directly isolate the ratio |m_J|/|m_Δ| for a fixed eigenvector; it therefore does not fully close the gap left by the absence of a controlled expansion. Since Eq. (15) is the only place where the J dependence drops out of the mode speed, this is the most load-bearing step in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies magnon excitations in the classical Heisenberg antiferromagnet on the pyrochlore lattice with weak Gaussian exchange disorder, in the regime 0 < Δ ≪ J. Using the mapping of the clean ground-state manifold to divergenceless emergent U(1) gauge fields, and a decomposition of fluctuations into O(Δ) soft gauge-field modes and O(J) cluster-magnetization modes, the authors argue that the low-frequency Goldstone modes are smooth rotations of the frozen gauge-field configuration. Their central claim is that the mode speed is c ∼ aΔ, independent of J, in contrast to the conventional Halperin–Saslow result c ∼ a√(JΔ). The paper derives the scaling laws N(λ) ∼ (λ/Δ)^{3/2} and D(ω) ∼ (ω/Δ)^3, and supports them by numerical diagonalization on clusters with L = 3 to 7, including scaling collapses of eigenvector correlators and of the smooth magnetisation amplitude in Fig. 2(c).","tokens_in":9463,"tokens_out":11235,"duration_ms":120424,"significance":"If correct, the result identifies a genuinely new regime in which the emergent gauge-field stiffness is controlled by disorder alone, and it yields a concrete experimental fingerprint: a large magnetic heat capacity C ∼ T^3/Δ^3 at T ≪ T_F. The paper's strengths are that the central scaling laws are parameter-free predictions, that the numerical analysis tests eigenvector structure rather than only spectra, and that the prediction c ∼ aΔ is explicitly falsifiable. The contrast with the naïve extension of Halperin–Saslow theory, which gives c ∼ a√(JΔ), makes the claim physically significant for the frustrated-magnet community.","major_comments":[{"comment":"The coefficient a²Δ²/J in Eq. (15), and therefore the central prediction c ∼ aΔ, rests entirely on the estimate |m_{r,J}| ∼ (Δ/J)|m_{r,Δ}| for the amplitudes of the hard- and soft-subspace components of the magnetisation fluctuation. This estimate is asserted without derivation. The earlier result |M_α| ∼ Δ/J quoted in Section II concerns ground-state tetrahedron magnetisations, not the eigenvector decomposition used here, so it does not by itself justify Eq. (15). The numerical collapse in Fig. 2(c) is consistent with the estimate, but it plots (J/Δ)⟨m_{J+}²⟩^{1/2} against ω/Δ; it does not directly measure the ratio |m_J|/|m_Δ| for individual low-lying eigenvectors. If the true ratio were O(1) rather than O(Δ/J), Eq. (15) would have coefficient a²Δ and the mode speed would revert to a√(JΔ), invalidating the paper's headline claim. I ask the authors to derive this ratio from a controlled expansion in Δ/J, or to present a direct numerical measurement of |m_J|/|m_Δ| for fixed low-lying modes over the same range of Δ/J used in Fig. 2(c).","section":"Section IV, paragraph following Eq. (13), and Eq. (15)"},{"comment":"There is an apparent tension between the stiffness O(Δ) obtained for smooth rotations in Eq. (9) and the coefficient a²Δ²/J used in the equation of motion (15). A reader who combines the kinetic term θ̇_Δ ∼ J m_J from Eq. (14) with the natural canonical equation ṁ_J = −τ θ_Δ and τ ∼ −Δ a²∇² would obtain c ∼ a√(JΔ), which is precisely the conventional result the paper sets out to disprove. The paper should explain more explicitly how the elimination of m_Δ through the relation |m_J| ∼ (Δ/J)|m_Δ| renormalises the effective stiffness in the coarse-grained dynamics; as written, this is the single point where the J-dependence cancels and the argument is not fully transparent.","section":"Section IV, Eqs. (9) and (15)"}],"minor_comments":[{"comment":"The dashed line in Fig. 1(d) is labelled N(ω) ∝ (ω/Δ)^3, but the integrated dynamical density was defined as D(ω) in Section III; please make the notation uniform.","section":"Section V, caption of Fig. 1(d)"},{"comment":"The factor 2 and the proportionality constant in Eq. (11) are not defined; please state the disorder average and normalization explicitly, and comment briefly on the validity of replacing the pyrochlore ground-state average by a Gaussian ensemble of divergenceless fields.","section":"Section IV, Eq. (11)"},{"comment":"The ranges '2 −10 ≤ Δ/J ≤ 2−6' appear with missing superscripts in the text; please ensure the exponents are typeset consistently as 2^{-10} and 2^{-6}.","section":"Sections II and V"},{"comment":"The statement that D(ω) is convex and hence the measured exponent α decreases with increasing T toward T_F would benefit from one additional sentence, since the heat capacity involves an integral of D(ω) weighted by the Bose factor and the connection is not immediate.","section":"Section VI, discussion of heat capacity"},{"comment":"The notation O(Δ) after the first term of Eq. (6) should be clarified: it would help to state explicitly that this term collects corrections from the randomness to the Hessian in the basis of clean-system ground-state coordinates.","section":"Section III, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"This is a thought-provoking paper with a clear central claim and partly convincing numerical support. My main concern is the unproven amplitude estimate that controls Eq. (15); it is fixable within the scope of the manuscript, but I do not think the current text fully closes the gap. I would not reject on this basis, because the scaling collapse in Fig. 2(c) points in the right direction, but I would like the authors to either supply a derivation or add the direct eigenvector-ratio test."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper predicts that Goldstone modes in the pyrochlore Heisenberg spin glass have speed c ~ aΔ, independent of the mean exchange J, rather than the naive Halperin-Saslow result c ~ a√(JΔ). That is a genuinely new and testable claim, and the paper gives a coherent physical mechanism: the soft modes are excitations of emergent U(1) gauge fields with stiffness Δ, and the smooth magnetization that couples to twists is suppressed by Δ/J. The numerical data on L=3..7 clusters support the predicted scaling forms N(λ) ~ (λ/Δ)^{3/2} and D(ω) ~ (ω/Δ)^3, plus the eigenvector scaling collapse. For a scaling-level theory that is decent evidence.\n\nThe soft spot is exactly where the stress test points: Eq. (15) hangs on the estimate |m_{r,J}| ~ (Δ/J)|m_{r,Δ}|. The paper says this correction alone survives coarse-graining and asserts its size, but it does not derive it from a controlled expansion. Fig. 2(c) shows (J/Δ)<m^2_{J+}>^{1/2} collapsing as a function of ω/Δ, which is consistent with |m_J| ~ (Δ/J)|m_Δ|, but it does not isolate the ratio for a fixed eigenvector. So the central speed prediction rests on an order-of-magnitude estimate. That is not fatal — the paper is upfront about the level of rigor — but it is the one load-bearing step I would want a direct numerical ratio check for. The experimental heat-capacity and neutron predictions are clearly stated and would distinguish the theory.\n\nI think the paper deserves a serious referee. The derivation is heuristic but the physical argument is clear, the numerics are appropriate for the claim, and the prediction contradicts the standard theory in a way worth resolving. Read it if you work on frustrated magnets or spin-glass hydrodynamics; the gauge-field framing is useful even if the coefficient question remains open.","headline":"A scaling-level argument with supporting numerics that plausibly corrects the Halperin-Saslow mode speed for weakly disordered pyrochlore magnets; the load-bearing estimate is heuristic, not a fatal flaw.","tokens_in":10143,"tokens_out":1415,"would_cite":true,"duration_ms":15294,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper predicts that spin-glass Goldstone modes in a pyrochlore antiferromagnet travel at speed $c \\sim a\\Delta$, set by exchange disorder alone, not by the mean exchange $J$.","keywords":["frustrated magnetism","pyrochlore lattice","Heisenberg antiferromagnet","spin glass","emergent gauge field","Goldstone modes","exchange disorder","magnon dispersion"],"falsifier":"In a numerical simulation of the model, compute the lowest dynamical-mode frequency $\\omega$ at the smallest wavevector for several system sizes $L$ and disorder strengths $\\Delta/J$: the claim requires $\\omega/q\\to c\\propto\\Delta$, while a scaling $\\propto\\sqrt{J\\Delta}$ at the smallest $\\Delta/J$ would refute it. A thermodynamic measurement of the magnetic heat capacity below $T_F$, giving an exponent $2$ rather than $3$ on a three-dimensional pyrochlore, would also contradict the predicted density of states $D(\\omega)\\propto\\omega^3$.","tokens_in":9045,"feed_emoji":"🧲","tokens_out":6034,"duration_ms":58591,"temperature":0.7,"pith_summary":"Spin-glass order in a weakly disordered pyrochlore Heisenberg antiferromagnet has low-energy Goldstone modes that are best described as excitations of three emergent U(1) gauge fields. The paper derives a modified hydrodynamic theory in which smooth rotations of the frozen gauge-field configuration are conjugate to a smooth, small magnetization density, giving a linear dispersion with speed $c \\sim a\\Delta$ set only by the exchange disorder strength $\\Delta$. This replaces the conventional Halperin-Saslow prediction $c \\sim a\\sqrt{J\\Delta}$, which mixed the disorder scale with the exchange scale. The result ties a large low-temperature magnetic heat capacity to the smallness of $\\Delta$, and gives concrete neutron-scattering signatures.","feed_headline":"Spin-glass waves travel at a disorder-set speed","feed_subtitle":"In a pyrochlore antiferromagnet, Goldstone modes are emergent-gauge-field excitations, modifying Halperin-Saslow theory.","key_machinery":"The central object is the emergent divergenceless tensor field $B^{ai}_r = S^a_r e^i_r$, with spin index $a$ and spatial index $i$, whose flux through each tetrahedron is the cluster magnetization; in the frozen state it supplies three copies of a U(1) gauge field related by global spin rotations. The argument runs on a mode decomposition of fluctuations into a soft subspace with Hessian eigenvalues $O(\\Delta)$ (gauge-field twists) and a hard subspace with eigenvalues $O(J)$ (cluster-magnetization fluctuations). The load-bearing identities are the two coarse-grained equations of motion, $\\dot{\\theta}_\\Delta\\sim J m_J$ and $\\dot{m}_J\\sim a^2(\\Delta^2/J)\\nabla^2\\theta_\\Delta$, which combine a smooth rotation $\\theta_\\Delta$ with the surviving smooth component $m_J$ of the hard magnetization, whose magnitude is suppressed by $\\Delta/J$.","core_discovery":"The central claim is that in the regime $0<\\Delta\\ll J$, the lowest-energy magnons are smooth twists of the spin-glass state and are simultaneously excitations of the emergent gauge fields. The equations of motion are $\\dot{\\theta}_\\Delta \\sim J m_J$ and $\\dot{m}_J \\sim a^2 (\\Delta^2/J)\\nabla^2\\theta_\\Delta$, so their combination yields linearly dispersing Goldstone excitations with speed $c \\sim a\\Delta$, independent of $J$. Because the stiffness of the gauge field is of order $\\Delta$ and the smooth magnetization that couples to it is suppressed by $\\Delta/J$, the conventional hydrodynamic speed $\\sqrt{\\rho/\\chi_0}\\sim a\\sqrt{J\\Delta}$ is not realized; the paper argues that a naive extension mixes in high-energy degrees of freedom. Numerical data on Hessian eigenvalues, eigenvector correlators, and magnetization scaling support the description.","pith_inferences":["If the mechanism is generic, similar disorder-selected gauge-field states on other frustrated lattices should show the same $c\\sim a\\Delta$ scaling whenever the Hessian splits into $O(\\Delta)$ and $O(J)$ subspaces; the paper's own footnote indicates the two-dimensional case may acquire logarithmic corrections to the stiffness.","The predicted $C_M\\propto T^3$ could be tested in pyrochlore materials such as NaCaNi2F7, Y2Mo2O7, and Lu2Mo2O7 at temperatures well below their freezing transitions, where reported exponents are closer to 2; a crossover or sample-dependent exponent would clarify the role of disorder.","A direct test is to measure the dispersion slope $c$ as $\\Delta/J$ is tuned, for example by chemical pressure or strain: the paper's claim requires $c/\\Delta$ to be constant, while the conventional theory requires $c/\\sqrt{J\\Delta}$ to be constant."],"forward_implications":["The magnon speed vanishes linearly with disorder strength, so weakly disordered frustrated magnets should show very soft, gapless low-temperature excitations even when the exchange scale $J$ is large.","The magnetic heat capacity from these modes scales as $(T/\\Delta)^3$ below the freezing temperature, explaining large $C_M$ and predicting the exponent in clean samples.","Inelastic neutron scattering should show pinch-point correlations and, below $T_F$, a triple-peaked energy spectrum; above $T_F$ the same wavevector dependence persists with a Lorentzian lineshape.","The standard Halperin-Saslow formula $c=\\sqrt{\\rho/\\chi_0}$ does not apply here; the gauge-field stiffness is $\\rho\\sim\\Delta$ and the conjugate magnetization amplitude is what suppresses the speed.","The low-frequency sector is purely gauge-field-like as $\\Delta/J\\to0$, so the Goldstone modes coincide with the emergent gauge degrees of freedom."],"supporting_citations":[{"why":"Supplies the mapping of ground states to divergenceless emergent gauge fields on the pyrochlore lattice.","marker":"[2]"},{"why":"Establishes ground-state degeneracy and precessional dynamics of the clean pyrochlore Heisenberg antiferromagnet, giving the dynamically conjugate zero modes.","marker":"[14]"},{"why":"Provides the normal-mode formulation of the clean model that the weak-disorder treatment extends.","marker":"[15]"},{"why":"Shows that weak exchange randomness stabilizes a spin-glass state at $T_F\\sim\\Delta$ in this setting.","marker":"[19]"},{"why":"Analyzes the disordered pyrochlore Heisenberg model, supporting the spin-glass phase and disorder-selected flux configurations.","marker":"[20]"},{"why":"Halperin and Saslow: the conventional hydrodynamic theory of spin-glass Goldstone modes that the paper modifies.","marker":"[21]"},{"why":"Supplies the numerical energy-minimization protocol used to generate low-lying frozen states in the simulations.","marker":"[28]"}],"fun_headline_variants":["Spin-glass wave speed comes from disorder, not exchange","Emergent gauge fields drive spin-glass Goldstone modes","Frustrated magnet's spin waves rewrite Halperin-Saslow theory","Pyrochlore antiferromagnet: spin waves as gauge-field ripples","Spin-glass magnons tie to emergent gauge fields, not standard theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction $c\\sim a\\Delta$ rests on the estimate that the smooth part of the magnetization in a dynamical mode is of order $\\Delta/J$ rather than order one; if that amplitude were order one instead, the conventional $a\\sqrt{J\\Delta}$ speed would return.","fun_headline_variants_meta":{"raw":{"variants":["Spin-glass wave speed comes from disorder, not exchange","Emergent gauge fields drive spin-glass Goldstone modes","Frustrated magnet's spin waves rewrite Halperin-Saslow theory","Pyrochlore antiferromagnet: spin waves as gauge-field ripples","Spin-glass magnons tie to emergent gauge fields, not standard theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000362,"raw_usage":{"total_tokens":1899,"prompt_tokens":834,"completion_tokens":1065,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":450,"completion_tokens_details":{"reasoning_tokens":975}},"tokens_in":450,"tokens_out":1065,"duration_ms":10598,"temperature":1.0,"reasoning_tokens":975,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:49:29.216233+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a numerical simulation of the model, compute the lowest dynamical-mode frequency $\\omega$ at the smallest wavevector for several system sizes $L$ and disorder strengths $\\Delta/J$: the claim requires $\\omega/q\\to c\\propto\\Delta$, while a scaling $\\propto\\sqrt{J\\Delta}$ at the smallest $\\Delta/J$ would refute it. A thermodynamic measurement of the magnetic heat capacity below $T_F$, giving an exponent $2$ rather than $3$ on a three-dimensional pyrochlore, would also contradict the predicted density of states $D(\\omega)\\propto\\omega^3$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Halperin and Saslow: the conventional hydrodynamic theory of spin-glass Goldstone modes that the paper modifies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the numerical energy-minimization protocol used to generate low-lying frozen states in the simulations."}],"review_version":1}