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.

We report on the production of inclusive $\Upsilon$(1S) and $\Upsilon$(2S) in p-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV at the LHC. The measurement is performed with the ALICE detector at backward ($-4.46< y_{{\rm cms}}<-2.96$) and forward ($2.03< y_{{\rm cms}}<3.53$) rapidity down to zero transverse momentum. The production cross sections of the $\Upsilon$(1S) and $\Upsilon$(2S) are presented, as well as the nuclear modification factor and the ratio of the forward to backward yields of $\Upsilon$(1S). A suppression of the inclusive $\Upsilon$(1S) yield in p-Pb collisions with respect to the yield from pp collisions scaled by the number of binary nucleon-nucleon collisions is observed at forward rapidity but not at backward rapidity. The results are compared to theoretical model calculations including nuclear shadowing or partonic energy loss effects.

The production of Z$^0$ bosons at large rapidities in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV is reported. Z$^0$ candidates are reconstructed in the dimuon decay channel (${\rm Z}^0 \rightarrow \mu^+\mu^-$), based on muons selected with pseudo-rapidity $-4.0<\eta<-2.5$ and $p_{\rm T}>20$ GeV/$c$. The invariant yield and the nuclear modification factor, $R_{\rm AA}$, are presented as a function of rapidity and collision centrality. The value of $R_{\rm AA}$ for the 0-20% central Pb-Pb collisions is $0.67 \pm 0.11 \, \mbox{(stat.)} \, \pm 0.03 \, \mbox{(syst.)} \, \pm 0.06 \, \mbox{(corr. syst.)}$, exhibiting a deviation of $2.6 \sigma$ from unity. The results are well-described by calculations that include nuclear modifications of the parton distribution functions, while the predictions using vacuum PDFs deviate from data by $2.3\sigma$ in the 0-90% centrality class and by $3\sigma$ in the 0-20% central collisions.

We report the measured transverse momentum ($p_{\rm T}$) spectra of primary charged particles from pp, p-Pb and Pb-Pb collisions at a center-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$ TeV in the kinematic range of $0.15<p_{\rm T}<50$ GeV/$c$ and $|\eta|< 0.8$. A significant improvement of systematic uncertainties motivated the reanalysis of data in pp and Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV, as well as in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, which is also presented. Spectra from Pb-Pb collisions are presented in nine centrality intervals and are compared to a reference spectrum from pp collisions scaled by the number of binary nucleon-nucleon collisions. For central collisions, the $p_{\rm T}$ spectra are suppressed by more than a factor of 7 around 6-7 GeV/$c$ with a significant reduction in suppression towards higher momenta up to 30 GeV/$c$. The nuclear modification factor $R_{\rm pPb}$, constructed from the pp and p-Pb spectra measured at the same collision energy, is consistent with unity above 8 GeV/$c$. While the spectra in both pp and Pb-Pb collisions are substantially harder at $\sqrt{s_{\rm NN}} = 5.02$ TeV compared to 2.76 TeV, the nuclear modification factors show no significant collision energy dependence. The obtained results should provide further constraints on the parton energy loss calculations to determine the transport properties of the hot and dense QCD matter.

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.

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.