{"id":"628beb91-3d34-4c27-939b-73c85c7462f7","arxiv_id":"2501.01201","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Inelastic neutron scattering in the field-polarized state of GdRu2Si2 yields a magnon dispersion fitted by bilinear Heisenberg exchange interactions J1 through J8, with no sign of higher-order or anisotropic terms.","lead":"Neutron scattering on the skyrmion candidate GdRu2Si2 in its high-field state reveals a single magnon band described by eight bilinear Heisenberg exchange couplings, with no need for anisotropic or higher-order exchange terms. The result constrains the magnetic Hamiltonian of a material that forms zero-field multi-Q order, which is hard to reconcile with purely bilinear exchange.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Biquadratic and longitudinal anisotropic exchange are degenerate with bilinear Heisenberg couplings in the linear spin-wave dispersion, so the fitted magnon band cannot support the paper's bound of ≲1 µeV on such terms.","rationale":"The reader's CONDITIONAL verdict is well-founded, but the stress-test identifies a stronger, structural reason the null result cannot be concluded from the magnon-band dispersion alone. In linear spin-wave theory around a fully polarized ferromagnet, the quadratic Hamiltonian is sensitive only to the Fourier transform of the transverse exchange couplings plus a constant term. Biquadratic exchange (S_i·S_j)^2 expands to a term proportional to the same boson operator as bilinear exchange, so it merely renormalizes the fitted J's; the longitudinal part of anisotropic exchange and single-ion anisotropy enter only through the same fitted constant C that already absorbs -gµ_BB and unaccounted exchange. Thus the statement 'no evidence for anisotropic or higher-order-exchange terms' is not a falsifiable outcome of the dispersion fit—any such terms can be absorbed. The paper's bound '≲1 µeV' is therefore an assertion, not a result of the analysis; it could in principle be obtained by comparing the extracted J_eff with the ab initio J_ij of Bouaziz et al., but the paper does not perform that comparison and instead attributes the differences to omitted long-range RKKY interactions. The proposed concrete test—demonstrating the analytic degeneracy and its consequence for the parameter covariance—would settle whether the data can support the bound, and it would force a reframing of the central claim as 'the dispersion is consistent with an effective bilinear model' rather than 'higher-order exchange is absent.' The verdict remains CONDITIONAL, since the effective model and the interaction-invariant method are still valuable contributions, but the advertised negative finding should be revised.","tokens_in":16592,"tokens_out":13740,"duration_ms":137997,"concrete_test":"Analytical check: derive the linear spin-wave expansion of H = Σ_ij [J_ij S_i·S_j + B_ij (S_i·S_j)^2] for a collinear ferromagnet and verify that the quadratic boson Hamiltonian depends only on the combination J_ij + 2S^2 B_ij. Then, as a numerical demonstration, take the fitted J1→8 from Table I, add a biquadratic term with B = 5 µeV on the J2 bonds, and refit an effective bilinear model; show that the same dispersion is recovered with a renormalized J2. More directly, fit the experimental dispersion points with a model that includes both J_ij and B_ij on all bonds and examine the covariance matrix: if J and B are perfectly collinear (a zero eigenvalue along dJ_ij = -2S^2 dB_ij), the data cannot bound B. This test settles whether the claim that biquadratic terms are 'limited to ≲1 µeV' is derivable from the fitting procedure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Linear spin-wave theory around the fully polarized state makes the central null result structurally untestable. Expanding a bond-biquadratic term B_ij (S_i·S_j)^2 to quadratic order in Holstein-Primakoff bosons gives -2B_ij S^3 (a_i†a_i + a_j†a_j - a_i†a_j - a_j†a_i), which is exactly the operator produced by a bilinear Heisenberg exchange with renormalized coupling J_ij + 2S^2 B_ij (S = 7/2, so 2S^2 = 24.5). Thus the measured dispersion ℏω(q) = S[J(q) - J(0)] + C determines only the effective transverse couplings; a model with no biquadratic exchange and one with arbitrarily large B_ij (compensated by different J_ij) produce identical spectra. A similar degeneracy holds for the longitudinal part of anisotropic exchange and for single-ion anisotropy, which enter only through the fitted constant C and are explicitly acknowledged as inseparable from -gµ_BB and unaccounted exchange. Consequently, the paper's statement that higher-order exchange terms are 'not required' is trivially true, and the inference that such terms are 'likely small and limited to ≲1 µeV' does not follow from the magnon-band fit. The bound could only come from an independent determination of the bare J_ij (e.g., from the ab initio values), but the paper does not convert the J_eff - J_abinitio differences into biquadratic bounds. The interaction-invariant method does not resolve this degeneracy, since it operates on the same dispersion formula.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports inelastic neutron scattering measurements of the magnon dispersion in the field-polarized phase of GdRu2Si2 at T = 25 K and B = 8.8 T. The authors fit the dispersion using linear spin-wave theory with an effective bilinear Heisenberg exchange model including interactions J1 through J8 (with J7 set to zero) plus a constant term C. They introduce a 'method of interaction invariant path analysis' to select reciprocal-space cuts that decorrelate the exchange parameters. The fitted model reproduces the main features of the measured dispersion and is qualitatively consistent with earlier ab initio calculations. The paper concludes that no anisotropic or higher-order exchange terms are required, and estimates that such terms are likely limited to about 1 micro-eV. The implications for the zero-field multi-Q ground state are discussed.","tokens_in":17027,"tokens_out":9246,"duration_ms":79005,"significance":"If the central null result were established, the paper would be significant: it would show that bilinear RKKY exchange alone describes the spin dynamics in the field-polarized state, and that higher-order exchange is negligible, which bears directly on the mechanism stabilizing multi-Q order in this material. The experimental dataset is valuable, and the interaction invariant path analysis provides a systematic way to reduce parameter correlations in magnon fits. The paper is commendably transparent about its limitations, including low signal-to-noise, median kernel filtering, and the unresolved bottom of the magnon band. However, as detailed in the major comments, the central null result is not actually testable with the present linear spin-wave analysis, so the paper's headline conclusion overreaches the data.","major_comments":[{"comment":"The statement in the abstract and Discussion that 'No evidence for anisotropic or higher-order-exchange terms beyond bilinear Heisenberg exchange is found,' and the quantitative bound that such terms are 'likely small and limited to ≲ 1 µeV,' are not supported by the analysis. In linear spin-wave theory around a field-polarized state, a biquadratic term B_ij (S_i·S_j)^2 expands to quadratic order as -2 B S^3 (a_i†a_i + a_j†a_j - a_i†a_j - a_j†a_i), which is exactly the same operator form as the quadratic part of a bilinear Heisenberg coupling J_ij with the replacement J_eff = J_ij + 2S^2 B_ij. With S = 7/2, 2S^2 = 24.5, so the fitted exchange constants in Table I are effective transverse couplings that cannot distinguish biquadratic exchange from bilinear exchange. A similar degeneracy holds for the longitudinal component of anisotropic exchange and for single-ion anisotropy, which enter only through the fitted constant C in Eq. (2), as the paper itself notes. Consequently, the observation that the dispersion can be fitted without higher-order terms is trivially true, and the bound of ≲1 µeV does not follow from the data. To support such a bound, the authors would need to compare the fitted effective J values with independently known bare J values and quantitatively estimate how much of the difference could be attributed to biquadratic terms; the paper does not do this.","section":"Discussion; Eq. (2)"},{"comment":"The inclusion of J8 is explicitly motivated by the need to place the global minima of the dispersion at q = [0.15, 0, 0], i.e., in the experimentally observed propagation direction. Therefore the agreement between the model's minimum and the lower-field propagation vector is a fitting constraint, not an independent prediction. The paper also notes that the minimum is at 0.15 r.l.u. rather than the observed q_e = [0.22, 0, 0]. This limits the strength of the claim that the model 'captures the salient features of the magnon dispersion, including global minima at incommensurate positions which are characteristic of the lower-field magnetism.'","section":"Analysis; Table I"},{"comment":"The paper does not provide a sensitivity analysis or a model comparison that would justify the quantitative bound on higher-order terms. The data have low signal-to-noise, are processed with a nonlinear median kernel filter, and the bottom of the band below ~0.8 meV is not resolved. Without fitting models that explicitly include biquadratic or anisotropic terms and assessing their statistical significance (e.g., via a likelihood-ratio test or error propagation), the estimate that such terms are ≲1 µeV is an assertion rather than a derived result. The paper should either remove this bound or replace it with a statement that the data are consistent with zero within the sensitivity of the measurement, and provide an estimate of that sensitivity.","section":"Analysis; Supplemental Material"}],"minor_comments":[{"comment":"The name of the spin-wave code is written inconsistently: 'Sun(n)y' in the main text and 'Su(n)ny' in the Supplemental Material. It is likely 'Sunny'; please standardize.","section":"Supplemental Material"},{"comment":"There is a typo in Section II: 'consistent the the S(Q,ω ) calculations' should be 'consistent with the S(Q,ω ) calculations'.","section":"Supplemental Material"},{"comment":"The phrase 'Gd 3+ and Eu 2+ intermetallics' has spacing issues; please correct the formatting.","section":"Introduction"},{"comment":"The sentence 'A Hamiltonian accounting for the dispersion relation has been derived with linear spin wave theory' could be phrased more clearly: the Hamiltonian was fitted to the dispersion, not derived from it.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The paper is from an established experimental group and contains a unique dataset that is difficult to obtain. The main concern is that the central conclusion about higher-order exchange is not supported by the linear spin-wave analysis due to the bilinear-biquadratic degeneracy. The authors should revise the paper to explicitly acknowledge this degeneracy, remove the quantitative bound, and reframe the conclusion as consistency with an effective bilinear model. With these changes, the paper would be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper reports field-polarized single-crystal INS magnon data for GdRu2Si2 and fits a J1–J8 bilinear Heisenberg exchange model. The data itself is the new thing — no prior single-crystal magnon measurement in the polarized state. The interaction invariant path analysis is a genuine methodological contribution; choosing cuts where a given interaction is non-dispersive is a clever way to decorrelate exchange parameters, and the supplementary material explains it carefully. The ab initio comparison with Bouaziz et al. is useful and gives an external anchor. The authors are also commendably explicit about the low signal-to-noise, median filtering, unresolved band bottom, and the C degeneracy.\n\nThe soft spot is central. The paper's headline null result — no anisotropic or higher-order exchange, with an estimated upper bound of ≲1 µeV — is not something the magnon band fit can actually constrain. Linear spin wave theory around the fully polarized state has a symmetry: a biquadratic term (S_i·S_j)^2 produces, at quadratic order, exactly the same operator as a bilinear Heisenberg coupling with renormalized J_ij. With S = 7/2, the renormalization factor is 24.5, so any biquadratic B can be traded against J. The dispersion ω(q) = S[J(q) − J(0)] only determines the effective transverse couplings; it cannot distinguish the source. The same is true for the longitudinal part of anisotropic exchange and for single-ion anisotropy, which enter only through the constant C and are acknowledged as inseparable from the Zeeman term. So the statement \"no evidence for higher-order terms\" is trivially true in this fitting scheme, and the ≲1 µeV bound does not follow from the data. A real bound would need an independent determination of bare J_ij, e.g. from ab initio, and then a comparison of J_eff − J_ab initio converted into B_ij. The paper does not do that. This doesn't invalidate the measured dispersion or the effective Hamiltonian as a parametrization, but it does mean the paper's key interpretive claim about the absence of higher-order exchange is unsupported.\n\nAlso worth noting: the model is fitted, not predictive; the bottom of the band is unresolved below ~0.8 meV; the model accounts for 72% of the Curie-Weiss temperature; and J8 is added precisely to put the minima in the observed propagation direction. None of these are fatal, but together they mean the Hamiltonian is a sensible effective model, not a uniquely established one.\n\nWho is this for? Experimentalists and theorists working on GdRu2Si2 and other centrosymmetric skyrmion hosts. The data and method deserve a serious referee. My recommendation: send it out, but the referees should require the authors to temper the higher-order exchange claim — either drop the ≲1 µeV bound or support it with an explicit calculation.","headline":"Valuable new magnon data and a useful fitting method, but the paper's central null result on higher-order exchange is not actually testable by the magnon band fit.","tokens_in":17560,"tokens_out":2433,"would_cite":true,"duration_ms":24005,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.30.Ds","75.40.Gb","78.70.Nx"],"model":"deepseek-v4-flash","headline":"The spin-wave spectrum of the skyrmion candidate GdRu2Si2 in its field-polarized phase is reproduced by eight bilinear Heisenberg exchange interactions, with no evidence for higher-order or anisotropic exchange.","keywords":["GdRu2Si2","magnon dispersion","linear spin wave theory","Heisenberg exchange","RKKY interaction","multi-Q magnetic order","inelastic neutron scattering","skyrmion candidate"],"falsifier":"A measurement of the magnon dispersion below 0.8 meV—for example with a cold-neutron triple-axis spectrometer as the authors propose—would settle the claim: if the bottom of the band deviates from the $J_1$-$J_8$ bilinear prediction, or if a fit with a biquadratic exchange term significantly outperforms the bilinear model, the central null claim is refuted. A simpler check on existing data is to add an $(\\mathbf{S}_i\\cdot\\mathbf{S}_j)^2$ term to the fit and test whether the residuals drop meaningfully.","tokens_in":16405,"feed_emoji":"🧲","tokens_out":8176,"duration_ms":64981,"temperature":0.7,"pith_summary":"This paper measures the magnon dispersion of the centrosymmetric skyrmion candidate GdRu2Si2 in the field-polarized phase using inelastic neutron scattering, and shows that the spectrum can be accounted for by a Hamiltonian containing eight bilinear Heisenberg exchange interactions, the Zeeman term, and negligible single-ion anisotropy. The authors find no evidence for anisotropic or higher-order (e.g., biquadratic) exchange, and argue that if such terms exist they are limited to roughly 1 micro-eV. This result matters because anisotropic and higher-order exchange are usually invoked to stabilize the multi-Q and skyrmion states that this material hosts in zero field; the paper's finding sharpens the puzzle of what actually stabilizes those states. The fitted exchange parameters are consistent with long-ranged RKKY interactions and with ab initio calculations, and the dispersion's incommensurate minima point toward the propagation direction of the lower-field multi-Q order.","feed_headline":"No exotic exchange terms needed for GdRu2Si2 spin waves","feed_subtitle":"Inelastic neutron scattering shows that eight bilinear Heisenberg interactions fully account for the magnon band.","key_machinery":"The central object is the linear spin-wave dispersion $\\hbar\\omega(\\mathbf{q}) = S[J(\\mathbf{q})-J(0)] - g\\mu_B B - 2KS$, obtained from a Holstein-Primakoff transformation of a Hamiltonian with bilinear Heisenberg exchange, Zeeman, and single-ion anisotropy terms. The paper introduces 'interaction invariant path analysis': reciprocal-space cuts are chosen along which a particular interaction, such as $J_2$ with $\\delta=[1/2,1/2,1/2]$, is non-dispersive (its cosine terms cancel), so that other weaker interactions can be refined independently. This method guides the selection of fits, and the final $J_1$-$J_8$ model is fitted to the extracted band by least-squares minimization of $\\hbar\\omega(\\mathbf{q})$.","core_discovery":"Using time-of-flight inelastic neutron scattering on an isotopically enriched single crystal, the spin excitations of GdRu2Si2 in the field-polarized (forced ferromagnetic) phase were measured along a 90-degree arc of reciprocal space. After median-kernel filtering to extract the magnon band from strong absorption-limited noise, the dispersion was fitted with linear spin wave theory. A model with eight bilinear Heisenberg exchange interactions, $J_1$ through $J_8$ (with $J_2$ dominant at $-65.1$ $\\mu$eV), plus a constant term consistent with the applied field, reproduces the measured dispersion. The key result is a null result: no anisotropic or higher-order exchange terms are required to fit the data, and the authors estimate that if present these terms are smaller than about $1$ $\\mu$eV. The fitted model places the global dispersion minima at incommensurate wavevectors along the $[1,0,0]$ direction, consistent in direction with the magnetic propagation vector of the lower-field multi-Q state, although the exact minimum position is uncertain because the bottom of the band lies below the energy resolution.","pith_inferences":["A direct test of the null claim is to fit existing or future higher-resolution data with models that include a biquadratic term $(\\mathbf{S}_i\\cdot\\mathbf{S}_j)^2$ or anisotropic exchange and compare the goodness of fit; a significant improvement would falsify the 'bilinear-only' conclusion.","The paper's own Curie-Weiss comparison shows the fitted model accounts for only about 72% of the net ferromagnetic molecular field, so the $J_1$-$J_8$ set is an effective truncation of a longer-range exchange tail; the 'no higher-order terms' conclusion should be read as applying to the resolvable portion of the spectrum, not as a complete Hamiltonian for the material.","The method of interaction invariant path analysis could be transferred to other strongly absorbing or low-signal magnets where full $\\mathbf{S}(\\mathbf{Q},\\omega)$ fitting is impractical, and the authors note plans to extend it to non-linear invariant loops."],"forward_implications":["If the Hamiltonian is correct, the spin dynamics of the field-polarized phase of GdRu2Si2 is governed by long-ranged bilinear RKKY exchange, and higher-order or anisotropic terms play no measurable role there.","The fitted exchange parameters support the Fermi-surface-based ab initio model cited in the paper, lending weight to the view that the exchange is RKKY-mediated.","The incommensurate minima of the dispersion along the [1,0,0] direction show that the tendency toward the lower-field multi-Q order is already present in the bilinear exchange; however, the exact minimum position differs from the observed propagation vector, suggesting that missing long-range exchange terms are important.","The absence of observable anisotropic or higher-order terms implies that the zero-field multi-Q state must be stabilized either by field- or temperature-dependent RKKY interactions tied to the Fermi surface, or by coupling to charge-density-wave modulations discussed in the paper."],"supporting_citations":[{"why":"Supplies the ab initio RKKY exchange parameters and Fermi-surface mechanism that the fitted model is compared against.","marker":"[29]"},{"why":"Determines the zero-field double-Q ground state and the magnetic propagation vector used to assess the model's minima.","marker":"[18]"},{"why":"Zero-field powder spin dynamics study that corroborates the ab initio model and the absence of higher-order exchange.","marker":"[54]"},{"why":"Documents the multi-Q phase zoology and field-direction anisotropy that motivate the search for anisotropic terms.","marker":"[17]"},{"why":"Supplies the theorem that bilinear Heisenberg exchange alone stabilizes single-Q order, motivating the need for additional terms in multi-Q states.","marker":"[32]"},{"why":"Extends the theorem on ground-state spin configurations relevant to the implication about multi-Q stability.","marker":"[33]"},{"why":"Quantum oscillations evidence for Fermi-surface reconstruction between field-polarized and multi-Q states, supporting the field-dependent RKKY proposal.","marker":"[41]"},{"why":"Shows coupling between itinerant electrons and local moments and charge density waves, used to argue for exchange modulation in zero-field.","marker":"[15]"}],"fun_headline_variants":["GdRu2Si2 magnons need only bilinear exchange","Spin waves in GdRu2Si2 fit with Heisenberg exchange","No higher-order exchange in GdRu2Si2 spin waves","Simplest exchange model explains GdRu2Si2 magnons","Magnon data rules out exotic exchange in GdRu2Si2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The null result for anisotropic and higher-order exchange relies on the measurement being sensitive enough to reveal such terms; the data are low signal-to-noise, median-filtered, and the dispersion is not resolved below about 0.8 meV, so terms below roughly 1 micro-eV would go undetected.","fun_headline_variants_meta":{"raw":{"variants":["GdRu2Si2 magnons need only bilinear exchange","Spin waves in GdRu2Si2 fit with Heisenberg exchange","No higher-order exchange in GdRu2Si2 spin waves","Simplest exchange model explains GdRu2Si2 magnons","Magnon data rules out exotic exchange in GdRu2Si2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00017,"raw_usage":{"total_tokens":1256,"prompt_tokens":923,"completion_tokens":333,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":244}},"tokens_in":539,"tokens_out":333,"duration_ms":3224,"temperature":1.0,"reasoning_tokens":244,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:33:11.011918+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the magnon dispersion below 0.8 meV—for example with a cold-neutron triple-axis spectrometer as the authors propose—would settle the claim: if the bottom of the band deviates from the $J_1$-$J_8$ bilinear prediction, or if a fit with a biquadratic exchange term significantly outperforms the bilinear model, the central null claim is refuted. A simpler check on existing data is to add an $(\\mathbf{S}_i\\cdot\\mathbf{S}_j)^2$ term to the fit and test whether the residuals drop meaningfully.","supporting_citations":[{"cited_title":"Bouaziz, E","cited_arxiv_id":null,"evidence_quote":"Supplies the ab initio RKKY exchange parameters and Fermi-surface mechanism that the fitted model is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Determines the zero-field double-Q ground state and the magnetic propagation vector used to assess the model's minima."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Zero-field powder spin dynamics study that corroborates the ab initio model and the absence of higher-order exchange."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the multi-Q phase zoology and field-direction anisotropy that motivate the search for anisotropic terms."},{"cited_title":"Luttinger and L","cited_arxiv_id":null,"evidence_quote":"Supplies the theorem that bilinear Heisenberg exchange alone stabilizes single-Q order, motivating the need for additional terms in multi-Q states."},{"cited_title":"Lyons and T","cited_arxiv_id":null,"evidence_quote":"Extends the theorem on ground-state spin configurations relevant to the implication about multi-Q stability."},{"cited_title":"Matsuyama, T","cited_arxiv_id":null,"evidence_quote":"Quantum oscillations evidence for Fermi-surface reconstruction between field-polarized and multi-Q states, supporting the field-dependent RKKY proposal."},{"cited_title":"Yasui, C","cited_arxiv_id":null,"evidence_quote":"Shows coupling between itinerant electrons and local moments and charge density waves, used to argue for exchange modulation in zero-field."}],"review_version":1}