---
title: Transverse Momentum Distributions in QCD
url: https://www.emergentmind.com/topics/transverse-momentum-distributions
type: topic
---

# Transverse Momentum Distributions in QCD

Transverse momentum distributions (TMDs) are quantitative descriptions of the probability density for finding a parton, hadron, or other process-specific particle with a particular value of transverse momentum relative to a well-defined high-energy axis of reference within a hadron or jet. TMDs are central to the theoretical and phenomenological understanding of multidimensional partonic structure, QCD factorization, nonperturbative dynamics, and the resummation of soft and collinear radiation. They underpin a wide range of processes, from low-mass Drell–Yan, Higgs and electroweak boson $p_T$ spectra, through semi-inclusive deep-inelastic scattering (SIDIS) and $e^+e^-$ annihilation, to jet substructure and hadron fragmentation. TMDs are operator-defined correlators, admitting systematic factorization theorems, perturbative expansions, all-order resummation, and nonperturbative modeling.

## 1. Operator Definitions, Factorization, and Evolution

A TMD is defined as a QCD matrix element with explicit dependence on transverse spatial separation or transverse momentum. For quarks in a hadron $h$, the unpolarized TMD PDF in impact-parameter space is
\[
F_{q/h}(x,\mathbf{b};\mu,\zeta) = \int \frac{d\lambda}{2\pi} e^{-ix\lambda p^+}
\langle p| \bar q(\lambda n + \mathbf{b}) W[\lambda n + \mathbf{b},\infty]\Gamma W[\infty,0]q(0) |p\rangle_{(\mu,\zeta)}
\]
where $x$ is the longitudinal momentum fraction, $\mathbf{b}$ is the transverse separation (Fourier-conjugate to $k_T$), $\mu$ and $\zeta$ are UV and rapidity renormalization scales, $W$ are Wilson lines, and $\Gamma$ denotes Dirac structure (e.g., $\gamma^+$ for $f_1$, $\gamma^+\gamma_5$ for $g_1$) [1209.3892, 2205.05997, 2402.01836, 2501.17274].

TMDs enter factorization theorems for cross sections differential in transverse momentum, as in Drell–Yan,
\[
\frac{d\sigma}{dQ^2 dY d^2q_T} = \sum_{ij} \sigma_{ij}^{(0)}(Q)
\int \frac{d^2 \mathbf{b}}{(2\pi)^2} e^{i\mathbf{b}\cdot \mathbf{q}_T}
F_{1,i\leftarrow h_1}(x_1,b;\mu,\zeta_1)\, F_{1,j\leftarrow h_2}(x_2,b;\mu,\zeta_2)
\]
with $Q$ the hard scale, $Y$ rapidity, $q_T$ the observed transverse momentum, and a mandate that $\zeta_1 \zeta_2 = Q^4$ [2205.05997, 1108.3609].

The Collins–Soper (CS) and renormalization group (RG) equations for TMDs are
\[
\frac{\partial}{\partial\ln\sqrt{\zeta}} F(x,b;\mu,\zeta) = -\mathcal{D}(b;\mu)F(x,b;\mu,\zeta),\qquad
\frac{\partial}{\partial\ln\mu} F(x,b;\mu,\zeta) = \gamma_F(\mu,\zeta) F(x,b;\mu,\zeta)
\]
$\mathcal{D}$ is the rapidity anomalous dimension; $\gamma_F$ involves the cusp anomalous dimension and noncusp terms. Solutions exponentiate the scale dependence, with scales sometimes fixed by the $\zeta$-prescription to minimize spurious logarithms [1209.3892, 2402.01836].

At small $b$, TMDs match onto collinear PDFs via a perturbative operator product expansion (OPE):
\[
F_{i/h}(x,b;\mu,\zeta) = \sum_j \int_x^1 \frac{dy}{y}\, C_{i\leftarrow j}(y,b;\mu,\zeta)\, f_{j/h}\left(\frac{x}{y};\mu\right) + \mathcal O(b^2)
\]
The Wilson coefficients $C_{i\leftarrow j}$ contain singular threshold terms at $y\to1$ which require special all-orders treatment at large $x$ [2501.17274].

## 2. Perturbative Expansion, Resummation, and Matching

At fixed order in $\alpha_s$, the $p_T$ spectrum for processes such as Drell–Yan or electroweak boson production is computed as
\[
\frac{d^3 \sigma}{dM^2\, dy\, dp_T^2} = \sum_{a,b} \int dx_a dx_b \, f_{a/A}(x_a, \mu_F) f_{b/B}(x_b, \mu_F) \frac{d\hat{\sigma}_{ab}}{dM^2\,dy\, dp_T^2}
\]
where $d\hat{\sigma}_{ab}$ includes real and virtual corrections up to NNLO and $f_{i/H}$ are collinear PDFs [2110.15839]. The first non-vanishing $p_T$ contributions arise from real emission beyond Born level.

For $p_T \ll Q$, logarithmically enhanced terms $\alpha_s^n \ln^m(Q^2/p_T^2)$ appear, with $m\leq 2n-1$. All-order resummation is required, performed either in $b$-space (CSS, SCET) or directly in momentum space. Resummation is achieved via exponentiation of the Sudakov form factor and matching to fixed-order for $p_T \sim Q$ [1108.3609, 1007.3773]. Resummation accuracy attains N$^3$LL in current state-of-the-art analyses [2501.17274].

A summary of benchmark results for low-mass Drell–Yan $p_T$ spectra at RHIC (PHENIX) and fixed-target (NuSea) is:

| Regime         | Order | K-factor     | Scale Unc. | Data/Theory | 
| -------------- | ----- | ----------- | -----------| --------------|
| RHIC (PHENIX)  | LO    | –           | ±50%       | $\sim$2–3 low | 
|                | NLO   | +30%        | ±25%       | ~1 for $p_T>1.5$ GeV |
|                | NNLO  | +25%        | ±15%       | within errors |
| NuSea (E866)   | NNLO  | +35–70%     | ±35–45%    | $\sim$1.5–10 |

For $p_T \lesssim 1.5$ GeV, fixed-order diverges and explicit resummation is mandatory [2110.15839].

## 3. Nonperturbative Corrections and Intrinsic $k_T$

Comparison with experimental data, particularly in the fixed-target regime and at high Feynman-$x_F$, reveals that even NNLO perturbative calculations severely underestimate the observed $p_T$ spectra. The discrepancy is nearly $p_T$-independent and can be largely resolved by introducing a nonperturbative intrinsic momentum shift,
\[
\Delta p_T \simeq 1~\text{GeV}
\]
which amounts to Gaussian smearing in the incoming partons' $k_T$,
\[
F_{\rm non\text{-}pert}(k_T) \propto \exp[-k_T^2/\langle k_T^2 \rangle], \qquad \sqrt{\langle k_T^2 \rangle} \in [0.5,1.0]~\text{GeV}
\]
Earlier TMD fits confirm that nonperturbative input is essential for a quantitative description at low $p_T$, especially in proton and pion-induced Drell–Yan [2110.15839, 1709.03324, 2210.01733].

In the TMD framework, the nonperturbative sector is parametrized as
\[
f^{\text{NP}}_{f\leftarrow h}(x,b) = \exp\left(-\frac{(1-x)\lambda_1^f + x\lambda_2^f}{b^2\sqrt{1+\lambda_0 x^2 b^2}}\right)
\]
or similar variants, whose parameters are fitted globally. PDF uncertainties and "PDF bias" due to imperfect knowledge of collinear distributions remain dominant theoretical errors at high precision [2205.05997].

## 4. Large-$x$ Structure and Threshold Resummation

The $x\to 1$ limit probes the endpoint of the parton distribution and is dominated by soft and collinear emissions. The matching coefficient $C_{q\leftarrow q}(x,b)$ has an all-orders exponentiation structure,
\[
C_{q\leftarrow q}(x,b) \sim (1-x)^{\alpha_s A_1} [1+\mathcal O(1-x)],
\]
where $A_1$ is set by the cusp anomalous dimension. This leads to a resummed form for the coefficient functions,
\[
\Sigma(b,\mu;x) = \exp[L g_1(\alpha_s L) + g_2(\alpha_s L) + \ldots], \quad L = \ln(1-x)
\]
The resummation is process-independent for all TMDs matching onto leading-twist, subject only to the cusp and soft anomalous dimensions. This structure improves perturbative convergence and restricts nonperturbative model forms by constraining the behavior of $b^*(b)$ prescriptions to avoid non-integrable endpoint singularities [2501.17274].

At NNLO and N$^3$LL, this resummation reduces the fixed-order dependence in the threshold region to the percent level, enabling robust extrapolation and theoretical uncertainty control.

## 5. Phenomenological Extraction and Universality

Dedicated global analyses using SIDIS, Drell–Yan, $Z$-boson production, and jet substructure provide simultaneous fits of unpolarized TMD PDFs, TMD FFs, and pion TMDs [1709.03324, 2210.01733, 1705.08443]. The intrinsic widths extracted typically satisfy $\langle k_T^2\rangle \sim 0.3$–$0.6~\text{GeV}^2$ at $Q_0 = 1$ GeV for protons and are broader for pions ($\sim$0.47 GeV$^2$). Flavour dependence is found to be moderate but non-negligible, and model uncertainties remain substantial at small $x$ and for gluons.

Jet-based observables and groomed-jet measurements admit factorization theorems where TMD PDFs and FFs can be accessed in cleaner configurations, allowing for more direct probes of gluon TMDs and fragmentation at the LHC and EIC [1705.08443, 1907.05896].

Thermal models using the Tsallis distribution also provide excellent fits to $p_T$ spectra over several orders of magnitude at LHC and RHIC energies, consistently yielding a non-extensivity parameter $q=1.1$–$1.2$ and effective temperatures $T=0.07$–$0.1$ GeV [1203.4343, 1310.0217]. While phenomenologically successful, these statistical approaches do not capture the QCD factorization and evolution encoded in field-theoretic TMDs.

## 6. Future Challenges and Outlook

Current limitations and future directions include:

- The need for combined NNLO (or higher) fixed-order calculations, N$^3$LL resummation, and flexible nonperturbative models ("NNLO+N$^3$LL+np") to describe the full $p_T$ range with precision [2110.15839, 2501.17274].
- Inclusion of flavor-dependent nonperturbative functions to mitigate PDF bias and understand universality breaking between protons, pions, and nuclei [2205.05997, 2210.01733].
- Extraction and evolution of polarized TMDs, especially the helicity and orbital angular momentum content, and their phenomenological impact [2409.08110].
- Explicit operator definitions for heavy-quark TMDs, spin-asymmetry sum rules, and the behavior of TMDs in the heavy quark limit [2305.15461].
- Lattice QCD calculations of TMD moments and their conversion to standard schemes using perturbative matching at up to three loops [2402.01836].
- Ongoing and forthcoming experimental input from Jefferson Lab 12 GeV, COMPASS, LHC, and the future EIC, aiming to reduce uncertainties, constrain the $k_T$-dependence, and test theoretical frameworks across flavors, scales, and processes.

The unification of higher-order resummation, robust nonperturbative modeling, and precise global data analysis is the central challenge for the next decade of TMD phenomenology and theory.

Source: https://www.emergentmind.com/topics/transverse-momentum-distributions