We present a calculation of the perturbative quark-to-quark transverse parton distribution function at next-to-next-to-leading order based on a gauge invariant operator definition. We demonstrate for the first time that such a definition works beyond the first nontrivial order. We extract from our calculation the coefficient functions relevant for a next-to-next-to-next-to-leading logarithmic $Q_T$ resummation in a large class of processes at hadron colliders.

• Received 7 September 2012

DOI:https://doi.org/10.1103/PhysRevLett.109.242003