---
title: Collins–Soper Operator Definition
url: https://www.emergentmind.com/topics/collins-soper-operator-definition
type: topic
---

# Collins–Soper Operator Definition

The Collins–Soper operator definition establishes the nonperturbative kernel governing the rapidity evolution of transverse-momentum-dependent (TMD) distributions in QCD. Through a precise construction involving gauge-invariant Wilson-line operators, the Collins–Soper (CS) kernel isolates the universal, process-independent rapidity anomalous dimension that mediates TMD evolution. It is central to factorization theorems for processes such as Drell–Yan and SIDIS, connects perturbative and nonperturbative scales, and facilitates direct determination via lattice QCD. Below, the technical, operator-level foundation, practical lattice extraction, regulator schemes, and renormalization systematics are summarized based on explicit arXiv data.

## 1. Operator Definition of the Collins–Soper Kernel

The foundational object is the TMD soft function, a vacuum matrix element of Wilson lines forming a staple-shaped contour in coordinate space. For transverse separation $b_\perp$ and light-cone vectors $n$, $\bar n$ ($n^2 = \bar n^2 = 0$, $n \cdot \bar n = 2$), the unsubtracted soft function is
\[
S(b_\perp; \mu, \zeta) = \frac{1}{N_c} \langle 0 | U_{\bar n}^\dagger(\infty \to 0; 0)\, U_{n}(0 \to \infty; b_\perp)\, U_n^\dagger(\infty \to 0; b_\perp)\, U_{\bar n}(0 \to \infty; 0) |0\rangle_{\mu,\,\text{rap.\,reg.}(\zeta)},
\]
where the path-ordered Wilson lines $U_{n}$ and $U_{\bar n}$ extend to infinity along $n$ and $\bar n$, then connect transversely by $b_\perp$.

The Collins–Soper kernel $K(b_\perp; \mu)$ is defined as the rapidity (log-$\zeta$) derivative of the logarithm of the soft function:
\[
K(b_\perp; \mu) = \frac{d}{d\ln\zeta}\, \ln S(b_\perp; \mu, \zeta).
\]
This operator-level definition is fully gauge invariant and universal for all leading-twist TMD factorization applications [2306.06488, 2003.06063, 2003.02288].

## 2. Rapidity Regulation and Renormalization

Purely lightlike Wilson lines introduce rapidity (light-cone) divergences. Regulator schemes include tilting staple directions off the light cone by introducing $v^2 < 0$ or parameterizing the staple lengths $L$ (Euclidean) or using exponential damping factors $e^{-\delta |s|}$, where the regulator $\delta\to 0$ after renormalization [2306.06488, 2502.15033].

Ultraviolet logarithmic divergences from Wilson-line cusps and quark self-energies are removed by explicit renormalization constants $Z_\text{soft}(\mu)$, nonperturbatively determined in schemes such as RI/MOM, with perturbative matching to $\overline{\mathrm{MS}}$ [1910.08569, 2307.12359].

On the lattice, linear (self-energy) divergences are canceled by dividing by nonperturbative Wilson loops of matching geometry, e.g., $Z_E(2L, b_\perp)$ [2204.00200, 2509.26316].

## 3. Practical Extraction in Lattice QCD

To enable direct determination, the staple geometry is mapped to quasi-TMD operators in the LaMET approach. The ratios
\[
K(b_\perp; \mu) = \frac{1}{\ln(P_1^z/P_2^z)}\,\ln \frac{H(\dots,P_2^z)\,\tilde \Psi^+(x, b_\perp, \mu, P_1^z)}{H(\dots,P_1^z)\,\tilde \Psi^+(x, b_\perp, \mu, P_2^z)} + O((P^z)^{-2})
\]
is employed, where $\tilde\Psi^+$ denotes the (Euclidean) quasi-TMD wave function extracted at two large hadron momenta $P_1^z$ and $P_2^z$, and $H$ is a perturbative matching kernel [2306.06488, 2402.06725].

Complete nonperturbative cancellation of rapidity and UV divergences is achieved by taking ratios of fully renormalized operators. In the Coulomb-gauge-fixed framework, the operator reduces to a simple quark bilinear with no Wilson line, and the CS kernel is obtained from the $\ln P^z$ dependence of the matrix element ratio. Soft factor dependence cancels in these ratios, and the resulting kernel is process- and $x$-independent [2407.10739, 2504.04625].

## 4. Collateral Structures: Soft Functions and Matching

The soft function $S(b_\perp; \mu, \zeta)$ can be factorized into a rapidity-independent (intrinsic soft) function $S_I(b_\perp; \mu)$ and a rapidity-dependent (Collins–Soper) kernel part. In the lattice context, $S_I$ is isolated using ratios of high-momentum form factors to quasi-TMD wave functions:
\[
S_I(b_\perp; \mu, P_1, P_2) = \frac{F(b_\perp; P_1, P_2)}{\int dx_1\, dx_2\, H(x_1, x_2)\, \tilde{\Psi}^{\pm*}(x_2, b_\perp, P^z)\, \tilde{\Psi}^{\pm}(x_1, b_\perp, P^z)}
\]
where $F$ is a form factor normalized by appropriate local matrix elements and $H$ is the matching kernel [2306.06488, 2509.26316].

Leading-order and higher-order matching kernels for both quasi-TMDPDFs and quasi-TMDWFs are supplied at one-loop and beyond, constituting essential ingredients for the continuum extrapolation of lattice results [2509.26316, 2402.06725].

## 5. Evolution Equations and Renormalization-Group Consistency

The TMD distributions $F(x, b_\perp; \mu, \zeta)$ obey a pair of coupled evolution equations:
\[
\frac{d}{d\ln \mu}\,\ln F = \gamma_F(\mu, \zeta), \qquad \frac{d}{d\ln \zeta}\,\ln F = K(b_\perp, \mu)
\]
where $K$ is the Collins–Soper kernel. The RG consistency condition,
\[
\frac{d}{d\ln \mu} K(b_\perp, \mu) = -\Gamma_\mathrm{cusp}(\mu),
\]
enforces a nontrivial constraint linking the kernel to the universal cusp anomalous dimension, guaranteeing two-dimensional evolution consistency in the $(\mu, \zeta)$ plane [2003.02288, 2003.06063].

## 6. Distinctions in Implementation Schemes

Different implementations yield equivalent Collins–Soper kernels but offer distinct technical advantages:

- **Gauge-invariant staple-link (Euclidean/lattice) schemes**: Based on nonlocal quark bilinears connected by staple-shaped Wilson lines, regulated via staple length, with full nonperturbative renormalization [2307.12359, 1910.08569, 2103.16991].
- **Coulomb-gauge-fixed operators**: Avoid Wilson lines entirely, dramatically suppress noise and operator mixing, and are compatible with chiral-preserving discretizations. Ratios at different quark separations yield the CS kernel with minimal contamination [2407.10739, 2504.04625].
- **Off-light-cone factorization in DIS and Mellin moments**: The same operator structure, after soft-subtraction, allows extraction of $K$ as the rapidity derivative of a Wilson loop with tilted paths, applicable at threshold and beyond [2502.15033].

## 7. Systematic Uncertainties and Physical Implications

Comprehensive lattice studies have controlled for discretization $(a/b_\perp)$, operator mixing (full $16\times16$ Dirac mixing matrices), staple-length extrapolation, quark-mass dependence, and higher-twist contributions—especially via Fierz rearrangement of Dirac structures to project leading-twist channels [2306.06488, 2402.06725, 2307.12359, 2204.00200]. Systematic uncertainties from NNLL matching and nonperturbative renormalization matrix elements are included.

The operator-level definition of the Collins–Soper kernel establishes a process-independent, gauge-invariant anchor for the rapidity evolution of all leading-twist TMDs. First-principles QCD determinations of $K(b_\perp, \mu)$ at nonperturbative $b_\perp$ are now available, discriminating between phenomenological parameterizations and providing robust inputs to TMD evolution in collider phenomenology [2504.04625, 2407.10739, 2307.12359].

Source: https://www.emergentmind.com/topics/collins-soper-operator-definition