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.

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 K$^*$(892)$^\pm$ meson resonance is measured at midrapidity ($|y|<0.5$) in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV using the ALICE detector at the LHC. The resonance is reconstructed via its hadronic decay channel K$^*$(892)$^\pm \rightarrow \rm{K^0_S \pi^\pm}$. The transverse momentum distributions are obtained for various centrality intervals in the $p_{\rm T}$ range of 0.4-16 GeV/$c$. The reported measurements of integrated yields, mean transverse momenta, and particle yield ratios are consistent with previous ALICE measurements for K$^*$(892)$^0$. The $p_{\rm T}$-integrated yield ratio 2K$^*$(892)$^\pm$/($\rm{K^+ + K^-}$) in central Pb-Pb collisions shows a significant suppression (9.3$\sigma$) relative to pp collisions. Thermal model calculations overpredict the particle yield ratio. Although both simulations consider the hadronic phase, only HRG-PCE accurately represents the measurements, whereas MUSIC+SMASH tends to overpredict them. These observations, along with the kinetic freeze-out temperatures extracted from the yields of light-flavored hadrons using the HRG-PCE model, indicate a finite hadronic phase lifetime, which increases towards central collisions. The $p_{\rm T}$-differential yield ratios 2K$^*$(892)$^\pm$/($\rm{K^+ + K^-}$) and 2K$^*$(892)$^\pm$/($\rm{\pi^+ + \pi^-}$) are suppressed by up to a factor of five at $p_{\rm T}<2$ GeV/$c$ in central Pb-Pb collisions compared to pp collisions at $\sqrt{s} =$ 5.02 TeV. Both particle ratios and are qualitatively consistent with expectations for rescattering effects in the hadronic phase. The nuclear modification factor shows a smooth evolution with centrality and is below unity at $p_{\rm T}>8$ GeV/$c$, consistent with measurements for other light-flavored hadrons. The smallest values are observed in most central collisions, indicating larger energy loss of partons traversing the dense medium.

An excess of J/$\psi$ yield at very low transverse momentum ($p_{\rm T} < 0.3$ GeV/$c$), originating from coherent photoproduction, is observed in peripheral and semicentral hadronic Pb$-$Pb collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}} = 5.02$ TeV. The measurement is performed with the ALICE detector via the dimuon decay channel at forward rapidity ($2.5<y<4$). The nuclear modification factor at very low $p_{\rm T}$ and the coherent photoproduction cross section are measured as a function of centrality down to the 10% most central collisions. These results extend the previous study at $\sqrt{s_{\rm NN}} = 2.76$ TeV, confirming the clear excess over hadronic production in the $p_{\rm T}$ range 0$-$0.3 GeV/$c$ and the centrality range 70$-$90%, and establishing an excess with a significance greater than 5$\sigma$ also in the 50$-$70% and 30$-$50% centrality ranges. The results are compared with earlier measurements at $\sqrt{s_{\rm NN}} = 2.76$ TeV and with different theoretical predictions aiming at describing how coherent photoproduction occurs in hadronic interactions with nuclear overlap.

The production of electrons from beauty-hadron decays was measured at midrapidity in proton-proton (pp) and central Pb-Pb collisions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV, using the ALICE detector at the LHC. The cross section measured in pp collisions in the transverse momentum interval $2 < p_{\rm T} < 8$ GeV/$c$ was compared with models based on perturbative quantum chromodynamics calculations. The yield in the 10% most central Pb-Pb collisions, measured in the interval $2 < p_{\rm T} < 26$ GeV/$c$, was used to compute the nuclear modification factor $R_{\rm AA}$, extrapolating the pp reference cross section to $p_{\rm T}$ larger than 8 GeV/$c$. The measured $R_{\rm AA}$ shows significant suppression of the yield of electrons from beauty-hadron decays at high $p_{\rm T}$ and does not show a significant dependence on $p_{\rm T}$ above 8 GeV/$c$ within uncertainties. The results are described by several theoretical models based on different implementations of the interaction of heavy quarks with a quark-gluon plasma, which predict a smaller energy loss for beauty quarks compared to light and charm quarks.

The transverse momentum ($p_{\rm T}$) and centrality dependence of the nuclear modification factor $R_{\rm AA}$ of prompt and non-prompt J$/\psi$, the latter originating from the weak decays of beauty hadrons, have been measured by the ALICE collaboration in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The measurements are carried out through the ${\rm e}^{+}{\rm e}^{-}$ decay channel at midrapidity ($|y| < 0.9$) in the transverse momentum region $1.5 < p_{\rm T} < 10$ GeV/$c$. Both prompt and non-prompt J$/\psi$ measurements indicate a significant suppression for $p_{\rm T} >$ 5 GeV/$c$, which becomes stronger with increasing collision centrality. The results are consistent with similar LHC measurements in the overlapping $p_{\rm T}$ intervals, and cover the kinematic region down to $p_{\rm T}$ = 1.5 GeV/$c$ at midrapidity, not accessible by other LHC experiments. The suppression of prompt J$/\psi$ in central and semicentral collisions exhibits a decreasing trend towards lower transverse momentum, described within uncertainties by models implementing J$/\psi$ production from recombination of c and $\overline{\rm c}$ quarks produced independently in different partonic scatterings. At high transverse momentum, transport models including quarkonium dissociation are able to describe the suppression for prompt J$/\psi$. For non-prompt J$/\psi$, the suppression predicted by models including both collisional and radiative processes for the computation of the beauty-quark energy loss inside the quark$-$gluon plasma is consistent with measurements within uncertainties.

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.

Transverse momentum ($p_{\rm T}$) spectra of charged particles at mid-pseudorapidity in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV measured with the ALICE apparatus at the Large Hadron Collider are reported. The kinematic range $0.15 < p_{\rm T} < 50$ GeV/$c$ and $|\eta| < 0.8$ is covered. Results are presented in nine classes of collision centrality in the 0-80% range. For comparison, a pp reference at the collision energy of $\sqrt{s}$ = 5.44 TeV is obtained by interpolating between existing \pp measurements at $\sqrt{s}$ = 5.02 and 7 TeV. The nuclear modification factors in central Xe-Xe collisions and Pb-Pb collisions at a similar center-of-mass energy of $\sqrt{s_{\rm NN}}$ = 5.02 TeV, and in addition at 2.76 TeV, at analogous ranges of charged particle multiplicity density $\left\langle\rm{d}N_{\rm ch}/\rm{d}\eta\right\rangle$ show a remarkable similarity at $p_{\rm T}> 10$ GeV/$c$. The comparison of the measured $R_{\rm AA}$ values in the two colliding systems could provide insight on the path length dependence of medium-induced parton energy loss. The centrality dependence of the ratio of the average transverse momentum $\left\langle p_{\rm{T}}\right\rangle$ in Xe-Xe collisions over Pb-Pb collision at $\sqrt{s}$ = 5.02 TeV is compared to hydrodynamical model calculations.

The production of non-prompt ${\rm D}^{0}$ mesons from beauty-hadron decays was measured at midrapidity ($\left| y \right| < 0.5$) in Pb-Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{\rm NN}}=5.02~\mathrm{TeV}$ with the ALICE experiment at the LHC. Their nuclear modification factor ($R_{\rm AA}$), measured for the first time down to $p_{\rm T}=1~\mathrm{GeV}/c$ in the $0-10$% and $30-50$% centrality classes, indicates a significant suppression, up to a factor of about three, for $p_{\rm T} > 5~\mathrm{GeV}/c$ in the $0-10$% central Pb-Pb collisions. The data are described by models that include both collisional and radiative processes in the calculation of beauty-quark energy loss in the quark-gluon plasma, and quark recombination in addition to fragmentation as a hadronisation mechanism. The ratio of the non-prompt to prompt ${\rm D}^{0}$-meson $R_{\rm AA}$ is larger than unity for $p_{\rm T} > 4~\mathrm{GeV}/c$ in the $0-10$% central Pb-Pb collisions, as predicted by models in which beauty quarks lose less energy than charm quarks in the quark-gluon plasma because of their larger mass.

The production yield and angular anisotropy of prompt ${D_s^+}$ mesons were measured as a function of transverse momentum ($p_{ T}$) in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{ NN}} = 5.02$ TeV collected with the ALICE detector at the LHC. ${D_s^+}$ mesons and their charge conjugates were reconstructed at midrapidity ($|y|<0.5$) from their hadronic decay channel ${D_s^+ \to \phi \pi^+}$, with ${\phi \to K^-K^+}$, in the $p_{ T}$ intervals $2<p_{ T}<50$ GeV/$c$ and $2<p_{ T}<36$ GeV/$c$ for the 0-10% and 30-50% centrality intervals. For $p_{ T}>10$ GeV/$c$, the measured ${D_s^+}$-meson nuclear modification factor $R_{ AA}$ is consistent with the one of non-strange D mesons within uncertainties, while at lower $p_{ T}$ a hint for a ${D_s^+}$-meson $R_{ AA}$ larger than that of non-strange D mesons is seen. The enhanced production of ${D_s^+}$ relative to non-strange D mesons is also studied by comparing the $p_{ T}$-dependent ${D_s^+/D^0}$ production yield ratios in Pb-Pb and in pp collisions. The ratio measured in Pb-Pb collisions is found to be on average higher than that in pp collisions in the interval $2<p_{ T} <8$ GeV/$c$ with a significance of 2.3$\sigma$ and 2.4$\sigma$ for the 0-10% and 30-50% centrality intervals. The azimuthal anisotropy coefficient $v_2$ of prompt ${D_s^+}$ mesons was measured in Pb-Pb collisions in the 30-50% centrality interval and is found to be compatible with that of non-strange D mesons. The main features of the measured $R_{ AA}$, ${D_s^+/D^0}$ ratio, and $v_2$ as a function of $p_{ T}$ are described by theoretical calculations of charm-quark transport in a hydrodynamically expanding quark-gluon plasma including hadronisation via charm-quark recombination with light quarks from the medium. The $p_{ T}$-integrated production yield of ${D_s^+}$ mesons is compatible with the prediction of the statistical hadronisation model.