Within NRQCD factorization framework, in this work we compute, at the lowest order in velocity expansion, the next-to-leading-order (NLO) perturbative corrections to the short-distance coefficients associated with heavy quark fragmentation into the $${}^1S_0^{(1,8)}$$
1
S
0
(
1
,
8
)
components of a heavy quarkonium. Starting from the Collins and Soper’s operator definition of the quark fragmentation function, we apply the sector decomposition method to facilitate the numerical manipulation. It is found that the NLO QCD corrections have a significant impact.