pith. machine review for the scientific record. sign in

arxiv: 2604.26731 · v1 · submitted 2026-04-29 · ⚛️ nucl-th · hep-ph· nucl-ex

Recognition: unknown

Thermal and geometric normal modes of spectral fluctuations in heavy-ion collisions

Authors on Pith no claims yet

Pith reviewed 2026-05-07 11:42 UTC · model grok-4.3

classification ⚛️ nucl-th hep-phnucl-ex
keywords heavy-ion collisionsspectral fluctuationsprincipal component analysisthermal fluctuationsgeometric fluctuationselliptic flowquark-gluon plasmanormal modes
0
0 comments X

The pith

Principal component analysis decomposes heavy-ion spectral fluctuations into two dominant thermal and geometric modes.

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

The paper establishes that event-by-event fluctuations in the transverse momentum spectrum of particles from heavy-ion collisions arise from a combination of thermal and geometric effects in the initial state. By performing principal component analysis on the joint covariance of the normalized spectrum, mean transverse momentum, and elliptic flow squared, the authors find that the first two modes capture 99.5 percent of the variance. These modes are then rotated using constraints tied to physical initial-state responses, yielding direct correspondences to observed flow coefficients: the thermal mode accounts for all of the measured v0(pT), while the geometric mode drives much of the v02(pT) behavior including its low-pT sign change. A reader would care because this supplies a concrete, data-driven separation of fluctuation sources that links microscopic initial conditions to macroscopic experimental signals in the quark-gluon plasma.

Core claim

Principal component analysis applied to the joint covariance structure of the normalized transverse-momentum spectrum, mean transverse momentum, and elliptic flow squared isolates two leading modes that together explain 99.5 percent of the total variance. Orthogonal rotation of these modes, performed under constraints motivated by the initial-state thermal and geometric responses, produces normal modes directly analogous to the vibrational modes of a linear triatomic molecule. The resulting thermal mode fully accounts for the experimentally measured v0(pT), while the geometric mode supplies the dominant contribution to v02(pT) in non-central collisions and explains its characteristic low-pT

What carries the argument

Principal component analysis on the joint covariance matrix of normalized spectrum, mean pT and v2 squared, followed by orthogonal rotation constrained by initial-state thermal and geometric responses.

If this is right

  • The thermal mode accounts entirely for the pT dependence of the measured v0 coefficient.
  • The geometric mode supplies the dominant part of v02(pT) in non-central collisions and explains its low-pT sign change.
  • The two-mode decomposition provides a transparent experimental window into the separate thermal and geometric structure of the initial state.
  • The normal-mode analogy to a linear triatomic molecule supplies a simple physical picture for how thermal and geometric fluctuations combine.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same covariance-based decomposition could be extended to higher-order flow harmonics to separate additional fluctuation sources.
  • If the separation holds, hydrodynamic simulations could be validated by checking whether their predicted thermal and geometric mode contributions match the data-driven ones.
  • This framework suggests a route to constrain initial-state models by demanding consistency with both the variance captured and the flow observables reproduced.

Load-bearing premise

The physical constraints used to rotate the PCA modes correctly isolate independent thermal and geometric responses without circular dependence on the fluctuation data itself.

What would settle it

Experimental data in which the rotated thermal mode fails to reproduce the full measured pT dependence of v0 or the geometric mode fails to reproduce the sign change and magnitude of v02 in non-central collisions.

Figures

Figures reproduced from arXiv: 2604.26731 by Rupam Samanta.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Scatter plot between event-by-event mean trans view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Cartoon of the symmetric stretching (coherent) mode of vibration in a linear triatomic molecule, where the middle view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Results for the thermal (red) and geometric (blue) view at source ↗
read the original abstract

The transverse momentum spectrum of charged particles in ultra-relativistic heavy-ion collisions fluctuates event-by-event, encoding signatures of underlying collective dynamics. Such fluctuations originate from a combined effect of thermal and geometric fluctuations in the initial state. We present a direct decomposition of these spectral fluctuations through principal component analysis performed on the joint covariance structure of normalized spectrum, mean transverse momentum and elliptic flow squared. The first two leading modes explain 99.5\% of the total variance, and are orthogonally rotated by imposing physical constraints motivated by the initial state thermal and geometric response. The resulting thermal and geometric modes bear direct analogy with the vibrational normal modes of a linear triatomic molecule. The thermal mode entirely drives the experimentally measured $v_0(p_T)$, while the geometric mode contributes substantially to $v_{02}(p_T)$ in non-central collisions, providing a transparent explanation of its characteristic low-$p_T$ sign change. The study establishes the first physically motivated interpretation of principal component modes in the field of heavy-ion collisions and provides an experimental window into the thermo-geometric structure of the QGP initial state.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript applies principal component analysis to the joint covariance matrix constructed from binned normalized p_T spectra, event-wise mean p_T, and v_2^2 in heavy-ion collisions. It reports that the first two eigenvectors capture 99.5% of the total variance; these are then subjected to an orthogonal rotation whose angle is determined by two physical constraints motivated by initial-state thermal (isotropic) and geometric (anisotropic) responses. The resulting rotated modes are interpreted as thermal and geometric normal modes analogous to the vibrational modes of a linear triatomic molecule. The thermal mode is stated to fully account for the measured v_0(p_T), while the geometric mode contributes substantially to v_{02}(p_T) in non-central collisions and explains its characteristic low-p_T sign change.

Significance. If the rotation constraints can be shown to be independent of the analyzed covariance data, the work would provide the first physically motivated interpretation of PCA modes in heavy-ion collisions. It would establish a direct link between spectral fluctuations and the thermo-geometric structure of the initial state, furnish a transparent explanation for observed features in v_0 and v_{02}, and introduce a useful molecular-vibration analogy that could guide further modeling and experimental tests.

major comments (1)
  1. [Section on orthogonal rotation of PCA modes] The section describing the orthogonal rotation (following the PCA decomposition): the two physical constraints used to fix the rotation angle are motivated by expected thermal and geometric initial-state responses, yet the manuscript does not demonstrate that these constraints can be derived from independent first-principles initial-state models without reference to the same covariance structure. Consequently, it remains possible that an alternative rotation angle or constraint set would produce equally valid but differently labeled modes, undermining the uniqueness of the thermal/geometric identification.
minor comments (1)
  1. [Abstract and results discussion] The analogy to the vibrational normal modes of a linear triatomic molecule is stated but would be strengthened by an explicit mapping (e.g., a figure or table) showing the correspondence between the rotated eigenvectors and the molecular degrees of freedom.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the major comment below.

read point-by-point responses
  1. Referee: The section describing the orthogonal rotation (following the PCA decomposition): the two physical constraints used to fix the rotation angle are motivated by expected thermal and geometric initial-state responses, yet the manuscript does not demonstrate that these constraints can be derived from independent first-principles initial-state models without reference to the same covariance structure. Consequently, it remains possible that an alternative rotation angle or constraint set would produce equally valid but differently labeled modes, undermining the uniqueness of the thermal/geometric identification.

    Authors: We thank the referee for this observation. The rotation angle is fixed by two constraints that encode general physical expectations: the thermal mode must be isotropic (vanishing elliptic anisotropy by symmetry, as it arises from uniform temperature fluctuations without preferred direction), while the geometric mode encodes the anisotropic shape response. These conditions follow from standard hydrodynamic response theory and initial-state symmetry considerations, without any dependence on the numerical entries or eigenvectors of the covariance matrix. They are therefore independent of the analyzed data structure. We will revise the relevant section to explicitly demonstrate this independence via symmetry arguments and consistency with hydrodynamic expectations, thereby confirming that alternative rotations would violate these physical requirements and that the thermal/geometric identification is unique under the stated constraints. revision: yes

Circularity Check

0 steps flagged

No significant circularity: PCA decomposition and physical interpretation remain independent.

full rationale

The paper applies standard PCA to the joint covariance matrix of normalized p_T spectra, event-wise mean p_T, and v_2^2, directly yielding the leading eigenvectors that capture 99.5% variance as a data-driven result. The subsequent orthogonal rotation is described as imposed by physical constraints motivated by independent initial-state thermal and geometric responses (e.g., isotropic vs. anisotropic fluctuations), not derived from or fitted to the same covariance structure. No equation reduces the rotated modes to the input eigenvectors by construction, nor does any self-citation chain or ansatz smuggle in the target interpretation. The assignment of thermal mode to v_0(p_T) and geometric mode to v_{02}(p_T) features is presented as an interpretive outcome supported by the rotated basis, not a tautological re-labeling of fitted quantities. The derivation chain is therefore self-contained against external initial-state benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract; no explicit free parameters, invented entities, or non-standard axioms are identifiable. The physical constraints for mode rotation constitute an unstated modeling choice whose independence cannot be assessed without the full text.

axioms (1)
  • standard math Principal component analysis applied to the joint covariance matrix of normalized transverse momentum spectrum, mean pT, and elliptic flow squared yields interpretable modes after rotation.
    Standard linear algebra technique; invoked implicitly in the abstract description of the decomposition.

pith-pipeline@v0.9.0 · 5487 in / 1453 out tokens · 77235 ms · 2026-05-07T11:42:02.636599+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

68 extracted references · 57 canonical work pages · 3 internal anchors

  1. [1]

    on the other hand is primarily driven by fluctuations of event-by-event participant eccentricity squared (ϵ 2

  2. [2]

    of the fireball [31– 37] which is strongly correlated with the collision impact parameter. Pearson correlation coefficient between [p T ] and the elliptic flow squared (v 2 2), called the Bo˙ zek’s cor- relator [38], can therefore serve as an excellent probe of ∗ rupam.samanta@ifj.edu.pl the thermo-geometric coupling in the initial state at fixed initial ...

  3. [3]

    Thermal and geometric normal modes of spectral fluctuations in heavy-ion collisions

    As shown in reference [47],v 02 can be decomposed into two components: one part cor- responding to the fluctuations ofN(p T ) associated with [pT ] fluctuations at fixedv 2 2, while the other part cap- tures the fluctuation of spectra associated withv 2 2 at fixed [p T ]. These two components can be regarded as response coefficients or modes, and could be...

  4. [4]

    The reason for such choice of augmented structure is to connect the principal modes obtained from PCA to the experimen- tally accessible observablesv 0(pT ) andv 02(pT ). Without augmenting, by performing a PCA on the spectral covari- ance alone, the physical modes can still be approximately obtained, but no connection to the experimentally mea- surable o...

  5. [5]

    To gain better in- sights about these physical modes we draw an analogy to thenormalmodes of a coupled oscillator, in particular with the case of a linear triatomic molecule. Further- more, we show that thermal and geometric modes can be experimentally accessed via the observablesv 0(pT ) and v02(pT ) by quantifying the mixing coefficients based on our an...

  6. [6]

    The events at fixed entropy density are shown in dark red andrdenotes the Pearson correlation coefficient between [pT ] andϵ 2 2 for those events correlations [48–50, 53–56]

    (b) Similar scat- ter plot between mean transverse momentum and eccentric- ity squared. The events at fixed entropy density are shown in dark red andrdenotes the Pearson correlation coefficient between [pT ] andϵ 2 2 for those events correlations [48–50, 53–56]. In reference [49], the princi- pal components of the flow fluctuations were related to the ini...

  7. [7]

    We choose normalizations in such a way that theN(p T )−[p T ] andN(p T )−v 2 2 blocks match the precise definitions ofv0(pT ) andv 02(pT ) as used in [22, 47]. In particular, we normalize the pT -spectrum and elliptic flow squared by their event- averages, whereas the mean transverse momentum is nor- malized by its standard deviation over events. With thi...

  8. [8]

    (e) Thermal normal mode of spectral fluctuation in Pb+Pb collision at 5.02 TeV for 30-40 % centrality, plotted as a function ofp T having a single node. (f) Geometric normal mode having characteristic double sign change in same centrality class 5 original modes such that the event-by-event coherent changes in the spectra is maximized, hence called the max...

  9. [9]

    Such rotation is expected to play most consequential role in cen- tral collisions where the mixing is substantial

    The effect 4 It should be noted that such orthogonal rotation is a general requirement according to the rigorous relationship between final and initial state fluctuations shown in Appendix A, irrespective of the mixing between the physical modes into principal modes. Such rotation is expected to play most consequential role in cen- tral collisions where t...

  10. [10]

    Hydrodynamic Models for Heavy Ion Collisions,

    P. Huovinen and P. V. Ruuskanen, Hydrodynamic Mod- els for Heavy Ion Collisions, Ann. Rev. Nucl. Part. Sci. 56, 163 (2006), arXiv:nucl-th/0605008

  11. [11]

    Ollitrault, Relativistic hydrodynamics for heavy- ion collisions, European Journal of Physics29, 275 (2008), arXiv:0708.2433 [nucl-th]

    J.-Y. Ollitrault, Relativistic hydrodynamics for heavy-ion collisions, Eur. J. Phys.29, 275 (2008), arXiv:0708.2433 [nucl-th]

  12. [12]

    Muller, J

    B. Muller, J. Schukraft, and B. Wyslouch, First Results from Pb+Pb collisions at the LHC, Ann. Rev. Nucl. Part. Sci.62, 361 (2012), arXiv:1202.3233 [hep-ex]

  13. [13]

    Collective flow and viscosity in relativistic heavy-ion collisions

    U. Heinz and R. Snellings, Collective flow and viscosity in relativistic heavy-ion collisions, Ann. Rev. Nucl. Part. Sci.63, 123 (2013), arXiv:1301.2826 [nucl-th]

  14. [14]

    Heavy Ion Collisions: The Big Picture, and the Big Questions

    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, 339 (2018), arXiv:1802.04801 [hep-ph]

  15. [15]

    Broniowski, M

    W. Broniowski, M. Chojnacki, and L. Obara, Size fluctuations of the initial source and the event-by- event transverse momentum fluctuations in relativistic heavy-ion collisions, Phys. Rev. C80, 051902 (2009), arXiv:0907.3216 [nucl-th]

  16. [16]

    Bozek and W

    P. Bozek and W. Broniowski, Transverse-momentum fluctuations in relativistic heavy-ion collisions from event-by-event viscous hydrodynamics, Phys. Rev. C85, 044910 (2012), arXiv:1203.1810 [nucl-th]

  17. [17]

    Bo˙ zek and W

    P. Bo˙ zek and W. Broniowski, Transverse momentum fluc- tuations in ultrarelativistic Pb + Pb and p + Pb colli- sions with “wounded” quarks, Phys. Rev. C96, 014904 (2017), arXiv:1701.09105 [nucl-th]

  18. [18]

    F. G. Gardim, G. Giacalone, M. Luzum, and J.-Y. Olli- trault, Thermodynamics of hot strong-interaction matter from ultrarelativistic nuclear collisions, Nature Phys.16, 615 (2020), arXiv:1908.09728 [nucl-th]

  19. [19]

    Samanta, S

    R. Samanta, S. Bhatta, J. Jia, M. Luzum, and J.-Y. Ol- litrault, Thermalization at the femtoscale seen in high- energy Pb+Pb collisions, Phys. Rev. C109, L051902 (2024), arXiv:2303.15323 [nucl-th]

  20. [20]

    Ollitrault, Anisotropy as a signature of transverse collective flow, Phys

    J.-Y. Ollitrault, Anisotropy as a signature of transverse collective flow, Phys. Rev. D46, 229 (1992)

  21. [21]

    Voloshin and Y

    S. Voloshin and Y. Zhang, Flow study in relativistic nu- clear collisions by Fourier expansion of Azimuthal parti- cle distributions, Z. Phys. C70, 665 (1996), arXiv:hep- ph/9407282

  22. [22]

    Sorge, Highly sensitive centrality dependence of el- liptic flow: A novel signature of the phase transition in QCD, Phys

    H. Sorge, Highly sensitive centrality dependence of el- liptic flow: A novel signature of the phase transition in QCD, Phys. Rev. Lett.82, 2048 (1999), arXiv:nucl- th/9812057

  23. [23]

    Ollitrault, Measures of azimuthal anisotropy in high-energy collisions, Eur

    J.-Y. Ollitrault, Measures of azimuthal anisotropy in high-energy collisions, Eur. Phys. J. A59, 236 (2023), arXiv:2308.11674 [nucl-ex]

  24. [24]

    Adams et al

    J. Adams et al. (STAR), Azimuthal anisotropy in Au+Au collisions at s(NN)**(1/2) = 200-GeV, Phys. Rev. C72, 014904 (2005), arXiv:nucl-ex/0409033

  25. [25]

    Aamodt et al

    K. Aamodt et al. (ALICE), Elliptic flow of charged par- ticles in Pb-Pb collisions at 2.76 TeV, Phys. Rev. Lett. 105, 252302 (2010), arXiv:1011.3914 [nucl-ex]

  26. [26]

    Aad et al

    G. Aad et al. (ATLAS), Measurement of the azimuthal anisotropy for charged particle production in √sN N = 2.76 TeV lead-lead collisions with the ATLAS detector, Phys. Rev. C86, 014907 (2012), arXiv:1203.3087 [hep- ex]

  27. [27]

    F. G. Gardim, G. Giacalone, M. Luzum, and J.-Y. Ollitrault, Effects of initial state fluctuations on the mean transverse momentum, Nucl. Phys. A1005, 121999 (2021), arXiv:2002.07008 [nucl-th]

  28. [28]

    Samanta, J

    R. Samanta, J. P. Picchetti, M. Luzum, and J.-Y. Olli- trault, Non-Gaussian transverse momentum fluctuations from impact parameter fluctuations, Phys. Rev. C108, 024908 (2023), arXiv:2306.09294 [nucl-th]

  29. [29]

    Alqahtani, T

    M. Alqahtani, T. Parida, and J.-Y. Ollitrault, Extracting the speed of sound of QCD from transverse momentum fluctuations, (2026), arXiv:2603.09647 [hep-ph]

  30. [30]

    Schenke, C

    B. Schenke, C. Shen, and D. Teaney, Transverse momen- tum fluctuations and their correlation with elliptic flow in nuclear collision, Phys. Rev. C102, 034905 (2020), arXiv:2004.00690 [nucl-th]

  31. [31]

    Parida, R

    T. Parida, R. Samanta, and J.-Y. Ollitrault, Probing collectivity in heavy-ion collisions with fluctuations of the pT spectrum, Phys. Lett. B857, 138985 (2024), arXiv:2407.17313 [nucl-th]

  32. [32]

    The shape of differential radial flow $v_0(p_T)$, not its zero-crossing, carries physical information

    S. Bhatta, A. Dimri, and J. Jia, Disentangling the global multiplicity and spectral shape fluctuations in radial flow, (2025), arXiv:2504.20008 [nucl-th]

  33. [33]

    S. Saha, R. Singh, and B. Mohanty, pT-differential radial flow in a blast-wave model, Phys. Rev. C112, 024902 (2025), arXiv:2505.19697 [nucl-ex]

  34. [34]

    Jia, Sources of Radial-Flow Fluctuations in the Quark- Gluon Plasma, Phys

    J. Jia, Sources of Radial-Flow Fluctuations in the Quark- Gluon Plasma, Phys. Rev. Lett.136, 112301 (2026), arXiv:2507.14399 [nucl-th]

  35. [35]

    Du, Phys

    L. Du, Characterizing radial flow fluctuations in rela- tivistic heavy-ion collisions at top RHIC and LHC ener- gies, Phys. Rev. C113, 014901 (2026), arXiv:2508.07184 [nucl-ex]

  36. [36]

    Agarwala, D

    R. Agarwala, D. Basak, and K. Dey, Observation of par- tonic collectivity viap T-differential radial flow fluctua- tions in Au+Au collisions at √sNN = 200 GeV, (2025), arXiv:2512.15026 [nucl-ex]

  37. [37]

    Rapidity dependence of mean transverse momentum fluctuation and decorrelation in baryon-dense medium

    T. Parida, Rapidity dependence of mean transverse mo- mentum fluctuation and decorrelation in baryon-dense medium, (2026), arXiv:2602.16369 [nucl-th]

  38. [38]

    Aadet al.(ATLAS), Phys

    G. Aad et al. (ATLAS), Evidence for the Collective Na- ture of Radial Flow in Pb+Pb Collisions with the AT- LAS Detector, Phys. Rev. Lett.136, 032301 (2026), arXiv:2503.24125 [nucl-ex]

  39. [39]

    Acharyaet al.(ALICE), Phys

    S. Acharya et al. (ALICE), Long-range transverse mo- mentum correlations and radial flow in Pb−Pb colli- sions at the LHC, Phys. Rev. Lett.136, 032302 (2026), arXiv:2504.04796 [nucl-ex]

  40. [40]

    R. S. Bhalerao and J.-Y. Ollitrault, Eccentricity fluctua- tions and elliptic flow at RHIC, Phys. Lett. B641, 260 (2006), arXiv:nucl-th/0607009

  41. [41]

    Broniowski, P

    W. Broniowski, P. Bozek, and M. Rybczynski, Fluctu- ating initial conditions in heavy-ion collisions from the Glauber approach, Phys. Rev. C76, 054905 (2007), arXiv:0706.4266 [nucl-th]

  42. [42]

    Alver and G

    B. Alver and G. Roland, Collision geometry fluctuations and triangular flow in heavy-ion collisions, Phys. Rev. C81, 054905 (2010), [Erratum: Phys.Rev.C 82, 039903 (2010)], arXiv:1003.0194 [nucl-th]

  43. [43]

    Holopainen, H

    H. Holopainen, H. Niemi, and K. J. Eskola, Event- by-event hydrodynamics and elliptic flow from fluctu- ating initial state, Phys. Rev. C83, 034901 (2011), 8 arXiv:1007.0368 [hep-ph]

  44. [44]

    G.-Y. Qin, H. Petersen, S. A. Bass, and B. Muller, Trans- lation of collision geometry fluctuations into momen- tum anisotropies in relativistic heavy-ion collisions, Phys. Rev. C82, 064903 (2010), arXiv:1009.1847 [nucl-th]

  45. [45]

    R. S. Bhalerao, M. Luzum, and J.-Y. Ollitrault, Deter- mining initial-state fluctuations from flow measurements in heavy-ion collisions, Phys. Rev. C84, 034910 (2011), arXiv:1104.4740 [nucl-th]

  46. [46]

    Initial State Fluctuations and Final State Correlations in Relativistic Heavy-Ion Collisions

    M. Luzum and H. Petersen, Initial State Fluctuations and Final State Correlations in Relativistic Heavy-Ion Collisions, J. Phys. G41, 063102 (2014), arXiv:1312.5503 [nucl-th]

  47. [47]

    Bozek, Transverse-momentum–flow correlations in rel- ativistic heavy-ion collisions, Phys

    P. Bozek, Transverse-momentum–flow correlations in rel- ativistic heavy-ion collisions, Phys. Rev. C93, 044908 (2016), arXiv:1601.04513 [nucl-th]

  48. [48]

    Bozek and H

    P. Bozek and H. Mehrabpour, Correlation coefficient between harmonic flow and transverse momentum in heavy-ion collisions, Phys. Rev. C101, 064902 (2020), arXiv:2002.08832 [nucl-th]

  49. [49]

    Correlation between mean transverse momen- tum and anisotropic flow in heavy-ion collisions

    G. Giacalone, B. Schenke, and C. Shen, Observable sig- natures of initial state momentum anisotropies in nu- clear collisions, Phys. Rev. Lett.125, 192301 (2020), arXiv:2006.15721 [nucl-th]

  50. [50]

    Giacalone, F

    G. Giacalone, F. G. Gardim, J. Noronha-Hostler, and J.-Y. Ollitrault, Correlation between mean transverse momentum and anisotropic flow in heavy-ion collisions, Phys. Rev. C103, 024909 (2021), arXiv:2004.01765 [nucl-th]

  51. [51]

    Bozek and R

    P. Bozek and R. Samanta, Higher order cumulants of transverse momentum and harmonic flow in relativistic heavy ion collisions, Phys. Rev. C104, 014905 (2021), arXiv:2103.15338 [nucl-th]

  52. [52]

    Samanta and P

    R. Samanta and P. Bozek, Momentum-dependent mea- sures of correlations between mean transverse momentum and harmonic flow in heavy-ion collisions, Phys. Rev. C 109, 064910 (2024), arXiv:2308.11565 [nucl-th]

  53. [53]

    Acharya et al

    S. Acharya et al. (ALICE), Characterizing the initial con- ditions of heavy-ion collisions at the LHC with mean transverse momentum and anisotropic flow correlations, Phys. Lett. B834, 137393 (2022), arXiv:2111.06106 [nucl-ex]

  54. [54]

    Aad et al

    G. Aad et al. (ATLAS), Correlations between flow and transverse momentum in Xe+Xe and Pb+Pb collisions at the LHC with the ATLAS detector: A probe of the heavy-ion initial state and nuclear deformation, Phys. Rev. C107, 054910 (2023), arXiv:2205.00039 [nucl-ex]

  55. [55]

    M. I. Abdulhamid et al. (STAR), Imaging shapes of atomic nuclei in high-energy nuclear collisions, Nature 635, 67 (2024), arXiv:2401.06625 [nucl-ex]

  56. [56]

    Parida, R

    T. Parida, R. Samanta, and J.-Y. Ollitrault, Correlation between particle spectra and elliptic flow, Phys. Lett. B 868, 139729 (2025), arXiv:2506.18690 [nucl-th]

  57. [57]

    R. S. Bhalerao, J.-Y. Ollitrault, S. Pal, and D. Teaney, Principal component analysis of event-by-event fluc- tuations, Phys. Rev. Lett.114, 152301 (2015), arXiv:1410.7739 [nucl-th]

  58. [58]

    Mazeliauskas and D

    A. Mazeliauskas and D. Teaney, Fluctuations of har- monic and radial flow in heavy ion collisions with prin- cipal components, Phys. Rev. C93, 024913 (2016), arXiv:1509.07492 [nucl-th]

  59. [59]

    A. M. Sirunyan et al. (CMS), Principal-component anal- ysis of two-particle azimuthal correlations in PbPb and pPb collisions at CMS, Phys. Rev. C96, 064902 (2017), arXiv:1708.07113 [nucl-ex]

  60. [60]

    M. B. Richman, Rotation of principal components, Jour- nal of Climatology6, 293

  61. [61]

    I. T. Jolliffe, Rotation of ill-defined principal components, Journal of the Royal Statistical Society. Series C (Applied Statistics)38, 139 (1989)

  62. [62]

    Mazeliauskas and D

    A. Mazeliauskas and D. Teaney, Subleading harmonic flows in hydrodynamic simulations of heavy ion collisions, Phys. Rev. C91, 044902 (2015), arXiv:1501.03138 [nucl- th]

  63. [63]

    Bo˙ zek, Angle and magnitude decorrelation in the fac- torization breaking of collective flow, Phys

    P. Bo˙ zek, Angle and magnitude decorrelation in the fac- torization breaking of collective flow, Phys. Rev. C98, 064906 (2018), arXiv:1808.04248 [nucl-th]

  64. [64]

    Z. Liu, W. Zhao, and H. Song, Principal Component Analysis of collective flow in Relativistic Heavy-Ion Col- lisions, Eur. Phys. J. C79, 870 (2019), arXiv:1903.09833 [nucl-th]

  65. [65]

    Z. Liu, A. Behera, H. Song, and J. Jia, Robustness of principal component analysis of harmonic flow in heavy ion collisions, Phys. Rev. C102, 024911 (2020), arXiv:2002.06061 [nucl-ex]

  66. [66]

    J. S. Moreland, J. E. Bernhard, and S. A. Bass, Alter- native ansatz to wounded nucleon and binary collision scaling in high-energy nuclear collisions, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708 [nucl-th]

  67. [67]

    Schenke, S

    B. Schenke, S. Jeon, and C. Gale, (3+1)D hydrodynamic simulation of relativistic heavy-ion collisions, Phys. Rev. C82, 014903 (2010), arXiv:1004.1408 [hep-ph]

  68. [68]

    Bozek and R

    P. Bozek and R. Samanta, Factorization breaking for higher moments of harmonic flow, Phys. Rev. C105, 034904 (2022), arXiv:2109.07781 [nucl-th]. Appendix A: Relation between final and initial state In our analysis the event-by-event fluctuation of final state observables can be written in matrix form: δX= δN(p T1). . . δN(p Ti). . . δp T δv 2 2 (A1) where...