Gluon TMD fragmentation function into quarkonium

We compute the gluon transverse-momentum-dependent fragmentation function (TMDFF) at next-to-leading order (NLO) into heavy quarkonium in the color-octet $^3S_1^{[8]}$ channel, based on the NRQCD factorization approach. The spurious rapidity divergences are explicitly shown to cancel in a well-defined TMDFF, which incorporates the needed soft factor. We also compute the integrated gluon FF at NLO in the same $^3S_1^{[8]}$ channel, and show that the matching coefficient of the TMDFF onto the FF at large transverse momentum is the expected one. These results are relevant to perform precise and sensible phenomenological studies of transverse-momentum spectra of quarkonium production, for which the production mechanism through fragmentation plays a relevant role, like in the future Electron-Ion Collider.