{"id":"d715e71c-9562-4189-8f6a-cbb13c628c1f","arxiv_id":"2412.06350","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Spectroscopic modes in BaCo2V2O8 are identified as repulsively bound two- and three-magnon states by matching terahertz measurements to numerical simulations of a Heisenberg-Ising spin chain.","lead":"Terahertz spectroscopy of the antiferromagnet BaCo2V2O8 reveals high-energy excitations that match numerical simulations of repulsively bound two- and three-magnon states in a spin chain. These are the first solid-state analogues of repulsively bound atom pairs seen in optical lattices, and they may become building blocks for magnonic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bound-state label for D/T rests on a self-cited real-space distance threshold rather than on a direct demonstration that the peaks are discrete states outside the multi-magnon continuum; an exact-diagonalization check would settle this.","rationale":"The reader's weakest assumption points mainly at the quantitative model and the newly introduced staggered g-factor; I agree that is a risk. I want to isolate a second, equally load-bearing condition: the theoretical assignment of the D/T modes as repulsively bound states. The tMPS spectra are almost certainly reliable for the stated Hamiltonian, and the seven-mode agreement plus magnetization is strong independent support for the model. However, the 'bound' label is not directly observable and is justified only by the average-distance criterion from Ref. 30, which is threshold-based and not validated in this paper. The experimental observation of peaks at the predicted positions is consistent with bound states, but it is also consistent with two-/three-magnon continuum features that are kinematically enhanced by the staggered field; the paper does not show that D and T are discrete poles. This is the soft spot in the transition from 'the model reproduces the modes' to 'the modes are repulsively bound magnons.' A small exact-diagonalization calculation would settle it without changing the experimental work.","tokens_in":15704,"tokens_out":8813,"duration_ms":100460,"concrete_test":"Perform exact diagonalization of Eq. (1) on finite chains (L=16–20) for the fitted parameters at B=30 T and B=40.3 T, with the same g-tensor couplings used in S(q,ω). For the eigenstates carrying the Dπ, Dπ/2, and Tπ/2 peaks, compute the two- and three-magnon content relative to the polarized state, the distances between the flipped spins (full distribution, not just the mean), and the energy separation from the exact two-/three-magnon continuum computed in the same staggered model. If the peaks are discrete states outside the continuum with sharply peaked short-distance distributions, the bound-state label is confirmed; if they lie inside the continuum or have broad distance distributions, the modes are continuum enhancements and the paper's central claim should be weakened.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Within the theoretical comparison, the experimental peak positions are matched to features of S(q,ω) computed by tMPS for Eq. (1), but the statement that D and T are repulsively bound states is an additional interpretive step. The only evidence given in this paper for that step is the average-distance criterion 'the distance should be small (e.g. <2.5 sites [and] <3.5 sites)' taken from Ref. 30. This criterion is not derived, no distribution or variance is reported, and it is not shown to distinguish a discrete two- or three-magnon bound state from the upper edge of the corresponding continuum in the staggered-field model. A continuum edge can also have enhanced spectral weight and a smaller typical inter-magnon distance, especially under the staggered field that is essential for visibility. Consequently, even if Eq. (1) is the correct effective Hamiltonian and the observed modes are faithfully reproduced, the central claim that these modes are repulsively bound magnons is not settled by the presented comparison; it requires a direct classification of the eigenstates contributing to the D and T peaks.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports terahertz and electron-spin-resonance spectroscopy of the quasi-one-dimensional Ising-like antiferromagnet BaCo2V2O8 in transverse magnetic fields. The authors observe several field-dependent modes below and above the quantum critical field Bc ≈ 40 T and compare their frequencies with the dynamical spin structure factor computed by time-dependent matrix-product-state (tMPS) simulations of the XXZ chain in Eq. (1), including a uniform and a staggered transverse g-factor. On the basis of this comparison, they assign the modes M to single-magnon excitations and the modes D and T to repulsively bound two- and three-magnon states, respectively. They argue that the staggered effective field is essential for making these bound states spectroscopically visible and estimate a lifetime of about 8 ps from the linewidths.","tokens_in":15937,"tokens_out":4229,"duration_ms":46871,"significance":"If the identification is correct, this would be the first observation of repulsively bound magnon composites in a solid-state magnet, extending the physics of repulsively bound pairs from ultracold atoms in optical lattices to condensed matter. The paper has clear strengths: the experiments combine pulsed-field ESR, static-field Fourier-transform THz spectroscopy, and magnetization measurements; the tMPS calculations report convergence parameters (bond dimension, truncation errors, system size); and the dependence of the spectra on the staggered field is systematically explored in Extended Data. The magnetization curve provides an additional constraint on the model parameters. However, the central claim that D and T are repulsively bound states depends on two load-bearing points that need stronger support: the model parameters are partly fitted to the same experimental modes used for the identification, and the bound-state classification rests on an average-distance criterion whose discriminatory power is not demonstrated.","major_comments":[{"comment":"The five parameters J, Delta, g_u^x, g_s^x, and g_z are introduced after stating 'Very good agreement is achieved for all seven detected modes, using the parameters ...' and are also constrained by the magnetization. Because the theoretical frequencies in Fig. 1h are generated with parameters that were adjusted to reproduce the same experimental modes, the agreement is in part postdiction rather than independent prediction. This weakens the claim of unambiguous identification in the section 'Identification of repulsively bound magnons'. Please report the fitting procedure, the residuals per mode, and at least one out-of-sample check (for example, a prediction for the field dependence of a mode not used in the fit, or an independent constraint on g_s^x from the g-tensor analysis in Eq. (2)).","section":"Comparison between experiment and theory (Fig. 1g,h)"},{"comment":"The classification of the D and T modes as repulsively bound two- and three-magnon states is based on the statement that the average distance between flipped spins is small ('< 2.5 sites' or '< 3.5 sites', taken from Ref. 30). This criterion is not derived in the present manuscript, no distribution or variance of the distance is reported, and it is not shown to distinguish a discrete bound state from the upper edge of the corresponding multi-magnon continuum. In the presence of the staggered field, the continuum edge can also carry enhanced spectral weight and involve predominantly nearby flipped spins. Please provide a direct many-body diagnosis, for example an exact-diagonalization spectrum of a finite chain showing discrete bound-state branches outside the continuum with significant weight in S(q, omega), or an analysis of the dominant eigenstates contributing to the D and T peaks; and show explicitly that the states at the continuum edge do not satisfy the same distance criterion.","section":"Close to the critical field; Identification of repulsively bound magnons"},{"comment":"The value g_s^x = 0.66 is a new parameter that was set to zero in previous literature (Extended Data Table 2), and the Methods state that without the staggered field the D and T features would have too little spectral weight to be detected. The observability and hence the experimental identification of the bound states therefore hinges on a fitted parameter whose independent justification is not given. Please provide a microscopic estimate of g_s^x from the reported g-tensor components and tilt angle theta = 5 deg using Eq. (2), or an independent experimental determination, and discuss how the mode assignment would change if g_s^x were substantially smaller.","section":"Methods: Anisotropic g-factors and effective magnetic field; Extended Data Table 2"}],"minor_comments":[{"comment":"The definition of S(q, omega) in Eq. (3) is unusual because the absolute values are taken before summing over the four sites l, which is not the standard total transverse structure factor. Please define the normalization and the relation of Eq. (3) to the usual S^{yy}(q,omega) + S^{zz}(q,omega) more explicitly, and spell out the definitions of G_l^alpha.","section":"Methods, Eq. (3)-(4)"},{"comment":"There are several typos: 'eigentstates' in the introduction, 'correspoding' in the main text, and 'respulsively bound' in the Methods. These should be corrected.","section":"Introduction and Methods"},{"comment":"The caption states that experimental and theoretical linewidths are indicated by bars, but no scale or definition of the linewidth extraction is given. Please specify how the linewidths were obtained and what the bars represent.","section":"Fig. 1h caption"},{"comment":"The values of J and Delta from the cited literature differ by about 10% (e.g., J from 2.58 to 3.07 meV and Delta from 1.89 to 2.17). The statement that the present parameters are 'consistent with previously reported values' should be quantified with the appropriate comparison and error estimates.","section":"Extended Data Table 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the readership if the authors can supply a direct demonstration that the D and T modes are discrete many-body bound states rather than continuum-edge features, and if they are transparent about the role of the fitted staggered g-factor in the identification. The reliance on a self-cited distance criterion from Ref. 30 is not inappropriate per se, but it is the only quantitative diagnostic offered and should be strengthened before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis paper reports the first experimental observation of repulsively bound two- and three-magnon states in a solid, extending a phenomenon previously seen only in cold atoms. That claim is new and, on the evidence here, likely correct. The terahertz and ESR measurements are carefully done, the tMPS calculations are described with convergence parameters, and the fit to seven modes plus the magnetization is genuinely good. The parameter set is mostly consistent with previous work; the one new ingredient, a staggered transverse g-factor g_s^xx = 0.66, is fitted here and is what makes the bound-state features visible in the calculated spectrum. That is a legitimate but important caveat: the 'prediction' is not fully independent, because the same parameters were chosen to match the modes being explained.\n\nThe real soft spot is the classification of the D and T peaks as bound states. The paper uses a real-space average-distance criterion (<2.5 or <3.5 sites) taken from the authors' own prior theory paper (Ref. 30), and it does not re-derive it or show the distribution. The stress-test note is right that a continuum edge can also have enhanced spectral weight and small typical distance, especially with a staggered field. That does not mean the interpretation is wrong—the sharpness of the modes and their separation from the continuum are good evidence—but 'unambiguous identification' is stronger than what is actually shown. A referee should ask for either an exact diagonalization of the relevant sector or a direct computation of the eigenstates contributing to S(q,ω) at the D and T peaks, with the distance distribution reported.\n\nI also would like to see error bars on the experimental frequencies and a sensitivity analysis for g_s^xx. These are requests, not refutations. The paper is well within the range of a serious journal; it deserves peer review and, after those clarifications, probably publication. I would cite it and bring it to a reading group on quantum magnetism.\n\nBest,\n\n[Your name]","headline":"First solid-state repulsively bound magnons: careful experiments, good numerics, but the bound-state label rests on a self-cited distance threshold that deserves scrutiny.","tokens_in":16516,"tokens_out":2839,"would_cite":true,"duration_ms":29309,"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":"Terahertz spectroscopy in BaCo2V2O8 reveals repulsively bound two- and three-magnon states, matching predictions for a Heisenberg-Ising chain.","keywords":["repulsively bound magnons","Heisenberg-Ising chain","XXZ model","BaCo2V2O8","terahertz spectroscopy","quantum criticality","matrix product states","spin-chain dynamics"],"falsifier":"A decisive check would be momentum-resolved inelastic neutron scattering on BaCo2V2O8 at 30 T: the model predicts flat, discrete bound-state bands at $q = \\pi/2$ (for $D_{\\pi/2}$ and $T_{\\pi/2}$) and at $q = \\pi$ (for $D_\\pi$) above the single-magnon continuum, whereas scattering spread at the sum of two or three single-magnon energies would show the modes are not bound. Alternatively, measuring the same high-frequency modes in a closely related Ising-like chain with negligible staggered $g$-factor should make the bound-state spectral weight almost vanish, since the paper's own simulations show the $D$ and $T$ features become indiscernible at $g_s^{xx} = 0$.","tokens_in":15524,"feed_emoji":"🧲","tokens_out":16404,"duration_ms":140759,"temperature":0.7,"pith_summary":"Composite particles are normally held together by attraction, yet this paper reports spectroscopic signatures of magnons bound by repulsive antiferromagnetic exchange in a solid. In BaCo2V2O8 under a strong transverse field below its quantum critical point, high-frequency terahertz absorption modes match numerically exact predictions for a Heisenberg-Ising spin chain, and the paper identifies two modes as repulsively bound two-magnon states and one as a repulsively bound three-magnon state. The bound states sit above the unbound-magnon continuum, owe their visibility to a staggered effective field, and live about eight picoseconds. If the identification holds, repulsively bound composites are no longer confined to engineered optical lattices but exist in bulk quantum magnets, where they could influence spin transport and magnonics.","feed_headline":"Repulsively bound magnons observed in a real solid","feed_subtitle":"Antiferromagnetic repulsion, not attraction, binds magnon pairs and triples in BaCo2V2O8.","key_machinery":"The load-bearing object is the repulsively bound magnon in the spin-$1/2$ Heisenberg-Ising (XXZ) chain, defined by the single-chain Hamiltonian in Eq. (1) with antiferromagnetic $J > 0$ and easy-axis anisotropy $\\Delta > 1$. In a strong transverse field near the polarized regime, two or three flipped spins sitting on adjacent sites cost more energy than separated magnons, so the composite is held together by repulsion and propagates as a unit. The mechanism that makes these fragile states observable is a staggered transverse magnetic field produced by site-dependent $g$-factors ($g_s^{xx} = 0.66$), which splits the single-magnon band and concentrates dynamical spectral weight in the bound-state peaks. The paper classifies a feature as bound by computing the average distance between flipped spins in the contributing states, taking distances below 2.5 sites for pairs and below 3.5 sites for triples as the bound-state criterion.","core_discovery":"The central discovery claimed is that the terahertz modes labeled $D_\\pi$, $D_{\\pi/2}$, and $T_{\\pi/2}$ in BaCo2V2O8 are repulsively bound magnon excitations of a one-dimensional transverse-field Heisenberg-Ising antiferromagnet. $D_\\pi$ and $D_{\\pi/2}$ are said to be two-magnon bound pairs at quasi-momenta $\\pi$ and $\\pi/2$, and $T_{\\pi/2}$ a three-magnon bound state at $\\pi/2$. The identification is made by matching seven measured field-dependent resonance frequencies to the dynamical spin structure factor computed with time-dependent matrix product states for a single chain with $J = 2.82\\,\\mathrm{meV}$, $\\Delta = 1.92$, $g_u^{x} = 3.06$, $g_s^{xx} = 0.66$, and $g_z = 0.21$. A staggered transverse field $g_s^{xx}$, previously neglected in the literature, is the ingredient that gives the bound states enough spectral weight to be detected. The paper emphasizes that these are not merely theoretical features: the resonances are well separated from continua and have a measured lifetime of roughly $8 \\pm 2\\,\\mathrm{ps}$.","pith_inferences":["Extending beyond the paper, a natural search would target four-magnon repulsively bound states at higher frequencies, where the same average-distance criterion should produce a distinct band above the three-magnon state.","Because visibility is controlled by the staggered $g$-factor, applying uniaxial pressure to change the CoO6 octahedron tilt angle in BaCo2V2O8 would offer a direct experimental test of the assignment.","The bound states should contribute to dynamical spin or thermal transport in a characteristic frequency window, so heat-current or spin-current correlation measurements could separate that contribution from single-magnon and continuum channels.","A practical consequence the authors leave implicit is a screening recipe for other materials: strong antiferromagnetic intrachain coupling, weak interchain coupling, low dissipation, and a sizable staggered transverse field."],"forward_implications":["Other quasi-one-dimensional Ising-like antiferromagnets with antiferromagnetic exchange and a staggered transverse $g$-factor should show similar $D$ and $T$ resonances in high-field terahertz spectra.","The staggered field acts as a control parameter: increasing it should enhance the bound-state spectral weight, while reducing it should make the features fade into the continuum.","Because the bound states are well separated in energy and live for about $8\\pm2$ ps, they provide a spectroscopic window on multi-magnon correlations in a paradigmatic quantum spin chain.","Above the quantum critical field the same model predicts the bound-state response becomes negligibly weak, consistent with the modes disappearing from the experimental spectra for $B > B_c$."],"supporting_citations":[{"why":"Supplies the anisotropic $g$-factors and exchange parameters that set the single-chain Hamiltonian for BaCo2V2O8.","marker":"[13]"},{"why":"Provides the magnetic phase diagram and single-crystal details used to place the experiment in the transverse-field regime.","marker":"[14]"},{"why":"Establishes the quantum critical field $B_c = 40$ T and the suppression of interchain order, fixing the field window below $B_c$.","marker":"[15]"},{"why":"Predicts repulsively bound magnons in the XXZ chain with a staggered field and supplies the average-distance criterion used to classify $D$ and $T$.","marker":"[30]"},{"why":"Supplies the real-time density-matrix renormalization group evolution used in the tMPS computation of the dynamical spin structure factor.","marker":"[32]"},{"why":"Supplies the matrix product state formalism and accuracy benchmarks for the dynamical simulations.","marker":"[33]"},{"why":"Describes the terahertz free-electron-laser electron spin resonance technique used for the pulsed high-field absorption measurements.","marker":"[41]"}],"fun_headline_variants":["Repulsive forces bind magnons in BaCo2V2O8","Bound magnon pairs and triples seen in magnet","Repulsively bound magnons: solid-state evidence","Terahertz reveals repulsively bound magnon states","Antiferromagnet hosts repulsively bound magnons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that BaCo2V2O8 is quantitatively described by a single Heisenberg-Ising chain with a staggered transverse $g$-factor $g_s^{xx} = 0.66$, and that without this staggered field the bound-state peaks would be too weak to explain the observed $D$ and $T$ modes; if the true staggered field were much smaller or interchain couplings reshaped the spectrum, those modes could have a different origin.","fun_headline_variants_meta":{"raw":{"variants":["Repulsive forces bind magnons in BaCo2V2O8","Bound magnon pairs and triples seen in magnet","Repulsively bound magnons: solid-state evidence","Terahertz reveals repulsively bound magnon states","Antiferromagnet hosts repulsively bound magnons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000836,"raw_usage":{"total_tokens":3700,"prompt_tokens":1054,"completion_tokens":2646,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":2564}},"tokens_in":670,"tokens_out":2646,"duration_ms":17209,"temperature":1.0,"reasoning_tokens":2564,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:44:30.306643+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be momentum-resolved inelastic neutron scattering on BaCo2V2O8 at 30 T: the model predicts flat, discrete bound-state bands at $q = \\pi/2$ (for $D_{\\pi/2}$ and $T_{\\pi/2}$) and at $q = \\pi$ (for $D_\\pi$) above the single-magnon continuum, whereas scattering spread at the sum of two or three single-magnon energies would show the modes are not bound. Alternatively, measuring the same high-frequency modes in a closely related Ising-like chain with negligible staggered $g$-factor should make the bound-state spectral weight almost vanish, since the paper's own simulations show the $D$ and $T$ features become indiscernible at $g_s^{xx} = 0$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the quantum critical field $B_c = 40$ T and the suppression of interchain order, fixing the field window below $B_c$."},{"cited_title":"& Bernier, J.-S","cited_arxiv_id":null,"evidence_quote":"Predicts repulsively bound magnons in the XXZ chain with a staggered field and supplies the average-distance criterion used to classify $D$ and $T$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the real-time density-matrix renormalization group evolution used in the tMPS computation of the dynamical spin structure factor."},{"cited_title":"The density-matrix renormalization group in the age of matrix product states","cited_arxiv_id":null,"evidence_quote":"Supplies the matrix product state formalism and accuracy benchmarks for the dynamical simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the terahertz free-electron-laser electron spin resonance technique used for the pulsed high-field absorption measurements."}],"review_version":1}