Probing nuclear structure with the Balitsky-Kovchegov equation in full impact-parameter dependence

This paper extends the solution of the Balitsky-Kovchegov equation with full impact-parameter dependence to nuclear targets, providing predictions for key processes like deep-inelastic scattering and diffractive vector meson production to investigate gluon saturation and nuclear structure effects relevant for future facilities like the EIC and current LHC studies.

The Gist

The Big Picture: The "Crowded Room" Problem

Imagine you are trying to understand how a nucleus (the core of an atom) is built. Inside this nucleus, there are tiny particles called gluons that act like the "glue" holding everything together.

When you smash a proton (a single particle) or a heavy nucleus (like Lead or Gold) with a high-energy beam, you are essentially throwing a ball into a room full of people.

• The Proton is like a small, empty room with just a few people.
• The Heavy Nucleus is like a massive, packed concert hall.

The scientists in this paper are trying to predict what happens when you throw that ball into these rooms. Specifically, they are looking for a phenomenon called Gluon Saturation.

The Analogy:
Imagine the gluons are people in a room.

• Low Energy: If the room is empty, people can move freely. If you add more people (increase energy), the crowd grows linearly.
• High Energy (Saturation): Eventually, the room gets so packed that people start bumping into each other. They can't move freely anymore; they start merging or pushing back. The crowd density stops growing because there is simply no room left. This "traffic jam" is Gluon Saturation.

The Tool: The "Traffic Simulator" (The BK Equation)

To predict this, the authors use a mathematical tool called the Balitsky-Kovchegov (BK) equation. Think of this as a super-advanced traffic simulator.

Previously, this simulator could only handle the "Proton" (the small room). It was a bit like a 2D map that ignored where people were standing, only counting how many were there.

What this paper does:

1. Full 3D Simulation: They upgraded the simulator to include impact-parameter dependence. This means the simulator now knows exactly where in the room the people are standing, not just how many. It maps the density of the crowd from the center of the nucleus to the very edge.
2. From Proton to Nucleus: They took this upgraded simulator and applied it to heavy nuclei (like Oxygen, Gold, and Lead) to see how the "traffic jam" forms in a crowded concert hall compared to a small room.

The Experiment: Two Versions of Reality

To prove that "saturation" (the traffic jam) is real, the authors ran two different simulations:

1. The "Real World" Model (Non-linear BK): This model includes the rule that "people bump into each other and stop moving." This represents the physics of gluon saturation.
2. The "Fantasy World" Model (Linearized BK): This model removes the "bumping" rule. It assumes people can magically overlap and the crowd can grow infinitely without stopping. This is like a fantasy where the room expands as people enter.

The Result:
When they compared the two models against real data from the Large Hadron Collider (LHC):

• The "Real World" model matched the data perfectly.
• The "Fantasy World" model failed miserably, predicting that the crowd would grow too fast and behave differently than what we actually see.

This proves that the "traffic jam" (gluon saturation) is a real physical phenomenon, especially in heavy nuclei.

The Special Case: The Oxygen Tetrahedron

One of the most interesting parts of the paper is how they modeled Oxygen.

Usually, scientists model the nucleus as a smooth, fuzzy ball (like a cloud of gas). This is called the Woods-Saxon distribution.

However, the authors asked: What if Oxygen isn't a smooth ball, but a specific shape made of smaller blocks?

• Oxygen has 16 protons/neutrons.
• They modeled it as four Helium nuclei (alpha particles) arranged in a tetrahedron (a pyramid shape with a triangular base).

The Metaphor:
Imagine building a house.

• Standard Model: You use a giant bag of sand and pour it into a mold. It's smooth and round.
• Tetrahedral Model: You use four distinct bricks arranged in a pyramid.

They ran the simulation with this "pyramid" shape.

• The Finding: Surprisingly, for most measurements, the "pyramid" and the "smooth ball" look almost identical. The differences are too small to see with current tools.
• The Exception: If you look at the very edges of the collision (high momentum transfer), the "pyramid" creates a slightly different pattern of ripples (diffraction dips) than the smooth ball. It's like the difference between a smooth wave hitting a round rock vs. a jagged rock.

Why Does This Matter?

1. Future Factories (EIC): The authors are preparing for the Electron-Ion Collider (EIC), a new machine being built to study these collisions. They have provided a "menu" of predictions for what the EIC should see when it smashes different types of nuclei (Carbon, Oxygen, Gold, Lead).
2. Current Data (LHC): Their model already explains data coming out of the LHC right now, specifically regarding how heavy nuclei produce particles like the J/psi meson (a heavy particle made of a charm quark and an anti-charm quark).
3. The Smoking Gun: They identified a specific way to prove saturation exists: looking at how the production of these particles changes as you increase the energy. The "Real World" model predicts a drop in production at certain angles, while the "Fantasy" model predicts a continuous rise. This is the "smoking gun" experiment future facilities can perform.

Summary in a Nutshell

The authors took a complex math equation used to describe particle collisions, upgraded it to see the "shape" of the nucleus in 3D, and applied it to heavy atoms. They proved that inside heavy atoms, gluons get so crowded they form a "traffic jam" (saturation). They also tested if Oxygen is a smooth ball or a pyramid of smaller blocks, finding that while the shape is interesting, the "traffic jam" effect is the dominant story. This work helps physicists know exactly what to look for when they turn on the next generation of particle accelerators.

Technical Summary

1. Problem Statement

The paper addresses the need to understand gluon saturation and the internal structure of atomic nuclei, particularly in the context of future facilities like the Electron-Ion Collider (EIC) and current Large Hadron Collider (LHC) programs.

• The Gap: While the Balitsky-Kovchegov (BK) equation has been successfully used to describe gluon dynamics in proton targets with full impact-parameter dependence, its application to nuclear targets has been limited. Previous nuclear studies often relied on the Glauber model or simplified impact-parameter independent frameworks.
• The Challenge: Accurately modeling the transition from proton to nuclear targets requires handling the complex geometry of nuclei (e.g., light nuclei like Oxygen vs. heavy nuclei like Lead) and distinguishing between linear gluon evolution (BFKL) and non-linear saturation effects (BK) in a regime where experimental data is sparse or ambiguous.
• Specific Goal: To extend the recently solved, full impact-parameter dependent BK equation (originally for protons) to various nuclear targets, predict key observables, and identify experimental signatures that can distinguish saturation physics from linear evolution.

2. Methodology

The authors employ a phenomenological approach based on the Color Dipole Picture and the BK evolution equation.

• The BK Equation: They utilize the target-rapidity, collinearly-improved BK equation with full dependence on dipole size ($r$), impact parameter ($b$), and the relative angle ($\theta$).
  • The equation includes a non-linear term ($-\kappa N^2$) representing gluon recombination/saturation.
  • Linearization: To isolate saturation effects, they introduce a scalar regulator $\kappa$.
    • $\kappa = 1$: Standard non-linear BK equation (saturation active).
    • $\kappa = 0$: Linearized version (equivalent to BFKL evolution, no saturation).
• Nuclear Initial Conditions: The transition from proton to nucleus is achieved by modifying the initial conditions of the evolution:
  1. Saturation Scale: Scaled by the mass number $A$: $Q_{s0}^2 \to Q_{s0,A}^2 = 0.1 \text{ GeV}^2 A^{1/3}$.
  2. Profile Function: The Gaussian proton profile is replaced by nuclear thickness functions $T_A(b, r)$.
     • Standard Models: Woods-Saxon (2-point Fermi) for heavy nuclei; 3-point Fermi for Oxygen; Modified Harmonic Oscillator for Carbon.
     • Tetrahedral Oxygen Model: A novel approach for Oxygen ($^{16}\text{O}$) treating it as a cluster of four $\alpha$-particles (helium nuclei) arranged in a tetrahedron. This model accounts for the non-central density distribution of light nuclei.
• Observables Calculated:
  • Deep-Inelastic Scattering (DIS): Nuclear structure functions ($F_2$) and the nuclear modification factor ($R_{pA}$).
  • Vector Meson Production: Coherent diffractive photoproduction of $J/\psi$ mesons, calculating both total cross-sections and differential cross-sections ($d\sigma/d|t|$).

3. Key Contributions

1. Extension to Nuclear Targets: The first application of the full 3D (r, b, $\theta$) BK solution to a wide palette of nuclear targets (C, O, Ca, Fe, Cu, Au, Pb) without introducing new free parameters, relying solely on modified initial conditions.
2. Tetrahedral Oxygen Model: Implementation of an $\alpha$-cluster model for Oxygen, providing a Python library for generating these profiles. This allows for testing deviations from standard isotropic Woods-Saxon distributions.
3. Saturation Discrimination: A systematic comparison between the non-linear (saturation) and linearized (no saturation) BK equations to identify "smoking gun" observables for gluon saturation in nuclei.
4. Public Software: Release of the nuclear helper software (PyPI library) for generating tetrahedral and standard nuclear profiles.

4. Key Results

• Dipole Amplitudes:
  • In the non-linear case ($\kappa=1$), the amplitude evolution slows down (travelling wave behavior) and does not grow indefinitely.
  • In the linear case ($\kappa=0$), the amplitude grows steadily, leading to unphysically high cross-sections at high energies.
• Structure Functions ($F_2$) & $R_{pA}$:
  • Both models fit proton data well.
  • For nuclei, the non-linear model predicts significant suppression (shadowing) in $R_{pA}$, especially for heavy nuclei, due to gluon recombination.
  • The linear model fails to reproduce this suppression, showing enhancement instead. The difference grows with the mass number $A$.
• Vector Meson Production ($J/\psi$):
  • ALICE Data: The non-linear BK model successfully describes ALICE data for $J/\psi$ photoproduction on Lead (Pb). The linearized model shows a sizeable discrepancy.
  • Differential Cross-Sections ($d\sigma/d|t|$): The position of diffractive dips correlates with nuclear size.
  • Energy Dependence ($W$): A crucial finding is the behavior of the cross-section at small momentum transfer ($|t| \to 0$).
    • Non-linear: Predicts a drop in cross-section at small $|t|$ as energy increases (due to saturation).
    • Linear: Predicts a continuous rise at all $|t|$.
    • This differential behavior is identified as the most promising experimental signature to distinguish saturation from linear evolution.
• Oxygen Models:
  • The tetrahedral $\alpha$-cluster model shifts the transverse density maximum to the outer shell compared to the Woods-Saxon model.
  • However, this structural difference has a negligible effect on integrated observables. The only measurable difference appears in the differential $J/\psi$ cross-sections at large $|t|$, where the dip positions differ slightly.

5. Significance

• For Future Facilities (EIC/LHeC): The paper provides robust theoretical predictions for a wide range of nuclear targets, essential for planning experiments at the Electron-Ion Collider. It highlights that differential measurements (specifically $W$-dependence at fixed $|t|$) are critical for proving the existence of gluon saturation.
• For Current Physics (LHC): The results offer a framework to interpret current LHC data on ultra-peripheral collisions (UPC) and vector meson production, confirming that saturation effects are necessary to describe the data.
• Theoretical Advancement: By demonstrating that nuclear effects can be modeled purely through the dynamics of the BK equation with modified initial conditions (and no ad-hoc nuclear parameters), the work strengthens the universality of the saturation framework.
• Light Nuclei Structure: The study of the tetrahedral Oxygen model suggests that while exotic nuclear geometries exist, their impact on high-energy scattering observables is subtle, primarily visible only in high-momentum-transfer differential distributions.

In conclusion, this work establishes a comprehensive, parameter-free framework for studying nuclear gluon saturation, identifying specific differential observables in vector meson production as the key to experimentally validating the non-linear dynamics of QCD in nuclei.