This article presents the first measurement of the interaction between charm hadrons and nucleons. The two-particle momentum correlations of $\mathrm{pD^-}$ and $\mathrm{\overline{p}D}^+$ pairs are measured by the ALICE Collaboration in high-multiplicity pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$. The data are compatible with the Coulomb-only interaction hypothesis within (1.1-1.5)$\sigma$. The level of agreement slightly improves if an attractive nucleon(N)$\overline{\mathrm{D}}$ strong interaction is considered, in contrast to most model predictions which suggest an overall repulsive interaction. This measurement allows for the first time an estimation of the 68% confidence level interval for the isospin $\mathrm{I}=0$ inverse scattering length of the $\mathrm{N\overline{D}}$ state ${f_{0,~\mathrm{I}=0}^{-1} \in [-0.4,0.9]~\mathrm{fm^{-1}}}$, assuming negligible interaction for the isospin $\mathrm{I}=1$ channel.

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.

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 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 yields of non-prompt $\mathrm{D_s}^{+}$ mesons, namely $\mathrm{D_s}^{+}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_{\rm T}$) at midrapidity ($|y|<0.5$) in central and semi-central Pb$-$Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE experiment at the LHC. The $\mathrm{D_s}^{+}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel $\mathrm{D_s}^{+} \to \phi\pi^+$, with $\phi \to {\rm K}^-{\rm K}^+$, in the $4 < p_{\rm T}<36$ GeV/$c$ and $2 < p_{\rm T}<24$ GeV/$c$ intervals for the 0$-$10% and 30$-$50% centrality classes, respectively. The measured yields of non-prompt $\mathrm{D_s}^{+}$ mesons are compared to those of prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons by calculating the ratios of the production yields in Pb$-$Pb collisions and the nuclear modification factor $R_\mathrm{AA}$. The ratio between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and prompt $\mathrm{D_s}^{+}$ mesons, and that between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons in central Pb$-$Pb collisions are found to be on average higher than unity in the $4< p_{\rm T}<12$ GeV/$c$ interval with a statistical significance of about $1.6\,\sigma$ and $1.7\,\sigma$, respectively. The measured $R_\mathrm{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

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 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 elliptic flow ($v_2$) of $\rm D^0$ mesons from beauty-hadron decays (non-prompt $\rm D^0$) was measured in midcentral (30-50%) Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The $\rm D^0$ mesons were reconstructed at midrapidity ($|y|<0.8$) from their hadronic decay $\mathrm{D^0 \to K^-\pi^+}$, in the transverse momentum interval $2 < p_{\rm T} < 12$ GeV/$c$. The result indicates a positive $v_2$ for non-prompt $\rm D^0$ mesons with a significance of 2.7$\sigma$. The non-prompt $\rm D^0$-meson $v_2$ is lower than that of prompt non-strange D mesons with 3.2$\sigma$ significance in $2 < p_{\rm T} < 8$ GeV/$c$, and compatible with the $v_2$ of beauty-decay electrons. Theoretical calculations of beauty-quark transport in a hydrodynamically expanding medium describe the measurement within uncertainties.

Three-body nuclear forces play an important role in the structure of nuclei and hypernuclei and are also incorporated in models to describe the dynamics of dense baryonic matter, such as in neutron stars. So far, only indirect measurements anchored to the binding energies of nuclei can be used to constrain the three-nucleon force, and if hyperons are considered, the scarce data on hypernuclei impose only weak constraints on the three-body forces. In this work, we present the first direct measurement of the p$-$p$-$p and p$-$p$-\Lambda$ systems in terms of three-particle correlation functions carried out for pp collisions at $\sqrt{s} = 13$ TeV. Three-particle cumulants are extracted from the correlation functions by applying the Kubo formalism, where the three-particle interaction contribution to these correlations can be isolated after subtracting the known two-body interaction terms. A negative cumulant is found for the p$-$p$-$p system, hinting to the presence of a residual three-body effect while for p$-$p$-\Lambda$ the cumulant is consistent with zero. This measurement demonstrates the accessibility of three-baryon correlations at the LHC.

The production cross sections of $\mathrm {D^0}$, $\mathrm {D^+}$, and $\mathrm {\Lambda_{c}^{+}}$ hadrons originating from beauty-hadron decays (i.e. non-prompt) were measured for the first time at midrapidity in proton$-$lead (p$-$Pb) collisions at the center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV. Nuclear modification factors ($R_{\mathrm {pPb}}$) of non-prompt $\mathrm {D^0}$, $\mathrm {D^+}$, and $\mathrm {\Lambda_{c}^{+}}$ are calculated as a function of the transverse momentum ($p_{\mathrm T}$) to investigate the modification of the momentum spectra measured in p$-$Pb collisions with respect to those measured in proton$-$proton (pp) collisions at the same energy. The $R_{\mathrm {pPb}}$ measurements are compatible with unity and with the measurements in the prompt charm sector, and do not show a significant $p_{\mathrm T}$ dependence. The $p_{\mathrm T}$-integrated cross sections and $p_{\mathrm T}$-integrated $R_{\mathrm {pPb}}$ of non-prompt $\mathrm {D^0}$ and $\mathrm {D^+}$ mesons are also computed by extrapolating the visible cross sections down to $p_{\mathrm T}$ = 0. The non-prompt D-meson $R_{\mathrm {pPb}}$ integrated over $p_{\mathrm T}$ is compatible with unity and with model calculations implementing modification of the parton distribution functions of nucleons bound in nuclei with respect to free nucleons. The non-prompt $\mathrm {\Lambda_{c}^{+}/D^{0}}$ and $\mathrm{D^+/D^0}$ production ratios are computed to investigate hadronisation mechanisms of beauty quarks into mesons and baryons. The measured ratios as a function of $p_{\mathrm T}$ display a similar trend to that measured for charm hadrons in the same collision system.