Color-singlet relativistic correction to inclusive $J/ψ$ production associated with light hadrons at $B$ factories
Yu Jia
TL;DR
The paper develops a novel NRQCD matching framework that uses physical hadronic kinematics and the Gremm-Kapustin relation to compute the first-order relativistic ($v^2$) correction to inclusive $J/\u001bpsi$ production in $e^+e^-$ annihilation at $B$ factories, focusing on $e^+e^-\to J/\u001bpsi+gg$ in the color-singlet channel. It provides analytic and numerical results for the integrated cross section and differential distributions (energy, polarization, and angular) at $\sqrt{s}=10.58$ GeV, finding a modest ~30% enhancement from the $O(v^2)$ term and only mild changes to distributions, with endpoint logarithms signaling the need for resummation or higher-order effects. The work also clarifies the relative merits of the proposed matching strategy, demonstrates consistency with orthodox approaches at the integrated level, and discusses implications for color-octet contributions and future refinements. Overall, it advances a practical method to incorporate relativistic corrections in complex quarkonium production kinematics and highlights theoretical uncertainties and endpoint issues that guide future improvements.
Abstract
We study the first-order relativistic correction to the associated production of $J/ψ$ with light hadrons at $B$ factory experiments at $\sqrt{s}=10.58$ GeV, in the context of NRQCD factorization. We employ a strategy for NRQCD expansion that slightly deviates from the orthodox doctrine, in that the matching coefficients are not truly of ``short-distance" nature, but explicitly depend upon physical kinematic variables rather than partonic ones. Our matching method, with validity guaranteed by the Gremm-Kapustin relation, is particularly suited for the inclusive quarkonium production and decay processes with involved kinematics, exemplified by the process $e^+e^-\to J/ψ+gg$ considered in this work. Despite some intrinsic ambiguity affiliated with the order-$v^2$ NRQCD matrix element, if we choose its value as what has been extracted from a recent Cornell-potential-model-based analysis, including the relative order-$v^2$ effect is found to increase the lowest-order prediction for the integrated $J/ψ$ cross section by about 30\%, and exert a modest impact on $J/ψ$ energy, angular and polarization distributions except near the very upper end of the $J/ψ$ energy. The order-$v^2$ contribution to the energy spectrum becomes logarithmically divergent at the maximum of $J/ψ$ energy. A consistent analysis may require that these large end-point logarithms be resummed to all orders in $α_s$.
