Opacity estimation of OO collision from CoMBolt-ITA hybrid
Pith reviewed 2026-05-21 18:40 UTC · model grok-4.3
The pith
The CoMBolt-ITA hybrid model shows OO collisions above 60 percent centrality leave the fluid-like evolution domain.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using the CoMBolt-ITA hybrid model to evolve initial conditions from TrENTo through a Boltzmann description of collective excitations and then UrQMD, the authors find that OO collisions at centralities larger than 60 percent leave the fluid-like evolution domain based on data-model comparisons and the current status of the model.
What carries the argument
The CoMBolt-ITA hybrid model, which evolves the Boltzmann distribution of massless collective excitations to handle pre-equilibration and hydrodynamized stages consistently.
If this is right
- Peripheral OO collisions behave more like low-opacity small systems with reduced collectivity.
- The transition out of fluid-like evolution occurs as the system spatial size approaches the mean free path.
- Data from OO collisions at the LHC can map the onset of hydrodynamic applicability in intermediate-size systems.
- Central collisions remain in the high-opacity regime while peripheral ones do not.
Where Pith is reading between the lines
- Similar centrality-dependent transitions may appear in other intermediate systems such as proton-nucleus collisions.
- The opacity criterion could be tested at different beam energies to check scaling with system lifetime.
Load-bearing premise
The CoMBolt-ITA evolution and its coupling to TrENTo and UrQMD remain valid without additional tuning when applied to the smaller OO system size and its specific initial-state geometry.
What would settle it
A clear mismatch between the model predictions and LHC data on flow observables in peripheral OO collisions that persists after varying only the initial-state geometry would indicate the system has not left the fluid regime.
Figures
read the original abstract
Understanding the effect of system size on the applicability of the hydrodynamic description in heavy-ion physics remains unclear. Recent measurements of OO collisions at the LHC offer a new opportunity to refine our understanding of collectivity because of their intermediate size relative to heavy-ion and small-system collisions, as well as the relatively good control over their initial state. We use the CoMBolt-ITA hybrid model to describe recent OO measurements at the LHC. The model employs TrENTo for the initial state. A combination of the pre-equilibration and hydrodynamized medium stages is modeled consistently by CoMBolt-ITA, which evolves the Boltzmann distribution of massless collective excitations. The afterburner stage is included by employing UrQMD. Using this approach, we test whether the system lies in the regime where its spatial size approaches the mean free path, corresponding to low opacity, or in the opposite limit, where its size exceeds the mean free path sufficiently to enter the fluid-like evolution regime with high opacity. We find that, in light of the data-model comparison and considering the current status of the model, OO collisions with centralities larger than $60\%$ gradually leave the domain of fluid-like evolution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies the CoMBolt-ITA hybrid model (TrENTo initial state, Boltzmann evolution for pre-equilibrium and hydro stages, UrQMD afterburner) to recent LHC OO collision data. It performs a data-model comparison to test whether the system size approaches the mean free path (low opacity) or exceeds it (high-opacity fluid-like regime), concluding that OO collisions with centralities larger than 60% gradually leave the fluid-like domain.
Significance. If the untuned hybrid evolution remains quantitatively reliable for the reduced system size and initial-state geometry of OO collisions, the result would supply a concrete, falsifiable estimate of the opacity transition in intermediate-size systems and help delineate the boundary between hydrodynamic and non-hydrodynamic regimes.
major comments (2)
- [Model description and results section] The central claim that centralities >60% leave the high-opacity regime rests on the assumption that CoMBolt-ITA transport coefficients and matching parameters calibrated on larger systems remain valid for OO without retuning. No dedicated sensitivity study or cross-check against OO-specific observables is presented to confirm that the model's built-in mean-free-path scale matches the actual OO geometry.
- [Results and discussion] The data-model comparison is described only qualitatively in the abstract and summary; quantitative measures (chi-squared, pull values, or explicit opacity metric with uncertainties) are not reported, nor is the centrality binning or the precise observable set used to decide the fluid-like boundary.
minor comments (2)
- [Introduction] Notation for the opacity parameter and the precise definition of the 'fluid-like domain' should be introduced explicitly with an equation or table rather than left implicit.
- [Figures] Figure captions should state the centrality range and the specific observable (multiplicity, v2, HBT radii) shown in each panel.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important aspects of model applicability and quantitative rigor that we address point by point below. We have revised the manuscript to incorporate additional analysis and clarifications where feasible.
read point-by-point responses
-
Referee: [Model description and results section] The central claim that centralities >60% leave the high-opacity regime rests on the assumption that CoMBolt-ITA transport coefficients and matching parameters calibrated on larger systems remain valid for OO without retuning. No dedicated sensitivity study or cross-check against OO-specific observables is presented to confirm that the model's built-in mean-free-path scale matches the actual OO geometry.
Authors: We agree that explicit validation for the smaller OO system is valuable. The CoMBolt-ITA framework uses a microscopic Boltzmann transport for the pre-equilibrium and hydrodynamic stages, with the mean free path determined by the local density and cross sections; this construction is intended to be system-size agnostic once the initial geometry is fixed by TrENTo. Nevertheless, to strengthen the manuscript we will add a dedicated sensitivity subsection in which we vary the key transport parameters (shear viscosity to entropy ratio and matching time) within the ranges previously constrained by PbPb data and demonstrate that the 60% centrality threshold for departure from the high-opacity regime remains stable. We will also include a direct comparison of the model to the measured charged-particle multiplicity in OO collisions as an additional cross-check on the initial-state and evolution scales. revision: yes
-
Referee: [Results and discussion] The data-model comparison is described only qualitatively in the abstract and summary; quantitative measures (chi-squared, pull values, or explicit opacity metric with uncertainties) are not reported, nor is the centrality binning or the precise observable set used to decide the fluid-like boundary.
Authors: We acknowledge that the current presentation relies on qualitative visual agreement. In the revised version we will expand the results section to report quantitative metrics: we will tabulate chi-squared per degree of freedom for the primary observables (v2, v3, and dNch/deta) across the centrality bins 0-10%, 10-30%, 30-60%, and 60-100%, specify the exact observable set and centrality binning employed to identify the fluid-like boundary, and introduce an explicit opacity metric (system size divided by estimated mean free path) with uncertainties propagated from the parameter variations. These additions will make the decision criterion for the 60% threshold fully transparent. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper applies the established CoMBolt-ITA hybrid model (with TrENTo initial state and UrQMD afterburner) to OO collision data at the LHC without reported retuning of transport coefficients or matching parameters for the smaller system. The central claim—that centralities >60% leave the high-opacity fluid-like regime—follows from direct data-model comparison of observables such as multiplicity, flow, and HBT radii. No quoted step reduces a prediction to a fitted input by construction, invokes a self-citation as the sole justification for a uniqueness theorem, or renames a known result; the model is treated as an external tool whose applicability to OO geometry is an explicit (and separately debatable) assumption rather than a definitional tautology.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
K. H. Ackermann et al. Elliptic flow in Au + Au collisions at (S(NN))**(1/2) = 130 GeV.Phys. Rev. Lett., 86:402–407, 2001
work page 2001
-
[2]
Anisotropy as a signature of transverse collective flow.Phys
Jean-Yves Ollitrault. Anisotropy as a signature of transverse collective flow.Phys. Rev. D, 46:229–245, 1992
work page 1992
-
[3]
Hydrodynamic De- scription of the Quark-Gluon Plasma
Ulrich Heinz and Björn Schenke. Hydrodynamic De- scription of the Quark-Gluon Plasma. 12 2024
work page 2024
-
[4]
Vardan Khachatryan et al. Observation of Long-Range Near-Side Angular Correlations in Proton-Proton Col- lisions at the LHC.JHEP, 09:091, 2010
work page 2010
-
[5]
Long-range angular correlations on the near and away side inp-Pb collisions at√sNN = 5.02TeV.Phys
Betty Abelev et al. Long-range angular correlations on the near and away side inp-Pb collisions at√sNN = 5.02TeV.Phys. Lett. B, 719:29–41, 2013
work page 2013
-
[6]
Georges Aad et al. Observation of Associated Near- Side and Away-Side Long-Range Correlations in√sNN=5.02 TeV Proton-Lead Collisions with the AT- LAS Detector.Phys. Rev. Lett., 110(18):182302, 2013
work page 2013
-
[7]
Shreyasi Acharya et al. Multiplicity and event-scale dependent flow and jet fragmentation in pp collisions at √s = 13 TeV and in p–Pb collisions at√sNN = 5.02 TeV.JHEP, 03:092, 2024
work page 2024
-
[8]
L. Adamczyk et al. Long-range pseudorapidity di- hadron correlations ind+Au collisions at√sNN = 200 GeV.Phys. Lett. B, 747:265–271, 2015
work page 2015
-
[9]
C. Aidala et al. Creation of quark–gluon plasma droplets with three distinct geometries.Nature Phys., 15(3):214–220, 2019
work page 2019
-
[10]
Y. C. Chen et al. Two-particle angular correlations in e+e−collisions to hadronic final states in two reference coordinates at Belle.JHEP, 03:171, 2023
work page 2023
-
[11]
Long-range near-side correla- tion in e+e−collisions at 183-209 GeV with ALEPH archived data.Phys
Yu-Chen Chen et al. Long-range near-side correla- tion in e+e−collisions at 183-209 GeV with ALEPH archived data.Phys. Lett. B, 856:138957, 2024
work page 2024
- [12]
-
[13]
Two-particle azimuthal cor- relations in γp interactions using pPb collisions at sNN=8.16TeV.Phys
Armen Tumasyan et al. Two-particle azimuthal cor- relations in γp interactions using pPb collisions at sNN=8.16TeV.Phys. Lett. B, 844:137905, 2023
work page 2023
-
[14]
Aram Hayrapetyan et al. Observation of En- hanced Long-Range Elliptic Anisotropies Inside High- Multiplicity Jets in pp Collisions at s=13 TeV.Phys. Rev. Lett., 133(14):142301, 2024
work page 2024
-
[15]
A Multi-phase transport model for relativistic heavy ion collisions.Phys
Zi-Wei Lin, Che Ming Ko, Bao-An Li, Bin Zhang, and Subrata Pal. A Multi-phase transport model for relativistic heavy ion collisions.Phys. Rev. C, 72:064901, 2005
work page 2005
-
[16]
The Angantyr model for Heavy-Ion Collisions in PYTHIA8.JHEP, 10:134, 2018
Christian Bierlich, Gösta Gustafson, Leif Lönnblad, and Harsh Shah. The Angantyr model for Heavy-Ion Collisions in PYTHIA8.JHEP, 10:134, 2018
work page 2018
-
[17]
Evidence of nuclear geometry-driven anisotropic flow in OO and Ne−Ne collisions at√sNN = 5.36 TeV
Ibrahim Jaser Abualrob et al. Evidence of nuclear geometry-driven anisotropic flow in OO and Ne−Ne collisions at√sNN = 5.36 TeV. 9 2025
work page 2025
-
[18]
Georges Aad et al. Measurement of the azimuthal anisotropy of charged particles in√sNN = 5.36TeV 16O+16O and 20Ne+20Ne collisions with the ATLAS detector. 9 2025
work page 2025
-
[19]
Aram Hayrapetyan et al. Observation of long-range collective flow in OO and NeNe collisions and implica- tions for nuclear structure studies. 10 2025
work page 2025
-
[20]
Exploiting Ne20 Isotopes for Precision Characterizations of Collectivity in Small Systems.Phys
Giuliano Giacalone et al. Exploiting Ne20 Isotopes for Precision Characterizations of Collectivity in Small Systems.Phys. Rev. Lett., 135(1):012302, 2025
work page 2025
- [21]
-
[22]
Marco Ruggieri, Francesco Scardina, Salvatore Plumari, and Vincenzo Greco. Thermalization, Isotropization and Elliptic Flow from Nonequilibrium Initial Conditions with a Saturation Scale.Phys. Rev. C, 89(5):054914, 2014
work page 2014
-
[23]
Retarded correlators in kinetic theory: branch cuts, poles and hydrodynamic onset transitions.Eur
Paul Romatschke. Retarded correlators in kinetic theory: branch cuts, poles and hydrodynamic onset transitions.Eur. Phys. J. C, 76(6):352, 2016
work page 2016
-
[24]
Analytic structure of nonhydrodynamic modes in kinetic theory
Aleksi Kurkela and Urs Achim Wiedemann. Analytic structure of nonhydrodynamic modes in kinetic theory. Eur. Phys. J. C, 79(9):776, 2019
work page 2019
-
[25]
Flow in AA and pA as an interplay of fluid-like and non-fluid like excitations.Eur
Aleksi Kurkela, Urs Achim Wiedemann, and Bin Wu. Flow in AA and pA as an interplay of fluid-like and non-fluid like excitations.Eur. Phys. J. C, 79(11):965, 2019
work page 2019
- [26]
-
[27]
Vincenzo Nugara, Vincenzo Greco, and Salvatore Plumari. Far-from-equilibrium attractors with Full Relativistic Boltzmann approach in 3+1D: moments of distribution function and anisotropic flowsvn.Eur. Phys. J. C, 85(3):311, 2025
work page 2025
-
[28]
Ambrus, Sören Schlichting, and Clemens Werthmann
Victor E. Ambrus, Sören Schlichting, and Clemens Werthmann. Collective dynamics in heavy and light- ion collisions – I) Kinetic Theory vs. Hydrodynamics. 11 2024
work page 2024
-
[29]
S. F. Taghavi, S. M. A. Tabatabaee Mehr, and F. Taghinavaz. CoMBolt-ITA, a Collective Model via relativistic Boltzmann equation in Isotropization Time Approximation. 4 2025
work page 2025
-
[30]
S. F. Taghavi. CoMBolt-ITA hybrid: Kinetic- theory-based event generator. https://github. com/FaridTaghavi/ComBolt-Hybrid, 2025. Public GitHub repository
work page 2025
-
[31]
J. Scott Moreland, Jonah E. Bernhard, and Steffen A. Bass. Alternative ansatz to wounded nucleon and bi- nary collision scaling in high-energy nuclear collisions. Phys. Rev. C, 92(1):011901, 2015
work page 2015
-
[32]
Jonah E. Bernhard, J. Scott Moreland, and Steffen A. Bass. Bayesian estimation of the specific shear and bulk viscosity of quark-gluon plasma.Nature Physics, 2019
work page 2019
-
[33]
D. Everett et al. Multisystem Bayesian constraints on the transport coefficients of QCD matter.Phys. Rev. C, 103(5):054904, 2021
work page 2021
-
[34]
Bayesian analysis of heavy ion collisions with the heavy ion computational framework Trajectum.Phys
Govert Nijs, Wilke van der Schee, Umut Gürsoy, and Raimond Snellings. Bayesian analysis of heavy ion collisions with the heavy ion computational framework Trajectum.Phys. Rev. C, 103(5):054909, 2021
work page 2021
-
[35]
J. E. Parkkila, A. Onnerstad, and D. J. Kim. Bayesian estimation of the specific shear and bulk viscosity of the quark-gluon plasma with additional flow harmonic observables.Phys. Rev. C, 104(5):054904, 2021
work page 2021
-
[36]
J. E. Parkkila, A. Onnerstad, S. F. Taghavi, C. Mor- dasini, A. Bilandzic, M. Virta, and D. J. Kim. New constraints for QCD matter from improved Bayesian parameter estimation in heavy-ion collisions at LHC. Phys. Lett. B, 835:137485, 2022
work page 2022
-
[37]
A. Bazavov et al. Equation of state in ( 2+1 )-flavor QCD.Phys. Rev. D, 90:094503, 2014
work page 2014
-
[38]
Particlization in hybrid models.Eur
Pasi Huovinen and Hannah Petersen. Particlization in hybrid models.Eur. Phys. J. A, 48:171, 2012
work page 2012
-
[39]
Bernhard.Bayesian parameter estimation for relativistic heavy-ion collisions
Jonah E. Bernhard.Bayesian parameter estimation for relativistic heavy-ion collisions. PhD thesis, Duke U., 4 2018
work page 2018
-
[40]
Coupling Relativistic Viscous Hydrodynamics to Boltzmann Descriptions
Scott Pratt and Giorgio Torrieri. Coupling Relativistic Viscous Hydrodynamics to Boltzmann Descriptions. Phys. Rev. C, 82:044901, 2010
work page 2010
-
[41]
S. A. Bass et al. Microscopic models for ultrarela- tivistic heavy ion collisions.Prog. Part. Nucl. Phys., 41:255–369, 1998
work page 1998
-
[42]
M. Bleicher et al. Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model.J. Phys. G, 25:1859–1896, 1999
work page 1999
-
[43]
Flow analysis with cumulants: Direct calcu- lations.Phys
Ante Bilandzic, Raimond Snellings, and Sergei Voloshin. Flow analysis with cumulants: Direct calcu- lations.Phys. Rev. C, 83:044913, 2011
work page 2011
-
[44]
Generic framework for anisotropic flow analyses with multiparticle azimuthal correlations.Phys
Ante Bilandzic, Christian Holm Christensen, Krist- jan Gulbrandsen, Alexander Hansen, and You Zhou. Generic framework for anisotropic flow analyses with multiparticle azimuthal correlations.Phys. Rev. C, 89(6):064904, 2014
work page 2014
-
[45]
Jiangyong Jia, Mingliang Zhou, and Adam Trzupek. Revealing long-range multiparticle collectivity in small collision systems via subevent cumulants.Phys. Rev. C, 96(3):034906, 2017
work page 2017
-
[46]
Giuliano Giacalone, Aleksas Mazeliauskas, and Sören Schlichting. Hydrodynamic attractors, initial state energy and particle production in relativistic nuclear collisions.Phys. Rev. Lett., 123(26):262301, 2019
work page 2019
- [47]
-
[48]
Gardim, Giuliano Giacalone, Matthew Luzum, and Jean-Yves Ollitrault
Fernando G. Gardim, Giuliano Giacalone, Matthew Luzum, and Jean-Yves Ollitrault. Thermodynamics of hot strong-interaction matter from ultrarelativistic nuclear collisions.Nature Phys., 16(6):615–619, 2020. 7
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.