Quark Transverse Parton Distribution at the Next-to-Next-to-Next-to-Leading Order

We report a calculation of the perturbative matching coefficients for the transverse-momentum-dependent parton distribution functions for quark at the next-to-next-to-next-to-leading order in QCD, which involves calculation of nonstandard Feynman integrals with rapidity divergence. We introduce a set of generalized integration-by-parts equations, which allows an algorithmic evaluation of such integrals using the machinery of modern Feynman integral calculation.