The measurements of the inclusive J/$\psi$ yield at midrapidity ($\left | y \right | < 0.9$) and forward rapidity (2.5 $< y <$ 4) in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ALICE detector at the LHC are reported. The inclusive J/$\psi$ production yields and nuclear modification factors, $R_{\rm AA}$, are measured as a function of the collision centrality, J/$\psi$ transverse momentum ($p_{\rm T}$), and rapidity. The J/$\psi$ average transverse momentum and squared transverse momentum ($\langle p_{\mathrm{T}}\rangle$ and $\langle p_{\mathrm{T}}^{\mathrm{2}}\rangle$) are evaluated as a function of the centrality at midrapidity. Compared to the previous ALICE publications, here the entire Pb$-$Pb collisions dataset collected during the LHC Run 2 is used, which improves the precision of the measurements and extends the $p_{\rm T}$ coverage. The $p_{\rm T}$-integrated $R_{\rm AA}$ shows a hint of an increasing trend towards unity from semicentral to central collisions at midrapidity, while it is flat at forward rapidity. The $p_{\rm T}$-differential $R_{\rm AA}$ shows a strong suppression at high $p_{\rm T}$ with less suppression at low $p_{\rm T}$ where it reaches a larger value at midrapidity compared to forward rapidity. The ratio of the $p_{\rm T}$-integrated yields of J/$\psi$ to those of D$^{0}$ mesons is reported for the first time for the central and semicentral event classes at midrapidity. Model calculations implementing charmonium production via the coalescence of charm quarks and antiquarks during the fireball evolution (transport models) or in a statistical approach with thermal weights are in good agreement with the data at low $p_{\rm T}$. At higher $p_{\rm T}$, the data are well described by transport models and a model based on energy loss in the strongly-interacting medium produced in nuclear collisions at the LHC.

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 azimuthal anisotropy of particles associated with jets (jet particles) at midrapidity is measured for the first time in p-Pb and Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV down to transverse momentum ($p_{\rm T}$) of 0.5 GeV/$c$ and 2 GeV/$c$, respectively, with ALICE. The results obtained in p-Pb collisions are based on a novel three-particle correlation technique. The azimuthal anisotropy coefficient $v_2$ in high-multiplicity p-Pb collisions is positive, with a significance reaching 6.8$\sigma$ at low $p_{\rm T}$, and its magnitude is smaller than in semicentral Pb-Pb collisions. In contrast to the measurements in Pb-Pb collisions, the $v_2$ coefficient is also found independent of $p_{\rm T}$ within uncertainties. Comparisons with the inclusive charged-particle $v_2$ and with AMPT calculations are discussed. The predictions suggest that parton interactions play an important role in generating a non-zero jet-particle $v_2$ in p-Pb collisions, even though they overestimate the reported measurement. These observations shed new insights on the understanding of the origin of the collective behaviour of jet particles in small systems such as p-Pb collisions, and provide significant stringent new constraints to models.

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.

Two-particle transverse momentum differential correlators, recently measured in Pb--Pb collisions at energies available at the CERN Large Hadron Collider (LHC), provide an additional tool to gain insights into particle production mechanisms and infer transport properties, such as the ratio of shear viscosity to entropy density, of the medium created in Pb-Pb collisions. The longitudinal long-range correlations and the large azimuthal anisotropy measured at low transverse momenta in small collision systems, namely pp and p-Pb, at LHC energies resemble manifestations of collective behaviour. This suggests that locally equilibrated matter may be produced in these small collision systems, similar to what is observed in Pb-Pb collisions. In this work, the same two-particle transverse momentum differential correlators are exploited in pp and p-Pb collisions at $\sqrt{s} = 7$ TeV and $\sqrt{s_{\rm NN}} = 5.02$ TeV, respectively, to seek evidence for viscous effects. Specifically, the strength and shape of the correlators are studied as a function of the produced particle multiplicity to identify evidence for longitudinal broadening that might reveal the presence of viscous effects in these smaller systems. The measured correlators and their evolution from pp and p--Pb to Pb--Pb collisions are additionally compared to predictions from Monte Carlo event generators, and the potential presence of viscous effects is discussed.

The ALICE Collaboration reports a differential measurement of inclusive jet suppression using pp and Pb$-$Pb collision data at a center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}} = 5.02$ TeV. Charged-particle jets are reconstructed using the anti-$k_{\rm T}$ algorithm with resolution parameters $R =$ 0.2, 0.3, 0.4, 0.5, and 0.6 in pp collisions and $R =$ 0.2, 0.4, 0.6 in central (0$-$10%), semi-central (30$-$50%), and peripheral (60$-$80%) Pb$-$Pb collisions. A novel approach based on machine learning is employed to mitigate the influence of jet background. This enables measurements of inclusive jet suppression in new regions of phase space, including down to the lowest jet $p_{\rm T} \geq 40$ GeV/$c$ at $R = 0.6$ in central Pb$-$Pb collisions. This is an important step for discriminating different models of jet quenching in the quark-gluon plasma. The transverse momentum spectra, nuclear modification factors, derived cross section, and nuclear modification factor ratios for different jet resolution parameters of charged-particle jets are presented and compared to model predictions. A mild dependence of the nuclear modification factor ratios on collision centrality and resolution parameter is observed. The results are compared to a variety of jet-quenching models with varying levels of agreement.

High-energy nuclear collisions create a quark-gluon plasma, whose initial condition and subsequent expansion vary from event to event, impacting the distribution of the event-wise average transverse momentum ($P([p_{\mathrm{T}}])$). Distinguishing between contributions from fluctuations in the size of the nuclear overlap area (geometrical component) and other sources at fixed size (intrinsic component) presents a challenge. Here, these two components are distinguished by measuring the mean, variance, and skewness of $P([p_{\mathrm{T}}])$ in $^{208}$Pb+$^{208}$Pb and $^{129}$Xe+$^{129}$Xe collisions at $\sqrt{s_{{\mathrm{NN}}}} = 5.02$ and 5.44 TeV, respectively, using the ATLAS detector at the LHC. All observables show distinct changes in behavior in ultra-central collisions, where the geometrical variations are suppressed as the overlap area reaches its maximum. These results demonstrate a new technique to disentangle geometrical and intrinsic fluctuations, enabling constraints on initial condition and properties of the quark-gluon plasma, such as the speed of sound.

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.

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