We report a calculation of the perturbative matching coefficients for the transverse-momentumdependent parton distribution functions for quark at the next-to-next-to-next-to-leading order in QCD, which involves calculation of non-standard 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.

We calculate in this paper the perturbative gluon transverse momentum dependent parton distribution functions (TMDPDFs) and fragmentation functions (TMDFFs) using the exponential regulator for rapidity divergences. We obtain results for both unpolarized and linearly polarized distributions through next-to-next-to leading order in strong coupling constant, and through O( 2 ) in dimensional regulator (finding discrepancy for the linearly polarized gluon TMDPDFs with a previous result in the literature). We find a nontrivial momentum conservation sum rule for the linearly polarized component for both TMDPDFs and TMDFFs in the N = 1 super-Yang-Mills theory. The TMDFFs are used to calculate the two-loop gluon jet function for the energy-energy correlator in Higgs gluonic decay in the back-to-back limit.

In this paper we calculate analytically the perturbative matching coefficients for unpolarized quark and gluon Transverse-Momentum-Dependent (TMD) Parton Distribution Functions (PDFs) and Fragmentation Functions (FFs) through Next-to-Next-to-Next-to-Leading Order (N3LO) in QCD. The N3LO TMD PDFs are calculated by solving a system of differential equation of Feynman and phase space integrals. The TMD FFs are obtained by analytic continuation from space-like quantities to time-like quantities, taking into account the probability interpretation of TMD PDFs and FFs properly. The coefficient functions for TMD FFs exhibit double logarithmic enhancement at small momentum fraction z. We resum such logarithmic terms to the third order in the expansion of αs. Our results constitute important ingredients for precision determination of TMD PDFs and FFs in current and future experiments.