⚛️ phenomenology

Scheme transformations as the gauge group of DGLAP: sum rules, classification and solution

This paper establishes that DGLAP evolution kernels can be treated as gauge fields under scheme transformations, providing a complete classification of symmetry-preserving transformations, a method to repair sum-rule violations, and a novel solution technique that reduces the evolution equations to a single second-order differential equation while exactly respecting physical sum rules.

Original authors: Tommaso Rainaldi

Published 2026-09-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Tommaso Rainaldi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Inside every proton, a chaotic storm of smaller particles called quarks and gluons whips around, carrying the particle's mass and momentum. To understand how these particles behave when they smash into one another at high speeds, physicists rely on a set of rules known as the DGLAP equations. These rules describe how the density of quarks and gluons changes as the energy of the interaction shifts. However, there is a catch: the way we define these densities depends on the mathematical "ruler" or scheme we choose to measure them. Just as a map can be drawn in different projections, each preserving some features while distorting others, physicists can choose different schemes to describe the same physical reality. For decades, researchers have used a standard ruler, but they have also known that infinitely many other choices are possible. The critical question has always been: which of these alternative rulers still tell the truth about the proton's structure? Specifically, which ones ensure that the total momentum of the proton remains conserved and that the number of quarks stays fixed, regardless of how we choose to measure them?

A physicist at Stony Brook University has now mapped the entire landscape of these possible rulers, identifying exactly which ones preserve the fundamental laws of the proton. The researcher approached the problem by treating the mathematical rules that govern the evolution of these particles not as static numbers, but as a kind of flexible framework, similar to how a gauge transformation works in electromagnetism. By viewing the problem through this lens, the study reveals that while there are many ways to redefine the densities, the universe is surprisingly strict about what is allowed. The only transformations that keep the proton's total momentum and the count of its valence quarks intact are a very specific, limited set. These allowed changes essentially boil down to two things: a reshuffling of momentum between the quarks and the gluons, and a uniform scaling of certain particle combinations. Any other change would break the fundamental accounting of the proton, causing the total momentum to drift or the number of quarks to fluctuate with energy, which contradicts what we know about the physical world.

The study goes further by showing that this limited set of allowed changes is not just a collection of random options, but a complete and exhaustive classification. The researcher proved that if a transformation respects the symmetries of the proton and keeps the total momentum and quark numbers constant, it must belong to this specific group. Conversely, any transformation outside this group inevitably violates these conservation laws. This finding is significant because it removes the guesswork from theoretical calculations. Previously, physicists had to check each new scheme individually to see if it broke the rules. Now, they have a definitive checklist. The research also clarifies the relationship between different types of particle densities, such as those that depend on transverse momentum. It demonstrates that by applying a simple, fixed correction factor derived from the conservation laws, one can convert a density defined by a transverse-momentum integral into a standard collinear density without losing the physical meaning of the total momentum or quark count.

Beyond simply classifying the rules, the paper uses this new understanding to solve the equations that describe how particles evolve. By treating the evolution process as a gauge transformation, the researcher found a way to simplify the complex system of equations into a single, manageable second-order equation. This simplification works at any level of precision, from the most basic approximations to the most detailed calculations. This approach provides a powerful tool for estimating the uncertainty in theoretical predictions. When physicists stop calculating at a certain level of detail, they are left with "missing" higher-order effects that could change the result. The new method allows them to estimate these missing pieces by varying the scheme within the allowed group. Because this variation respects the conservation laws by design, the resulting uncertainty bands are more reliable than previous methods, which sometimes inadvertently violated the fundamental rules of momentum and quark number.

The study also sheds light on the stability of these particle densities. It shows that while the total momentum and quark numbers are rigidly protected, the positivity of the densities—meaning the fact that the probability of finding a particle is never negative—is much more fragile. The research identifies a specific subset of the allowed transformations that preserves this positivity, but it finds that this subset is smaller than the full group of allowed changes. In fact, the standard evolution process itself, when run backward in energy, can generate negative densities, a known issue that this work helps to contextualize. The findings confirm that the standard way of calculating particle evolution is robust at the most basic level, but as calculations become more precise, the choice of scheme matters more. The work provides a rigorous framework for making those choices, ensuring that the theoretical descriptions of the proton remain consistent with the fundamental conservation laws that govern our universe.

In the end, this work unifies two seemingly different aspects of particle physics: the classification of mathematical schemes and the solution of the equations that describe particle evolution. By recognizing that the freedom to choose a scheme is a form of gauge symmetry, the researcher has shown that the rules governing the proton are far more constrained than previously thought. The classification is complete, meaning no other valid schemes exist outside the ones identified. This clarity allows for more precise predictions in high-energy physics experiments, where understanding the exact behavior of quarks and gluons is essential for interpreting the results of particle colliders. The study does not claim to have discovered new particles or forces, but rather to have provided a clearer, more complete map of the mathematical territory we already inhabit, ensuring that our calculations of the proton's inner workings are as solid as the conservation laws they are built upon.

Technical Summary: Scheme Transformations as the Gauge Group of DGLAP

Problem Statement
Parton distribution functions (PDFs) are not scheme-independent; their definition relies on a renormalization prescription. While the standard $\overline{\text{MS}}$ scheme is widely used, infinitely many other schemes are theoretically allowed. The central problem addressed is characterizing the entire set of scheme transformations that preserve the fundamental physical constraints of the parton model: the DGLAP evolution form, charge-conjugation and flavor symmetries, and the momentum and valence number sum rules. Furthermore, the paper seeks to utilize this gauge freedom to solve the DGLAP evolution equations and to organize estimates of missing higher-order uncertainties in a way that strictly respects these sum rules.

Methodology
The analysis is conducted in Mellin space, where convolutions become matrix products. The author treats the DGLAP evolution kernels as gauge connections and scheme transformations as gauge transformations acting on the flavor space of parton densities.

Symmetry Constraints: Using Schur's lemma, the paper decomposes the flavor space into irreducible representations of the symmetry group $G \equiv \langle C \rangle \times S_{N_f}$ (charge conjugation and flavor permutations). This reduces the general transformation matrix to a block-diagonal form consisting of a $2 \times 2$ singlet-gluon block and scalar blocks for non-singlet sectors.
Sum Rule Constraints: The momentum sum rule (at Mellin moment $z=2$ ) and valence number sum rules (at $z=1$ ) are imposed as conditions on the transformation matrix. These act as covector conditions, requiring specific eigenvectors of the transformation to remain invariant.
Boundary Conditions: Asymptotic freedom is used as a boundary condition, requiring transformations to approach the identity matrix as the coupling $\alpha_s \to 0$ (or $\mu \to \infty$ ).
Gauge Solving: In the second part of the paper, the boundary condition is relaxed to "gauge away" the evolution connection. This reduces the DGLAP system to a single second-order linear differential equation (a Riccati equation in disguise) plus a quadrature.

Key Contributions and Results

Complete Classification of Admissible Schemes:
The paper provides a complete classification of scheme transformations that preserve symmetries and sum rules to all orders. The admissible transformations form a group isomorphic to:
$\text{Aff}(1, \mathbb{R})_+ \times \mathbb{R}_{>0} \times \{1\}^2$
(for $N_f > 1$ ).
- The $\text{Aff}(1, \mathbb{R})_+$ factor represents the freedom to reshuffle momentum between the quark singlet and the gluon at $z=2$ . This confirms that the split of hadron momentum between quarks and gluons is a convention.
- The $\mathbb{R}_{>0}$ factor represents a positive rescaling of C-even non-singlet combinations.
- The $\{1\}^2$ factors indicate that the C-odd sectors (valence numbers) are completely frozen at $z=1$ ; valence numbers are scheme invariants.
  The classification is proven to be exhaustive (Theorem 2) because the map from a transformation to its effect on the kernels is injective under the asymptotic freedom boundary condition.
Positivity and Semigroup Structure:
The paper investigates which transformations preserve the positivity of PDFs in $x$ -space. It finds that the set of positivity-preserving transformations forms a semigroup, not a group, meaning they cannot be inverted while maintaining positivity.
- For moment-independent transformations, positivity restricts the freedom to a single parameter: the gluon can only transfer momentum to quarks, never receive it, and non-singlet rescalings are frozen.
- The evolution operator itself belongs to this positivity semigroup only at leading order (LO) and only for upward evolution. Downward evolution or higher-order kernels can generate negative densities.
Collinear Densities from TMD Integrals:
The paper addresses the definition of collinear densities as integrals of Transverse-Momentum-Dependent (TMD) densities. It demonstrates that the unnormalized integral (unnormalized cutoff scheme) generally violates sum rules. The paper defines a normalized cutoff scheme by applying a specific, moment-independent renormalization matrix (determined by three numbers: two for momentum, one for valence) to the unnormalized integral. This restores the sum rules exactly while maintaining the integral relation between TMDs and collinear densities.
Solution of DGLAP via Gauge Transformation:
By treating the evolution as a gauge connection, the author shows that DGLAP evolution can be solved by transforming to a "no-evolution" scheme. This reduces the problem to solving a single second-order linear ODE (Eq. 52) for the singlet-gluon block, with the solution expressible as a Frobenius series.
- This approach explicitly identifies the radius of convergence and the "resonant moments" where the series solution becomes degenerate (where indices differ by integers).
- It separates the known coefficients from the unknown ones at any perturbative order.
Uncertainty Estimation for Missing Higher Orders:
The paper proposes a method to estimate missing higher-order (MHO) uncertainties by varying the scheme within the classified admissible transformations.
- Unlike standard scale variation, which is a specific subset of scheme changes, this method varies the kernels directly using the admissible group structure.
- This ensures that the uncertainty band respects the sum rules exactly.
- The paper derives a closed-form expression (Eq. 61) for the shift in the evolution coefficients due to missing orders, linking the scheme variation directly to the unknown anomalous dimensions.

Significance and Claims
The paper claims that viewing scheme transformations as the gauge group of DGLAP provides a unifying language for several known results (such as the Adler sum rule constraints and small- $x$ resummation techniques) and yields new insights:

Completeness: The classification of sum-rule-preserving schemes is proven to be complete; no other transformations exist that satisfy the physical constraints and asymptotic freedom.
Sum Rule Repair: Sum rules are never an obstruction to a transformation; any transformation with the correct symmetries that fails the sum rules can be "repaired" by a simple normalization factor determined by the sum rules themselves.
Exact Sum Rule Preservation in Uncertainty: The proposed scheme variation method for estimating MHO uncertainties is unique in that it respects the sum rules exactly, whereas standard scale variation (while preserving them) does not explore the full space of scheme dependence, and other variations might violate them.
Structural Solution: The reduction of DGLAP to a single second-order equation provides a structural understanding of the solution, clarifying the origin of resonant moments and the convergence properties of the perturbative series.

The author notes that while the mathematical structure is rigorous, the paper does not claim immediate phenomenological superiority over existing tools (like Pegasus) but rather offers a structural framework to understand and organize the evolution and its uncertainties. The work suggests that future global fits could vary schemes within this classified group to test if sum-rule-preserving uncertainty bands differ from standard scale-variation bands.

More like this