Pith. sign in

REVIEW 4 major objections 5 minor 69 references

Topological analysis of scale-invariant spatial fluctuations in ultrarelativistic heavy-ion collisions

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read A two-stage topological machine-learning pipeline can recover the critical intermittency signal that standard factorial moments wash out in heavy-ion collisions.

desk verdict A useful new benchmark for TDA-ML intermittency recovery in (eta, phi), but the headline recovery is tuned to the reference and needs reframing. read the letter →

arxiv 2608.08259 v1 pith:ANWMHRCP submitted 2026-08-08 hep-ph hep-ex

classification hep-phhep-ex
keywords intermittencycriticalfluctuationstopologicaldataanalysispersistenthomologymachinelearningheavy-ioncollisionsfactorialmomentsQCDpoint
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper argues that scale-invariant critical fluctuations, even when they contribute only a few percent of the final-state tracks, can be separated from an overwhelming thermal background by a two-stage pipeline that first classifies whole events and then prunes individual tracks. In the first stage, each event's $(\eta,\varphi)$ particle distribution is treated as a point cloud, and Betti curves derived from a Delaunay sub-level-set filtration are corrected for multiplicity bias by azimuthal randomization; a convolutional network and a boosted-tree ensemble then separate background events from events containing an injected 5% critical signal. The second stage keeps only tracks with nearest-neighbour distance below 0.02 rad, stripping away the diffuse background and leaving the dense critical clusters. On this filtered sample the normalized factorial moments resume their power-law rise, and the recovered intermittency index $\varphi_2\approx0.73$ matches the pure critical reference value $\varphi_2\approx0.75$. If this works on real collision data, it would provide a tool for probing the QCD critical point that conventional factorial-moment analysis cannot.

What carries the argument

The load-bearing object is the $\Delta$-Betti curve: for each event, a Delaunay triangulation of the $(\eta,\varphi)$ particle cloud is filtered by nearest-neighbour distance, the counts of connected components ($\beta_0$) and loops ($\beta_1$) are tracked against filtration scale $\varepsilon$, and the same curves are recomputed after randomly shuffling azimuthal angles. Subtracting the randomized baseline leaves a fingerprint of genuine spatial clustering, and the concatenated 600-dimensional $\Delta$-Betti vector is the sole input to the classifiers. The second-stage selector is the particle-level nearest-neighbour cut $\varepsilon_{\rm cut}=0.02$ rad, which physically isolates the dense Lévy clusters from the thermal bulk.

What would settle it

Run the full two-stage pipeline on pure EPOS events with no injected CMC signal, applying the same $\varepsilon_{\rm cut}=0.02$ rad particle filter, and measure $\varphi_2$ from the surviving tracks; if $\varphi_2$ rises substantially above zero, the density cut is manufacturing the power law it is supposed to recover.

Watch

Extended reading notes

Core claim

The central discovery is that the dilution that hides critical intermittency can be undone geometrically. Standard normalized factorial moments $F_2(M)$ of a 5% CMC-injected mixture collapse to the inert background value $\varphi_2\approx0$, because the non-critical tracks outnumber the critical clusters. After an event-level topological classifier retains only high-score events, and after a particle-level density cut $d_{\rm NN}\le0.02$ rad removes the diffuse thermal bulk, the same factorial moments recover a power law and give $\varphi_2\approx0.73$, consistent with the pure CMC reference $\varphi_2\approx0.75$. The paper also establishes a calibration point for this geometry: the pure CMC value in $(\eta,\varphi)$ space is about 0.75, shifted from the theoretical $2/3$, which the paper notes is a prediction for $(p_x,p_y)$ space, not angular space.

Load-bearing premise

The recovery rests on the assumption that the 0.02-radian density threshold can be fixed from the pure CMC reference without peeking at the injected signal, and that a Lévy walk with exponent 1/6 performed directly in (η, φ) is the right model for critical angular fluctuations.

Editorial extensions

If this is right

  • At a 5% injected signal fraction, event-level topological classification reaches an AUC of about 0.99 with both a convolutional network and a boosted-tree ensemble, while performance degrades to about 0.95 at 3% and about 0.8 at 1%.
  • Factorial moments computed on event-filtered samples stay flat, so the particle-level density cut is both necessary and sufficient to restore the critical power-law scaling.
  • Applying the density cut to the pure CMC sample does not distort its scaling, indicating that the recovered index reflects the injected signal rather than an artifact of the filter.
  • The pure CMC baseline in $(\eta,\varphi)$ space is $\varphi_2\approx0.75$, so future intermittency searches in angular coordinates should compare against this shifted reference rather than the $(p_x,p_y)$ value of $2/3$.
  • The same pipeline can be recalibrated with dynamically generated critical configurations, such as the successive-contraction-and-randomization model, to test whether the recovery survives more realistic phase evolution.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the choice of $\varepsilon_{\rm cut}=0.02$ rad was validated against the known pure CMC reference, so applying this pipeline to real data, where no pure reference exists, will require an independent calibration of the density threshold or the method risks circularity.
  • Editorial inference: because the paper itself notes that the theoretical $\varphi_2=2/3$ applies to $(p_x,p_y)$ and not $(\eta,\varphi)$, the recovered values 0.73 and 0.75 are geometric calibration numbers; testing the Ising prediction directly would require running the same $\Delta$-Betti pipeline in a cumulatively flattened momentum space.
  • Editorial inference: a sharp testable extension is to run the second-stage density filter on pure EPOS events and on a Poissonian toy background with no injected clusters; if $\varphi_2$ rises substantially above zero, the filter itself would be creating apparent intermittency.
  • Editorial inference: the nearest-neighbour threshold sits exactly in the angular-separation regime where detector track merging and splitting occur, so simulating a realistic detector response is the immediate next step before this can be applied to experimental data.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript presents the first Critical Monte Carlo (CMC) based intermittency analysis in the two-dimensional angular (η,φ) phase space, using EPOS as the background model for Pb–Pb collisions at √sNN = 5.02 TeV. The authors embed CMC signal tracks at fractions λ = 1%, 3%, and 5%, show that standard normalized factorial moments (NFMs) are diluted by the background, and then propose a two-stage topological machine learning pipeline: event-level classification with TopoPointNet and XGBoost on Δ-Betti curves, followed by a particle-level nearest-neighbor density filter with threshold ε_cut. They report that this pipeline restores power-law scaling of NFMs and recovers an intermittency index φ2 ≈ 0.73 for the filtered 5% mixture, in close agreement with their pure CMC reference of φ2 ≈ 0.75. The paper also reports classification AUCs ~0.99, SHAP/gain feature attributions, and a discussion of limitations for experimental deployment.

Significance. If the recovery claim were validated out of sample, the paper would be a useful methodological contribution: it extends CMC-based intermittency studies from (p_x,p_y) to the experimentally more accessible (η,φ) plane, and it combines persistent homology, two independent classifiers, and explicit multiplicity-bias corrections. The periodic Delaunay construction with phantom points, the Euler-characteristic shortcut for β1, the use of azimuthal randomization, and the candid limitation paragraph are all strengths. However, the central quantitative claim is currently a benchmark fit rather than an out-of-sample measurement: the particle-level cut ε_cut is selected after inspecting which value best reproduces the pure CMC reference. Because the surviving tracks after such a density cut are mostly the injected CMC particles, the NFM recovery is to a substantial degree built into the construction. The paper needs a pre-specified threshold selection rule or an independent validation set, plus a consistent statement of the pure CMC reference value and its uncertainty.

major comments (4)
  1. [§VI.3 and Figs. 12–13] The central recovery claim rests on a post hoc choice of ε_cut. The text states that among the tested thresholds, ‘the ε_cut ≥ 0.02 sample lies closest to the pure CMC reference,’ and Figs. 12–13 show that this is how the operating point was identified. Because the threshold is selected by closeness to the benchmark it is supposed to recover, the reported φ2 ≈ 0.73 is a fit on the benchmark, not an out-of-sample measurement. Please specify the threshold before any comparison with the pure CMC reference, or validate the threshold on a held-out set of mixtures and then apply it to the 5% sample without reference to the outcome.
  2. [§II, Fig. 4 vs. §VI.3 and Conclusions] The pure CMC reference value is quoted inconsistently: Fig. 4 gives φ2 ≈ 0.73, while §VI.3 and the Conclusions give φ2 ≈ 0.75, with no uncertainty on either. If the two values come from different fit ranges or different CMC samples, that must be stated; otherwise this is an internal inconsistency that directly affects the claimed 1% residual discrepancy between the recovered index and the reference.
  3. [§VI.3 and Sec. IV.1] The particle-level filter uses the same nearest-neighbor distance d_NN that defines the CMC clusters: after retaining tracks with d_NN ≤ 0.02 rad, the surviving sample is dominated by the injected CMC tracks by construction. This makes the NFM recovery partly mechanical, independent of the ML stage. Please quantify the post-filter signal purity and test a null control in which the same number of tracks is selected by a random or density-mismatched rule, to show that the restored power law is not an artifact of the cut itself.
  4. [§II and Sec. VI.3] The paper correctly notes that the theoretical value φ2 = 2/3 is derived for (p_x,p_y) space, not (η,φ), yet the CMC Lévy walk with μ = 1/6 is used directly in (η,φ) as a critical benchmark. The recovered index is therefore calibrated against a model-specific reference whose relation to QCD criticality in angular space is not independently established. This is a correctness-risk concern rather than a circularity claim; a concrete test would be to compare against the SCR model mentioned in the Conclusions, or to derive the expected angular-plane intermittency index for a critical system before using the CMC value as ground truth.
minor comments (5)
  1. [Abstract and Introduction] There are several typographical errors, including ‘a extremely hot’ and ‘behaviuor’; these should be corrected before publication.
  2. [§II, Eq. (1)] The event-averaging convention in Eq. (1) is somewhat ambiguous as written; please clarify that the factorial moments are computed per event and then averaged, or provide the standard Bialas–Peschanski expression explicitly.
  3. [Fig. 11 caption] The caption contains a grammatical issue: ‘which exhibits a secondary peaks’ should read ‘which exhibit secondary peaks.’
  4. [References] Reference [12] contains an unresolved ‘?’ placeholder, and §VI.1 refers to ‘Section VI 2’ without a space; both should be fixed.
  5. [§VI.3] No statistical uncertainties are quoted for the recovered φ2 values, and Figs. 4 and 13 show no error bars; even a bootstrap or fit-uncertainty estimate would materially strengthen the comparison between the filtered mixture and the pure CMC reference.

Circularity Check

2 steps flagged · score 6.0 of 10

The headline recovery of φ2≈0.73 is a benchmark fit rather than an out-of-sample prediction: the particle-level cut ε_cut=0.02 is selected after comparison with the pure-CMC reference, and the filter is defined by the same local-density metric that distinguishes the injected CMC clusters.

  1. fitted input called prediction [Sec. VI.3 (Recovery of the NFMs), Figs. 12–13]
    "Re-evaluating the NFMs on this filtered track sample yields a clear restoration of the power-law scaling. Fig. 12 demonstrate the scaling behaviour of ln F2(M) as a function of ln M2 across the EPOS background, the pure CMC reference and the unfiltered 5% signal mixture, along with the fully filtered pipeline samples at different filtration thresholds, ε_cut. Among these, the ε_cut ≥ 0.02 sample lies closest to the pure CMC reference, where a looser cut of 0.05 retains an excess of background tracks."

    The central quantitative claim — recovered φ2≈0.73 in close agreement with pure CMC φ2≈0.75 — is established by selecting the cut that reproduces the reference, not by a pre-specified rule fixed before seeing the reference. The text reports that the 0.02 rad sample 'lies closest to the pure CMC reference'; Figs. 12–13 are then used as the demonstration of success. The SHAP/gain motivation cited for the cut is broad (0.02–0.05 rad) and does not determine the operating point. Thus the 'accurate recovery' is a fit to the benchmark, and the agreement is forced by the selection.

  2. self definitional [Sec. VI.3; Sec. IV.1; Conclusions (limitations paragraph)]
    "Within these signal enriched events, individual tracks are retained exclusively if their local nearest-neighbour distance satisfies d_NN ≤ 0.02 rad. This geometric threshold explicitly strips away the diffuse thermal bulk, isolating the densely packed particles comprising the embedded critical clusters."

    By construction, the CMC signal is a 'localized CMC Lévy walk' whose 'tightly packed critical clusters yield small d_NN values', while the EPOS background is a 'uniform distribution' with large d_NN. The particle-level filter is therefore a selector for exactly the injected CMC particles, and the paper's own limitations paragraph concedes that the classifier and the filter 'are not strictly statistically independent' because both use 'the identical d_NN metric'. Computing NFMs on the surviving tracks is equivalent to recomputing the NFM of the input pure-CMC cluster; the restored power law is inherited from the injection rather than independently derived.

full rationale

The first stage of the pipeline (event-level classification of Δ-Betti curves) is evaluated on a held-out test set with AUC≈0.99 and is not circular in itself. The circularity is concentrated in the second stage, which carries the paper's central claim. The particle-level cut ε_cut=0.02 rad is selected after comparing different cuts with the pure-CMC reference (Fig. 12: 'the ε_cut ≥ 0.02 sample lies closest to the pure CMC reference'), and the same Sec. VI.3 then uses that cut to report φ2≈0.73 'in close agreement with the pure CMC reference' (φ2≈0.75). This is benchmark fitting: no pre-registered rule fixes the operating point, and the quoted agreement is the selection criterion. Moreover, the filter is defined by nearest-neighbour distance, which is the same property used to characterize the CMC Lévy clusters as 'tightly packed' and the EPOS background as uniform; the paper's limitations paragraph admits both stages rest on 'the identical d_NN metric'. The surviving tracks are therefore predominantly the injected CMC particles, so the 'restored' NFM power law is the input signal's power law by construction. The pure-CMC reference is also generated by the same CMC model injected into the mixture, so the recovery target is not an external first-principles prediction; the paper itself notes the theoretical φ2=2/3 applies to (p_x,p_y), not (η,φ). These points do not invalidate the classification study, but they mean the headline 'accurate recovery' reduces to a self-selected benchmark plus a density selector matched to the injected signal. Score 6 reflects partial circularity in the central recovery claim; no load-bearing self-citation or uniqueness import is present.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

The benchmark is built on several chosen parameters and physical assumptions. The most important free parameter is the density cut epsilon_cut = 0.02 rad, which is tuned to recover the known CMC reference. The physical axioms concern the validity of CMC geometry in angular space and the inertness of EPOS background. No invented physical entities are introduced.

free parameters (6)
  • lambda (signal replacement fraction) = 0.01, 0.03, 0.05 (operational 0.05)
    Controls the injected signal fraction; the paper identifies 5% as the threshold for effective event-level discrimination.
  • epsilon_cut (particle-level density filter) = 0.02 rad
    Selected because it produces the closest recovery of the pure CMC phi2; a looser cut of 0.05 degrades recovery.
  • S threshold (event-level score threshold) = 0.90
    Strict operating point chosen for high signal purity; efficiencies are quoted at this threshold.
  • N_epsilon (number of filtration steps) = 300
    Chosen because classification AUC saturates beyond 300 steps.
  • NFM fit range = ln M^2 in [6.9, 9.1]
    Linear fit range for the intermittency index, selected in the higher-M^2 region without independent justification.
  • Levy walk step-size bounds = r_min = r_max x 10^-7, r_max = 2 pi
    Boundary values of the Levy walk; the lower bound is a numerical regularization.
assumptions (5)
  • domain assumption CMC Levy walk with exponent mu = 1/6 in (eta, phi) represents critical fluctuations of the 3D Ising universality class.
    Section II; the theoretical value phi2 = 2/3 is derived for (p_x, p_y), not angular space, and the paper does not derive the angular-space analogue.
  • domain assumption EPOS events contain no critical fluctuations and serve as an inert thermal background.
    Section III; based on earlier EPOS intermittency studies cited as [34, 35].
  • domain assumption Subtracting azimuthally randomized Betti curves removes all trivial multiplicity and flow effects, isolating dynamical correlations.
    Section IV.3; assumes the randomized control preserves the same distribution shape except for azimuthal correlations.
  • standard math The Delaunay filtration and Euler-Poincare formula exactly capture the homology of the event point cloud.
    Section IV; standard TDA result for 2D simplicial complexes.
  • ad hoc to paper Track-level nearest-neighbor distance maps inversely to local density and can separate signal clusters from background.
    Section VI.3; the cut is applied at 0.02 rad with no independent calibration.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Topological analysis of scale-invariant spatial fluctuations in ultrarelativistic heavy-ion collisions." pith.science (2026). https://pith.science/paper/ANWMHRCP

@misc{pith2026260808259,
  author       = {Pith},
  title        = {Pith review of: Topological analysis of scale-invariant spatial fluctuations in ultrarelativistic heavy-ion collisions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ANWMHRCP}},
  note         = {Machine review of arXiv:2608.08259}
}
abstract

The QGP-to-hadronic matter phase transition and QCD critical point in heavy-ion collisions can be identified by studying spatial fluctuations among final-state particles using intermittency analysis. First CMC-based intermittency analysis in the two-dimensional angular ($\eta$, $\varphi$) phase space, using EPOS as the background model is presented. Critical fluctuation signals are extremely weak, constituting only a few percent of the total event sample and are severely diluted by the overwhelming non-critical background, rendering traditional intermittency analyses insufficient for reliable signal extraction. To extract the weak critical signal, we employ a two-stage topological machine learning framework combining Topological Data Analysis (TDA) with deep learning. In the first stage, particle events are represented as two-dimensional point clouds and a Delaunay-based sub-level set filtration is constructed to extract Betti curves as multiscale topological invariants, corrected for multiplicity bias via azimuthal randomisation and classified by two complementary architectures, a TopoPointNet (TPN) and Boosted Decision Trees (BDT). Since event-level classification alone is insufficient to restore the critical scaling, a second stage applies a particle-level density filter, explicitly stripping away the diffuse thermal background and isolating the densely packed critical clusters. The two stage pipeline successfully restores the power-law scaling of the normalized factorial moments, enabling accurate recovery of the intermittency index in ($\eta$, $\varphi$) space and establishing topological machine learning as a robust data driven tool for probing the QCD critical point and the phase structure of strongly interacting matter in heavy-ion collisions at LHC energies.

Figures

Figures reproduced from arXiv: 2608.08259 by the authors.

Figure 1
Figure 1. FIG. 1. Two-dimensional spatial track distribution in the [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Log-log dependence of the second-order NFM, [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Scaling dilution of second-order NFM, [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (8 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Evolution of the 2D periodic Delaunay filtration in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Multiplicity-normalized raw Betti curves [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Azimuthally randomized differential Betti curves [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Output classifier score distributions [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. ROC performance curves for TopoPointNet (red) and BDT (green) across signal replacement fractions 1% (left, AUC [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 10
Figure 10. Figure 10: shows the classification performance quantified by the signal efficiency 𝜀sig and background rejection 1 − 𝜀bkg as a function of the classifier score threshold. Both architec￾tures exhibit strong separation over the full threshold range. At a loose threshold of S ≥ 0.…
Figure 11
Figure 11. Figure 11: FIG. 11. Topological feature importance across Delaunay filtration scale, [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12. Restoration of intermittency scaling, [PITH_FULL_IMAGE:figures/full_fig_p013_12.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

69 extracted references · 37 canonical work pages

  1. [1]

    NFMs are calculated for the event samples and shown in Figure 2 as a log-log dependence of𝐹 2(M) on the number of phase-space bins𝑀 2

    Dilution in the embedded signal The section details the sensitivity of intermittency to the embedded signal. NFMs are calculated for the event samples and shown in Figure 2 as a log-log dependence of𝐹 2(M) on the number of phase-space bins𝑀 2. For the pure EPOS background (left),𝐹 2(M) exhibit no significant dependence on the increasing phase space bins. ...

  2. [2]

    Delaunay filtration In order to analyse the multi-particle correlations within an event, topology relies on the construction of a simplicial com- plex. A simplicial complex generalizes a standard network graph by incorporating higher-dimensional geometric struc- tures, for a two-dimensional point cloud, this comprises ver- tices (0-simplices), edges (1-si...

  3. [3]

    The ranks of these homology groups are the Betti num- bers, which quantify the distinct topological features present at a given spatial resolution

    Extraction of Betti curves The sequence of subcomplexes generated by the sub-level- set filtration is topologically characterized by its homology groups. The ranks of these homology groups are the Betti num- bers, which quantify the distinct topological features present at a given spatial resolution. For the two-dimensional(𝜂,𝜑) phase space, the relevant ...

  4. [4]

    A machine learning classifier trained on these raw Betti curves could exploit this amplitude scaling to distin- guish classes by multiplicity rather than by cluster geometry

    Azimuthal randomisation in betti curves 𝛽0(𝜀)and𝛽 1(𝜀)curves scale with the total number of ac- tive particles and hence, with the event track multiplicity𝑁, which fluctuates across collisions even within a narrow cen- trality class. A machine learning classifier trained on these raw Betti curves could exploit this amplitude scaling to distin- guish class...

  5. [5]

    By treating the discrete scale parameter𝜀as an ordered sequence, the topological features are naturally processed by a one- dimensional convolutional neural network (1D CNN) [54]

    TopoPointNet (TPN) The delta Betti curves𝛽 Δ 0(𝜀)and𝛽 Δ 1(𝜀)constitute ordered sequences of topological information, tracking the structural evolution of the event across the spatial filtration scale. By treating the discrete scale parameter𝜀as an ordered sequence, the topological features are naturally processed by a one- dimensional convolutional neural...

  6. [6]

    Boosted Decision Trees (BDT) As an alternative architecture to the deep convolutional network, the eXtreme Gradient Boosting (XGBoost) algo- rithm [53] is used. The foundational element of this model is the decision tree, which partitions the topological feature space into discrete regions through a sequence of binary threshold splits, assigning a constan...

  7. [7]

    Figure 8 shows the resulting score distributions for signal (𝜆=0.05) and back- ground, normalized to unit area (probability density)

    Classification performance A continuous classification score bounded within the in- terval is evaluated, representing the likelihood of an event containing the embedded critical signal. Figure 8 shows the resulting score distributions for signal (𝜆=0.05) and back- ground, normalized to unit area (probability density). For both the TPN and the BDT ensemble...

  8. [8]

    Classifier feature attribution To interpret the physical mechanisms driving the BDT clas- sification decisions, the global feature importance is evaluated as a function of the filtration parameter𝜀. Two complemen- tary interpretability metrics are employed, the mean absolute Shapley value [55] averaged across the test ensemble for each discreteΔ-Betti fea...

Show all 69 references
  1. [9]

    Indian participation in the ALICE experiment at CERN,

    Recovery of the NFMs To quantify the recovery of the critical signal, the bench- marks by computing have been established in Figure 3,4. Ap- plying the primary stage of the topological pipeline filters the event ensemble, retaining only those collision configurations that sati...

  2. [10]

    E. V. Shuryak, Quark-Gluon Plasma and Hadronic Production of Leptons, Photons and Psions, Phys. Lett. B78, 150 (1978)

  3. [11]

    Niida and Y

    T. Niida and Y. Miake, Signatures of QGP at RHIC and the LHC, AAPPS Bull.31, 12 (2021), arXiv:2104.11406 [nucl-ex]

  4. [12]

    Woithe, G

    J. Woithe, G. J. Wiener, and F. F. Van der Veken, Let’s have a coffee with the Standard Model of particle physics!, Phys. Educ. 52, 034001 (2017)

  5. [13]

    D. P. Barberet al., Discovery of Three Jet Events and a Test of Quantum Chromodynamics at PETRA Energies, Phys. Rev. Lett.43, 830 (1979)

  6. [14]

    V. N. Gribov, The Theory of quark confinement, Eur. Phys. J. C 10, 91 (1999), arXiv:hep-ph/9902279

  7. [15]

    D. J. Gross and F. Wilczek, Ultraviolet Behavior of Nonabelian Gauge Theories, Phys. Rev. Lett.30, 1343 (1973)

  8. [16]

    D. J. Gross and F. Wilczek, Asymptotically Free Gauge Theories - I, Phys. Rev. D8, 3633 (1973)

  9. [17]

    D. J. Gross and F. Wilczek, Asymptotically Free Gauge Theo- ries 2, Phys. Rev. D9, 980 (1974)

  10. [18]

    S. S. Adleret al.(PHENIX), J / psi production in Au Au col- lisions at s(NN)**(1/2) = 200-GeV at the Relativistic Heavy Ion Collider, Phys. Rev. C69, 014901 (2004), arXiv:nucl- ex/0305030

  11. [19]

    Koch and V

    V. Koch and V. Vovchenko, Exploring the QCD phase diagram through correlations and fluctuations, Eur. Phys. J. Spec. Top. 10.1140/epjs/s11734-026-02307-w (2025), arXiv:2512.04288 [nucl-th]

  12. [20]

    Chenget al., The Transition temperature in QCD, Phys

    M. Chenget al., The Transition temperature in QCD, Phys. Rev. D74, 054507 (2006), arXiv:hep-lat/0608013

  13. [21]

    Odyniec (STAR), Beam Energy Scan Program at RHIC (BES I and BES II) – Probing QCD Phase Diagram with Heavy-Ion Collisions, PoSCORFU2018, 151 (2019)

    G. Odyniec (STAR), Beam Energy Scan Program at RHIC (BES I and BES II) – Probing QCD Phase Diagram with Heavy-Ion Collisions, PoSCORFU2018, 151 (2019)

  14. [22]

    Baym and H

    G. Baym and H. Heiselberg, Event-by-event fluctuations in ul- trarelativistic heavy ion collisions, Phys. Lett. B469, 7 (1999), arXiv:nucl-th/9905022

  15. [23]

    V. Koch, M. Bleicher, and S. Jeon, Event-by-event fluctua- tions and the QGP, Nucl. Phys. A698, 261 (2002), arXiv:nucl- th/0103084

  16. [24]

    Bialas and R

    A. Bialas and R. B. Peschanski, Intermittency in Multiparticle Production at High-Energy, Nucl. Phys. B308, 857 (1988)

  17. [25]

    E. A. De Wolf, I. M. Dremin, and W. Kittel, Scaling laws for density correlations and fluctuations in multiparticle dynamics, Phys. Rept.270, 1 (1996), arXiv:hep-ph/9508325

  18. [26]

    Bialas and R

    A. Bialas and R. B. Peschanski, Moments of Rapidity Distribu- tions as a Measure of Short Range Fluctuations in High-Energy Collisions, Nucl. Phys. B273, 703 (1986)

  19. [27]

    R. C. Hwa and C. B. Yang, Observable Properties of Quark- Hadron Phase Transition at the Large Hadron Collider, Acta Phys. Polon. B48, 23 (2017), arXiv:1601.04671 [nucl-th]

  20. [28]

    N. G. Antoniou, F. K. Diakonos, A. S. Kapoyannis, and K. S. Kousouris, Critical opalescence in baryonic QCD matter, Phys. Rev. Lett.97, 032002 (2006), arXiv:hep-ph/0602051

  21. [29]

    N. G. Antoniou, Y. F. Contoyiannis, F. K. Diakonos, and G. Mavromanolakis, Critical QCD in nuclear collisions, Nucl. Phys. A761, 149 (2005), arXiv:hep-ph/0505185

  22. [30]

    J. Wu, Z. Li, X. Luo, M. Xu, and Y. Wu, Intermittency of charged particles in the hybrid UrQMD+CMC model at energies available at the BNL Relativistic Heavy Ion Collider, Phys. Rev. C106, 054905 (2022), arXiv:2209.07135 [nucl-th]

  23. [31]

    V. Z. Reyna Ortiz (NA61/SHINE), Search for the Critical Point via Intermittency Analysis in NA61/SHINE, Ukr. J. Phys.69, 858 (2024)

  24. [32]

    Podlaski, NA61/SHINE Overview, EPJ Web Conf.296, 01008 (2024), arXiv:2402.10973 [nucl-ex]

    P. Podlaski, NA61/SHINE Overview, EPJ Web Conf.296, 01008 (2024), arXiv:2402.10973 [nucl-ex]

  25. [33]

    Adhikaryet al.(NA61/SHINE), Search for a critical point of strongly-interacting matter in central 40Ar + 45Sc collisions at 13 A–75 A GeV/c beam momentum, Eur

    H. Adhikaryet al.(NA61/SHINE), Search for a critical point of strongly-interacting matter in central 40Ar + 45Sc collisions at 13 A–75 A GeV/c beam momentum, Eur. Phys. J. C84, 741 (2024), arXiv:2401.03445 [nucl-ex]

  26. [34]

    Abdulhamidet al.(STAR), Energy dependence of intermit- tency for charged hadrons in Au+Au collisions at RHIC, Phys

    M. Abdulhamidet al.(STAR), Energy dependence of intermit- tency for charged hadrons in Au+Au collisions at RHIC, Phys. Lett. B845, 138165 (2023), arXiv:2301.11062 [nucl-ex]

  27. [35]

    B. B. Abelevet al.(ALICE), Performance of the ALICE Ex- periment at the CERN LHC, Int. J. Mod. Phys. A29, 1430044 (2014), arXiv:1402.4476 [nucl-ex]

  28. [36]

    Sharma and R

    S. Sharma and R. Gupta (ALICE), Local Multiplicity Fluc- tuations in Pb−Pb Collisions at √𝑠NN = 2.76 TeV with AL- ICE at the LHC, Springer Proc. Phys.304, 860 (2024), arXiv:2307.14407 [nucl-ex]

  29. [37]

    S. K. Malik, S. Sharma, and R. Gupta (ALICE), Event-by-event local multiplicity fluctuations in charged particle production at the LHC energies with ALICE, DAE Symp. Nucl. Phys.67, 1003 (2024)

  30. [38]

    R. C. Hwa and M. T. Nazirov, Intermittency in second order phase transition, Phys. Rev. Lett.69, 741 (1992)

  31. [39]

    N. G. Antoniou, Y. F. Contoyiannis, F. K. Diakonos, and C. G. Papadopoulos, Fractals at T = T(c) due to instanton - like configurations, Phys. Rev. Lett.81, 4289 (1998), arXiv:hep- ph/9810383

  32. [40]

    N. G. Antoniou, Y. F. Contoyiannis, F. K. Diakonos, A. I. Karanikas, and C. N. Ktorides, Pion production from a crit- ical QCD phase, Nucl. Phys. A693, 799 (2001), arXiv:hep- ph/0012164

  33. [41]

    R. C. Hwa and C. B. Yang, Local Multiplicity Fluctuations as a Signature of Critical Hadronization at LHC, Phys. Rev. C85, 044914 (2012), arXiv:1111.6651 [nucl-th]

  34. [42]

    Sharma and R

    R. Sharma and R. Gupta, Scaling Properties of Multiplicity Fluctuations in the AMPT Model, Adv. High Energy Phys.2018, 6283801 (2018), arXiv:1806.10854 [hep-ph]

  35. [43]

    Gupta and S

    R. Gupta and S. K. Malik, Intermittency study of charged particles generated in Pb-Pb collisions at √𝑠NN= 2.76 TeV us- ing EPOS3, Adv. High Energy Phys.2020, 5073042 (2020), arXiv:1911.13111 [hep-ex]

  36. [44]

    F. U. Haider, S. K. Malik, R. Gupta, and B. Singh, Study of the charged particle multiplicity fluctutions in EPOS4 for Pb-Pb 𝑠𝑞𝑟𝑡𝑠 𝑁𝑁 = 5.02 TeV, DAE Symp. Nucl. Phys.69, 1081 (2026)

  37. [45]

    Sarma and B

    P. Sarma and B. Bhattacharjee, Color reconnection as a possible mechanism of intermittency in the emission spectra of charged particles in PYTHIA-generated high-multiplicity𝑝𝑝collisions at energies available at the CERN Large Hadron Collider, Phys. Rev. C99, 034901 (2019), arX...

  38. [46]

    Singh, A

    A. Singh, A. Kumar, A. Chandra, S. Singh, S. Ahmad, and B. K. Singh, Scaling properties of particle density fluctuations at LHC energies, EPL148, 14001 (2024)

  39. [47]

    Sharma, S

    S. Sharma, S. K. Malik, Z. Banoo, and R. Gupta, Normalized factorial moments of spatial distributions of particles in high multiplicity events: A Toy model study, Nucl. Phys. A1053, 122963 (2025), arXiv:2309.07712 [hep-ph]

  40. [48]

    K. Werner, Parallel scattering, saturation, and generalized Abramovskii-Gribov-Kancheli (AGK) theorem in the EPOS4 framework, with applications for heavy-ion collisions at sNN of 5.02 TeV and 200 GeV, Phys. Rev. C109, 034918 (2024), arXiv:2310.09380 [hep-ph]

  41. [49]

    Werner, Epos4: New theoretical concepts for modeling 16 proton-proton and ion-ion scattering at very high energies (2024), arXiv:2410.09955 [hep-ph]

    K. Werner, Epos4: New theoretical concepts for modeling 16 proton-proton and ion-ion scattering at very high energies (2024), arXiv:2410.09955 [hep-ph]

  42. [50]

    K. Werner, Revealing a deep connection between factorization and saturation: New insight into modeling high-energy proton- proton and nucleus-nucleus scattering in the EPOS4 framework, Phys. Rev. C108, 064903 (2023), arXiv:2301.12517 [hep-ph]

  43. [51]

    S. A. Basset al., Microscopic models for ultrarelativistic heavy ion collisions, Prog. Part. Nucl. Phys.41, 255 (1998), arXiv:nucl-th/9803035

  44. [52]

    L. D. McLerran and R. Venugopalan, Gluon distribution func- tions for very large nuclei at small transverse momentum, Phys. Rev. D49, 3352 (1994), arXiv:hep-ph/9311205

  45. [53]

    J. Wu, Y. Lin, Z. Li, X. Luo, and Y. Wu, Intermittency analysis of proton numbers in heavy-ion collisions at energies available at the BNL Relativistic Heavy Ion Collider, Phys. Rev. C104, 034902 (2021), arXiv:2104.11524 [nucl-th]

  46. [54]

    Abbaset al.(ALICE), Performance of the ALICE VZERO system, JINST8, P10016, arXiv:1306.3130 [nucl-ex]

    E. Abbaset al.(ALICE), Performance of the ALICE VZERO system, JINST8, P10016, arXiv:1306.3130 [nucl-ex]

  47. [55]

    Centrality determination in heavy ion collisions, (2018)

  48. [56]

    Anticicet al.(NA49), Critical fluctuations of the proton den- sity in A+A collisions at 158𝐴GeV, Eur

    T. Anticicet al.(NA49), Critical fluctuations of the proton den- sity in A+A collisions at 158𝐴GeV, Eur. Phys. J. C75, 587 (2015), arXiv:1208.5292 [nucl-ex]

  49. [57]

    Leykam and D

    D. Leykam and D. G. Angelakis, Topological data analy- sis and machine learning, Adv. Phys. X8, 2202331 (2023), arXiv:2206.15075 [cond-mat.mes-hall]

  50. [58]

    Hamilton, T

    G. Hamilton, T. Dore, and C. Plumberg, Applications of persis- tent homology in nuclear collisions, Phys. Rev. C106, 064912 (2022), arXiv:2209.15480 [nucl-th]

  51. [59]

    Capellino, A

    F. Capellino, A. Dubla, S. Masciocchi, G. Nijs, and D. Spitz, To- ward a topological data analysis for heavy-ion collisions, Phys. Rev. C112, 054909 (2025), arXiv:2509.02339 [nucl-th]

  52. [60]

    R. Wang, C. Qiu, C.-S. Hu, Z. Li, and Y. Wu, Identifying weak critical fluctuations of intermittency in heavy-ion collisions with topological machine learning, Phys. Lett. B864, 139405 (2025), arXiv:2412.06151 [nucl-th]

  53. [61]

    T. K. Duy, Y. Hiraoka, and T. Shirai, Limit theorems for persis- tence diagrams (2016), arXiv:1612.08371 [math.PR]

  54. [62]

    Chen and C

    T. Chen and C. Guestrin, Xgboost: A scalable tree boosting system, inProceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, KDD ’16 (Association for Computing Machinery, New York, NY, USA, 2016) p. 785–794

  55. [63]

    Huang, L.-G

    Y. Huang, L.-G. Pang, X. Luo, and X.-N. Wang, Probing crit- icality with deep learning in relativistic heavy-ion collisions, Phys. Lett. B827, 137001 (2022), arXiv:2107.11828 [nucl-th]

  56. [64]

    Lundberg and S.-I

    S. Lundberg and S.-I. Lee, A unified approach to interpreting model predictions (2017), arXiv:1705.07874 [cs.AI]

  57. [65]

    Lippmann (ALICE), Upgrade of the ALICE Time Projection Chamber, (2014)

    C. Lippmann (ALICE), Upgrade of the ALICE Time Projection Chamber, (2014)

  58. [66]

    Bialas and M

    A. Bialas and M. Gazdzicki, A New variable to study intermit- tency, Phys. Lett. B252, 483 (1990)

  59. [67]

    Weinberg,The First Three Minutes

    S. Weinberg,The First Three Minutes. A Modern View of the Origin of the Universe(Basic Books, 1977)

  60. [68]

    Wan, C.-Z

    J. Wan, C.-Z. Wang, Q.-Y. Shou, and Y.-G. Ma, Tracing pT- differential radial flow from blast-wave analytics to quark coa- lescence, Phys. Rev. C113, 064902 (2026), arXiv:2509.24889 [nucl-th]

  61. [69]

    J. Wu, X. Luo, Z. Li, and Y. Wu, [article title], Nuclear Physics Review41, 580 (2024)

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.