pith. sign in

arxiv: 2606.17827 · v1 · pith:AQ7RSRQRnew · submitted 2026-06-16 · ✦ hep-ph · hep-ex

Semi-analytical results for e^+e^-to J/psi + X_{{rm non\,}cbar{c}} up to mathcal{O}(α_s v²) at B factories

Pith reviewed 2026-06-27 00:17 UTC · model grok-4.3

classification ✦ hep-ph hep-ex
keywords NRQCD factorizationJ/ψ productionB factoriescolor-singlethigher-order correctionsangular distributionfeeddown
0
0 comments X

The pith

The O(α_s v²) correction to color-singlet J/ψ production in e⁺e⁻ collisions is obtained for the first time, verifying NRQCD factorization at B factories.

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

The paper calculates the strong-coupling, velocity, and mixed corrections up to O(α_s v²) for the color-singlet part of e⁺e⁻ → J/ψ + X without charm quarks. Short-distance coefficients are extracted as high-order asymptotic series in the mass ratio r = m_c/√s. The new mixed term turns out small, the factorization structure holds, and after adding ψ(2S) feeddown the total rate matches the Belle measurement while the predicted angular shape does not.

Core claim

Within NRQCD factorization the O(α_s), O(v²) and O(α_s v²) corrections to both the unpolarized cross section and the J/ψ angular distribution are computed for the color-singlet channel. The O(α_s v²) short-distance coefficients are obtained for the first time via the differential equation method as expansions in r = m_c/√s up to r^40; these reproduce exact results to relative precision 10^{-14} for the rate and 10^{-7} for the angular parameter at B-factory energies. With the scale choice μ_R = √s/2 the O(α_s) term reaches roughly 50 % of the leading-order rate while the two velocity corrections remain accidentally small. After including feeddown from ψ(2S) the predicted cross section is 0.5

What carries the argument

Short-distance coefficients for color-singlet J/ψ production computed to O(α_s v²) as asymptotic expansions in r = m_c/√s up to order r^40 via the differential equation method.

If this is right

  • The O(α_s) correction is large (~50 % of LO) while O(v²) and O(α_s v²) are small, so the perturbative and velocity expansions can be kept separate at present precision.
  • Feeddown from ψ(2S) is required to bring the predicted cross section into agreement with Belle data.
  • The angular-distribution parameter remains discrepant with experiment by >2σ, indicating that further contributions or measurements are needed.
  • The asymptotic expansions reproduce exact results to 10^{-14} relative accuracy for the rate at B-factory energies.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same expansion technique could be used for other charmonium states or different production channels to test factorization at this order.
  • The persistent angular-distribution mismatch may point to the importance of color-octet channels or still-higher-order terms once the LDMEs are allowed to adjust at O(α_s v²).
  • Re-extracting the LDMEs from data that now includes the new coefficients would provide a direct test of whether the present scale choice and matrix-element values stay stable.

Load-bearing premise

The long-distance matrix elements and renormalization-scale choice taken from earlier lower-order fits remain appropriate once the new O(α_s v²) coefficients are added.

What would settle it

An exact (non-asymptotic) evaluation of any O(α_s v²) short-distance coefficient at √s ≈ 10.6 GeV that differs from the r^40 expansion result by more than the stated 10^{-14} relative error on the cross section.

read the original abstract

Within the NRQCD factorization framework, we investigate the color-singlet contribution to $e^+e^- \to J/\psi + X_{{\rm non\,}c\bar{c}}$ at B factories, computing the $\mathcal{O}(\alpha_s)$, $\mathcal{O}(v^2)$, and $\mathcal{O}(\alpha_s v^2)$ corrections to both the unpolarized cross section and the $J/\psi$ angular distribution. The $\mathcal{O}(\alpha_s v^2)$ correction is obtained for the first time, and the validity of NRQCD factorization at this order is explicitly verified. Using the differential equation method, the short-distance coefficients are obtained as asymptotic expansions in $r = m_c/\sqrt{s}$ up to $r^{40}$, which reproduce exact results with high precision at B factory energies, achieving relative errors around $10^{-14}$ for the cross section and around $10^{-7}$ for the angular distribution. Notably, with the same input parameters, our $\mathcal{O}(\alpha_s)$ and $\mathcal{O}(v^2)$ corrections are consistent with those reported in the literature. Phenomenologically, the $\mathcal{O}(\alpha_s)$ correction (with $\mu_R=\sqrt{s}/2$) reaches about $50\%$ of the leading-order cross section, while the $\mathcal{O}(v^2)$ and $\mathcal{O}(\alpha_s v^2)$ corrections are accidentally small. After including feeddown contributions from $\psi(2S)$, the predicted cross section $0.530_{-0.113}^{+0.122}$ pb agrees with the {\tt Belle} measurement within uncertainties. However, the predicted angular distribution parameter $0.120_{-0.019}^{+0.027}$ deviates from the experimental value $5.71\pm 2.51$ by more than $2\sigma$, calling for further experimental and theoretical investigations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript computes the color-singlet contributions to e⁺e⁻ → J/ψ + X_non-c c̄ at B-factory energies in NRQCD factorization, obtaining the O(α_s), O(v²), and first-time O(α_s v²) corrections to both the unpolarized cross section and the J/ψ angular distribution. Short-distance coefficients are derived via the differential-equation method as asymptotic expansions in r = m_c/√s up to r⁴⁰, reproducing exact results to relative precision ~10^{-14} (cross section) and ~10^{-7} (angular). Lower-order terms match literature values with the same inputs; after ψ(2S) feeddown the predicted cross section is 0.530 pb (agreeing with Belle within uncertainties) while the angular parameter 0.120 deviates by >2σ.

Significance. If the results hold, the work supplies the first O(α_s v²) short-distance coefficients for this process together with high-order semi-analytical expansions that serve as a practical tool at B-factory energies. Explicit verification of NRQCD factorization at the new order and consistency of the O(α_s) and O(v²) pieces with prior literature are concrete strengths. The phenomenological section illustrates both the size of perturbative corrections and a tension in the angular observable that warrants further study.

major comments (1)
  1. [Abstract] Abstract and phenomenological discussion: the quoted agreement of the total cross section (0.530 pb) with Belle data and the assessment that O(v²) and O(α_s v²) corrections are “accidentally small” both rest on LDMEs and the choice μ_R = √s/2 taken unchanged from lower-order literature; because the O(α_s) term already reaches ~50 % of LO with this scale, the numerical impact of the new O(α_s v²) coefficient on the final prediction is not independently tested by a refit or by a scale-variation study at the enlarged perturbative order.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive evaluation and the detailed comment. We respond to the single major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract and phenomenological discussion: the quoted agreement of the total cross section (0.530 pb) with Belle data and the assessment that O(v²) and O(α_s v²) corrections are “accidentally small” both rest on LDMEs and the choice μ_R = √s/2 taken unchanged from lower-order literature; because the O(α_s) term already reaches ~50 % of LO with this scale, the numerical impact of the new O(α_s v²) coefficient on the final prediction is not independently tested by a refit or by a scale-variation study at the enlarged perturbative order.

    Authors: We agree that the quoted numerical results and the characterization of the higher-order corrections as accidentally small are obtained with the LDMEs and scale choice μ_R=√s/2 taken from the existing literature. These inputs were deliberately retained to permit a direct, apples-to-apples comparison of the newly computed O(α_s v²) short-distance coefficients with the lower-order terms already published. A global refit of the LDMEs that incorporates the O(α_s v²) corrections, or a dedicated scale-variation study at the enlarged perturbative order, would constitute a separate phenomenological analysis and lies outside the scope of the present work, whose primary goal is the derivation and verification of the perturbative coefficients themselves. We will add a brief clarifying sentence in the abstract and in the phenomenological discussion to make this choice of inputs explicit. revision: partial

Circularity Check

0 steps flagged

No circularity; short-distance coefficients derived independently via differential equations

full rationale

The paper computes the new O(α_s v²) short-distance coefficients from first principles using the differential equation method with asymptotic expansions in r = m_c/√s up to r^40. Lower-order results are reproduced for consistency with external literature but are not used to define or fit the new term. LDMEs and μ_R = √s/2 are adopted as standard external inputs for phenomenology without refitting or redefinition at the new order; the factorization verification applies only to the independently computed coefficients. No step reduces the central result to a fitted quantity, self-citation chain, or input by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The calculation rests on the NRQCD factorization framework with standard inputs (LDMEs, m_c, α_s) taken from prior literature; the new content is the short-distance coefficients at O(α_s v²).

free parameters (2)
  • renormalization scale μ_R = √s/2
    Chosen as √s/2 so that the O(α_s) correction reaches ~50% of the leading-order cross section.
  • long-distance matrix elements
    Standard NRQCD inputs; values taken from literature or previous fits and not re-determined at the new order.
axioms (1)
  • domain assumption NRQCD factorization remains valid at O(α_s v²)
    The paper states that this validity is explicitly verified, yet the separation of short- and long-distance physics is the foundational premise of the entire framework.

pith-pipeline@v0.9.1-grok · 5917 in / 1626 out tokens · 48480 ms · 2026-06-27T00:17:53.078209+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

59 extracted references · 39 linked inside Pith

  1. [1]

    Keung,Off Resonance Production of Heavy Vector Quarkonium States ine+e− Annihilation,Phys

    W.-Y. Keung,Off Resonance Production of Heavy Vector Quarkonium States ine+e− Annihilation,Phys. Rev. D23(1981) 2072

  2. [2]

    Kiselev, A.K

    V.V. Kiselev, A.K. Likhoded and M.V. Shevlyagin,Double charmed baryon production at B factory,Phys. Lett. B332(1994) 411 [hep-ph/9408407]

  3. [3]

    Braaten and Y.-Q

    E. Braaten and Y.-Q. Chen,Signature for color octet production ofJ/ Psi ine +e− annihilation,Phys. Rev. Lett.76(1996) 730 [hep-ph/9508373]

  4. [4]

    Yuan, C.-F

    F. Yuan, C.-F. Qiao and K.-T. Chao,PromptJ/ψproduction ate +e− colliders,Phys. Rev. D56(1997) 321 [hep-ph/9703438]

  5. [5]

    Cho and A.K

    P.L. Cho and A.K. Leibovich,Color singlet psi(Q) production at e+ e- colliders,Phys. Rev. D54(1996) 6690 [hep-ph/9606229]

  6. [6]

    Beneke, I.Z

    M. Beneke, I.Z. Rothstein and M.B. Wise,Kinematic enhancement of nonperturbative corrections to quarkonium production,Phys. Lett. B408(1997) 373 [hep-ph/9705286]

  7. [7]

    Schuler,Testing factorization of charmonium production,Eur

    G.A. Schuler,Testing factorization of charmonium production,Eur. Phys. J. C8(1999) 273 [hep-ph/9804349]. – 18 –

  8. [8]

    Liu, Z.-G

    K.-Y. Liu, Z.-G. He and K.-T. Chao,Production of J / psi + c anti-c through two photons in e+ e- annihilation,Phys. Rev. D68(2003) 031501 [hep-ph/0305084]

  9. [9]

    Liu, Z.-G

    K.-Y. Liu, Z.-G. He and K.-T. Chao,Inclusive charmonium production via doublec¯cine +e− annihilation,Phys. Rev. D69(2004) 094027 [hep-ph/0301218]

  10. [10]

    Fleming, A.K

    S. Fleming, A.K. Leibovich and T. Mehen,Resumming the Color Octet Contribution to e+e− →J/ψ+X,Phys. Rev. D68(2003) 094011 [hep-ph/0306139]

  11. [11]

    Zhang and K.-T

    Y.-J. Zhang and K.-T. Chao,Double charm production e+ e- —>J / psi + c anti-c at B factories with next-to-leading order QCD correction,Phys. Rev. Lett.98(2007) 092003 [hep-ph/0611086]

  12. [12]

    Ma, Y.-J

    Y.-Q. Ma, Y.-J. Zhang and K.-T. Chao,QCD correction toe +e− →J/ψ+ggat B Factories, Phys. Rev. Lett.102(2009) 162002 [0812.5106]

  13. [13]

    Gong and J.-X

    B. Gong and J.-X. Wang,Next-to-leading-order QCD corrections to e+e- –>J/psi(cc) at the B factories,Phys. Rev. D80(2009) 054015 [0904.1103]

  14. [14]

    Gong and J.-X

    B. Gong and J.-X. Wang,Next-to-Leading-Order QCD Corrections to e+ e- —>J/psi gg at the B Factories,Phys. Rev. Lett.102(2009) 162003 [0901.0117]

  15. [15]

    Z.-G. He, Y. Fan and K.-T. Chao,Relativistic correction to e+ e- —>J/psi + gg at B factories and constraint on color-octet matrix elements,Phys. Rev. D81(2010) 054036 [0910.3636]

  16. [16]

    Jia,Color-singlet relativistic correction to inclusive J/psi production associated with light hadrons at B factories,Phys

    Y. Jia,Color-singlet relativistic correction to inclusive J/psi production associated with light hadrons at B factories,Phys. Rev. D82(2010) 034017 [0912.5498]

  17. [17]

    Zhang, Y.-Q

    Y.-J. Zhang, Y.-Q. Ma, K. Wang and K.-T. Chao,QCD radiative correction to color-octet J/ψinclusive production at B Factories,Phys. Rev. D81(2010) 034015 [0911.2166]

  18. [18]

    Shao,Initial state radiation effects in inclusiveJ/ψproduction at B factories,JHEP 04(2014) 182 [1402.5840]

    H.-S. Shao,Initial state radiation effects in inclusiveJ/ψproduction at B factories,JHEP 04(2014) 182 [1402.5840]

  19. [19]

    Lee and G

    K. Lee and G. Sterman,Power expansion for heavy quarkonium production at next-to-leading order in e+e− annihilation,JHEP09(2020) 046 [2006.07375]

  20. [20]

    Brambilla et al.,Heavy Quarkonium: Progress, Puzzles, and Opportunities,Eur

    N. Brambilla et al.,Heavy Quarkonium: Progress, Puzzles, and Opportunities,Eur. Phys. J. C71(2011) 1534 [1010.5827]

  21. [21]

    Chen, Y.-Q

    A.-P. Chen, Y.-Q. Ma and H. Zhang,A Short Theoretical Review of Charmonium Production,Adv. High Energy Phys.2022(2022) 7475923 [2109.04028]

  22. [22]

    Lansberg,New Observables in Inclusive Production of Quarkonia,Phys

    J.-P. Lansberg,New Observables in Inclusive Production of Quarkonia,Phys. Rept.889 (2020) 1 [1903.09185]

  23. [23]

    Kuhn and H

    J.H. Kuhn and H. Schneider,Testing QCD Through InclusiveJ/ψProduction ine +e− Annihilations,Z. Phys. C11(1981) 263

  24. [24]

    Kuhn and H

    J.H. Kuhn and H. Schneider,InclusiveJ/ψ ′s ine +e− Annihilations,Phys. Rev. D24(1981) 2996

  25. [25]

    Driesen, J.H

    V.M. Driesen, J.H. Kuhn and E. Mirkes,Testing J / psi production and decay properties in e+ e- annihilation,Phys. Rev. D49(1994) 3197

  26. [26]

    S. Baek, P. Ko, J. Lee and H.S. Song,Polarized J / psi production at CLEO,J. Korean Phys. Soc.33(1998) 97 [hep-ph/9804455]. – 19 –

  27. [27]

    Hagiwara, E

    K. Hagiwara, E. Kou, Z.H. Lin, C.F. Qiao and G.H. Zhu,InclusiveJ/ψproductions ate +e− colliders,Phys. Rev. D70(2004) 034013 [hep-ph/0401246]

  28. [28]

    Zhang, B.-Q

    Y.-J. Zhang, B.-Q. Li and K.-Y. Liu,J/ψelectromagnetic production associated with light hadrons atBfactories,1003.5566

  29. [29]

    He and J.-X

    Z.-G. He and J.-X. Wang,QED contribution toJ/ψproduction with light hadrons at B factories,Phys. Rev. D100(2019) 114025 [1301.4455]

  30. [30]

    Gong, Y.-D

    B. Gong, Y.-D. Wang and J.-X. Wang,Perturbative QCD forJ/ψinclusive production via initial state radiation ine+e− collisions,Chin. Phys. C43(2019) 083104 [1903.08161]

  31. [31]

    Wang,The Color octet effect frome +e− →J/ψ+X+γatBfactory,AIP Conf

    J.-X. Wang,The Color octet effect frome +e− →J/ψ+X+γatBfactory,AIP Conf. Proc.1092(2009) 1 [hep-ph/0311292]

  32. [32]

    Li, G.-Z

    Y.-J. Li, G.-Z. Xu, P.-P. Zhang, Y.-J. Zhang and K.-Y. Liu,Study of Color Octet Matrix Elements ThroughJ/ψProduction ine +e− Annihilation,Eur. Phys. J. C77(2017) 597 [1409.2293]

  33. [33]

    Chen, X.-B

    A.-P. Chen, X.-B. Jin, Y.-Q. Ma and C. Meng,Color-octet contributions for J/ψinclusive production at B factories in soft gluon factorization,JHEP03(2022) 202 [2201.04492]

  34. [34]

    F. Feng, Y. Jia and W.-L. Sang,Reconciling the nonrelativistic QCD prediction and the J/ψ→3γdata,Phys. Rev. D87(2013) 051501 [1210.6337]

  35. [35]

    W.-L. Sang, F. Feng and Y. Jia,O(α sv2)corrections to the hadronic decay of vector quarkonia,Phys. Rev. D102(2020) 094021 [2010.11896]

  36. [36]

    Bodwin, E

    G.T. Bodwin, E. Braaten and G.P. Lepage,Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium,Phys. Rev. D51(1995) 1125 [hep-ph/9407339]

  37. [37]

    Bodwin and A

    G.T. Bodwin and A. Petrelli,Order-v 4 corrections toS-wave quarkonium decay,Phys. Rev. D66(2002) 094011 [hep-ph/0205210]

  38. [38]

    Luke and M.J

    M.E. Luke and M.J. Savage,Power counting in dimensionally regularized NRQCD,Phys. Rev. D57(1998) 413 [hep-ph/9707313]

  39. [39]

    Bodwin, H.S

    G.T. Bodwin, H.S. Chung, J. Lee and C. Yu,Order-alpha(s) corrections to the quarkonium electromagnetic current at all orders in the heavy-quark velocity,Phys. Rev. D79(2009) 014007 [0807.2634]

  40. [40]

    H.-R. Dong, F. Feng and Y. Jia,O(α sv2)correction toe +e− →J/ψ+η c atBfactories, Phys. Rev. D85(2012) 114018 [1204.4128]

  41. [41]

    Li and J.-X

    X.-H. Li and J.-X. Wang,O(α sυ2)correction toJ/ψplusη c production ine +e− annihilation at √s=10.6 GeV,Chin. Phys. C38(2014) 043101 [1301.0376]

  42. [42]

    Xu, Y.-J

    G.-Z. Xu, Y.-J. Li, K.-Y. Liu and Y.-J. Zhang,α sv2 corrections toηc andχ cJ production recoiled with a photon ate+e− colliders,JHEP10(2014) 071 [1407.3783]

  43. [43]

    Gremm and A

    M. Gremm and A. Kapustin,Annihilation of S wave quarkonia and the measurement of alpha-s,Phys. Lett. B407(1997) 323 [hep-ph/9701353]

  44. [44]

    Hahn,Generating Feynman diagrams and amplitudes with FeynArts 3,Comput

    T. Hahn,Generating Feynman diagrams and amplitudes with FeynArts 3,Comput. Phys. Commun.140(2001) 418 [hep-ph/0012260]

  45. [45]

    Mertig, M

    R. Mertig, M. Bohm and A. Denner,FEYN CALC: Computer algebraic calculation of Feynman amplitudes,Comput. Phys. Commun.64(1991) 345. – 20 –

  46. [46]

    Feng and R

    F. Feng and R. Mertig,FormLink/FeynCalcFormLink : Embedding FORM in Mathematica and FeynCalc,1212.3522

  47. [47]

    Binosi, J

    D. Binosi, J. Collins, C. Kaufhold and L. Theussl,JaxoDraw: A Graphical user interface for drawing Feynman diagrams. Version 2.0 release notes,Comput. Phys. Commun.180(2009) 1709 [0811.4113]

  48. [48]

    X. Guan, X. Liu, Y.-Q. Ma and W.-H. Wu,Blade: A package for block-triangular form improved Feynman integrals decomposition,Comput. Phys. Commun.310(2025) 109538 [2405.14621]

  49. [49]

    Liu, Y.-Q

    X. Liu, Y.-Q. Ma and C.-Y. Wang,A Systematic and Efficient Method to Compute Multi-loop Master Integrals,Phys. Lett. B779(2018) 353 [1711.09572]

  50. [50]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma,AMFlow: A Mathematica package for Feynman integrals computation via auxiliary mass flow,Comput. Phys. Commun.283(2023) 108565 [2201.11669]

  51. [51]

    Liu and Y.-Q

    Z.-F. Liu and Y.-Q. Ma,Determining Feynman Integrals with Only Input from Linear Algebra,Phys. Rev. Lett.129(2022) 222001 [2201.11637]

  52. [52]

    Liu, Y.-Q

    X. Liu, Y.-Q. Ma, W. Tao and P. Zhang,Calculation of Feynman loop integration and phase-space integration via auxiliary mass flow,Chin. Phys. C45(2021) 013115 [2009.07987]

  53. [53]

    Li, X.-D

    C. Li, X.-D. Huang and W.-L. Sang,Two loop QCD corrections to e+e−→J/ψ+ηc in asymptotic expansion,Phys. Lett. B873(2026) 140150 [2506.16317]

  54. [54]

    Li, W.-L

    C. Li, W.-L. Sang and H.-F. Zhang,The next-to-next-to-leading-order QCD corrections to e+e−→η c/χcJ +γat B factories,JHEP05(2026) 163 [2512.04758]

  55. [55]

    H.-Y. Liu, C. Li and W.-L. Sang,Semi-analytical two-loop QCD corrections to e+e− →J/ψ+χ cJ at B factories,2606.02182

  56. [56]

    X. Chen, X. Guan, C.-Q. He, Y.-Q. Ma, J. Wang and D.-J. Zhang,Analytical two-loop amplitudes of e+e-→J/ψ+ηc at B factories,Phys. Rev. D113(2026) 074023 [2508.20777]

  57. [57]

    X. Chen, X. Guan, C.-Q. He and Y.-Q. Ma,Two-loop QCD corrections to e+e− →Z ∗ →J/ψ+J/ψ,Phys. Rev. D113(2026) 094017

  58. [58]

    Herren and M

    F. Herren and M. Steinhauser,Version 3 of RunDec and CRunDec,Comput. Phys. Commun.224(2018) 333 [1703.03751]

  59. [59]

    Bodwin, H.S

    G.T. Bodwin, H.S. Chung, D. Kang, J. Lee and C. Yu,Improved determination of color-singlet nonrelativistic QCD matrix elements for S-wave charmonium,Phys. Rev. D77 (2008) 094017 [0710.0994]. [64]Particle Data Groupcollaboration,Review of particle physics,Phys. Rev. D110(2024) 030001. – 21 –