Algorithm to find an all-order in the running coupling solution to an equation of the DGLAP type

This paper proposes an algorithm utilizing complex analysis to solve DGLAP-type integro-differential equations to all orders in the running coupling $\alpha$, offering a simplified method for calculating the necessary contour integrals compared to existing techniques.

The Gist

The Big Picture: Predicting the Unpredictable

Imagine you are trying to predict the weather inside a hurricane. You know the laws of physics (the "equations"), but the storm is so chaotic and changes so fast that solving the math directly is like trying to count every single raindrop while the wind is blowing them away.

In the world of particle physics, scientists study protons (the building blocks of atoms) by smashing them together at high speeds. Inside a proton, there are smaller particles called partons (quarks and gluons) zooming around. The paper by Igor Kondrashuk is about a new, clever way to calculate how these particles behave and change as the energy of the collision changes.

The Problem: A Mathematical Maze

The scientists use a famous set of rules called the DGLAP equations to describe how these particles evolve. Think of these equations as a complex recipe for a cake.

• The Ingredients: The "splitting functions" (how one particle breaks into two).
• The Cooking Time: The "running coupling" (how the strength of the force changes as you cook).

The problem is that the recipe is written in a language called Integro-Differential Equations. This is a mix of calculus that is incredibly hard to solve directly. It's like trying to bake a cake where the instructions tell you to "mix the batter while simultaneously measuring the heat of the oven and the humidity of the room, all while the oven temperature keeps changing."

Traditionally, to solve this, scientists use a mathematical trick called the Mellin Transform. Imagine taking a photo of the cake, turning it into a black-and-white negative (the Mellin space), solving the math there, and then trying to print the photo back into color (the real world).

• The Old Way: To get the photo back, they had to add up an infinite number of tiny layers (a series expansion). It was slow, messy, and required massive computer power.

The Solution: A Secret Shortcut (The "Magic Door")

The author proposes a new algorithm. Instead of struggling to add up all those layers, he suggests opening a "secret door" in the math.

The Analogy: The Shape-Shifting Map
Imagine you are lost in a dense, foggy forest (the complex math problem).

1. The Old Way: You try to walk through the forest, counting every tree and stepping over every root. It takes forever.
2. The New Way: You realize that if you look at the forest from a specific angle (a complex diffeomorphism), the trees suddenly line up perfectly into a straight, paved highway.

In this paper, the author uses Complex Analysis (math involving imaginary numbers) to "reshape" the problem. He takes the messy, twisted math and stretches and twists it (like a rubber sheet) until it looks like a standard, simple problem that everyone already knows how to solve.

The Key Trick: The "Jacobian" and the "Laplace"

Here is the core of his method, broken down:

1. The Transformation: He changes the variables in the equation. He swaps the "Mellin moment" (the weird negative photo) for a new variable.
2. The Jacobian: When you stretch a rubber sheet, the area changes. The "Jacobian" is just a number that tells you how much the area stretched. The author realized that if you calculate this stretch factor correctly, the messy equation turns into a Laplace Transform.
3. The Table Trick: Laplace Transforms are like the "Cheat Codes" of math. There are huge tables (like a dictionary) where you can look up the answer for almost any standard Laplace problem.
   • Before: "I have to derive this from scratch using infinite series."
   • Now: "I have transformed my problem into a shape that matches entry #42 in the table. The answer is right there!"

Why This Matters: The "Running" Engine

The paper specifically tackles the case where the "engine" of the universe (the strong nuclear force) changes its strength as things move faster. This is called the running coupling.

• Fixed Coupling: Imagine driving a car at a constant speed on a flat road. Easy to predict.
• Running Coupling: Imagine driving a car where the engine power changes every second depending on the terrain. Much harder to predict.

The author's algorithm works even when the engine is changing. By using these complex "shape-shifters," he can find the solution for all levels of complexity (all orders) without getting bogged down in infinite calculations.

The "Duality" Connection

The paper also mentions a connection between two different theories: DGLAP (how particles split) and BFKL (how particles interact at high speeds).

• Think of DGLAP and BFKL as two different languages describing the same story.
• The author found a "Rosetta Stone" (the complex mapping) that translates between them instantly. This proves that the messy math of one theory is actually just a distorted version of the clean math of the other.

The Conclusion: Less Math, More Physics

The main takeaway is this:
Instead of forcing the computer to do millions of calculations to solve a messy equation, the author found a way to rearrange the furniture in the room so the problem solves itself.

By turning the problem into a Laplace Transform (a standard table lookup) and then converting that into a Barnes Integral (a known mathematical shape), he bypasses the need for heavy numerical simulations.

In short: He didn't build a faster car to get to the destination; he found a shortcut tunnel that makes the trip instant. This allows physicists to predict how protons behave with much greater precision and less computing power.

Technical Summary

1. Problem Statement

The paper addresses the challenge of solving the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) integro-differential equations (IDEs) for parton distribution functions (PDFs) in Quantum Chromodynamics (QCD).

• Context: While the evolution of PDFs with respect to the momentum transfer scale ($Q^2$) is well understood, obtaining an analytic solution for all orders in the running coupling $\alpha(u)$ remains difficult.
• Current Limitations:
  • Standard methods involve taking Mellin moments to convert IDEs into differential equations, solving them, and then performing an inverse Mellin transformation back to Bjorken-$x$ space.
  • The inverse transformation typically requires expanding the exponential function in terms of the coupling and calculating residues on the complex plane. This results in an infinite series where higher-order calculations become increasingly complex.
  • Existing numerical packages (e.g., PartonEvolution, QCDNUM) or advanced analytical tools based on algebraic geometry (shuffle products, harmonic polylogarithms) are often computationally heavy or limited to specific orders.
  • A common philosophy in the field is to parameterize PDFs in Mellin space and perform the final inversion numerically. The author argues this is unnecessary if a more efficient analytic path exists.

2. Methodology

The author proposes a novel algorithm based on complex analysis and complex diffeomorphisms (variable transformations in the complex plane) to simplify the inverse Mellin transformation.

• Core Strategy: Instead of directly calculating residues for the inverse Mellin transform, the method transforms the integral into an inverse Laplace transformation of the Jacobian corresponding to a specific complex mapping.
• Key Steps:
  1. Complex Diffeomorphism: A change of variables is performed in the complex plane of the Mellin moment ($N$). A new variable $M$ is introduced such that the structure of the integrand becomes uniform.
  2. Duality Relation: The method leverages the duality between the DGLAP and BFKL (Balitsky-Fadin-Kuraev-Lipatov) equations. By defining a mapping where the anomalous dimension $\gamma(N)$ relates to the new variable $M$ via $\gamma(N) = M$, the author constructs an inverse function $\chi(M) = N$.
  3. Transformation to Laplace Form: The original contour integral (vertical line in the $N$-plane) is mapped to a new contour. Through this mapping, the integral is recast as the inverse Laplace transformation of the Jacobian of the diffeomorphism.
  4. Contour Deformation to Barnes Integrals:
     • The vertical line of the inverse Laplace transform is deformed into a Hankel contour.
     • This deformation allows the extraction of the Euler Gamma function in the denominator.
     • The resulting integral is identified as a Barnes integral, which is the standard representation for generalized hypergeometric functions.
• Handling the Running Coupling: The algorithm extends to the case of a running coupling $\alpha(u)$ by treating the coupling as the evolution variable. A new variable $M$ is chosen such that $M\alpha = F(\alpha, N)$, where $F$ is the integral of the anomalous dimension over the coupling, ensuring the integrand remains simple.

3. Key Contributions

• New Analytic Algorithm: Proposes a systematic way to find all-order solutions in the running coupling by reducing the inverse Mellin transformation to the calculation of inverse Laplace transformations of Jacobians.
• Simplification of Integrals: Demonstrates that complex mappings can reduce complex integrands to standard table integrals (specifically Barnes integrals), bypassing the need for infinite series expansions of residues.
• Unification of DGLAP and BFKL: Provides a mathematical framework showing how complex diffeomorphisms relate the DGLAP and BFKL equations, offering a deeper understanding of their duality in the complex plane.
• Avoidance of Numerical Inversion: Challenges the standard practice of performing the final inversion numerically, suggesting that the entire process can be handled analytically (or semi-analytically via standard functions) up to the final convolution with initial conditions.

4. Results

• Toy Model Verification: The author validates the method using a simplified toy model where the splitting function is $P_{GG}(z) = 2z$.
  • The standard DGLAP solution yields a Bessel function: $x I_0(2\sqrt{\frac{\alpha}{\pi} \ln u \ln \frac{1}{x}})$.
  • Using the proposed complex diffeomorphism (mapping $N \to M$), the author successfully reproduces the exact same Bessel function result, confirming the validity of the transformation.
• Generalization to Real QCD:
  • The paper acknowledges that in real QCD, the inverse function $\chi(M)$ involves the Euler $\psi$ function (digamma function), which does not have a simple explicit representation.
  • However, the author notes that computer algorithms exist to evaluate this inverse function, and the structure of the resulting integrals (Barnes integrals) remains valid.
  • The method suggests that the "new special functions" often required for analytic inversions in QCD originate from these specific complex diffeomorphisms.

5. Significance

• Computational Efficiency: By converting complex contour integrals into Barnes integrals (standard hypergeometric functions), the method offers a potentially more efficient route for calculating PDF evolution at higher orders (NNLO and beyond) compared to residue calculus or numerical integration.
• Theoretical Insight: The work provides a geometric interpretation of the DGLAP/BFKL duality, suggesting that the complexity of QCD evolution equations can be managed through appropriate coordinate transformations in the complex Mellin plane.
• Practical Application: The algorithm reduces the problem of finding all-order solutions to the evaluation of Jacobians for specific complex maps, which can be tabulated or computed systematically, potentially improving the speed and accuracy of parton distribution function evolution codes used in high-energy physics experiments.

In summary, Kondrashuk presents a sophisticated mathematical technique that uses complex analysis to bypass the traditional bottlenecks in solving DGLAP equations with a running coupling, offering a pathway to fully analytic solutions via Barnes integrals.