---
title: 'Angantyr: Heavy-Ion Extension of PYTHIA8'
url: https://www.emergentmind.com/topics/angantyr
type: topic
---

# Angantyr: Heavy-Ion Extension of PYTHIA8

Searching arXiv for Angantyr papers and related applications to ground the article in the literature.
arxiv_search(query="Angantyr PYTHIA heavy-ion model", max_results=10, sort_by="relevance")
arxiv_search(query="Angantyr model heavy-ion PYTHIA8 applications", max_results=10, sort_by="relevance")
Angantyr is the heavy-ion extension of PYTHIA8: a model for generating complete exclusive hadronic final states in proton–nucleus and nucleus–nucleus collisions by extrapolating the pp event-generation framework of PYTHIA8 to pA and AA systems with as few new assumptions and tunings as possible [1806.10820]. It is built around a Glauber-model description of nuclear geometry, wounded-nucleon and Fritiof-inspired ideas, and explicit treatment of multi-parton interactions and diffractive excitation in each nucleon–nucleon sub-collision [1806.10820]. Across later phenomenology, Angantyr is repeatedly used as a baseline model without a hot thermalized medium, hydrodynamic flow, or jet quenching, so that deviations between data and simulation can be interpreted as signatures of collectivity, QGP-medium effects, critical fluctuations, or other dynamics absent from the model [2309.12995].

## 1. Origins and conceptual role

Angantyr was introduced to bridge a large part of the gap between high-energy hadron phenomenology and heavy-ion phenomenology by constructing heavy-ion events as a direct extrapolation of high-energy pp collisions as described by PYTHIA [1806.10820]. The core motivation is twofold: to provide a realistic non-collective baseline for heavy-ion observables, and to test how far a pp-based description can go before genuinely collective physics must be added [1806.10820].

In the model’s own framing, heavy-ion collisions are not generated through an explicit QGP stage or hydrodynamic evolution. Instead, Angantyr generalizes the usual pp framework to pA and AA collisions by combining many nucleon–nucleon-like sub-collisions into one nuclear event [2403.01833]. This places it in a distinct methodological category relative to hybrid approaches such as SMASH-vHLLE or transport-based models such as AMPT: Angantyr is designed to capture geometry, fluctuating nucleon substructure, multiple scattering, parton showers, and Lund-string hadronization, while omitting a de-confined medium and collective expansion [2509.05613].

The baseline status of Angantyr is central to its later use. The model is described as suitable for understanding fluctuations expected from geometry, hadronic production, and acceptance effects alone, and as a reference model without deconfined QGP formation [2403.01833]. This suggests that Angantyr is best interpreted not as a complete theory of heavy-ion dynamics, but as a controlled null hypothesis for observables sensitive to collectivity, medium response, or criticality.

## 2. Event construction and nuclear geometry

Angantyr constructs a pA or AA event as a superposition of multiple nucleon–nucleon sub-collisions determined by Glauber geometry and event-by-event nucleon fluctuations [1806.10820]. In practical descriptions across the literature, nucleons are distributed inside nuclei according to Woods–Saxon profiles, the collision geometry is handled with a Glauber-model picture, and the number of wounded nucleons is determined with Gribov-corrected Glauber fluctuations [2403.01833].

At the formal level, the nuclear S-matrix is written as
$$
S^{AB}(b) = \prod_{\mu=1}^A \prod_{\nu=1}^B S^{(N_\mu,N_\nu)}(b_{\mu\nu}),
$$
with nucleon positions generated event by event and propagated on straight-line trajectories through the opposing nucleus [1806.10820]. Angantyr is explicitly inspired by the wounded-nucleon model and the old Fritiof model. In the wounded-nucleon picture, soft production scales with wounded projectile and target nucleons through
$$
\frac{dN_{ch}}{d\eta} = w_p F(\eta) + w_t F(-\eta),
$$
while the Fritiof-inspired component supplies a mechanism for longitudinal excitation and string-like particle production [1806.10820].

A distinctive feature of Angantyr is its treatment of projectile and target fluctuations. Diffractive excitation is interpreted in the Good–Walker picture as a consequence of fluctuations in the nucleon’s partonic substructure, with
$$
\frac{d\sigma_D}{d^2b} = \langle T^2\rangle - \langle T\rangle^2,
$$
so diffraction is identified with the variance of the interaction amplitude [1806.10820]. The model includes fluctuations in both projectile and target nucleons, which are represented through fluctuating effective radii drawn from a Gamma distribution and kept frozen during passage through the opposite nucleus [1806.10820]. A plausible implication is that Angantyr’s fluctuation structure is not reducible to purely geometric participant counting; it also embeds event-by-event fluctuations in the underlying nucleon state.

## 3. Sub-collision classes and PYTHIA machinery

Each nucleon–nucleon encounter in Angantyr is classified as elastic, diffractive, or absorptive, and the full heavy-ion final state is obtained by stacking the corresponding sub-events [1806.10820]. Primary absorptive collisions are generated as ordinary PYTHIA minimum-bias events, while additional non-diffractive interactions of already wounded nucleons are treated as secondary non-diffractive collisions modeled through a diffractive-like prescription [2303.11747]. This primary/secondary distinction is one of the model’s essential algorithmic devices for extrapolating pp-like event generation to nuclear collisions [2303.11747].

Within each sub-collision, Angantyr uses standard PYTHIA machinery: hard and soft interactions, initial- and final-state radiation, multi-parton interactions, parton showers, fragmentation, and decays [2108.08626]. Final hadrons are produced via Lund string fragmentation, and the resulting sub-events are merged into one nuclear event [2403.01833]. The standard PYTHIA infrared regularization of partonic scatterings,
$$
\frac{d\sigma_{2\rightarrow 2}}{dp_\perp^2} \propto \frac{\alpha_s^2(p_\perp^2)}{p_\perp^4} \rightarrow \frac{\alpha_s^2(p_\perp^2+p_{\perp 0}^2)}{(p_\perp^2+p_{\perp 0}^2)^2},
$$
is retained as part of the MPI framework that Angantyr inherits from PYTHIA8 [1806.10820].

The model’s treatment of diffractive excitation is not peripheral bookkeeping. It is tied to fluctuations in nucleon substructure and to the generation of secondary wounded-nucleon contributions [1806.10820]. In later work, Angantyr’s heavy-ion event construction is summarized as a stack of multiple PYTHIA-like nucleon–nucleon sub-collisions determined by a Glauber model with Good–Walker fluctuations, where one interaction is treated as primary and further interactions are handled as secondary non-diffractive events [2303.11747]. This suggests that Angantyr’s predictive power for multiplicity and rapidity systematics depends strongly on how these secondary excitations are modeled rather than on hydrodynamic response.

## 4. MPI, color reconnection, and locality

Two essential ingredients in Angantyr phenomenology are multi-parton interactions and color reconnection. MPI are repeatedly identified as central to multiplicity distributions and underlying-event activity, while CR reduces string length and can generate flow-like patterns even in a model without hydrodynamics [2108.08626]. In Pb–Pb studies, multiplicity distributions and the trend of $\langle p_T\rangle$ with multiplicity are described best when MPI and CR are both active, whereas turning MPI off strongly suppresses particle production; if MPI is off, turning CR on or off makes no difference [2108.08626].

This MPI–CR dependence also appears in light-nuclei studies. In Au+Au collisions, the light-nuclei yield ratio $N_tN_p/N_d^2$ is affected by CR only when MPI is present; if MPI is turned off, CR has no visible effect on the ratio [2309.12995]. A closely related analysis including three-nucleon correlations reaches the same conclusion and states it explicitly: color reconnection depends entirely on MPI; if MPI is turned off, CR has no effect [2504.17394]. This establishes a model-internal constraint: CR is not an autonomous mechanism in Angantyr, but a reorganization of the color topology generated by multiple partonic scatterings.

Originally, Angantyr effectively restricted CR to partons within individual nucleon–nucleon sub-collisions. A later extension introduced a spatially constrained QCD-based color reconnection model, allowing reconnections across different sub-collisions only if dipoles are sufficiently close in transverse space [2303.11747]. The new parameter,
$$
\texttt{ColourReconnection:dipoleMaxDist},
$$
was tuned to
$$
\delta b_c = 0.5\ \mathrm{fm},
$$
so that CR becomes local in impact parameter rather than either fully intra-subcollision or unrealistically all-to-all [2303.11747]. The same study reports a pp retune with
$\texttt{pT0Ref}=2.37$,
$\texttt{m0}=1.05$,
$\texttt{junctionCorrection}=1.37$,
and $\texttt{dipoleMaxDist}=0.5$ fm, followed by Angantyr-specific adjustments for pA and AA via $\texttt{HISigmaDiffractive:PomFluxEpsilon}=-0.04$ [2303.11747].

A common misconception is that color reconnection in Angantyr is intrinsically a heavy-ion collective mechanism. The literature instead treats it as a hadronization-level QCD effect that can mimic some collective-looking trends, such as baryon enhancement or a bump in $p/\pi$ around intermediate $p_T$, without implying a thermalized medium [2108.08626]. Later p–Pb and d–Au studies reinforce this point by showing that CR can generate flow-like peaks or baryon enhancement while leaving other medium-sensitive observables, such as high-$p_T$ suppression, absent [2407.07724].

## 5. Phenomenological domains and typical observables

Angantyr has been used in a wide range of phenomenological contexts, typically as a non-collective reference model.

| Domain | Representative observable(s) | Role of Angantyr |
|---|---|---|
| Heavy-ion fluctuations | $\omega_{ch}$, net-proton cumulants, $N_tN_p/N_d^2$ | Hadronic or non-collective baseline |
| Small-system collectivity tests | $v_n\{2\}$, $v_n\{4\}$, $R_{\mathrm{AA}}$, particle ratios | Null hypothesis for no collective effects |
| Nuclear-structure studies | $dN_{ch}/d\eta$, $\langle p_T\rangle$, $v_2$ | Propagation of initial geometry without hydrodynamics |
| Air-shower modeling | longitudinal profiles, muon distributions | High-energy hadronic interaction model in CORSIKA 8 |

In fluctuation studies, Angantyr is used to quantify what can be explained by geometry, hadronic production, and acceptance effects alone. For charged-particle multiplicity fluctuations in Au–Au and Pb–Pb, the scaled variance
$$
\omega_{ch} = \frac{\langle N_{ch}^2\rangle - \langle N_{ch}\rangle^2}{\langle N_{ch}\rangle}
$$
is found to be systematically below the expectations of a simple participant superposition model, with weak beam-energy dependence and only mild centrality dependence [2403.01833]. In net-proton studies at top RHIC and LHC energies, Angantyr provides a baseline for cumulants and their ratios in a framework where a thermalized medium is not assumed, while separate afterburner studies show that radial flow and weak decays can substantially modify these cumulants relative to the pure Angantyr baseline [1909.12113].

In light-nuclei production, Angantyr has been used to study neutron density fluctuation, neutron–proton correlation, and three-nucleon correlation effects through the yield ratio
$$
\frac{N_tN_p}{N_d^2}.
$$
The ratio is not a direct fluctuation measure unless the neutron–proton correlation term is negligible; in the model, that correlation is generally not negligible [2309.12995]. Angantyr predicts that the ratio is nearly unchanged by rapidity coverage and centrality, rises only weakly with beam energy, and cannot reproduce the non-monotonic STAR energy dependence with a peak around $\sqrt{s_{NN}}\sim 20$–$30$ GeV [2309.12995].

In small-system and light-ion studies, Angantyr frequently functions as a control model for the absence of collective expansion. In O–O collisions at $\sqrt{s_{NN}}=5.36$ TeV, Angantyr serves as a baseline for no collective effects in comparisons to SMASH and a hydrodynamic hybrid model; apparent nonzero $v_2\{2\}$ and $v_3\{2\}$ trends are interpreted as dominated by non-flow, while 4-particle cumulants remain mostly close to zero [2509.05613]. In O–O, p–O, and Ne–Ne structural studies, Angantyr has also been used to propagate initial-state $\alpha$-cluster or elongated nuclear configurations into final-state multiplicity, $\langle p_T\rangle$, and elliptic-flow-like observables without invoking hydrodynamics [2503.20339].

Beyond collider heavy-ion phenomenology, Pythia 8/Angantyr has been embedded in CORSIKA 8 as the high-energy hadronic interaction model for air-shower simulations. A combined tuning to LEP, LHC, and NA22 data adjusted the parameters `MultipartonInteractions:pT0Ref`, `StringZ:aLund`, and `StringPT:sigma`, with the study concluding that collider and fixed-target data prefer different parameter values, especially for `pT0Ref` [2508.11458]. This demonstrates that Angantyr has also become relevant outside heavy-ion phenomenology, in contexts where hadron–nucleus-like multiple scattering is essential.

## 6. Empirical performance, limitations, and extensions

Angantyr’s empirical performance is strongest for inclusive baselines and weakest for observables dominated by dense-medium dynamics. The original model paper states that it gives a good description of general final-state properties such as multiplicity and transverse-momentum distributions in pA and AA collisions, and that many observables are reproduced at the level of roughly 5–10% without tuning to heavy-ion data [1806.10820]. Later studies confirm that Angantyr often reproduces minimum-bias multiplicity and pseudorapidity distributions reasonably well, including in d–Au and p–Pb, while failing more visibly in the most central events or in observables that are known to be medium sensitive [2503.23019].

Several recurring limitations are explicit in the literature. Angantyr does not include hydrodynamic flow, a thermalized QGP medium, jet quenching, critical fluctuations near the QCD critical point, or, in the baseline setups emphasized by multiple comparison papers, collective behavior such as string shoving and ropes [2509.05613]. These omissions have direct phenomenological consequences. In Pb–Pb angular correlations, PYTHIA 8.235/Angantyr does not reproduce the ALICE near-side jet-peak widths well, especially in $\Delta\eta$, and performs worse than AMPT [1904.00304]. In identified-particle balance functions in Pb–Pb, the Monash 2013 tune describes peripheral collisions reasonably well but does not quantitatively reproduce central data, leading the authors to suggest that a dedicated heavy-ion tuning of Angantyr may be needed [2604.19585].

The model also has documented fluctuation mismatches. In p+Pb collisions at zero impact parameter, a Bayesian reconstruction of multiplicity fluctuations and rapidity correlations from ATLAS data finds that Angantyr overpredicts the relative variances and covariance of central multiplicity observables, which is interpreted as overly large proton-size fluctuations in the model [2208.12175]. In jet-background studies, the width of random-cone fluctuations in Angantyr exceeds gamma-based analytic expectations by approximately 13% unless leading jets are removed, indicating that its background is not equivalent to a random soft bath and is sensitive to the detailed single-particle $p_T$ spectrum [2005.02320].

A major extension path is to couple Angantyr to later-stage dynamics rather than to alter its QGP-free core. The Angantyr+UrQMD framework adds hadronic rescattering after the PYTHIA event generation, enabled by a hadron vertex model that supplies space-time production points to UrQMD [2002.10236]. In Pb–Pb at 2.76 TeV, this hadronic cascade leads to substantial suppression of mid- to high-$p_\perp$ yields, away-side jet suppression, and modifications of $v_2\{2\}$ even without any partonic energy loss [2002.10236]. This does not convert Angantyr into a hydrodynamic model, but it demonstrates that the non-QGP baseline can be systematically extended to include a realistic hadronic phase.

## 7. Interpretive significance

The principal significance of Angantyr lies in its role as a falsifiable baseline. Because it extrapolates PYTHIA’s pp dynamics to nuclear collisions while omitting QGP-medium formation, agreement with Angantyr indicates that a given observable may be explainable through geometry, fluctuating nucleon substructure, multiple scattering, Lund-string hadronization, MPI, CR, resonance decays, and acceptance effects alone [2403.01833]. Disagreement, by contrast, can isolate the sector of heavy-ion phenomenology that requires collective expansion, rescattering, critical dynamics, or medium-induced energy loss.

This role is especially clear in studies that contrast Angantyr with models containing explicit collectivity. In O–O collisions, the hybrid model produces negative and sizable 4-particle cumulants, larger $v_2$ and $v_3$ in central collisions, and strong intermediate-$p_T$ enhancement in $R_{\mathrm{AA}}$, whereas these signatures are absent in Angantyr [2509.05613]. In Au+Au light-nuclei observables, Angantyr underestimates STAR data and cannot reproduce the non-monotonic peak near 20–30 GeV, reinforcing the interpretation that the observed structure may be connected to physics beyond standard hadronic transport and string-based modeling [2504.17394].

At the same time, Angantyr also shows that some collective-looking signals do not uniquely imply a QGP. Spatially constrained CR can alter multiplicity systematics and baryon production across pp, pA, and AA [2303.11747]. In central d–Au, SC CR can mimic baryon enhancement without producing high-$p_T$ suppression [2503.23019]. In Pb–Pb, MPI+CR can reproduce a rising $\langle p_T\rangle$ with multiplicity and baryon/meson-ratio bumps usually associated with radial flow [2108.08626]. The appropriate conclusion is therefore not that Angantyr invalidates medium interpretations, but that it sharpens them by quantifying which effects survive in a QGP-free construction.

In this sense, Angantyr occupies a specific and durable place in contemporary nuclear-collision phenomenology: it is a PYTHIA-based, Glauber-driven, fluctuation-aware generator of complete hadronic final states for pA and AA systems, and its scientific utility derives precisely from what it excludes as much as from what it includes [1806.10820].

Source: https://www.emergentmind.com/topics/angantyr