The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

The first measurements of anisotropic flow coefficients $v_{\rm{n}}$ for mid-rapidity charged particles in Xe-Xe collisions at $\sqrt{s_{\rm{NN}}} = 5.44$ TeV are presented. Comparing these measurements to those from Pb-Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV, $v_{2}$ is found to be suppressed for mid-central collisions at the same centrality, and enhanced for central collisions. The values of $v_{3}$ are generally larger in Xe-Xe than in Pb-Pb at a given centrality. These observations are consistent with expectations from hydrodynamic predictions. When both $v_{2}$ and $v_{3}$ are divided by their corresponding eccentricities for a variety of initial state models, they generally scale with transverse density when comparing Xe-Xe and Pb-Pb, with some deviations observed in central Xe-Xe and Pb-Pb collisions. These results assist in placing strong constraints on both the initial state geometry and medium response for relativistic heavy-ion collisions.

This letter presents the first measurement of jet mass in Pb-Pb and p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV and 5.02 TeV, respectively. Both the jet energy and the jet mass are expected to be sensitive to jet quenching in the hot Quantum Chromodynamics (QCD) matter created in nuclear collisions at collider energies. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm and resolution parameter $R = 0.4$. The jets are measured in the pseudorapidity range $|\eta_{\rm{jet}}|<0.5$ and in three intervals of transverse momentum between 60 GeV/$c$ and 120 GeV/$c$. The measurement of the jet mass in central Pb-Pb collisions is compared to the jet mass as measured in p-Pb reference collisions, to vacuum event generators, and to models including jet quenching. It is observed that the jet mass in central Pb-Pb collisions is consistent within uncertainties with p-Pb reference measurements. Furthermore, the measured jet mass in Pb-Pb collisions is not reproduced by the quenching models considered in this letter and is found to be consistent with PYTHIA expectations within systematic uncertainties.

The production of the charm-strange baryon $\Xi_{\rm c}^0$ is measured for the first time at the LHC via its semileptonic decay into e$^+\Xi^-\nu_{\rm e}$ in pp collisions at $\sqrt{s}=7$ TeV with the ALICE detector. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval 1 $<$ $p_{\rm T}$ $<$ 8 GeV/$c$ at mid-rapidity, $|y|$ $<$ 0.5. The transverse momentum dependence of the $\Xi_{\rm c}^0$ baryon production relative to the D$^0$ meson production is compared to predictions of event generators with various tunes of the hadronisation mechanism, which are found to underestimate the measured cross-section ratio.

$\phi$ meson measurements provide insight into strangeness production, which is one of the key observables for the hot medium formed in high-energy heavy-ion collisions. ALICE measured $\phi$ production through its decay in muon pairs in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 2.76 TeV in the intermediate transverse momentum range $2 < p_{\rm T} < 5$ GeV/$c$ and in the rapidity interval $2.5<y<4$. The $\phi$ yield was measured as a function of the transverse momentum and collision centrality. The nuclear modification factor was obtained as a function of the average number of participating nucleons. Results were compared with the ones obtained via the kaon decay channel in the same $p_{\rm T}$ range at midrapidity. The values of the nuclear modification factor in the two rapidity regions are in agreement within uncertainties.

We report on the measurement of the inclusive J/$\psi$ polarization parameters in pp collisions at a center of mass energy $\sqrt{s} = 8$ TeV with the ALICE detector at the LHC. The analysis is based on a data sample corresponding to an integrated luminosity of 1.23 pb$^{-1}$. J/$\psi$ resonances are reconstructed in their di-muon decay channel in the rapidity interval $2.5 < y < 4.0$ and over the transverse-momentum interval $2 < p_{\rm T} < 15$ GeV/$c$. The three polarization parameters ($\lambda_\theta$, $\lambda_\varphi$, $\lambda_{\theta\varphi}$) are measured as a function of $p_{\rm T}$ both in the helicity and Collins-Soper reference frames. The measured J/$\psi$ polarization parameters are found to be compatible with zero within uncertainties, contrary to expectations from all available predictions. The results are compared with the measurement in pp collisions at $\sqrt{s} = 7$ TeV.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.

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.

We report measurements of the production of prompt D$^0$, D$^+$, D$^{*+}$ and D$^+_{\rm s}$ mesons in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV, in the centrality classes 0-10%, 30-50% and 60-80%. The D-meson production yields are measured at mid-rapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm T}$). The $p_{\rm T}$ intervals covered in central collisions are: $1<p_{\rm T}<50$ Gev/$c$ for D$^0$, $2<p_{\rm T}<50$ GeV/$c$ for D$^+$, $3<p_{\rm T}<50$ GeV/$c$ for D$^{*+}$, and $4<p_{\rm T}<16$ GeV/$c$ for D$^+_{\rm s}$ mesons. The nuclear modification factors ($R_{\rm AA}$) for non-strange D mesons (D$^0$, D$^+$, D$^{*+}$) show minimum values of about 0.2 for $p_{\rm T}$ = 6-10 GeV/$c$ in the most central collisions and are compatible within uncertainties with those measured at $\sqrt{s_{\rm NN}}=2.76$ TeV. For D$^+_{\rm s}$ mesons, the values of $R_{\rm AA}$ are larger than those of non-strange D mesons, but compatible within uncertainties. In central collisions the average $R_{\rm AA}$ of non-strange D mesons is compatible with that of charged particles for $p_{\rm T} > 8$ GeV/$c$, while it is larger at lower $p_{\rm T}$. The nuclear modification factors for strange and non-strange D mesons are also compared to theoretical models with different implementations of in-medium energy loss.

Measurements of anisotropic flow coefficients with two- and multi-particle cumulants for inclusive charged particles in Pb-Pb collisions at $\sqrt{{\textit s}_\text{NN}} = 5.02$ and 2.76 TeV are reported in the pseudorapidity range $|\eta| < 0.8$ and transverse momentum $0.2 < p_\text{T} < 50$ GeV/$c$. The full data sample collected by the ALICE detector in 2015 (2010), corresponding to an integrated luminosity of 12.7 (2.0) $\mu$b$^{-1}$ in the centrality range 0-80%, is analysed. Flow coefficients up to the sixth flow harmonic ($v_6$) are reported and a detailed comparison among results at the two energies is carried out. The $p_\text{T}$ dependence of anisotropic flow coefficients and its evolution with respect to centrality and harmonic number $n$ are investigated. An approximate power-law scaling of the form $v_n(p_\text{T}) \sim p_\text{T}^{n/3}$ is observed for all flow harmonics at low $p_\text{T}$ ($0.2 < p_\text{T} < 3$ GeV/$c$). At the same time, the ratios $v_n/v_m^{n/m}$ are observed to be essentially independent of $p_\text{T}$ for most centralities up to about $p_\text{T} = 10$ GeV/$c$. Analysing the differences among higher-order cumulants of elliptic flow ($v_2$), which have different sensitivities to flow fluctuations, a measurement of the standardised skewness of the event-by-event $v_2$ distribution $P(v_2)$ is reported and constraints on its higher moments are provided. The Elliptic Power distribution is used to parametrise $P(v_2)$, extracting its parameters from fits to cumulants. The measurements are compared to different model predictions in order to discriminate among initial-state models and to constrain the temperature dependence of the shear viscosity to entropy-density ratio.

We present a measurement of azimuthal correlations between inclusive J/$\psi$ and charged hadrons in p-Pb collisions recorded with the ALICE detector at the CERN LHC. The J/$\psi$ are reconstructed at forward (p-going, 2.03 $<$ y $<$ 3.53) and backward (Pb-going, $-$4.46 $<$ y $<$ $-$2.96) rapidity via their $\mu^+\mu^-$ decay channel, while the charged hadrons are reconstructed at mid-rapidity ($|\eta|$ $<$ 1.8). The correlations are expressed in terms of associated charged-hadron yields per J/$\psi$ trigger. A rapidity gap of at least 1.5 units is required between the trigger J/$\psi$ and the associated charged hadrons. Possible correlations due to collective effects are assessed by subtracting the associated per-trigger yields in the low-multiplicity collisions from those in the high-multiplicity collisions. After the subtraction, we observe a strong indication of remaining symmetric structures at $\Delta\varphi$ $\approx$ 0 and $\Delta\varphi$ $\approx$ $\pi$, similar to those previously found in two-particle correlations at middle and forward rapidity. The corresponding second-order Fourier coefficient ($v_2$) in the transverse momentum interval between 3 and 6 GeV/$c$ is found to be positive with a significance of about 5$\sigma$. The obtained results are similar to the J/$\psi$ $v_2$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV, suggesting a common mechanism at the origin of the J/$\psi$ $v_2$.

Azimuthally-differential femtoscopic measurements, being sensitive to spatio-temporal characteristics of the source as well as to the collective velocity fields at freeze out, provide very important information on the nature and dynamics of the system evolution. While the HBT radii oscillations relative to the second harmonic event plane measured recently reflect mostly the spatial geometry of the source, model studies have shown that the HBT radii oscillations relative to the third harmonic event plane are predominantly defined by the velocity fields. In this Letter, we present the first results on azimuthally-differential pion femtoscopy relative to the third harmonic event plane as a function of the pion pair transverse momentum $k_{\rm T}$ for different collision centralities in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. We find that the $R_{\rm side}$ and $R_{\rm out}$ radii, which characterize the pion source size in the directions perpendicular and parallel to the pion transverse momentum, oscillate in phase relative to the third harmonic event plane, similar to the results from 3+1D hydrodynamical calculations. The observed radii oscillations unambiguously signal a collective expansion and anisotropy in the velocity fields. A comparison of the measured radii oscillations with the Blast-Wave model calculations indicate that the initial state triangularity is washed-out at freeze out.

Ultrarelativistic heavy ion collisions recreate in the laboratory the thermodynamical conditions prevailing in the early universe up to 10$^{-6}$ seconds, thereby allowing the study of the quark-gluon plasma (QGP), a state of quantum chromodynamics (QCD) matter with deconfined partons. The top quark, the heaviest elementary particle known, is accessible in nucleus-nucleus collisions at the CERN LHC, and constitutes a novel probe of the QGP. Here, we report the first-ever evidence for the production of top quarks in nucleus-nucleus collisions, using lead-lead collision data at a nucleon-nucleon centre-of-mass energy of 5.02 TeV recorded by the CMS experiment. Two methods are used to measure the cross section for top quark pair production ($\sigma_\mathrm{t\bar{t}}$) via the decay into charged leptons (electrons or muons) and bottom quarks. One method relies on the leptonic information alone, and the second one exploits, in addition, the presence of bottom quarks. The measured cross sections, $\sigma_\mathrm{t\bar{t}} = $ 2.54 $^{+0.84}_{-0.74}$ and 2.03 $^{+0.71}_{-0.64}$ $\mu$b, respectively, are compatible with expectations from scaled proton-proton data and QCD predictions.

The production of the W$^\pm$ bosons measured in p$-$Pb collisions at a centre-of-mass energy per nucleon$-$nucleon collision $\sqrt{s_{NN}} = 8.16$ TeV and Pb$-$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV with ALICE at the LHC is presented. The W$^\pm$ bosons are measured via their muonic decay channel, with the muon reconstructed in the pseudorapidity region $-4 < \eta^\mu_{\rm lab} < -2.5$ with transverse momentum $p_{\rm T}^\mu > 10$ GeV/$c$. While in Pb$-$Pb collisions the measurements are performed in the forward ($2.5 < y^\mu_{\rm cms} < 4$) rapidity region, in p$-$Pb collisions, where the centre-of-mass frame is boosted with respect to the laboratory frame, the measurements are performed in the backward ($-4.46 < y^\mu_{\rm cms} < -2.96$) and forward ($2.03 < y^\mu_{\rm cms} < 3.53$) rapidity regions. The W$^{-}$ and W$^{+}$ production cross sections, lepton-charge asymmetry, and nuclear modification factors are evaluated as a function of the muon rapidity. In order to study the production as a function of the p$-$Pb collision centrality, the production cross sections of the W$^{-}$ and W$^{+}$ bosons are combined and normalised to the average number of binary nucleon$-$nucleon collision $\langle N_\mathrm{coll} \rangle$. In Pb$-$Pb collisions, the same measurements are presented as a function of the collision centrality. Study of the binary scaling of the W$^\pm$-boson cross sections in p$-$Pb and Pb$-$Pb collisions is also reported. The results are compared with perturbative QCD (pQCD) calculations, with and without nuclear modifications of the Parton Distribution Functions (PDFs), as well as with available data at the LHC. Significant deviations from the theory expectations are found in the two collision systems, indicating that the measurements can provide additional constraints for the determination of nuclear PDF (nPDFs) and in particular of the light-quark distributions.

This article reports measurements characterizing the Underlying Event (UE) associated with hard scatterings at midrapidity in pp collisions at $\sqrt{s}=13$ TeV. The hard scatterings are identified by the leading particle, the charged particle with the highest transverse momentum ($p_{\rm T}^{\rm leading}$) in the event. Charged-particle number and summed transverse-momentum densities are measured in different azimuthal regions defined with respect to the leading particle direction: Toward, Transverse, and Away. The Toward and Away regions contain the fragmentation products of the hard scatterings in addition to the UE contribution, whereas particles in the Transverse region are expected to originate predominantly from the UE. The study is performed as a function of $p_{\rm T}^{\rm leading}$ with three different $p_{\rm T}$ thresholds for the associated particles, $p_{\rm T}^{\rm min} >$ 0.15, 0.5, and 1.0 GeV/$c$. The charged-particle density in the Transverse region rises steeply for low values of $p_{\rm T}^{\rm leading}$ and reaches a plateau. The results confirm the trend that the charged-particle density in the Transverse region shows a stronger increase with $\sqrt{s}$ than the inclusive charged-particle density at midrapidity. The UE activity is increased by approximately 20% when going from 7 to 13 TeV. The plateau in the Transverse region ($5 < p_{\rm T}^{\rm leading} < ~ 40$ GeV/$c$ ) is further characterized by the probability distribution of its charged-particle multiplicity normalized to its average value (relative transverse activity, $R_{T}$) and the mean transverse momentum as a function of $R_{T}$. Experimental results are compared to model calculations using PYTHIA 8 and EPOS LHC. The overall agreement between models and data is within 30%. These measurements provide new insights on the interplay between hard scatterings and the associated UE in pp collisions.

The polarization of inclusive J/$\psi$ and $\Upsilon(1{\rm S})$ produced in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV at the LHC is measured with the ALICE detector. The study is carried out by reconstructing the quarkonium through its decay to muon pairs in the rapidity region $2.5<y<4$ and measuring the polar and azimuthal angular distributions of the muons. The polarization parameters $\lambda_{\theta}$, $\lambda_{\phi}$ and $\lambda_{\theta\phi}$ are measured in the helicity and Collins-Soper reference frames, in the transverse momentum interval $2<p_{\rm T}<10$ GeV/$c$ and $p_{\rm T}<15$ GeV/$c$ for the J/$\psi$ and $\Upsilon(1{\rm S})$, respectively. The polarization parameters for the J/$\psi$ are found to be compatible with zero, within a maximum of about two standard deviations at low $p_{\rm T}$, for both reference frames and over the whole $p_{\rm T}$ range. The values are compared with the corresponding results obtained for pp collisions at $\sqrt{s}=7$ and 8 TeV in a similar kinematic region by the ALICE and LHCb experiments. Although with much larger uncertainties, the polarization parameters for $\Upsilon(1{\rm S})$ production in Pb-Pb collisions are also consistent with zero.

The first measurement of e$^+$e$^-$ pair production at mid-rapidity ($|\eta_{{\rm e}}|$ $<$ 0.8) in pp collisions at $\sqrt{s} = 7$ TeV with ALICE at the LHC is presented. The dielectron production is studied as a function of the invariant mass ($m_{\rm ee}$ $<$ 3.3 GeV/$c^{2}$), the pair transverse momentum ($p_{\rm T,ee}$ $<$ 8 GeV/$c$), and the pair transverse impact parameter (DCA$_{{\rm ee}}$), i.e., the average distance of closest approach of the reconstructed electron and positron tracks to the collision vertex, normalised to its resolution. The results are compared with the expectations from a cocktail of known hadronic sources and are well described when PYTHIA is used to generate the heavy-flavour contributions. In the low-mass region (0.14 $<$ $m_{\rm ee}$ $<$ 1.1 GeV/$c^{2}$), prompt and non-prompt e$^+$e$^-$ sources can be separated via the DCA$_{\rm ee}$. In the intermediate-mass region (1.1 $<$ $m_{\rm ee}$ $<$ 2.7 GeV/$c^{2}$), a double-differential fit to the data in $m_{\rm ee}$ and $p_{\rm T,ee}$ and a fit of the DCA$_{\rm ee}$ distribution allow the total ${\rm c\overline c}$ and ${\rm b\overline b}$ cross sections to be extracted. Two different event generators, PYTHIA and POWHEG, can reproduce the shape of the two-dimensional $m_{\rm ee}$ and $p_{\rm T,ee}$ spectra, as well as the shape of the DCA$_{\rm ee}$ distribution, reasonably well. However, differences in the ${\rm c\overline c}$ and ${\rm b\overline b}$ cross sections are observed when using the generators to extrapolate to full phase space. Finally, the ratio of inclusive to decay photons is studied via the measurement of virtual direct photons in the transverse-momentum range 1 $<$ $p_{\rm T}$ $<$ 8 GeV/$c$. This is found to be unity within the statistical and systematic uncertainties and consistent with expectations from next-to-leading order perturbative quantum chromodynamic calculations.

The $p_{\rm T}$-differential production cross section of prompt $\Lambda_{\rm c}^+$ charmed baryons was measured with the ALICE detector at the Large Hadron Collider (LHC) in pp collisions at $\sqrt{s} = 7$ TeV and in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV at midrapidity. The $\Lambda_{\rm c}^+$ and ${\overline{\Lambda}}_{\rm c}^-$ were reconstructed in the hadronic decay modes $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K^-}\pi^+$, $\Lambda_{\rm c}^+\rightarrow {\rm p}{\rm K_{\rm S}^0}$ and in the semileptonic channel $\Lambda_{\rm c}^+\rightarrow {\rm e^+}\nu_{\rm e}\Lambda$ (and charge conjugates). The measured values of the $\Lambda_{\rm c}^+/{\rm D_0}$ ratio, which is sensitive to the c-quark hadronisation mechanism, and in particular to the production of baryons, are presented and are larger than those measured previously in different colliding systems, centre-of-mass energies, rapidity and $p_{\rm T}$ intervals, where the $\Lambda_{\rm c}^+$ production process may differ. The results are compared with the expectations obtained from perturbative Quantum Chromodynamics calculations and Monte Carlo event generators. Neither perturbative QCD calculations nor Monte Carlo models reproduce the data, indicating that the fragmentation of heavy-flavour baryons is not well understood. The first measurement at the LHC of the $\Lambda_{\rm c}^+$ nuclear modification factor, $R_{\rm pPb}$, is also presented. The $R_{\rm pPb}$ is found to be consistent with unity and with that of D mesons within the uncertainties, and consistent with a theoretical calculation that includes cold nuclear matter effects and a calculation that includes charm quark interactions with a deconfined medium.

Measurement of Z-boson production in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=8.16$ TeV and Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV is reported. It is performed in the dimuon decay channel, through the detection of muons with pseudorapidity $-4 < \eta_{\mu} < -2.5$ and transverse momentum $p_{\rm T}^{\mu} > 20$ GeV/$c$ in the laboratory frame. The invariant yield and nuclear modification factor are measured for opposite-sign dimuons with invariant mass $60 < m^{\mu\mu} < 120$ GeV$c^2$ and rapidity $2.5 < y_{cms}^{\mu\mu} < 4$. They are presented as a function of rapidity and, for the Pb-Pb collisions, of centrality as well. The results are compared with theoretical calculations, both with and without nuclear modifications to the Parton Distribution Functions (PDFs). In p-Pb collisions the center-of-mass frame is boosted with respect to the laboratory frame, and the measurements cover the backward ($-4.46< y_{cms}^{\mu\mu}<-2.96$) and forward ($2.03< y_{cms}^{\mu\mu}<3.53$) rapidity regions. For the p-Pb collisions, the results are consistent within experimental and theoretical uncertainties with calculations that include both free-nucleon and nuclear-modified PDFs. For the Pb-Pb collisions, a $3.4\sigma$ deviation is seen in the integrated yield between the data and calculations based on the free-nucleon PDFs, while good agreement is found once nuclear modifications are considered.

This Letter presents the first direct investigation of the p-$\Sigma^{0}$ interaction, using the femtoscopy technique in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV measured by the ALICE detector. The $\Sigma^{0}$ is reconstructed via the decay channel to $\Lambda \gamma$, and the subsequent decay of $\Lambda$ to p$\pi^-$. The photon is detected via the conversion in material to e$^{+}$e$^{-}$ pairs exploiting the unique capability of the ALICE detector to measure electrons at low transverse momenta. The measured p-$\Sigma^{0}$ correlation indicates a shallow strong interaction. The comparison of the data to several theoretical predictions obtained employing the $Correlation~Analysis~Tool~using~the~Schr\"odinger~Equation$ (CATS) and the Lednick\'y-Lyuboshits approach shows that the current experimental precision does not yet allow to discriminate between different models, as it is the case for the available scattering and hypernuclei data. Nevertheless, the p-$\Sigma^{0}$ correlation function is found to be sensitive to the strong interaction, and driven by the interplay of the different spin and isospin channels. This pioneering study demonstrates the feasibility of a femtoscopic measurement in the p-$\Sigma^{0}$ channel and with the expected larger data samples in LHC Run 3 and Run 4, the p-$\Sigma^{0}$ interaction will be constrained with high precision.

Deuteron production in high-energy collisions is sensitive to the space-time evolution of the collision system, and is typically described by a coalescence mechanism. For the first time, we present results on jet-associated deuteron production in \pp collisions at $\sqrt{s}=13$ TeV, providing an opportunity to test the established picture for deuteron production in events with a hard scattering. Using a trigger particle with high transverse-momentum ($p_{\rm T}>5$ GeV/$c$) as a proxy for the presence of a jet at midrapidity, we observe a measurable population of deuterons being produced around the jet proxy. The associated deuteron yield measured in a narrow angular range around the trigger particle differs by 2.4-4.8 standard deviations from the uncorrelated background. The data are described by PYTHIA model calculations featuring baryon coalescence.

The first evidence of spin alignment of vector mesons ($K^{*0}$ and $\phi$) in heavy-ion collisions at the Large Hadron Collider (LHC) is reported. The spin density matrix element $\rho_{00}$ is measured at midrapidity ($|y| <$ 0.5) in Pb-Pb collisions at a center-of-mass energy ($\sqrt{s_{\rm NN}}$) of 2.76 TeV with the ALICE detector. $\rho_{00}$ values are found to be less than 1/3 (1/3 implies no spin alignment) at low transverse momentum ($p_{\rm T} <$ 2 GeV/$c$) for $K^{*0}$ and $\phi$ at a level of 3$\sigma$ and 2$\sigma$, respectively. No significant spin alignment is observed for the $K^0_S$ meson (spin = 0) in Pb-Pb collisions and for the vector mesons in $pp$ collisions. The measured spin alignment is unexpectedly large but qualitatively consistent with the expectation from models which attribute it to a polarization of quarks in the presence of angular momentum in heavy-ion collisions and a subsequent hadronization by the process of recombination.

The invariant differential cross section of inclusive $\omega(782)$ meson production at midrapidity ($|y|<0.5$) in pp collisions at $\sqrt{s}$ = 7 TeV was measured with the ALICE detector at the LHC over a transverse momentum range of 2 < $p_{\rm{T}}$ < 17 GeV/$c$. The $\omega$ meson was reconstructed via its $\omega\rightarrow\pi^+\pi^-\pi^0$ decay channel. The measured $\omega$ production cross section is compared to various calculations: PYTHIA 8.2 Monash 2013 describes the data, while PYTHIA 8.2 Tune 4C overestimates the data by about 50%. A recent NLO calculation, which includes a model describing the fragmentation of the whole vector-meson nonet, describes the data within uncertainties below 6 GeV/$c$, while it overestimates the data by up to 50% for higher $p_{\rm{T}}$. The $\omega/\pi^0$ ratio is in agreement with previous measurements at lower collision energies and the PYTHIA calculations. In addition, the measurement is compatible with transverse mass scaling within the measured $p_{\rm{T}}$ range and the ratio is constant with $C^{\omega/\pi^{0}}$ = 0.67 $\pm$ 0.03 (stat) $\pm$ 0.04 (sys) above a transverse momentum of 2.5 GeV/$c$.

The second- and third-order azimuthal anisotropy Fourier harmonics of charged particles produced in pPb collisions, at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV, are studied over a wide range of event multiplicities. Multiparticle correlations are used to isolate global properties stemming from the collision overlap geometry. The second-order "elliptic" harmonic moment is obtained with high precision through four-, six-, and eight-particle correlations and, for the first time, the third-order "triangular" harmonic moment is studied using four-particle correlations. A sample of peripheral PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV that covers a similar range of event multiplicities as the pPb results is also analyzed. Model calculations of initial-state fluctuations in pPb and PbPb collisions can be directly compared to the high precision experimental results. This work provides new insight into the fluctuation-driven origin of the $v_3$ coefficients in pPb and PbPb collisions, and into the dominating overall collision geometry in PbPb collisions at the earliest stages of heavy ion interactions.