Extended DGLAP Evolution Equation

• Extended DGLAP-Type Evolution Equations are generalizations that include mass effects, QCD×QED mixing, and non-linear multi-parton correlations.
• They unify collinear (DGLAP) and small-x (BFKL) resummations by employing novel angular-ordering techniques and subtraction schemes.
• The framework directly evolves physical observables such as structure functions, generalized parton distributions, and diffractive distributions for improved precision.

An extended DGLAP-type evolution equation refers to any generalization or augmentation of the standard Dokshitzer–Gribov–Lipatov–Altarelli–Parisi (DGLAP) equations, designed to account for additional physical regimes, degrees of freedom, or theoretical constraints beyond the conventional leading-twist, massless, collinear parton evolution formalism. Modern developments have produced several classes of such extended equations, addressing key challenges in QCD phenomenology: unified evolutions bridging collinear (DGLAP) and low-x (BFKL) resummations, mass effects, multi-parton correlations, QCD×QED mixing, generalized parton distributions, diffraction, and even extensions to accommodate quantum anomalies.

1. Motivation for Extending DGLAP Evolution

The original DGLAP equations describe the scale dependence of parton distribution functions (PDFs) under collinear factorization, resumming leading logarithms in $Q^2$. However, a wide range of theoretical and phenomenological issues call for more general frameworks:

• At small $x$, power corrections and large $\ln(1/x)$ terms become significant, motivating hybrid BFKL+DGLAP schemes.
• Physical observables (structure functions, energy flows, diffractive distributions) sometimes demand evolution equations directly for those quantities rather than for PDFs.
• Massive quarks induce kinematic thresholds and flavor-number discontinuities not present in massless evolution.
• Mixed QCD-QED corrections and electroweak anomalies modify both the splitting kernels and the structure of sum rules.
• Multi-particle and correlation-sensitive observables require evolution equations including non-linear mixing of multi-fragmentation functions.
• Generalized parton distributions (GPDs), relevant for exclusive and semi-exclusive processes, necessitate evolution kernels depending on additional kinematic variables.

These realities motivate systematic generalizations—“extended DGLAP-type equations”—incorporating these effects within a consistent perturbative QCD (and QED) formalism.

2. Unified BFKL+DGLAP Approaches and Evolution Variable Choice

At small $x$, the BFKL equation resums leading $\ln(1/x)$ terms, resulting in evolution more naturally ordered in rapidity or emission angle, while DGLAP resums $\ln(Q^2)$. The need for simultaneous resummation led to unified evolution equations for the integrated or unintegrated gluon density.

A key result is the so-called “hybrid” or “unified” equation. For the unintegrated gluon density $\mathcal{F}(x,k^2)$, it reads (Toton, 2014): $$\mathcal{F}(x,k^2) = \mathcal{F}_0(x,k^2) + \bar{\alpha}_s \!\int_x^1\! \frac{dz}{z}\int_{k_0^2}^{\infty} \frac{dl^2}{l^2} \left[ \frac{\ldots}{l^2-k^2} + \frac{k^2\mathcal{F}(x/z,k^2)}{\sqrt{4l^4 + k^4}} \right] + \frac{\bar{\alpha}_s}{k^2} \int_x^1\!dz \left( \frac{zP_{gg}(z)}{2N_c} - 1 \right) \int_{k_0^2}^{k^2} dl'^2\, \mathcal{F}(x/z, l'^2)$$ with a subtraction scheme to avoid double counting collinear logs.

Changing the evolution variable to gluon emission angle $\theta = k/(xp)$ (with $k$ the gluon $k_t$, $x$ the longitudinal momentum fraction, $p$ the parent’s momentum) provides numerical and conceptual advantages (Oliveira et al., 2014). In angle-ordered evolution, the collinear and small-$x$ contributions are treated in a single variable, and angular ordering directly enforces coherence: $$\partial_{\ln \theta^2} [xg(x, \theta)] = f_0(x, \theta_0) + \frac{\alpha_s}{2\pi} \left[ \text{BFKL kernel} + \text{DGLAP term} + \text{subtractions} \right]$$ The $\theta$-evolution framework accelerates numerical convergence and naturally handles both $\ln Q^2$ and $\ln(1/x)$ resummation.

3. Structure-Function and Physical-Basis Evolution

Conventional DGLAP equations evolve PDFs, but physical cross sections in DIS or hadron-hadron collisions depend on observable structure functions. An extended DGLAP-type approach directly evolves a complete set of structure functions in a “physical basis” (Lappi et al., 2023): $$\frac{d}{d\ln Q^2} F_i(x, Q^2) = \frac{\alpha_s}{2\pi} \sum_{j=1}^6 \left[ P_{ij} \otimes F_j \right] (x, Q^2)$$ The evolution kernel $P_{ij}(z)$ forms a $6 \times 6$ matrix constructed from the usual $P_{qq}$, $P_{qg}$, $P_{gq}$, $P_{gg}$ functions, but rotated into the structure-function basis by invertible linear maps. This formalism guarantees manifest independence from factorization schemes and unobservable intermediate PDFs, while allowing evolution to NLO or NNLO order with kernels built from the appropriate higher-order splitting and coefficient functions.

Numerical comparisons confirm the physical-basis and standard PDF-evolution results agree to within $10^{-4}$ over $10$ orders of magnitude in $Q^2$ except at extreme $x$ (Lappi et al., 2023).

4. Mass Effects and Generalization to Heavy Partons

Treatment of massive quarks in DGLAP evolution introduces flavor thresholds and kinematical phase-space boundaries. The “continuous-flavor-number scheme” (cfns) introduces activity factors $\ell_h(Q^2) = 1/(1+4m_h^2/Q^2)$ for each heavy quark flavor $h$ (Pascaud, 2011). Splitting kernels and $\alpha_s$ running are modified: $$\widetilde{P}_{oi}(x,Q^2) = P_{oi}\Bigl( \frac{x}{\ell_o(Q^2)} \Bigr) \frac{1}{\ell_o(Q^2)}\, \theta\bigl(x < \ell_o(Q^2)\bigr)$$

$N_f(Q^2) = \sum_{h=1}^6 \ell_h(Q^2)$

This produces fully smooth parton distributions and observables through heavy-quark thresholds, preserving momentum and flavor sum rules, and avoiding discontinuities or ad-hoc matching scales.

5. Correlations and Collinear Evolution Beyond Single-Particle DGLAP

Jet substructure and multi-hadron observables depend on correlations beyond single-hadron fragmentation, requiring an extended RG evolution for fragmentation (track) functions $T_i(x,\mu)$ that describe the probability a parton $i$ deposits energy fraction $x$ into a subset of hadrons (Chen et al., 2022): $$\frac{d}{d\ln \mu^2} T_i(x, \mu) = a_s(\mu) [ K^{(0)}_{i\rightarrow i} T_i + K^{(0)}_{i\rightarrow i_1 i_2}\otimes T_{i_1} T_{i_2} ](x) + a_s^2(\mu)[ K^{(1)}_{i\rightarrow i} T_i + K^{(1)}_{i\rightarrow i_1 i_2}\otimes T_{i_1} T_{i_2} + K^{(1)}_{i\rightarrow i_1 i_2 i_3}\otimes T_{i_1}T_{i_2}T_{i_3} ](x) + \mathcal{O}(a_s^3)$$ Here, the non-linear structure encodes multi-particle correlations, with explicit $1\rightarrow 3$ kernels at NLO, going beyond the strictly linear evolution of DGLAP. This formalism encompasses standard single-hadron evolution as a limit and, upon repeated delta-function substitution, reproduces the tower of multi-hadron fragmentation function evolutions. It enables NNLO-accurate predictions for observables such as the energy fraction in charged particles in $e^+e^- \rightarrow$ hadrons (Chen et al., 2022).

6. QCD×QED Mixing and Anomalous Extensions

Inclusion of QED corrections alters splitting kernels and induces charge-separated evolution equations. The generalized system, augmented to $\mathcal{O}(\alpha_s\alpha)$, features

$$P_{ij}(x) = a_s P_{ij}^{(1,0)}(x) + a P_{ij}^{(0,1)}(x) + a_s a P_{ij}^{(1,1)}(x) + \ldots$$

with coupled evolution for quark, gluon, and photon distributions (Sborlini et al., 2016, Florian et al., 6 May 2025). Mixed QCD$\otimes$QED kernels with explicit charge dependence quantify the corrections to photon PDFs and break degeneracy between $u$ and $d$ channels, which is critical for percent-level LHC precision analyses. Analytic solutions via Mellin space further enhance computational efficiency and provide exact control over mixed corrections (Florian et al., 6 May 2025).

In scenarios with a $U(1)_B$ (baryon-number) anomaly, the DGLAP equations receive inhomogeneous, scale-localized contributions, modifying the standard conservation sum rules (Sainapha, 2018): $$\frac{\partial f_q(x,Q^2)}{\partial\ln Q^2} = \frac{\alpha_s}{2\pi}\ldots + \mathcal{A}_q(Q^2) \delta(1-z)\frac{d}{d\ln Q^2} Q_\text{CS}$$ where $Q_\text{CS}$ is the Chern–Simons number. This ensures self-consistency with anomaly-induced baryon-number violation and directly encodes Sakharov's requirements for baryogenesis within the evolution structure.

7. Applications to GPDs, Diffraction, and Low-x Physics

Extended DGLAP equations have been formulated for generalized parton distributions (GPDs), yielding kernels $K(x,z,\xi)$ with explicit dependence on skewness $\xi$ (Bertone et al., 2022). These kernels reproduce both standard DGLAP and ERBL equations in their respective regions, and continuity at $x=\pm\xi$ is maintained by construction. Open-source evolution codes (APFEL++, PARTONS) now implement these kernels, including heavy-quark thresholds.

For diffractive DIS, DGLAP-type equations are derived for diffractive gluon and sea-quark distributions $\beta G^D(\beta, Y_0, \xi)$, with physical input fixed by Color Glass Condensate (CGC) dynamics (Contreras et al., 2018). The evolution structure is identical to DGLAP but includes a nontrivial, perturbative initial condition and relates directly to the underlying dilute gluon emission regime in diffraction.

At very low $x$ and high $Q^2$, JIMWLK/BK evolution equations receive a DGLAP-resummed addition in $\ln Q^2$, producing a combined evolution for the dressed gluon $S$-matrix (Armesto et al., 31 Jan 2025): $$\Bigl(\frac{\partial}{\partial Y} + \frac{\partial}{\partial\ln Q^2} \Bigr) \mathbb{S}_Q = - (H_\text{JIMWLK} + H_\text{DGLAP})\,\mathbb{S}_Q$$ This simultaneously resums high-energy (small-$x$) and collinear (large $Q^2$) logarithms, leading to universal, coupling-independent patterns of unitarity violation.

---

The spectrum of “extended DGLAP-type evolution equations” thus subsumes methods addressing multi-scale resummation, mass thresholds, higher-order and mixed (QCD/QED) corrections, multi-hadron and GPD evolution, non-linear and correlation-sensitive observables, and anomalies, all within the RG paradigm. These advances underpin the current standard of precision QCD phenomenology and are critical for future collider and cosmic-ray analyses.