⚛️ lattice

Unitarity dressing of the dynamical gluon mass scale

This paper demonstrates that enforcing $s$ -channel unitarity in high-energy elastic $pp$ scattering models significantly increases the extracted dynamical gluon mass scale from approximately 300–420 MeV to 750–1100 MeV, revealing that the physical infrared scale is a unitarity-dressed quantity determined jointly by the nonperturbative gluon propagator and multiple-exchange dynamics.

Original authors: T. V. Iser, V. Li, E. G. S. Luna

Published 2026-09-11

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: T. V. Iser, V. Li, E. G. S. Luna

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

High-energy physics often feels like trying to understand a storm by looking at a single raindrop. At the heart of this storm lies the strong nuclear force, the invisible glue that binds the fundamental particles of matter together. This force is carried by particles called gluons. In the simplest mathematical descriptions of particle collisions, these gluons are treated as weightless, massless entities. However, when scientists look at how these particles behave in the real world, particularly at very low energies or over short distances, the picture changes. Theoretical work and massive computer simulations suggest that gluons acquire a kind of effective mass, a resistance to change that emerges from the complex interactions of the force itself. This "dynamical mass" is not a fixed weight like a brick, but a property that arises from the environment of the collision. Determining exactly how heavy this mass is, and how it influences the way particles scatter off one another, is a crucial puzzle for understanding the structure of matter.

A team of researchers in Brazil has taken a fresh look at this puzzle by examining data from the world's most powerful particle accelerator. They focused on the elastic scattering of protons, a process where two protons collide and bounce off each other without breaking apart. In these collisions, the protons exchange a force carrier known as a Pomeron, which in this specific theoretical model is represented by a pair of gluons working together. For years, scientists analyzing these collisions at the "Born level"—a term for the simplest, first-pass calculation—have consistently found that the gluon mass needed to match the data falls in a range of a few hundred million electron volts. This value, roughly between 300 and 400 MeV, has become the standard expectation, fitting comfortably with other measurements of how the strong force behaves at low energies.

The researchers in this study asked a critical question: is this familiar mass value the true, intrinsic weight of the gluon pair, or is it an illusion created by the way we analyze the data? They realized that at the incredibly high energies of the Large Hadron Collider, the simple picture is incomplete. When protons collide at these speeds, they do not just exchange a single pair of gluons; they undergo a complex series of multiple exchanges and interactions. These additional interactions, governed by a fundamental rule of physics called unitarity, act like a dressing or a coating on the basic collision process. The team wanted to see what happened to the estimated gluon mass when they stopped looking at the collision as a single event and instead modeled it as a full, multi-step interaction where these extra exchanges were included.

To test this, the team took the same basic input—the theoretical model of the two-gluon exchange—and ran it through two different, sophisticated mathematical frameworks designed to account for these multiple interactions. They compared the results against real experimental data collected by the ATLAS and TOTEM collaborations at the Large Hadron Collider. The data covered proton collisions at three different energy levels: 7, 8, and 13 tera-electron volts. The researchers were careful to treat the data from these two different experiments separately, as they are known to have slight differences in how they measure the overall scale of the collisions, ensuring that their findings were not just an artifact of averaging two slightly different datasets.

The results were striking and consistent. When the researchers applied the complex unitarity corrections to their model, the value of the gluon mass required to fit the data jumped dramatically. Instead of the familiar range of 300 to 400 MeV, the new, more complete models demanded a mass scale between roughly 750 and 1000 MeV. This represents an increase by a factor of about 2.5. This shift was not a fluke; it appeared regardless of which of the two mathematical frameworks they used, regardless of which specific energy dataset they analyzed, and even when they tried two different mathematical formulas for how the gluon mass changes with energy. The consistency of this jump across all scenarios suggests that the effect is robust and real.

The physical reason for this shift lies in how the collision is visualized. In the simple, single-exchange model, the gluon pair must do all the heavy lifting to explain the observed spread and shape of the collision data. To match the wide, gentle curve seen in the experiments, the model needs a lighter, more spread-out gluon mass. However, once the researchers added the multiple exchanges, the picture changed. These additional interactions naturally contribute to the spread and shape of the collision profile. Because the multiple exchanges help build the physical shape of the collision, the initial, basic gluon pair no longer needs to be as spread out to match the data. It can be more compact. In the language of this theory, a more compact gluon pair corresponds to a heavier mass. Therefore, the mass value that fits the data is not just the mass of the gluon in isolation, but a "dressed" mass that includes the effects of the surrounding quantum environment.

This finding clarifies a long-standing ambiguity in how we interpret particle physics data. The mass scale extracted from high-energy elastic scattering is not a direct measurement of an isolated particle property, but rather a value that is fixed jointly by the fundamental nature of the gluon and the complex dynamics of the collision itself. The study confirms that to truly understand the infrared scale of the strong force, one cannot ignore the nonlinear mapping between the elementary building blocks and the final physical outcome. The "unitarity-dressed" scale of roughly 750 to 1000 MeV appears to be the more accurate reflection of the physics at play in these high-energy collisions, bridging the gap between theoretical predictions of the gluon mass gap and the experimental observations of how protons scatter.

Technical Summary: Unitarity Dressing of the Dynamical Gluon Mass Scale

Problem Statement
The paper addresses a discrepancy in the determination of the dynamical gluon mass scale ( $m_g$ ) within Quantum Chromodynamics (QCD). Nonperturbative analyses, such as Schwinger-Dyson equations and lattice simulations, suggest a finite infrared gluon propagator characterized by a momentum-dependent dynamical mass. Phenomenological applications in hadronic scattering, specifically using Landshoff–Nachtmann (LN) models where the soft Pomeron is a color-singlet pair of massive gluons, typically extract $m_g$ values in the range of a few hundred MeV (e.g., 299–422 MeV) when fitting high-energy elastic $pp$ scattering data at the Born (tree-level) approximation.

However, at LHC energies, the identification of the elementary one-Pomeron input with the physical elastic amplitude is problematic. $s$ -channel unitarity requires nonlinear corrections and multiple exchanges, which relate the elementary input to the physical impact-parameter profile in a complex manner. The central question is whether the "few-hundred-MeV" scale is an intrinsic property of the elementary two-gluon exchange or an effective scale resulting from unitarity dressing (rescattering effects).

Methodology
The authors investigate this by embedding a specific Reggeized LN two-gluon exchange amplitude into two standard one-channel unitarization schemes: the eikonal and the U-matrix formalisms.

Elementary Input: The Born amplitude ( $M_B$ ) is constructed using a Reggeized two-gluon exchange. The soft Pomeron trajectory is defined as $\alpha_P(t) = 1 + \epsilon + \alpha'_P t$ . The nonperturbative QCD content is introduced via an infrared-finite gluon propagator $D(k^2)$ and an effective charge $\bar{\alpha}(k^2)$ .
Mass Solutions: Two distinct solutions for the running gluon mass $m(Q^2)$ $m (Q^{2})$ derived from Schwinger-Dyson analyses are tested:
- A logarithmic solution ( $m_{log}$ ).
- A power-law solution ( $m_{pl}$ ).
  In both cases, the infrared scale $m_g = m(0)$ is the primary fitted parameter.
Unitarization: The elementary amplitude is converted into impact-parameter space ( $\chi(s, b)$ $χ (s, b)$ ) and then unitarized:
- Eikonal scheme: $A_E(s, b) = i[1 - \exp(i\chi)]$ .
- U-matrix scheme: $A_U(s, b) = \chi / (1 - i\chi/2)$ .
  Both schemes reproduce the Born term and the leading nonlinear rescattering correction but differ in higher-order multiple exchanges and high-opacity limits.
Data Sets: The analysis fits forward elastic $pp$ differential cross-section data ( $d\sigma/dt$ ) at $\sqrt{s} = 7, 8,$ and $13$ TeV. To avoid normalization tensions between experiments, two separate ensembles are used: Ensemble A (ATLAS data) and Ensemble T (TOTEM data). The fits are restricted to the forward region ( $|t| \le 0.11$ GeV $^2$ ) to ensure the validity of the minimal proton form factor.

Key Results
The study yields a clear and robust shift in the extracted gluon mass scale depending on the treatment of unitarity:

Born Level: When the Reggeized two-gluon exchange is fitted without unitarization, the extracted mass scales fall within the conventional phenomenological window:
- Logarithmic mass: $m_g \approx 299–355$ MeV.
- Power-law mass: $m_g \approx 363–422$ MeV.
Unitarized Level: When the same input is embedded in the eikonal or U-matrix schemes, the preferred mass scale increases significantly:
- Eikonal: $m_g \approx 779–1101$ MeV.
- U-matrix: $m_g \approx 749–1033$ MeV.
Enhancement Factor: The ratio of the unitarized mass to the Born mass ( $R = m_g^{unit} / m_g^{Born}$ ) is consistently close to 2.5 across all data ensembles (ATLAS and TOTEM) and both mass running solutions.
Stability: This enhancement is stable against:
- The choice of experimental data set (ATLAS vs. TOTEM).
- The analytic form of the running mass (logarithmic vs. power-law).
- The specific unitarization prescription (eikonal vs. U-matrix).
- The variation of other fit parameters ( $\epsilon$ and the proton form factor slope $a_1$ ), which adjust naturally without undergoing a comparable qualitative shift.

Physical Interpretation
The authors attribute this shift to the nonlinear mapping between the elementary two-gluon kernel and the physical impact-parameter profile.

At the Born level, the single exchange must independently reproduce the observed transverse extension and normalization of the data. A smaller $m_g$ provides a longer-ranged, stronger infrared kernel, which mimics the physical profile's width.
With unitarization, multiple exchanges (rescattering) contribute to the transverse structure and curvature of the amplitude. Consequently, the elementary kernel required to fit the data can be more compact (shorter-ranged). Since the correlation length of the kernel scales as $\ell_g \sim 1/m_g$ , a more compact kernel corresponds to a larger fitted $m_g$ .

Significance and Claims
The paper concludes that the dynamical gluon mass scale extracted from high-energy elastic scattering should not be viewed as the mass of an isolated one-Pomeron exchange. Instead, it is a "unitarity-dressed infrared scale."

Nature of the Scale: The value is fixed jointly by the nonperturbative gluon propagator and the many-body dynamics required by $s$ -channel unitarity.
Comparison with Lattice/Functional Methods: The authors note that the unitarity-dressed interval ( $\sim 0.75–1.1$ GeV) overlaps with the lower end of the "gluon mass gap" ( $\sim 0.7–0.8$ GeV) found in recent functional and lattice analyses of the Landau-gauge propagator. They suggest that elastic diffraction, once unitarity is implemented, may probe the same intermediate nonperturbative momentum domain associated with the mass gap, rather than just the deep infrared saturation scale.
Limitations: The authors explicitly state that their one-channel model is not intended to describe the entire pre-dip region or the full $t$ -dependence beyond the forward cone. The rapid deterioration of fit quality when extending $|t|$ beyond 0.11 GeV $^2$ indicates the need for additional transverse structure or diffractive-channel effects in that region. However, within the forward cone, the conclusion regarding the necessity of unitarity dressing for quantitative extraction of $m_g$ remains robust.

Technical Summary: Unitarity Dressing of the Dynamical Gluon Mass Scale

More like this