{"id":"3063133d-431b-42d2-9d9d-8d8375b35bd0","arxiv_id":"2509.04672","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using CLAS threshold photoproduction data and the JM model, the rho-proton scattering length is extracted as 0.23 +/- 0.03 fm, about four times smaller than the omega-proton value.","lead":"This paper uses Jefferson Lab CLAS data to estimate how strongly a rho meson interacts with a proton, reporting a scattering length of 0.23 fm. The value is about four times smaller than the omega-proton interaction, a difference that could test models of hadron structure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Broad rho width undermines the stable-particle threshold expansion: Eq. (2) omits a q^0 term, so the fitted a and hence |alpha_rho p| = 0.23 fm may not be a scattering length.","rationale":"The reader's weakest_assumption identified the quasi-stable treatment of the rho as the key issue. I agree and sharpen it: the problem is not only the possibility of bound states, but the finite width itself, which removes the clean two-body threshold on which Eq. (2) depends. Because the first data point lies only ~24 MeV above the nominal threshold while Gamma_rho ~ 147 MeV, the odd-power expansion in the nominal q is not justified without a quantitative test. The proposed JM model evaluation at q=0 would directly show whether a q^0 term is present, and the c0 re-fit would check whether the published quasi-data require such a term. The paper currently provides neither, so the central numerical claim is conditional on an untested assumption. This supports the reader's CONDITIONAL verdict rather than changing it.","tokens_in":6877,"tokens_out":7235,"duration_ms":84489,"concrete_test":"Use the JM model to compute the integrated gamma p -> rho p cross section at W = m_p + m_rho (q=0) with the full rho Breit-Wigner width. If this value is non-negligible compared with the fitted a q term at qmin (a*143 MeV/c ~ 5.7 microbarn), then Eq. (2) is missing a q^0 term and the extracted a is not the stable-particle threshold coefficient. As a complementary data-level check, re-fit the 11 CLAS/JM cross sections to sigma_t = c0 + a q + b q^3 and test whether c0 is consistent with zero; if not, the linear coefficient is contaminated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on Eq. (2), sigma_t = a q + b q^3 + c q^5, and the VMD identification of a with the S-wave scattering-length amplitude via Eqs. (3)-(5). This expansion assumes a narrow, quasi-stable rho with a genuine two-body threshold at q=0. But the rho width is 147 MeV, and the CLAS/JM cross sections are extracted from gamma p -> pi+pi-p through a Breit-Wigner propagator. The first fitted point has W=1.737 GeV, only about 24 MeV above m_p+m_rho, much smaller than Gamma_rho; the physical pi+pi-p threshold is near 1.08 GeV. At q=0 the rho-p contribution need not vanish, because the rho spectral function extends below the nominal mass. Equation (2) therefore omits a possible q-independent term, and the fitted a can absorb low-mass rho strength. The paper's caveat 'we treat the rho meson as a quasi-stable object' is the assumption itself, not a justification. Agreement with Wang et al. [9] only shows that the two analyses share the same assumption. The bound-state caveat does not address the width problem.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extracts the rho-proton scattering length from CLAS/JM-model integrated gamma p -> rho p photoproduction cross sections near threshold. The cross sections (Table I) are fitted to the odd-power expansion sigma_t = a q + b q^3 + c q^5, and the linear coefficient a is converted to |alpha_rho p| through a VMD-motivated factor B_V, Eqs. (3)-(5). The central result is |alpha_rho p| = 0.23 +/- 0.03 fm, about a factor of four smaller than the omega-proton value. The paper explicitly assumes that the broad rho can be treated as a quasi-stable particle, that no rho-N bound states exist below q_min = 143 MeV/c, and that VMD applies.","tokens_in":1592,"tokens_out":1630,"duration_ms":58842,"significance":"If the extraction is valid, the result provides a first modern determination of the rho-proton scattering length from CLAS data and a sharp phenomenological contrast with the omega-nucleon scattering length, which is relevant for vector-meson-nucleon interactions and hadron-structure theory. The paper is transparent about its main assumptions and uses an existing amplitude analysis of CLAS data. However, the central number rests on a threshold expansion that is questionable for a resonance of width 147 MeV, and no model-systematic uncertainty is quantified. The comparison with the earlier SAPHIR-based result [9] is useful but does not independently validate the quasi-stable-rho assumption because that analysis uses the same expansion.","major_comments":[{"comment":"The central extraction uses sigma_t = a q + b q^3 + c q^5, which assumes a genuine two-body threshold at q=0 where the cross section vanishes linearly. For the rho, Gamma_rho = 147 MeV and the physical pi+pi-p threshold is near 1.08 GeV, while the first fitted point has W=1.737 GeV, only about 24 MeV above m_p + m_rho. The rho spectral function contributes below the nominal threshold, so sigma_t(q=0) need not vanish. If a q^0 term is present, the fitted coefficient a can absorb low-mass rho strength and is not simply related to the on-shell S-wave scattering length. This is the load-bearing assumption of the paper; the sentence 'we treat the rho-meson as a quasi-stable object similar to omega' is the assumption itself, not a justification. Please quantify this effect, for example by refitting Table I with sigma_t = c0 + a q + b q^3 + c q^5 and reporting the change in |alpha_rho p|, or by","section":"Section III, Eq. (2) and Table I"},{"comment":"The quoted uncertainty of 0.03 fm is derived only from the fit to the JM-model cross sections and from JM-model parameter uncertainties. The cross sections themselves are quasi-data, obtained through a specific meson-baryon model with a particular non-resonant ansatz (exponential t-channel propagator) and a fixed resonance set. No systematic uncertainty is assigned to this model choice, so the error bar likely underestimates the total uncertainty. An independent extraction from the same or related data, or at least an estimate of model-form sensitivity, is needed to support the claimed precision of the central result.","section":"Section II and Table I"},{"comment":"The paper argues that working at q_min = 143 MeV/c avoids the effect of the broad rho width. This is not evident: 143 MeV/c corresponds to W = 1.737 GeV, and the momentum is defined with respect to the central rho mass. A 147 MeV-wide resonance has substantial strength at invariant masses far below m_rho, so the distance in q from the nominal threshold does not by itself suppress off-shell effects. The statement 'we assume that there are no rho-N bound states below the experimental q_min' addresses a different question from the width problem; even without bound states, the spectral function below q_min can contribute. Please provide a quantitative argument or a data-driven test that the fitted linear term is insensitive to the rho width.","section":"Section III, paragraph after Eq. (5)"}],"minor_comments":[{"comment":"'a factor of 4 smaller than the size of the hadron' is imprecise; specify that the comparison is to the omega-nucleon scattering length and to the hadronic size scale.","section":"Abstract and Section IV"},{"comment":"The caption uses 'showed' in all three columns; use 'shows' for a present-tense caption.","section":"Table I caption"},{"comment":"The photon CM momentum k is not defined at which W; near threshold it should be specified, since the relation may depend on the choice.","section":"Eq. (4)"},{"comment":"The fitted parameters b and c are not reported. For reproducibility and for judging the quality of the polynomial truncation, the full fit results and chi^2 should be given.","section":"Section III, Eq. (2)"},{"comment":"The figure caption does not describe the data symbols for the omega points or the meaning of the dash-dotted curves beyond 'polynomial fits'; please make the legend complete.","section":"Fig. 3"},{"comment":"The sentence 'unfortunately, we cannot consider Quarkonium beyond the Upsilon' is informal for a journal article; consider rewording.","section":"Introduction"},{"comment":"Reference [8] is a PhD thesis and workshop proceedings; provide the published SAPHIR cross-section reference where available for easier verification.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is short and the central claim is interesting, but the broad-rho threshold expansion is a serious correctness risk that should be tested quantitatively. The comparison with Wang et al. [9] is not an independent check because it shares the same expansion. I would ask the authors to add a c0-term fit and to state the resulting change in |alpha_rho p|; if the change is small, this would substantially strengthen the paper. The lack of model-systematic uncertainty in the JM extraction is also a barrier to claiming a 0.03 fm precision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the new result is |alpha_rho p| = 0.23 +/- 0.03 fm from CLAS/JM, about a factor of 4 below the omega-proton value. It agrees with the earlier SAPHIR-based extraction, so the genuinely new content is an independent cross-check plus a threshold comparison of rho and omega cross sections. The paper is straightforward about its assumptions, which is good, because the central number leans on them.\n\nWhat's actually new: the integrated gamma p -> rho p cross sections are extracted from the CLAS nine-fold differential data using the JM reaction model, and the polynomial fit to sigma_t(q) is done on those points. The fit and the VMD conversion are standard. The paper does not overclaim: it says it determines the absolute value, not the sign, and doesn't separate partial waves. That's honest.\n\nSoft spots, in order. First, the quasi-stable rho assumption is load-bearing. The stress-test note is right: Eq. (2) starts at q^1, and with Gamma=147 MeV the rho spectral function extends below nominal threshold, so a q^0 piece is not obviously absent. The paper's defense is 'we work far from threshold' at qmin=143 MeV/c, but that's only ~24 MeV above nominal threshold in W and far below the width. The absence of bound-state evidence doesn't address the width-induced q^0 term. That said, the same assumption underlies the prior SAPHIR extraction, so agreement with Wang et al. is a consistency check, not independent validation.\n\nSecond, no model systematic uncertainty is reported for the JM extraction. The cross sections carry fit uncertainties but not a model band from, say, alternative resonance parameters or the non-resonant ansatz. The paper states the ansatz works for tpp'<1 GeV^2, but doesn't say what happens if you vary it. This is a real gap, though not disqualifying—the JM model is published and the data are public.\n\nThird, the central result is algebraically tied to the fit coefficient a via Eqs. (3)-(5); it's a rescaled VMD parameter, not a direct amplitude extraction. The paper is clear about this, but readers should not mistake 0.23 fm for a model-independent observable.\n\nWho is this for? People working on vector meson photoproduction and VMD phenomenology. It's a niche but legitimate result. A serious referee would ask for the model-systematic estimate and a more careful discussion of the width effect, but the paper is coherent and worth engaging. I'd send it out.","headline":"An honest, internally coherent extraction of the rho-proton scattering length from CLAS data; the number is VMD-model-dependent and rests on a quasi-stable rho assumption, but the paper is clear about it and deserves a serious referee.","tokens_in":7721,"tokens_out":1700,"would_cite":true,"duration_ms":15700,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The rho-proton scattering length is about 0.23 fm, roughly a quarter of the omega-proton value.","keywords":["rho meson","scattering length","photoproduction","vector meson dominance","CLAS","threshold cross section","nucleon resonance"],"falsifier":"Find the gamma p -> rho p cross section at center-of-mass momenta below 143 MeV/c: a resonance-like bump or bound-state cusp there would break the smooth a q + b q^3 + c q^5 extrapolation and change the linear coefficient; a precise lattice QCD calculation of the rho-nucleon S-wave scattering length that disagrees with 0.23 fm would also settle the question.","tokens_in":6826,"feed_emoji":"⚛️","tokens_out":7087,"duration_ms":62152,"temperature":0.7,"pith_summary":"This paper extends the authors' vector-meson-nucleon scattering-length studies to the rho meson, which had been set aside because its 147 MeV width makes threshold studies ambiguous. The authors use CLAS photoproduction data on gamma p -> pi+ pi- p, analyzed with the Jefferson Lab–Moscow State (JM) meson-baryon reaction model, to isolate the gamma p -> rho p cross section near threshold, then apply a vector-meson-dominance-inspired expansion of the total cross section in odd powers of the rho momentum. They obtain |alpha_rho p| = 0.23 ± 0.03 fm, in agreement with an earlier SAPHIR-based extraction, and about a factor of four smaller than the omega-proton scattering length (~0.82 fm). If correct, the rho meson scatters far more weakly from the proton than its isoscalar partner, a difference that bears on how vector mesons couple to hadronic matter.","feed_headline":"Rho-proton scattering length: 0.23 fm, one-fourth of omega's","feed_subtitle":"A factor-of-four gap between rho-N and omega-N couplings that hadron-structure models will have to explain.","key_machinery":"The extraction chain is: CLAS nine one-fold differential cross sections are fed into the JM meson-baryon reaction model, which isolates the integrated gamma p -> rho p cross section as a function of the rho center-of-mass momentum q. Near-threshold total cross sections are expanded as sigma_t = a q + b q^3 + c q^5, where the linear coefficient a is set by two independent S-waves. The absolute scattering length is then |alpha_V p| = sqrt(a) * B_V, with B_V = alpha m_V k / (12 pi Gamma(V -> e+e-)), the VMD-motivated kinematic factor built from the fine-structure constant, meson mass, photon momentum, and leptonic width.","core_discovery":"The central claim is that the absolute value of the rho-proton S-wave scattering length is |alpha_rho p| = 0.23 ± 0.03 fm, extracted from the CLAS/JM integrated gamma p -> rho p cross sections. The total cross sections for gamma p -> rho p and gamma p -> omega p are nearly equal below a center-of-mass momentum q = 0.2 GeV/c, but the rho cross section then rises much faster. Fitting the rho cross sections to sigma_t = a q + b q^3 + c q^5 gives a = (3.99 ± 0.52) × 10^-2 microbarns/(MeV/c), and converting the linear coefficient to a scattering length through the VMD kinematic factor yields 0.23 fm. The authors emphasize this is a factor of four smaller than the omega-proton value of roughly 0.8","pith_inferences":["If the 0.23 fm value holds, combined rho/omega analyses may need a dynamical mechanism—such as the rho's large width or its two-pion decay—to break the near-degeneracy of the two vector mesons' couplings to nucleons, since simple VMD universality would not generate a factor-of-four difference.","A direct check would come from measuring gamma p -> pi+ pi- p at lower W, where q < 143 MeV/c; the polynomial extrapolation predicts a smooth continuation, and any resonance-like structure there would shift the linear term and invalidate the quasi-stable-rho treatment.","Lattice QCD calculations of the rho-nucleon S-wave scattering length, which do not require the quasi-stable assumption, could test whether the photoproduction extrapolation is reliable."],"forward_implications":["The rho-proton and omega-proton scattering lengths differ by about a factor of four, so vector-meson-nucleon interactions are not universal across the nonet despite the similar masses of rho and omega.","The CLAS/JM extraction agrees with the older SAPHIR-based value of 0.24 fm, giving two independent photoproduction determinations of the rho-proton scattering length.","The near-threshold equality of rho and omega total cross sections below q = 0.2 GeV/c, followed by a rapid rho rise, means the rho-nucleon interaction is only visible clearly through the higher-momentum part of the data.","Because the extracted length is well below the hadron size, the result is consistent with the assumption of no rho-nucleon bound states near threshold.","The same VMD-based expansion connects photoproduction data to scattering lengths for other vector mesons, extending the systematics from narrow states like omega, phi, J/psi, and Upsilon to a broad state like the rho."],"supporting_citations":[{"why":"CLAS nine one-fold differential pi+pi- p cross sections and the JM-model integrated gamma p -> rho p cross sections; the central dataset.","marker":"[2]"},{"why":"A2/MAMI omega threshold photoproduction data used for the omega scattering-length comparison and for the omega fit coefficient.","marker":"[15]"},{"why":"The authors' prior survey of VN scattering lengths for narrow vector mesons, giving the VMD-based method and systematics.","marker":"[1]"},{"why":"SAPHIR-based extraction of |alpha_rho0 p| = 0.24 fm; the comparison value for the rho result.","marker":"[9]"},{"why":"Study establishing VMD applicability for light vector mesons, justifying the method for rho and omega.","marker":"[5]"},{"why":"PDG values for rho/omega masses, widths, and Gamma(V -> e+e-) used in the kinematic factor.","marker":"[7]"},{"why":"Soding ansatz for the non-resonant t-channel amplitudes in the JM model, used to isolate the rho-p channel.","marker":"[11]"}],"fun_headline_variants":["Rho-proton scattering length: 0.23 fm, a quarter of omega's","CLAS data pin rho-nucleon length at 0.23 fm, fourfold below omega","Rho-N length 0.23 fm: factor-of-four gap with omega-N","Rho-proton length 0.23 fm: four times smaller than omega's"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The result rests on treating the 147 MeV-wide rho as a quasi-stable particle and assuming no rho-nucleon bound states or resonance structures below q = 143 MeV/c, so the linear term of the momentum expansion is the true on-shell rho-proton scattering length.","fun_headline_variants_meta":{"raw":{"variants":["Rho-proton scattering length: 0.23 fm, a quarter of omega's","CLAS data pin rho-nucleon length at 0.23 fm, fourfold below omega","Rho-N length 0.23 fm: factor-of-four gap with omega-N","Rho-proton length 0.23 fm: four times smaller than omega's"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001175,"raw_usage":{"total_tokens":4714,"prompt_tokens":782,"completion_tokens":3932,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":3836}},"tokens_in":526,"tokens_out":3932,"duration_ms":27090,"temperature":1.0,"reasoning_tokens":3836,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T05:56:28.044468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find the gamma p -> rho p cross section at center-of-mass momenta below 143 MeV/c: a resonance-like bump or bound-state cusp there would break the smooth a q + b q^3 + c q^5 extrapolation and change the linear coefficient; a precise lattice QCD calculation of the rho-nucleon S-wave scattering length that disagrees with 0.23 fm would also settle the question.","supporting_citations":[{"cited_title":"First results on nucleon resonance photocouplings from the γp → π+π−p reaction,","cited_arxiv_id":null,"evidence_quote":"CLAS nine one-fold differential pi+pi- p cross sections and the JM-model integrated gamma p -> rho p cross sections; the central dataset."},{"cited_title":"Pho- toproduction of the ω meson on the proton near thresh- old,","cited_arxiv_id":null,"evidence_quote":"A2/MAMI omega threshold photoproduction data used for the omega scattering-length comparison and for the omega fit coefficient."},{"cited_title":"Threshold Upsilon-meson photoproduction at the EIC and EicC,","cited_arxiv_id":null,"evidence_quote":"The authors' prior survey of VN scattering lengths for narrow vector mesons, giving the VMD-based method and systematics."},{"cited_title":"First ex- traction of the proton mass radius and scattering length |αρ0p| from ρ0 photoproduction,","cited_arxiv_id":null,"evidence_quote":"SAPHIR-based extraction of |alpha_rho0 p| = 0.24 fm; the comparison value for the rho result."},{"cited_title":"Vector-meson production and vector meson dominance,","cited_arxiv_id":null,"evidence_quote":"Study establishing VMD applicability for light vector mesons, justifying the method for rho and omega."},{"cited_title":"On the apparent shift of the ρ meson mass in photoproduction,","cited_arxiv_id":null,"evidence_quote":"Soding ansatz for the non-resonant t-channel amplitudes in the JM model, used to isolate the rho-p channel."}],"review_version":1}