Fragmentation functions of polarized heavy quarkonium

Investigating the production of polarized heavy quarkonia in terms of recently proposed QCD factorization formalism requires the knowledge of a large number of input fragmentation functions (FFs) from a single parton or a heavy quark-antiquark pair to a polarized heavy quarkonium. We study these universal FFs at the input factorization scale mu(0) greater than or similar to 2m(Q), with heavy quark mass m(Q), in the framework of nonrelativistic QCD (NRQCD) factorization. We express these FFs in terms of perturbatively calculable coefficients for producing a heavy quark-antiquark pair in all possible NRQCD states, multiplied by corresponding NRQCD long-distance matrix elements for the pair to transmute into a polarized heavy quarkonium. We derive all relevant NRQCD operators for the long-distance matrix elements based on symmetries, and introduce a self-consistent scheme to define them in arbitrary d-dimensions. We compute, up to the first non-trivial order in alpha(s), the perturbative coefficients for producing a heavy quark pair in all possible S-wave and P-wave NRQCD states. We also discuss the role of the polarized FFs in generating QCD predictions for the polarization of J/produced at collider energies.