{"id":"cbfa6389-d57d-4102-b5b6-f5435eb84b0b","arxiv_id":"2501.00592","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"New QED-based predictions show that e+e- → V l+ l- (V = J/ψ, ρ0, ω, φ) production rates at BESIII and Belle II are large enough to observe, with J/ψ μ+ μ- forming a significant background to J/ψ pair production.","lead":"This paper computes the cross sections for producing a neutral vector meson (J/ψ, ρ0, ω, φ) together with a muon or tau pair in electron-positron annihilation at BESIII and Belle II energies. It predicts that these rates are large enough to be observed, and that one channel, J/ψ plus a muon pair, is a major background to searches for double charmonium production.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ρ finite-width treatment is the main quantitative weak point; the zero-width Table II row for ρ is not robust, though narrow-meson predictions stand.","rationale":"The reader's weakest_assumption correctly flags the neglect of the vector meson width, and that is the concern I regard as load-bearing. However, I do not fully agree with framing the VMD vertex itself as a shaky assumption: for a stable on-shell vector meson, the coupling gV is fixed by the exact matrix element ⟨V|J_EM|0⟩ = gV ε*, and the photon that converts to V carries q^2 = M_V^2, so the contact vertex (4πα gV/M_V^2)γμ is a valid leading-order QED+VMD representation. The genuine weak point is the zero-width treatment of the broad ρ0, where 'on-shell' is not a well-defined observable without an invariant-mass cut. This concern affects the largest numerical entries in Table II, but it does not overturn the central observability claim because the J/ψ, ω, and φ channels already provide thousands of events at BESIII and tens of thousands at Belle II. The reader's CONDITIONAL verdict already accommodates this caveat, so I recommend no change to the verdict. The proposed Breit-Wigner integration test would settle whether the ρ row shifts materially.","tokens_in":8616,"tokens_out":36875,"duration_ms":399231,"concrete_test":"Recompute σ(e+e−→ρ μ+μ−) by replacing the stable-ρ amplitude with a Breit-Wigner integral over the ρ invariant mass m: σ_ρ = ∫ dm^2 (1/π) [m Γ_ρ(m)/((m^2−M_ρ^2)^2 + m^2 Γ_ρ^2(m))] σ_0(m), where Γ_ρ(m) uses the P-wave ππ phase-space factor Γ_ρ(m) = Γ_ρ (M_ρ/m) ((m^2−4m_π^2)/(M_ρ^2−4m_π^2))^{3/2}, and σ_0(m) is the same formula as Eq. (2)/(4) with M_V→m and g_V scaled according to the VMD off-shell prescription, e.g. g_ρ(m) = g_ρ(M_ρ) (M_ρ/m)^n for standard VMD choices n = 0, 1, 2. Repeat the integration for Belle II and BESIII energies and compare the result with the ρ rows in Table II. If the integrated rate changes by more than about 20% or depends strongly on the choice of n, the zero-width ρ prediction is not robust; if it changes by less, the finite-width concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative claim that observation prospects are bright rests on the cross sections in Table II, all of which are computed under a zero-width, on-shell treatment of the vector meson. For the narrow J/ψ, ω, and φ this is harmless (Γ/M below a few percent). For the ρ0, however, Γ ≈ 150 MeV is about 20% of M_ρ, and an on-shell 'ρ0' is not an asymptotic state. Experimentally ρ0 is reconstructed through the π+π− invariant mass in some signal window, and the production amplitude should be integrated over the ρ lineshape with an energy-dependent width and an off-shell γ*→ρ transition form factor. The zero-width prediction of 1558 fb at BESIII and 505 fb at Belle II is the largest row in Table II, so this is not a negligible caveat. Because the paper provides no uncertainty estimates, the ρ row could shift by O(Γ/M) or more, and the shift is model dependent through the off-shell ρππ coupling. The J/ψ μ+μ− channel alone (4458 events at BESIII, 83k at Belle II) is enough to support the qualitative observability claim, so this concern does not invalidate the paper's main conclusion; it does mean the ρ-specific rates in Table II should not be used as quantitative predictions until the lineshape integration is performed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the reactions e+e− → V l+l− (V = J/ψ, ρ0, ω, φ; l = μ, τ) at BESIII and Belle II energies. Within a VMD framework, the neutral vector meson is emitted from a lepton line through an effective llV vertex containing g_V, whose values are taken from measured leptonic widths. The four lowest-order diagrams are split into ISE and FSE classes; the interference is argued to vanish by Furry's theorem. Analytical expressions (Eqs. (2)-(5)) and numerical cross sections (Table II) are given, showing that ISE dominates via a ln^2 m_l^2 double-log, and event yields at BESIII and Belle II are estimated. The paper concludes that several channels are observable.","tokens_in":8850,"tokens_out":12227,"duration_ms":123833,"significance":"The paper provides a transparent, falsifiable set of predictions for a class of processes that have not been considered before. The analytic treatment is a strength: the Furry-theorem cancellation of ISE-FSE interference is correct, and the identification of the double-log enhancement from the triple-collinear limit is convincing. The inputs g_V are not fitted to the target cross sections, so the numbers are genuine predictions. The appendix's QED fracture function (Eq. (A4)) is a useful byproduct. However, the quantitative reliability is limited by the zero-width treatment of the ρ0 and by the absence of uncertainty estimates, and the claim of model independence is overstated.","major_comments":[{"comment":"For the ρ0, the finite width (Γ ≈ 150 MeV, about 19% of M_ρ) cannot be neglected for an on-shell particle in the final state. The quoted cross sections of 1558.6 fb (BESIII) and 505.5 fb (Belle II) are obtained with a zero-width, on-shell treatment, and are the largest entries in Table II. A realistic estimate requires convolution of the production amplitude over the π+π− invariant mass with an energy-dependent width and an off-shell γ*→ρ transition form factor; the result will depend on model assumptions and can shift by O(Γ/M). Please either perform this lineshape integration or present the ρ0 numbers only with an explicit caveat, and exclude them from quantitative observability claims.","section":"§II, Table II (ρ0 row)"},{"comment":"The statements that the processes 'can be precisely accounted within QED' and that the framework is 'model independent' overstate the content of the calculation. The calculation uses the VMD contact vertex (4πα g_V/M_V^2)γ_μ for the llV coupling; g_V from the leptonic width fixes the on-shell γ*→V amplitude, but off-shell effects and corrections from higher resonances or continuum are not quantified. Table II reports cross sections to six significant figures without any uncertainty. Please soften the QED/model-independence wording, quote the input uncertainties from the PDG leptonic widths, and give a theoretical-error estimate (for example, a range obtained by varying the VMD scale or the cutoff in Eq. (3)).","section":"§I (Abstract) and §II first paragraph"},{"comment":"The text states that BESIII and Belle II 'have already achieved about 20 fb−1 and 1500 fb−1 integrated luminosity.' The Belle II figure is not the currently recorded luminosity (which is of order 500 fb−1 as of 2024); it appears to be a planned value. The quoted event yields for Belle II in Section II are therefore overestimated by about a factor of three. Please correct the luminosity and recompute the yields, or clearly label the yields as projections for 1500 fb−1. The qualitative conclusion remains unchanged.","section":"§II (event yields after Table II)"}],"minor_comments":[{"comment":"The HCubature package is not the same as the CUBA library; reference [6] describes CUBA, not HCubature. Please cite the actual package (for example, S. G. Johnson, HCubature.jl) or state that CUBA was used.","section":"References [6]"},{"comment":"There are typos: 'fraction function' should read 'fracture function' in several places (for example, in the text below Eq. (A1) and before Eq. (A5)).","section":"Appendix A"},{"comment":"The notation 'ln^2 m_l^2' is shorthand; it should be written explicitly, for example as ln^2(s/m_l^2), to avoid confusion.","section":"§III.A, Eq. (3)"},{"comment":"There is a typo 'is is supported' in the acknowledgment of X. X.'s grant.","section":"Acknowledgments"},{"comment":"The phrase 'There is no enough phase space' should read 'There is not enough phase space.'","section":"§II"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the core derivation is sound. The main issues requiring revision are the explicit treatment (or clear caveat) of the broad ρ0 width and the overclaim of model independence; both are fixable. The novelty appears adequate, though a quick check for prior work on e+e− → J/ψμ+μ− as a background would strengthen the introduction. The self-citation of reference [9] is appropriate and does not affect the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe paper is a straightforward but complete LO QED calculation of e+e- -> V l+l- for V = J/psi, rho0, omega, phi and l = mu, tau, using VMD couplings. The genuinely new pieces are the full ISE/FSE decomposition, the double-log analysis showing why ISE dominates, and the QED fracture function in Appendix A. The derivations in Sec. III are solid: the Furry-theorem argument that ISE/FSE interference cancels is correct, and the asymptotic formulas match the numerics. The g_V inputs come from measured leptonic widths, not from a fit to the target cross sections, so the predictions are real observables.\n\nThe main soft spot is the rho. Table II gives 1558 fb for rho0 mu+mu- at BESIII, the largest row, but the calculation treats the rho as a zero-width on-shell state. With Gamma_rho ~ 150 MeV, that approximation is not obviously safe, and the paper gives no uncertainty estimate. The narrow mesons (J/psi, omega, phi) should be robust. The J/psi mu+mu- rate alone (228 fb at BESIII, 55 fb at Belle II) supports the qualitative observability claim, so the rho issue doesn't sink the paper, but the rho-specific numbers shouldn't be quoted as quantitative predictions until a lineshape integration is done.\n\nAlso, calling the framework \"model independent\" overstates things. It is VMD with a fixed on-shell coupling; the off-shell continuation is assumed. That is worth a sentence of hedging, and an estimate of the theoretical uncertainty.\n\nThe fracture function appendix is a nice byproduct, though the derivation is compressed. It reproduces the leading double log, which is enough for the paper's purpose.\n\nCitation pattern looks fine; the comparison to the authors' own NRQCD prediction for J/psi J/psi is relevant, not self-promotional.\n\nBottom line: this is a clean, useful phenomenological calculation. It deserves a serious referee. My recommendation would be to send it out, with a request that the authors quantify the rho-width effect and the VMD uncertainty before publication.","headline":"A clean, standard LO QED calculation of e+e- -> V l+l- that makes a defensible observability claim; the rho row needs a width treatment before use, but the J/psi channel stands.","tokens_in":9446,"tokens_out":1952,"would_cite":true,"duration_ms":19032,"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 paper predicts that $e^+e^-\\to V l^+l^-$ for $V=J/\\psi,\\rho^0,\\omega,\\phi$ and $l=\\mu,\\tau$ is a QED-controlled production mechanism with rates large enough to observe at BESIII and Belle II, dominated by initial-state emission with a…","keywords":["vector meson production","e+e- annihilation","QED","vector meson dominance","initial-state emission","double logarithm","BESIII","Belle II"],"falsifier":"The cleanest test is to measure $\\sigma(e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-)$ at BESIII at $\\sqrt{s}=3.77$ GeV; the paper predicts about 228 fb, almost all from initial-state emission, with the $J/\\psi$ energy spectrum peaking near $E_{J/\\psi}\\approx(s+M_{J/\\psi}^2-4m_\\mu^2)/(2\\sqrt{s})$ and the angular distribution peaked at $\\theta\\approx 0,\\pi$. An observed rate much smaller than the prediction, or a distribution lacking those collinear peaks, would falsify the VMD-based QED description.","tokens_in":8390,"feed_emoji":"⚛️","tokens_out":10138,"duration_ms":95783,"temperature":0.7,"pith_summary":"This paper argues that neutral vector mesons can be produced in $e^+e^-$ annihilation together with a lepton pair through a purely QED mechanism: a virtual photon converts into the vector meson while the leptons are radiated. The authors compute leading-order cross sections for $e^+e^-\\to V l^+l^-$ at BESIII ($\\sqrt{s}=3.77$ GeV) and Belle II ($\\sqrt{s}=10.58$ GeV) and find sizable rates, for example $\\sigma(e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-)\\approx 228$ fb at BESIII and $\\approx 55.6$ fb at Belle II. They show that the dominant contribution comes from diagrams where the vector meson is emitted from an incident electron or positron, and that this initial-state emission carries a $\\ln^2 m_l^2$ double-logarithmic enhancement arising from the triple collinear limit in which the final leptons fly nearly parallel to the beam. The central message is that these channels should be observable with existing integrated luminosities, and that $e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-$ is an important background for double charmonium searches.","feed_headline":"New QED channel puts vector mesons within reach at BESIII and Belle II","feed_subtitle":"Double-log QED enhancement makes J/ψ, ρ, ω and φ production visible in current data.","key_machinery":"The argument turns on the effective lepton-lepton-vector vertex $(4\\pi\\alpha g_V/M_V^2)\\gamma_\\mu$, built from vector meson dominance, which lets each neutral vector meson be emitted directly from a lepton line; the coupling $g_V$ is fixed by the measured leptonic width of $V$. The analytic control comes from expanding the ISE cross section in powers of the final lepton-pair invariant mass $s_{ll}$, which isolates the leading $\\ln^2 m_l^2$ term, and from a QED fracture function---a light-cone correlator giving the probability to find an electron parton accompanied by two forward leptons---that reproduces the same double logarithm in the triple collinear limit. The fracture function construction is the mechanism that makes the leading enhancement calculable and universal rather than process-specific.","core_discovery":"The paper's central claim is that the reactions $e^+e^-\\to V l^+l^-$ proceed at leading order through four QED diagrams, split into initial-state emission (ISE), where $V$ is radiated from the incoming electron or positron, and final-state emission (FSE), where it is radiated from one of the produced leptons. The integrated cross section is dominated by the ISE diagrams, whose asymptotic behavior is controlled by a double logarithm $\\ln^2 m_l^2$ generated by the triple collinear limit in which the two final leptons move nearly parallel to the beam electron. The FSE contribution instead carries a Sudakov double logarithm $\\ln^2(s/M_V^2)$, and the interference between the two classes vanishes after integration by Furry's theorem. Numerically, the paper predicts $\\sigma(e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-)\\approx 227.9$ fb at BESIII and $55.6$ fb at Belle II, with event yields in the thousands to hundreds of thousands at the quoted integrated luminosities.","pith_inferences":["The paper does not say this, but the same QED fracture function should apply to $e^+e^-\\to Z^0\\,\\mu^+\\mu^-$ at future colliders, where the triple-collinear double logarithm would appear in the $\\mu^+\\mu^-$ invariant-mass distribution.","A natural test of the core assumption that the paper does not perform is to include the finite width of the $\\rho^0$; a line-shape analysis of the $\\rho^0\\mu^+\\mu^-$ channel would either confirm the VMD contact vertex or expose its breakdown.","One could also suppress the $J/\\psi\\,\\mu^+\\mu^-$ background to double charmonium by cutting on the $\\mu^+\\mu^-$ invariant mass or the $J/\\psi$ emission angle, since ISE kinematics concentrate the signal near the beam direction.","At center-of-mass energies far above the $B$ factory, the FSE Sudakov logarithm $\\ln^2(s/M_V^2)$ could catch up with the ISE $\\ln^2 m_l^2$ term, so the paper's hierarchy ISE $\\gg$ FSE is tied to the specific energies considered and may not persist."],"forward_implications":["Existing BESIII data ($\\sim20$ fb$^{-1}$) should contain about 4,500 $J/\\psi\\,\\mu^+\\mu^-$ events and 31,000 $\\rho^0\\mu^+\\mu^-$ events; Belle II data ($\\sim1{,}500$ fb$^{-1}$) should contain about 83,000 and 758,000 events, respectively.","The $e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-$ rate is roughly 25 times the NRQCD prediction for $e^+e^-\\to J/\\psi J/\\psi$, so this channel must be subtracted as a leading background in double charmonium searches.","Because ISE dominates, the vector meson's energy spectrum peaks near its kinematic maximum and its angular distribution peaks along the beam directions, a distinctive signature that separates this mechanism from other production processes.","At Belle II, the $V\\tau^+\\tau^-$ channels, which have no phase space at BESIII, should also be observable, with predicted cross sections between 6 and 144 fb.","Predicted $V l^+l^-$ rates are considerably larger than those for $e^+e^-\\to V_1V_2$ through double vector meson dominance."],"supporting_citations":[{"why":"supplies the particle masses used in the numerics","marker":"[4]"},{"why":"supplies the running QED coupling values at the two collider energies","marker":"[5]"},{"why":"the four lowest-order diagrams are obtained by crossing the $J/\\psi\\to\\mu^+\\mu^-e^+e^-$ decay diagrams","marker":"[3]"},{"why":"the HCubature package performs the phase-space integration for the cross-section values","marker":"[6]"},{"why":"provide the double-VMD $e^+e^-\\to V_1V_2$ rates used as a comparison baseline","marker":"[7, 8]"},{"why":"gives the NRQCD prediction for $e^+e^-\\to J/\\psi J/\\psi$ that the $J/\\psi\\,\\mu^+\\mu^-$ channel backgrounds","marker":"[9]"},{"why":"gives an independent NRQCD prediction for the same double charmonium background comparison","marker":"[10]"},{"why":"supplies the fragmentation treatment of $Z^0\\to J/\\psi l^+l^-$ used to reproduce the FSE double logarithm","marker":"[11]"},{"why":"introduces the fracture function that underlies the QED factorization derivation of the ISE double logarithm","marker":"[12]"},{"why":"provides the light-cone correlator technology used to define and compute the QED fracture function","marker":"[13]"}],"fun_headline_variants":["Double-log QED boost reveals vector mesons in e+e- collisions","QED double log lifts vector meson yields at BESIII and Belle II","Vector mesons via lepton-pair QED in e+e- collisions","New QED route to J/psi and friends at BESIII and Belle II","Lepton-pair channel exposes vector mesons in e+e- data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything quantitative rests on vector meson dominance: that a photon emitted from a lepton line converts into an on-shell vector meson with the coupling $g_V$ extracted from the measured leptonic width, and that the vector meson's finite width, especially the broad $\\rho^0$, can be neglected.","fun_headline_variants_meta":{"raw":{"variants":["Double-log QED boost reveals vector mesons in e+e- collisions","QED double log lifts vector meson yields at BESIII and Belle II","Vector mesons via lepton-pair QED in e+e- collisions","New QED route to J/psi and friends at BESIII and Belle II","Lepton-pair channel exposes vector mesons in e+e- data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000612,"raw_usage":{"total_tokens":2844,"prompt_tokens":943,"completion_tokens":1901,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":1799}},"tokens_in":559,"tokens_out":1901,"duration_ms":12164,"temperature":1.0,"reasoning_tokens":1799,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:47:49.687243+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The cleanest test is to measure $\\sigma(e^+e^-\\to J/\\psi\\,\\mu^+\\mu^-)$ at BESIII at $\\sqrt{s}=3.77$ GeV; the paper predicts about 228 fb, almost all from initial-state emission, with the $J/\\psi$ energy spectrum peaking near $E_{J/\\psi}\\approx(s+M_{J/\\psi}^2-4m_\\mu^2)/(2\\sqrt{s})$ and the angular distribution peaked at $\\theta\\approx 0,\\pi$. An observed rate much smaller than the prediction, or a distribution lacking those collinear peaks, would falsify the VMD-based QED description.","supporting_citations":[{"cited_title":"Four-lepton decays of neutral vector mesons","cited_arxiv_id":"2009.12363","evidence_quote":"supplies the particle masses used in the numerics"},{"cited_title":"Exclusive Two-Vector-Meson Production from e^+ e^- Annihilation","cited_arxiv_id":"hep-ph/0608200","evidence_quote":"gives the NRQCD prediction for $e^+e^-\\to J/\\psi J/\\psi$ that the $J/\\psi\\,\\mu^+\\mu^-$ channel backgrounds"},{"cited_title":"J/Psi Production from Electromagnetic Fragmentation in Z decay","cited_arxiv_id":"hep-ph/9403396","evidence_quote":"introduces the fracture function that underlies the QED factorization derivation of the ISE double logarithm"},{"cited_title":"Trentadue and G","cited_arxiv_id":null,"evidence_quote":"provides the light-cone correlator technology used to define and compute the QED fracture function"}],"review_version":1}