Abstract:The failure of observing the $e^+ e^- \to J/\psi+J/\psi$ events at $B$ factories to date is often attributed to the significant negative order-$\alpha_s$ correction. In this work we compute the ${\cal O}(\alpha^2_s)$ 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 NRQCD 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/\psi$ decay constant. We thus implement an improved NRQCD factorization framework, by decomposing the amplitude into the photon-fragmentation piece and the non-fragmentation piece. With the measured $J/\psi$ 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 non-fragmentation type of the amplitude is then computed through NNLO in $\alpha_s$ and at lowest order in velocity. Both the ${\cal O}(\alpha_s)$ and ${\cal O}(\alpha^2_s)$ corrections in the interference term become positive and exhibit a decent convergence behavior. Our finest prediction is $\sigma(e^+ e^- \to J/\psi+J/\psi)= 2.13^{+0.30}_{-0.06}$ fb at $\sqrt{s}=10.58$ GeV. With the projected integrated luminosity of 50 ${\rm ab}^{-1}$, the prospect to observe this exclusive process at \texttt{Belle} 2 experiment appears to be bright.

Calculation of Electroproduction to NNLO and Precision Determination of $α_s$

Mixed NNLO QCD $\times$ electroweak corrections to single-Z production in pole approximation: differential distributions and forward-backward asymmetry

Determination of the Strong Coupling \boldmath{\as} from hadronic Event Shapes and NNLO QCD predictions using JADE Data

Calculations for deep inelastic scattering using fast interpolation grid techniques at NNLO in QCD and the extraction of $α_s$ from HERA data

A first determination of the strong coupling $α_S$ at approximate N$^{3}$LO order in a global PDF fit

$W^\pm Z$ production at NNLO QCD and NLO EW matched to parton showers with MiNNLO$_{\rm PS}$

Two loop operator matrix elements calculated up to finite terms for polarized deep inelastic lepton - hadron scattering

NLO QCD + NLO EW Corrections to WZZ Productions with Leptonic Decays at the LHC

QCD analysis of valence structure functions using deep inelastic lepton-nucleon scattering

The NN scattering 3S1-3D1 mixing angle at NNLO

Optimized ${\cal O}(α_s^2)$ correction to exclusive double $J/ψ$ production at $B$ factories

Complete Next-to-Leading Order Calculation of Single Inclusive $π^0$ Production in Forward Proton-Nucleus Collisions

A complete NLO calculation of the $J/ψ$ and $ψ'$ production at hadron colliders

NLO corrections to triple vector-boson production in final states with three charged leptons and two jets

Next-to-next-to-leading Order QCD⊗EW Corrections to Z-boson Pair Production at Electron-Positron Colliders

Mixed NNLO QCD × electroweak corrections to single-Z production in pole approximation: differential distributions and forward-backward asymmetry

Matching N3LO QCD calculations to parton showers

Next-to-next-to-leading order $\text{QCD} \otimes \text{EW}$ corrections to $Z$-boson pair production at electron-positron colliders

N$^{\mathbf{3}}$LL + $\mathcal{O}(α_s^2)$ predictions of lepton-jet azimuthal angular distribution in deep-inelastic scattering

s-channel single top quark production and decay at next-to-next-to-leading-order in QCD