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.

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.

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.

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 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).

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.