Understanding the production mechanism of light (anti)nuclei is one of the key challenges of nuclear physics and has important consequences for astrophysics, since it provides an input for indirect dark-matter searches in space. In this paper, the latest results about the production of light (anti)nuclei in pp collisions at $\sqrt{s} = 13$ TeV are presented, focusing on the comparison with the predictions of coalescence and thermal models. For the first time, the coalescence parameters $B_2$ for deuterons and $B_3$ for helions are compared with parameter-free theoretical predictions that are directly constrained by the femtoscopic measurement of the source radius in the same event class. A fair description of the data with a Gaussian wave function is observed for both deuteron and helion, supporting the coalescence mechanism for the production of light (anti)nuclei in pp collisions. This method paves the way for future investigations of the internal structure of more complex nuclear clusters, including the hypertriton.

The production of prompt $D^{0}$, $D^{+}_{\rm s}$, and $\Lambda_{\rm c}^{+}$ hadrons, and their ratios, $D^{+}_{\rm s}$/$D^{0}$ and $\Lambda_{\rm c}^{+}$/$D^{0}$, are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity ($|y| <0.5$) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum ($p_{\rm T}$) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange $D^{+}_{\rm s}$/$D^{0}$ ratio indicates no significant multiplicity dependence, the baryon-to-meson $p_{\rm T}$-differential $\Lambda_{\rm c}^{+}$/$D^{0}$ ratio shows a multiplicity-dependent enhancement, with a significance of 5.3$\sigma$ for $1< p_{\rm T} < 12$ GeV/$c$, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The $\Lambda_{\rm c}^{+}$/$D^{0}$ ratios as a function of $p_{\rm T}$ present a similar shape and magnitude as the $\Lambda/K^{0}_{s}$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The $p_{\rm T}$-integrated ratios, extrapolated down to $p_{\rm T}$=0, do not show a significant dependence on multiplicity within the uncertainties.

First measurements of balance functions (BFs) of all combinations of identified charged hadron $(\pi,\rm K,\rm p)$ pairs in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV recorded by the ALICE detector are presented. The BF measurements are carried out as two-dimensional differential correlators versus the relative rapidity ($\Delta y$) and azimuthal angle ($\Delta\varphi$) of hadron pairs, and studied as a function of collision centrality. The $\Delta\varphi$ dependence of BFs is expected to be sensitive to the light quark diffusivity in the quark$-$gluon plasma. While the BF azimuthal widths of all pairs substantially decrease from peripheral to central collisions, the longitudinal widths exhibit mixed behaviors: BFs of $\pi\pi$ and cross-species pairs narrow significantly in more central collisions, whereas those of $\rm KK$ and $\rm pp$ are found to be independent of collision centrality. This dichotomy is qualitatively consistent with the presence of strong radial flow effects and the existence of two stages of quark production in relativistic heavy-ion collisions. Finally, the first measurements of the collision centrality evolution of BF integrals are presented, with the observation that charge balancing fractions are nearly independent of collision centrality in Pb$-$Pb collisions. Overall, the results presented provide new and challenging constraints for theoretical models of hadron production and transport in relativistic heavy-ion collisions.

The measurement of the production of deuterons, tritons and $^{3}\mathrm{He}$ and their antiparticles in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV is presented in this article. The measurements are carried out at midrapidity ($|y| < $ 0.5) as a function of collision centrality using the ALICE detector. The $p_{\rm T}$-integrated yields, the coalescence parameters and the ratios to protons and antiprotons are reported and compared with nucleosynthesis models. The comparison of these results in different collision systems at different centre-of-mass collision energies reveals a suppression of nucleus production in small systems. In the Statistical Hadronisation Model framework, this can be explained by a small correlation volume where the baryon number is conserved, as already shown in previous fluctuation analyses. However, a different size of the correlation volume is required to describe the proton yields in the same data sets. The coalescence model can describe this suppression by the fact that the wave functions of the nuclei are large and the fireball size starts to become comparable and even much smaller than the actual nucleus at low multiplicities.

The study of the production of nuclei and antinuclei 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 protons, deuterons and $^3$He and their charge conjugates at midrapidity is studied as a function of the charged-particle multiplicity in inelastic pp collisions at $\sqrt{s}=5.02$ TeV using the ALICE detector. Within the uncertainties, the yields of nuclei in pp collisions at $\sqrt{s}=5.02$ TeV are compatible with those in pp collisions at different energies and to those in p-Pb collisions when compared at similar multiplicities. The measurements are compared with the expectations of coalescence and Statistical Hadronisation Models. The results suggest a common formation mechanism behind the production of light nuclei in hadronic interactions and confirm that they do not depend on the collision energy but on the number of produced particles.

The first measurements of elliptic flow of $\pi^\pm$, ${\rm K}^\pm$, p+$\overline{\rm p}$, ${\rm K_{S}^0}$, $\Lambda$+$\overline{\Lambda}$, $\phi$, $\Xi^-$+$\Xi^+$, and $\Omega^-$+$\Omega^+$ using multiparticle cumulants in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented. Results obtained with two- ($v_2\{2\}$) and four-particle cumulants ($v_2\{4\}$) are shown as a function of transverse momentum, $p_{\rm T}$, for various collision centrality intervals. Combining the data for both $v_2\{2\}$ and $v_2\{4\}$ also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark-gluon plasma. The characteristic features observed in previous $p_{\rm T}$-differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low $p_{\rm T}$ and the approximate scaling with the number of constituent quarks at intermediate $p_{\rm T}$, are similarly present in the four-particle correlations and the combinations of $v_2\{2\}$ and $v_2\{4\}$. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide $p_{\rm T}$ range.

The study of the azimuthal anisotropy of inclusive muons produced in p-Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV, using the ALICE detector at the LHC is reported. The measurement of the second-order Fourier coefficient of the particle azimuthal distribution, $v_2$, is performed as a function of transverse momentum $p_{\rm T}$ in the 0-20% high-multiplicity interval at both forward ($2.03 < y_{\rm CMS} < 3.53$) and backward ($-4.46 < y_{\rm CMS} < -2.96$) rapidities over a wide $p_{\rm T}$ range, $0.5 < p_{\rm T} < 10$ GeV/$c$, in which a dominant contribution of muons from heavy-flavour hadron decays is expected at $p_{\rm T} > 2$ GeV/$c$. The $v_2$ coefficient of inclusive muons is extracted using two different techniques, namely two-particle cumulants, used for the first time for heavy-flavour measurements, and forward-central two-particle correlations. Both techniques give compatible results. A positive $v_2$ is measured at both forward and backward rapidities with a significance larger than $4.7\sigma$ and $7.6\sigma$, respectively, in the interval $2 < p_{\rm T} < 6$ GeV/$c$. Comparisons with previous measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, and with AMPT and CGC-based theoretical calculations are discussed. The findings impose new constraints on the theoretical interpretations of the origin of the collective behaviour in small collision systems.

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 production of $\pi^\pm$, ${\rm K}^\pm$, and $(\overline{\rm p})$p is measured in pp collisions at $\sqrt{s}=13$ TeV in different topological regions. Particle transverse momentum ($p_{\rm T}$) spectra are measured in the ``toward'', ``transverse'', and ``away'' angular regions defined with respect to the direction of the leading particle in the event. While the toward and away regions contain the fragmentation products of the near-side and away-side jets, respectively, the transverse region is dominated by particles from the Underlying Event (UE). The relative transverse activity classifier, $R_{\rm T}=N_{\rm T}/\langle N_{\rm T}\rangle$, is used to group events according to their UE activity, where $N_{\rm T}$ is the measured charged-particle multiplicity per event in the transverse region and $\langle N_{\rm T}\rangle$ is the mean value over all the analysed events. The first measurements of identified particle $p_{\rm T}$ spectra as a function of $R_{\rm T}$ in the three topological regions are reported. The yield of high transverse momentum particles relative to the $R_{\rm T}$-integrated measurement decreases with increasing $R_{\rm T}$ in both the toward and away regions, indicating that the softer UE dominates particle production as $R_{\rm T}$ increases and validating that $R_{\rm T}$ can be used to control the magnitude of the UE. Conversely, the spectral shapes in the transverse region harden significantly with increasing $R_{\rm T}$. This hardening follows a mass ordering, being more significant for heavier particles. The $p_{\rm T}$-differential particle ratios $({\rm p+\overline{p}})/(\pi^+ +\pi^-)$ and $({\rm K^+ +K^-})/(\pi^+ +\pi^-)$ in the low UE limit $(R_{\rm T}\rightarrow 0)$ approach expectations from Monte Carlo generators such as PYTHIA 8 with Monash 2013 tune and EPOS LHC, where the jet-fragmentation models have been tuned to reproduce ${\rm e^+ e^-}$ results.

The fractions of non-prompt (i.e. originating from beauty-hadron decays) D$^0$ and D$^+$ mesons with respect to the inclusive yield are measured as a function of the charged-particle multiplicity in proton$-$proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV with the ALICE detector at the LHC. The results are reported in intervals of transverse momentum ($p_{\rm T}$) and integrated in the range $1 < p_{\rm T} < 24$ GeV/$c$. The fraction of non-prompt D$^0$ and D$^+$ mesons is found to increase slightly as a function of $p_{\rm T}$ in all the measured multiplicity intervals, while no significant dependence on the charged-particle multiplicity is observed. In order to investigate the production and hadronisation mechanisms of charm and beauty quarks, the results are compared to PYTHIA 8 as well as EPOS 3 and EPOS 4 Monte Carlo simulations, and to calculations based on the colour glass condensate including three-pomeron fusion.