Analytical solution to DGLAP integro-differential equation via complex maps in domains of contour integrals

This paper details the analytical solution to the DGLAP integro-differential equation within a simplified QCD model by employing complex maps and contour integrals to relate the solution to Bessel functions, establish a connection to the BFKL equation, and express the inverse Laplace transformation as a Barnes contour integral.

The Gist

Imagine you are trying to predict how a crowd of people (particles called partons) moves and interacts inside a tiny, invisible ball (a proton) when you smash it with a high-energy beam. Physicists have a specific rulebook for this movement called the DGLAP equation.

Think of the DGLAP equation as a massive, complicated recipe for a soup. It tells you how the ingredients (quarks and gluons) change as you cook them at different temperatures (energy levels). However, this recipe is written in a very strange language: it's an "integro-differential equation." In plain English, this means the recipe doesn't just tell you how to mix ingredients; it tells you how the entire history of the soup affects the current taste, involving infinite loops and complex math that is incredibly hard to solve directly.

The Problem: A Maze of Infinite Paths

Usually, to solve this recipe, physicists use a technique called "residue calculus." Imagine trying to find a specific person in a giant, dark maze by calculating the exact coordinates of every single turn. It works, but it's tedious, prone to errors, and gets exponentially harder if you add more ingredients (higher levels of precision).

The authors of this paper, Gustavo and Igor, asked: "Is there a shortcut? Can we map this maze to a straight road?"

The Solution: The "Magic Map" (Complex Maps)

The authors propose a clever trick using Complex Maps.

The Analogy:
Imagine you are looking at a distorted, warped photograph of a city. The streets are twisted, and the buildings are stretched. Trying to navigate this photo is a nightmare.

• The Old Way: You try to walk through the twisted streets, calculating every step carefully.
• The Authors' Way: They take the photo and run it through a "magic lens" (a complex mathematical transformation). Suddenly, the warped city straightens out. The twisted streets become a perfect grid, and the distorted buildings look like standard, easy-to-read blocks.

In the paper, they take the messy DGLAP equation and apply this "magic lens" (a complex map) to the variables. This transforms the difficult, looping math into something much simpler: a Laplace Transform.

The "Jacobian" and the "Barnes Integral"

Here is where the paper gets technical, but we can keep it simple with two more analogies:

1. The Jacobian (The Stretch Factor):
   When you stretch a rubber sheet (the complex map), some parts stretch more than others. The "Jacobian" is just a measure of how much the sheet stretched at any given point. The authors realized that the messy solution to the particle soup recipe could be described simply as the "stretch factor" of this magic map.

2. The Barnes Integral (The Universal Translator):
   Once they found the stretch factor, they had a new kind of math problem. Instead of a twisted maze, they had a specific type of loop integral. They realized this new loop looked exactly like a famous, pre-solved puzzle known as a Barnes Integral.

   • The Metaphor: Imagine you have a weird, unique puzzle piece. Instead of trying to figure out how it fits from scratch, you realize it's actually a standard piece from a "Universal Puzzle Set" (the Barnes Integrals). Because this set is well-studied, you can immediately look up the answer in a giant encyclopedia (the tables of special functions) without doing the hard work yourself.

Why Does This Matter?

The authors show that by using these "magic maps," the solution to the particle soup recipe turns out to be a Bessel Function.

• What is a Bessel Function? Think of it as a specific, well-known wave pattern (like the ripples in a pond). It's a "famous face" in the world of math.
• The Benefit: Before this, physicists had to do heavy lifting to find the answer. Now, they can say, "Oh, this is just a Bessel function!" and instantly know its properties.

The Big Picture

The authors are essentially building a translation tool for physicists.

• Old Method: "Let's calculate the infinite loops manually. It will take us 100 pages of math."
• New Method: "Let's apply a complex map. It turns the infinite loops into a standard Bessel function. We can look it up in the dictionary."

This is crucial for the future of computing. If you want to build a computer program (like a neural network) to predict how particles behave in the Large Hadron Collider, you don't want the computer to struggle with infinite loops every time. You want it to recognize patterns. By converting these complex equations into standard "Barnes Integral" forms, the authors are making it much easier for computers to classify, understand, and predict the behavior of the universe's smallest building blocks.

In summary: They found a way to turn a tangled ball of yarn (the DGLAP equation) into a straight, labeled string (a Bessel function) by using a mathematical magic trick (complex maps), making it much easier for scientists and computers to understand how the universe works at its most fundamental level.

Technical Summary

1. Problem Statement

The paper addresses the challenge of solving the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) integro-differential equation, which governs the evolution of parton distribution functions (PDFs) in Quantum Chromodynamics (QCD).

• The Difficulty: While the DGLAP equation can be converted into a first-order differential equation via Mellin moments, the inverse Mellin transformation required to return to Bjorken-$x$ space typically involves complex contour integrals.
• Current Limitations: Standard methods rely on calculating residues (Cauchy integral formula), which often leads to infinite series that are difficult to classify or sum, especially in realistic QCD scenarios involving running couplings and higher-order corrections.
• The Goal: The authors aim to develop a systematic strategy to transform these contour integrals into a form that allows for classification via special functions (specifically generalized hypergeometric functions) and facilitates the construction of computational algorithms. They focus on a simplified QCD model where the splitting function contains only one term, yielding a solution expressible via Bessel functions.

2. Methodology

The authors propose a novel approach based on complex diffeomorphisms (complex maps) applied to the integration variable in the Mellin moment plane. The methodology proceeds in three main stages:

A. Complex Mapping to Jacobian Form

Instead of evaluating the contour integral directly via residues, the authors introduce a change of variables (a complex map) $N \to M$ in the Mellin plane.

• They start with the contour integral solution for the simplified DGLAP model:
  $$\phi(x, u) = \int_{-1+\delta-i\infty}^{-1+\delta+i\infty} dN \, x^{-N} \frac{1}{N+1} u^{1/(N+1)}$$
• By defining a specific map $M(N)$, they transform the integrand. The Jacobian of this transformation, $\frac{dN}{dM}$, becomes part of the new integrand.
• This transformation converts the original integral into a Laplace transform of the Jacobian. In the specific model considered, this results in an integral representation of the Bessel function $I_0$.

B. Connection to Standard Tables

The authors demonstrate that the resulting Laplace transform of the Jacobian corresponds to standard integrals found in tables (e.g., Gradshteyn and Ryzhik). This step validates that the complex map successfully simplifies the structure of the integrand into a known, solvable form.

C. Transformation to Barnes Integrals

The core innovation lies in the second transformation. The authors convert the Laplace transform of the Jacobian into a Barnes contour integral.

• Barnes integrals are defined as contour integrals where the integrand is a ratio of products of Euler Gamma functions.
• The authors utilize the inverse Laplace transformation and properties of the Gamma function to rewrite the Jacobian form (which often involves multivalued functions like square roots and requires integration over branch cuts) into a Barnes integral form.
• This transformation eliminates the need for integration over cuts in the complex plane, as the Barnes representation relies on the analytic properties of Gamma functions.

3. Key Contributions

1. Systematic Transformation Strategy: The paper provides a detailed derivation showing how to map a DGLAP contour integral solution to a Barnes integral via an intermediate Jacobian/Laplace form. This offers a unified framework for classifying solutions.
2. Duality Verification: The work reinforces the connection between the DGLAP and BFKL equations. It confirms that a complex map in the Mellin moment plane can relate the DGLAP solution to the BFKL equation, treating them as duals.
3. Handling Multivalued Functions: The authors show how to avoid the complexities of integrating over branch cuts (required when dealing with the square-root Jacobian of the complex map) by converting the expression into a Barnes integral, which is single-valued in its standard representation.
4. Algorithmic Foundation: By reducing the problem to ratios of Gamma functions, the authors lay the groundwork for constructing computer algorithms. These algorithms can systematically classify solutions in terms of generalized hypergeometric functions, which is crucial for handling the high complexity of anomalous dimensions at higher loop orders (NLO, NNLO).

4. Key Results

• Analytical Solution: For the simplified model (single-term splitting function), the authors explicitly derive that the solution is the Bessel function $I_0(2\sqrt{\ln u \ln(1/x)})$.
• Integral Identity: They prove the identity relating the Jacobian form to the Barnes integral:
  $$\int_{C''} dM \frac{e^{Mx}}{\sqrt{M^2-1}} = \left(\frac{2}{x}\right)^2 \int_{C_2} du \frac{\Gamma(1-u)}{\Gamma(u)} \left(-\frac{x^2}{4}\right)^u$$
  The left side represents the Jacobian/Laplace form, and the right side is the Barnes integral.
• Series Expansion: The Barnes integral is shown to expand directly into the standard series definition of the Bessel function $I_0$, confirming the validity of the transformation.

5. Significance

• Computational Utility: The primary significance is the potential for automated calculation. Current QCD calculations often rely on numerical packages or complex algebraic geometry tools (shuffle products). This method suggests a path where contour integrals can be systematically converted into Barnes integrals, allowing for the automatic classification of results using known special functions.
• Pedagogical and Theoretical Insight: The paper clarifies the mathematical structure underlying the DGLAP solution in the small-$x$ regime, capturing the "Bessel-like" behavior of gluon distributions.
• Neural Network Training: The authors suggest that these simple, analytically solvable models can serve as training data or consistency checks for neural networks used in global PDF fits (e.g., NNPDF), helping to constrain parameters in regimes where numerical solutions might behave poorly (low momentum transfer).
• Future Development: The strategy is designed to be scalable. While demonstrated on a simple model, the authors argue that the method of mapping to Barnes integrals is applicable to the full, multi-term DGLAP equations, offering a robust alternative to current residue-based summation techniques.