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.

In this Letter, the first evidence of the ${}^4_{\bar{\Lambda}}\overline{\mathrm{He}}$ antihypernucleus is presented, along with the first measurement at the LHC of the production of (anti)hypernuclei with mass number $A=4$, specifically (anti)${}^4_{\Lambda}\mathrm{H}$ and (anti)${}^4_{\Lambda}\mathrm{He}$. In addition, the antiparticle-to-particle ratios for both hypernuclei (${}^4_{\bar{\Lambda}}\overline{\mathrm{H}}$ / ${}^4_{\Lambda}\mathrm{H}$~and ${}^4_{\bar{\Lambda}}\overline{\mathrm{He}}$ / ${}^4_{\Lambda}\mathrm{He}$) are shown, which are sensitive to the baryochemical potential of the strongly-interacting matter created in heavy-ion collisions. The results are obtained from a data sample of central Pb--Pb collisions, collected during the 2018 LHC data-taking at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} = $ 5.02 TeV. The yields measured for the average of the charge-conjugated states are found to be $[0.78 \; \pm \; 0.19 \; \mathrm{(stat.)} \; \pm \; 0.17 \; \mathrm{(syst.)}] \times 10^{-6}$ for the (anti)${}^4_{\Lambda}\mathrm{H}$ and $[1.08 \; \pm \; 0.34 \; \mathrm{(stat.)} \; \pm \; 0.20 \; \mathrm{(syst.)}] \times 10^{-6}$ for the (anti)${}^4_{\Lambda}\mathrm{He}$, and the measured antiparticle-to-particle ratios are in agreement with unity. The presence of (anti)${}^4_{\Lambda}\mathrm{H}$ and (anti)${}^4_{\Lambda}\mathrm{He}$ excited states is expected to strongly enhance the production yield of these hypernuclei. The yield values exhibit a combined deviation of 3.3$\sigma$ from the theoretical ground-state-only expectation, while the inclusion of the excited states in the calculations leads to an agreement within 0.6$\sigma$ with the present measurements. Additionally, the measured (anti)${}^4_{\Lambda}\mathrm{H}$ and (anti)${}^4_{\Lambda}\mathrm{He}$ masses are compatible with the world-average values within the uncertainties.

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.

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 interactions of kaons (K) and antikaons ($\mathrm{\overline{K}}$) with few nucleons (N) were studied so far using kaonic atom data and measurements of kaon production and interaction yields in nuclei. Some details of the three-body KNN and $\mathrm{\overline{K}}$NN dynamics are still not well understood, mainly due to the overlap with multi-nucleon interactions in nuclei. An alternative method to probe the dynamics of three-body systems with kaons is to study the final state interaction within triplet of particles emitted in pp collisions at the Large Hadron Collider, which are free from effects due to the presence of bound nucleons. This Letter reports the first femtoscopic study of p$-$p$-$K$^+$ and p$-$p$-$K$^-$ correlations measured in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV by the ALICE Collaboration. The analysis shows that the measured p$-$p$-$K$^+$ and p$-$p$-$K$^-$ correlation functions can be interpreted in terms of pairwise interactions in the triplets, indicating that the dynamics of such systems is dominated by the two-body interactions without significant contributions from three-body effects or bound states.

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 nuclei and antinuclei production has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high-energy hadronic collisions. The first measurement of the production of ${\rm ^{3}_{\Lambda}\rm H}$ in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV is presented in this Letter. Its production yield measured in the rapidity interval $-1 < y < 0$ for the 40% highest multiplicity p-Pb collisions is ${\rm d} N /{\rm d} y =[\mathrm{6.3 \pm 1.8 (stat.) \pm 1.2 (syst.) ] \times 10^{-7}}$. The measurement is compared with the expectations of statistical hadronisation and coalescence models, which describe the nucleosynthesis in hadronic collisions. These two models predict very different yields of the hypertriton in charged particle multiplicity environments relevant to small collision systems such as p-Pb and therefore the measurement of ${\rm d} N /{\rm d} y$ is crucial to distinguish between them. The precision of this measurement leads to the exclusion with a significance larger than 6.9$\sigma$ of some configurations of the statistical hadronization model, thus constraining the theory behind the production of loosely bound states at hadron colliders.

Results on the transverse spherocity dependence of light-flavor particle production ($\pi$, K, p, $\phi$, ${\rm K^{*0}}$, ${\rm K}^{0}_{\rm{S}}$, $\Lambda$, $\Xi$) at midrapidity in high-multiplicity pp collisions at $\sqrt{s} = 13$ TeV were obtained with the ALICE apparatus. The transverse spherocity estimator ($S_{{\rm O}}^{{\it p}_{\rm T}=1}$) categorizes events by their azimuthal topology. Utilizing narrow selections on $S_{\text{O}}^{{\it p}_{\rm T}=1}$, it is possible to contrast particle production in collisions dominated by many soft initial interactions with that observed in collisions dominated by one or more hard scatterings. Results are reported for two multiplicity estimators covering different pseudorapidity regions. The $S_{{\rm O}}^{{\it p}_{\rm T}=1}$ estimator is found to effectively constrain the hardness of the events when the midrapidity ($\left | \eta \right |< 0.8$) estimator is used. The production rates of strange particles are found to be slightly higher for soft isotropic topologies, and severely suppressed in hard jet-like topologies. These effects are more pronounced for hadrons with larger mass and strangeness content, and observed when the topological selection is done within a narrow multiplicity interval. This demonstrates that an important aspect of the universal scaling of strangeness enhancement with final-state multiplicity is that high-multiplicity collisions are dominated by soft, isotropic processes. On the contrary, strangeness production in events with jet-like processes is significantly reduced. The results presented in this article are compared with several QCD-inspired Monte Carlo event generators. Models that incorporate a two-component phenomenology, either through mechanisms accounting for string density, or thermal production, are able to describe the observed strangeness enhancement as a function of $S_{{\rm O}}^{{\it p}_{\rm T}=1}$.

The production yields of the $\Sigma(1385)^{\pm}$ and $\Xi(1530)^{0}$ resonances are measured in pp collisions at $\sqrt{s}=13$ TeV with ALICE. The measurements are performed as a function of the charged-particle multiplicity $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$, which is related to the energy density produced in the collision. The results include transverse momentum ($p_{\rm T}$) distributions, $p_{\rm T}$-integrated yields, mean transverse momenta of $\Sigma(1385)^{\pm}$ and $\Xi(1530)^{0}$, as well as ratios of the $p_{\rm T}$-integrated resonance yields relative to yields of other hadron species. The $\Sigma(1385)^{\pm}/\pi^{\pm}$ and $\Xi(1530)^{0}/\pi^{\pm}$ yield ratios are consistent with the trend of the enhancement of strangeness production from low to high multiplicity pp collisions, which was previously observed for strange and multi-strange baryons. The yield ratio between the measured resonances and the long-lived baryons with the same strangeness content exhibits a hint of a mild increasing trend at low multiplicity, despite too large uncertainties to exclude the flat behaviour. The results are compared with predictions from models such as EPOS-LHC and PYTHIA 8 with Rope shoving. The latter provides the best description of the multiplicity dependence of the $\Sigma(1385)^{\pm}$ and $\Xi(1530)^{0}$ production in pp collisions at $\sqrt{s}=13$ TeV.