Semi-analytical two-loop QCD corrections to $e^{+}{e}^{-}\to J/\ddot{I}ˆ+\ddot{I}{‡}_{cJ}$ at B factories

Abstract

In this work, we compute the next-to-next-to-leading-order (NNLO) QCD corrections to the process $e^{+}{e}^{-}\to J/\ddot{I}ˆ+\ddot{I}{‡}_{cJ}$ at B factories within the NRQCD factorization framework. The helicity amplitudes are obtained via asymptotic expansions around $r=0$ and $r=1$, with $r=16m_c^2/s$. Our asymptotic expressions reproduce the exact numerical results with high accuracy across the entire range $0\le r\le 1$, achieving a relative error below $10^{-5}$, which is sufficient for phenomenological applications. Notably, the large logarithmic terms are obtained analytically. We compute the unpolarized cross sections. The $\mathcal{O}(\hat{I}{\pm }_s)$ correction is found to be large, while the $\mathcal{O}(\hat{I}{\pm }_s^2)$ correction for $\ddot{I}{‡}_{c0}$ production amounts to $33\%$ of the leading-order (LO) cross section, significantly reducing the scale uncertainties. For $\ddot{I}{‡}_{c1}$, the $\mathcal{O}(\hat{I}{\pm }_s)$ and $\mathcal{O}(\hat{I}{\pm }_s^2)$ corrections correspond to $35\%$ and $-15\%$, respectively. For $\ddot{I}{‡}_{c2}$, the corresponding corrections are $25\%$ and $-38\%$. The large cancellation between the corrections for $\ddot{I}{‡}_{c2}$ brings the NNLO cross section close to the LO prediction. Our prediction for $\ddot{I}{‡}_{c0}$ is consistent with the {\tt Belle} measurement and agrees with the {\tt BaBar} data within $2\ddot{I}ƒ$. We also predict the angular distribution parameters $\hat{I}{\pm }_{\hat{I}}^J¸$, which are independent of nonperturbative inputs. A sharp discrepancy between the theory and the {\tt Belle} measurement is observed for $\hat{I}{\pm }_{\hat{I}}^0¸$, calling for further experimental and theoretical investigations. Moreover, future measurements of the angular distribution parameters for $\ddot{I}{‡}_{c1}$ and $\ddot{I}{‡}_{c2}$ will provide important tests of the theoretical framework.

30 pages, 5 figures, 2 tables