Collective behavior has been observed in high-energy heavy-ion collisions for several decades. Collectivity is driven by the high particle multiplicities that are produced in these collisions. At the Large Hadron Collider (LHC), features of collectivity have also been seen in high-multiplicity proton-proton collisions that can attain particle multiplicities comparable to peripheral Pb-Pb collisions. One of the possible signatures of collective behavior is the decrease of femtoscopic radii extracted from pion and kaon pairs emitted from high-multiplicity collisions with increasing pair transverse momentum. This decrease can be described in terms of an approximate transverse mass scaling. In the present work, femtoscopic analyses are carried out by the ALICE collaboration on charged pion and kaon pairs produced in pp collisions at $\sqrt{s}=13$ TeV from the LHC to study possible collectivity in pp collisions. The event-shape analysis method based on transverse sphericity is used to select for spherical versus jet-like events, and the effects of this selection on the femtoscopic radii for both charged pion and kaon pairs are studied. This is the first time this selection method has been applied to charged kaon pairs. An approximate transverse-mass scaling of the radii is found in all multiplicity ranges studied when the difference in the Lorentz boost for pions and kaons is taken into account. This observation does not support the hypothesis of collective expansion of hot and dense matter that should only occur in high-multiplicity events. A possible alternate explanation of the present results is based on a scenario of common emission conditions for pions and kaons in pp collisions for the multiplicity ranges studied.

Measurements of charged-particle production in pp, p$-$Pb, and Pb$-$Pb collisions in the toward, away, and transverse regions with the ALICE detector are discussed. These regions are defined event-by-event relative to the azimuthal direction of the charged trigger particle, which is the reconstructed particle with the largest transverse momentum ($p_{\mathrm{T}}^{\rm trig}$) in the range $8<p_{\mathrm{T}}^{\rm trig}<15$ GeV$/c$. The toward and away regions contain the primary and recoil jets, respectively; both regions are accompanied by the underlying event (UE). In contrast, the transverse region perpendicular to the direction of the trigger particle is dominated by the so-called UE dynamics, and includes also contributions from initial- and final-state radiation. The relative transverse activity classifier, $R_{\mathrm{T}}=N_{\mathrm{ch}}^{\mathrm{T}}/\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$, is used to group events according to their UE activity, where $N_{\mathrm{ch}}^{\mathrm{T}}$ is the charged-particle multiplicity per event in the transverse region and $\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$ is the mean value over the whole analysed sample. The energy dependence of the $R_{\mathrm{T}}$ distributions in pp collisions at $\sqrt{s}=2.76$, 5.02, 7, and 13 TeV is reported, exploring the Koba-Nielsen-Olesen (KNO) scaling properties of the multiplicity distributions. The first measurements of charged-particle $p_{\rm T}$ spectra as a function of $R_{\mathrm{T}}$ in the three azimuthal regions in pp, p$-$Pb, and Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV are also reported. Data are compared with predictions obtained from the event generators PYTHIA 8 and EPOS LHC. This set of measurements is expected to contribute to the understanding of the origin of collective-like effects in small collision systems (pp and p$-$Pb).

The Chiral Magnetic Wave (CMW) phenomenon is essential to provide insights into the strong interaction in QCD, the properties of the quark-gluon plasma, and the topological characteristics of the early universe, offering a deeper understanding of fundamental physics in high-energy collisions. Measurements of the charge-dependent anisotropic flow coefficients are studied in Pb-Pb collisions at center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\mathrm{NN}}}=$ 5.02 TeV to probe the CMW. In particular, the slope of the normalized difference in elliptic ($v_{2}$) and triangular ($v_{3}$) flow coefficients of positively and negatively charged particles as a function of their event-wise normalized number difference, is reported for inclusive and identified particles. The slope $r_{3}^{\rm Norm}$ is found to be larger than zero and to have a magnitude similar to $r_{2}^{\rm Norm}$, thus pointing to a large background contribution for these measurements. Furthermore, $r_{2}^{\rm Norm}$ can be described by a blast wave model calculation that incorporates local charge conservation. In addition, using the event shape engineering technique yields a fraction of CMW ($f_{\rm CMW}$) contribution to this measurement which is compatible with zero. This measurement provides the very first upper limit for $f_{\rm CMW}$, and in the 10-60% centrality interval it is found to be 26% (38%) at 95% (99.7%) confidence level.

The transverse momentum ($p_{\rm T}$) and centrality dependence of the nuclear modification factor $R_{\rm AA}$ of prompt and non-prompt J$/\psi$, the latter originating from the weak decays of beauty hadrons, have been measured by the ALICE collaboration in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The measurements are carried out through the ${\rm e}^{+}{\rm e}^{-}$ decay channel at midrapidity ($|y| < 0.9$) in the transverse momentum region $1.5 < p_{\rm T} < 10$ GeV/$c$. Both prompt and non-prompt J$/\psi$ measurements indicate a significant suppression for $p_{\rm T} >$ 5 GeV/$c$, which becomes stronger with increasing collision centrality. The results are consistent with similar LHC measurements in the overlapping $p_{\rm T}$ intervals, and cover the kinematic region down to $p_{\rm T}$ = 1.5 GeV/$c$ at midrapidity, not accessible by other LHC experiments. The suppression of prompt J$/\psi$ in central and semicentral collisions exhibits a decreasing trend towards lower transverse momentum, described within uncertainties by models implementing J$/\psi$ production from recombination of c and $\overline{\rm c}$ quarks produced independently in different partonic scatterings. At high transverse momentum, transport models including quarkonium dissociation are able to describe the suppression for prompt J$/\psi$. For non-prompt J$/\psi$, the suppression predicted by models including both collisional and radiative processes for the computation of the beauty-quark energy loss inside the quark$-$gluon plasma is consistent with measurements within uncertainties.

The production of K$^*$(892)$^\pm$ meson resonance is measured at midrapidity ($|y|<0.5$) in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV using the ALICE detector at the LHC. The resonance is reconstructed via its hadronic decay channel K$^*$(892)$^\pm \rightarrow \rm{K^0_S \pi^\pm}$. The transverse momentum distributions are obtained for various centrality intervals in the $p_{\rm T}$ range of 0.4-16 GeV/$c$. The reported measurements of integrated yields, mean transverse momenta, and particle yield ratios are consistent with previous ALICE measurements for K$^*$(892)$^0$. The $p_{\rm T}$-integrated yield ratio 2K$^*$(892)$^\pm$/($\rm{K^+ + K^-}$) in central Pb-Pb collisions shows a significant suppression (9.3$\sigma$) relative to pp collisions. Thermal model calculations overpredict the particle yield ratio. Although both simulations consider the hadronic phase, only HRG-PCE accurately represents the measurements, whereas MUSIC+SMASH tends to overpredict them. These observations, along with the kinetic freeze-out temperatures extracted from the yields of light-flavored hadrons using the HRG-PCE model, indicate a finite hadronic phase lifetime, which increases towards central collisions. The $p_{\rm T}$-differential yield ratios 2K$^*$(892)$^\pm$/($\rm{K^+ + K^-}$) and 2K$^*$(892)$^\pm$/($\rm{\pi^+ + \pi^-}$) are suppressed by up to a factor of five at $p_{\rm T}<2$ GeV/$c$ in central Pb-Pb collisions compared to pp collisions at $\sqrt{s} =$ 5.02 TeV. Both particle ratios and are qualitatively consistent with expectations for rescattering effects in the hadronic phase. The nuclear modification factor shows a smooth evolution with centrality and is below unity at $p_{\rm T}>8$ GeV/$c$, consistent with measurements for other light-flavored hadrons. The smallest values are observed in most central collisions, indicating larger energy loss of partons traversing the dense medium.

The production cross sections of ${\rm D^0}$ and $\Lambda^+_{\rm c}$ hadrons originating from beauty-hadron decays (i.e. non-prompt) were measured for the first time at midrapidity ($|y|<0.5$) by the ALICE Collaboration in proton-proton collisions at a center-of-mass energy $\sqrt{s}=13$ TeV. They are described within uncertainties by perturbative QCD calculations employing the fragmentation fractions of beauty quarks to baryons measured at forward rapidity by the LHCb Collaboration. The ${\rm b\overline{b}}$ production cross section per unit of rapidity at midrapidity, estimated from these measurements, is ${\rm d}\sigma_{\rm b\overline{b}}/{\rm d}y|_{|y|<0.5} = 83.1 \pm 3.5 (\mathrm{stat.}) \pm 5.4(\mathrm{syst.}) ^{+12.3}_{-3.2} (\mathrm{extrap.})\,\mu$b. The baryon-to-meson ratios are computed to investigate the hadronization mechanism of beauty quarks. The non-prompt $\Lambda^+_{\rm c}/{\rm D^0}$ production ratio has a similar trend to the one measured for the promptly produced charmed particles and to the p$/\pi^+$ and $\Lambda/{\rm K^0_S}$ ratios, suggesting a similar baryon-formation mechanism among light, strange, charm, and beauty hadrons. The $p_{\rm T}$-integrated non-prompt $\Lambda_{\rm c}/{\rm D^0}$ ratio is found to be significantly higher than the one measured in e$^+$e$^-$ collisions.

Measurements of the production cross sections of prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D^{\ast +}}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, and ${\rm \Xi_{c}^{+}}$ charm hadrons at midrapidity in proton$-$proton collisions at $\sqrt{s}=13$ TeV with the ALICE detector are presented. The D-meson cross sections as a function of transverse momentum ($p_{\rm T}$) are provided with improved precision and granularity. The ratios of $p_{\rm T}$-differential meson production cross sections based on this publication and on measurements at different rapidity and collision energy provide a constraint on gluon parton distribution functions at low values of Bjorken-$x$ ($10^{-5}-10^{-4}$). The measurements of ${\rm \Lambda_{c}^{+}}$ (${\rm \Xi_{c}^{+}}$) baryon production extend the measured $p_{\rm T}$ intervals down to $p_{\rm T}=0(3)$~GeV$/c$. These measurements are used to determine the charm-quark fragmentation fractions and the ${\rm c\overline{c}}$ production cross section at midrapidity ($|y|<0.5$) based on the sum of the cross sections of the weakly-decaying ground-state charm hadrons ${\rm D^0}$, ${\rm D^+}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, ${\rm \Xi_{c}^{0}}$ and, for the first time, ${\rm \Xi_{c}^{+}}$, and of the strongly-decaying J/$psi$ mesons. The first measurements of ${\rm \Xi_{c}^{+}}$ and ${\rm \Sigma_{c}^{0,++}}$ fragmentation fractions at midrapidity are also reported. A significantly larger fraction of charm quarks hadronising to baryons is found compared to e$^+$e$^-$ and ep collisions. The ${\rm c\overline{c}}$ production cross section at midrapidity is found to be at the upper bound of state-of-the-art perturbative QCD calculations.

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.

Charged-particle jet yields have been measured in semicentral Pb$-$Pb collisions at center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}} = 5$ TeV with the ALICE detector at the LHC. These yields are reported as a function of the jet transverse momentum, and further classified by their angle with respect to the event plane and the event shape, characterized by ellipticity, in an effort to study the path-length dependence of jet quenching. Jets were reconstructed at midrapidity from charged-particle tracks using the anti-$k_{\rm T}$ algorithm with resolution parameters $R =$ 0.2 and 0.4, with event-plane angle and event-shape values determined using information from forward scintillating detectors. The results presented in this letter show that, in semicentral Pb$-$Pb collisions, there is no significant difference between jet yields in predominantly isotropic and elliptical events. However, out-of-plane jets are observed to be more suppressed than in-plane jets. Further, this relative suppression is greater for low transverse momentum ($<$ 50 GeV/$c$) $R =$ 0.2 jets produced in elliptical events, with out-of-plane to in-plane jet-yield ratios varying up to 5.2$\sigma$ between different event-shape classes. These results agree with previous studies indicating that jets experience azimuthally anisotropic suppression when traversing the QGP medium, and can provide additional constraints on the path-length dependence of jet energy loss.

Measurements of inclusive charged-particle jet production in pp and p-Pb collisions at center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}} = 5.02$ TeV and the corresponding nuclear modification factor $R_{\rm pPb}^{\rm ch\,jet}$ are presented, using data collected with the ALICE detector at the LHC. Jets are reconstructed in the central rapidity region $|\eta_{\rm jet}| < 0.5$ from charged particles using the anti-$k_{\rm T}$ algorithm with resolution parameters $R = 0.2$, 0.3, and 0.4. The $p_{\rm T}$-differential inclusive production cross section of charged-particle jets, as well as the corresponding cross-section ratios, are reported for pp and p-Pb collisions in the transverse momentum range $10 < p^{\rm ch}_{\rm T,jet} < 140$ GeV/$c$ and $10 < p^{\rm ch}_{\rm T,jet} < 160$ GeV/$c$, respectively, together with the nuclear modification factor $R_{\rm pPb}^{\rm ch\,jet}$ in the range $10 < p^{\rm ch}_{\rm T,jet} < 140$ GeV/$c$. The analysis extends the $p_{\rm T}$ range of the previously-reported charged-particle jet measurements by the ALICE Collaboration. The nuclear modification factor is found to be consistent with one and independent of the jet resolution parameter with the improved precision of this study, indicating that the possible influence of cold nuclear matter effects on the production cross section of charged-particle jets in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV is smaller than the current precision. The obtained results are in agreement with other minimum bias jet measurements available for RHIC and LHC energies, and are well reproduced by the NLO perturbative QCD POWHEG calculations with parton shower provided by PYTHIA8 as well as by JETSCAPE simulations.

The pseudorapidity dependence of elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of charged particles measured in Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}=5.02$ TeV and in Xe-Xe collisions at $\sqrt{s_{\rm NN}}=5.44$ TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range $-3.5 < \eta < 5$ for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient $v_2$ calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to $3+1$ dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.

The cross section for coherent photonuclear production of J/$\psi$ is presented as a function of the electromagnetic dissociation (EMD) of Pb. The measurement is performed with the ALICE detector in ultra-peripheral Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV. Cross sections are presented in five different J/$\psi$ rapidity ranges within $|y|<4$, with the J/$\psi$ reconstructed via its dilepton decay channels. In some events the J/$\psi$ is not accompanied by EMD, while other events do produce neutrons from EMD at beam rapidities either in one or the other beam direction, or in both. The cross sections in a given rapidity range and for different configurations of neutrons from EMD allow for the extraction of the energy dependence of this process in the range $17 < W_{\gamma\, \mathrm{Pb, n}} <920$ GeV, where $W_{\gamma\, \mathrm{Pb, n}}$ is the centre-of-mass energy per nucleon of the $\gamma\,\mathrm{Pb}$ system. This range corresponds to a Bjorken-$x$ interval spanning about three orders of magnitude: $ 1.1\times10^{-5}<x<3.3\times 10^{-2}$. In addition to the ultra-peripheral and photonuclear cross sections, the nuclear suppression factor is obtained. These measurements point to a strong depletion of the gluon distribution in Pb nuclei over a broad, previously unexplored, energy range. These results, together with previous ALICE measurements, provide unprecedented information to probe quantum chromodynamics at high energies.

The ALICE Collaboration reports three measurements in ultra-peripheral proton$-$lead collisions at forward rapidity. The exclusive two-photon process \ggmm and the exclusive photoproduction of J/$\psi$ are studied. J/$\psi$ photoproduction with proton dissociation is measured for the first time at a hadron collider. The cross section for the two-photon process of dimuons in the invariant mass range from 1 to 2.5 GeV/$c^2$ agrees with leading order quantum electrodynamics calculations. The exclusive and dissociative cross sections for J/$\psi$ photoproductions are measured for photon$-$proton centre-of-mass energies from 27 to 57 GeV. They are in good agreement with HERA results.

Measurements of the production of electrons from heavy-flavour hadron decays in pp collisions at $\sqrt{s} = 13$ TeV at midrapidity with the ALICE detector are presented down to a transverse momentum ($p_{\rm T}$) of 0.2 GeV$/c$ and up to $p_{\rm T} = 35$ GeV$/c$, which is the largest momentum range probed for inclusive electron measurements in ALICE. In p$-$Pb collisions, the production cross section and the nuclear modification factor of electrons from heavy-flavour hadron decays are measured in the $p_{\rm T}$ range $0.5 < p_{\rm T} < 26$ GeV$/c$ at $\sqrt{s_{\rm NN}} = 8.16$ TeV. The nuclear modification factor is found to be consistent with unity within the statistical and systematic uncertainties. In both collision systems, first measurements of the yields of electrons from heavy-flavour hadron decays in different multiplicity intervals normalised to the multiplicity-integrated yield (self-normalised yield) at midrapidity are reported as a function of the self-normalised charged-particle multiplicity estimated at midrapidity. The self-normalised yields in pp and p$-$Pb collisions grow faster than linear with the self-normalised multiplicity. A strong $p_{\rm T}$ dependence is observed in pp collisions, where the yield of high-$p_{\rm T}$ electrons increases faster as a function of multiplicity than the one of low-$p_{\rm T}$ electrons. The measurement in p$-$Pb collisions shows no $p_{\rm T}$ dependence within uncertainties. The self-normalised yields in pp and p$-$Pb collisions are compared with measurements of other heavy-flavour, light-flavour, and strange particles, and with Monte Carlo simulations.

The production cross section of inclusive J/$\psi$ pairs in pp collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV is measured with ALICE. The measurement is performed for J/$\psi$ in the rapidity interval $2.5 < y < 4.0$ and for transverse momentum $p_{\rm T} > 0$. The production cross section of inclusive J/$\psi$ pairs is reported to be $10.3 \pm 2.3 {\rm (stat.)} \pm 1.3 {\rm (syst.)}$ nb in this kinematic interval. The contribution from non-prompt J/$\psi$ (i.e. originated from beauty-hadron decays) to the inclusive sample is evaluated. The results are discussed and compared with data.

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 correlations between different moments of two flow amplitudes, extracted with the recently developed asymmetric cumulants, are measured in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV recorded by the ALICE detector at the CERN Large Hadron Collider. The magnitudes of the measured observables show a dependence on the different moments as well as on the collision centrality, indicating the presence of non-linear response in all even moments up to the eighth. Furthermore, the higher-order asymmetric cumulants show different signatures than the symmetric and lower-order asymmetric cumulants. Comparisons with state-of-the-art event generators using two different parametrizations obtained from Bayesian optimization show differences between data and simulations in many of the studied observables, indicating a need for further tuning of the models behind those event generators. These results provide new and independent constraints on the initial conditions and transport properties of the system created in heavy-ion collisions.

The azimuthal ($\Delta\varphi$) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p$-$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV. Results are reported for electrons with transverse momentum $4<p_{\rm T}<16$ GeV/$c$ and pseudorapidity $|\eta|<0.6$. The associated charged particles are selected with transverse momentum $1<p_{\rm T}<7$ GeV/$c$, and relative pseudorapidity separation with the leading electron $|\Delta\eta| < 1$. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p$-$Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The $\Delta\varphi$ distribution and the peak observables in pp and p$-$Pb collisions are compared with calculations from various Monte Carlo event generators.

A study of multiplicity and pseudorapidity distributions of inclusive photons measured in pp and p$-$Pb collisions at a center-of-mass energy per nucleon$-$nucleon collision of $\sqrt{s_{\rm NN}} = 5.02$ TeV using the ALICE detector in the forward pseudorapidity region $2.3 < \eta_{\rm lab} < 3.9$ is presented. Measurements in p$-$Pb collisions are reported for two beam configurations in which the directions of the proton and lead ion beam were reversed. The pseudorapidity distributions in p$-$Pb collisions are obtained for seven centrality classes which are defined based on different event activity estimators, i.e., the charged-particle multiplicity measured at midrapidity as well as the energy deposited in a calorimeter at beam rapidity. The inclusive photon multiplicity distributions for both pp and p$-$Pb collisions are described by double negative binomial distributions. The pseudorapidity distributions of inclusive photons are compared to those of charged particles at midrapidity in \pp collisions and for different centrality classes in p$-$Pb collisions. The results are compared to predictions from various Monte Carlo event generators. None of the generators considered in this paper reproduces the inclusive photon multiplicity distributions in the reported multiplicity range. The pseudorapidity distributions are, however, better described by the same generators.

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.

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.

A newly developed observable for correlations between symmetry planes, which characterize the direction of the anisotropic emission of produced particles, is measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with ALICE. This so-called Gaussian Estimator allows for the first time the study of these quantities without the influence of correlations between different flow amplitudes. The centrality dependence of various correlations between two, three and four symmetry planes is presented. The ordering of magnitude between these symmetry plane correlations is discussed and the results of the Gaussian Estimator are compared with measurements of previously used estimators. The results utilizing the new estimator lead to significantly smaller correlations than reported by studies using the Scalar Product method. Furthermore, the obtained symmetry plane correlations are compared to state-of-the-art hydrodynamic model calculations for the evolution of heavy-ion collisions. While the model predictions provide a qualitative description of the data, quantitative agreement is not always observed, particularly for correlators with significant non-linear response of the medium to initial state anisotropies of the collision system. As these results provide unique and independent information, their usage in future Bayesian analysis can further constrain our knowledge on the properties of the QCD matter produced in ultrarelativistic heavy-ion collisions.

Recent measurements of charm-baryon production in hadronic collisions have questioned the universality of charm-quark fragmentation across different collision systems. In this work the fragmentation of charm quarks into charm baryons is probed, by presenting the first measurement of the longitudinal jet momentum fraction carried by $\Lambda_{\rm c}^{+}$ baryons, $z^{\mathrm{ch}}_\mathrm{||}$, in hadronic collisions. The results are obtained in proton-proton (pp) collisions at $\sqrt{s}=13$ TeV at the LHC, with $\Lambda_{\rm c}^{+}$ baryons and charged (track-based) jets reconstructed in the transverse momentum intervals of $3 \leq p_{\rm T}^{\Lambda_{\rm c}^{+}} < 15$ GeV/$c$ and $7 \leq p_{\rm T}^{\rm jet\;ch} < 15$ GeV/$c$, respectively. The $z^{\mathrm{ch}}_\mathrm{||}$ distribution is compared to a measurement of ${\rm D}^0$-tagged charged jets in pp collisions as well as to PYTHIA 8 simulations. The data hints that the fragmentation of charm quarks into charm baryons is softer with respect to charm mesons, in the measured kinematic interval, as predicted by hadronisation models which include colour correlations beyond leading-colour in the string formation.

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 pseudorapidity density of charged particles with minimum transverse momentum ($p_{\rm T}$) thresholds of 0.15, 0.5, 1, and 2 GeV$/c$ is measured in pp collisions at the centre of mass energies of $\sqrt{s} =$ 5.02 and 13 TeV with the ALICE detector. The study is carried out for inelastic collisions with at least one primary charged particle having a pseudorapidity ($\eta$) within $\pm0.8$ and $p_{\rm T}$ larger than the corresponding threshold. In addition, measurements without $p_{\rm T}$-thresholds are performed for inelastic and non-single-diffractive events as well as for inelastic events with at least one charged particle having $|\eta|<1$ in pp collisions at $\sqrt{s} =$ 5.02 TeV for the first time at the LHC. These measurements are compared to the PYTHIA 6, PYTHIA 8, and EPOS-LHC models. In general, the models describe the $\eta$ dependence of particle production well. However, discrepancies are observed for the highest transverse momentum threshold ($p_{\rm T}>2 {\rm\ GeV}/c$), highlighting the importance of such measurements for tuning event generators. The new measurements agree within uncertainties with results from the ATLAS and CMS experiments obtained at $\sqrt{s} = 13$ TeV.

Multiplicity ($N_{\rm ch}$) distributions and transverse momentum ($p_{\rm T}$) spectra of inclusive primary charged particles in the kinematic range of $|\eta| < 0.8$ and 0.15 GeV/$c$$< p_{T} <$ 10 GeV/$c$ are reported for pp, p-Pb, Xe-Xe and Pb-Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{\rm NN}} = 2.76$ TeV up to $13$ TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the $p_{\rm T}$ spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

The production of inclusive, prompt and non-prompt J/$\psi$ was studied for the first time at midrapidity ($ -1.37 < y_{\rm cms} < 0.43$) in p$-$Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV with the ALICE detector at the LHC. The inclusive J/$\psi$ mesons were reconstructed in the dielectron decay channel in the transverse momentum ($p_{\rm T}$) interval $0 < p_{\rm T} < 14$ GeV/$c$ and the prompt and non-prompt contributions were separated on a statistical basis for $p_{\rm T} > 2$ GeV/$c$. The study of the J/$\psi$ mesons in the dielectron channel used for the first time in ALICE online single-electron triggers from the Transition Radiation Detector, providing a data sample corresponding to an integrated luminosity of $689 \pm 13 \mu{\rm b}^{-1}$. The proton$-$proton reference cross section for inclusive J/$\psi$ was obtained based on interpolations of measured data at different centre-of-mass energies and a universal function describing the $p_{\rm T}$-differential J/$\psi$ production cross sections. The $p_{\rm T}$-differential nuclear modification factors $R_{\rm pPb}$ of inclusive, prompt, and non-prompt J/$\psi$ are consistent with unity and described by theoretical models implementing only nuclear shadowing.

The production of electrons from beauty-hadron decays was measured at midrapidity in proton-proton (pp) and central Pb-Pb collisions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV, using the ALICE detector at the LHC. The cross section measured in pp collisions in the transverse momentum interval $2 < p_{\rm T} < 8$ GeV/$c$ was compared with models based on perturbative quantum chromodynamics calculations. The yield in the 10% most central Pb-Pb collisions, measured in the interval $2 < p_{\rm T} < 26$ GeV/$c$, was used to compute the nuclear modification factor $R_{\rm AA}$, extrapolating the pp reference cross section to $p_{\rm T}$ larger than 8 GeV/$c$. The measured $R_{\rm AA}$ shows significant suppression of the yield of electrons from beauty-hadron decays at high $p_{\rm T}$ and does not show a significant dependence on $p_{\rm T}$ above 8 GeV/$c$ within uncertainties. The results are described by several theoretical models based on different implementations of the interaction of heavy quarks with a quark-gluon plasma, which predict a smaller energy loss for beauty quarks compared to light and charm quarks.

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.

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.

This article reports measurements of the angle between differently defined jet axes in pp collisions at $\sqrt{s} = 5.02$ TeV carried out by the ALICE Collaboration. Charged particles at midrapidity are clustered into jets with resolution parameters $R=0.2$ and 0.4. The jet axis, before and after Soft Drop grooming, is compared to the jet axis from the Winner-Takes-All (WTA) recombination scheme. The angle between these axes, $\Delta R_{\mathrm{axis}}$, probes a wide phase space of the jet formation and evolution, ranging from the initial high-momentum-transfer scattering to the hadronization process. The $\Delta R_{\mathrm{axis}}$ observable is presented for $20 < {p_{\mathrm{T}}^{\mathrm{ch\; jet}}}< 100$ GeV/$c$, and compared to predictions from the PYTHIA 8 and Herwig 7 event generators. The distributions can also be calculated analytically with a leading hadronization correction related to the non-perturbative component of the Collins$-$Soper$-$Sterman (CSS) evolution kernel. Comparisons to analytical predictions at next-to-leading-logarithmic accuracy with leading hadronization correction implemented from experimental extractions of the CSS kernel in Drell$-$Yan measurements are presented. The analytical predictions describe the measured data within 20% in the perturbative regime, with surprising agreement in the non-perturbative regime as well. These results are compatible with the universality of the CSS kernel in the context of jet substructure.

The production of strange hadrons (K$^{0}_{\rm S}$, $\Lambda$, $\Xi^{\pm}$, and $\Omega^{\pm}$), baryon-to-meson ratios ($\Lambda/{\rm K}^0_{\rm S}$, $\Xi/{\rm K}^0_{\rm S }$, and $\Omega/{\rm K}^0_{\rm S}$), and baryon-to-baryon ratios ($\Xi/\Lambda$, $\Omega/\Lambda$, and $\Omega/\Xi$) associated with jets and the underlying event were measured as a function of transverse momentum ($p_{\rm T}$) in pp collisions at $\sqrt{s} = 13$ TeV and p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The inclusive production of the same particle species and the corresponding ratios are also reported. The production of multi-strange hadrons, $\Xi^{\pm}$ and $\Omega^{\pm}$, and their associated particle ratios in jets and in the underlying event are measured for the first time. In both pp and p-Pb collisions, the baryon-to-meson and baryon-to-baryon yield ratios measured in jets differ from the inclusive particle production for low and intermediate hadron $p_{\rm T}$ (0.6$-$6 GeV/$c$). Ratios measured in the underlying event are in turn similar to those measured for inclusive particle production. In pp collisions, the particle production in jets is compared with PYTHIA 8 predictions with three colour-reconnection implementation modes. None of them fully reproduces the data in the measured hadron $p_{\rm T}$ region. The maximum deviation is observed for $\Xi^{\pm}$ and $\Omega^{\pm}$, which reaches a factor of about six. In p-Pb collisions, there is no significant event-multiplicity dependence for particle production in jets, in contrast to what is observed in the underlying event. The presented measurements provide novel constraints on hadronisation and its Monte Carlo description. In particular, they demonstrate that the fragmentation of jets alone is insufficient to describe the strange and multi-strange particle production in hadronic collisions at LHC energies.

Two-particle correlations with ${\rm K}^{0}_{\rm{S}}$, $\Lambda$/$\bar{\Lambda}$, and charged hadrons as trigger particles in the transverse momentum range $8 < p_\mathrm{T,trig}<16$ GeV/$c$, and associated charged particles within $1 < p_\mathrm{T,assoc}<8$ GeV/$c$, are studied at mid-rapidity in pp and central Pb-Pb collisions at a centre-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV with the ALICE detector at the LHC. After subtracting the contributions of the flow background, the per-trigger yields are extracted on both the near and away sides, and the ratio in Pb-Pb collisions with respect to pp collisions ($I_{\mathrm {AA}}$) is computed. The per-trigger yield in Pb-Pb collisions on the away side is strongly suppressed to the level of $I_{\mathrm {AA}} \approx 0.6$ for $p_\mathrm{T,assoc}>3$ GeV/$c$ as expected from strong in-medium energy loss, while an enhancement develops at low $p_\mathrm{T,assoc}$ on both the near and away sides, reaching $I_{\mathrm {AA}} \approx 1.8$ and $2.7$ respectively. These findings are in good agreement with previous ALICE measurements from two-particle correlations triggered by neutral pions ($\pi^{0}$-h) and charged hadrons (h-h) in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 2.76$ TeV. Moreover, the correlations with ${\rm K}^{0}_{\rm{S}}$ mesons and $\Lambda$/$\bar{\Lambda}$ baryons as trigger particles are compared to those of inclusive charged hadrons. The results are compared with the predictions of Monte Carlo models.

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 of the $\psi(2S)$ charmonium state was measured with ALICE in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the dimuon decay channel. A significant signal was observed for the first time at LHC energies down to zero transverse momentum, at forward rapidity ($2.5<y<4$). The measurement of the ratio of the inclusive production cross sections of the $\psi(2S)$ and J/$\psi$ resonances is reported as a function of the centrality of the collisions and of transverse momentum, in the region $p_{\rm T}<12$ GeV/$c$. The results are compared with the corresponding measurements in pp collisions, by forming the double ratio $[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{Pb-Pb}}/[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{pp}}$. It is found that in Pb-Pb collisions the $\psi(2S)$ is suppressed by a factor of $\sim 2$ with respect to the J/$\psi$. The $\psi(2S)$ nuclear modification factor $R_{\rm AA}$ was also obtained as a function of both centrality and $p_{\rm T}$. The results show that the $\psi(2S)$ resonance yield is strongly suppressed in Pb-Pb collisions, by a factor up to $\sim 3$ with respect to pp. Comparisons of cross section ratios with previous SPS findings by the NA50 experiment and of $R_{\rm AA}$ with higher-$p_{\rm T}$ results at LHC energy are also reported. These results and the corresponding comparisons with calculations of transport and statistical models address questions on the presence and properties of charmonium states in the quark-gluon plasma formed in nuclear collisions at the LHC.

Understanding the role of parton mass and Casimir colour factors in the quantum chromodynamics parton shower represents an important step in characterising the emission properties of heavy quarks. Recent experimental advances in jet substructure techniques have provided the opportunity to isolate and characterise gluon emissions from heavy quarks. In this work, the first direct experimental constraint on the charm-quark splitting function is presented, obtained via the measurement of the groomed shared momentum fraction of the first splitting in charm jets, tagged by a reconstructed ${\rm D}^0$ meson. The measurement is made in proton-proton collisions at $\sqrt{s}$ = 13 TeV, in the low jet transverse-momentum interval of $15 \leq p_{\rm T}^{\rm jet~ch} < 30$ GeV/$c$ where the emission properties are sensitive to parton mass effects. In addition, the opening angle of the first perturbative emission of the charm quark, as well as the number of perturbative emissions it undergoes, is reported. Comparisons to measurements of an inclusive-jet sample show a steeper splitting function for charm quarks compared with gluons and light quarks. Charm quarks also undergo fewer perturbative emissions in the parton shower, with a reduced probability of large-angle emissions.

The measurement of the production of ${\rm f}_{0}(980)$ in inelastic pp collisions at $\sqrt{s} = 5.02$ TeV is presented. This is the first reported measurement of inclusive ${\rm f}_{0}(980)$ yield at LHC energies. The production is measured at midrapidity, $|y| < 0.5$, in a wide transverse momentum range, $0 < p_{\rm T} < 16$ GeV/$c$, by reconstructing the resonance in the ${\rm f}_{0}(980) \rightarrow \pi^{+}\pi^{-}$ hadronic decay channel using the ALICE detector. The $p_{\rm T}$-differential yields are compared to those of pions, protons and $\phi$ mesons as well as to predictions from the HERWIG 7.2 QCD-inspired Monte Carlo event generator and calculations from a coalescence model that uses the AMPT model as an input. The ratio of the $p_{\rm T}$-integrated yield of ${\rm f}_{0}(980)$ relative to pions is compared to measurements in ${\rm e}^{+}{\rm e}^{-}$ and pp collisions at lower energies and predictions from statistical hadronisation models and HERWIG 7.2. A mild collision energy dependence of the ${\rm f}_{0}(980)$ to pion production is observed in pp collisions from SPS to LHC energies. All considered models underpredict the $p_{\rm T}$-integrated $2{\rm f}_{0}(980)/(\pi^{+}+\pi^{-})$ ratio. The prediction from the canonical statistical hadronisation model assuming a zero total strangeness content of ${\rm f}_{0}(980)$ is consistent with the data within 1.9$\sigma$ and is the closest to the data. The results provide an essential reference for future measurements of the particle yield and nuclear modification in p$-$Pb and Pb$-$Pb collisions, which have been proposed to be instrumental to probe the elusive nature and quark composition of the ${\rm f}_{0}(980)$ scalar meson.

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.

Fluctuation measurements are important sources of information on the mechanism of particle production at LHC energies. This article reports the first experimental results on third-order cumulants of the net-proton distributions in Pb$-$Pb collisions at a center-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$ TeV recorded by the ALICE detector. The results on the second-order cumulants of net-proton distributions at $\sqrt{s_{\rm NN}} = 2.76$ and $5.02$ TeV are also discussed in view of effects due to the global and local baryon number conservation. The results demonstrate the presence of long-range rapidity correlations between protons and antiprotons. Such correlations originate from the early phase of the collision. The experimental results are compared with HIJING and EPOS model calculations, and the dependence of the fluctuation measurements on the phase-space coverage is examined in the context of lattice quantum chromodynamics (LQCD) and hadron resonance gas (HRG) model estimations. The measured third-order cumulants are consistent with zero within experimental uncertainties of about 4% and are described well by LQCD and HRG predictions.

The inclusive production of the charm-strange baryon $\rm \Omega_c^0$ is measured for the first time via its hadronic decay into $\rm \Omega^-\pi^+$ at midrapidity ($|y|<0.5$) in proton-proton (pp) collisions at the centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval $2<p_{\rm T}<12~{\rm GeV}/c$. The $p_{\rm T}$ dependence of the $\rm \Omega_c^0$-baryon production relative to the prompt $\rm D^0$-meson and to the prompt $\rm \Xi_c^0$-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured $p_{\rm T}$ interval, the ratio of the $p_{\rm T}$-integrated cross sections of $\rm \Omega_c^0$ and prompt $\Lambda_{\rm c}^{+}$ baryons multiplied by the $\rm \Omega^-\pi^+$ branching ratio is found to be larger by a factor of about 20 with a significance of about 4$\sigma$ when compared to $\rm e^+e^-$ collisions.

Hadronic resonances are used to probe the hadron gas produced in the late stage of heavy-ion collisions since they decay on the same timescale, of the order of 1 to 10 fm/$c$, as the decoupling time of the system. In the hadron gas, (pseudo)elastic scatterings among the products of resonances that decayed before the kinetic freeze-out and regeneration processes counteract each other, the net effect depending on the resonance lifetime, the duration of the hadronic phase, and the hadronic cross sections at play. In this context, the $\Sigma(1385)^{\pm}$ particle is of particular interest as models predict that regeneration dominates over rescattering despite its relatively short lifetime of about 5.5 fm/$c$. The first measurement of the $\Sigma(1385)^{\pm}$ resonance production at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}= 5.02$ TeV with the ALICE detector is presented in this Letter. The resonances are reconstructed via their hadronic decay channel, $\Lambda\pi$, as a function of the transverse momentum ($p_{\rm T}$) and the collision centrality. The results are discussed in comparison with the measured yield of pions and with expectations from the statistical hadronization model as well as commonly employed event generators, including PYTHIA8/Angantyr and EPOS3 coupled to the UrQMD hadronic cascade afterburner. None of the models can describe the data. For $\Sigma(1385)^{\pm}$, a similar behaviour as ${\rm K}^{*} (892)^{0}$ is observed in data unlike the predictions of EPOS3 with afterburner.

The production of the W$^\pm$ bosons measured in p$-$Pb collisions at a centre-of-mass energy per nucleon$-$nucleon collision $\sqrt{s_{NN}} = 8.16$ TeV and Pb$-$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV with ALICE at the LHC is presented. The W$^\pm$ bosons are measured via their muonic decay channel, with the muon reconstructed in the pseudorapidity region $-4 < \eta^\mu_{\rm lab} < -2.5$ with transverse momentum $p_{\rm T}^\mu > 10$ GeV/$c$. While in Pb$-$Pb collisions the measurements are performed in the forward ($2.5 < y^\mu_{\rm cms} < 4$) rapidity region, in p$-$Pb collisions, where the centre-of-mass frame is boosted with respect to the laboratory frame, the measurements are performed in the backward ($-4.46 < y^\mu_{\rm cms} < -2.96$) and forward ($2.03 < y^\mu_{\rm cms} < 3.53$) rapidity regions. The W$^{-}$ and W$^{+}$ production cross sections, lepton-charge asymmetry, and nuclear modification factors are evaluated as a function of the muon rapidity. In order to study the production as a function of the p$-$Pb collision centrality, the production cross sections of the W$^{-}$ and W$^{+}$ bosons are combined and normalised to the average number of binary nucleon$-$nucleon collision $\langle N_\mathrm{coll} \rangle$. In Pb$-$Pb collisions, the same measurements are presented as a function of the collision centrality. Study of the binary scaling of the W$^\pm$-boson cross sections in p$-$Pb and Pb$-$Pb collisions is also reported. The results are compared with perturbative QCD (pQCD) calculations, with and without nuclear modifications of the Parton Distribution Functions (PDFs), as well as with available data at the LHC. Significant deviations from the theory expectations are found in the two collision systems, indicating that the measurements can provide additional constraints for the determination of nuclear PDF (nPDFs) and in particular of the light-quark distributions.

An excess of J/$\psi$ yield at very low transverse momentum ($p_{\rm T} < 0.3$ GeV/$c$), originating from coherent photoproduction, is observed in peripheral and semicentral hadronic Pb$-$Pb collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}} = 5.02$ TeV. The measurement is performed with the ALICE detector via the dimuon decay channel at forward rapidity ($2.5<y<4$). The nuclear modification factor at very low $p_{\rm T}$ and the coherent photoproduction cross section are measured as a function of centrality down to the 10% most central collisions. These results extend the previous study at $\sqrt{s_{\rm NN}} = 2.76$ TeV, confirming the clear excess over hadronic production in the $p_{\rm T}$ range 0$-$0.3 GeV/$c$ and the centrality range 70$-$90%, and establishing an excess with a significance greater than 5$\sigma$ also in the 50$-$70% and 30$-$50% centrality ranges. The results are compared with earlier measurements at $\sqrt{s_{\rm NN}} = 2.76$ TeV and with different theoretical predictions aiming at describing how coherent photoproduction occurs in hadronic interactions with nuclear overlap.

The first experimental information on the strong interaction between $\Lambda$ and $\Xi^-$ strange baryons is presented in this Letter. The correlation function of $\Lambda-\Xi^-$ and $\overline{\Lambda}-\overline{\Xi}^{+}$ pairs produced in high-multiplicity proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the LHC is measured as a function of the relative momentum of the pair. The femtoscopy method is used to calculate the correlation function, which is then compared with theoretical expectations obtained using a meson exchange model, chiral effective field theory, and Lattice QCD calculations close to the physical point. Data support predictions of small scattering parameters while discarding versions with large ones, thus suggesting a weak $\Lambda-\Xi^{-}$ interaction. The limited statistical significance of the data does not yet allow one to constrain the effects of coupled channels like $\Sigma-\Xi$ and N$-\Omega$.

The measurement of the production of charm jets, identified by the presence of a ${\rm D^0}$ meson in the jet constituents, is presented in proton-proton collisions at centre-of-mass energies of $\sqrt{s}$ = 5.02 and 13 TeV with the ALICE detector at the CERN LHC. The ${\rm D^0}$ mesons were reconstructed from their hadronic decay ${\rm D^0} \rightarrow {\rm K^-}\pi^+$ and the respective charge conjugate. Jets were reconstructed from ${\rm D^0}$-meson candidates and charged particles using the anti-$k_{\rm T}$ algorithm, in the jet transverse momentum range $5<p_{\rm T;chjet}<50$ GeV/$c$, pseudorapidity $|\eta_{\rm jet}| <0.9-R$, and with the jet resolution parameters $R$ = 0.2, 0.4, 0.6. The distribution of the jet momentum fraction carried by a ${\rm D^0}$ meson along the jet axis ($z^{\rm ch}_{||}$) was measured in the range $0.4 < z^{\rm ch}_{||} < 1.0$ in four ranges of the jet transverse momentum. Comparisons of results for different collision energies and jet resolution parameters are also presented. The measurements are compared to predictions from Monte Carlo event generators based on leading-order and next-to-leading-order perturbative quantum chromodynamics calculations. A generally good description of the main features of the data is obtained in spite of a few discrepancies at low $p_{\rm T;chjet}$. Measurements were also done for $R = 0.3$ at $\sqrt{s}$ = 5.02 TeV and are shown along with their comparisons to theoretical predictions in an appendix to this paper.

This article presents measurements of the groomed jet radius and momentum splitting fraction in pp collisions at $\sqrt{s}=5.02$ TeV with the ALICE detector at the Large Hadron Collider. Inclusive charged-particle jets are reconstructed at midrapidity using the anti-$k_{\rm{T}}$ algorithm for transverse momentum $60< p_{\mathrm{T}}^{\rm{ch\; jet}}<80$ GeV/$c$. We report results using two different grooming algorithms: soft drop and, for the first time, dynamical grooming. For each grooming algorithm, a variety of grooming settings are used in order to explore the impact of collinear radiation on these jet substructure observables. These results are compared to perturbative calculations that include resummation of large logarithms at all orders in the strong coupling constant. We find good agreement of the theoretical predictions with the data for all grooming settings considered.

We present the first systematic comparison of the charged-particle pseudorapidity densities for three widely different collision systems, pp, p-Pb, and Pb-Pb, at the top energy of the Large Hadron Collider ($\sqrt{s_{\rm NN}} = 5.02$ TeV) measured over a wide pseudorapidity range (${-3.5 <\eta <5}$), the widest possible among the four experiments at that facility. The systematic uncertainties are minimised since the measurements are recorded by the same experimental apparatus (ALICE). The distributions for p-Pb and Pb-Pb collisions are determined as a function of the centrality of the collisions, while results from pp collisions are reported for inelastic events with at least one charged particle at midrapidity. The charged-particle pseudorapidity densities are, under simple and robust assumptions, transformed to charged-particle rapidity densities. This allows for the calculation and the presentation of the evolution of the width of the rapidity distributions and of a lower bound on the Bjorken energy density, as a function of the number of participants in all three collision systems. We find a decreasing width of the particle production, and roughly a smooth ten fold increase in the energy density, as the system size grows, which is consistent with a gradually higher dense phase of matter.

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.

Production of inclusive charmonia in pp collisions at center-of-mass energy of $\sqrt{s}$ = 13 TeV and p-Pb collisions at center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 8.16 TeV is studied as a function of charged-particle pseudorapidity density with ALICE. Ground and excited charmonium states (J/$\psi$, $\psi$(2S)) are measured from their dimuon decays in the interval of rapidity in the center-of-mass frame $2.5 < y_{\rm cms} < 4.0$ for pp collisions, and $2.03 < y_{\rm cms} < 3.53$ and $-4.46 < y_{\rm cms} < -2.96$ for p-Pb collisions. The charged-particle pseudorapidity density is measured around midrapidity ($|\eta|<1.0$). In pp collisions, the measured charged-particle multiplicity extends to about six times the average value, while in p-Pb collisions at forward (backward) rapidity a multiplicity corresponding to about three (four) times the average is reached. The $\psi$(2S) yield increases with the charged-particle pseudorapidity density. The ratio of $\psi$(2S) over J/$\psi$ yield does not show a significant multiplicity dependence in either colliding system, suggesting a similar behavior of J/$\psi$ and $\psi$(2S) yields with respect to charged-particle pseudorapidity density. Results for the $\psi$(2S) yield and its ratio with respect to J/$\psi$ agree with available model calculations.

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons $v_2${$\Psi_{\rm SP}$} in collisions of Pb ions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$=2.76 TeV and Xe ions at $\sqrt{s_{\rm NN}}$=5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum. The ratio between $v_2${$\Psi_{\rm SP}$} and the elliptic flow coefficient relative to the participant plane $v_2$4, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the TRENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models of initial state fluctuations. A significant transverse momentum dependence of the ratio $v_2${$\Psi_{\rm SP}$}/$v_2${4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of $v_2${$\Psi_{\rm SP}$} and $v_2${4} to the corresponding initial state eccentricities for Xe-Xe and Pb-Pb collisions at similar initial entropy density show a difference of $(7.0 \pm 0.9)$% with an additional variation of +1.8% when including RHIC data in the TRENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark-gluon plasma.