We systematically calculate the relativistic corrections to the polarization
variables of prompt J/ψ photoproduction and hadroproduction using the
factorization formalism of nonrelativistic QCD. Specifically, we include the
3S1[1] and 3PJ[1] color-singlet and the
3S1[8], 1S0[8], and 3PJ[8] color-octet
channels as well as the effects due to the mixing between the 3S1[8]
and 3D1[8] channels. We provide all the squared hard-scattering
amplitudes in analytic form. Assuming the nonrelativistic-QCD long-distance
matrix elements to satisfy the velocity scaling rules, we find the relativistic
corrections to be appreciable, especially at small transverse momentum pT
and large inelasticity z. The results obtained here and in our previous work
on the unpolarized yield [Phys. Rev. D 90, 014045 (2014)] will help to render
global analyses of prompt J/ψ production data more complete and hopefully
to shed light on the J/ψ polarization puzzle.Comment: 42 pages, 4 figures, matches journal version, PRD Editors' Suggestio