{"id":"82bb425d-aa8d-4640-999f-2f548a52a29b","arxiv_id":"1909.00146","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Full Brillouin-zone dispersions of two- and three-magnon Bethe string states are observed in SrCo2V2O8 and match Bethe Ansatz calculations.","lead":"Neutron scattering on the quantum magnet SrCo2V2O8 reveals the full energy-momentum dispersions of two- and three-magnon bound states called Bethe strings, matching Bethe Ansatz predictions. The result provides the first complete experimental map of these long-predicted many-body excitations across the Brillouin zone.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mode identification rests on a string-ansatz DSF computed at N=100 with no error check; exact finite-size comparison would decide whether the assigned 2- and 3-string dispersions are secure.","rationale":"Stress-testing the paper, I found no reason to doubt the quality of the experimental data or the qualitative existence of sharp high-energy features. The reader's weakest assumption captures the main logical hinge: the theoretical spectra are the sole tool for assigning the observed continuum to χ(2)/χ(3), and these spectra are produced by an approximation with no stated error estimate. I agree with the reader's identification. I recommend a conditional verdict rather than unconditional acceptance because the needed error check is cheap and would close the main gap; if the exact-BAE energies coincide with the ideal-string boundaries, the concern disappears. The discrepancies already visible in Table S1 make this check non-rhetorical. I do not see a basis for reject; the experiment appears carefully done and the mode assignment is plausible.","tokens_in":20886,"tokens_out":12269,"duration_ms":205603,"concrete_test":"Solve the full Bethe ansatz equations for N=100 without the ideal-string simplification for the actual string-containing eigenstates used in Fig. 2 at 2m=0.08 and 0.12, obtaining the exact complex rapidities and eigenenergies for the two- and three-string branches, and compare these energies to the string-ansatz boundaries. Repeat at N=200; if the exact-boundary positions change by less than ~0.2 meV (the approximate experimental energy resolution in the string region), the string-ansatz approximation is validated. If the shift exceeds this, the theoretical spectra and hence the mode identification are not reliable.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the observed neutron spectra unambiguously identify the two- and three-string states, but this identification is only as strong as the Bethe-ansatz spectra used to draw the black/red dispersion boundaries and to label every peak. Those spectra are computed in the string ansatz, in which the complex-rapidity deviations δ_j^n are neglected because N is assumed large (SI, 'Bethe ansatz formalism'), and the numerical results use N=100. The SI gives no test of this approximation at the actual magnetizations (2m=0.08–0.60) and fields, and the determinant formulas for string form factors must be regularized, so there are two unvalidated steps before any experimental comparison. Independent support is limited: the theoretical curves come from the authors' own preprint (ref 26), and the displayed first-moment sum-rule saturation (Fig. S6) is a necessary but not sufficient check of the individual line shapes and intensities used for mode labeling. The absolute-energy offsets in Table S1 (e.g. χπ(2): 2.5 vs 3.2 meV at 6 T) show the benchmark is not exact, so the error budget is too tight for an unchecked O(1/N) or string-deviation correction to be ignored. Without a finite-size validation, the 'unambiguous' identification is conditional on an approximation whose numerical accuracy is not demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports inelastic neutron scattering measurements of the quasi-one-dimensional spin-chain compound SrCo2V2O8 in longitudinal magnetic fields up to 25.9 T, covering the quantum critical regime. The authors observe dispersive magnetic excitations at intermediate and high energies, which they assign to two-magnon and three-magnon Bethe-string states, together with psinon-psinon and psinon-antipsinon modes. The experimental dispersions, intensities, and field-dependent slopes are compared mode-by-mode with Bethe ansatz calculations of the dynamical structure factor using the Hamiltonian parameters J = 3.55 meV, Delta = 2, g_parallel = 5.5, and g_perp = 2.79, with account taken of the screw-chain structure factor, zone folding, the Co2+ form factor, polarization factors, and instrumental resolution. The central claim is that the excellent agreement between experiment and theory permits the unambiguous identification and full characterization of the two- and three-string states, completing the experimental verification of many-body Bethe strings.","tokens_in":21212,"tokens_out":3259,"duration_ms":33694,"significance":"If the identification holds, this is the first measurement of the full Brillouin-zone dispersion relations and magnetic-field dependence of two- and three-magnon Bethe-string states, a long-standing prediction of the Bethe ansatz. The paper's strengths include the high quality of the INS data, the broad field range (6-25.9 T), the explicit mode-by-mode comparison of energies, intensities, and slopes dE/dB, and the use of parameter-free predictions in the sense that J, Delta, and the g-factors are fixed by earlier work rather than fitted to the present data. The comparison of experimental and theoretical slopes in Table S1 for five modes is a particularly persuasive falsifiable test, and the first-moment sum-rule saturation checks in Fig. S6 provide a necessary global validation of the computed spectral decomposition. The work is likely to be influential for the study of bound states in integrable and near-integrable spin chains.","major_comments":[{"comment":"The theoretical spectra used to label every observed peak are computed within the string ansatz, in which the string deviations are neglected 'when N is large' (SI), and the numerical results use N = 100. The SI states that the determinant formulas for string form factors 'need to be regularized' but gives no test of the accuracy of this approximation at the magnetizations (2m = 0.08-0.60) and fields studied. Because the 'unambiguous identification' claim rests on these theoretical spectra, a finite-size or string-deviation validation is load-bearing. Please add a convergence check with respect to system size (e.g., N = 200, 400) and/or a comparison against an independent method such as exact diagonalization or MPS/DMRG dynamical structure factors at representative fields, and report how much the string-mode boundaries and intensities shift under this check.","section":"SI, 'Bethe ansatz formalism', paragraph on string states"},{"comment":"The absolute eigenenergies show systematic deviations between experiment and theory, for example chi_pi^(2) at 6 T: 2.5 meV (exp.) vs 3.2 meV (theory), and chi_(pi/2)^(3) at 6 T: 10.9 vs 11.5 meV. The text attributes these differences to interchain couplings, but the SI estimates J_perp_eff/J < 10^-2, i.e., J_perp < 0.04 meV for J = 3.55 meV, which is an order of magnitude too small to explain a 0.7 meV shift. Please provide a quantitative estimate of the expected interchain contribution to the mode energies, or discuss alternative sources of the offset (parameter uncertainties in J, Delta, g; finite-size effects; the neglected string deviations), and state how the claimed 'excellent agreement' is affected by these systematic offsets.","section":"Table S1"},{"comment":"The conclusion that the Bethe strings are 'unambiguously identified' is stronger than the evidence presented, because the identification is made by overlaying theoretical boundaries and intensities from a single calculation (the authors' own preprint, Ref. 26) and no alternative assignment is considered or excluded. A decisive cross-check would be to compare not only the peak positions but also the relative integrated intensities of the two- and three-string modes as a function of field, or to compute the DSF with an independent method at one or two fields and show that the string features survive without the string-ansatz approximation. Without such a test, the wording 'unambiguous' should be softened or additional evidence provided.","section":"Main text, last paragraph of Results and Discussion"}],"minor_comments":[{"comment":"The caption contains a typo, 'Bethe Ansalz', which should read 'Bethe Ansatz'.","section":"Figure 3 caption"},{"comment":"Reference 14 is listed as 'Supplimentary information'; the correct spelling is 'Supplementary information'.","section":"Reference 14"},{"comment":"The text claims agreement in linewidth, but no quantitative linewidth fits or extracted widths are shown; consider either presenting fitted linewidths or phrasing the statement in terms of the observed line shapes shown in Fig. 3.","section":"Main text, paragraph after Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The experimental work is impressive and the comparison with Bethe ansatz calculations is extensive. The main risk is the unvalidated string-ansatz approximation at N = 100 and the fact that the theoretical spectra used for mode labeling come from the authors' own preprint; the requested finite-size/systematic cross-check would substantially strengthen the 'unambiguous identification' claim. I do not see this as a reason to reject, but it is a load-bearing point that needs to be addressed in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the full Brillouin-zone, field-resolved neutron scattering map of the 2- and 3-string states in SrCo2V2O8, including first observations of the q=pi/2 two-string and the longitudinal R_PAP modes that THz cannot access. The zone-folding analysis is careful, and the mode-by-mode comparison with Bethe ansatz over energies, slopes, and linewidths is unusually detailed. I think the central claim is very likely correct: these are the two- and three-magnon Bethe string dispersions, measured over the whole zone.\n\nCredit where due: this is a hard experiment. The data come from three different instruments with different resolution and field capabilities, and the authors handle that honestly—for example, they state when intensities are not comparable. Table S1 gives experimental and theoretical eigenenergies for five modes across six fields, and the slopes dE/dB match within error bars for most entries. The SI provides the determinant formulas, the sum-rule checks, and the full theoretical formalism, even though no code or raw data are deposited. The first-moment sum-rule saturation in Fig. S6 is a necessary check and it passes.\n\nThe soft spot is the theoretical benchmark. The DSFs are computed in the string ansatz, in which the deviations of complex rapidities from ideal string values are neglected, and the numerics use N=100 with no test of convergence or of that approximation at the actual magnetizations (2m=0.08-0.60). The determinant formulas for string form factors need regularization, which is mentioned but not validated. That matters because every peak label and every dispersion boundary is drawn from those calculated spectra. The authors' own preprint (ref 26) is the source, so independent confirmation is limited. The absolute-energy offsets in Table S1—up to about 1 meV—are attributed to interchain couplings, but a finite-size or string-deviation correction could also produce shifts of that size. I would not call the identification 'unambiguous' in the strictest sense; I would call it well-supported but conditional on the string ansatz at finite N.\n\nThat said, this is not a fatal flaw. The theory was not fitted to the present data, the parameters come from earlier work, and the predicted field dependences and dispersion shapes are genuinely matched. A referee with Bethe-ansatz expertise can push on the N-dependence, and the authors should add a finite-size comparison or exact diagonalization for small systems before publication.\n\nVerdict: this paper deserves serious peer review. It is a significant experimental advance and the main caveat is addressable. I would bring it to the reading group and would cite the experimental result in my own work.","headline":"Strong experimental case for full-BZ Bethe string dispersions in SrCo2V2O8, with the mode identification resting on an N=100 string-ansatz calculation that lacks a finite-size check.","tokens_in":21772,"tokens_out":3000,"would_cite":true,"duration_ms":29566,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q05","82B23"],"pacs":["75.10.Jm","75.40.Gb","78.70.Nx"],"model":"deepseek-v4-flash","headline":"Neutron scattering reveals the full Brillouin-zone dispersions of two- and three-magnon Bethe strings in SrCo2V2O8.","keywords":["Bethe strings","XXZ spin chain","inelastic neutron scattering","quantum bound states","Brillouin zone","SrCo2V2O8","magnetic field dependence"],"falsifier":"A measurement at higher energy resolution of the 3-string continuum at q = π/2 and q = 0, comparing the observed line shape with the Bethe Ansatz prediction: if the line shape or its field dependence deviates from the string-ansatz calculation while the psinon modes still match, the string assignment would be called into question. Alternatively, a direct comparison of measured field-dependent intensities with the determinant-formula predictions in the D^+− channel at multiple wavevectors, beyond the q-points already reported, would test whether the string assignment holds.","tokens_in":20714,"feed_emoji":"🧲","tokens_out":8441,"duration_ms":70878,"temperature":0.7,"pith_summary":"This paper reports inelastic neutron scattering experiments on the one-dimensional spin-chain compound SrCo2V2O8 in high longitudinal magnetic fields, and claims to observe the dispersion relations of two- and three-magnon Bethe string states across the full Brillouin zone, together with their magnetic-field dependence. The measured energies, scattering intensities, and linewidths agree with Bethe Ansatz calculations of the XXZ model for the material. If correct, this completes the experimental verification of the long-predicted many-body Bethe strings and shows that these bound states carry a measurable part of the spin dynamics of one-dimensional quantum magnets, not merely a theoretical curiosity.","feed_headline":"Neutrons reveal the full dispersion of Bethe strings","feed_subtitle":"Inelastic scattering on SrCo2V2O8 traces two- and three-magnon bound states over the whole Brillouin zone, matching Bethe Ansatz theory.","key_machinery":"The central object is the Bethe string: a many-body bound state of magnons in an integrable spin-1/2 XXZ chain, whose Bethe rapidities form complex-valued patterns of length n. The calculations use the algebraic Bethe Ansatz with determinant formulas for the dynamical spin structure factors, including the screw-chain structure factor, Co2+ form factor, anisotropic g-factors, and polarization factors for comparison with the neutron data. The measured spectra are simulated by zone-folding four copies of the single-chain spectrum, and the theoretical intensities are superimposed on the data. The string ansatz neglects string deviations, which is valid when the system size is large; here the numerical system size is N = 100.","core_discovery":"The central claim is that inelastic neutron scattering on SrCo2V2O8, a spin-1/2 Heisenberg-Ising (XXZ) antiferromagnetic chain with easy-axis anisotropy, reveals the complete dispersion relations of two-string and three-string Bethe states over the full Brillouin zone, along with their field dependences up to 25.9 T. The observed modes include the 2-string states at q = π and q = π/2, the 3-string state at q = π/2, and the psinon-psinon and psinon-antipsinon modes; the comparison with Bethe Ansatz calculations of energies, intensities, and lineshapes is reported as excellent. The authors conclude that the many-body Bethe strings are unambiguously identified and fully characterised.","pith_inferences":["A testable extension would be to measure the momentum-resolved linewidths of the string states at higher resolution and compare them with the Bethe Ansatz lineshapes, since the present comparison is concentrated at selected wavevectors.","The small systematic shifts between experimental and theoretical mode energies, attributed by the authors to interchain couplings, could be modelled quantitatively to extract the effective interchain exchange J⊥ rather than treating it as a residual discrepancy.","Similar string-state dispersions should appear in related XXZ chain compounds with different Ising anisotropy Δ; mapping how the two- and three-string spectra evolve with Δ would test the generality of the string-ansatz description.","If string states contribute to the thermal conductivity or spin transport in the quantum critical regime, the field-dependent spectra reported here could inform calculations of dynamical transport in 1D integrable magnets."],"forward_implications":["If the identification is correct, neutron scattering can be used as a direct probe of string-state dispersions in other quasi-one-dimensional integrable magnets, extending beyond the zero-momentum information accessible by terahertz spectroscopy.","The measured field-dependent slopes dE/dB of the string and psinon modes provide a quantitative test of the XXZ model parameters J = 3.55 meV, Δ = 2, and the anisotropic g-factors.","The observation of the 3-string state, though weak in intensity and spread over all wavevectors, demonstrates that even high-order bound states contribute visibly to the dynamical spin structure factor in the field-induced critical regime.","The agreement justifies using the Bethe Ansatz as a predictive tool for the spin dynamics of SrCo2V2O8, including the relative intensities of the transverse and longitudinal response channels.","The full Brillouin-zone mapping, enabled by the fourfold screw-chain zone folding, establishes SrCo2V2O8 as a benchmark material for studying quantum many-body bound states beyond the lowest-energy excitations."],"supporting_citations":[{"why":"Bethe's original 1931 prediction of two-magnon bound states in the one-dimensional Heisenberg chain, which the present work claims to verify.","marker":"[1]"},{"why":"The earlier terahertz spectroscopy study that observed string states only at the zone centre; the present work extends the observation to the full Brillouin zone.","marker":"[13]"},{"why":"The Bethe Ansatz calculation of quantum spin dynamics of the axial XXZ chain, supplying the theoretical spectra used for comparison.","marker":"[26]"},{"why":"Yang and Yang's proof of Bethe's hypothesis for the ground state of the anisotropic chain, underpinning the Bethe Ansatz formalism used.","marker":"[28]"},{"why":"Korepin's determinant formula for norms of Bethe wavefunctions, used in the form-factor calculations of the dynamical structure factors.","marker":"[38]"},{"why":"Slavnov's calculation of scalar products and form factors in the algebraic Bethe Ansatz framework, used for the intensity calculations.","marker":"[39]"},{"why":"Caux, Hagemans, and Maillet's computation of dynamical correlation functions for Heisenberg chains, supplying the method for evaluating the dynamical structure factors of string states.","marker":"[48]"},{"why":"The supplementary information of this paper, providing the crystal structure, zone-folding, and parameter details needed for the data-theory comparison.","marker":"[14]"},{"why":"Earlier study of spinon confinement in SrCo2V2O8 that determined J = 3.55 meV and Δ = 2, the model parameters used in the Bethe Ansatz calculations.","marker":"[16]"}],"fun_headline_variants":["Full dispersion of Bethe strings mapped in neutron experiment","Bethe strings' complete dispersion revealed by neutrons","Neutron scattering pins down Bethe string dispersions","SrCo2V2O8 neutron data expose full Bethe string dispersion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole identification rests on the assumption that the complex mathematical solutions for the bound states can be approximated as perfectly rigid 'strings'; if that approximation is not accurate for the chain of 100 spins used in the calculations, the predicted dispersions and intensities could be wrong.","fun_headline_variants_meta":{"raw":{"variants":["Full dispersion of Bethe strings mapped in neutron experiment","Bethe strings' complete dispersion revealed by neutrons","Neutron scattering pins down Bethe string dispersions","SrCo2V2O8 neutron data expose full Bethe string dispersion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000174,"raw_usage":{"total_tokens":1239,"prompt_tokens":860,"completion_tokens":379,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":311}},"tokens_in":476,"tokens_out":379,"duration_ms":4090,"temperature":1.0,"reasoning_tokens":311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T06:00:20.120266+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement at higher energy resolution of the 3-string continuum at q = π/2 and q = 0, comparing the observed line shape with the Bethe Ansatz prediction: if the line shape or its field dependence deviates from the string-ansatz calculation while the psinon modes still match, the string assignment would be called into question. Alternatively, a direct comparison of measured field-dependent intensities with the determinant-formula predictions in the D^+− channel at multiple wavevectors, beyond the q-points already reported, would test whether the string assignment holds.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Bethe's original 1931 prediction of two-magnon bound states in the one-dimensional Heisenberg chain, which the present work claims to verify."}],"review_version":1}