NLO QCD corrections to pseudoscalar quarkonium production with two heavy flavors in photon-photon collision

We calculate the next-to-leading order (NLO) quantum chromodynamics (QCD) corrections to $\gamma+\gamma\to \eta_c+c+\bar{c}$, $\gamma+\gamma\to \eta_b+b+\bar{b}$ and $\gamma+\gamma\to B_c+b+\bar{c}$ processes in the framework of non-relativistic QCD (NRQCD) factorization formalism. The cross sections at the SuperKEKB electron-positron collider, as well as the future collider like the Circular Electron Positron Collider (CEPC), are evaluated. Numerical results indicate that the NLO corrections are significant, and the uncertainties in theoretical predictions with NLO corrections are reduced as expected. Due to the high luminosity of the SuperKEKB collider, the $\eta_c+c+\bar{c}$ production is hopefully observable in the near future.