The ALICE experiment has measured the inclusive J/$\psi$ production in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}} } = 2.76$ TeV down to zero transverse momentum in the rapidity range $2.5 < y < 4$. A suppression of the inclusive J/$\psi$ yield in Pb-Pb is observed with respect to the one measured in pp collisions scaled by the number of binary nucleon-nucleon collisions. The nuclear modification factor, integrated over the 0-80% most central collisions, is $0.545 \pm 0.032 \rm{(stat.)} \pm 0.083 \rm{(syst.)}$ and does not exhibit a significant dependence on the collision centrality. These features appear significantly different from measurements at lower collision energies. Models including J/$\psi$ production from charm quarks in a deconfined partonic phase can describe our data.

The first measurements of the scattering parameters of $\Lambda$K pairs in all three charge combinations ($\Lambda$K$^{+}$, $\Lambda$K$^{-}$, and $\Lambda\mathrm{K^{0}_{S}}$) are presented. The results are achieved through a femtoscopic analysis of $\Lambda$K correlations in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 2.76 TeV recorded by ALICE at the LHC. The femtoscopic correlations result from strong final-state interactions, and are fit with a parametrization allowing for both the characterization of the pair emission source and the measurement of the scattering parameters for the particle pairs. Extensive studies with the THERMINATOR 2 event generator provide a good description of the non-femtoscopic background, which results mainly from collective effects, with unprecedented precision. Furthermore, together with HIJING simulations, this model is used to account for contributions from residual correlations induced by feed-down from particle decays. The extracted scattering parameters indicate that the strong force is repulsive in the $\Lambda\rm{K}^{+}$ interaction and attractive in the $\Lambda\rm{K}^{-}$ interaction. The data hint that the and $\Lambda\rm{K}^{0}_{S}$ interaction is attractive, however the uncertainty of the result does not permit such a decisive conclusion. The results suggest an effect arising either from different quark-antiquark interactions between the pairs ($\rm s\overline{s}$ in $\Lambda$K$^{+}$ and $\rm u\overline{u}$ in $\Lambda$K$^{-}$) or from different net strangeness for each system (S = 0 for $\Lambda$K$^{+}$, and S = $-2$ for $\Lambda$K$^{-}$). Finally, the $\Lambda$K systems exhibit source radii larger than expected from extrapolation from identical particle femtoscopic studies. This effect is interpreted as resulting from the separation in space-time of the single-particle $\Lambda$ and K source distributions.

A measurement of the production of prompt $\Lambda_{\rm c}^{+}$ baryons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02 $ TeV with the ALICE detector at the LHC is reported. The $\Lambda_{\rm c}^{+}$ and $\overline{\Lambda}_{\rm c}^{-}$ were reconstructed at midrapidity ($|y| < 0.5$) via the hadronic decay channel $\Lambda_{\rm c}^{+}\rightarrow {\rm p} {\rm K}_{\rm S}^{0}$ (and charge conjugate) in the transverse momentum and centrality intervals $6 < p_{\rm T} <12$ GeV/$c$ and 0-80%. The $\Lambda_{\rm c}^{+}$/D$^0$ ratio, which is sensitive to the charm quark hadronisation mechanisms in the medium, is measured and found to be larger than the ratio measured in minimum-bias pp collisions at $\sqrt{s} = 7$ TeV and in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02 $ TeV. In particular, the values in p-Pb and Pb-Pb collisions differ by about two standard deviations of the combined statistical and systematic uncertainties in the common $p_{\rm T}$ interval covered by the measurements in the two collision system. The $\Lambda_{\rm c}^{+}$/D$^0$ ratio is also compared with model calculations including different implementations of charm quark hadronisation. The measured ratio is reproduced by models implementing a pure coalescence scenario, while adding a fragmentation contribution leads to an underestimation. The $\Lambda_{\rm c}^{+}$ nuclear modification factor, $R_{\rm AA}$, is also presented. The measured values of the $R_{\rm AA}$ of $\Lambda_{\rm c}^{+}$, D$_{\rm s}$ and non-strange D mesons are compatible within the combined statistical and systematic uncertainties. They show, however, a hint of a hierarchy $(R_{\rm AA}^{{\rm D}^{0}}<R_{\rm AA}^{{\rm D}_{\rm s}}<R_{\rm AA}^{\Lambda_{\rm c}^{+}})$, conceivable with a contribution of recombination mechanisms to charm hadron formation in the medium.

The $p_{\rm T}$-differential production cross section of prompt $\Lambda_{\rm c}^+$ charmed baryons was measured with the ALICE detector at the Large Hadron Collider (LHC) in pp collisions at $\sqrt{s} = 7$ TeV and in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV at midrapidity. The $\Lambda_{\rm c}^+$ and ${\overline{\Lambda}}_{\rm c}^-$ were reconstructed in the hadronic decay modes $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K^-}\pi^+$, $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K_{\rm S}^0}$ and in the semileptonic channel $\Lambda_{\rm c}^+\rightarrow {\rm e^+}\nu_{\rm e}\Lambda$ (and charge conjugates). The measured values of the $\Lambda_{\rm c}^+/{\rm D_0}$ ratio, which is sensitive to the c-quark hadronisation mechanism, and in particular to the production of baryons, are presented and are larger than those measured previously in different colliding systems, centre-of-mass energies, rapidity and $p_{\rm T}$ intervals, where the $\Lambda_{\rm c}^+$ production process may differ. The results are compared with the expectations obtained from perturbative Quantum Chromodynamics calculations and Monte Carlo event generators. Neither perturbative QCD calculations nor Monte Carlo models reproduce the data, indicating that the fragmentation of heavy-flavour baryons is not well understood. The first measurement at the LHC of the $\Lambda_{\rm c}^+$ nuclear modification factor, $R_{\rm pPb}$, is also presented. The $R_{\rm pPb}$ is found to be consistent with unity and with that of D mesons within the uncertainties, and consistent with a theoretical calculation that includes cold nuclear matter effects and a calculation that includes charm quark interactions with a deconfined medium.

Hadronic resonances are used to probe the hadron gas produced in the late stage of heavy-ion collisions since they decay on the same timescale, of the order of 1 to 10 fm/$c$, as the decoupling time of the system. In the hadron gas, (pseudo)elastic scatterings among the products of resonances that decayed before the kinetic freeze-out and regeneration processes counteract each other, the net effect depending on the resonance lifetime, the duration of the hadronic phase, and the hadronic cross sections at play. In this context, the $\Sigma(1385)^{\pm}$ particle is of particular interest as models predict that regeneration dominates over rescattering despite its relatively short lifetime of about 5.5 fm/$c$. The first measurement of the $\Sigma(1385)^{\pm}$ resonance production at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}= 5.02$ TeV with the ALICE detector is presented in this Letter. The resonances are reconstructed via their hadronic decay channel, $\Lambda\pi$, as a function of the transverse momentum ($p_{\rm T}$) and the collision centrality. The results are discussed in comparison with the measured yield of pions and with expectations from the statistical hadronization model as well as commonly employed event generators, including PYTHIA8/Angantyr and EPOS3 coupled to the UrQMD hadronic cascade afterburner. None of the models can describe the data. For $\Sigma(1385)^{\pm}$, a similar behaviour as ${\rm K}^{*} (892)^{0}$ is observed in data unlike the predictions of EPOS3 with afterburner.

The production of $\Upsilon$ mesons in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5 TeV is measured with the muon spectrometer of the ALICE detector at the LHC. The yields as well as the nuclear modification factors are determined in the forward rapidity region $2.5<y<4.0$, as a function of rapidity, transverse momentum and collision centrality. The results show that the production of the $\Upsilon$(1S) meson is suppressed by a factor of about three with respect to the production in proton-proton collisions. For the first time, a significant signal for the $\Upsilon$(2S) meson is observed at forward rapidity, indicating a suppression stronger by about a factor 2-3 with respect to the ground state. The measurements are compared with transport, hydrodynamic, comover and statistical hadronisation model calculations.

Inclusive $\Upsilon$(1S) and $\Upsilon$(2S) production have been measured in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{_{\rm NN}}}=5.02$ TeV, using the ALICE detector at the CERN LHC. The $\Upsilon$ mesons are reconstructed in the centre-of-mass rapidity interval $2.5<y<4$ and in the transverse-momentum range $p_{\rm T}<15$ GeV/$c$, via their decays to muon pairs. In this Letter, we present results on the inclusive $\Upsilon$(1S) nuclear modification factor $R_{\rm AA}$ as a function of collision centrality, transverse momentum and rapidity. The $\Upsilon$(1S) and $\Upsilon$(2S) $R_{\rm AA}$, integrated over the centrality range 0-90%, are $0.37 \pm 0.02 {\rm{(stat)}}\pm 0.03 {\rm{(syst)}}$ and $0.10 \pm 0.04 {\rm{(stat)}}\pm 0.02 {\rm{(syst)}}$, respectively, leading to a ratio $R_{\rm{AA}}^{\Upsilon(\rm2S)}/R_{\rm{AA}}^{\Upsilon(\rm1S)}$ of $0.28\pm0.12\text{(stat)}\pm0.06\text{(syst)}$. The observed $\Upsilon$(1S) suppression increases with the centrality of the collision and no significant variation is observed as a function of transverse momentum and rapidity.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

$\phi$ meson measurements provide insight into strangeness production, which is one of the key observables for the hot medium formed in high-energy heavy-ion collisions. ALICE measured $\phi$ production through its decay in muon pairs in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 2.76 TeV in the intermediate transverse momentum range $2 < p_{\rm T} < 5$ GeV/$c$ and in the rapidity interval $2.5<y<4$. The $\phi$ yield was measured as a function of the transverse momentum and collision centrality. The nuclear modification factor was obtained as a function of the average number of participating nucleons. Results were compared with the ones obtained via the kaon decay channel in the same $p_{\rm T}$ range at midrapidity. The values of the nuclear modification factor in the two rapidity regions are in agreement within uncertainties.

The first study of $\phi$-meson production in p-Pb collisions at forward and backward rapidity, at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$~TeV, has been performed with the ALICE apparatus at the LHC. The $\phi$-mesons have been identified in the dimuon decay channel in the transverse momentum ($p_{\rm T}$) range $1 < p_{\rm T} < 7$ GeV/$c$, both in the p-going ($2.03 < y < 3.53$) and the Pb-going ($-4.46 < y < -2.96$) directions, where $y$ stands for the rapidity in the nucleon-nucleon centre-of-mass, the integrated luminosity amounting to $5.01 \pm 0.19$~nb$^{-1}$ and $5.81 \pm 0.20$~nb$^{-1}$, respectively, for the two data samples. Differential cross sections as a function of transverse momentum and rapidity are presented. The forward-backward ratio for $\phi$-meson production is measured for $2.96<|y|<3.53$, resulting in a ratio $\sim 0.5$ with no significant $p_{\rm T}$ dependence within the uncertainties. The $p_{\rm T}$ dependence of the $\phi$ nuclear modification factor $R_{\rm pPb}$ exhibits an enhancement up to a factor 1.6 at $p_{\rm T}$ = 3-4 GeV/$c$ in the Pb-going direction. The $p_{\rm T}$ dependence of the $\phi$-meson cross section in pp collisions at $\sqrt{s}$ = 2.76 TeV, which is used to determine a reference for the p-Pb results, is also presented here for $1 < p_{\rm T} < 5$ GeV/$c$ and $2.5 <y < 4$ for a $78 \pm 3$~nb$^{-1}$ integrated luminosity sample.

An invariant differential cross section measurement of inclusive $\pi^{0}$ and $\eta$ meson production at mid-rapidity in pp collisions at $\sqrt{s}=8$ TeV was carried out by the ALICE experiment at the LHC. The spectra of $\pi^{0}$ and $\eta$ mesons were measured in transverse momentum ranges of $0.3<p_{\rm T}<35$ GeV/$c$ and $0.5<p_{\rm T}<35$ GeV/$c$, respectively. Next-to-leading order perturbative QCD calculations using fragmentation functions DSS14 for the $\pi^{0}$ and AESSS for the $\eta$ overestimate the cross sections of both neutral mesons, although such calculations agree with the measured $\eta/\pi^{0}$ ratio within uncertainties. The results were also compared with PYTHIA~8.2 predictions for which the Monash~2013 tune yields the best agreement with the measured neutral meson spectra. The measurements confirm a universal behavior of the $\eta/\pi^{0}$ ratio seen for NA27, PHENIX and ALICE data for pp collisions from $\sqrt{s}=27.5$ GeV to $\sqrt{s}=8$ TeV within experimental uncertainties. A relation between the $\pi^{0}$ and $\eta$ production cross sections for pp collisions at $\sqrt{s}=8$ TeV is given by $m_{\rm T}$ scaling for $p_{\rm T}>3.5$ GeV/$c$. However, a deviation from this empirical scaling rule is observed for transverse momenta below $p_{\rm T}<3.5$ GeV/$c$ in the $\eta/\pi^0$ ratio with a significance of $6.2\sigma$.

The production of the $\psi(2S)$ charmonium state was measured with ALICE in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the dimuon decay channel. A significant signal was observed for the first time at LHC energies down to zero transverse momentum, at forward rapidity ($2.5<y<4$). The measurement of the ratio of the inclusive production cross sections of the $\psi(2S)$ and J/$\psi$ resonances is reported as a function of the centrality of the collisions and of transverse momentum, in the region $p_{\rm T}<12$ GeV/$c$. The results are compared with the corresponding measurements in pp collisions, by forming the double ratio $[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{Pb-Pb}}/[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{pp}}$. It is found that in Pb-Pb collisions the $\psi(2S)$ is suppressed by a factor of $\sim 2$ with respect to the J/$\psi$. The $\psi(2S)$ nuclear modification factor $R_{\rm AA}$ was also obtained as a function of both centrality and $p_{\rm T}$. The results show that the $\psi(2S)$ resonance yield is strongly suppressed in Pb-Pb collisions, by a factor up to $\sim 3$ with respect to pp. Comparisons of cross section ratios with previous SPS findings by the NA50 experiment and of $R_{\rm AA}$ with higher-$p_{\rm T}$ results at LHC energy are also reported. These results and the corresponding comparisons with calculations of transport and statistical models address questions on the presence and properties of charmonium states in the quark-gluon plasma formed in nuclear collisions at the LHC.

The prompt production of the charm baryon $\rm \Lambda_{c}^{+}$ and the $\rm \Lambda_{c}^{+}/\mathrm {D^0}$ production ratios were measured at midrapidity with the ALICE detector in pp and p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$TeV. These new measurements show a clear decrease of the $\rm \Lambda_{c}^{+}/\mathrm {D^0}$ ratio with increasing transverse momentum ($p_{\rm T}$) in both collision systems in the range $2<p_{\rm T}<12$ GeV/$c$, exhibiting similarities with the light-flavour baryon-to-meson ratios ${\rm p}/\pi$ and $\Lambda/\mathrm {K^0_S}$. At low $p_{\rm T}$, predictions that include additional colour-reconnection mechanisms beyond the leading-colour approximation; assume the existence of additional higher-mass charm-baryon states; or include hadronisation via coalescence can describe the data, while predictions driven by charm-quark fragmentation processes measured in $\mathrm {e^+e^-}$ and $\mathrm {e^-p}$ collisions significantly underestimate the data. The results presented in this letter provide significant evidence that the established assumption of universality (colliding-system independence) of parton-to-hadron fragmentation is not sufficient to describe charm-baryon production in hadronic collisions at LHC energies.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ in pp collisions at $\sqrt{s}$ = 8 TeV was measured using Run 1 data collected by the ALICE collaboration at the LHC. The $p_{\rm{T}}$-differential yields d$^{\rm 2}N$/d$y$d$p_{\rm{T}}$ in the range $0 < p_{\rm{T}} < 20$ GeV/$c$ for $\rm{K}^{*0}$ and $0.4 < p_{\rm{T}} < 16$ GeV/$c$ for $\phi$ have been measured at midrapidity, $|y| < 0.5$. Moreover, improved measurements of the $\rm{K}^{*}(892)^{0}$ and $\phi(1020)$ at $\sqrt{s} = 7$TeV are presented. The collision energy dependence of $p_{\rm{T}}$ distributions, $p_{\rm{T}}$-integrated yields and particle ratios in inelastic pp collisions are examined. The results are also compared with different collision systems. The values of the particle ratios are found to be similar to those measured at other LHC energies. In pp collisions a hardening of the particle spectra is observed with increasing energy, but at the same time it is also observed that the relative particle abundances are independent of the collision energy. The $p_{\rm{T}}$-differential yields of $\rm{K}^{*0}$ and $\phi$ in pp collisions at $\sqrt{s} = 8$ TeV are compared with the expectations of different Monte Carlo event generators.

The production of the hypertriton nuclei $^{3}_{\Lambda}\mathrm H$ and $^{3}_{\bar{\Lambda}} \overline{\mathrm H}$ has been measured for the first time in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV with the ALICE experiment at LHC energies. The total yield, d$N$/d$y$ $\times \mathrm{B.R.}_{\left( ^{3}_{\Lambda}\mathrm H \rightarrow ^{3}\mathrm{He},\pi^{-} \right)} = \left( 3.86 \pm 0.77 (\mathrm{stat.}) \pm 0.68 (\mathrm{syst.})\right) \times 10^{-5}$ in the 0-10% most central collisions, is consistent with the predictions from a statistical thermal model using the same temperature as for the light hadrons. The coalescence parameter $B_3$ shows a dependence on the transverse momentum, similar to the $B_2$ of deuterons and the $B_3$ of $^{3}\mathrm{He}$ nuclei. The ratio of yields $S_3$ = $^{3}_{\Lambda}\mathrm H$/($^{3}\mathrm{He}$ $\times \Lambda/\mathrm{p}$) was measured to be $S_3$ = 0.60 $\pm$ 0.13 (stat.) $\pm$ 0.21 (syst.) in 0-10% centrality events; this value is compared to different theoretical models. The measured $S_3$ is fully compatible with thermal model predictions. The measured $^{3}_{\Lambda}\mathrm H$ lifetime, $ \tau = 181^{+54}_{-39} (\mathrm{stat.}) \pm 33 (\mathrm{syst.})\ \mathrm{ps}$ is compatible within 1$\sigma$ with the world average value.

The measurement of the production of ${\rm f}_{0}(980)$ in inelastic pp collisions at $\sqrt{s} = 5.02$ TeV is presented. This is the first reported measurement of inclusive ${\rm f}_{0}(980)$ yield at LHC energies. The production is measured at midrapidity, $|y| < 0.5$, in a wide transverse momentum range, $0 < p_{\rm T} < 16$ GeV/$c$, by reconstructing the resonance in the ${\rm f}_{0}(980) \rightarrow \pi^{+}\pi^{-}$ hadronic decay channel using the ALICE detector. The $p_{\rm T}$-differential yields are compared to those of pions, protons and $\phi$ mesons as well as to predictions from the HERWIG 7.2 QCD-inspired Monte Carlo event generator and calculations from a coalescence model that uses the AMPT model as an input. The ratio of the $p_{\rm T}$-integrated yield of ${\rm f}_{0}(980)$ relative to pions is compared to measurements in ${\rm e}^{+}{\rm e}^{-}$ and pp collisions at lower energies and predictions from statistical hadronisation models and HERWIG 7.2. A mild collision energy dependence of the ${\rm f}_{0}(980)$ to pion production is observed in pp collisions from SPS to LHC energies. All considered models underpredict the $p_{\rm T}$-integrated $2{\rm f}_{0}(980)/(\pi^{+}+\pi^{-})$ ratio. The prediction from the canonical statistical hadronisation model assuming a zero total strangeness content of ${\rm f}_{0}(980)$ is consistent with the data within 1.9$\sigma$ and is the closest to the data. The results provide an essential reference for future measurements of the particle yield and nuclear modification in p$-$Pb and Pb$-$Pb collisions, which have been proposed to be instrumental to probe the elusive nature and quark composition of the ${\rm f}_{0}(980)$ scalar meson.

The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

Charged-particle spectra at midrapidity are measured in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV and presented in centrality classes ranging from most central (0-5%) to most peripheral (95-100%) collisions. Possible medium effects are quantified using the nuclear modification factor ($R_{\rm AA}$) by comparing the measured spectra with those from proton-proton collisions, scaled by the number of independent nucleon-nucleon collisions obtained from a Glauber model. At large transverse momenta ($8<p_{\rm T}<20$ GeV/$c$), the average $R_{\rm AA}$ is found to increase from about $0.15$ in 0-5% central to a maximum value of about $0.8$ in 75-85% peripheral collisions, beyond which it falls off strongly to below $0.2$ for the most peripheral collisions. Furthermore, $R_{\rm AA}$ initially exhibits a positive slope as a function of $p_{\rm T}$ in the $8$-$20$ GeV/$c$ interval, while for collisions beyond the 80% class the slope is negative. To reduce uncertainties related to event selection and normalization, we also provide the ratio of $R_{\rm AA}$ in adjacent centrality intervals. Our results in peripheral collisions are consistent with a PYTHIA-based model without nuclear modification, demonstrating that biases caused by the event selection and collision geometry can lead to the apparent suppression in peripheral collisions. This explains the unintuitive observation that $R_{\rm AA}$ is below unity in peripheral Pb-Pb, but equal to unity in minimum-bias p-Pb collisions despite similar charged-particle multiplicities.

The first measurements of elliptic flow of $\pi^\pm$, ${\rm K}^\pm$, p+$\overline{\rm p}$, ${\rm K_{S}^0}$, $\Lambda$+$\overline{\Lambda}$, $\phi$, $\Xi^-$+$\Xi^+$, and $\Omega^-$+$\Omega^+$ using multiparticle cumulants in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented. Results obtained with two- ($v_2\{2\}$) and four-particle cumulants ($v_2\{4\}$) are shown as a function of transverse momentum, $p_{\rm T}$, for various collision centrality intervals. Combining the data for both $v_2\{2\}$ and $v_2\{4\}$ also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark-gluon plasma. The characteristic features observed in previous $p_{\rm T}$-differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low $p_{\rm T}$ and the approximate scaling with the number of constituent quarks at intermediate $p_{\rm T}$, are similarly present in the four-particle correlations and the combinations of $v_2\{2\}$ and $v_2\{4\}$. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide $p_{\rm T}$ range.

Measurements of elliptic ($v_2$) and triangular ($v_3$) flow coefficients of $\pi^{\pm}$, K$^{\pm}$, p+$\rm \overline{p}$, K$^0_{\rm S}$, and $\Lambda + \overline{\Lambda}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range $\left | y \right |<0.5$ and reported as a function of transverse momentum, $p_{\rm T}$, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T}<3$ GeV/$c$, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3< $p_{\rm T}$ <8 GeV/$c$). The magnitude of the baryon $v_{2}$ is larger than that of mesons up to $p_{\rm T}$ = 6 GeV/$c$. The centrality dependence of the shape evolution of the $p_{\rm T}$-differential $v_2$ is studied for the various hadron species. The $v_2$ coefficients of $\pi^{\pm}$, K$^{\pm}$, and p+$\rm \overline{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for $p_{\rm T} <1$ GeV/$c$. A comparison with $v_{\rm n}$ measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV yields an enhanced $v_2$ in central collisions and diminished value in semicentral collisions.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.

The production of the prompt charmed mesons $D^0$, $D^+$ and $D^{*+}$ relative to the reaction plane was measured in Pb-Pb collisions at a centre-of-mass energy per nucleon-nucleon collision of $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the LHC. D mesons were reconstructed via their hadronic decays at central rapidity in the transverse momentum ($p_{\rm T}$) interval of 2-16 GeV/$c$. The azimuthal anisotropy is quantified in terms of the second coefficient $v_2$ in a Fourier expansion of the D meson azimuthal distribution, and in terms of the nuclear modification factor $R_{\rm AA}$, measured in the direction of the reaction plane and orthogonal to it. The $v_2$ coefficient was measured with three different methods and in three centrality classes in the interval 0-50%. A positive $v_2$ is observed in mid-central collisions (30-50% centrality class), with an mean value of $0.204_{-0.036}^{+0.099}$ (tot.unc.) in the interval $2 < p_{\rm T} < 6$ GeV/$c$, which decreases towards more central collisions (10-30% and 0-10% classes). The positive $v_2$ is also reflected in the nuclear modification factor, which shows a stronger suppression in the direction orthogonal to the reaction plane for mid-central collisions. The measurements are compared to theoretical calculations of charm quark transport and energy loss in high-density strongly-interacting matter at high temperature. The models that include substantial elastic interactions with an expanding medium provide a good description of the observed anisotropy. However, they are challenged to simultaneously describe the strong suppression of high-$p_{\rm T}$ yield of D mesons in central collisions and their azimuthal anisotropy in non-central collisions.

We present measurements of the azimuthal dependence of charged jet production in central and semi-central $\sqrt{s_{\mathrm{NN}}}$ = 2.76 TeV Pb-Pb collisions with respect to the second harmonic event plane, quantified as $v_{2}^{\mathrm{ch~jet}}$. Jet finding is performed employing the anti-$k_{\mathrm{T}}$ algorithm with a resolution parameter $R$ = 0.2 using charged tracks from the ALICE tracking system. The contribution of the azimuthal anisotropy of the underlying event is taken into account event-by-event. The remaining (statistical) region-to-region fluctuations are removed on an ensemble basis by unfolding the jet spectra for different event plane orientations independently. Significant non-zero $v_{2}^{\mathrm{ch~jet}}$ is observed in semi-central collisions (30-50\% centrality) for 20 $<$ $p_{\mathrm{T}}^{\rm ch~jet}$ $<$ 90 ${\mathrm{GeV}\kern-0.05em/\kern-0.02em c}$. The azimuthal dependence of the charged jet production is similar to the dependence observed for jets comprising both charged and neutral fragments, and compatible with measurements of the $v_2$ of single charged particles at high $p_{\mathrm{T}}$. Good agreement between the data and predictions from JEWEL, an event generator simulating parton shower evolution in the presence of a dense QCD medium, is found in semi-central collisions.

The measurement of the azimuthal-correlation function of prompt D mesons with charged particles in pp collisions at $\sqrt{s}$ = 5.02 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC is reported. The D$^{\rm 0}$, D$^{\rm +}$, and D$^{\rm *+}$ mesons, together with their charge conjugates, were reconstructed at midrapidity in the transverse momentum interval 3 < $p_{\rm T}$ < 24 GeV/c and correlated with charged particles having $p_{\rm T}$ > 0.3 GeV/c and pseudorapidity $|\eta| <$ 0.8. The properties of the correlation peaks appearing in the near- and away-side regions (for $\Delta \varphi \approx$ 0 and $\Delta \varphi \approx \pi$, respectively) were extracted via a fit to the azimuthal correlation functions. The shape of the correlation functions and the near- and away-side peak features are found to be consistent in pp and p-Pb collisions, showing no modifications due to nuclear effects within uncertainties. The results are compared with predictions from Monte Carlo simulations performed with the PYTHIA, POWHEG+PYTHIA, HERWIG, and EPOS 3 event generators.

We present the first azimuthally differential measurements of the pion source size relative to the second harmonic event plane in Pb-Pb collisions at a center-of-mass energy per nucleon-nucleon pair of $\sqrt{s_{\rm NN}}=2.76$ TeV. The measurements have been performed in the centrality range 0-50% and for pion pair transverse momenta $0.2 < k_{\rm T} < 0.7$ GeV/$c$. We find that the $R_{\rm side}$ and $R_{\rm out}$ radii, which characterize the pion source size in the directions perpendicular and parallel to the pion transverse momentum, oscillate out of phase, similar to what was observed at the Relativistic Heavy Ion Collider (RHIC). The final-state source eccentricity, estimated via $R_{\rm side}$ oscillations, is found to be significantly smaller than the initial-state source eccentricity, but remains positive; indicating that even after a stronger expansion in the in-plane direction, the pion source at the freeze-out is still elongated in the out-of-plane direction. The 3+1D hydrodynamic calculations are in qualitative agreement with observed centrality and transverse momentum $R_{\rm side}$ oscillations, but systematically underestimate the oscillation magnitude.