Multiplicity distribution of produced gluons in deep inelastic scattering: main equations and their homotopy solutions for heavy nuclei

This paper derives new equations for the multiplicity distribution of produced gluons in deep inelastic scattering on heavy nuclei using a dipole approach and AGK cutting rules, develops a homotopy method to solve them, and provides an analytical solution for large multiplicities to calculate the resulting gluon entropy.

The Gist

The Big Picture: Shattering a Glass Ball

Imagine you have a very fast-moving glass ball (a proton) and you smash it with a tiny, high-energy hammer (a virtual photon). This is what happens in Deep Inelastic Scattering (DIS).

When they collide, the glass ball doesn't just break into two pieces; it shatters into a chaotic cloud of tiny shards (gluons and quarks). Physicists want to know: How many shards do we get? Is it usually 10? 100? 1,000? And how does the number of shards change depending on how hard we hit it?

This paper tries to answer that question for heavy nuclei (like a big chunk of lead) using the rules of Quantum Chromodynamics (QCD), the physics of how particles stick together.

---

The Three Main Achievements

The authors claim they did three big things to solve this puzzle.

1. Rewriting the Rulebook (The "Dipole" Approach)

The Problem: Traditionally, to calculate how many shards fly out, physicists used a set of complex rules called the AGK cutting rules. These rules are like a specific instruction manual for counting broken pieces in a chaotic mess.

The Solution: The authors said, "We don't need to rely on that old manual directly." Instead, they looked at the collision as a game of billiard balls splitting apart.

• The Analogy: Imagine a fast-moving cue ball (the proton) hitting a target. As it moves, it splits into smaller and smaller balls (dipoles).
• They derived a new set of equations based on how these "billiard balls" split and interact. Surprisingly, their new math leads to the exact same results as the old AGK rules, but it feels more natural because it's built from the ground up using the "splitting balls" picture.

2. The "Homotopy" Ladder (Solving the Math)

The Problem: The equations they wrote down are non-linear. In plain English, this means they are incredibly messy. If you try to solve them like a standard algebra problem, you get stuck. It's like trying to untangle a knot that keeps tightening itself the more you pull.

The Solution: They used a technique called the Homotopy Approach.

• The Analogy: Imagine you are trying to climb a steep, foggy mountain (the solution). You can't see the top.
  • Step 1: You build a small, easy ramp (the first iteration) that gets you partway up. You can calculate this easily.
  • Step 2: You use a computer to build a slightly steeper ramp based on the first one.
  • Step 3: You keep adding ramps, getting closer and closer to the top.
• They found that after just three or four steps, their "ramp" was accurate enough (within 0.2% error) to give a reliable answer for small numbers of shards.

3. The "Big Number" Secret (Entropy and Chaos)

The Problem: What happens when the collision is so violent that it produces a massive number of shards (large $n$)? The "ramp" method above gets too slow and messy for huge numbers.

The Solution: They found a special analytical formula (a clean, direct equation) that works specifically when the number of shards is huge.

• The Analogy: Think of a crowd of people. If you have 5 people, you can track everyone individually. If you have 10,000 people, you stop counting individuals and start looking at the "density" of the crowd.
• The Discovery: They found that when the number of shards is huge, the distribution follows a specific pattern called KNO scaling. It's like a universal law of shattering: no matter how big the target is, the shape of the distribution of shards looks the same if you scale it correctly.

The "Entropy" Surprise:
Using this formula, they calculated the Entropy (a measure of disorder or information) of the produced gluons.

• They found a beautiful, simple result: Entropy = $\ln(\text{Number of Shards})$.
• Why it matters: This confirms a previous theory that the "messiness" of the collision is directly tied to the logarithm of how many particles are created. It's a fundamental link between the chaos of the collision and the information it contains.

---

The Catch (The "Matching" Problem)

The paper admits one thing they haven't fully solved yet: The Middle Ground.

• They have a great method for small numbers of shards (using the "ramp" method).
• They have a great method for huge numbers of shards (using the "Big Number" formula).
• The Gap: They don't have a perfect way to smoothly connect the two in the middle. It's like having a map of the bottom of the mountain and the top of the mountain, but the middle section is still foggy.

However, they argue that for the most important questions (like the total entropy), the "huge number" part dominates, so their answer is still correct for the big picture.

Summary for the Everyday Person

Imagine you are trying to predict how a specific type of glass shatters when hit by a laser.

1. New Rules: You invented a new way to describe the glass breaking that matches the old rules but makes more sense physically.
2. Step-by-Step Math: You built a computer program that solves the breaking pattern step-by-step, getting more accurate with each step.
3. The Big Picture: You discovered that when the glass shatters into thousands of pieces, there is a simple, universal law governing the chaos.
4. The Result: You proved that the "disorder" of the shattered glass is exactly equal to the natural log of the number of pieces.

This paper helps physicists understand the fundamental "messiness" of the universe at its smallest scales, confirming that even in chaos, there is a mathematical order.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenge of calculating the multiplicity distribution of produced gluons in Deep Inelastic Scattering (DIS) processes at high energies within the framework of Quantum Chromodynamics (QCD). Specifically, the authors focus on:

• Determining the cross-sections ($\sigma_n$) for the production of $n$-cut Pomerons in the final state.
• Solving the non-linear evolution equations governing these cross-sections for heavy nuclei targets.
• Calculating the entropy of the produced gluons ($S_E$) and verifying its relationship to the multiplicity distribution, a topic of recent interest regarding entanglement entropy in DIS.
• Bridging the gap between the standard AGK (Abramovsky-Gribov-Kancheli) cutting rules and the dipole approach to QCD without explicitly invoking the cutting rules in the derivation.

2. Methodology

The authors employ a multi-faceted approach combining dipole formalism, evolution equations, and a specific mathematical technique for solving non-linear differential equations.

A. Theoretical Framework: Dipole Approach & AGK Equations

• Dipole Formalism: The scattering amplitude is treated using the color dipole model, where a fast hadron is represented as a wave function of dipoles.
• Derivation without Explicit AGK Rules: Instead of applying AGK cutting rules as an external input, the authors derive the evolution equations for $\sigma_n$ (cross-sections of $n$-cut Pomerons) directly from the dipole wave function and the $S$-matrix unitarity. They demonstrate that this derivation naturally reproduces the coefficients found in the standard AGK rules.
• Equation Structure: The evolution equation for $\sigma_n$ involves a linear BFKL (Balitsky-Fadin-Kuraev-Lipatov) kernel and non-linear terms representing dipole interactions and shadowing effects. The equation is derived by summing contributions from $k$-Pomeron exchanges and projecting onto fixed multiplicities.

B. The Homotopy Approach

To solve the resulting non-linear partial differential equations, the authors utilize the Homotopy Perturbation Method (HPM):

• Decomposition: The non-linear operator equation $L[u] + NL[u] = 0$ is split into a linear part $L$ (solvable analytically) and a non-linear part $NL$.
• Iterative Solution: A homotopy parameter $p$ is introduced ($H(p, u) = L[u_p] + p NL[u_p] = 0$). The solution is expanded as a series in $p$: $u = u_0 + p u_1 + p^2 u_2 + \dots$.
• Convergence: The authors calculate the first few iterations ($p=0, 1, 2, \dots$) numerically and analytically. They find that the series converges rapidly, with the fourth iteration providing accuracy better than 2-3% for small $n$.

C. Kinematic Regions

The analysis distinguishes between two kinematic regions based on the variable $z = \ln(r^2 Q_s^2)$:

• Region I: Geometric scaling region where the amplitude depends only on $z$.
• Region II: The region where the initial conditions for heavy nuclei (McLerran-Venugopalan formula) apply, requiring dependence on both rapidity and dipole size.

3. Key Contributions

1. New Derivation of $\sigma_n$ Equations: The paper provides a rigorous derivation of the evolution equations for $n$-cut Pomeron cross-sections using the dipole approach, showing they are equivalent to those derived via AGK rules but grounded in the wave function evolution.
2. Homotopy Solution for Small $n$: The authors successfully apply the homotopy approach to solve the evolution equations for small multiplicities ($n=1, 2, 3, \dots$). They demonstrate that the method is stable and convergent, allowing for the calculation of moments of the multiplicity distribution.
3. Analytical Solution for Large $n$: For large multiplicities ($n \gtrsim N(z)$), where $N(z)$ is the typical multiplicity, the authors derive an analytical solution. They show that in this limit, the distribution follows KNO (Koba-Nielsen-Olesen) scaling.
4. Entropy Calculation: Using the large-$n$ solution, they calculate the von Neumann entropy of the produced gluons.

4. Key Results

• Multiplicity Distribution ($\sigma_n$):

  • For small $n$, the distribution is calculated numerically using the homotopy method. The corrections from higher iterations decay rapidly.
  • For large $n$, the solution exhibits the form:
    $$\sigma_n(z) \propto \frac{1}{N(z)} \Psi\left(\frac{n}{N(z)}\right)$$
    where $\Psi(\xi) = \xi e^{-\xi}$ is the KNO scaling function, and $N(z)$ is the characteristic multiplicity scaling as $N(z) \sim z \exp(z^2/2\kappa)$.
  • The distribution decreases as $\exp(-z^2/2\kappa)$ at large $z$.

• Entropy of Produced Gluons ($S_E$):

  • The authors calculate the entropy defined as $S_E = -\sum P_n \ln P_n$.
  • Result: $S_E = \ln(N(z))$.
  • This result confirms the hypothesis (from Ref. [8]) that the entropy in DIS is closely related to the entanglement entropy between the probed spatial region and the rest of the proton.

• Inelastic Cross Section:

  • The sum of all $\sigma_n$ yields an inelastic cross-section $\sigma_{in}$.
  • The analytical large-$n$ solution suggests $\sigma_{in} \to 2$ at large $z$, which violates unitarity (which requires $\sigma_{in} \to 1$).
  • Resolution: The authors argue that the large-$n$ solution is valid only for $n \ge N(z)$. The contribution from small $n$ (where the homotopy solution applies) is negligible for the total cross-section at high energies, but the matching between the small-$n$ and large-$n$ regimes remains an open technical challenge.

5. Significance and Conclusion

• Validation of Entropy-Entanglement Link: The calculation $S_E = \ln(N(z))$ provides strong theoretical support for the idea that the entropy measured in high-energy scattering is a manifestation of quantum entanglement.
• Methodological Advance: The successful application of the homotopy approach to high-energy QCD evolution equations demonstrates a viable path for solving complex non-linear problems where exact analytical solutions are unavailable.
• KNO Scaling: The derivation of KNO scaling directly from the non-linear evolution equations in the saturation regime reinforces the universality of this scaling behavior in high-energy QCD.
• Future Work: The paper highlights the need for a rigorous matching procedure between the small-$n$ (homotopy) solution and the large-$n$ (analytical) solution to fully satisfy unitarity constraints across all multiplicities.

In summary, this paper advances the understanding of gluon multiplicity distributions in DIS by providing a unified framework that connects dipole evolution, AGK rules, and entropy calculations, offering both numerical precision for low multiplicities and analytical insight for the high-multiplicity tail.