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.

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.

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

The study of the production of nuclei and antinuclei in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high-energy hadronic collisions. In this paper, the production of protons, deuterons and $^3$He and their charge conjugates at midrapidity is studied as a function of the charged-particle multiplicity in inelastic pp collisions at $\sqrt{s}=5.02$ TeV using the ALICE detector. Within the uncertainties, the yields of nuclei in pp collisions at $\sqrt{s}=5.02$ TeV are compatible with those in pp collisions at different energies and to those in p-Pb collisions when compared at similar multiplicities. The measurements are compared with the expectations of coalescence and Statistical Hadronisation Models. The results suggest a common formation mechanism behind the production of light nuclei in hadronic interactions and confirm that they do not depend on the collision energy but on the number of produced particles.

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 production of $\Lambda$ baryons and ${\rm K}^{0}_{\rm S}$ mesons (${\rm V}^{0}$ particles) was measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and pp collisions at $\sqrt{s} = 7$ TeV with ALICE at the LHC. The production of these strange particles is studied separately for particles associated with hard scatterings and the underlying event to shed light on the baryon-to-meson ratio enhancement observed at intermediate transverse momentum ($p_{\rm T}$) in high multiplicity pp and p-Pb collisions. Hard scatterings are selected on an event-by-event basis with jets reconstructed with the anti-$k_{\rm T}$ algorithm using charged particles. The production of strange particles associated with jets $p_{\rm T,\;jet}^{\rm ch}>10$ and $p_{\rm T,\;jet}^{\rm ch}>20$ GeV/$c$ in p-Pb collisions, and with jet $p_{\rm T,\;jet}^{\rm ch}>10$ GeV/$c$ in pp collisions is reported as a function of $p_{\rm T}$. Its dependence on angular distance from the jet axis, $R({\rm V}^{0},\;{\rm jet})$, for jets with $p_{\rm T,\;jet}^{\rm ch}>10$ GeV/$c$ in p-Pb collisions is reported as well. The $p_{\rm T}$-differential production spectra of strange particles associated with jets are found to be harder compared to that in the underlying event and both differ from the inclusive measurements. In events containing a jet, the density of the ${\rm V}^{0}$ particles in the underlying event is found to be larger than the density in the minimum bias events. The $\Lambda/{\rm K}^{0}_{\rm S}$ ratio associated with jets in p-Pb collisions is consistent with the ratio in pp collisions and follows the expectation of jets fragmenting in vacuum. On the other hand, this ratio within jets is consistently lower than the one obtained in the underlying event and it does not show the characteristic enhancement of baryons at intermediate $p_{\rm T}$ often referred to as "baryon anomaly" in the inclusive measurements.

The production of $\phi$ mesons has been studied in pp collisions at LHC energies with the ALICE detector via the dimuon decay channel in the rapidity region $2.5 < y < 4$. Measurements of the differential cross section ${\rm d}^2\sigma/{\rm d}y {\rm d}p_{\rm T}$ are presented as a function of the transverse momentum ($p_{\rm T}$) at the center-of-mass energies $\sqrt{s}=5.02$, 8 and 13 TeV and compared with the ALICE results at midrapidity. The differential cross sections at $\sqrt{s}=5.02$ and 13 TeV are also studied in several rapidity intervals as a function of $p_{\rm T}$, and as a function of rapidity in three $p_{\rm T}$ intervals. A hardening of the $p_{\rm T}$-differential cross section with the collision energy is observed, while, for a given energy, $p_{\rm T}$ spectra soften with increasing rapidity and, conversely, rapidity distributions get slightly narrower at increasing $p_{\rm T}$. The new results, complementing the published measurements at $\sqrt{s}=2.76$ and 7 TeV, allow one to establish the energy dependence of $\phi$ meson production and to compare the measured cross sections with phenomenological models. None of the considered models manages to describe the evolution of the cross section with $p_{\rm T}$ and rapidity at all the energies.

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.

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

The jet angularities are a class of jet substructure observables which characterize the angular and momentum distribution of particles within jets. These observables are sensitive to momentum scales ranging from perturbative hard scatterings to nonperturbative fragmentation into final-state hadrons. We report measurements of several groomed and ungroomed jet angularities in pp collisions at $\sqrt{s}=5.02$ TeV with the ALICE detector. Jets are reconstructed using charged particle tracks at midrapidity ($|\eta| < 0.9$). The anti-$k_{\rm T}$ algorithm is used with jet resolution parameters $R=0.2$ and $R=0.4$ for several transverse momentum $p_{\rm T}^{\text{ch jet}}$ intervals in the 20$-$100 GeV/$c$ range. Using the jet grooming algorithm Soft Drop, the sensitivity to softer, wide-angle processes, as well as the underlying event, can be reduced in a way which is well-controlled in theoretical calculations. We report the ungroomed jet angularities, $\lambda_{\alpha}$, and groomed jet angularities, $\lambda_{\alpha\text{,g}}$, to investigate the interplay between perturbative and nonperturbative effects at low jet momenta. Various angular exponent parameters $\alpha = 1$, 1.5, 2, and 3 are used to systematically vary the sensitivity of the observable to collinear and soft radiation. Results are compared to analytical predictions at next-to-leading-logarithmic accuracy, which provide a generally good description of the data in the perturbative regime but exhibit discrepancies in the nonperturbative regime. Moreover, these measurements serve as a baseline for future ones in heavy-ion collisions by providing new insight into the interplay between perturbative and nonperturbative effects in the angular and momentum substructure of jets. They supply crucial guidance on the selection of jet resolution parameter, jet transverse momentum, and angular scaling variable for jet quenching studies.

Recent $p_{\rm T}$-integrated cross section measurements of the ground-state charm mesons and baryons, D$^{\rm 0}$, D$^+$, D$_{\rm s}^{+}$, $\Lambda_{\rm c}^{+}$, and $\Xi_{\rm c}^0$, are used to evaluate the charm fragmentation fractions and production cross section per unit of rapidity at midrapidity ($|y|<0.5$), in pp collisions at $\sqrt{s} = 5.02$ TeV at the LHC. The latter is ${\rm d} \sigma^{\rm c \overline{c}}/{\rm d} y|_{|y| < 0.5}$ =1165 $\pm 44(\rm{stat})^{+134}_{-101}(\rm{syst})$$\mu b$. These measurements were obtained for the first time in hadronic collisions at the LHC including the charm baryon states, recently measured by ALICE at midrapidity. The charm fragmentation fractions differ significantly from the values measured in e$^+$e$^-$ and ep collisions, providing evidence of the dependence of the parton-to-hadron fragmentation fractions on the collision system, indicating that the assumption of their universality is not supported by the measured cross sections. An increase of a factor of about 3.3 for the fragmentation fraction for the $\Lambda_{\rm c}^{+}$ with a significance of $5\,\sigma$ between the values obtained in pp collisions and those obtained in e$^+$e$^-$ (ep) collisions is reported. The fragmentation fraction for the $\Xi_{\rm c}^0$ was obtained for the first time in any collision system. The measured fragmentation fractions were used to update the $\rm c \overline{c}$ cross sections per unit of rapidity at $|y|<0.5$ at $\sqrt{s} = 2.76$ and 7 TeV, which are about 40% higher than the previously published results. The data were compared with perturbative-QCD calculations and lie at the upper edge of the theoretical bands.

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.

In quantum scattering processes between two particles, aspects characterizing the strong and Coulomb forces can be observed in kinematic distributions of the particle pairs. The sensitivity to the interaction potential reaches a maximum at low relative momentum and vanishing distance between the two particles. Ultrarelativistic heavy-ion collisions at the LHC provide an abundant source of many hadron species and can be employed as a measurement method of scattering parameters that is complementary to scattering experiments. This study confirms that momentum correlations of particles produced in Pb-Pb collisions at the LHC provide an accurate measurement of kaon-proton scattering parameters at low relative momentum, allowing precise access to the $ {K}^{-} p\rightarrow {K}^{-} p$ process. This work also validates the femtoscopic measurement in ultrarelativistic heavy-ion collisions as an alternative to scattering experiments and a complementary tool to the study of exotic atoms with comparable precision. In this work, the first femtoscopic measurement of momentum correlations of ${K}^{-} p\ ({K}^{+}\overline{p})$ and ${K}^{+}p ({K}^{-}\overline{p})$ pairs in Pb-Pb collisions at centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}} = 5.02$ TeV registered by the ALICE experiment is reported. The components of the ${K}^{-} p$ complex scattering length are extracted and found to be $\Re f_0=-0.91\pm~{0.03}$(stat)$^{+0.17}_{-0.03}$(syst) and $\Im f_0 = 0.92\pm~{0.05}$(stat)$^{+0.12}_{-0.33}$(syst). The results are compared with chiral effective field theory predictions as well as with existing data from dedicated scattering and exotic kaonic atom experiments.

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.

Correlations between moments of different flow coefficients are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV recorded with the ALICE detector. These new measurements are based on multiparticle mixed harmonic cumulants calculated using charged particles in the pseudorapidity region $|\eta|<0.8$ with the transverse momentum range $0.2 < p_{\rm T} < 5.0$ GeV/$c$. The centrality dependence of correlations between two flow coefficients as well as the correlations between three flow coefficients, both in terms of their second moments, are shown. In addition, a collection of mixed harmonic cumulants involving higher moments of $v_2$ and $v_3$ is measured for the first time, where the characteristic signature of negative, positive and negative signs of four-, six- and eight-particle cumulants are observed, respectively. The measurements are compared to the hydrodynamic calculations using iEBE-VISHNU with AMPT and TRENTo initial conditions. It is shown that the measurements carried out using the LHC Run 2 data in 2015 have the precision to explore the details of initial-state fluctuations and probe the nonlinear hydrodynamic response of $v_2$ and $v_3$ to their corresponding initial anisotropy coefficients $\varepsilon_2$ and $\varepsilon_3$. These new studies on correlations between three flow coefficients as well as correlations between higher moments of two different flow coefficients will pave the way to tighten constraints on initial-state models and help to extract precise information on the dynamic evolution of the hot and dense matter created in heavy-ion collisions at the LHC.

Jet fragmentation transverse momentum ($j_{\rm T}$) distributions are measured in proton-proton (pp) and proton-lead (p-Pb) collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE experiment at the LHC. Jets are reconstructed with the ALICE tracking detectors and electromagnetic calorimeter using the anti-$k_{\rm T}$ algorithm with resolution parameter $R=0.4$ in the pseudorapidity range $|\eta|<0.25$. The $j_{\rm T}$ values are calculated for charged particles inside a fixed cone with a radius $R = 0.4$ around the reconstructed jet axis. The measured $j_{\rm T}$ distributions are compared with a variety of parton-shower models. Herwig and PYTHIA 8 based models describe the data well for the higher $j_{\rm T}$ region, while they underestimate the lower $j_{\rm T}$ region. The $j_{\rm T}$ distributions are further characterised by fitting them with a function composed of an inverse gamma function for higher $j_{\rm T}$ values (called the "wide component"), related to the perturbative component of the fragmentation process, and with a Gaussian for lower $j_{\rm T}$ values (called the "narrow component"), predominantly connected to the hadronisation process. The width of the Gaussian has only a weak dependence on jet transverse momentum, while that of the inverse gamma function increases with increasing jet transverse momentum. For the narrow component, the measured trends are successfully described by all models except for Herwig. For the wide component, Herwig and PYTHIA 8 based models slightly underestimate the data for the higher jet transverse momentum region. These measurements set constraints on models of jet fragmentation and hadronisation.

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 production of prompt D$^{0}$, D$^{+}$, and D$^{*+}$ mesons was measured at midrapidity (|y|< 0.5) in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decay channels and their production yields were measured in central (0-10%) and semicentral (30-50%) collisions. The measurement was performed up to a transverse momentum ($p_{\rm T}$) of 36 or 50 GeV/$c$ depending on the D meson species and the centrality interval. For the first time in Pb-Pb collisions at the LHC, the yield of D$^0$ mesons was measured down to $p_{\rm T}$ = 0, which allowed a model-independent determination of the $p_{\rm T}$-integrated yield per unit of rapidity (d$N$/d$y$). A maximum suppression by a factor 5 and 2.5 was observed with the nuclear modification factor ($R_{\rm AA}$) of prompt D mesons at $p_{\rm T}$ = 6-8 GeV/$c$ for the 0-10% and 30-50% centrality classes, respectively. The D-meson $R_{\rm AA}$ is compared with that of charged pions, charged hadrons, and J/$\psi$ mesons as well as with theoretical predictions. The analysis of the agreement between the measured $R_{\rm AA}$, elliptic ($v_2$) and triangular ($v_3$) flow, and the model predictions allowed us to constrain the charm spatial diffusion coefficient $D_s$. Furthermore the comparison of $R_{\rm AA}$ and $v_2$ with different implementations of the same models provides an important insight into the role of radiative energy loss as well as charm quark recombination in the hadronisation mechanisms.

Correlations between mean transverse momentum $[p_{\rm T}]$ and anisotropic flow coefficients $v_{\rm 2}$ or $v_{\rm 3}$ are measured as a function of centrality in Pb$-$Pb and Xe$-$Xe collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between $[p_{\rm T}]$, $v_{\rm 2}$, and $v_{\rm 3}$ is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and $\rm T_{R}ENTo$ initial-state shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the $\rm T_{R}ENTo$ based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the $\rm T_{R}ENTo$ initial state model but with different parameter settings, fail to describe the measurements. As the correlations between $[p_{\rm T}]$ and $v_{\rm n}$ are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark$-$gluon plasma recreated in relativistic heavy-ion collisions.

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.

The energy deposited at very forward rapidities (very forward energy) is a powerful tool for characterising proton fragmentation in pp and p$-$Pb collisions. The correlation of very forward energy with particle production at midrapidity provides direct insights into the initial stages and the subsequent evolution of the collision. Furthermore, the correlation with the production of particles with large transverse momenta at midrapidity provides information complementary to the measurements of the underlying event, which are usually interpreted in the framework of models implementing centrality-dependent multiple parton interactions. Results about very forward energy, measured by the ALICE zero degree calorimeters (ZDCs), and its dependence on the activity measured at midrapidity in pp collisions at $\sqrt{s}=13$ TeV and in p$-$Pb collisions at $\sqrt{s_{\rm{NN}}}=8.16$ TeV are discussed. The measurements performed in pp collisions are compared with the expectations of three hadronic interaction event generators: PYTHIA 6 (Perugia 2011 tune), PYTHIA 8 (Monash tune), and EPOS LHC. These results provide new constraints on the validity of models in describing the beam remnants at very forward rapidities, where perturbative QCD cannot be used.

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.

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.

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.

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.

This article presents new measurements of the fragmentation properties of jets in both proton-proton (pp) and heavy-ion collisions with the ALICE experiment at the LHC. We report distributions of the fraction $z_r$ of transverse momentum carried by subjets of radius $r$ within jets of radius $R$. Charged-particle jets are reconstructed at midrapidity using the anti-$k_{\rm{T}}$ algorithm with jet radius $R=0.4$, and subjets are reconstructed by reclustering the jet constituents using the anti-$k_{\rm{T}}$ algorithm with radii $r=0.1$ and $r=0.2$. In pp collisions, we measure both the inclusive and leading subjet distributions. We compare these measurements to perturbative calculations at next-to-leading logarithmic accuracy, which suggest a large impact of threshold resummation and hadronization effects on the $z_r$ distribution. In heavy-ion collisions, we measure the leading subjet distributions, which allow access to a region of harder jet fragmentation than has been probed by previous measurements of jet quenching via hadron fragmentation distributions. The $z_r$ distributions enable extraction of the parton-to-subjet fragmentation function and allow for tests of the universality of jet fragmentation functions in the quark-gluon plasma (QGP). We find no significant modification of $z_r$ distributions in Pb-Pb compared to pp collisions. However, the distributions are also consistent with a hardening trend for $z_r<0.95$, as predicted by several jet quenching models. As $z_r \rightarrow 1$ our results indicate that any such hardening effects cease, exposing qualitatively new possibilities to disentangle competing jet quenching mechanisms. By comparing our results to theoretical calculations based on an independent extraction of the parton-to-jet fragmentation function, we find consistency with the universality of jet fragmentation and no indication of factorization breaking in the QGP.

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.

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.

A measurement of the inclusive b-jet production cross section is presented in pp and p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, using data collected with the ALICE detector at the LHC. The jets were reconstructed in the central rapidity region $|\eta|<0.5$ from charged particles using the anti-$k_{\rm T}$ algorithm with resolution parameter $R=0.4$. Identification of b jets exploits the long lifetime of b hadrons, using the properties of secondary vertices and impact parameter distributions. The $p_{\rm T}$-differential inclusive production cross section of b jets, as well as the corresponding inclusive b-jet fraction, are reported for pp and p-Pb collisions in the jet transverse momentum range $10 \le p_{\text{T, ch jet}} \le 100$ GeV/$c$, together with the nuclear modification factor, $R_{\rm pPb}^{\text{b-jet}}$. The analysis thus extends the lower $p_{\rm T}$ limit of b-jet measurements at the LHC. The nuclear modification factor is found to be consistent with unity, indicating that the production of b jets in p-Pb at $\sqrt{s_{\rm NN}} = 5.02$ TeV is not affected by cold nuclear matter effects within the current precision. The measurements are well reproduced by POWHEG NLO pQCD calculations with PYTHIA fragmentation.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

The production of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons in proton-proton (pp) and lead-lead (Pb-Pb) collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV has been measured using the ALICE detector at the Large Hadron Collider (LHC). The transverse momentum ($p_{\mathrm{T}}$) distributions of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons have been measured at midrapidity $(|y|<0.5)$ up to $p_{\mathrm{T}} = 20$ GeV$/c$ in inelastic pp collisions and for several Pb-Pb collision centralities. The collision centrality and collision energy dependence of the average transverse momenta agree with the radial flow scenario observed with stable hadrons, showing that the effect is stronger for more central collisions and higher collision energies. The $\mathrm{K^{*0}/K}$ ratio is found to be suppressed in Pb-Pb collisions relative to pp collisions: this indicates a loss of the measured K$^{*}(892)^{0}$ signal due to rescattering of its decay products in the hadronic phase. In contrast, for the longer-lived $\phi(1020)$ mesons, no such suppression is observed. The nuclear modification factors ($R_{\rm AA}$) of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons are calculated using pp reference spectra at the same collision energy. In central Pb-Pb collisions for $p_{\rm T} > 8$ GeV$/c$, the $R_{\rm AA}$ values of K$^{*}(892)^{0}$ and $\phi(1020)$ are below unity and observed to be similar to those of pions, kaons, and (anti)protons. The $R_{\rm AA}$ values at high $p_{\mathrm T}$ ($>$~8 GeV$/c$) for K$^{*}(892)^{0}$ and $\phi(1020)$ mesons are in agreement within uncertainties for $\sqrt{s_\mathrm{NN}} = 5.02$ and 2.76 TeV.

The coherent photoproduction of $\rm{J/\psi}$ and $\rm{\psi'}$ mesons was measured in ultra-peripheral Pb-Pb collisions at a center-of-mass energy $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV with the ALICE detector. Charmonia are detected in the central rapidity region for events where the hadronic interactions are strongly suppressed. The $\rm{J/\psi}$ is reconstructed using the dilepton ($l^{+} l^{-}$) and proton-antiproton decay channels, while for the $\rm{\psi'}$, the dilepton and the $l^{+} l^{-} \pi^{+} \pi^{-}$ decay channels are studied. The analysis is based on an event sample corresponding to an integrated luminosity of about 233 ${\mu b}^{-1}$. The results are compared with theoretical models for coherent $\rm{J/\psi}$ and $\rm{\psi'}$ photoproduction. The coherent cross section is found to be in a good agreement with models incorporating moderate nuclear gluon shadowing of about 0.65 at a Bjorken-$x$ of around $6\times 10^{-4}$, such as the EPS09 parametrization, however none of the models is able to fully describe the rapidity dependence of the coherent $\rm{J/\psi}$ cross section including ALICE measurements at forward rapidity. The ratio of $\rm{\psi'}$ to $\rm{J/\psi}$ coherent photoproduction cross sections was also measured and found to be consistent with the one for photoproduction off protons.

The production of $\Upsilon$ mesons in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5 TeV is measured with the muon spectrometer of the ALICE detector at the LHC. The yields as well as the nuclear modification factors are determined in the forward rapidity region $2.5<y<4.0$, as a function of rapidity, transverse momentum and collision centrality. The results show that the production of the $\Upsilon$(1S) meson is suppressed by a factor of about three with respect to the production in proton-proton collisions. For the first time, a significant signal for the $\Upsilon$(2S) meson is observed at forward rapidity, indicating a suppression stronger by about a factor 2-3 with respect to the ground state. The measurements are compared with transport, hydrodynamic, comover and statistical hadronisation model 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.

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.

This article presents groomed jet substructure measurements in pp and Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector. The Soft Drop grooming algorithm provides access to the hard parton splittings inside a jet by removing soft wide-angle radiation. We report the groomed jet momentum splitting fraction, $z_{\rm g}$, and the (scaled) groomed jet radius, $\theta_{\rm g}$. Charged-particle jets are reconstructed at midrapidity using the anti-kT algorithm with resolution parameters $R = 0.2$ and $R = 0.4$. In heavy-ion collisions, the large underlying event poses a challenge for the reconstruction of groomed jet observables, since fluctuations in the background can cause groomed parton splittings to be misidentified. By using strong grooming conditions to reduce this background, we report these observables fully corrected for detector effects and background fluctuations for the first time. A narrowing of the $\theta_{\rm g}$ distribution in Pb$-$Pb collisions compared to pp collisions is seen, which provides direct evidence of the modification of the angular structure of jets in the quark$-$gluon plasma. No significant modification of the $z_{\rm g}$ distribution in Pb$-$Pb collisions compared to pp collisions is observed. These results are compared with a variety of theoretical models of jet quenching, and provide constraints on jet energy-loss mechanisms and coherence effects in the quark$-$gluon plasma.

A measurement of inclusive, prompt, and non-prompt ${\rm J}/\psi$ production in p$-$Pb collisions at a nucleon$-$nucleon centre-of-mass energy $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV is presented. The inclusive ${\rm J}/\psi$ mesons are reconstructed in the dielectron decay channel at midrapidity down to a transverse momentum $p_{\rm T} = 0$. The inclusive ${\rm J}/\psi$ nuclear modification factor $R_{\rm pPb}$ is calculated by comparing the new results in p$-$Pb collisions to a recently measured proton$-$proton reference at the same centre-of-mass energy. Non-prompt ${\rm J}/\psi$ mesons, which originate from the decay of beauty hadrons, are separated from promptly produced ${\rm J}/\psi$ on a statistical basis for $p_{\rm T}$ larger than 1.0 GeV/$c$. These results are based on the data sample collected by the ALICE detector during the 2016 LHC p$-$Pb run, corresponding to an integrated luminosity ${\cal L}_{\rm int} = 292 \pm 11 \; {\rm \mu b}^{-1}$, which is six times larger than the previous publications. The total uncertainty on the $p_{\rm T}$-integrated inclusive ${\rm J}/\psi$ and non-prompt ${\rm J}/\psi$ cross section are reduced by a factor 1.7 and 2.2, respectively. The measured cross sections and $R_{\rm pPb}$ are compared with theoretical models that include various combinations of cold nuclear matter effects. From the non-prompt ${\rm J}/\psi$ production cross section, the ${\rm b\overline{b}}$ production cross section at midrapidity, $\mathrm{d}\sigma_{\rm b\overline{b}}/\mathrm{d}y$, and the total cross section extrapolated over full phase space, $\sigma_{\rm b\overline{b}}$, are derived.

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 of J/$\psi$ is measured at midrapidity ($|y|<0.9$) in proton-proton collisions at $\sqrt{s}$ = 5.02 and 13 TeV, through the dielectron decay channel, using the ALICE detector at the Large Hadron Collider. The data sets used for the analyses correspond to integrated luminosities of $\mathcal{L}_{\rm int}$ = 19.4 $\pm$ 0.4 nb$^{-1}$ and $\mathcal{L}_{\rm int}$ = 32.2 $\pm$ 0.5 nb$^{-1}$ at $\sqrt{s}$ = 5.02 and 13 TeV, respectively. The fraction of non-prompt J/$\psi$ mesons, i.e. those originating from the decay of beauty hadrons, is measured down to a transverse momentum $p_{\rm T}$ = 2 GeV/$c$ (1 GeV/$c$) at $\sqrt{s}$ = 5.02 TeV (13 TeV). The $p_{\rm T}$ and rapidity ($y$) differential cross sections, as well as the corresponding values integrated over $p_{\rm T}$ and $y$, are carried out separately for prompt and non-prompt J/$\psi$ mesons. The results are compared with measurements from other experiments and theoretical calculations based on quantum chromodynamics (QCD). The shapes of the $p_{\rm T}$ and $y$ distributions of beauty quarks predicted by state-of-the-art perturbative QCD models are used to extrapolate an estimate of the $\rm b\overline{b}$ pair cross section at midrapidity and in the total phase space. The total $\rm b\overline{b}$ cross sections are found to be $\sigma_{\rm b \overline{\rm b}} = 541 \pm 45 (\rm stat.) \pm 69 (\rm syst.)_{-12}^{+10} (\rm extr.)~{\rm \mu b}$ and $\sigma_{\rm b \overline{\rm b}}~=~218 \pm 37 (\rm stat.) \pm 31 (\rm syst.)_{-9.1}^{+8.2} (\rm extr.)~{\rm \mu b}$ at $\sqrt{s}$ = 13 and 5.02 TeV, respectively. The value obtained from the combination of ALICE and LHCb measurements in pp collisions at $\sqrt{s}$ = 13 TeV is also provided.

The polarization of the $\Lambda$ and $\overline\Lambda$ hyperons along the beam ($z$) direction, $P_{\rm z}$, has been measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02TeV recorded with ALICE at the Large Hadron Collider (LHC). The main contribution to $P_{\rm z}$ comes from elliptic flow induced vorticity and can be characterized by the second Fourier sine coefficient $P_{\rm z,s2} = \langle P_{\rm z} \sin(2\varphi - 2 \Psi_{\rm 2}) \rangle$, where $\varphi$ is the hyperon azimuthal emission angle, and $\Psi_{\rm 2}$ is the elliptic flow plane angle. We report the measurement of $P_{\rm z,\,{\rm s2}}$ for different collision centralities, and in the 30-50% centrality interval as a function of the hyperon transverse momentum and rapidity. The $P_{\rm z,\,{\rm s2}}$ is positive similarly as measured by the STAR Collaboration in Au-Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV, with somewhat smaller amplitude in the semi-central collisions. This is the first experimental evidence of a non-zero hyperon $P_{\rm z}$ in Pb-Pb collisions at the LHC. The comparison of the measured $P_{\rm z,\,{\rm s2}}$ with the hydrodynamic model calculations shows sensitivity to the competing contributions from thermal and the recently found shear induced vorticity, as well as to whether the polarization is acquired at the quark-gluon plasma or the hadronic phase.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

The production yield and angular anisotropy of prompt ${D_s^+}$ mesons were measured as a function of transverse momentum ($p_{ T}$) in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{ NN}} = 5.02$ TeV collected with the ALICE detector at the LHC. ${D_s^+}$ mesons and their charge conjugates were reconstructed at midrapidity ($|y|<0.5$) from their hadronic decay channel ${D_s^+ \to \phi \pi^+}$, with ${\phi \to K^-K^+}$, in the $p_{ T}$ intervals $2<p_{ T}<50$ GeV/$c$ and $2<p_{ T}<36$ GeV/$c$ for the 0-10% and 30-50% centrality intervals. For $p_{ T}>10$ GeV/$c$, the measured ${D_s^+}$-meson nuclear modification factor $R_{ AA}$ is consistent with the one of non-strange D mesons within uncertainties, while at lower $p_{ T}$ a hint for a ${D_s^+}$-meson $R_{ AA}$ larger than that of non-strange D mesons is seen. The enhanced production of ${D_s^+}$ relative to non-strange D mesons is also studied by comparing the $p_{ T}$-dependent ${D_s^+/D^0}$ production yield ratios in Pb-Pb and in pp collisions. The ratio measured in Pb-Pb collisions is found to be on average higher than that in pp collisions in the interval $2<p_{ T} <8$ GeV/$c$ with a significance of 2.3$\sigma$ and 2.4$\sigma$ for the 0-10% and 30-50% centrality intervals. The azimuthal anisotropy coefficient $v_2$ of prompt ${D_s^+}$ mesons was measured in Pb-Pb collisions in the 30-50% centrality interval and is found to be compatible with that of non-strange D mesons. The main features of the measured $R_{ AA}$, ${D_s^+/D^0}$ ratio, and $v_2$ as a function of $p_{ T}$ are described by theoretical calculations of charm-quark transport in a hydrodynamically expanding quark-gluon plasma including hadronisation via charm-quark recombination with light quarks from the medium. The $p_{ T}$-integrated production yield of ${D_s^+}$ mesons is compatible with the prediction of the statistical hadronisation model.

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

The production of prompt $\rm \Lambda_{\rm c}^{+}$ baryons at midrapidity ($|y|<0.5$) was measured in central (0-10%) and mid-central (30-50%) Pb-Pb collisions at the center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector. The results are more precise, more differential in centrality, and reach much lower transverse momentum ($p_{\rm T}=1$ GeV/$c$) with respect to previous measurements performed by the ALICE, STAR, and CMS Collaborations in nucleus-nucleus collisions, allowing for an extrapolation down to $p_{\rm T}=0$. The $p_{\rm T}$-differential $\rm \Lambda_{\rm c}^{+}$/D$^0$ production ratio is enhanced with respect to the pp measurement for $4<p_{\rm T}<8$ GeV/$c$ by 3.7 standard deviations ($\sigma$), while the $p_{\rm T}$-integrated ratios are compatible within 1$\sigma$. The observed trend is similar to that observed in the strange sector for the $\Lambda/$K$^0_{\rm S}$ ratio. Model calculations including coalescence or statistical hadronization for charm-hadron formation are compared with the data.

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 transverse momentum ($p_{\rm T}$) differential cross section of the charm-strange baryon $\Xi^0_{\rm c}$ is measured at midrapidity ($|y|<$ 0.5) via its semileptonic decay into ${\rm e^{+}}\Xi^{-}\nu_{\rm e}$ in pp collisions at $\sqrt{s}$ = 5.02 TeV with the ALICE detector at the LHC. The ratio of the $p_{\rm T}$-differential $\Xi^0_{\rm c}$-baryon and ${\rm D^0}$-meson production cross sections is also reported. The measurements are compared with simulations with different tunes of the PYTHIA 8 event generator, with predictions from a statistical hadronisation model (SHM) with a largely augmented set of charm-baryon states beyond the current lists of the Particle Data Group, and with models including hadronisation via quark coalescence. The $p_{\rm T}$-integrated cross section of prompt $\Xi^0_{\rm c}$-baryon production at midrapidity is also reported, which is used to calculate the baryon-to-meson ratio $\Xi^0_{\rm c}/{\rm D^0} = 0.20 \pm 0.04~{\rm (stat.)} ^{+0.08}_{-0.07}~{\rm (syst.)}$. These results provide an additional indication of a modification of the charm fragmentation from $\rm e^+e^-$ and $\rm e^{-}p$ collisions to pp collisions.

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

The transverse-momentum ($p_{\rm T}$) spectra of ${\rm K}^*(892)^0$ and $\phi(1020)$ measured with the ALICE detector up to $p_{\rm T}$ = 16 GeV/$c$ in the rapidity range $-1.2 < y < 0.3$, in p-Pb collisions at the center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $p_{\rm T}$ distributions show a dependence on both multiplicity and rapidity at low $p_{\rm T}$ whereas no significant dependence is observed at high $p_{\rm T}$. A rapidity dependence is observed in the $p_{\rm T}$-integrated yield (d$N$/d$y$), whereas the mean transverse momentum ($\left< p_{\rm T} \right>$) shows a flat behavior as a function of rapidity. The rapidity asymmetry ($Y_{\rm asym}$) at low $p_{\rm T}$ ( < 5 GeV/$c$) is more significant for higher multiplicity classes. At high $p_{\rm T}$, no significant rapidity asymmetry is observed in any of the multiplicity classes. Both ${\rm K}^*(892)^0$ and $\phi(1020)$ show similar $Y_{\rm asym}$. The nuclear modification factor ($Q_{\rm CP}$) as a function of $p_{\rm T}$ shows a Cronin-like enhancement at intermediate $p_{\rm T}$, which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $p_{\rm T}$ (> 5 GeV/$c$), the $Q_{\rm CP}$ values are greater than unity and no significant rapidity dependence is observed.