DGLAP Evolution and Scaling Violation

DGLAP evolution predicts how renormalized parton distributions change with the factorization scale. It resums collinear logarithms, turns naive Bjorken scaling into calculable logarithmic scaling violation, and preserves flavor number and total momentum when kernels, indices, thresholds, and distributions are implemented consistently.

Required background. Collinear factorization and operator-defined PDFs supplies the distributions, convolution, and scheme cancellation. Coefficient evolution supplies the compensating evolution of short-distance coefficients.

Helpful background. Mellin transforms and scaling asymptotics supplies the moment-space diagonalization used below.

Spacelike evolution in momentum fraction

For a PDF vector $f_i(x,\mu)$,

$$\frac{\partial f_i(x,\mu)}{\partial\ln\mu^2} =\sum_j\left[P_{ij}(\alpha_s(\mu))\otimes f_j(\mu)\right](x),$$ $$[P\otimes f](x)=\int_x^1\frac{dz}{z}\,P(z)f\!\left(\frac{x}{z}\right), \qquad P_{ij}=\frac{\alpha_s}{2\pi}P_{ij}^{(0)}+ \left(\frac{\alpha_s}{2\pi}\right)^2P_{ij}^{(1)}+\cdots.$$

$P_{ij}$ describes parton $j$ feeding parton $i$ in this index convention. Reversing that convention transposes the singlet matrix, so the convention must travel with any numerical kernel. The original leading-logarithmic equations and kernels are given by Altarelli and Parisi 1977, §§2–4, pp. 298–318.

For a quark nonsinglet combination, gluon mixing cancels. One useful leading-order form is

$$P_{\mathrm{NS}}^{(0)}(z)=C_F \left[\frac{2}{(1-z)_+}-1-z+\frac32\delta(1-z)\right].$$

The plus distribution is defined by

$$\int_0^1dz\,\frac{\phi(z)}{(1-z)_+} =\int_0^1dz\,\frac{\phi(z)-\phi(1)}{1-z}.$$

It represents the real-emission singularity together with the endpoint subtraction needed for a finite distribution. The $\delta(1-z)$ term contains the virtual contribution. Dropping either piece breaks conserved moments even if the pointwise kernel looks plausible.

Moment-space derivation and exact checks

Define Mellin moments

$$f_N(\mu)=\int_0^1dx\,x^{N-1}f(x,\mu), \qquad \gamma_N=\int_0^1dz\,z^{N-1}P^{(0)}(z).$$

Changing variables $x=zy$ in the convolution gives

$$\int_0^1dx\,x^{N-1}[P\otimes f](x)=\gamma_Nf_N.$$

The integro-differential equation therefore becomes an ordinary one. For a fixed coupling and $t=\ln(\mu^2/\mu_0^2)$,

$$f_N(t)=\exp\!\left[\frac{\alpha_s}{2\pi}\gamma_Nt\right]f_N(0).$$

With the one-loop running coupling in a fixed-$n_f$ interval,

$$f_N(\mu)=f_N(\mu_0) \left[\frac{\alpha_s(\mu)}{\alpha_s(\mu_0)}\right]^{-2\gamma_N/\beta_0}.$$

For the nonsinglet kernel above,

$$\gamma_1=C_F\left[0-\frac32+\frac32\right]=0,$$

so the net nonsinglet quark number is scale independent. The next moment is an especially useful sign check:

$$\gamma_2=C_F\left[-2-\frac56+\frac32\right] =-\frac{16}{9}.$$

Thus this quark moment decreases as the resolution scale grows: radiation redistributes longitudinal momentum toward smaller $x$.

Singlet–gluon mixing and momentum conservation

Let $\Sigma=\sum_f(q_f+\bar q_f)$. Its evolution mixes with the gluon:

$$\frac{d}{d\ln\mu^2} \begin{pmatrix}\Sigma\\ g\end{pmatrix} =\frac{\alpha_s}{2\pi} \begin{pmatrix}P_{qq}&P_{qg}\\P_{gq}&P_{gg}\end{pmatrix} \otimes \begin{pmatrix}\Sigma\\ g\end{pmatrix}.$$

At moment $N=2$, momentum conservation requires the column sums of the moment matrix to vanish. For $n_f=3$ in the leading-order normalization used here, this check is

$$\Gamma_2= \begin{pmatrix} -16/9&1\\[2pt] 16/9&-1 \end{pmatrix}, \qquad (1,1)\Gamma_2=(0,0).$$

Therefore $\Sigma_2+g_2$ remains constant under this evolution. This condition is stronger than checking either component independently and catches an accidental transpose immediately.

The physical pattern is a redistribution, not creation of hadron momentum. Higher resolution resolves more gluons and sea quarks at small $x$, while valence-like large-$x$ moments tend to decrease. The derivation from collinear emission and its application to DIS scaling violation are developed in Schwartz 2014, §32.2, pp. 677–81.

Thresholds, accuracy, and solution methods

A complete evolution record states:

Item	Required statement
boundary data	distributions and covariance at $\mu_0$ in a named scheme
kernels	spacelike order and index convention
coupling	running order, reference input, and active flavors
thresholds	matching scales, PDF matching matrices, and heavy-mass scheme
solver	$x$-space grid, Mellin inversion, or another method with tolerances
validation	sum-rule residuals, convergence, and round-trip error

When the scale crosses a heavy flavor, coupling evolution, PDFs, and coefficient functions must all be matched. A variable-flavor description is not obtained by merely adding a new component to the PDF vector.

Logarithmic accuracy also must be consistent: an evolution kernel at one order, boundary distributions extracted with another, and coefficient functions at a third do not automatically combine into the highest of those labels.

Failure modes and limits

Space- versus timelike evolution. Fragmentation functions obey timelike equations. Leading kernels are related by a transpose under a suitable convention, but the equality does not persist unchanged beyond leading order.

Small $x$. When $\alpha_s\ln(1/x)$ is order one, fixed-order DGLAP kernels may need high-energy resummation and matching.

Large $x$. Threshold logarithms and power corrections are enhanced as $x\to1$; evolution alone is not an endpoint prediction.

Numerical drift. Positivity-looking curves can still violate exact moments. Monitor valence and momentum sum rules at every scale and make grid refinement part of the validation record.

Handoff

The evolution map passed onward is

$$U^{\mathrm S}_{ij}(\mu,\mu_0; P^{\mathrm S},\alpha_s,n_f,\mathcal S,\text{thresholds}),$$

together with boundary distributions, their covariance, and numerical tolerances. Use the superscript $\mathrm S$ to prevent accidental reuse as the timelike map for fragmentation functions. Use high-energy QCD when the small-$x$ logarithm changes the counting.

References

• Altarelli, Guido, and Giorgio Parisi. “Asymptotic Freedom in Parton Language.” Nuclear Physics B 126, no. 2 (1977): 298–318. DOI.
• Schwartz, Matthew D. Quantum Field Theory and the Standard Model. Cambridge University Press, 2014, §32.2, pp. 677–81. DOI.