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.

The differential invariant yield as a function of transverse momentum ($p_\mathrm{T}$) of electrons from semileptonic heavy-flavour hadron decays was measured at midrapidity in central (0-10%), semi-central (30-50%) and peripheral (60-80%) lead-lead (Pb-Pb) collisions at $\sqrt{s_{\mathrm{NN}}}=5.02\text{ TeV}$ in the $p_{\mathrm{T}}$ intervals 0.5-26 GeV/$c$ (0-10% and 30-50%) and 0.5-10 GeV/$c$ (60-80%). The production cross section in proton-proton (pp) collisions at $\sqrt{s}=5.02$ TeV was measured as well in $0.5<p_\mathrm{T}<10$ GeV/$c$ and it lies close to the upper band of perturbative QCD calculation uncertainties up to $p_\mathrm{T}=5$ GeV/$c$ and close to the mean value for larger $p_\mathrm{T}$. The modification of the electron yield with respect to what is expected for an incoherent superposition of nucleon-nucleon collisions is evaluated by measuring the nuclear modification factor $R_{\mathrm{AA}}$. The measurement of the $R_{\mathrm{AA}}$ in different centrality classes allows in-medium energy loss of charm and beauty quarks to be investigated. The $R_{\mathrm{AA}}$ shows a suppression with respect to unity at intermediate $p_\mathrm{T}$, which increases while moving towards more central collisions. Moreover, the measured $R_{\mathrm{AA}}$ is sensitive to the modification of the parton distribution functions (PDF) in nuclei, like nuclear shadowing, which causes a suppression of the heavy-quark production at low $p_\mathrm{T}$ in heavy-ion collisions at LHC.

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.

Inclusive J/$\psi$ production is studied in Xe-Xe interactions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}= 5.44$ TeV, using the ALICE detector at the CERN LHC. The J/$\psi$ meson is reconstructed via its decay into a muon pair, in the centre-of-mass rapidity interval $2.5<y<4$ and down to zero transverse momentum. In this Letter, the nuclear modification factors $R_{\rm AA}$ for inclusive J/$\psi$, measured in the centrality range 0-90% as well as in the centrality intervals 0-20% and 20-90% are presented. The $R_{\rm AA}$ values are compared to previously published results for Pb-Pb collisions at $\sqrt{s_{\rm NN}}= 5.02$ TeV and to the calculation of a transport model. A good agreement is found between Xe-Xe and Pb-Pb results as well as between data and the model.

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 $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.