Mapping jet substructure in heavy-ion collisions with track functions
Pith reviewed 2026-06-26 07:39 UTC · model grok-4.3
The pith
Track functions show that medium-induced energy loss modifies jet fragmentation differently across quenching models.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In both the JEWEL and HYBRID Monte-Carlo frameworks, the higher moments and cumulants of track-function distributions exhibit sizable deviations from vacuum baselines, demonstrating that medium-induced energy loss imprints itself on the in-medium jet fragmentation pattern; the quantitative magnitude of these modifications differs significantly between the two models, identifying track functions as observables capable of discriminating between competing microscopic pictures of jet-medium interactions. The in-medium renormalization-group flows for the leading moments remain in qualitative agreement with the in-vacuum evolution.
What carries the argument
Track functions, non-perturbative objects that encode the flow of energy from an initiating parton to all charged hadrons and obey non-linear renormalization-group evolution.
If this is right
- Medium-induced energy loss affects the full jet fragmentation pattern captured by track functions.
- Track functions can distinguish between different microscopic models of jet quenching.
- In-medium renormalization-group flows for leading moments can be compared directly with vacuum evolution.
- Extraction of track functions from experimental heavy-ion data is feasible though challenging.
Where Pith is reading between the lines
- Track functions could supply new experimental handles on energy-loss mechanisms that are not accessible with standard jet substructure observables.
- The approach offers a route to confront basic quantum-field-theory predictions for partonic branching inside a strongly interacting medium.
- Similar non-linear evolution objects might be constructed for other classes of hadrons or for different collision systems.
Load-bearing premise
The JEWEL and HYBRID Monte Carlo implementations faithfully capture distinct microscopic mechanisms of jet quenching.
What would settle it
A measurement of track-function moments in heavy-ion data that shows either no deviation from vacuum values or identical quantitative modifications in both JEWEL and HYBRID frameworks.
read the original abstract
We investigate the dynamics of high-energy QCD cascades in heavy-ion collisions, focusing on modifications to jets' substructure induced by the presence of a quark-gluon plasma (QGP). To this end, we study the properties and scale evolution of track functions, non-perturbative objects encoding the flow of energy from an initiating parton to all charged hadrons. In contrast to the more standard fragmentation functions, these objects have a non-linear renormalization group (RG) evolution, being sensitive to the entire jet fragmentation process. Using the JEWEL and HYBRID Monte-Carlo models of jet quenching, we find that in both frameworks the higher moments and cumulants of the track functions' distributions exhibit sizable deviations from their vacuum baselines, demonstrating that medium-induced energy loss imprints itself in the in-medium jet fragmentation pattern. We find that the quantitative magnitude of these modifications differs significantly between the two models. This identifies track functions as a class of observables capable of discriminating between competing microscopic pictures of jet-medium interactions. We further examined the (in-medium) RG flows for the leading moments, which remain in qualitative agreement with the in-vacuum evolution. This provides an avenue to directly test the RG flows inside a QGP, linking heavy-ion jet measurements to basic QFT understanding of partonic branching. Finally, we comment on the feasibility and challenges associated with extracting track functions from experimental data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that track functions—non-perturbative objects encoding energy flow from an initiating parton to charged hadrons—exhibit sizable deviations in their higher moments and cumulants from vacuum baselines in both the JEWEL and HYBRID jet-quenching Monte Carlo models. These modifications differ quantitatively between the two frameworks, positioning track functions as observables that can discriminate between competing microscopic pictures of jet-medium interactions. The work also reports that the in-medium renormalization-group flows for the leading moments remain in qualitative agreement with vacuum evolution and comments on experimental extraction challenges.
Significance. If substantiated, the results would establish track functions as a new class of jet-substructure observables sensitive to medium-induced modifications of the full fragmentation process, with potential to link heavy-ion data to fundamental QFT aspects of partonic branching via RG evolution. The use of two independent event generators provides a basic cross-check, though the discrimination power hinges on whether inter-model differences faithfully reflect distinct physics rather than shared modeling approximations.
major comments (2)
- [Abstract] Abstract: the central discrimination claim—that quantitatively different magnitudes of deviations between JEWEL and HYBRID identify track functions as model-discriminating observables—rests on the assumption that observed spreads arise purely from distinct microscopic energy-loss mechanisms. No quantitative breakdown of how the two codes differ in their in-medium track-function evolution, nor cross-validation against a third framework, is provided to support this interpretation over possible shared biases (e.g., underlying parton shower or recoil handling).
- [Abstract] Abstract: the reported 'sizable deviations' and 'significantly different' magnitudes lack any mention of the extraction procedure for track functions from the Monte Carlo events, the event statistics employed, or uncertainty quantification. Without these elements the support for the central claim that medium-induced energy loss imprints on the fragmentation pattern cannot be fully assessed.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each major comment below and indicate planned revisions where appropriate.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central discrimination claim—that quantitatively different magnitudes of deviations between JEWEL and HYBRID identify track functions as model-discriminating observables—rests on the assumption that observed spreads arise purely from distinct microscopic energy-loss mechanisms. No quantitative breakdown of how the two codes differ in their in-medium track-function evolution, nor cross-validation against a third framework, is provided to support this interpretation over possible shared biases (e.g., underlying parton shower or recoil handling).
Authors: We agree that additional discussion of model differences would strengthen the interpretation. We will add a paragraph in the revised manuscript that qualitatively compares the in-medium implementations in JEWEL and HYBRID, focusing on differences in energy-loss modeling and recoil treatment that underlie the distinct track-function modifications. This supports that the quantitative differences reflect distinct physics rather than shared biases. A cross-validation with a third framework lies outside the present scope. revision: partial
-
Referee: [Abstract] Abstract: the reported 'sizable deviations' and 'significantly different' magnitudes lack any mention of the extraction procedure for track functions from the Monte Carlo events, the event statistics employed, or uncertainty quantification. Without these elements the support for the central claim that medium-induced energy loss imprints on the fragmentation pattern cannot be fully assessed.
Authors: Details of the track-function extraction from Monte Carlo events, the event statistics used, and the uncertainty estimation are provided in Sections 2 and 3 of the manuscript. To improve the abstract's self-contained support for the claims, we will revise the abstract to include a concise reference to the Monte Carlo extraction procedure and statistical uncertainties. revision: yes
- Cross-validation of the discrimination claim against a third independent jet-quenching framework
Circularity Check
No circularity: results from direct Monte Carlo sampling in independent generators
full rationale
The paper computes track-function moments and cumulants via explicit sampling inside JEWEL and HYBRID event generators, then compares the in-medium distributions to their vacuum baselines. These quantities are obtained by running the codes on the same parton-level events with and without medium; no parameter is fitted to a subset of the data and then re-labeled as a prediction, nor is any observable defined in terms of itself. The reported inter-model differences therefore constitute an external numerical test rather than a self-referential reduction. The RG-flow remarks are presented as qualitative observations, not as load-bearing derivations. Consequently the central claim rests on independent simulation output and receives the default non-circularity score.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The Monte Carlo event generators JEWEL and HYBRID accurately represent the physics of jet quenching in the QGP.
Reference graph
Works this paper leans on
-
[1]
W. Busza, K. Rajagopal, and W. van der Schee,Heavy Ion Collisions: The Big Picture, and the Big Questions,Ann. Rev. Nucl. Part. Sci.68(2018) 339–376, [arXiv:1802.04801]
Pith/arXiv arXiv 2018
-
[2]
Y. Mehtar-Tani, J. G. Milhano, and K. Tywoniuk,Jet physics in heavy-ion collisions,Int. J. Mod. Phys. A28(2013) 1340013, [arXiv:1302.2579]
Pith/arXiv arXiv 2013
-
[3]
L. Apolinário, Y.-J. Lee, and M. Winn,Heavy quarks and jets as probes of the QGP,Prog. Part. Nucl. Phys.127(2022) 103990, [arXiv:2203.16352]
arXiv 2022
-
[4]
Mehtar-Tani,The Physics of Jet Quenching in Perturbative QCD,arXiv:2509.26394
Y. Mehtar-Tani,The Physics of Jet Quenching in Perturbative QCD,arXiv:2509.26394
-
[5]
S. Cao and X.-N. Wang,Jet quenching and medium response in high-energy heavy-ion collisions: a review,Rept. Prog. Phys.84(2021), no. 2 024301, [arXiv:2002.04028]
arXiv 2021
-
[6]
S. Cao, A. Majumder, R. Modarresi-Yazdi, I. Soudi, and Y. Tachibana,Jet quenching: From theory to simulation,Int. J. Mod. Phys. E33(2024), no. 08 2430002, [arXiv:2401.10026]
arXiv 2024
-
[7]
P. Caucal, E. Iancu, A. H. Mueller, and G. Soyez,Vacuum-like jet fragmentation in a dense QCD medium,Phys. Rev. Lett.120(2018) 232001, [arXiv:1801.09703]. – 15 –
Pith/arXiv arXiv 2018
-
[8]
L. Cunqueiro, D. Pablos, A. Soto-Ontoso, M. Spousta, A. Takacs, and M. Verweij,Isolating perturbative QCD splittings in heavy-ion collisions,Phys. Rev. D110(2024), no. 1 014015, [arXiv:2311.07643]
arXiv 2024
-
[9]
Y. Mehtar-Tani and K. Tywoniuk,Radiative energy loss of neighboring subjets,Nucl. Phys. A979(2018) 165–203, [arXiv:1706.06047]
Pith/arXiv arXiv 2018
-
[10]
Y. Mehtar-Tani and K. Tywoniuk,Sudakov suppression of jets in QCD media,Phys. Rev. D 98(2018), no. 5 051501, [arXiv:1707.07361]
Pith/arXiv arXiv 2018
- [11]
- [12]
-
[13]
M. V. Kuzmin, X. Mayo López, J. Reiten, and A. V. Sadofyev,Jet quenching in anisotropic flowing matter,Phys. Rev. D109(2024), no. 1 014036, [arXiv:2309.00683]
arXiv 2024
-
[14]
K. Boguslavski, A. Kurkela, T. Lappi, F. Lindenbauer, and J. Peuron,Jet momentum broadening during initial stages in heavy-ion collisions,Phys. Lett. B850(2024) 138525, [arXiv:2303.12595]
arXiv 2024
-
[15]
M. E. Carrington, A. Czajka, and S. Mrowczynski,Jet quenching in glasma,Phys. Lett. B 834(2022) 137464, [arXiv:2112.06812]
arXiv 2022
-
[16]
S. P. Adhya, C. A. Salgado, M. Spousta, and K. Tywoniuk,Jet quenching and scaling properties of medium-evolved gluon cascade in expanding media,PoSHardProbes2020 (2021) 129, [arXiv:2009.03122]
arXiv 2021
-
[17]
J. M. Silva and A. Soto-Ontoso,Collinear spin correlations of final-state radiation in dense QCD matter,arXiv:2511.17737
- [18]
-
[19]
S. Hauksson and E. Iancu,Jet polarisation in an anisotropic medium,JHEP08(2023) 027, [arXiv:2303.03914]
arXiv 2023
-
[20]
E. Carrió and D. Pablos,Elliptic Anisotropy from Quantum Diffraction,arXiv:2603.12346
- [21]
- [22]
- [23]
- [24]
-
[25]
J. Casalderrey-Solana, Y. Mehtar-Tani, C. A. Salgado, and K. Tywoniuk,New picture of jet quenching dictated by color coherence,Phys. Lett. B725(2013) 357–360, [arXiv:1210.7765]. – 16 –
Pith/arXiv arXiv 2013
-
[26]
Y. Mehtar-Tani, C. A. Salgado, and K. Tywoniuk,The Radiation pattern of a QCD antenna in a dense medium,JHEP10(2012) 197, [arXiv:1205.5739]
Pith/arXiv arXiv 2012
-
[27]
J.-P. Blaizot, F. Dominguez, E. Iancu, and Y. Mehtar-Tani,Probabilistic picture for medium-induced jet evolution,JHEP06(2014) 075, [arXiv:1311.5823]
Pith/arXiv arXiv 2014
-
[28]
M. V. Kuzmin and J. M. Silva,QCD antenna radiative spectrum in dense media within the Improved Opacity Expansion,JHEP12(2025) 022, [arXiv:2507.00143]. [29]CMSCollaboration, A. M. Sirunyan et al.,Measurement of the Splitting Function in pp and Pb-Pb Collisions at √sNN = 5.02TeV,Phys. Rev. Lett.120(2018) 142302, [arXiv:1708.09429]. [30]ALICECollaboration, S...
arXiv 2025
-
[29]
J. Casalderrey-Solana, Z. Hulcher, G. Milhano, D. Pablos, and K. Rajagopal,Simultaneous description of hadron and jet suppression in heavy-ion collisions,Phys. Rev. C99(2019), no. 5 051901, [arXiv:1808.07386]
Pith/arXiv arXiv 2019
- [30]
-
[31]
W. Chen, S. Cao, T. Luo, L.-G. Pang, and X.-N. Wang,Medium modification ofγ-jet fragmentation functions in Pb+Pb collisions at LHC,Phys. Lett. B810(2020) 135783, [arXiv:2005.09678]. [41]JETSCAPECollaboration, Y. Tachibana et al.,Hard jet substructure in a multistage approach,Phys. Rev. C110(2024), no. 4 044907, [arXiv:2301.02485]
arXiv 2020
-
[32]
Y. L. Dokshitzer,Calculation of the structure functions for deep inelastic scattering and – 17 – e+e− annihilation by perturbation theory in quantum chromodynamics,Sov. Phys. JETP46 (1977) 641–653
1977
-
[33]
Altarelli and G
G. Altarelli and G. Parisi,Asymptotic freedom in parton language,Nucl. Phys. B126(1977) 298–318
1977
-
[34]
H.-M. Chang, M. Procura, J. Thaler, and W. J. Waalewijn,Calculating Track-Based Observables for the LHC,Phys. Rev. Lett.111(2013) 102002, [arXiv:1303.6637]
Pith/arXiv arXiv 2013
-
[35]
H.-M. Chang, M. Procura, J. Thaler, and W. J. Waalewijn,Calculating Track Thrust with Track Functions,Phys. Rev. D88(2013) 034030, [arXiv:1306.6630]
Pith/arXiv arXiv 2013
-
[36]
Y. Li, I. Moult, S. S. van Velzen, W. J. Waalewijn, and H. X. Zhu,Extending Precision Perturbative QCD with Track Functions,Phys. Rev. Lett.128(2022), no. 18 182001, [arXiv:2108.01674]
arXiv 2022
-
[37]
M. Jaarsma, Y. Li, I. Moult, W. Waalewijn, and H. X. Zhu,Renormalization group flows for track function moments,JHEP06(2022) 139, [arXiv:2201.05166]
arXiv 2022
-
[38]
H. Chen, M. Jaarsma, Y. Li, I. Moult, W. J. Waalewijn, and H. X. Zhu,Multi-collinear splitting kernels for track function evolution,JHEP07(2023) 185, [arXiv:2210.10058]
arXiv 2023
-
[39]
H. Chen, M. Jaarsma, Y. Li, I. Moult, W. J. Waalewijn, and H. X. Zhu,Collinear parton dynamics beyond Dokshitzer-Gribov-Lipatov-Altarelli-Parisi framework,Phys. Rev. D111 (2025), no. 7 076021, [arXiv:2210.10061]
arXiv 2025
-
[40]
I. Moult and H. X. Zhu,Energy Correlators: A Journey From Theory to Experiment, arXiv:2506.09119
-
[41]
M. Jaarsma, Y. Li, I. Moult, W. J. Waalewijn, and H. X. Zhu,Energy correlators on tracks: resummation and non-perturbative effects,JHEP12(2023) 087, [arXiv:2307.15739]
arXiv 2023
-
[42]
H. Chen, I. Moult, X. Zhang, and H. X. Zhu,Rethinking jets with energy correlators: Tracks, resummation, and analytic continuation,Phys. Rev. D102(2020), no. 5 054012, [arXiv:2004.11381]
arXiv 2020
- [43]
- [44]
- [45]
- [46]
-
[47]
J. Barata and R. Szafron,Leading order track functions in a hot and dense QGP,Phys. Rev. D110(2024), no. 3 L031501, [arXiv:2401.04164]
arXiv 2024
-
[48]
Mehtar-Tani,Non-linear dynamics of jet quenching,JHEP04(2025) 163, [arXiv:2411.11992]
Y. Mehtar-Tani,Non-linear dynamics of jet quenching,JHEP04(2025) 163, [arXiv:2411.11992]
arXiv 2025
-
[49]
J. Casalderrey-Solana, D. C. Gulhan, J. G. Milhano, D. Pablos, and K. Rajagopal,A Hybrid Strong/Weak Coupling Approach to Jet Quenching,JHEP10(2014) 019, [arXiv:1405.3864]. [Erratum: JHEP 09, 175 (2015)]
Pith/arXiv arXiv 2014
-
[50]
J. Casalderrey-Solana, D. C. Gulhan, J. G. Milhano, D. Pablos, and K. Rajagopal, Predictions for Boson-Jet Observables and Fragmentation Function Ratios from a Hybrid Strong/Weak Coupling Model for Jet Quenching,JHEP03(2016) 053, [arXiv:1508.00815]. – 18 –
Pith/arXiv arXiv 2016
-
[51]
J. Casalderrey-Solana, D. Gulhan, G. Milhano, D. Pablos, and K. Rajagopal,Angular Structure of Jet Quenching Within a Hybrid Strong/Weak Coupling Model,JHEP03(2017) 135, [arXiv:1609.05842]
Pith/arXiv arXiv 2017
-
[52]
Bierlich et al.,A comprehensive guide to the physics and usage of PYTHIA 8.3,SciPost Phys
C. Bierlich et al.,A comprehensive guide to the physics and usage of PYTHIA 8.3,SciPost Phys. Codeb.2022(2022) 8, [arXiv:2203.11601]
Pith/arXiv arXiv 2022
-
[53]
J. Casalderrey-Solana, J. G. Milhano, and P. Quiroga-Arias,Out of Medium Fragmentation from Long-Lived Jet Showers,Phys. Lett. B710(2012) 175–181, [arXiv:1111.0310]
Pith/arXiv arXiv 2012
-
[54]
P. M. Chesler and K. Rajagopal,Jet quenching in strongly coupled plasma,Phys. Rev. D90 (2014), no. 2 025033, [arXiv:1402.6756]
Pith/arXiv arXiv 2014
-
[55]
P. M. Chesler and K. Rajagopal,On the Evolution of Jet Energy and Opening Angle in Strongly Coupled Plasma,JHEP05(2016) 098, [arXiv:1511.07567]
Pith/arXiv arXiv 2016
-
[56]
P. M. Chesler and L. G. Yaffe,The Stress-energy tensor of a quark moving through a strongly-coupled N=4 supersymmetric Yang-Mills plasma: Comparing hydrodynamics and AdS/CFT,Phys. Rev. D78(2008) 045013, [arXiv:0712.0050]
Pith/arXiv arXiv 2008
-
[57]
Z. Hulcher, A. S. Kudinoor, D. Pablos, and K. Rajagopal,Sensitivity of Jet Observables to Molière Scattering Off Quasiparticles in Quark-Gluon Plasma,arXiv:2603.08776
-
[58]
Z. Hulcher, D. Pablos, and K. Rajagopal,Resolution Effects in the Hybrid Strong/Weak Coupling Model,JHEP03(2018) 010, [arXiv:1707.05245]
Pith/arXiv arXiv 2018
-
[59]
A. B. Migdal,Bremsstrahlung and pair production in condensed media at high energies,Phys. Rev.103(1956) 1811–1820
1956
-
[60]
L. D. Landau and I. Pomeranchuk,Electron cascade process at very high energies,Dokl. Akad. Nauk Ser. Fiz.92(1953) 735–738
1953
-
[61]
L. D. Landau and I. Pomeranchuk,Limits of applicability of the theory of bremsstrahlung electrons and pair production at high energies,Dokl. Akad. Nauk Ser. Fiz.92(1953) 535
1953
-
[62]
G. Milhano, U. A. Wiedemann, and K. C. Zapp,Sensitivity of jet substructure to jet-induced medium response,Phys. Lett. B779(2018) 409–413, [arXiv:1707.04142]
Pith/arXiv arXiv 2018
-
[63]
Zapp,A dynamical implementation of colour coherence for quenched jets in JEWEL, arXiv:2604.13932
K. Zapp,A dynamical implementation of colour coherence for quenched jets in JEWEL, arXiv:2604.13932
-
[64]
D. Pablos and S. Sanjurjo,Color coherence effects in dipole-quark scattering in the soft limit, Phys. Rev. D110(2024), no. 11 L111502, [arXiv:2406.08550]
arXiv 2024
-
[65]
R. Kunnawalkam Elayavalli and K. C. Zapp,Simulating V+jet processes in heavy ion collisions with JEWEL,Eur. Phys. J. C76(2016), no. 12 695, [arXiv:1608.03099]
Pith/arXiv arXiv 2016
-
[66]
M. Cacciari, G. P. Salam, and G. Soyez,The anti-kt jet clustering algorithm,JHEP04 (2008) 063, [arXiv:0802.1189]
Pith/arXiv arXiv 2008
-
[67]
M. Cacciari, G. P. Salam, and G. Soyez,FastJet User Manual,Eur. Phys. J. C72(2012) 1896, [arXiv:1111.6097]. [78]A TLASCollaboration, ATLAS,Measurement of jet track functions inppcollisions at√s= 13TeV with the ATLAS detector,arXiv:2502.02062. [79]CMSCollaboration, A. Belyaev et al.,Probing early parton emissions in heavy ion collisions using the Lund jet ...
Pith/arXiv arXiv 2012
-
[68]
K. Lee, I. Moult, F. Ringer, and W. J. Waalewijn,A formalism for extracting track functions from jet measurements,JHEP01(2024) 194, [arXiv:2308.00028]
arXiv 2024
-
[69]
B. Schenke, C. Gale, and S. Jeon,MARTINI: An Event generator for relativistic heavy-ion collisions,Phys. Rev. C80(2009) 054913, [arXiv:0909.2037]
Pith/arXiv arXiv 2009
-
[70]
Y. He, T. Luo, X.-N. Wang, and Y. Zhu,Linear Boltzmann Transport for Jet Propagation in the Quark-Gluon Plasma: Elastic Processes and Medium Recoil,Phys. Rev. C91(2015) 054908, [arXiv:1503.03313]. [Erratum: Phys.Rev.C 97, 019902 (2018)]
Pith/arXiv arXiv 2015
- [71]
-
[72]
Y. Mehtar-Tani, D. Pablos, and K. Tywoniuk,Cone-Size Dependence of Jet Suppression in Heavy-Ion Collisions,Phys. Rev. Lett.127(2021), no. 25 252301, [arXiv:2101.01742]
arXiv 2021
-
[73]
A. Takacs and K. Tywoniuk,Quenching effects in the cumulative jet spectrum,JHEP10 (2021) 038, [arXiv:2103.14676]
arXiv 2021
-
[74]
L. Apolinário, D. Costa, and A. Soto-Ontoso,Bottom-up approach to describe groomed jet data in heavy-ion collisions,arXiv:2601.13310
-
[75]
J. F. Grosse-Oetringhaus and U. A. Wiedemann,A Decade of Collectivity in Small Systems, arXiv:2407.07484. – 20 – Appendix The purpose of this Appendix is to complement the material provided in the main text by varying the scaleµ, the jet radiusR, and turning on and off medium response. A Impact of Jet Radius on Track Function Distribution Figure 8 shows t...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.