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.

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.

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 first measurements of the production of muons and electrons from heavy-flavour hadron decays in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV, using the ALICE detector at the LHC, are reported. The measurement of the nuclear modification factor $R_{\rm AA}$ is performed as a function of transverse momentum $p_{\rm T}$ in several centrality classes at forward rapidity ($2.5 < y <4$) and midrapidity ($\vert y \vert < 0.8$) for muons and electrons from heavy-flavour hadron decays, respectively. A suppression by a factor up to about 2.5 compared to the binary-scaled pp reference is observed in central collisions at both central and forward rapidities. The $R_{\rm AA}$ of muons from heavy-flavour hadron decays is compared to previous measurements in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. When the nuclear modification factors are compared in the centrality classes 0-10% for Xe-Xe collisions and 10-20% for Pb-Pb collisions, which have similar charged-particle multiplicity density, a similar suppression, with $R_{\rm AA} \sim 0.4$ in the $p_{\rm T}$ interval $4 < p_{\rm T} < 8$ GeV/$c$, is observed. The comparison of the measured $R_{\rm AA}$ values in the two collision systems brings new insights on the properties of the quark-gluon plasma by investigating the system-size and path-length dependence of medium-induced parton energy loss. The results of muons and electrons from heavy-flavour hadron decays provide new constraints to model calculations.

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 inclusive J/$\psi$ production in Pb-Pb collisions at the center-of-mass energy per nucleon pair $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV, measured with the ALICE detector at the CERN LHC, is reported. The J/$\psi$ meson is reconstructed via the dimuon decay channel at forward rapidity ($2.5<y<4$) down to zero transverse momentum. The suppression of the J/$\psi$ yield in Pb-Pb collisions with respect to binary-scaled pp collisions is quantified by the nuclear modification factor ($R_{\rm{AA}}$). The $R_{\rm{AA}}$ at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV is presented and compared with previous measurements at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV as a function of the centrality of the collision, and of the J/$\psi$ transverse momentum and rapidity. The inclusive J/$\psi$ $R_{\rm{AA}}$ shows a suppression increasing toward higher $p_{\rm{T}}$, with a steeper dependence for central collisions. The modification of the J/$\psi$ average $p_{\rm{T}}$ and $p_{\rm{T}}^{2}$ is also studied. Comparisons with the results of models based on a transport equation and on statistical hadronization are also carried out.

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

The inclusive J/$\psi$ meson production in Pb-Pb collisions at a center-of-mass energy per nucleon-nucleon collision of $\sqrt{s_{\rm NN}}$ = 5.02 TeV at midrapidity ($|y|$ < 0.9) is reported by the ALICE Collaboration. The measurements are performed in the dielectron decay channel, as a function of event centrality and J/$\psi$ transverse momentum $p_{\rm T}$, down to $p_{\rm T}$ = 0 GeV/$c$. The J/$\psi$ mean transverse momentum $\langle p_{\rm T} \rangle$ and $r_{\rm AA}$ ratio, defined as $\langle p^{\rm 2}_{\rm T} \rangle_{\rm PbPb}/\langle p^{\rm 2}_{\rm T} \rangle_{\rm pp}$, are evaluated. Both observables show a centrality dependence decreasing towards central (head-on) collisions. The J/$\psi$ nuclear modification factor $R_{\rm AA}$ exhibits a strong $p_{\rm T}$ dependence with a large suppression at high $p_{\rm T}$ and an increase to unity for decreasing $p_{\rm T}$. When integrating over the measured momentum range $p_{\rm T}$ < 10 GeV/$c$, the J/$\psi$ $R_{\rm AA}$ shows a weak centrality dependence. Each measurement is compared with results at lower center-of-mass energies and with ALICE measurements at forward rapidity, as well as to theory calculations. All reported features of the J/$\psi$ production at low $p_{\rm T}$ are consistent with a dominant contribution to the J/$\psi$ yield originating from charm quark (re)combination.