Processing math: 100%
CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIN-24-004
Characterizing the Initial State in XeXe and PbPb Collisions using Multiparticle Cumulants
Abstract: The nuclear shape and size dependences of Fourier flow harmonics (vn, n= 2, 3, 4) in XeXe collisions at sNN= 5.44 TeV and PbPb collisions at sNN= 5.36 TeV are studied using the CMS detector. For the first time, correlations between higher-power moments of two or three flow harmonics, as a function of collision centrality, are compared in XeXe and PbPb collisions. This is achieved by measuring multiparticle mixed harmonic cumulants (up to 8-particle cumulants) using charged particles in the pseudorapidity range |η|< 2.4 and transverse momentum range 0.5 <pT< 3.0 GeV/c. The results are compared to theoretical models to investigate the non-linear hydrodynamic response by examining correlations between v2, v3, and v4 and their corresponding eccentricities ε2, ε3, and ε4. This work provides valuable insights into the deformation parameters of the Xe nucleus and constrains initial-state model parameters influencing the evolution of the quark-gluon plasma in heavy-ion collisions at the LHC.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from the 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red solid lines show the prediction from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and the black line is for spherical Xe (Xe-(b)) when taken in ratio with Pb. The shaded bands show the different hydrodynamic predictions for different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png
Figure 1-a:
vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from the 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red solid lines show the prediction from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and the black line is for spherical Xe (Xe-(b)) when taken in ratio with Pb. The shaded bands show the different hydrodynamic predictions for different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png
Figure 1-b:
vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from the 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red solid lines show the prediction from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and the black line is for spherical Xe (Xe-(b)) when taken in ratio with Pb. The shaded bands show the different hydrodynamic predictions for different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png pdf
Figure 2:
vn{4}/vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model parameters for Xe-(a). The cyan line shows the theoretical prediction for v3{4}/v3{2,|Δη|>2} according to A1/4 scaling.

png
Figure 2-a:
vn{4}/vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model parameters for Xe-(a). The cyan line shows the theoretical prediction for v3{4}/v3{2,|Δη|>2} according to A1/4 scaling.

png
Figure 2-b:
vn{4}/vn{2,|Δη|>2} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of the initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model parameters for Xe-(a). The cyan line shows the theoretical prediction for v3{4}/v3{2,|Δη|>2} according to A1/4 scaling.

png pdf
Figure 3:
The top panel shows the normalized symmetric cumulants NSC(2,3) and NSC(2,4) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines show the prediction from IP-Glasma IC model having the same Xe and Pb nucleus parameters, while the shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model. The bottom panel shows the ratios of NSC(m,n) for XeXe/PbPb as a function of centrality. The red solid lines show the initial-state predictions from TRENTo-IC while the orange solid lines show those from IP-Glasma IC. The blue shaded band is the hydrodynamic prediction given by IP-Glasma+MUSIC+UrQMD.

png
Figure 3-a:
The top panel shows the normalized symmetric cumulants NSC(2,3) and NSC(2,4) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines show the prediction from IP-Glasma IC model having the same Xe and Pb nucleus parameters, while the shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model. The bottom panel shows the ratios of NSC(m,n) for XeXe/PbPb as a function of centrality. The red solid lines show the initial-state predictions from TRENTo-IC while the orange solid lines show those from IP-Glasma IC. The blue shaded band is the hydrodynamic prediction given by IP-Glasma+MUSIC+UrQMD.

png
Figure 3-b:
The top panel shows the normalized symmetric cumulants NSC(2,3) and NSC(2,4) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red and blue solid lines show the predictions from TRENTo-IC in terms of initial-state eccentricities εn for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines show the prediction from IP-Glasma IC model having the same Xe and Pb nucleus parameters, while the shaded band shows the hydrodynamic predictions for the IP-Glasma+MUSIC+UrQMD model. The bottom panel shows the ratios of NSC(m,n) for XeXe/PbPb as a function of centrality. The red solid lines show the initial-state predictions from TRENTo-IC while the orange solid lines show those from IP-Glasma IC. The blue shaded band is the hydrodynamic prediction given by IP-Glasma+MUSIC+UrQMD.

png pdf
Figure 4:
A selection of six-particle three-harmonic symmetric cumulants and normalized symmetric cumulants as a function of centrality in XeXe (left plot) and PbPb (right plot), calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png
Figure 4-a:
A selection of six-particle three-harmonic symmetric cumulants and normalized symmetric cumulants as a function of centrality in XeXe (left plot) and PbPb (right plot), calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png
Figure 4-b:
A selection of six-particle three-harmonic symmetric cumulants and normalized symmetric cumulants as a function of centrality in XeXe (left plot) and PbPb (right plot), calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png pdf
Figure 5:
Six-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png
Figure 5-a:
Six-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png
Figure 5-b:
Six-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei.

png pdf
Figure 6:
Eight-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines represent the initial-state model predictions.

png
Figure 6-a:
Eight-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines represent the initial-state model predictions.

png
Figure 6-b:
Eight-particle normalized mixed harmonic cumulants (nMHC) as a function of centrality in XeXe and PbPb, calculated from multiparticle cumulants. The bars represent the statistical and negligible systematic uncertainties, added in quadrature. The shaded bands show the corresponding hydrodynamic predictions from the IP-Glasma+MUSIC+UrQMD model for deformed Xe (Xe-(a)), and Pb nuclei. The solid lines represent the initial-state model predictions.

png pdf
Figure 7:
v4{2,|Δη|>2} as a function of centrality in XeXe and PbPb, calculated from 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red line shows the prediction from TRENTo-IC in terms in terms of ε4 for Xe-(a), while the shaded bands show the different hydrodynamic predictions for four different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png
Figure 7-a:
v4{2,|Δη|>2} as a function of centrality in XeXe and PbPb, calculated from 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red line shows the prediction from TRENTo-IC in terms in terms of ε4 for Xe-(a), while the shaded bands show the different hydrodynamic predictions for four different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png
Figure 7-b:
v4{2,|Δη|>2} as a function of centrality in XeXe and PbPb, calculated from 2-particle correlation method. The bars and the open boxes represent statistical and systematic uncertainties, respectively. The red line shows the prediction from TRENTo-IC in terms in terms of ε4 for Xe-(a), while the shaded bands show the different hydrodynamic predictions for four different sets of deformations from the IP-Glasma+MUSIC+UrQMD model.

png pdf
Figure 8:
vn{4} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated using 4-particle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively.

png
Figure 8-a:
vn{4} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated using 4-particle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively.

png
Figure 8-b:
vn{4} (n = 2, 3) as a function of centrality in XeXe and PbPb, calculated using 4-particle cumulants. The bars and the open boxes represent statistical and systematic uncertainties, respectively.
Tables

png pdf
Table 1:
Woods-Saxon deformation parameters for the four sets of Xe nucleus and the Pb nucleus.

png pdf
Table 2:
Signs of correlations between different moments of v2 and v3/v4.
Summary
For the first time, detailed comparisons of individual flow harmonics, two- and three-harmonic correlations, and higher-power mixed harmonic cumulants have been performed between XeXe collisions at sNN= 5.44 TeV and PbPb collisions at sNN= 5.36 TeV. The results have been compared with the theoretical predictions of the IP-Glasma+MUSIC+UrQMD hydrodynamic model and TRENTo-IC and IP-Glasma initial-state models, which were calculated with various sets of deformation parameters for the Xe nucleus. We find that the final-state IP-Glasma+MUSIC+UrQMD model calculation with R0= 5.601, a0= 0.492, β2= 0.207, β4=0.003 for Xe matches the data best. On the other hand, the difference in initial-state predictions from the TRENTo-IC and IP-Glasma IC models not only point to the sensitivity of these observables to the pre-equilibrium dynamics needed to model the experimental data, but also to the increasing non-linearity of the higher-power moments of the flow harmonics in peripheral collisions. The comparison of the higher-power moments in XeXe and PbPb collisions is a good way to study the impact of the system size on these observables, where it is seen that the observables containing higher-power moments of fluctuation-driven v3 and v4 have greater magnitude in XeXe, it being the smaller nucleus. Subtle differences exist between data and model calculations with the current sets of parameters. These differences highlight the importance of fine-tuning the model parameters, including not only the nuclear deformation parameters but also those intrinsic to hydrodynamic calculations such as transport coefficients and freeze-out criteria. Moreover, the degree of applicability of hydrodynamics to peripheral heavy-ion collisions remains an open question. Such analyses can significantly improve the overall understanding of initial conditions and transport properties of the Quark-Gluon Plasma (QGP) created at the LHC.
References
1 BRAHMS Collaboration Quark gluon plasma and color glass condensate at RHIC? The Perspective from the BRAHMS experiment Nucl. Phys. A 757 (2005) 1 nucl-ex/0410020
2 J.-Y. Ollitrault Anisotropy as a signature of transverse collective flow PRD 46 (1992) 229
3 C. Mordasini, A. Bilandzic, D. Karako ç , and S. F. Taghavi Higher order Symmetric Cumulants Phys. Rev. C 102 (2020) 024907 1901.06968
4 ALICE Collaboration Investigations of Anisotropic Flow Using Multiparticle Azimuthal Correlations in pp, p-Pb, Xe-Xe, and Pb-Pb Collisions at the LHC PRL 123 (2019) 142301 1903.01790
5 CMS Collaboration Observation of Correlated Azimuthal Anisotropy Fourier Harmonics in pp and p+Pb Collisions at the LHC PRL 120 (2018) 092301 CMS-HIN-16-022
1709.09189
6 J. Noronha-Hostler, L. Yan, F. G. Gardim, and J.-Y. Ollitrault Linear and cubic response to the initial eccentricity in heavy-ion collisions Phys. Rev. C 93 (2016) 014909 1511.03896
7 ALICE Collaboration Measurements of mixed harmonic cumulants in Pb-Pb collisions at sNN= 5.02 TeV PLB 818 (2021) 136354 2102.12180
8 CMS Collaboration Mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles in PbPb collisions at sNN= 2.76 and 5.02 TeV EPJC 80 (2020) 534 CMS-HIN-17-005
1910.08789
9 F. G. Gardim, F. Grassi, M. Luzum, and J.-Y. Ollitrault Breaking of factorization of two-particle correlations in hydrodynamics Phys. Rev. C 87 (2013) 031901 1211.0989
10 I. Kozlov et al. Signatures of collective behavior in small systems Nucl. Phys. A 931 (2014) 1045 1412.3147
11 ATLAS Collaboration Measurement of the correlation between flow harmonics of different order in lead-lead collisions at sNN= 2.76 TeV with the ATLAS detector Phys. Rev. C 92 (2015) 034903 1504.01289
12 ATLAS Collaboration Measurement of event-plane correlations in sNN= 2.76 TeV lead-lead collisions with the ATLAS detector Phys. Rev. C 90 (2014) 024905 1403.0489
13 ATLAS Collaboration Correlated long-range mixed-harmonic fluctuations measured in pp, p+Pb and low-multiplicity Pb+Pb collisions with the ATLAS detector PLB 789 (2019) 444 1807.02012
14 PHENIX Collaboration Systematic Study of Azimuthal Anisotropy in Cu+Cu and Au+Au Collisions at sNN= 62.4 and 200 GeV Phys. Rev. C 92 (2015) 034913 1412.1043
15 CMS Collaboration Description and performance of track and primary-vertex reconstruction with the CMS tracker JINST 9 (2014) P10009 CMS-TRK-11-001
1405.6569
16 CMS Collaboration The CMS Experiment at the CERN LHC JINST 3 (2008) S08004
17 GEANT4 Collaboration GEANT4--a simulation toolkit NIM A 506 (2003) 250
18 CMS Collaboration Charged-particle angular correlations in XeXe collisions at sNN= 5.44 TeV Phys. Rev. C 100 (2019) 044902 CMS-HIN-18-001
1901.07997
19 CMS Collaboration Azimuthal Anisotropy of Charged Particles at High Transverse Momenta in PbPb Collisions at sNN= 2.76 TeV PRL 109 (2012) 022301 CMS-HIN-11-012
1204.1850
20 CMS Collaboration Measurement of the Pseudorapidity and Centrality Dependence of the Transverse Energy Density in PbPb Collisions at sNN= 2.76 TeV PRL 109 (2012) 152303 CMS-HIN-11-003
1205.2488
21 A. Bilandzic, R. Snellings, and S. Voloshin Flow analysis with cumulants: Direct calculations Phys. Rev. C 83 (2011) 044913 1010.0233
22 CMS Collaboration Multiparticle correlation studies in pPb collisions at sNN= 8.16 TeV Phys. Rev. C 101 (2020) 014912 CMS-HIN-17-004
1904.11519
23 A. Bilandzic et al. Generic framework for anisotropic flow analyses with multiparticle azimuthal correlations Phys. Rev. C 89 (2014) 064904 1312.3572
24 CMS Collaboration Correlations of azimuthal anisotropy Fourier harmonics with subevent cumulants in pPb collisions at sNN= 8.16 TeV Phys. Rev. C 103 (2021) 014902 CMS-HIN-18-015
1905.09935
25 ALICE Collaboration Correlated event-by-event fluctuations of flow harmonics in Pb-Pb collisions at sNN= 2.76 TeV PRL 117 (2016) 182301 1604.07663
26 Z. Moravcova, K. Gulbrandsen, and Y. Zhou Generic algorithm for multiparticle cumulants of azimuthal correlations in high energy nucleus collisions Phys. Rev. C 103 (2021) 024913 2005.07974
27 B. Schenke, P. Tribedy, and R. Venugopalan Fluctuating Glasma initial conditions and flow in heavy ion collisions PRL 108 (2012) 252301 1202.6646
28 S. McDonald et al. A Detailed Study and Synthesis of Flow Observables in the IP-Glasma+MUSIC+UrQMD Framework Nucl. Phys. A 967 (2017) 393 1704.05362
29 C. Gale, S. Jeon, and B. Schenke Hydrodynamic Modeling of Heavy-Ion Collisions Int. J. Mod. Phys. A 28 (2013) 1340011 1301.5893
30 S. A. Bass et al. Microscopic models for ultrarelativistic heavy ion collisions Prog. Part. Nucl. Phys. 41 (1998) 255 nucl-th/9803035
31 J. S. Moreland, J. E. Bernhard, and S. A. Bass Alternative ansatz to wounded nucleon and binary collision scaling in high-energy nuclear collisions Phys. Rev. C 92 (2015) 011901 1412.4708
32 G. Giacalone, J. Noronha-Hostler, M. Luzum, and J.-Y. Ollitrault Hydrodynamic predictions for 5.44 TeV Xe+Xe collisions Phys. Rev. C 97 (2018) 034904 1711.08499
33 H. M ä ntysaari, B. Schenke, C. Shen, and W. Zhao Probing nuclear structure of heavy ions at energies available at the CERN Large Hadron Collider Phys. Rev. C 110 (2024) 054913 2409.19064
34 S. Zhao et al. Exploring the Nuclear Shape Phase Transition in Ultra-Relativistic 129Xe+129Xe Collisions at the LHC PRL 133 (2024) 192301 2403.07441
35 ALICE Collaboration Multiharmonic Correlations of Different Flow Amplitudes in Pb-Pb Collisions at sNN= 2.76 TeV PRL 127 (2021) 092302 2101.02579
36 W. Zhao, H.-j. Xu, and H. Song Collective flow in 2.76 A TeV and 5.02 A TeV Pb+Pb collisions EPJC 77 (2017) 645 1703.10792
Compact Muon Solenoid
LHC, CERN