Order alpha_s³ ln²(1/alpha_s) Corrections to Heavy-Quarkonium Creation and Annihilation

In the framework of nonrelativistic QCD, we compute the leading double-logarithmic corrections of order alpha_s^3 ln^2(1/alpha_s) to the heavy-quark-antiquark bound-state wave function at the origin, which determines the production and annihilation rates of heavy quarkonia. The phenomenological implications for the top-antitop and Upsilon systems are discussed.