IP-Glasma Framework

• IP-Glasma is a QCD-based framework that models the initial state in ultra-relativistic heavy-ion collisions using CGC theory and dynamic Yang–Mills evolution.
• It integrates IP-Sat constraints with event-by-event fluctuations to generate realistic energy density profiles and geometric eccentricities for hydrodynamic simulations.
• The framework accurately predicts observables such as gluon multiplicity and flow harmonics, enhancing quantitative QGP phenomenology.

The IP-Glasma framework is a first-principles initial-state model for ultra-relativistic heavy-ion collisions, combining the Color Glass Condensate (CGC) effective theory of high-density QCD matter with an impact parameter–dependent saturation (IP-Sat) model for parton distributions, and event-by-event classical Yang–Mills (CYM) dynamics of gluonic fields. It provides a comprehensive description of the early-time, pre-equilibrium evolution of strongly interacting QCD matter, generating lumpy, fluctuating distributions of energy, momentum, and spatial geometry that serve as event-by-event fluctuating initial conditions for viscous relativistic hydrodynamics. The framework has systematically evolved from a 2+1D (boost-invariant) realization to a full 3+1D (rapidity-dependent) formulation, incorporating state-of-the-art JIMWLK small-$x$ QCD evolution, and underpins much of the current quantitative phenomenology of QGP collectivity, flow harmonics, and early-time dynamics.

1. Theoretical Foundations: CGC, IP-Sat, and Wilson Lines

The IP-Glasma framework is rooted in the CGC effective theory, valid at small Bjorken-$x$, where high gluon occupancy allows the treatment of gluon fields as classical color fields generated by stochastic large-$x$ color charges $\rho^a(x^\mp, \mathbf{x}_\perp)$. The statistical distribution $W_Y[\rho]$ of these charges is evolved in rapidity $Y$ via the JIMWLK equation, implementing quantum small-$x$ evolution. The local color-charge squared fluctuation $\mu^2(\mathbf{x}_\perp)$ is tied to the local saturation scale $Q_s$.

The IP-Sat model provides constraints on $Q_s(x,b)$ from HERA DIS data, parameterizing the dipole–nucleus forward amplitude as

$x$0

with $x$1 the nucleon transverse profile and $x$2 the gluon density.

Wilson lines,

$x$3

encode the eikonal color rotation of high-energy partons traversing the nucleus. The IR regulator $x$4 ensures color neutrality.

2. Event-by-Event Fluctuations and Initial Condition Sampling

Each event's geometry is constructed by sampling nucleon positions from a Woods–Saxon distribution, followed by sampling local color charges $x$5 as independent Gaussian random variables with variance determined by the local $x$6 extracted from IP-Sat. Sub-nucleonic fluctuations (“hot spots”) and rapidity correlations are also implemented where required.

The sampled color charge densities are exponentiated to define Wilson lines for each nucleus, typically discretized into $x$7 longitudinal sheets (with $x$8 ranging from 10–50). The Gaussian two-point correlator,

$x$9

determines the quantum fluctuations at the sub-nucleonic scale.

For 3D (rapidity-dependent) implementations, Wilson lines are further evolved in rapidity using the JIMWLK equation, with stochastic Langevin updates incorporating Gaussian noise and the appropriate kernel regulated by $x$0.

3. Classical Yang–Mills Evolution and Glasma Field Dynamics

Following initialization, the classical (sourceless) Yang–Mills equations,

$x$1

are solved numerically on a transverse (2+1D) or full three-dimensional (3+1D, in $x$2, $x$3, $x$4, $x$5) lattice. Typical parameters are $x$6–$x$7, lattice spacing $x$8–$x$9 fm, and infrared regulator $\rho^a(x^\mp, \mathbf{x}_\perp)$0–$\rho^a(x^\mp, \mathbf{x}_\perp)$1 GeV.

In 2+1D (boost-invariant) simulations, the initial gauge fields $\rho^a(x^\mp, \mathbf{x}_\perp)$2 and $\rho^a(x^\mp, \mathbf{x}_\perp)$3 are evolved with vanishing $\rho^a(x^\mp, \mathbf{x}_\perp)$4 gradients. In 3+1D realizations, genuine rapidity dependence is implemented by constructing pure-gauge pre-collision fields in all three spatial directions. The Gauss law constraint $\rho^a(x^\mp, \mathbf{x}_\perp)$5 is solved numerically for the initial transverse electric field, typically via a covariant Poisson equation for a scalar potential $\rho^a(x^\mp, \mathbf{x}_\perp)$6.

Time stepping is accomplished via leapfrog integration, updating electric and magnetic fields and the link variables using the lattice Hamiltonian equations in temporal or Fock–Schwinger gauge.

4. Observables: Multiplicities, Geometry, and Event-by-Event Fluctuations

Key event-by-event observables computed in the IP-Glasma framework include:

• Gluon multiplicity: The gluon spectrum is extracted from equal-time field correlators in Coulomb gauge after CYM evolution,

$\rho^a(x^\mp, \mathbf{x}_\perp)$7

The event-by-event multiplicity distribution is well described by a negative binomial distribution, with parameters directly determined by the fluctuating initial Glasma fields. For instance, $\rho^a(x^\mp, \mathbf{x}_\perp)$8 in central collisions, decreasing in more peripheral events.

• Energy density: The lattice energy density at time $\rho^a(x^\mp, \mathbf{x}_\perp)$9 reads

$W_Y[\rho]$0

which exhibits spatial granularity on scales $W_Y[\rho]$1.

• Eccentricities: Energy-weighted eccentricities up to $W_Y[\rho]$2,

$W_Y[\rho]$3

quantify the initial geometry and drive final-state harmonic flow, with fluctuation-dominated odd moments $W_Y[\rho]$4 systematically larger than in MC-KLN or wounded-nucleon models.

• Longitudinal structure (3D): 3+1D implementations generate rapidity-dependent profiles, enabling computation of rapidity correlations, two-point $W_Y[\rho]$5, and harmonic decorrelation ratios $W_Y[\rho]$6.

5. Coupling to Hydrodynamics and Hadronic Cascade

The energy-momentum tensor $W_Y[\rho]$7 of the Glasma fields at a switching time $W_Y[\rho]$8 (typically $W_Y[\rho]$9–$Y$0 fm/c) is diagonalized (“Landau matching”) to extract the local energy density $Y$1 and flow velocity $Y$2. These serve as initial conditions for event-by-event relativistic viscous hydrodynamics, e.g., the MUSIC code, which is run with constant or temperature-dependent shear viscosity $Y$3 (e.g., $Y$4) and bulk viscosity $Y$5 tuned to single-particle spectra and flow harmonics.

Hydrodynamic evolution proceeds until the fluid reaches a switching temperature (e.g., $Y$6 MeV), after which a hadronization prescription (e.g., Cooper–Frye) generates particle distributions. The subsequent hadronic rescattering and resonance decays are simulated using a microscopic hadronic transport model such as UrQMD.

6. Key Phenomenological Insights

The IP-Glasma framework, when interfaced with hydrodynamics and hadronic afterburners, correctly predicts:

• The centrality and system dependence of charged-hadron multiplicities at both RHIC and LHC energies, with the inclusion of JIMWLK evolution necessary for reproducing the $Y$7-dependence from RHIC to LHC (Mäntysaari et al., 28 Aug 2025).
• Realistic distributions of flow harmonics $Y$8 up to $Y$9 and their fluctuations, due to large initial-state eccentricities. The ratio $x$0–$x$1 for mid-central Au+Au, close to experiment (Schenke et al., 2012).
• A more compact initial transverse energy density ($x$2) than wounded-nucleon (MC-Glauber) models, leading to stronger radial flow and requiring larger bulk viscosity $x$3 to match final $x$4 (Schenke et al., 2018).
• Nontrivial initial flow fields in both transverse and longitudinal directions, with pre-equilibrium longitudinal flow $x$5 of magnitude $x$6, not captured in purely boost-invariant (2+1D) treatments (McDonald et al., 2017).
• Genuine rapidity-dependent correlations and flow-vector decorrelation observables $x$7 in $x$8D, controlled directly by JIMWLK evolution, agreeing with CMS factorization-measurements (McDonald et al., 2020, Schenke et al., 2017).

7. Extensions and Ongoing Developments

Recent developments and extensions include:

• Efficient numerical methods for $x$9D CYM evolution in Milne coordinates with realistic treatment of color currents, boundary conditions, and Gauss-law enforcement (Matsuda et al., 2023).
• Fully event-by-event hydrodynamic simulations initialized by $\mu^2(\mathbf{x}_\perp)$0D IP-Glasma, facilitating systematic exploration of longitudinal observables and the impact of initial-state longitudinal fluctuations (McDonald et al., 2023).
• Incorporation of heavy-quark transport (e.g., charm) initialized and propagated through the Glasma and QGP using the MARTINI event generator, probing early-stage charm–medium interactions (Gale et al., 2024).

Compared to more parametric or geometric models (e.g., TRENTo, MC-Glauber, EKRT), the IP-Glasma framework is distinguished by its direct implementation of small-$\mu^2(\mathbf{x}_\perp)$1 QCD evolution, event-by-event geometry fluctuations, and a principled mapping to hydro and hadronic observables. The predictive power and flexibility of the framework have made it a standard tool for quantitative QGP phenomenology and for interpreting the initial conditions in relativistic heavy-ion collisions.