The first measurement of the $\Upsilon(1{\rm S})$ elliptic flow coefficient ($v_2$) is performed at forward rapidity (2.5 $<$ $y$ $<$ 4) in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The results are obtained with the scalar product method and are reported as a function of transverse momentum ($p_{\rm{T}}$) up to 15 GeV/$c$ in the 5-60% centrality interval. The measured $\Upsilon(1{\rm S})$ $v_2$ is consistent with zero and with the small positive values predicted by transport models within uncertainties. The $v_2$ coefficient in 2 $<$ $p_{\rm T}$ $<$ 15 GeV/$c$ is lower than that of inclusive J/$\psi$ mesons in the same $p_{\rm{T}}$ interval by 2.6 standard deviations. These results, combined with earlier suppression measurements, are in agreement with a scenario in which the $\Upsilon$(1S) production in Pb-Pb collisions at LHC energies is dominated by dissociation limited to the early stage of the collision whereas in the J/$\psi$ case there is substantial experimental evidence of an additional regeneration component.

Measurements of anisotropic flow coefficients ($v_n$) and their cross-correlations using two- and multi-particle cumulant methods are reported in collisions of pp at $\sqrt{s} = 13$ TeV, p-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV, Xe-Xe at $\sqrt{s_{_{\rm NN}}} = 5.44$ TeV, and Pb-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV recorded with the ALICE detector. The multiplicity dependence of $v_n$ is studied in a very wide range from 20 to 3000 particles produced in the mid-rapidity region $|\eta|<0.8$ for the transverse momentum range $0.2 < p_{\rm T} < 3.0$ GeV/$c$. An ordering of the coefficients $v_2 > v_3 > v_4$ is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak $v_2$ multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, $v_{2}$ measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between $v_n^2$ and $v_m^2$, evaluated with the symmetric cumulants SC$(m,n)$ is observed to be positive at all multiplicities for $v_2$ and $v_4$, while for $v_2$ and $v_3$ it is negative and changes sign for multiplicities below 100, which may indicate a different $v_n$ fluctuation pattern in this multiplicity range. The observed long-range multi-particle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by IP-Glasma+MUSIC+UrQMD model calculations, and hence provide new insights into the understanding of collective effects in small collision systems.

The measurement of the azimuthal-correlation function of prompt D mesons with charged particles in pp collisions at $\sqrt{s}$ = 5.02 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC is reported. The D$^{\rm 0}$, D$^{\rm +}$, and D$^{\rm *+}$ mesons, together with their charge conjugates, were reconstructed at midrapidity in the transverse momentum interval 3 < $p_{\rm T}$ < 24 GeV/c and correlated with charged particles having $p_{\rm T}$ > 0.3 GeV/c and pseudorapidity $|\eta| <$ 0.8. The properties of the correlation peaks appearing in the near- and away-side regions (for $\Delta \varphi \approx$ 0 and $\Delta \varphi \approx \pi$, respectively) were extracted via a fit to the azimuthal correlation functions. The shape of the correlation functions and the near- and away-side peak features are found to be consistent in pp and p-Pb collisions, showing no modifications due to nuclear effects within uncertainties. The results are compared with predictions from Monte Carlo simulations performed with the PYTHIA, POWHEG+PYTHIA, HERWIG, and EPOS 3 event generators.

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.

Inclusive J/$\psi$ yields and average transverse momenta in p-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 8.16 TeV are measured as a function of the charged-particle pseudorapidity density with ALICE. The J/$\psi$ mesons are reconstructed at forward $(2.03 < y_{\rm cms} < 3.53)$ and backward ($-4.46 < y_{\rm cms} < -2.96$) center-of-mass rapidity in their dimuon decay channel while the charged-particle pseudorapidity density is measured around midrapidity. The J/$\psi$ yields at forward and backward rapidity normalized to their respective average values increase with the normalized charged-particle pseudorapidity density, the former showing a weaker increase than the latter. The normalized average transverse momenta at forward and backward rapidity manifest a steady increase from low to high charged-particle pseudorapidity density with a saturation beyond the average value.

The production of J/$\psi$ is measured at midrapidity ($|y|<0.9$) in proton-proton collisions at $\sqrt{s}$ = 5.02 and 13 TeV, through the dielectron decay channel, using the ALICE detector at the Large Hadron Collider. The data sets used for the analyses correspond to integrated luminosities of $\mathcal{L}_{\rm int}$ = 19.4 $\pm$ 0.4 nb$^{-1}$ and $\mathcal{L}_{\rm int}$ = 32.2 $\pm$ 0.5 nb$^{-1}$ at $\sqrt{s}$ = 5.02 and 13 TeV, respectively. The fraction of non-prompt J/$\psi$ mesons, i.e. those originating from the decay of beauty hadrons, is measured down to a transverse momentum $p_{\rm T}$ = 2 GeV/$c$ (1 GeV/$c$) at $\sqrt{s}$ = 5.02 TeV (13 TeV). The $p_{\rm T}$ and rapidity ($y$) differential cross sections, as well as the corresponding values integrated over $p_{\rm T}$ and $y$, are carried out separately for prompt and non-prompt J/$\psi$ mesons. The results are compared with measurements from other experiments and theoretical calculations based on quantum chromodynamics (QCD). The shapes of the $p_{\rm T}$ and $y$ distributions of beauty quarks predicted by state-of-the-art perturbative QCD models are used to extrapolate an estimate of the $\rm b\overline{b}$ pair cross section at midrapidity and in the total phase space. The total $\rm b\overline{b}$ cross sections are found to be $\sigma_{\rm b \overline{\rm b}} = 541 \pm 45 (\rm stat.) \pm 69 (\rm syst.)_{-12}^{+10} (\rm extr.)~{\rm \mu b}$ and $\sigma_{\rm b \overline{\rm b}}~=~218 \pm 37 (\rm stat.) \pm 31 (\rm syst.)_{-9.1}^{+8.2} (\rm extr.)~{\rm \mu b}$ at $\sqrt{s}$ = 13 and 5.02 TeV, respectively. The value obtained from the combination of ALICE and LHCb measurements in pp collisions at $\sqrt{s}$ = 13 TeV is also provided.

The polarization of the $\Lambda$ and $\overline\Lambda$ hyperons along the beam ($z$) direction, $P_{\rm z}$, has been measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02TeV recorded with ALICE at the Large Hadron Collider (LHC). The main contribution to $P_{\rm z}$ comes from elliptic flow induced vorticity and can be characterized by the second Fourier sine coefficient $P_{\rm z,s2} = \langle P_{\rm z} \sin(2\varphi - 2 \Psi_{\rm 2}) \rangle$, where $\varphi$ is the hyperon azimuthal emission angle, and $\Psi_{\rm 2}$ is the elliptic flow plane angle. We report the measurement of $P_{\rm z,\,{\rm s2}}$ for different collision centralities, and in the 30-50% centrality interval as a function of the hyperon transverse momentum and rapidity. The $P_{\rm z,\,{\rm s2}}$ is positive similarly as measured by the STAR Collaboration in Au-Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV, with somewhat smaller amplitude in the semi-central collisions. This is the first experimental evidence of a non-zero hyperon $P_{\rm z}$ in Pb-Pb collisions at the LHC. The comparison of the measured $P_{\rm z,\,{\rm s2}}$ with the hydrodynamic model calculations shows sensitivity to the competing contributions from thermal and the recently found shear induced vorticity, as well as to whether the polarization is acquired at the quark-gluon plasma or the hadronic phase.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

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 transverse momentum ($p_{\rm T}$) spectra and elliptic flow coefficient ($v_2$) of deuterons and anti-deuterons at mid-rapidity ($|y|<0.5$) are measured with the ALICE detector at the LHC in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV. The measurement of the $p_{\rm T}$ spectra of (anti-)deuterons is done up to 8 GeV/$c$ in 0-10% centrality class and up to GeV/$c$ in 10-20% and 20-40% centrality classes. The $v_2$ is measured in the $0.8 <p_{\rm T} <5$ GeV/$c$ interval and in six different centrality intervals (0-5%, 5-10%, 10-20%, 20-30%, 30-40% and 40-50%) using the scalar product technique. Measured $\pi^{\pm}$, K$^{\pm}$ and p+$\overline{\mathrm{p}}$ transverse-momentum spectra and $v_2$ are used to predict the deuteron $p_{\rm T}$ spectra and $v_2$ within the Blast-Wave model. The predictions are able to reproduce the $v_2$ coefficient in the measured $p_{\rm T}$ range and the transverse-momentum spectra for $p_{\rm T}>1.8$ GeV/$c$ within the experimental uncertainties. The measurement of the coalescence parameter $B_2$ is performed, showing a $p_{\rm T}$ dependence in contrast with the simplest coalescence model, which fails to reproduce also the measured $v_2$ coefficient. In addition, the coalescence parameter $B_2$ and the elliptic flow coefficient in the 20-40% centrality interval are compared with the AMPT model which is able, in its version without string melting, to reproduce the measured $v_2$($p_{\rm T}$) and the $B_2$($p_{\rm T}$) trend.