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 revisit the calculation of perturbative quark transverse momentum dependent parton distribution functions and fragmentation functions using the exponential regulator for rapidity divergences. We show that the exponential regulator provides a consistent framework for the calculation of various ingredients in transverse momentum dependent factorization. Compared to existing regulators in the literature, the exponential regulator has a couple of advantages which we explain in detail. As a result, the calculation is greatly simplified and we are able to obtain the next-to-next-to-leading order results up to O( 2 ) in dimensional regularization. These terms are necessary for a higher order calculation which is made possible with the simplification brought by the new regulator. As a by-product, we have obtained the two-loop quark jet function for the Energy-Energy Correlator in the back-to-back limit, which is the last missing ingredient for its N 3 LL resummation.

Energy Correlators measure the energy deposited in multiple detectors as a function of the angles between the detectors. In this paper, we analytically compute the three particle correlator in the collinear limit in QCD for quark and gluon jets, and also in N = 4 super Yang-Mills theory. We find an intriguing duality between the integrals for the energy correlators and infrared finite Feynman parameter integrals, which maps the angles of the correlators to dual momentum variables. In N = 4, we use this duality to express our result as a rational sum of simple Feynman integrals (triangles and boxes). In QCD our result is expressed as a sum of the same transcendental functions, but with more complicated rational functions of cross ratio variables as coefficients. Our results represent the first analytic calculation of a three-prong jet substructure observable of phenomenological relevance for the LHC, revealing unexplored simplicity in the energy flow of QCD jets. They also provide valuable data for improving the understanding of the light-ray operator product expansion.