Extracting the speed of sound in the strongly interacting matter created in ultrarelativistic lead-lead collisions at the LHC

The CMS collaboration Hayrapetyan, Aram ; Tumasyan, Armen ; Adam, Wolfgang ; et al.
CMS-HIN-23-003, 2024.
Inspire Record 2747107 DOI 10.17182/hepdata.146016

Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark-gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we extracted the speed of sound in this medium created using lead-lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 $\pm$ 0.002 (stat) $\pm$ 0.016 (syst) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 $\pm$ 8 (syst) MeV. The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

3 data tables

The event fraction distribution as a function of the charged-particle multiplicity, $N_{\mathrm{ch}}$, within the kinematic range of $|\eta|<0.5$ and extrapolated to the full $p_{\mathrm{T}}$ range, in PbPb collisions at $\sqrt{s_{_{\mathrm{NN}}}}=5.02$ TeV. The $N_{\mathrm{ch}}$ value is normalized by its value in the $0-5\%$ centrality class ($N_{\mathrm{ch}}^{0}$).

The average transverse momentum of charged particles, $\langle p_{\mathrm{T}}\rangle$, as a function of the charged-particle multiplicity, $N_{\mathrm{ch}}$, within the kinematic range of $|\eta|<0.5$ and extrapolated to the full $p_{\mathrm{T}}$ range in PbPb collisions at $\sqrt{s_{_{\mathrm{NN}}}}=5.02$ TeV. Both $\langle p_{\mathrm{T}}\rangle$ and $N_{\mathrm{ch}}$ are normalized by their values in the $0-5\%$ centrality class (${\langle p_{\mathrm{T}}\rangle}^{0}$ and $N_{\mathrm{ch}}^{0}$).

The speed of sound, $c^2_{\mathrm{s}}$, as a function of the effective temperature, $T_{\mathrm{eff}}$, with the CMS data point obtained from ultra-central PbPb collision data at $\sqrt{s_{_{\mathrm{NN}}}}=5.02$ TeV.


Higher-order moments of the elliptic flow distribution in PbPb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

The CMS collaboration Tumasyan, Armen ; Adam, Wolfgang ; Andrejkovic, Janik Walter ; et al.
JHEP 2024 (2024) 106, 2024.
Inspire Record 2724506 DOI 10.17182/hepdata.132638

The hydrodynamic flow-like behavior of charged hadrons in high-energy lead-lead collisions is studied through multiparticle correlations. The elliptic anisotropy values based on different orders of multiparticle cumulants, $v_{2}\{2k\}$, are measured up to the tenth order ($k$ = 5) as functions of the collision centrality at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. A hierarchy is observed between the coefficients, with $v_{2}\{2\} > v_{2}\{4\} \gtrsim v_{2}\{6\} \gtrsim v_{2}\{8\} \gtrsim v_{2}\{10\}$. Based on these results, centrality-dependent moments for the fluctuation-driven event-by-event $v_{2}$ distribution are determined, including the skewness, kurtosis and, for the first time, superskewness. Assuming a hydrodynamic expansion of the produced medium, these moments directly probe the initial-state geometry in high-energy nucleus-nucleus collisions.

8 data tables

The $v_2\{2k\}$ ($k = 1, ..., 5$) values as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The relative differences $(v_2\{2k\}-v_2\{10\})/v_2\{10\}$ ($k = 1, ..., 4$) as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The $h_1$ and the $h_2$ hydrodynamic probes as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

More…

Collision-system and beam-energy dependence of anisotropic flow fluctuations

The STAR collaboration Abdallah, Mohamed ; Adam, Jaroslav ; Adamczyk, Leszek ; et al.
Phys.Rev.Lett. 129 (2022) 252301, 2022.
Inspire Record 2017211 DOI 10.17182/hepdata.116554

Elliptic flow measurements from two-, four- and six-particle correlations are used to investigate flow fluctuations in collisions of U+U at $\sqrt{s_{\rm NN}}$= 193 GeV, Cu+Au at $\sqrt{s_{\rm NN}}$= 200 GeV and Au+Au spanning the range $\sqrt{s_{\rm NN}}$= 11.5 - 200 GeV. The measurements show a strong dependence of the flow fluctuations on collision centrality, a modest dependence on system size, and very little if any, dependence on particle species and beam energy. The results, when compared to similar LHC measurements, viscous hydrodynamic calculations, and T$\mathrel{\protect\raisebox{-2.1pt}{R}}$ENTo model eccentricities, indicate that initial-state-driven fluctuations predominate the flow fluctuations generated in the collisions studied.

11 data tables

The Au+Au 200 GeV measurements of the two and four-particle elliptic flow and the elliptic flow fluctuations of the $\pi$ particle.

The Au+Au 200 GeV measurements of the two and four-particle elliptic flow and the elliptic flow fluctuations of the K particle.

The Au+Au 200 GeV measurements of the two and four-particle elliptic flow and the elliptic flow fluctuations of the p particle.

More…

Hyperon polarization along the beam direction relative to the second and third harmonic event planes in isobar collisions at $\sqrt{s_{NN}}$ = 200 GeV

The STAR collaboration Abdulhamid, Muhammad ; Aboona, Bassam ; Adam, Jaroslav ; et al.
Phys.Rev.Lett. 131 (2023) 202301, 2023.
Inspire Record 2643014 DOI 10.17182/hepdata.139081

The polarization of $\Lambda$ and $\bar{\Lambda}$ hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild $p_T$ dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagree with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and $p_T$ dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.

5 data tables

$sgn(\alpha_H)\langle\cos(\theta_p^{\ast})\rangle$ of $\Lambda$ and $\bar{\Lambda}$ as a function of hyperon azimuthal angle relative to the second-order event plane in isobar collisions at 200 GeV.

$sgn(\alpha_H)\langle\cos(\theta_p^{\ast})\rangle$ of $\Lambda$ and $\bar{\Lambda}$ as a function of hyperon azimuthal angle relative to the third-order event plane in isobar collisions at 200 GeV.

$P_z$ sine coefficients of $\Lambda+\bar{\Lambda}$ as a function of centrality in isobar collisions at 200 GeV.

More…

Measurement of electrons from open heavy-flavor hadron decays in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV with the STAR detector

The STAR collaboration Abdulhamid, M.I. ; Aboona, B.E. ; Adam, Jaroslav ; et al.
JHEP 06 (2023) 176, 2023.
Inspire Record 2641480 DOI 10.17182/hepdata.139080

We report a new measurement of the production of electrons from open heavy-flavor hadron decays (HFEs) at mid-rapidity ($|y|<$ 0.7) in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Invariant yields of HFEs are measured for the transverse momentum range of $3.5 < p_{\rm T} < 9$ GeV/$c$ in various configurations of the collision geometry. The HFE yields in head-on Au+Au collisions are suppressed by approximately a factor of 2 compared to that in $p$+$p$ collisions scaled by the average number of binary collisions, indicating strong interactions between heavy quarks and the hot and dense medium created in heavy-ion collisions. Comparison of these results with models provides additional tests of theoretical calculations of heavy quark energy loss in the quark-gluon plasma.

6 data tables

Ratios of NPE (non-photonic electron) to PHE (photonic electron) as a function of $p_{\rm T}$ in 0-10% central (yellow circles) and 40-80% peripheral (green squares) Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Vertical bars represent statistical uncertainties while boxes represent systematic uncertainties. Horizontal bars indicate the bin width.

Invariant yields of electrons from decays of prompt $J/\psi$ (dot-dashed line), $\Upsilon$ (dotted line), Drell-Yan (long dash-dotted line), light vector mesons (long dashed line) and the combined HDE (hadron decayed electron) contribution (solid line), estimated utilizing experimental measurements, theoretical calculations, and PYTHIA and $\rm E_{VT}G_{EN}$ event generators, in 0-10% central Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Color bands represent systematic uncertainties.

Invariant yields of electrons from decays of prompt $J/\psi$ (dot-dashed line), $\Upsilon$ (dotted line), Drell-Yan (long dash-dotted line), light vector mesons (long dashed line) and the combined HDE (hadron decayed electron) contribution (solid line), estimated utilizing experimental measurements, theoretical calculations, and PYTHIA and $\rm E_{VT}G_{EN}$ event generators, in 40-80% central Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Color bands represent systematic uncertainties.

More…

Study of charm hadronization with prompt $\Lambda^+_\mathrm{c}$ baryons in proton-proton and lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

The CMS collaboration Tumasyan, Armen ; Adam, Wolfgang ; Andrejkovic, Janik Walter ; et al.
JHEP 01 (2024) 128, 2024.
Inspire Record 2679262 DOI 10.17182/hepdata.135973

The production of prompt $\Lambda^+_\mathrm{c}$ baryons is measured via the exclusive decay channel $\Lambda^+_\mathrm{c}$$\to$ pK$^-\pi^+$ at a center-of-mass energy per nucleon pair of 5.02 TeV, using proton-proton (pp) and lead-lead (PbPb) collision data collected by the CMS experiment at the CERN LHC. The pp and PbPb data were obtained in 2017 and 2018 with integrated luminosities of 252 and 0.607 nb$^{-1}$, respectively. The measurements are performed within the $\Lambda^+_\mathrm{c}$ rapidity interval $\vert y \vert$$\lt$ 1 with transverse momentum ($p_\mathrm{T}$) ranges of 3-30 and 6-40 GeV/$c$ for pp and PbPb collisions, respectively. Compared to the yields in pp collisions scaled by the expected number of nucleon-nucleon interactions, the observed yields of $\Lambda^+_\mathrm{c}$ with $p_\mathrm{T}$$\gt$ 10 GeV/$c$ are strongly suppressed in PbPb collisions. The level of suppression depends significantly on the collision centrality. The $\Lambda^+_\mathrm{c}$ / D$^0$ production ratio is similar in PbPb and pp collisions at $p_\mathrm{T}$$\gt$ 10 GeV/$c$, suggesting that the coalescence process does not play a dominant role in prompt $\Lambda^+_\mathrm{c}$ baryon production at higher $p_\mathrm{T}$.

6 data tables

The product of acceptance and efficiency ($A \epsilon$) as a function of $p_{\mathrm{T}}$ for prompt $\Lambda^+_c$ in pp collisions and within centrality regions of 0-90, 0-10, 10-30, 30-50 and 50-90% in PbPb collisions.

The $p_{\mathrm{T}}$ differential cross sections for prompt $\Lambda^+_c$ production in pp collisions. The global fit uncertainty is 8.6%.

The $\mathrm{T_{AA}}$-scaled $\Lambda^+_c$ yields as a function of $p_{\mathrm{T}}$ for PbPb collisions within centrality regions of 0-90, 0-10, 10-30, 30-50 and 50-90%.

More…

Version 3
Centrality and transverse momentum dependence of $D^0$-meson production at mid-rapidity in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{200\,GeV}}$

The STAR collaboration Adam, Jaroslav ; Adamczyk, Leszek ; Adams, Joseph ; et al.
Phys.Rev.C 99 (2019) 034908, 2019.
Inspire Record 1711377 DOI 10.17182/hepdata.95750

We report a new measurement of $D^0$-meson production at mid-rapidity ($|y|$\,$<$\,1) in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{200\,GeV}}$ utilizing the Heavy Flavor Tracker, a high resolution silicon detector at the STAR experiment. Invariant yields of $D^0$-mesons with transverse momentum $p_{T}$ $\lesssim 9$\,GeV/$c$ are reported in various centrality bins (0--10\%, 10--20\%, 20--40\%, 40--60\% and 60--80\%). Blast-Wave thermal models are used to fit the $D^0$-meson $p_{T}$ spectra to study $D^0$ hadron kinetic freeze-out properties. The average radial flow velocity extracted from the fit is considerably smaller than that of light hadrons ($\pi,K$ and $p$), but comparable to that of hadrons containing multiple strange quarks ($\phi,\Xi^-$), indicating that $D^0$ mesons kinetically decouple from the system earlier than light hadrons. The calculated $D^0$ nuclear modification factors re-affirm that charm quarks suffer large amount of energy loss in the medium, similar to those of light quarks for $p_{T}$\,$>$\,4\,GeV/$c$ in central 0--10\% Au+Au collisions. At low $p_{T}$, the nuclear modification factors show a characteristic structure qualitatively consistent with the expectation from model predictions that charm quarks gain sizable collective motion during the medium evolution. The improved measurements are expected to offer new constraints to model calculations and help gain further insights into the hot and dense medium created in these collisions.

22 data tables

$D^0$ (in terms of (D0 +D0)/2)) invariant yield at mid-rapidity ($|y| < 1$) vs transverse momentum for different centrality classes. Error bars indicate statistical uncertainties and brackets depict systematic uncertainties. Global systematic uncertainties in B.R. are not plotted. Solid and dashed lines depict Levy function fits.

$D^0$ (in terms of (D0 +D0)/2)) spectra in pp collisions. Note, the $\sigma_{NSD}$ = 30 $m$b for p+p was used in the calculations.

Integrated $D^0$ cross section per nucleon-nucleon collision at mid-rapidity for $p_T >0$ (a) and $p_T >4$ GeV/c (b) as a function of centrality $N_{part}$. The statistical and systematic uncertainties are shown as error bars and brackets on the data points. The green boxes on the data points depict the overall normalization uncertainties in p+p and Au+Au data respectively.

More…

Two-particle Bose-Einstein correlations and their Lévy parameters in PbPb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

The CMS collaboration Tumasyan, Armen ; Adam, Wolfgang ; Andrejkovic, Janik Walter ; et al.
Phys.Rev.C 109 (2024) 024914, 2024.
Inspire Record 2670243 DOI 10.17182/hepdata.134676

Two-particle Bose-Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data sample, containing 4.27 $\times$$10^{9}$ minimum bias events corresponding to an integrated luminosity of 0.607 nb$^{-1}$, was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter ($\alpha$), the Lévy scale parameter ($R$), and the correlation strength parameter ($\lambda$). The source shape, characterized by $\alpha$, is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy $R$ parameter. The $\lambda$ parameter is studied in terms of the core-halo model.

12 data tables

The correlation function $C_2(q)$ of negatively charged hadron pairs with $1.00<K_{\mathrm{T}}<1.05$ GeV/$c$ in the 10-20% centrality bin.

The double-ratio correlation function $DR(q)$ of negatively charged hadron pairs with $1.30<K_{\mathrm{T}}<1.35$ GeV/$c$ in the 20-30% centrality bin.

The two-particle correlation function of negatively charged hadron pairs with $0.9<K_{\mathrm{T}}<1.3$ GeV/$c$ in the 0-100% centrality range, calculated using MC events with (Reco) and without (Gen) detector reconstruction.

More…

First Observation of Directed Flow of Hypernuclei $^3_{\Lambda}$H and $^4_{\Lambda}$H in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au Collisions at RHIC

The STAR collaboration Aboona, Bassam ; Adam, Jaroslav ; Adams, Joseph ; et al.
Phys.Rev.Lett. 130 (2023) 212301, 2023.
Inspire Record 2605845 DOI 10.17182/hepdata.136028

We report here the first observation of directed flow ($v_1$) of the hypernuclei $^3_{\Lambda}$H and $^4_{\Lambda}$H in mid-central Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165 $\times$ 10$^{6}$ events in 5%-40% centrality, about 8400 $^3_{\Lambda}$H and 5200 $^4_{\Lambda}$H candidates are reconstructed through two- and three-body decay channels. We observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity $v_1$ slopes of $^3_{\Lambda}$H and $^4_{\Lambda}$H follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in such collisions.

8 data tables

$\Lambda$ hyperon and hypernuclei directed flow $v_1$, shown as a function of rapidity, from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. In the case of $^{3}_{\Lambda}$H $v_1$, both two-body (dots) and three-body (triangles) decays are used. The linear terms of the fitting for $#Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H are shown as the yellow-red lines. The rapidity dependence of $v_1$ for $p$, $d$, $t$, $^3$He, and $^4$He are also shown as open markers (circles, diamonds, up-triangles, down-triangles and squares), and the linear terms of the fitting results are shown as dashed lines in the positive rapidity region.

$\Lambda$ hyperon and hypernuclei directed flow $v_1$, shown as a function of rapidity, from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. In the case of $^{3}_{\Lambda}$H $v_1$, both two-body (dots) and three-body (triangles) decays are used. The linear terms of the fitting for $#Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H are shown as the yellow-red lines. The rapidity dependence of $v_1$ for $p$, $d$, $t$, $^3$He, and $^4$He are also shown as open markers (circles, diamonds, up-triangles, down-triangles and squares), and the linear terms of the fitting results are shown as dashed lines in the positive rapidity region.

$\Lambda$ hyperon and hypernuclei directed flow $v_1$, shown as a function of rapidity, from the $\sqrt{s_{NN}}$ = 3 GeV 5-40% mid-central Au+Au collisions. In the case of $^{3}_{\Lambda}$H $v_1$, both two-body (dots) and three-body (triangles) decays are used. The linear terms of the fitting for $#Lambda$, $^{3}_{\Lambda}$H and $^{4}_{\Lambda}$H are shown as the yellow-red lines. The rapidity dependence of $v_1$ for $p$, $d$, $t$, $^3$He, and $^4$He are also shown as open markers (circles, diamonds, up-triangles, down-triangles and squares), and the linear terms of the fitting results are shown as dashed lines in the positive rapidity region.

More…

Beam Energy Dependence of Triton Production and Yield Ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$) in Au+Au Collisions at RHIC

The STAR collaboration Abdulhamid, Muhammad ; Aboona, Bassam ; Adam, Jaroslav ; et al.
Phys.Rev.Lett. 130 (2023) 202301, 2023.
Inspire Record 2152917 DOI 10.17182/hepdata.133992

We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

68 data tables

Invariant yields of tritons at 7.7 GeV, all centralities. The first uncertainty is statistical uncertainty, the second is systematic uncertainty.

Invariant yields of tritons at 11.5 GeV, all centralities. The first uncertainty is statistical uncertainty, the second is systematic uncertainty.

Invariant yields of tritons at 14.5 GeV, all centralities. The first uncertainty is statistical uncertainty, the second is systematic uncertainty.

More…