Within the framework of nonrelativistic QCD (NRQCD) factorization, we compute the three-loop QCD corrections to the decay constant of $B_c$. We reconstruct the analytical expressions for the three-loop renormalization constant and the corresponding anomalous dimension affiliated with the pseudoscalar current composed of two different heavy flavors in NRQCD, which are functions of the ratio between the charm and bottom quark masses. Meanwhile, the short-distance coefficient is obtained with very high numerical accuracy. The three-loop QCD correction turns out to be overwhelmingly large. The phenomenological implication of this new piece of radiative corrections for the $B_c$ leptonic decay is also addressed.
Comment: 9 pages, 4 figures and 2 tables