Pith. sign in

REVIEW 3 major objections 5 minor 14 references

Vector meson production associated with a lepton pair in $e^+$ $e^-$ annihilation

T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read 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…

desk verdict 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. read the letter →

arxiv 2501.00592 v1 pith:FXORTUXP submitted 2024-12-31 hep-ph hep-ex

classification hep-phhep-ex
keywords vectormesonproductione+e-annihilationQEDdominanceinitial-stateemissiondoublelogarithmBESIIIBelleII
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [§II, Table II (ρ0 row)] 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.
  2. [§I (Abstract) and §II first paragraph] 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)).
  3. [§II (event yields after Table II)] 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.
minor comments (5)
  1. [References [6]] 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.
  2. [Appendix A] There are typos: 'fraction function' should read 'fracture function' in several places (for example, in the text below Eq. (A1) and before Eq. (A5)).
  3. [§III.A, Eq. (3)] 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.
  4. [Acknowledgments] There is a typo 'is is supported' in the acknowledgment of X. X.'s grant.
  5. [§II] The phrase 'There is no enough phase space' should read 'There is not enough phase space.'

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the cross-section predictions use externally measured g_V inputs and direct QED phase-space integrals, with self-citations only in non-load-bearing context.

full rationale

The target rates σ(e+e−→V l+l−) are computed from the four LO QED diagrams using an effective llV vertex (4πα gV/MV^2)γμ, with gV taken from the measured V→e+e− width (Table I). That is an external input, not a fit to the predicted cross sections; the predicted observables are new and independent enough that a measurement would test the model. No parameter is adjusted to make the Table II numbers come out. The VMD ansatz is a modeling assumption, but it is not a circular reduction: the paper does not define the input in terms of the output. The ISE/FSE separation and Eqs. (2)–(5) are direct calculations; the fracture-function appendix independently reproduces the leading log term of Eq. (3), which is a consistency check rather than the source of the prediction. Self-citations [3] and [9] appear only in a crossing-relation remark and in a background-motivation comparison against NRQCD; neither is load-bearing in deriving σ(e+e−→Vl+l−). The finite-ρ-width caveat is a model-robustness concern, not a circularity. Therefore no circular step is identified.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on four g_V values fitted to measured leptonic widths, the VMD effective vertex, the narrow-width approximation for the vector mesons, the Furry-theorem cancellation of ISE-FSE interference, and a lowest-order factorization ansatz in the appendix. The g_V values are external inputs and are not fitted to the target cross sections, so the circularity burden is low. No new physical entities are introduced; the QED fracture function is a mathematical construct, not a new degree of freedom.

free parameters (4)
  • g_V for J/ψ = 0.8319 GeV^2
    Determined from the measured J/ψ → e+e- width; used as the VMD coupling input for the J/ψ predictions.
  • g_V for ρ0 = 0.1177 GeV^2
    Determined from the measured ρ0 → e+e- width; the ρ0 is broad, adding uncertainty to the narrow-width treatment.
  • g_V for ω = 0.0360 GeV^2
    Determined from the measured ω → e+e- width; used as the VMD coupling input for the ω predictions.
  • g_V for φ = 0.0753 GeV^2
    Determined from the measured φ → e+e- width; used as the VMD coupling input for the φ predictions.
assumptions (4)
  • domain assumption The neutral vector meson couples to a lepton line via the effective vertex (4πα g_V/M_V^2) γ_μ, with g_V taken from the measured leptonic width of V → e+e- (VMD hypothesis).
    Invoked in Sec. II, Fig. 1. The paper calls the framework 'model independent', but this vertex is the VMD model for the γ*→V conversion.
  • domain assumption Each vector meson is treated as a stable on-shell particle with fixed mass M_V; the finite width of the meson is neglected.
    Assumed in Sec. II, Table I. This is not discussed; for the broad ρ0 (Γ ≈ 150 MeV) this could affect the predicted rates.
  • standard math The interference between ISE and FSE amplitudes vanishes because the corresponding cut diagram contains a closed muon loop with an odd number of C-odd vector bosons (Furry's theorem).
    Sec. II, Fig. 5. Standard QED theorem; the extension to exclusive differential distributions is asserted.
  • domain assumption In the triple collinear limit, the ISE cross section factorizes into a universal QED fracture function times the partonic cross section, and the fracture function obeys the convolution formula of Eq. (A5).
    Appendix A. The factorization is assumed, not proven to higher orders; the fracture function itself is derived at lowest order.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Vector meson production associated with a lepton pair in $e^+$ $e^-$ annihilation." pith.science (2026). https://pith.science/paper/FXORTUXP

@misc{pith2026250100592,
  author       = {Pith},
  title        = {Pith review of: Vector meson production associated with a lepton pair in $e^+$ $e^-$ annihilation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FXORTUXP}},
  note         = {Machine review of arXiv:2501.00592}
}
abstract

In this work, we investigate a novel production mechanism of vector mesons, exemplified by the production of a neutral vector meson associated with a lepton pair in $e^+e^-$ annihilation, i.e., $e^+e^-\to V l^+l^-$ ($V=J/\psi, \rho^0, \omega, \phi$, and $l=\mu, \tau$). These vector meson production channels can be precisely accounted within QED. The production rates of these processes are dominated by those diagrams where the vector meson is emitted from either the incident electron or positron, which exhibit a $\ln^2 m_l^2$ enhancement stemming from the triple collinear limit of leptons. Our numerical analysis indicates that the corresponding production rates are substantial enough to warrant the observation of these novel vector meson production channels at BESIII and Belle II experiments in near future.

Figures

Figures reproduced from arXiv: 2501.00592 by the authors.

Figure 1
Figure 1. FIG. 1: An effective [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Four lowest-order diagrams for [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Energy and angular distributions of [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Energy and angular distributions of [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: A representative cut diagram by stitching one ISE dia [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: A schematic illustration of the factorization progr [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Leading-order diagram for the QED fracture function [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 8 canonical work pages

  1. [1]

    Brambilla, S

    N. Brambilla, S. Eidelman, B. K. Heltsley, R. Vogt, G. T. B odwin, E. Eichten, A. D. Frawley, A. B. Meyer, R. E. Mitchell and V. Papadimitriou, et al. Eur. Phys. J. C 71, 1534 (2011) doi:10.1140/epjc/s10052-010-1534-9 [arXiv :1010.5827 [hep- ph]]

  2. [2]

    W represents a light-like Wilson line, whose role is to ensure QED gauge inv ariance

    The QED fraction function related to e →e(x) +l+l− can be defined as the following light-cone correlator [13]: Ml+l− e/e ({x;y1, q1⊥;y2, q2⊥}) = ∫ dw− 4π e−ixP +w− ⟨e(P )| ¯ψ(0,w −, 0⊥)W(∞,w −)†|l+(q1)l−(q2)⟩ × ⟨l+(q1)l−(q2) ⏐ ⏐γ+W(∞, 0)ψ(0) ⏐ ⏐e(P )⟩, (A1) wherex,y1, 2 signify the light-cone momentum fraction carried by the electron pa rton and by the bea...

  3. [3]

    G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125-1171 (1995) [erratum: Phys. Rev. D 55, 5853 (1997)] doi:10.1103/PhysRevD.55.5853 [arXiv:hep-ph/9407339 [h ep-ph]]

  4. [4]

    W. Chen, Y. Jia, Z. Mo, J. Pan and X. Xiong, Phys. Rev. D 104, no.9, 094023 (2021) doi:10.1103/PhysRevD.104.094023 [arXiv:2009.12363 [hep-ph]]

  5. [5]

    Navas et al

    S. Navas et al. [Particle Data Group], Phys. Rev. D 110, no.3, 030001 (2024) doi:10.1103/PhysRevD.110.030001

  6. [6]

    Bierlich, S

    C. Bierlich, S. Chakraborty, N. Desai, L. Gellersen, I. H elenius, P. Ilten, L. L¨ onnblad, S. Mrenna, S. Prestel and C.T. Preuss, et al. SciPost Phys. Codeb. 2022, 8 (2022) doi:10.21468/SciPostPhysCodeb.8 [arXiv:2203. 11601 [hep-ph]]

  7. [7]

    Hahn, Comput

    T. Hahn, Comput. Phys. Commun. 168, 78-95 (2005) doi:10.1016/j.cpc.2005.01.010 [arXiv:hep -ph/0404043 [hep-ph]]

  8. [8]

    Davier, M

    M. Davier, M. E. Peskin and A. Snyder, [arXiv:hep-ph/060 6155 [hep-ph]]

Show all 14 references
  1. [9]

    G. T. Bodwin, E. Braaten, J. Lee and C. Yu, Phys. Rev. D 74, 074014 (2006) doi:10.1103/PhysRevD.74.074014 [arXiv:hep-ph/0608200 [hep-ph]]

  2. [10]

    W. L. Sang, F. Feng, Y. Jia, Z. Mo, J. Pan and J. Y. Zhang, Phy s. Rev. Lett. 131, no.16, 161904 (2023) doi:10.1103/PhysRevLett.131.161904 [arXiv:2306.11538 [hep-ph]]

  3. [11]

    X. D. Huang, B. Gong, R. C. Niu, H. M. Yu and J. X. Wang, JHEP 02, 055 (2024) doi:10.1007/JHEP02(2024)055 [arXiv:2311.04751 [hep-ph]]

  4. [12]

    Fleming, Phys

    S. Fleming, Phys. Rev. D 50, 5808-5815 (1994) doi:10.1103/PhysRevD.50.5808 [arXiv: hep-ph/9403396 [hep-ph]]

  5. [13]

    Trentadue and G

    L. Trentadue and G. Veneziano, Phys. Lett. B 323, 201-211 (1994) doi:10.1016/0370-2693(94)90292-5

  6. [14]

    Collins, Foundations of Perturbative QCD , Camb

    J. Collins, Foundations of Perturbative QCD , Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 32, 1-624 (2011) Cambridge University Press, 2011, ISBN 978-1-009-40184-5, 978-1-00 9-40183-8, 978-1-009-40182-1 doi:10.1017/978100940184 5

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.