⚛️ nuclear theory

Initial momentum anisotropies in the kT-factorization of the CGC I: Gradient Expansion

This paper demonstrates that $k_\perp$ -factorization in the Color Glass Condensate represents only the zeroth order of a systematic gradient expansion, revealing that higher-order corrections involving transverse derivatives of source distributions generate intrinsic initial-state momentum anisotropies and a more complete energy-momentum tensor that are missed by leading-order approximations.

Original authors: Oscar Garcia-Montero

Published 2026-09-09

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Oscar Garcia-Montero

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

When heavy atomic nuclei smash together at nearly the speed of light, they create a fleeting, super-hot soup of matter that existed only moments after the Big Bang. Scientists call this state of matter the quark-gluon plasma. To understand how this plasma behaves, researchers treat it like a fluid, using complex computer models to simulate how it flows and expands. However, these models need a starting point: a precise description of the energy and pressure inside the collision at the very first instant. For decades, the standard way to calculate this starting point has relied on a simplified assumption: that the energy is spread out evenly in all directions across the collision plane, like a flat pancake. This assumption has worked well enough to match many experimental results, but it ignores a crucial detail: the nuclei are not perfect, uniform spheres, and their internal structures vary from point to point.

A new study by Oscar Garcia-Monteroa challenges this long-standing simplification. The research focuses on the very first moments of a collision between two nuclei, specifically looking at how the tiny particles inside them—gluons—interact before the fluid even begins to form. The author demonstrates that the standard method, which assumes a smooth, even distribution of energy, is actually just the first, most basic step in a much more detailed calculation. By carefully accounting for the fact that the density of matter inside the nuclei changes from one spot to another, the study reveals that the initial state of the collision is not perfectly smooth. Instead, it possesses a built-in directional bias, or anisotropy, right from the start. This means the energy is already pushing harder in some directions than others before the fluid dynamics even begin to take over.

The paper achieves this by revisiting a theoretical framework known as the Color Glass Condensate, which describes the high-energy gluons inside fast-moving nuclei. The researcher performs a systematic expansion, a mathematical technique that adds layers of detail to a basic formula. In the simplest version, the nuclei are treated as if they are infinitely large and uniform, which leads to the standard result of even energy distribution. However, by including the next layers of detail—specifically, how the density of gluons changes across the face of the nucleus—the study uncovers a new set of effects. These effects arise because the gluons produced in the collision can "feel" the unevenness of the source they came from. The result is a more complete picture of the energy-momentum tensor, a mathematical object that describes how energy and pressure are distributed. This new picture includes terms that represent a flow of energy and a shear stress, which are forces that push matter sideways, creating an imbalance that the old, simpler models missed entirely.

To ensure these findings were not just an artifact of the calculation method, the author performed the analysis in two different ways. First, they looked at the collision as if it were a stream of individual particles, which is the traditional approach. Second, they treated the collision as a continuous, classical field, which is how the matter behaves in the earliest, most intense moments. Both methods led to the same conclusion: the standard formula is incomplete. The field-based calculation, in particular, revealed additional effects that arise from the interference between the waves of the collision. These interference patterns create a longitudinal flow of energy, a movement of energy along the direction of the collision, which was previously thought to be zero in the simplest models. This flow is a direct consequence of the real-time interaction between the two colliding nuclei, rather than a property of the nuclei themselves in isolation.

The significance of this work lies in its ability to explain how the geometry of the collision translates into the motion of the resulting matter. In the past, scientists assumed that any directional flow observed in the final debris was generated entirely by the fluid dynamics of the plasma as it expanded. This new study suggests that a portion of that flow is actually imprinted on the system at the very moment of creation. The effect is driven by the gradients, or slopes, in the density of the nuclei. Where the density changes sharply, the resulting pressure is not uniform. This mechanism is particularly important for collisions involving smaller nuclei, such as oxygen or neon, where fluctuations in density are more pronounced relative to the size of the system. In these smaller collisions, the initial directional bias is likely to be much stronger, potentially changing how quickly these small systems settle into a fluid-like state.

The study does not discard the old models but rather refines them. The new terms derived in the paper are local operators, meaning they can be added directly to existing computer simulations without requiring a complete overhaul of the code. This allows researchers to generate more realistic initial conditions for their simulations, incorporating a dynamically generated momentum anisotropy that was previously missing. The author notes that while the effect is present in all collisions, it is most significant in systems where the density profiles are steep and fluctuating. By including these corrections, future simulations will be able to test whether the observed flow in experiments is a true signature of the fluid's viscosity or if it is partly a remnant of the initial geometric imbalance.

Ultimately, this work provides a more faithful representation of the early universe's first moments. It shows that the transition from a collision of two nuclei to a flowing fluid is not a clean break between a static start and a dynamic evolution. Instead, the seeds of the fluid's motion are sown in the very first interactions, dictated by the uneven landscape of the colliding nuclei. The research bridges the gap between the microscopic physics of gluon production and the macroscopic behavior of the quark-gluon plasma, offering a clearer path to understanding the fundamental properties of matter under extreme conditions. The findings are presented as a systematic expansion, a ladder of corrections that can be climbed to achieve greater precision, with the first few rungs already revealing a richer, more complex reality than previously imagined.

Technical Summary: Initial Momentum Anisotropies in the kT-Factorization of the CGC I: Gradient Expansion

Problem Statement
Current frameworks for modeling heavy-ion collisions often delegate the conversion of spatial collision geometry into momentum-space anisotropy almost entirely to the hydrodynamic response. While models like IP-Glasma capture initial-state field dynamics, they are computationally expensive and face challenges in relating longitudinal current structures to small- $x$ nuclear structure. Conversely, tractable models based on the dilute-dilute limit of the Color Glass Condensate (CGC) utilize $k_\perp$ -factorization. However, the standard leading-order $k_\perp$ -factorized formula treats sources as locally homogeneous, resulting in an isotropic transverse momentum distribution at the initial stage. This obscures the extent to which momentum anisotropies are generated dynamically by the initial gluon production mechanism itself, before any hydrodynamic evolution. The paper addresses the question of how momentum anisotropies are handed over to hydrodynamics in the dilute-dilute limit of the CGC, specifically investigating whether the standard factorized formula misses intrinsic initial-state anisotropies.

Methodology
The author approaches the problem through two complementary methods to ensure robustness and capture different physical regimes:

Momentum-Space Gradient Expansion: Starting from the single-inclusive gluon production formula in the dilute-dilute limit, the author performs a systematic gradient expansion in the transverse positions of the two colliding sources. The expansion is organized in powers of $1/(Q_s R)$ , where $Q_s$ is the saturation scale and $R$ is the scale over which the source densities vary. By expanding the hard scattering kernel in powers of relative momenta (Fourier conjugates to spatial gradients), they derive an extended $k_\perp$ -factorized formula. This allows the construction of the energy-momentum tensor ( $T^{\mu\nu}$ ) directly from the gluon spectrum up to fourth order in gradients.
Real-Time Field Evolution: To capture infrared physics and coherent interference effects missed by the on-shell particle picture, the author performs a real-time computation of the classical Glasma fields in the future light cone. They solve the classical Yang-Mills equations perturbatively in the dilute-dilute limit and expand the resulting energy-momentum tensor around the "on-shell" limit (where amplitude and conjugate amplitudes share the same mode momenta). This method retains the coherent interference between field modes and their conjugates, revealing terms related to chromo-electric and chromo-magnetic interference.

Key Contributions and Results

Systematic Gradient Expansion: The paper demonstrates that the standard $k_\perp$ -factorization formula is merely the zeroth-order term of a systematic gradient expansion. Each higher order introduces tensor structures built from the local unintegrated gluon distributions (uGDFs) and their transverse derivatives.
Initial-State Anisotropies: The expansion reveals that momentum-space anisotropies (such as $T^{xx} - T^{yy} \neq 0$ and shear $T^{xy}$ ) are present in the initial state before hydrodynamic evolution. These arise from the coupling of the outgoing gluon momentum to the spatial gradients of the overlapping nuclear profiles. Even for isotropic uGDFs, finite impact parameters and non-vanishing density gradients generate these anisotropies.
Extended Factorization: The derived corrections retain the factorized form of the uGDFs but act as local differential operators on them. This means the corrections can be added directly to existing saturation-based initial-state models (e.g., McDipper) without abandoning the tractable factorized structure.
Field-Theoretic Corrections (Real-Time): The real-time computation recovers the momentum-space results in the eikonal limit but adds distinct terms:
- Longitudinal Energy Flux: Terms such as $T^{\tau\eta}$ and $T^{\eta i}$ , which strictly vanish in the on-shell yield computation, are found to be non-zero. These represent coherent longitudinal flow arising from the interference of the two nuclei's longitudinal fields.
- Chromo-Electric/Magnetic Interference: The expansion identifies terms driven by parity-odd invariants (related to $E_\perp \times B_\perp$ ) that are suppressed by powers of the large outgoing momentum $|p| \gg Q_s$ .
- Phase-Induced Flow: The interference phase between the amplitude and conjugate introduces a transverse, parity-odd term that generates a first-order momentum current (linear response to energy density gradients), a feature absent in the pure particle picture.

Significance and Claims
The paper claims that the local isotropy of the transverse stress in standard $k_\perp$ -factorized initial conditions is an artifact of the leading approximation, not a physical necessity. By retaining transverse gradients, the author shows that momentum-space anisotropy is dynamically generated at the level of single-gluon production.

Universality: These anisotropies are a property of the dynamical mechanism (the real-time interaction of colliding fields) and do not depend on the specific model chosen for the uGDFs, although the magnitude of the effect depends on the model's density profiles.
Phenomenological Relevance: The effect is controlled by density gradients and is expected to be most pronounced in fluctuation-dominated, dilute systems, such as light-ion collisions (e.g., O+O, Ne+Ne). The author suggests these initial anisotropies could influence how small systems approach a hydrodynamic regime and affect the extraction of transport coefficients from flow data.
Practical Implementation: A key practical contribution is that these corrections are local operators on the same unintegrated distributions used in current models. This allows for the generation of realistic event-by-event profiles with dynamically generated initial momentum anisotropies without requiring full real-time lattice simulations.

The author notes that while similar results regarding longitudinal flux were obtained in previous small- $\tau$ expansions of the Glasma, the origin of the gradients in this work is distinct: they arise from the real-time interactions resolving the spatial inhomogeneity of the incoming sources, rather than from the breaking of translation invariance within the nuclei itself. The phenomenological consequences of these findings are reserved for a companion paper.

More like this