IP-Glasma+MUSIC+UrQMD Framework

• The IP-Glasma+MUSIC+UrQMD framework is a hybrid, event-by-event simulation pipeline that integrates QCD-based initial conditions, relativistic viscous hydrodynamics, and hadronic transport.
• It employs dynamic Landau matching and advanced particlization techniques to accurately transfer the early Yang–Mills evolution into hydrodynamic fields and subsequent microscopic transport.
• The modular approach synthesizes diverse observables, including event-plane correlations and nonlinear mode-coupling coefficients, thereby enabling precise constraints on QGP shear and bulk viscosities.

Searching arXiv for the specific framework paper and closely related follow-up work to ground the article in current and historical literature. {"query":"IP-Glasma MUSIC UrQMD framework arXiv (McDonald et al., 2017) event-plane correlations nonlinear response coefficients", "max_results": 10} {"query":"arXiv IP-Glasma+MUSIC+UrQMD flow observables 5.02 TeV high-order anisotropic flow nonlinear hydrodynamic responses", "max_results": 10} The IP-Glasma+MUSIC+UrQMD framework is a hybrid, event-by-event pipeline designed to connect QCD-based initial-state dynamics to hydrodynamic evolution and late-stage hadronic transport, and to synthesize a broad set of flow observables for precision constraints on QGP transport properties and initial-state fluctuations. In the 2.76 TeV Pb–Pb study that established its detailed use as an overview framework, it was extended beyond single-particle spectra and anisotropic flow coefficients $v_n$ to include event-plane correlations, non-linear mode-coupling coefficients $\chi_{npq}$, and event-shape engineering, with the aim of constraining the fluctuation spectrum, linear and non-linear hydrodynamic response, and the sensitivity of the data to $\eta/s$ and $\zeta/s(T)$ (McDonald et al., 2017).

1. Architecture of the hybrid pipeline

The framework is organized as a sequence of dynamically matched stages. The initial energy-momentum tensor is constructed event-by-event from the IP-Glasma model, which combines realistic nuclear geometry with sub-nucleonic color charge fluctuations. Color sources are sampled and evolved in classical Yang–Mills fields, yielding gluon fields whose energy density and flow are computed at early time. These fields encode the saturation physics through the local saturation scale and produce fluctuating “Glasma” flux tubes whose transverse structure drives the initial eccentricities. The pre-equilibrium evolution is purely classical Yang–Mills up to the hydrodynamization time $\tau_0$, with no additional modeling inserted before hydrodynamics (McDonald et al., 2017).

At $\tau_0$, Landau matching converts the IP-Glasma energy-momentum tensor $T^{\mu\nu}(\tau_0)$ to hydrodynamic fields—energy density $e$, flow velocity $u^\mu$, and viscous stresses—preserving energy and momentum. The subsequent medium evolution is performed with MUSIC, a second-order relativistic viscous hydrodynamics code following Israel–Stewart-type equations with shear-stress tensor $\pi^{\mu\nu}$ and bulk pressure $\chi_{npq}$0, and using a lattice-QCD-based equation of state. When the hydrodynamic medium cools to the switching temperature $\chi_{npq}$1 MeV, particles are emitted from the isothermal hypersurface via the Cooper–Frye prescription, and the hadronic stage is initialized. UrQMD then simulates resonance decays and hadronic re-scatterings until kinetic freeze-out, encoding late-stage viscous damping and hadronic chemistry effects that modify spectra and anisotropies (McDonald et al., 2017).

In the 2.76 TeV implementation, the collision system is Pb–Pb at $\chi_{npq}$2 TeV. Centrality selection and detailed simulation parameters follow the implementation in the companion work (McDonald et al., 2016), with pre-equilibrium classical Yang–Mills until $\chi_{npq}$3. The calculations shown include results for 0–50% centrality bins and differential selections within those bins for event-shape engineering (McDonald et al., 2017).

2. Initial-state dynamics and hydrodynamic initialization

Within the IP-Glasma stage, the fundamental object passed downstream is the event-by-event stress tensor generated by classical gauge fields. The initial spatial anisotropies are encoded in the eccentricities $\chi_{npq}$4 and participant-plane angles $\chi_{npq}$5, defined from the transverse energy density $\chi_{npq}$6 as

$\chi_{npq}$7

The $\chi_{npq}$8 weight reflects the relevance of different harmonics to the hydrodynamic response. This definition is central because the later decomposition of final-state flow into linear and non-linear pieces is formulated relative to these initial eccentricities (McDonald et al., 2017).

Landau matching at $\chi_{npq}$9 converts the Yang–Mills tensor into hydrodynamic variables. In the framework summary, this matching is described as preserving energy and momentum while defining the initial conditions for viscous hydrodynamics. In the 2.76 TeV study, the hydrodynamic evolution uses a constant shear viscosity to entropy ratio $\eta/s$0 together with a temperature-dependent bulk viscosity $\eta/s$1 based on parametrizations motivated by Noronha-Hostler–Noronha–Greiner and Karsch–Kharzeev–Tuchin. The relaxation dynamics for $\eta/s$2 and $\eta/s$3 and the remaining second-order transport coefficients are those standardly used in MUSIC for phenomenology, although their explicit numerical values are not specified in the proceedings summary (McDonald et al., 2017).

A notable feature of the IP-Glasma initialization is that the early Yang–Mills evolution naturally generates pre-equilibrium flow. A separate 5.02 TeV study quantified the pre-equilibrium longitudinal component and found, for 100 events at $\eta/s$4 TeV in 0–5% central collisions at $\eta/s$5 fm/c, that the energy-density-weighted RMS value satisfies $\eta/s$6, emphasizing that the Glasma stress tensor generically contains non-zero longitudinal flow already before hydrodynamics (McDonald et al., 2017). This suggests that the matching surface is not merely a scalar-energy-density handoff, but a transfer of a tensorial pre-equilibrium state.

3. Hydrodynamics, particlization, and hadronic transport

The hydrodynamic stage evolves the viscous energy-momentum tensor using the lattice-QCD-based equation of state and the specified transport coefficients. In the 2.76 TeV synthesis study, the constant $\eta/s$7 was adopted throughout, while $\eta/s$8 peaks in the crossover region and was taken from earlier work in which bulk viscosity was shown to be essential for simultaneously describing yields and flow magnitudes (McDonald et al., 2017, Ryu et al., 2015). The switching hypersurface is found with the standard surface finder used in MUSIC, and details of the numerical implementation and hypersurface construction follow earlier work (McDonald et al., 2017).

Particlization is implemented through the Cooper–Frye formula,

$\eta/s$9

with viscous-corrected distributions $\zeta/s(T)$0. The shear correction is taken in the standard 14-moment form,

$\zeta/s(T)$1

and the bulk correction is taken in a commonly used ansatz,

$\zeta/s(T)$2

These forms, standard in MUSIC-based phenomenology, control the $\zeta/s(T)$3 dependence of identified spectra and $\zeta/s(T)$4 (McDonald et al., 2017).

The UrQMD stage simulates resonance decays and hadronic re-scatterings until kinetic freeze-out, using a broad set of hadronic species and cross sections. In the framework description, this late hadronic stage is not a small correction but an explicit transport sector that modifies both spectra and anisotropies through late-stage viscous damping and hadronic chemistry effects (McDonald et al., 2017). Earlier work on the MUSIC–UrQMD coupling emphasized the importance of a switching window in which both viscous hydrodynamics and hadronic kinetic theory are applicable, and established the technical infrastructure for passing hypersurfaces and sampled particles from MUSIC into UrQMD (Ryu et al., 2012).

4. Observable basis and analysis formalism

The framework is distinguished by the breadth of observables it synthesizes. Flow vectors are written event-by-event as

$\zeta/s(T)$5

with two-particle and four-particle cumulants defined by

$\zeta/s(T)$6

and

$\zeta/s(T)$7

Event-shape engineering is implemented through the reduced flow vector

$\zeta/s(T)$8

with $\zeta/s(T)$9 used to select event shape within fixed centrality (McDonald et al., 2017).

A central formal ingredient is the decomposition of higher harmonics into linear and non-linear response pieces. In schematic form,

$\tau_0$0

The study focuses on the standard set

$\tau_0$1

with event-averaged extractions

$\tau_0$2

$\tau_0$3

The RMS linear piece is then obtained through

$\tau_0$4

This decomposition supplies a direct language for separating geometry-driven and mode-coupled contributions to $\tau_0$5, $\tau_0$6, and $\tau_0$7 (McDonald et al., 2017).

Event-plane correlations are constructed from azimuthal constraints of the form $\tau_0$8. Two-plane correlators are written as

$\tau_0$9

The analysis includes correlators such as $\tau_0$0, $\tau_0$1, and $\tau_0$2, following ATLAS measurements at 2.76 TeV. Correlations of flow magnitudes are encoded by the symmetric cumulants

$\tau_0$3

Taken together, event-shape engineering, event-plane correlators, $\tau_0$4 coefficients, and symmetric cumulants provide a tightly interlocked description of fluctuation structure and multi-harmonic response (McDonald et al., 2017).

5. Pb–Pb results at $\tau_0$5 TeV

The 2.76 TeV study reported several qualitative and semi-quantitative findings. For event-shape engineering, using the reduced flow vector $\tau_0$6 within fixed centrality bins and the ATLAS selection $\tau_0$7, the framework reproduces the wide spread of $\tau_0$8 values and the robust anti-correlation between $\tau_0$9 and $T^{\mu\nu}(\tau_0)$0 in more peripheral collisions. At the same time, integrated $T^{\mu\nu}(\tau_0)$1 are somewhat overestimated in this kinematic window. The proceedings trace this to transport coefficients tuned to ALICE data with the wider range $T^{\mu\nu}(\tau_0)$2–$T^{\mu\nu}(\tau_0)$3 GeV, together with the model’s known $T^{\mu\nu}(\tau_0)$4 dependence of $T^{\mu\nu}(\tau_0)$5: slight underestimation below $T^{\mu\nu}(\tau_0)$6 GeV and overestimation above, so the narrower ATLAS cut biases the integrated $T^{\mu\nu}(\tau_0)$7 upward compared to the ALICE-tuned calibration (McDonald et al., 2017).

For event-plane correlations, the correlators involving $T^{\mu\nu}(\tau_0)$8 and $T^{\mu\nu}(\tau_0)$9 grow from central to peripheral collisions, while those involving $e$0 and $e$1 decrease. The transition point where the $e$2–$e$3 correlation overtakes $e$4–$e$5 occurs at roughly 20% centrality. This pattern is presented as consistent with the non-linear mode-coupling picture in which $e$6 draws strength from both $e$7 and $e$8 (McDonald et al., 2017).

The decomposition of $e$9 into linear and non-linear parts highlights the dominance hierarchy of mode coupling. Two quantitative results are emphasized. First, $u^\mu$0 is dominated by non-linear mode coupling: $u^\mu$1 is approximately two orders of magnitude smaller than the non-linear contributions $u^\mu$2 and $u^\mu$3. Second, with increasing centrality percentile, $u^\mu$4 rises and overtakes $u^\mu$5 around 20% centrality, tracking the growth of the global elliptic geometry and the diminishing relative role of fluctuations associated with triangularity (McDonald et al., 2017).

The extraction of $u^\mu$6 is validated event-by-event. For example, the scatter of $u^\mu$7 versus $u^\mu$8 shows that the slope of the best-fit line matches $u^\mu$9, while the vertical spread reflects the event-wise fluctuations of $\pi^{\mu\nu}$0. The proceedings note that the RMS nature of $\pi^{\mu\nu}$1 implies a slight overestimate of the intercept relative to the event-wise best-fit line (McDonald et al., 2017).

The physical interpretation of the observable set is explicit. The joint analysis of event-shape engineering, event-plane correlations, $\pi^{\mu\nu}$2 decompositions, and $\pi^{\mu\nu}$3 trends tightens constraints on both transport and initial-state parameters. The event-shape-engineering “boomerang” and the spread of $\pi^{\mu\nu}$4 within $\pi^{\mu\nu}$5 bins are described as highly sensitive to shear viscosity. Bulk viscosity primarily influences the $\pi^{\mu\nu}$6 slopes of spectra and the integrated $\pi^{\mu\nu}$7 via radial-flow competition, and can subtly shift the relative strengths $\pi^{\mu\nu}$8 and $\pi^{\mu\nu}$9 by altering $\chi_{npq}$00 and $\chi_{npq}$01 magnitudes. Taken together, these observables prefer a relatively small $\chi_{npq}$02 and a sizeable $\chi_{npq}$03 in the crossover region, corroborating earlier extractions from spectra and integrated $\chi_{npq}$04 while demonstrating consistency across more differential flow observables (McDonald et al., 2017).

6. Methodological scope, limitations, and later developments

The framework’s precision aspirations are accompanied by explicit caveats. Multi-particle and event-plane correlators demand large event statistics per centrality bin to stabilize RMS $\chi_{npq}$05 and $\chi_{npq}$06 extractions, and event-shape engineering further subdivides events by $\chi_{npq}$07, requiring additional statistics within each centrality bin. The choice of $\chi_{npq}$08 for shear and bulk corrections affects $\chi_{npq}$09 slopes and integrated $\chi_{npq}$10, particularly under restricted $\chi_{npq}$11 cuts. Sensitivity to $\chi_{npq}$12, $\chi_{npq}$13, and hypersurface construction can introduce variations, especially in spectra and low-$\chi_{npq}$14 $\chi_{npq}$15, although the qualitative trends in mode coupling and event-plane correlations are reported as robust. The proceedings do not quote explicit $\chi_{npq}$16 numerical values or detailed parameter scans of $\chi_{npq}$17 or $\chi_{npq}$18; a fuller Bayesian calibration against the expanded observable set is identified as a natural next step (McDonald et al., 2017).

Subsequent work has broadened the domain of the same hybrid logic. A later systematics study used one parameter set across p+p, p+Pb, Xe+Xe, Au+Au, and Pb+Pb, with $\chi_{npq}$19 fm/c, $\chi_{npq}$20, temperature-dependent $\chi_{npq}$21, $\chi_{npq}$22 MeV, and three subnucleon hot spots per nucleon, and concluded that the framework provides a unified description of bulk properties and multi-particle correlations from large to small systems, with caveats in pp (Schenke et al., 2020). A broader survey across Pb+Pb, Xe+Xe, O+O, Au+Au, U+U, Ru+Ru, Zr+Zr, p+Au, d+Au, $\chi_{npq}$23He+Au, p+Pb, and p+p likewise reported that many observables in the smaller systems are well described although they test the limits of the model (Schenke et al., 2020).

Several extensions target specific physics sectors. A 3+1D generalization coupled JIMWLK-evolved IP-Glasma initial conditions to MUSIC+UrQMD in order to describe longitudinal observables such as $\chi_{npq}$24, factorization ratios $\chi_{npq}$25, and $\chi_{npq}$26 in Pb–Pb at 2.76 TeV (McDonald et al., 2020). An O–O study at 7 TeV used IP-Glasma + MUSIC + iSS + UrQMD to compare Woods–Saxon and $\chi_{npq}$27-clustered nuclear geometries and found robust qualitative differences in $\chi_{npq}$28, $\chi_{npq}$29, $\chi_{npq}$30, and $\chi_{npq}$31 fluctuations versus multiplicity (Prasad et al., 2024). A high-statistics 5.02 and 5.36 TeV analysis extended the nonlinear-response program to charged-hadron $\chi_{npq}$32 up to $\chi_{npq}$33, extracting high-order mode couplings and showing that off-diagonal elements in the response matrices matter for selected coefficients at $\chi_{npq}$34 (Horecny et al., 3 Jun 2025). Heavy-flavor coupling has also been implemented: in Pb+Pb at 5.02 TeV, IP-Glasma+MUSIC+UrQMD was combined with PYTHIA, MARTINI, and heavy-quark hadronization, and the study concluded that even though there is significant momentum broadening in the earliest stage, D-meson $\chi_{npq}$35 and $\chi_{npq}$36 are only weakly sensitive to pre-equilibrium interactions (Singh et al., 23 Sep 2025).

The framework has also become a tool for nuclear-structure inference. In Xe–Xe at 5.44 TeV, comprehensive comparisons with IP-Glasma+MUSIC+UrQMD calculations using different nuclear parameters indicated that the choice $\chi_{npq}$37 and $\chi_{npq}$38 provides a better description of the presented flow measurements (Collaboration, 2024). Conversely, precision small-system data can expose limitations: ALICE measurements of $\chi_{npq}$39 in pp, p–Pb, and Pb–Pb showed that current IP-Glasma + MUSIC + UrQMD implementations fail to reproduce the positive sign and monotonic increase in the smallest systems, thereby imposing strong constraints on the existing theoretical models (Collaboration, 13 Mar 2026).

In this broader perspective, the IP-Glasma+MUSIC+UrQMD framework is best understood not as a single fixed parameterization, but as a modular event-by-event strategy in which CGC-based initial conditions, second-order viscous hydrodynamics, and hadronic transport are combined to test the consistency of QGP transport, initial-state fluctuation spectra, and late-stage hadronic dynamics across increasingly differential observables and increasingly diverse collision systems (McDonald et al., 2017, Schenke et al., 2020).