Optimized O(α_{s}^{2}) Correction to Exclusive Double-J/ψ Production at B Factories

The failure of observing the e^{+}e^{-}→J/ψJ/ψ events at B factories to date is often attributed to the significant negative order-α_{s} correction. In this work we compute the O(α_{s}^{2}) correction to this process for the first time. The magnitude of the next-to-next-to-leading order (NNLO) perturbative correction is substantially negative so that the standard nonrelativistic QCD prediction would suffer from an unphysical, negative cross section. This dilemma may be traced in the fact that the bulk contribution of the fixed-order radiative corrections stems from the perturbative corrections to the J/ψ decay constant. We thus implement an improved nonrelativistic QCD factorization framework, by decomposing the amplitude into the photon-fragmentation piece and the nonfragmentation piece. With the measured J/ψ decay constant as input, which amounts to resumming a specific class of radiative and relativistic corrections to all orders, the fragmentation-induced production rate can be predicted accurately and serves a benchmark prediction. The nonfragmentation type of the amplitude is then computed through NNLO in α_{s} and at lowest order in velocity. Both the O(α_{s}) and O(α_{s}^{2}) corrections in the interference term become positive and exhibit a decent convergence behavior. Our finest prediction is σ(e^{+}e^{-}→J/ψJ/ψ)=2.13_{-0.06}^{+0.30}  fb at sqrt[s]=10.58  GeV. With the projected integrated luminosity of 50  ab^{-1}, the prospect to observe this exclusive process at Belle 2 experiment appears to be bright.