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

The production of prompt $\rm \Lambda_{\rm c}^{+}$ baryons at midrapidity ($|y|<0.5$) was measured in central (0-10%) and mid-central (30-50%) Pb-Pb collisions at the center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector. The results are more precise, more differential in centrality, and reach much lower transverse momentum ($p_{\rm T}=1$ GeV/$c$) with respect to previous measurements performed by the ALICE, STAR, and CMS Collaborations in nucleus-nucleus collisions, allowing for an extrapolation down to $p_{\rm T}=0$. The $p_{\rm T}$-differential $\rm \Lambda_{\rm c}^{+}$/D$^0$ production ratio is enhanced with respect to the pp measurement for $4<p_{\rm T}<8$ GeV/$c$ by 3.7 standard deviations ($\sigma$), while the $p_{\rm T}$-integrated ratios are compatible within 1$\sigma$. The observed trend is similar to that observed in the strange sector for the $\Lambda/$K$^0_{\rm S}$ ratio. Model calculations including coalescence or statistical hadronization for charm-hadron formation are compared with the data.

Measurements of the production of muons from heavy-flavour hadron decays in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = $5.02$ and $2.76$ TeV using the ALICE detector at the LHC are reported. The nuclear modification factor $R_{\rm AA}$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV is measured at forward rapidity ($2.5 < y <4$) as a function of transverse momentum $p_{\rm T}$ in central, semi-central, and peripheral collisions over a wide $p_{\rm T}$ interval, $3 < p_{\rm T} < 20$ GeV/$c$, in which muons from beauty-hadron decays are expected to take over from charm as the dominant source at high $p_{\rm T}$ ($p_{\rm T} > 7$ GeV/$c$). The $R_{\rm AA}$ shows an increase of the suppression of the yields of muons from heavy-flavour hadron decays with increasing centrality. A suppression by a factor of about three is observed in the $10\%$ most central collisions. The $R_{\rm AA}$ at $\sqrt{s_{\rm NN}}$ = 5.02 is similar to that at 2.76 TeV. The precise $R_{\rm AA}$ measurements have the potential to distinguish between model predictions implementing different mechanisms of parton energy loss in the high-density medium formed in heavy-ion collisions. They place important constraints for the understanding of the heavy-quark interaction with the hot and dense QCD medium.

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 J/$\psi$ production has been studied with the ALICE detector in p-Pb collisions at the nucleon-nucleon center of mass energy $\sqrt{s_{\rm NN}}$ = 5.02 TeV at the CERN LHC. The measurement is performed in the center of mass rapidity domains $2.03<y_{\rm cms}<3.53$ and $-4.46<y_{\rm cms}<-2.96$, down to zero transverse momentum, studying the $\mu^+\mu^-$ decay mode. In this paper, the J/$\psi$ production cross section and the nuclear modification factor $R_{\rm pPb}$ for the rapidities under study are presented. While at forward rapidity, corresponding to the proton direction, a suppression of the J/$\psi$ yield with respect to binary-scaled pp collisions is observed, in the backward region no suppression is present. The ratio of the forward and backward yields is also measured differentially in rapidity and transverse momentum. Theoretical predictions based on nuclear shadowing, as well as on models including, in addition, a contribution from partonic energy loss, are in fair agreement with the experimental results.