When JIMWLK evolution really matters: the example of incoherent diffraction

This paper demonstrates that the Gaussian Approximation, while accurate for two-gluon exchange processes, fails to describe incoherent diffraction in photon-nucleus collisions because it significantly underestimates cross sections compared to the full JIMWLK evolution for correlators involving four-gluon exchanges.

The Gist

Imagine you are trying to predict how a crowd of people (representing subatomic particles called gluons) behaves when a single person (a photon) runs through them at nearly the speed of light. This happens in high-energy physics experiments, like those at the Large Hadron Collider.

Physicists have a very complex, precise rulebook for how this crowd moves and interacts, called the JIMWLK equation. It's like a super-accurate, computer-simulated weather forecast that tracks every single gust of wind and shift in the crowd's mood. However, running this simulation is incredibly difficult and slow, like trying to calculate the path of every single raindrop in a storm.

To make things easier, scientists often use a shortcut called the Gaussian Approximation (GA). Think of this as using a simple, smooth average to describe the crowd. Instead of tracking every individual, you just say, "On average, the crowd is moving this way." For many situations, this shortcut works amazingly well. It's like saying, "The average temperature is 70°F," which is a great guess for a sunny afternoon.

The Problem: When the Shortcut Fails

This paper asks a critical question: Does this shortcut always work?

The authors discovered that the shortcut fails spectacularly in a specific scenario: Incoherent Diffraction.

To understand this, imagine the crowd isn't just a smooth mass, but a group of people holding hands in a complex, shifting web.

• Coherent Diffraction (The Shortcut Works): If the crowd moves together as one big block, the "average" description works fine. The shortcut predicts the outcome correctly.
• Incoherent Diffraction (The Shortcut Fails): This happens when the crowd breaks apart into smaller, chaotic groups that move independently. The paper shows that in this chaotic state, the "average" description (the Gaussian Approximation) completely misses the mark. It's like trying to predict the behavior of a chaotic mosh pit by looking at the average movement of a calm line of people. The shortcut assumes the crowd is too smooth and orderly, ignoring the wild, individual fluctuations that actually drive the result.

The Analogy of the Four-Handshake

The paper explains that the shortcut works well when the interaction involves a simple "two-handshake" between particles. It's like two people shaking hands; the average rule covers it.

However, the "Incoherent Diffraction" scenario involves a complex "four-handshake" (a four-gluon exchange). Imagine four people trying to coordinate a dance move. The shortcut assumes they are just doing a simple, average dance. But in reality, they are doing a complex, synchronized routine that depends on their specific, individual positions. The shortcut misses these specific, complex connections, leading to a wrong prediction.

What the Authors Did

1. The Math Check: They did the math on a single step of the process and proved that the shortcut gives a different answer than the precise rulebook. Specifically, the shortcut predicted the result would be zero in certain geometric arrangements, while the precise rulebook showed it would be significant.
2. The Computer Simulation: They ran massive computer simulations using the precise rulebook (JIMWLK) and compared them to the shortcut (GA).
3. The Result: The precise rulebook consistently predicted much larger effects (cross-sections) than the shortcut did. In some cases, the shortcut was off by a factor of two.

The Bottom Line

The paper concludes that while the "average" shortcut (Gaussian Approximation) is a useful tool for many physics problems, it is dangerous to use when studying "incoherent diffraction" (where the target breaks apart or fluctuates wildly). In these cases, you cannot rely on the average; you must use the full, complex, and computationally expensive rulebook (JIMWLK) to get the right answer.

The authors emphasize that for these specific types of collisions, the "fluctuations" (the individual quirks of the crowd) are the most important part of the story, and the shortcut simply smooths them out too much, hiding the real physics.

Technical Summary

1. Problem Statement

The paper addresses a critical question in high-energy Quantum Chromodynamics (QCD) within the Color Glass Condensate (CGC) effective theory: Is the Gaussian Approximation (GA) a valid substitute for the full JIMWLK evolution equation for all observables?

• Context: The JIMWLK equation describes the energy evolution of Wilson line correlators, which encode the scattering of a dilute probe (like a photon) off a dense gluon target (a nucleus). Solving the full JIMWLK equation is computationally expensive.
• Current Practice: Phenomenological studies often use the Gaussian Approximation (GA), which assumes the target color fields follow a Gaussian distribution. The GA is known to be highly accurate for observables where the weak-scattering expansion begins with a two-gluon exchange (e.g., inclusive cross-sections, coherent diffraction).
• The Gap: It is unclear if the GA remains accurate for observables where the leading-order perturbative expansion involves four-gluon exchanges. The authors hypothesize that for such observables, the GA fails because it cannot capture the specific non-Gaussian correlations generated by JIMWLK evolution.
• Target Observable: The paper focuses on incoherent diffraction in photon-nucleus collisions (specifically diffractive dijet production). In this process, the target nucleus breaks up, and the cross-section is proportional to the variance of the scattering amplitude, $\langle T^2 \rangle - \langle T \rangle^2$. This quantity starts at order $\alpha_s^4$ (four-gluon exchange) in the weak scattering limit.

2. Methodology

The authors employ a dual approach, combining analytical derivations in the weak scattering limit with full numerical simulations.

A. Analytical Approach (Weak Scattering Limit)

• Setup: They analyze the evolution of the incoherent diffractive correlator $W_{xy\bar{y}\bar{x}} = \langle D_{xy}D_{\bar{y}\bar{x}} \rangle - \langle D_{xy} \rangle \langle D_{\bar{y}\bar{x}} \rangle$ over a single rapidity step ($\Delta Y$).
• Comparison:
  1. GA Prediction: They evolve the correlator assuming the Gaussian Hamiltonian ($H_{GA}$), which expresses higher-point functions purely in terms of the dipole amplitude $T$.
  2. JIMWLK Prediction: They use the full Balitsky hierarchy equations (specifically the equations for the dipole and the double-dipole) to evolve the same correlator. This includes contributions from "non-dipole" terms (sextupoles) that arise from color exchanges between the two dipoles.
• Model: They utilize the Golec-Biernat-Wusthoff (GBW) model for the initial dipole amplitude to perform explicit integrations over the transverse momentum of the emitted soft gluon.

B. Numerical Approach (Full Kinematics)

• Lattice Simulation: They solve the JIMWLK equation numerically using the Langevin formulation.
• Initial Condition: The simulation starts with the McLerran-Venugopalan (MV) model, which is Gaussian by construction (ensuring the GA is valid at $Y_0$).
• Evolution: They evolve the system by $\Delta Y \approx 5.18$ units using a running coupling prescription.
• Observables: They calculate the connected double-dipole correlator $W$ and the sextupole correlator $S$ for various geometric configurations (e.g., collinear, right-angled, square) and compare them directly to the GA results derived from the same evolved dipole amplitude.
• Cross-Section Calculation: Finally, they compute the ratio of incoherent to coherent diffractive cross-sections to quantify the phenomenological impact.

3. Key Contributions and Results

A. Analytical Breakdown of the GA

The authors prove analytically that the GA fails for incoherent diffraction even in the weak scattering limit:

• Functional Mismatch: The functional form of the evolution derived from the GA (Eq. 5.2) differs from the exact JIMWLK evolution (Eqs. 5.7 and 5.8).
• Angular Dependence: For a specific geometric configuration (three charges forming a triangle), the GA predicts a vanishing evolution rate when the angle $\theta = \pi/2$ (charges on a right triangle). In contrast, the full JIMWLK evolution yields a non-zero result.
• Magnitude: The JIMWLK evolution rate is systematically larger than the GA prediction. For averaged configurations, the JIMWLK result is approximately 50% larger than the GA result.
• Eigenvalue Analysis: They define an effective eigenvalue $\lambda_{eff}$ for the evolution. They find $\lambda_{eff} \approx 1.4 - 1.6 \times \lambda_{GA}$, indicating that the growth of incoherent fluctuations is faster than predicted by the Gaussian approximation.

B. Numerical Verification

The numerical lattice simulations confirm the analytical findings across a wide kinematic range:

• Systematic Underestimation: The GA consistently underestimates the incoherent correlator $W$ compared to the full JIMWLK solution.
• Magnitude of Deviation: The discrepancy grows with rapidity. In the weak scattering regime, the difference is large. Even in the saturation regime (distances $r \sim 1/Q_s$), the JIMWLK result remains roughly 50% larger than the GA result.
• Geometric Sensitivity: The deviation depends on the relative geometry of the Wilson lines. For "crossing" dipole configurations where the GA predicts zero, the JIMWLK evolution generates significant non-zero values ($\sim 0.02$).
• Sextupole Correlators: The GA also fails to accurately describe sextupole correlators (which appear in the evolution of $W$), showing deviations of tens of percent.

C. Phenomenological Impact

• Cross Sections: The ratio of incoherent to coherent diffractive cross-sections ($\sigma_{inc}/\sigma_{coh}$) calculated using JIMWLK is significantly higher than that calculated using GA.
• Conclusion: The GA underestimates the contribution of color charge fluctuations to the incoherent cross-section. This implies that previous phenomenological studies relying on GA for incoherent diffraction may have underestimated the magnitude of these events, particularly at high momentum transfer $|t|$.

4. Significance

• Theoretical Validity: The paper establishes a clear boundary for the validity of the Gaussian Approximation. While GA is excellent for two-gluon exchange processes, it is invalid for observables dominated by four-gluon exchanges (like incoherent diffraction).
• Necessity of Full JIMWLK: The authors demonstrate that for incoherent diffraction, one cannot rely on the computationally cheaper GA. Full JIMWLK evolution (via Langevin equations) is inevitable for accurate predictions.
• Future Experiments: The results are highly relevant for the upcoming Electron-Ion Collider (EIC) and LHC ultra-peripheral collisions. Incoherent diffraction is a key observable for probing proton/nucleus structure and fluctuations. Using the GA for these measurements would lead to incorrect interpretations of the data regarding the saturation scale and fluctuation dynamics.
• New Physics: The discrepancy highlights the importance of non-Gaussian correlations (specifically those involving color exchanges between different parts of the projectile) that are missed by mean-field approximations.

In summary, this paper provides a rigorous proof that the Gaussian Approximation is insufficient for describing incoherent diffraction, necessitating the use of full JIMWLK evolution to correctly capture the dynamics of color charge fluctuations in high-energy QCD.