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.

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.

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.

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 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.

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.

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.