{"id":"190be456-2a75-4535-873b-6b98dda85a6a","arxiv_id":"2608.08587","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":16,"one_line_summary":"Co3V2O8 shows sharp terahertz modes with field slopes two to four times the single-magnon slopes, interpreted as interacting multi-magnon quasiparticles stabilized by strong exchange anisotropy.","lead":"A terahertz spectroscopy study of the kagome-staircase magnet Co3V2O8 maps one-magnon modes and finds sharp high-energy branches whose field slopes are two to four times the one-magnon slopes. These branches are interpreted as interacting multi-magnon quasiparticles in a three-dimensional magnet, an unusual setting for such composite excitations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The integer magnon-number labels in Fig. 6 rest on an unverified field-slope-to-number mapping, so the central multi-magnon-quasiparticle claim lacks quantitative support.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the identification of multi-magnon quasiparticles relies on an unstated mapping from field-slope ratios to magnon number. My read agrees with that assessment. The experimental fingerprints—sharpness, polarization selectivity, and avoided crossings—are real, but they do not by themselves establish that the observed branches are interacting multi-magnon bound states. The fitted one-magnon model provides a useful noninteracting benchmark, but the multi-magnon sector is only treated through unquantified qualitative statements and noninteracting continua. The missing quantitative step is a calculation of the interacting multi-magnon spectrum from the fitted Hamiltonian. Because the paper explicitly states that a quantitative determination of binding energies and spectral weights requires calculations beyond linear spin-wave theory, this is a stated limitation of the manuscript itself. The CONDITIONAL verdict is appropriate: the central claim is plausible but not yet supported. I do not recommend rejection, because the proposed test could reasonably confirm the interpretation; I also do not recommend acceptance, because the key assignment is currently unjustified. No adjustment to the reader's verdict is needed.","tokens_in":15805,"tokens_out":7654,"duration_ms":93386,"concrete_test":"Solve the multi-magnon problem for the Hamiltonian fitted in Table I: for example, exact diagonalization of a finite cluster in the S^x sectors with the fitted exchanges (including the symmetry-allowed Jxy and Jxz terms, varied within the one-magnon fit constraints), or a two-magnon Bethe-Salpeter equation on the full lattice. Extract the field-dependent bound-state energies and the corresponding TDTS absorption, and overlay the predicted slopes on Fig. 6(c). If the lowest sharp computed branches occur at slopes of 2, 3, and 4 times the relevant one-magnon slopes and reproduce the 2m-1/3m-4 anticrossing, the assignment is supported; if they do not, the magnon-number labels and the interaction claim are unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central assertion—that the sharp branches above 1 THz are interacting n-magnon quasiparticles—rests on assigning integer n from ratios of field slopes (Sec. V.A, Fig. 6). This assignment is not derived. In the fitted model (Table I), the two Co sites have different effective g-factors, g_x^s = 4.80 and g_x^c = 3.41, and the measured one-magnon slopes differ by about 40%; the Zeeman slope of a multi-magnon state depends on which sites the spin flips occupy, and is not necessarily 2, 3, or 4 times any single observed one-magnon slope. Moreover, an n-magnon continuum edge or van Hove singularity also shifts with the sum of the constituent g-factors, so an integer slope ratio does not by itself imply a bound state. The paper's non-interacting continua (Fig. 5) are described only qualitatively as broad; no side-by-side quantitative comparison of these calculated continua with the measured spectra is given. The claimed even-odd hybridization (Sec. V.B) additionally invokes Jxy and Jxz terms that are absent from the fitted model in Eq. (2) and are unconstrained by the fit. Until the interacting multi-magnon spectrum of the fitted Hamiltonian is actually computed, the labels 2m/3m/4m and the conclusion that these branches are interaction-stabilized quasiparticles are not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports time-domain terahertz spectroscopy of the kagome-staircase antiferromagnet Co3V2O8. The authors identify low-energy one-magnon modes by their magnetic-dipole selection rules and follow their field and temperature evolution. Combining their TDTS zone-center data with previously published inelastic-neutron-scattering dispersions, they fit an effective spin-1/2 anisotropic-exchange Hamiltonian and show that it reproduces the one-magnon spectrum. Using this model as a noninteracting linear-spin-wave benchmark, they calculate multi-magnon continua. They then interpret sharp high-energy branches, with field slopes reported as two to four times the one-magnon slopes, as interacting two-, three-, and four-magnon quasiparticles, and interpret apparent avoided crossings as hybridization between even- and odd-magnon sectors. The central claim is that strong exchange anisotropy stabilizes multi-magnon quasiparticles in a three-dimensional magnet.","tokens_in":16229,"tokens_out":7966,"duration_ms":88591,"significance":"If established, this result would be significant: it would extend the small family of materials with spectroscopically identified multi-magnon quasiparticles to a three-dimensional anisotropic magnet, and it would demonstrate that TDTS can access these interaction-generated modes. The one-magnon part of the work is a genuine strength: the fitting strategy combines zone-center field dependence with published dispersions, reports parameter uncertainties, uses publicly available codes, and provides a parameter-free noninteracting benchmark once the Hamiltonian is fixed. The multi-magnon interpretation, however, is currently not supported to the same standard. The magnon-number assignments rest on an unproven slope-to-number mapping, the comparison with the calculated continua is qualitative, and there is an internal inconsistency between the stated parity selection rules and the observed polarizations of the labeled branches. For these reasons the central claim requires major additional quantitative work before it can be accepted.","major_comments":[{"comment":"The assignment of the high-field branches as 2m, 3m, and 4m quasiparticles rests entirely on the statement that their field slopes are two to four times the fundamental one-magnon slopes. This slope-to-number mapping is not derived. The fitted model has site-dependent g-factors, g_x^s = 4.80 and g_x^c = 3.41 (Table I), and the one-magnon modes themselves have different slopes (Fig. 2(b)). The Zeeman slope of an n-magnon state is the sum of the g-factors of the flipped sites, so the ratio of a multi-magnon slope to any particular one-magnon slope is not generally equal to n. Moreover, the edge of a noninteracting n-magnon continuum has the same slope property, so an approximately integer slope ratio does not by itself distinguish a bound state from a continuum edge. The labels should be justified by computing the field dependence of candidate n-magnon states in the fitted Hamiltonian, or at least by fitting the observed slopes with the site-dependent g-factors and showing that the inferred magnon numbers are unique.","section":"Sec. V.A, Fig. 6, Table I"},{"comment":"The central distinction between sharp quasiparticle branches and broad noninteracting continua is made by qualitative inspection. The calculated spectra in Fig. 5 are plotted on a logarithmic intensity scale and over a different field axis from the data in Fig. 6, with no overlay, no common frequency cut at fixed field, and no quantitative measure of linewidth or integrated weight. Without a quantitative comparison, the statement that the measured branches 'depart from the calculated multi-magnon continua' is not established. Please provide, for example, cuts of the measured absorption at selected fields with the calculated alpha_yy and alpha_zz superimposed, and specify a quantitative criterion for sharpness relative to the noninteracting continuum.","section":"Sec. V.A, Figs. 5 and 6"},{"comment":"There is an internal inconsistency in the magnon-number selection rules. The text states that the longitudinal response (h parallel to the ordered moment, along a) couples to even-magnon excitations, whereas transverse responses couple to odd-magnon excitations, and Fig. 5 indeed shows two-magnon continua only in alpha_xx and three-magnon continua in alpha_yy and alpha_zz. However, the experimental branches labeled 2m-1, 2m-2, and 4m-1 in Fig. 6 are measured with h parallel to b, which is a transverse geometry. Under the parity argument given in Sec. V.A, these even-magnon branches should be absent from h||b data. The later introduction of Jxy and Jxz terms to explain even-odd hybridization (Sec. V.B) contradicts the parity-conserving model used for the noninteracting benchmark. The authors must show, within one consistent Hamiltonian, how even-magnon spectral weight appears in the h||b response, and how the benchmark selection rules are modified.","section":"Sec. V.A and Sec. V.B, Fig. 6"},{"comment":"The conclusion that the sharp high-energy branches are interaction-stabilized multi-magnon quasiparticles is not supported by any calculation of the interacting multi-magnon spectrum of the fitted Hamiltonian. The manuscript acknowledges that 'a quantitative determination of their binding energies and spectral weights requires calculations beyond the LSWT scheme.' Since the paper's central claim is precisely that these modes are interacting quasiparticles, a calculation demonstrating that the fitted anisotropic exchange produces bound states at the observed energies and field slopes is needed. In the absence of such a calculation, the abstract and conclusion should be tempered, and the branches should be described as candidate multi-magnon excitations rather than established interacting quasiparticles.","section":"Sec. V.A and Sec. V.B"}],"minor_comments":[{"comment":"The statement that the sharp high-energy excitations are 'particularly pronounced in the h||b polarization configurations' cites Ref. [8] (Sala et al., Nat. Commun. 2021, on a honeycomb-lattice Van Hove singularity); this reference does not appear to support the claim. Please check and correct the citation.","section":"Sec. V.A, citation [8]"},{"comment":"The caption says the calculated INS spectra 'can be directly compared to the experimental report in Ref. [32]', but no side-by-side comparison or residual plot is shown. Displaying the experimental data alongside the calculated spectra would substantially strengthen the one-magnon validation.","section":"Fig. 4"},{"comment":"The Jxy and Jxz exchange terms are said to be 'set to zero' in the fits, but are later invoked as the mechanism for even-odd hybridization. Please clarify whether these terms are zero for symmetry reasons or merely unconstrained by the one-magnon data, and distinguish these two cases explicitly in the text.","section":"Eq. (2) and Sec. V.B"},{"comment":"The abstract and conclusion state that multi-magnon quasiparticles 'can be stabilized' and are 'identified' in CVO, while Sec. V.A notes that binding energies and spectral weights have not been calculated. The wording should be reconciled with the evidence presented in the paper.","section":"Sec. VI and Abstract"},{"comment":"There are several typographical errors: 'transition form the ferromagnetic phase' should be 'from'; 'the spectra weight exhibits' should be 'spectral weight'; 'Figures 5 present' should be 'Figure 5 presents'; and 'the high energy multi-magnon response' appears in Fig. 1(e) caption as an unspaced phrase.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The experimental dataset is valuable and the one-magnon modeling is a solid contribution. However, the multi-magnon interpretation currently overreaches: the slope-based magnon-number assignments are unjustified, the comparison to the noninteracting benchmark is qualitative, and the h||b observation of even-magnon-labeled branches conflicts with the paper's own parity selection rules. I would recommend requiring the authors to address these points quantitatively before publication, and to temper the abstract if the interacting calculation cannot be provided."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper delivers a genuinely new experimental dataset and a credible one-magnon Hamiltonian. The TDTS measurements resolve polarization-dependent modes, and the fit combining their field evolution with the prior INS dispersion produces an anisotropic-exchange model that plausibly describes the one-magnon sector. The noninteracting continua from that model are a legitimate, parameter-free benchmark.\n\nThe central multi-magnon claim, though, rests on an assumption the authors never justify. They label branches 2m, 3m, 4m from field slopes two to four times those of single magnons. That mapping fails to account for the site-dependent g-factors (4.80 vs 3.41), and an n-magnon continuum edge shifts with the same Zeeman sum. So integer slope ratios alone don't identify bound states. The paper also gives no quantitative comparison between the calculated continua and the measured spectra; 'broad and featureless' is asserted, not shown. The anticrossing in Fig. 6 is a real piece of data, but the explanation requires Jxy and Jxz exchanges that the fit deliberately sets to zero. The authors concede that a quantitative calculation is needed, which is honest, but it's also the missing piece that would make the claim.\n\nIf they computed the interacting multi-magnon spectrum of their own fitted model—two-magnon T-matrix or ED on a cluster—and showed branches below the continuum with Zeeman slopes from the fitted g-factors, I'd be convinced. As it stands, the paper is a strong experimental report with an overreaching interpretation. Citation pattern is fine; they engage the relevant CVO and bound-state literature.\n\nRecommendation: send it to peer review. It deserves referee time, and the right referee will ask for that calculation or at least for a slope analysis that uses the actual g-tensors. With that, the paper could be good; without it, the headline claim stays conditional.\n\nYours,\n[your name]","headline":"New THz data and a plausible one-magnon Hamiltonian, but the multi-magnon quasiparticle claim rests on an unverified slope-to-number mapping and needs an interacting calculation.","tokens_in":16747,"tokens_out":4368,"would_cite":true,"duration_ms":47832,"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":"Sharp terahertz branches in the kagome-staircase magnet Co3V2O8 are interacting multi-magnon quasiparticles, and a strongly anisotropic spin-1/2 model captures the one-magnon spectrum.","keywords":["terahertz spectroscopy","multi-magnon quasiparticles","magnon bound states","kagome staircase","Co3V2O8","anisotropic exchange","spin-1/2 model","spin dynamics"],"falsifier":"A beyond-linear-spin-wave calculation of the two- and three-magnon spectral functions using the fitted exchange matrices would settle the claim: if the sharp high-energy branches do not emerge as bound states below the noninteracting continua, or an exact calculation places them inside the continuum, the quasiparticle interpretation fails. A simpler check is to measure the field slopes of the labeled 2m and 3m branches with the field along a direction other than a and verify that the slope ratios track the appropriate site-resolved g-factors rather than simple integer multiples of the one-magnon slope.","tokens_in":15571,"feed_emoji":"🧲","tokens_out":7786,"duration_ms":74439,"temperature":0.7,"pith_summary":"The paper uses polarization-resolved, time-domain terahertz spectroscopy to map the magnetic excitations of Co3V2O8, a three-dimensional magnet whose cobalt ions form kagome-staircase layers. It argues that the low-energy modes are single magnons, elementary spin waves, obeying the magnetic-dipole selection rule, and that combining their field dependence with published neutron-scattering dispersions yields an effective spin-1/2 Hamiltonian with strongly anisotropic exchange that quantitatively reproduces the one-magnon spectrum. The same model, without interactions, predicts only weak broad multi-magnon continua above about 1 THz, but the experiment finds sharp branches whose field slopes are two to four times the one-magnon slopes, together with anticrossings between branches assigned different magnon numbers. The central claim is that these sharp branches are interacting multi-magnon quasiparticles, bound composites of two, three, or four spin flips, stabilized by strong exchange anisotropy rather than being continuum features. If true, Co3V2O8 is a three-dimensional magnet with well-defined composite magnon quasiparticles, a regime mostly seen before in one- and two-dimensional systems.","feed_headline":"Sharp THz modes in Co3V2O8 are multi-magnon quasiparticles","feed_subtitle":"A fitted spin-1/2 anisotropic-exchange model benchmarks the one-magnon spectrum and exposes bound magnon branches.","key_machinery":"The central object is an effective spin-1/2 Hamiltonian for the two inequivalent Co sites (spine and cross-tie) of the kagome-staircase lattice, with anisotropic exchange matrices J1, J3, J4, J6 and an isotropic J12, fitted to both the terahertz field dependence at the zone center and the neutron-scattering dispersions. The fitted model provides the noninteracting benchmark: it fixes the one-magnon spectrum, gives site-dependent g-factors along the ordered a axis ($g_s^x=4.80$, $g_c^x=3.41$), and produces the weak, broad multi-magnon continua in the longitudinal and transverse susceptibilities that the sharp experimental branches depart from. The interaction mechanism invoked for the bound states is the broken-bond attraction: because J3 and J4 are strongly ferromagnetic along the ordered direction, two spin flips on neighboring sites leave their shared bond ferromagnetically aligned, so the exchange cost is lower than for two separated flips. The parity argument is also load-bearing: the anisotropic exchange breaks continuous spin-rotation symmetry while preserving magnon-number parity, so even-magnon and odd-magnon sectors remain distinct, and the observed even-odd anticrossings point to Jxy and Jxz terms that change magnon number by odd integers but do not affect the harmonic one-magnon spectrum.","core_discovery":"The central discovery claim is that the high-energy terahertz response of Co3V2O8 contains sharp, magnetic-dipole-active branches whose magnetic-field slopes are two, three, and four times those of the elementary one-magnon modes, and that these branches are interacting multi-magnon quasiparticles rather than features of the noninteracting multi-magnon continuum. The support comes from a joint analysis of the zone-center field dependence measured by terahertz spectroscopy and the momentum-resolved one-magnon dispersions measured by inelastic neutron scattering; the resulting anisotropic-exchange spin-1/2 Hamiltonian reproduces the one-magnon spectrum quantitatively, and its linear-spin-wave multi-magnon continua are broad and weak, in contrast to the observed sharp branches. The observed avoided crossing between a two-magnon and a three-magnon branch, and the broadening and spectral-weight transfer at the three-magnon/four-magnon crossing, are read as hybridization among composite magnon states. The paper additionally attributes the interaction to a broken-bond mechanism: a neighboring pair of spin flips sharing a strongly ferromagnetic bond costs less exchange energy than two separated flips, so strong exchange anisotropy stabilizes bound composite magnons in three dimensions.","pith_inferences":["The slope-to-magnon-number assignment could be tested by computing the field dependence of two- and three-magnon bound states in the fitted model; if the site-dependent g-factors mix into the bound-state slopes, the integer-multiple labelling in the field-dependence plots may need revision.","The polarization-selective 3m-1 electromagnon suggests that magnetoelectric coupling acts differently on even- and odd-magnon sectors, so measuring its field dependence in both Faraday and Voigt geometries could separate electric- and magnetic-dipole matrix elements.","If strong exchange anisotropy is the stabilizing ingredient, then substituting Ni or Mg on the Co sites, which changes the spin-orbit-entangled crystal field, should systematically tune the binding energies and give a chemical series for testing the mechanism."],"forward_implications":["If the central claim is right, Co3V2O8 becomes a three-dimensional example where strong exchange anisotropy stabilizes well-defined composite magnon quasiparticles, extending the bound-magnon phenomenology beyond one- and two-dimensional settings.","The fitted anisotropic-exchange Hamiltonian, together with its one-magnon benchmark, gives future interaction-aware calculations a quantitative starting point for computing binding energies and spectral weights of the multi-magnon branches.","The parity argument predicts that even-magnon and odd-magnon excitations appear in different terahertz polarizations, so polarization-resolved spectra can be used to sort high-energy branches by magnon-number parity.","The observed even-odd anticrossings imply that Jxy and Jxz exchange components, though invisible in the one-magnon spectrum, are constrained by the multi-magnon hybridization data; including them should reproduce the avoided-crossing gaps.","Sharp multi-magnon branches with slopes two to four times the one-magnon slope provide a spectroscopic fingerprint that can be searched for in other strongly anisotropic three-dimensional magnets."],"supporting_citations":[{"why":"Supplies the inelastic-neutron-scattering one-magnon dispersions that the anisotropic-exchange Hamiltonian is fitted to.","marker":"[32]"},{"why":"Supplies the magnetic structure, ordering wavevector, and ordered-moment magnitudes used to constrain the ferromagnetic ground state and site g-factors.","marker":"[30]"},{"why":"Provides the magnetic-dipole absorption formula that connects the calculated response functions to the measured terahertz absorption coefficient.","marker":"[7]"},{"why":"Foundational theory of two-spin-wave bound states that motivates interpreting sharp branches as bound magnons rather than continuum features.","marker":"[4]"},{"why":"Provides the broken-bond mechanism by which neighboring spin flips attract in easy-axis magnets, the interaction invoked for the bound states.","marker":"[20]"},{"why":"Supports the possibility of two-magnon bound states in two and three dimensions, strengthening the case for a three-dimensional magnet.","marker":"[51]"},{"why":"Used to simulate the magnetic excitation spectra from the fitted Hamiltonian.","marker":"[42]"},{"why":"Used to calculate the noninteracting multi-magnon continua that serve as the benchmark from which the sharp experimental branches depart.","marker":"[43, 44]"}],"fun_headline_variants":["THz reveals interacting multi-magnon states in Co3V2O8","Sharp THz branches are bound multi-magnon quasiparticles","Co3V2O8 shows multi-magnon quasiparticles with steep slopes","Anisotropic model exposes multi-magnon bound states in Co3V2O8","Multi-magnon quasiparticles emerge in Co3V2O8 THz spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that a branch's field slope tells you its magnon number, twice as steep meaning two magnons and three times meaning three, even though the two cobalt sites have different couplings to the field (g=4.80 and g=3.41), so the slope-to-number mapping is not automatically guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["THz reveals interacting multi-magnon states in Co3V2O8","Sharp THz branches are bound multi-magnon quasiparticles","Co3V2O8 shows multi-magnon quasiparticles with steep slopes","Anisotropic model exposes multi-magnon bound states in Co3V2O8","Multi-magnon quasiparticles emerge in Co3V2O8 THz spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000708,"raw_usage":{"total_tokens":3241,"prompt_tokens":1048,"completion_tokens":2193,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":664,"completion_tokens_details":{"reasoning_tokens":2085}},"tokens_in":664,"tokens_out":2193,"duration_ms":15630,"temperature":1.0,"reasoning_tokens":2085,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:30:56.000981+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A beyond-linear-spin-wave calculation of the two- and three-magnon spectral functions using the fitted exchange matrices would settle the claim: if the sharp high-energy branches do not emerge as bound states below the noninteracting continua, or an exact calculation places them inside the continuum, the quasiparticle interpretation fails. A simpler check is to measure the field slopes of the labeled 2m and 3m branches with the field along a direction other than a and verify that the slope ratios track the appropriate site-resolved g-factors rather than simple integer multiples of the one-magnon slope.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the inelastic-neutron-scattering one-magnon dispersions that the anisotropic-exchange Hamiltonian is fitted to."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the magnetic structure, ordering wavevector, and ordered-moment magnitudes used to constrain the ferromagnetic ground state and site g-factors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the magnetic-dipole absorption formula that connects the calculated response functions to the measured terahertz absorption coefficient."},{"cited_title":"Wortis, Bound States of Two Spin Waves in the Heisenberg Ferromagnet, Phys","cited_arxiv_id":null,"evidence_quote":"Foundational theory of two-spin-wave bound states that motivates interpreting sharp branches as bound magnons rather than continuum features."},{"cited_title":"El Mendili, T","cited_arxiv_id":null,"evidence_quote":"Provides the broken-bond mechanism by which neighboring spin flips attract in easy-axis magnets, the interaction invoked for the bound states."},{"cited_title":"Dong and J.-X","cited_arxiv_id":null,"evidence_quote":"Supports the possibility of two-magnon bound states in two and three dimensions, strengthening the case for a three-dimensional magnet."}],"review_version":1}