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
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.
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
- 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.
Referee Report
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)
- [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
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
-
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
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
free parameters (2)
- renormalization scale μ_R =
√s/2
- long-distance matrix elements
axioms (1)
- domain assumption NRQCD factorization remains valid at O(α_s v²)
Reference graph
Works this paper leans on
-
[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
1981
-
[2]
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]
Pith/arXiv arXiv 1994
-
[3]
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]
Pith/arXiv arXiv 1996
-
[4]
F. Yuan, C.-F. Qiao and K.-T. Chao,PromptJ/ψproduction ate +e− colliders,Phys. Rev. D56(1997) 321 [hep-ph/9703438]
Pith/arXiv arXiv 1997
-
[5]
P.L. Cho and A.K. Leibovich,Color singlet psi(Q) production at e+ e- colliders,Phys. Rev. D54(1996) 6690 [hep-ph/9606229]
Pith/arXiv arXiv 1996
-
[6]
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]
Pith/arXiv arXiv 1997
-
[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 –
Pith/arXiv arXiv 1999
-
[8]
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]
Pith/arXiv arXiv 2003
-
[9]
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]
Pith/arXiv arXiv 2004
-
[10]
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]
Pith/arXiv arXiv 2003
-
[11]
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]
Pith/arXiv arXiv 2007
-
[12]
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]
Pith/arXiv arXiv 2009
-
[13]
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]
Pith/arXiv arXiv 2009
-
[14]
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]
Pith/arXiv arXiv 2009
-
[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]
Pith/arXiv arXiv 2010
-
[16]
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]
Pith/arXiv arXiv 2010
-
[17]
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]
Pith/arXiv arXiv 2010
-
[18]
H.-S. Shao,Initial state radiation effects in inclusiveJ/ψproduction at B factories,JHEP 04(2014) 182 [1402.5840]
Pith/arXiv arXiv 2014
- [19]
-
[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]
Pith/arXiv arXiv 2011
-
[21]
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]
arXiv 2022
-
[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]
arXiv 2020
-
[23]
Kuhn and H
J.H. Kuhn and H. Schneider,Testing QCD Through InclusiveJ/ψProduction ine +e− Annihilations,Z. Phys. C11(1981) 263
1981
-
[24]
Kuhn and H
J.H. Kuhn and H. Schneider,InclusiveJ/ψ ′s ine +e− Annihilations,Phys. Rev. D24(1981) 2996
1981
-
[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
1994
-
[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 –
Pith/arXiv arXiv 1998
-
[27]
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]
Pith/arXiv arXiv 2004
-
[28]
Y.-J. Zhang, B.-Q. Li and K.-Y. Liu,J/ψelectromagnetic production associated with light hadrons atBfactories,1003.5566
-
[29]
Z.-G. He and J.-X. Wang,QED contribution toJ/ψproduction with light hadrons at B factories,Phys. Rev. D100(2019) 114025 [1301.4455]
arXiv 2019
-
[30]
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]
arXiv 2019
-
[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]
Pith/arXiv arXiv 2009
-
[32]
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]
Pith/arXiv arXiv 2017
-
[33]
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]
arXiv 2022
-
[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]
Pith/arXiv arXiv 2013
-
[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]
arXiv 2020
-
[36]
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]
Pith/arXiv arXiv 1995
-
[37]
G.T. Bodwin and A. Petrelli,Order-v 4 corrections toS-wave quarkonium decay,Phys. Rev. D66(2002) 094011 [hep-ph/0205210]
Pith/arXiv arXiv 2002
-
[38]
M.E. Luke and M.J. Savage,Power counting in dimensionally regularized NRQCD,Phys. Rev. D57(1998) 413 [hep-ph/9707313]
Pith/arXiv arXiv 1998
-
[39]
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]
Pith/arXiv arXiv 2009
-
[40]
H.-R. Dong, F. Feng and Y. Jia,O(α sv2)correction toe +e− →J/ψ+η c atBfactories, Phys. Rev. D85(2012) 114018 [1204.4128]
Pith/arXiv arXiv 2012
-
[41]
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]
Pith/arXiv arXiv 2014
-
[42]
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]
Pith/arXiv arXiv 2014
-
[43]
M. Gremm and A. Kapustin,Annihilation of S wave quarkonia and the measurement of alpha-s,Phys. Lett. B407(1997) 323 [hep-ph/9701353]
Pith/arXiv arXiv 1997
-
[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]
Pith/arXiv arXiv 2001
-
[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 –
1991
-
[46]
F. Feng and R. Mertig,FormLink/FeynCalcFormLink : Embedding FORM in Mathematica and FeynCalc,1212.3522
-
[47]
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]
Pith/arXiv arXiv 2009
-
[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]
arXiv 2025
-
[49]
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]
Pith/arXiv arXiv 2018
-
[50]
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]
arXiv 2023
-
[51]
Z.-F. Liu and Y.-Q. Ma,Determining Feynman Integrals with Only Input from Linear Algebra,Phys. Rev. Lett.129(2022) 222001 [2201.11637]
arXiv 2022
- [52]
- [53]
- [54]
-
[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]
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]
arXiv 2026
-
[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
2026
-
[58]
F. Herren and M. Steinhauser,Version 3 of RunDec and CRunDec,Comput. Phys. Commun.224(2018) 333 [1703.03751]
Pith/arXiv arXiv 2018
-
[59]
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 –
Pith/arXiv arXiv 2008
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.