{"id":"1ded3675-1b48-422c-9922-8b6f0cc708ae","arxiv_id":"2608.08769","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":1.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":17,"one_line_summary":"This review argues that final-state hadron interactions fitted to cross-section data explain nearly all near-threshold resonances in e+e- annihilation, making new exotic particles unnecessary.","lead":"This review argues that the near-threshold resonances seen in electron-positron annihilation into hadron pairs are all explained by final-state interactions between the produced hadrons, not by new particles. It shows how simple potential models fitted to data reproduce the observed peaks, dips, and even the psi(3770) and Upsilon(4S) bumps.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim is not identifiable: the successful fits do not discriminate FSI from energy dependence placed in the production vertex, e.g., the ψ(3770) charmonium pole, which the model absorbs into a virtual level.","rationale":"The reader's weakest assumption — energy-independence of the production amplitude — is indeed the load-bearing point, and I agree it is not derived. Eq. (4) and Eq. (27) implement this factorization, while Eq. (28) already introduces an energy-dependent dipole form factor FD(s) with fixed s0 = 1 GeV^2; the model therefore does not assume a strictly constant production vertex but a particular smooth one. The decisive weakness is identifiability: the data cannot tell whether a near-threshold peak originates from the production amplitude (e.g., a charmonium pole) or from the FSI wave function, because the fitted rectangular-well potentials can absorb either. The paper's own treatment of ψ(3770) and Υ(4S) as virtual levels rather than quark-model states demonstrates this degeneracy. Since the central claim is a universal negative ('no new, unknown particles needed'), it can only be established if the alternative placement of energy dependence is ruled out; the paper does not perform such a test. Additional limitations noted in the text — exclusion of Belle Λc data due to contradiction with BESIII (Sec. V) and omission of superseded older experiments (Sec. VI) — reduce the force of the 'vast majority of cases' claim, but the primary concern remains the unidentifiability of the production/FSI split. A concrete refit with an explicit ψ(3770) Breit-Wigner term in the production amplitude would settle this, so the conditional verdict remains appropriate.","tokens_in":35979,"tokens_out":7835,"duration_ms":91565,"concrete_test":"Refit the 367 D(*)Dbar(*) cross-section points with the production constants g_i multiplied by a Breit-Wigner factor for ψ(3770) (M = 3773 MeV, Γ = 27 MeV) in the D0D0bar channel, leaving Γ free if needed, and compare χ²/Ndf and the extracted potentials with the FSI-only fit. If χ²/Ndf remains near 1.18, the data do not distinguish an FSI virtual level from a charmonium-pole production amplitude, and the central claim lacks support. Repeat at Υ(4S) in the B(*)Bbar(*) fit. As an analytical cross-check, replace the constant g in Eq. (27) with an energy-dependent B(E); the cross section becomes σ ∝ k^(2L+1)|B(E)|^2|F_L|^2, and since B(E) can be absorbed into a redefined potential, the production/FSI split is identifiable only if B(E) is fixed by independent data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is the factorization in Eqs. (4) and (27): the production amplitude is taken as g·FD(s) with FD fixed as a dipole and g energy-independent, so all near-threshold structure is forced into the FSI wave function at the origin. This premise is asserted from the scale separation, not derived. The paper itself registers a caveat by inserting FD(s) with s0 = 1 GeV^2 (Eq. (28)), which is strongly energy-dependent in the fitted region; its value is fixed a priori and no sensitivity study is shown. The central claim ('no new, unknown particles') is a universal negative, so it requires the model to be identifiable: the same data must not be reproducible with a different split between the production vertex and FSI. That condition is not met. Both known and candidate resonances can be moved into the bare amplitude: in Sec. VIII, ψ(3770) in e+e- -> D0D0bar is reproduced as a virtual level in the D-meson potential, and in Sec. IX, Υ(4S) is reproduced as a virtual level in the B-meson potential, with no charmonium/bottomonium pole in the model. Since the potentials are fitted to the cross-section data, the reported χ² values (e.g., 397/338 for D(*)Dbar(*)) only measure flexibility, not the absence of new particles. The predictive power of the claim therefore rests on the unverified assumption that the bare amplitude is exactly smooth apart from the chosen dipole; if the true amplitude contains a pole, the fitted potentials would absorb it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a phenomenological framework for final-state interactions (FSI) in hadronic pair production near threshold, based on solving the Schrödinger equation with simple potential wells and expressing production amplitudes through the wave function at the origin. The framework is applied to e+e- -> Lambda Lambda, Lambda_c Lambda_c, pbar p, nbar n, to multimeson processes via NbarN intermediate states, and to D(*)Dbar(*) and B(*)Bbar(*). The authors report good chi2/Ndf values (1.18-1.57) and claim that FSI can explain the observed near-threshold structures without invoking new particles. They also predict bound states of Lambda-Lambda_bar and Lambda_c-Lambda_c_bar and reproduce the psi(3770) and Upsilon(4S) bumps as virtual levels in meson-meson systems.","tokens_in":36602,"tokens_out":5508,"duration_ms":59688,"significance":"The paper is valuable as a coherent, broad phenomenological synthesis. Its strengths include a unified treatment of Coulomb, tensor, isospin-breaking, and multichannel effects; simultaneous fits to cross sections and form-factor ratios; and explicit, falsifiable predictions such as the exclusive charged/neutral B(*)Bbar(*) ratios R_ij. If the central 'no new particles' claim is to be sustained, however, the statistical identifiability of the production vertex versus the FSI wave function must be established; the present manuscript does not yet do so. The reported chi2 values are encouraging but not conclusive.","major_comments":[{"comment":"The load-bearing assumption is that the bare production amplitude is energy-independent up to a fixed dipole form factor FD(s)=1/(1-s/s0)^2 with s0=1 GeV^2. This assumption forces all energy dependence near threshold into the FSI enhancement factor. Since the potential parameters are fitted to the same cross sections, the good chi2/Ndf values (e.g., 376/29 in Sec. IV) only show that the model is flexible enough to absorb the data; they do not discriminate FSI from energy dependence in the production vertex. The manuscript needs a sensitivity study in s0 and, ideally, a test in which an explicit resonance pole is added to the bare amplitude (e.g., a charmonium pole in the D(*)Dbar(*) channel) and compared with the FSI-only fit using a model-selection criterion.","section":"Sec. II A and IV, Eqs. (4), (27), (28)"},{"comment":"The claim that these resonances arise as virtual levels in the D- and B-meson potentials is presented as support for the no-new-particles conclusion. However, the potentials are fitted to the cross-section data containing those very resonances, so the 'reproduction' is not independent evidence against a quarkonium pole. An identifiability check is required: for instance, compare the FSI-only model with a model containing an explicit charmonium/bottomonium pole and an energy-independent (or weakly energy-dependent) production vertex, using an information criterion or a chi2 difference test. Without such a test, the central conclusion overreaches.","section":"Sec. VIII and IX, virtual-level description of psi(3770) and Upsilon(4S)"},{"comment":"The Belle Lambda_c Lambda_c data are excluded because they 'contradict' BESIII data, and 'some results from older experiments' are excluded in the NbarN fits because they 'clearly contradict' newer data. The exclusion criteria are not quantified, and no fit including the excluded points is shown. Since the central claim rests on the quality of the quoted chi2 values, the selection must be documented transparently: list the excluded data sets, give the quantitative conflict criterion, and show the effect on the fit parameters and chi2.","section":"Sec. V and VI, data selection"},{"comment":"No parameter uncertainties or theory bands are reported for any fit; only chi2/Ndf values are quoted (e.g., 376/29, 105.6/89, 397/338, 50.2/32). Without uncertainties on V0, a, Uij, g, etc., it is impossible to judge whether the predicted bound states (e.g., Lambda-Lambda_bar at E0 about -30 MeV and Lambda_c-Lambda_c_bar at sqrt(s) about 4530 MeV) are robust or dominated by the arbitrary choice of rectangular-well parametrization. The authors should provide at least a profile-likelihood or bootstrap estimate for the key parameters and show representative theory bands in the figures.","section":"Throughout Secs. IV-IX, parameter uncertainties and theory bands"},{"comment":"The bound states are presented as 'predictions' although they are computed from potentials fitted to the same data that the model is said to explain. This is not circular in the logical sense, but it is a model extrapolation below threshold rather than an independent prediction. The wording should distinguish parameter-derived consequences from independent predictions, and the claims should be tested against observables not used in the fit, such as pbar p -> Lambda Lambda_bar scattering or line shapes in other final states.","section":"Sec. IV after Eq. (17) and Sec. V after Eq. (40), bound-state predictions"}],"minor_comments":[{"comment":"There is a typo in the text: 'We fixed the parameters 0 equal to 1 GeV2' should be 'We fixed the parameter s0 equal to 1 GeV^2'.","section":"Eq. (28)"},{"comment":"In both captions, sigma(1) and sigma(2) are both described as dashed lines, which makes the plots impossible to read; the line styles in the caption should match those in the figure.","section":"Figs. 4 and 5 captions"},{"comment":"Internal references such as 'Fig. II C' should be replaced with the actual figure numbers.","section":"Sec. II C"},{"comment":"The section refers to Ref. [147] for all details of the NbarN potential and fit; for a self-contained review, at least the potential parametrization and the fitted parameter values should be summarized in a table or in an appendix.","section":"Sec. VI"},{"comment":"The regularization parameter b=10 fm^-1 is said to give results 'practically independent' of b, but no supporting test is shown; a brief scan over b or a statement of the resulting variation would justify this claim.","section":"Sec. V, after Eq. (43)"},{"comment":"The caption cites 'Dong2020' without a corresponding entry in the reference list; the intended reference appears to be [168], which should be cited explicitly.","section":"Fig. 16 and reference list"},{"comment":"The text contains typographical errors such as 'multi-mezonic states' and 'mison'; these should be corrected in a final proofreading pass.","section":"Sec. VII"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is largely a synthesis of the authors' own previous publications, which is appropriate for a review-style paper. My main scientific concern is the identifiability of the universal negative claim 'no new particles are needed'; I would like to see the authors address this with concrete model-selection checks rather than only additional chi2 values. I do not see this as a permanent obstacle, but it is a load-bearing point that requires substantive revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you the short version: this is a review of the authors' own final-state-interaction program, not a new research result. The useful part is the collection: the framework is coherent, the fits cover a lot of channels, and the chi2 values are decent. The soft part is the punchline. They claim FSI explains virtually all cases and that no new particles are needed, but the model is not identifiable against the alternative of energy dependence in the production vertex. The stress-test note is right. The dipole form factor in Eq. (28) with s0 = 1 GeV2 is already strongly energy-dependent in the fitted region, so the bare amplitude is not smooth. With potentials fitted to the same data, a virtual level can mimic a resonance pole, and the predicted bound states below threshold are inferred from potentials fitted to the same cross sections, so they are not independent checks.\n\nCredit where due: the Schrodinger-equation approach with rectangular wells is transparent; the treatment of the Lambda_c Lambda_c process with Coulomb and tensor forces is genuinely nontrivial and describes the cross section and GE/GM ratio reasonably well; and the paper does not hide the fact that some data are excluded (Belle Lambda_c, older pbarp data). As a review it is well organized and will be useful for students coming into the area.\n\nThe soft spots are real but not fatal. No parameter uncertainties or theory bands are given, so you cannot judge how rigid the predictions are. The data exclusions are mentioned but not analyzed in detail. Reinterpreting psi(3770) and Upsilon(4S) as virtual levels in D and B meson channels is a provocative position, but it is not a disproof of quarkonium poles; it is a choice of how to distribute the energy dependence between the production vertex and the final-state interaction. The conclusion goes beyond what the current fits can establish.\n\nWho is this for? Experimentalists and students working on near-threshold hadroproduction who want a compact summary of the Novosibirsk approach. It is not a neutral review; it is an advocacy piece for their own program. Still, it is competent and the fits are worth seeing.\n\nRecommendation: send it to peer review. The identifiability issue should be a major referee concern, and the claim about no new particles should be substantially softened, but the paper deserves a serious referee rather than a desk rejection.","headline":"A readable review of the Novosibirsk FSI program, but the 'no new particles' conclusion is stronger than the fits can support.","tokens_in":37089,"tokens_out":2411,"would_cite":false,"duration_ms":30091,"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":"Final-state interactions between slow hadrons, not new exotic particles, explain the near-threshold peaks and dips seen in $e^+e^-$ pair-production cross sections.","keywords":["final-state interactions","near-threshold resonances","e+e- annihilation","hadron pair production","Schrödinger equation","coupled channels","electromagnetic form factors","exotic hadron candidates"],"falsifier":"A high-statistics measurement in a channel where the predicted FSI cross section from the fitted potentials is smooth, such as a narrowly peaked structure in $e^+e^-\\to \\Lambda\\bar\\Lambda$ or $e^+e^-\\to B\\bar B$ that would require an energy-dependent short-distance amplitude or an additional pole, would refute the claim that all near-threshold features are final-state effects. Concretely, the model predicts $\\Lambda\\bar\\Lambda$ and $\\Lambda_c\\bar\\Lambda_c$ bound states near $-30$ MeV and $-40$ MeV below their thresholds; if high-precision data exclude narrow structures at those energies, the bound-state predictions fail.","tokens_in":35796,"feed_emoji":"⚛️","tokens_out":8763,"duration_ms":89320,"temperature":0.7,"pith_summary":"This review argues that the nontrivial energy dependence of $e^+e^-$ annihilation into hadron pairs near threshold—sharp rises, resonance-like peaks, dips, and form-factor anomalies—is, in the vast majority of cases, produced by the strong interaction between the slow final-state hadrons rather than by new particles. The paper shows that one Schrödinger-equation treatment, with the production amplitude factored into an energy-independent short-distance part and a final-state enhancement factor, fits measured cross sections and form-factor ratios for $p\\bar p$, $n\\bar n$, $\\Lambda\\bar\\Lambda$, $\\Lambda_c\\bar\\Lambda_c$, $D^{(*)}\\bar D^{(*)}$, and $B^{(*)}\\bar B^{(*)}$ production with $\\chi^2/N_{\\rm df}$ between 1.2 and 1.6. It also predicts bound states of $\\Lambda\\bar\\Lambda$ and $\\Lambda_c\\bar\\Lambda_c$ below their thresholds. If this is right, many claimed exotic resonances near hadron-pair thresholds are reinterpreted as kinematic manifestations of hadron-hadron forces.","feed_headline":"No new particles needed for near-threshold hadron peaks","feed_subtitle":"One final-state interaction framework fits proton, Lambda, D, and B pair data, predicting two bound states.","key_machinery":"The carrying object is the final-state enhancement factor $F_L = \\frac{(2L+1)!!}{k^L L!}\\,\\frac{\\partial^L \\psi_L^{(R)}(0)}{\\partial r^L}$, which multiplies the bare production amplitude for a pair in partial wave $L$; the measured cross section is proportional to $k^{2L+1}|F_L|^2$. Here $\\psi_L^{(R)}(r)$ is the regular radial wave function of the nonrelativistic hadron pair obtained by solving the Schrödinger equation in a potential describing their interaction, and its value or derivative at $r=0$ captures the probability that the interacting pair is at short distances where it was formed. The paper computes this factor for single channels, optical potentials, coupled channels with nearby thresholds, tensor-force mixing of $S$ and $D$ waves, and Coulomb corrections, and uses it to convert fitted potential parameters into cross sections and form-factor ratios.","core_discovery":"The central claim is that the observed near-threshold structures in $e^+e^-$ annihilation into hadron pairs—including the rapid rise of $\\sigma(e^+e^-\\to p\\bar p)$ and $n\\bar n$, the $\\Lambda\\bar\\Lambda$ and $\\Lambda_c\\bar\\Lambda_c$ cross-section enhancements, the dips in $D^{(*)}\\bar D^{(*)}$ and $B^{(*)}\\bar B^{(*)}$ production, and the deviations of $|G_E/G_M|$ from unity—are all consequences of final-state interactions calculable from the Schrödinger equation. The paper exhibits a single framework in which the hadron pair's wave function at the origin, or its $L$-th derivative, multiplies an energy-independent bare amplitude, and with simple rectangular-well potentials, optical potentials for annihilation, Coulomb forces, tensor forces, and multichannel coupling it reproduces the data with $\\chi^2/N_{\\rm df}$ from 1.19 to 1.57. Consequently the paper concludes that it is not necessary to assume unknown intermediate particles to explain these features.","pith_inferences":["If the FSI explanation is universal, many resonances currently classified as exotic hadron candidates near two-body thresholds should be re-examined as threshold cusps or virtual-level enhancements; the burden of proof shifts to finding a peak whose shape cannot be produced by any reasonable hadronic potential.","The same $F_L$ enhancement applies to any production process with a short-distance source and slow outgoing hadrons, so the framework could be transferred to heavy-meson decays, two-photon collisions, and antiproton-proton annihilation, where similar threshold anomalies are observed.","A decisive test would be measuring the exclusive charged-to-neutral $B^{(*)}\\bar B^{(*)}$ ratios above $\\Upsilon(4S)$; the three parameter sets in the paper give different predictions for those ratios, so future data could select one set and rule out the others.","The fitted well parameters could be compared with first-principles strong-interaction calculations for the same channels; since the framework claims insensitivity to potential details, agreement would support the mechanism, while large discrepancies would indicate missing energy-dependent short-distance physics."],"forward_implications":["Near-threshold peaks in $e^+e^-$ cross sections, including the $\\psi(3770)$ and $\\Upsilon(4S)$ bumps, can arise from virtual levels in the interacting meson pairs, without treating them purely as $c\\bar c$ or $b\\bar b$ quarkonium states.","The model predicts a $\\Lambda\\bar\\Lambda$ bound state near $E_0\\approx -30$ MeV and a $\\Lambda_c\\bar\\Lambda_c$ bound state about 40 MeV below threshold; both should show up as sharp energy dependences in other production and decay processes.","Isospin-violating final-state interactions in $B^{(*)}\\bar B^{(*)}$ production can shift the measured $B^+$-$B^0$ mass difference by up to $\\delta M\\sim 0.4$ MeV, so FSI must be included in precision mass measurements.","Inelastic channels such as $e^+e^-\\to 6\\pi$ and $K^+K^-\\pi^+\\pi^-$ inherit a sharp drop at the $N\\bar N$ threshold from virtual $N\\bar N$ pairs in the intermediate state, and the fitted probability for an $I=1$ $N\\bar N$ pair to annihilate into six pions (62%) is close to the known annihilation fraction (56%)."],"supporting_citations":[{"why":"Provides the nucleon-antinucleon interaction model and fitted parameters used for $p\\bar p$ and $n\\bar n$ production and form factors.","marker":"[147]"},{"why":"Presents the $\\Lambda\\bar\\Lambda$ rectangular-well fit and parameter values reproduced in Section IV.","marker":"[125]"},{"why":"Presents the $\\Lambda_c\\bar\\Lambda_c$ coupled-channel potential fit including tensor forces and Coulomb interaction used in Section V.","marker":"[141]"},{"why":"Defines the six-channel $D^{(*)}\\bar D^{(*)}$ Schrödinger system with $L=1$ and the potential parametrization used in Section VIII.","marker":"[154]"},{"why":"Provides the three parameter sets for the $B^{(*)}$ meson interaction model whose predictions are compared in Section IX.","marker":"[167]"},{"why":"Supplies the partial-wave analysis of nucleon-antinucleon scattering data used to constrain the strong-interaction parameters of the model.","marker":"[14]"},{"why":"Provides the Born cross-section data near $\\psi(3770)$ used in the $D^{(*)}$ fit.","marker":"[105]"},{"why":"Gives the measurement of the charged-to-neutral $B$ production ratio used to test the isospin-violation prediction for the $B$ mass difference.","marker":"[165]"}],"fun_headline_variants":["Final-state interactions explain all near-threshold peaks","No new particles for hadron pair peaks near threshold","One framework fits proton, Lambda, D, B pair data","Near-threshold peaks from final-state forces only","One potential explains hadron resonances near threshold"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The short-distance quark-production amplitude is treated as a constant, apart from an empirical dipole form factor, in the final-state integral, so all near-threshold energy dependence is attributed to the hadronic wave function at the origin; if the bare amplitude carries additional energy dependence, for example from charmonium or bottomonium poles, the fitted potentials would absorb it and the explanation would become curve fitting.","fun_headline_variants_meta":{"raw":{"variants":["Final-state interactions explain all near-threshold peaks","No new particles for hadron pair peaks near threshold","One framework fits proton, Lambda, D, B pair data","Near-threshold peaks from final-state forces only","One potential explains hadron resonances near threshold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001508,"raw_usage":{"total_tokens":6002,"prompt_tokens":855,"completion_tokens":5147,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":5072}},"tokens_in":471,"tokens_out":5147,"duration_ms":39296,"temperature":1.0,"reasoning_tokens":5072,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:25:45.332842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics measurement in a channel where the predicted FSI cross section from the fitted potentials is smooth, such as a narrowly peaked structure in $e^+e^-\\to \\Lambda\\bar\\Lambda$ or $e^+e^-\\to B\\bar B$ that would require an energy-dependent short-distance amplitude or an additional pole, would refute the claim that all near-threshold features are final-state effects. Concretely, the model predicts $\\Lambda\\bar\\Lambda$ and $\\Lambda_c\\bar\\Lambda_c$ bound states near $-30$ MeV and $-40$ MeV below their thresholds; if high-precision data exclude narrow structures at those energies, the bound-state predictions fail.","supporting_citations":[{"cited_title":"Milstein and S.G","cited_arxiv_id":null,"evidence_quote":"Provides the nucleon-antinucleon interaction model and fitted parameters used for $p\\bar p$ and $n\\bar n$ production and form factors."},{"cited_title":"Feynman, Photon-Hadron Interactions, W.A","cited_arxiv_id":null,"evidence_quote":"Presents the $\\Lambda\\bar\\Lambda$ rectangular-well fit and parameter values reproduced in Section IV."},{"cited_title":"Dmitriev, A.I","cited_arxiv_id":null,"evidence_quote":"Presents the $\\Lambda_c\\bar\\Lambda_c$ coupled-channel potential fit including tensor forces and Coulomb interaction used in Section V."},{"cited_title":"Klempt, C","cited_arxiv_id":null,"evidence_quote":"Defines the six-channel $D^{(*)}\\bar D^{(*)}$ Schrödinger system with $L=1$ and the potential parametrization used in Section VIII."},{"cited_title":"Salnikov, A.E","cited_arxiv_id":null,"evidence_quote":"Provides the three parameter sets for the $B^{(*)}$ meson interaction model whose predictions are compared in Section IX."},{"cited_title":"H¨ usken, R.E","cited_arxiv_id":null,"evidence_quote":"Provides the Born cross-section data near $\\psi(3770)$ used in the $D^{(*)}$ fit."},{"cited_title":"Choudhury et al","cited_arxiv_id":null,"evidence_quote":"Gives the measurement of the charged-to-neutral $B$ production ratio used to test the isospin-violation prediction for the $B$ mass difference."}],"review_version":1}