We report measurements of the production of prompt D$^0$, D$^+$, D$^{*+}$ and D$^+_{\rm s}$ mesons in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV, in the centrality classes 0-10%, 30-50% and 60-80%. The D-meson production yields are measured at mid-rapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm T}$). The $p_{\rm T}$ intervals covered in central collisions are: $1<p_{\rm T}<50$ Gev/$c$ for D$^0$, $2<p_{\rm T}<50$ GeV/$c$ for D$^+$, $3<p_{\rm T}<50$ GeV/$c$ for D$^{*+}$, and $4<p_{\rm T}<16$ GeV/$c$ for D$^+_{\rm s}$ mesons. The nuclear modification factors ($R_{\rm AA}$) for non-strange D mesons (D$^0$, D$^+$, D$^{*+}$) show minimum values of about 0.2 for $p_{\rm T}$ = 6-10 GeV/$c$ in the most central collisions and are compatible within uncertainties with those measured at $\sqrt{s_{\rm NN}}=2.76$ TeV. For D$^+_{\rm s}$ mesons, the values of $R_{\rm AA}$ are larger than those of non-strange D mesons, but compatible within uncertainties. In central collisions the average $R_{\rm AA}$ of non-strange D mesons is compatible with that of charged particles for $p_{\rm T} > 8$ GeV/$c$, while it is larger at lower $p_{\rm T}$. The nuclear modification factors for strange and non-strange D mesons are also compared to theoretical models with different implementations of in-medium energy loss.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.

Measurements of anisotropic flow coefficients with two- and multi-particle cumulants for inclusive charged particles in Pb-Pb collisions at $\sqrt{{\textit s}_\text{NN}} = 5.02$ and 2.76 TeV are reported in the pseudorapidity range $|\eta| < 0.8$ and transverse momentum $0.2 < p_\text{T} < 50$ GeV/$c$. The full data sample collected by the ALICE detector in 2015 (2010), corresponding to an integrated luminosity of 12.7 (2.0) $\mu$b$^{-1}$ in the centrality range 0-80%, is analysed. Flow coefficients up to the sixth flow harmonic ($v_6$) are reported and a detailed comparison among results at the two energies is carried out. The $p_\text{T}$ dependence of anisotropic flow coefficients and its evolution with respect to centrality and harmonic number $n$ are investigated. An approximate power-law scaling of the form $v_n(p_\text{T}) \sim p_\text{T}^{n/3}$ is observed for all flow harmonics at low $p_\text{T}$ ($0.2 < p_\text{T} < 3$ GeV/$c$). At the same time, the ratios $v_n/v_m^{n/m}$ are observed to be essentially independent of $p_\text{T}$ for most centralities up to about $p_\text{T} = 10$ GeV/$c$. Analysing the differences among higher-order cumulants of elliptic flow ($v_2$), which have different sensitivities to flow fluctuations, a measurement of the standardised skewness of the event-by-event $v_2$ distribution $P(v_2)$ is reported and constraints on its higher moments are provided. The Elliptic Power distribution is used to parametrise $P(v_2)$, extracting its parameters from fits to cumulants. The measurements are compared to different model predictions in order to discriminate among initial-state models and to constrain the temperature dependence of the shear viscosity to entropy-density ratio.

The pseudorapidity density of charged particles, $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$, in p-Pb collisions has been measured at a centre-of-mass energy per nucleon-nucleon pair of $\sqrt{s_{\rm{NN}}}$ = 8.16 TeV at mid-pseudorapidity for non-single-diffractive events. The results cover 3.6 units of pseudorapidity, $|\eta|<1.8$. The $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ value is $19.1\pm0.7$ at $|\eta|<0.5$. This quantity divided by $\langle N_\rm{part} \rangle/2$, is $4.73\pm0.20$, which is 9.5% higher than the corresponding value for p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV. Measurements are compared with models based on different mechanisms for particle production. All models agree within uncertainties with data in the Pb-going side, while HIJING overestimates, showing a symmetric behaviour, and EPOS underestimates the p-going side of the $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ distribution. Saturation-based models reproduce the distributions well for $\eta>-1.3$. The $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ is also measured for different centrality estimators, based both on the charged-particle multiplicity and on the energy deposited in the Zero-Degree Calorimeters. A study of the implications of the large multiplicity fluctuations due to the small number of participants for systems like p-Pb in the centrality calculation for multiplicity-based estimators is discussed, demonstrating the advantages of determining the centrality with energy deposited near beam rapidity.

The transverse structure of jets was studied via jet fragmentation transverse momentum ($j_{\rm{T}}$) distributions, obtained using two-particle correlations in proton-proton and proton-lead collisions, measured with the ALICE experiment at the LHC. The highest transverse momentum particle in each event is used as the trigger particle and the region $3 < p_{\rm{Tt}} < 15$ GeV/$c$ is explored in this study. The measured distributions show a clear narrow Gaussian component and a wide non-Gaussian one. Based on Pythia simulations, the narrow component can be related to non-perturbative hadronization and the wide component to quantum chromodynamical splitting. The width of the narrow component shows a weak dependence on the transverse momentum of the trigger particle, in agreement with the expectation of universality of the hadronization process. On the other hand, the width of the wide component shows a rising trend suggesting increased branching for higher transverse momentum. The results obtained in pp collisions at $\sqrt{s}$ = 7 TeV and in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV are compatible within uncertainties and hence no significant cold nuclear matter effects are observed. The results are compared to previous measurements from CCOR and PHENIX as well as to Pythia 8 and Herwig 7 simulations.

Azimuthally-differential femtoscopic measurements, being sensitive to spatio-temporal characteristics of the source as well as to the collective velocity fields at freeze out, provide very important information on the nature and dynamics of the system evolution. While the HBT radii oscillations relative to the second harmonic event plane measured recently reflect mostly the spatial geometry of the source, model studies have shown that the HBT radii oscillations relative to the third harmonic event plane are predominantly defined by the velocity fields. In this Letter, we present the first results on azimuthally-differential pion femtoscopy relative to the third harmonic event plane as a function of the pion pair transverse momentum $k_{\rm T}$ for different collision centralities in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. We find that the $R_{\rm side}$ and $R_{\rm out}$ radii, which characterize the pion source size in the directions perpendicular and parallel to the pion transverse momentum, oscillate in phase relative to the third harmonic event plane, similar to the results from 3+1D hydrodynamical calculations. The observed radii oscillations unambiguously signal a collective expansion and anisotropy in the velocity fields. A comparison of the measured radii oscillations with the Blast-Wave model calculations indicate that the initial state triangularity is washed-out at freeze out.

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 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 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 correlations between event-by-event fluctuations of anisotropic flow harmonic amplitudes have been measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV with the ALICE detector at the LHC. The results are reported in terms of multiparticle correlation observables dubbed Symmetric Cumulants. These observables are robust against biases originating from nonflow effects. The centrality dependence of correlations between the higher order harmonics (the quadrangular $v_4$ and pentagonal $v_5$ flow) and the lower order harmonics (the elliptic $v_2$ and triangular $v_3$ flow) is presented. The transverse momentum dependence of correlations between $v_3$ and $v_2$ and between $v_4$ and $v_2$ is also reported. The results are compared to calculations from viscous hydrodynamics and A Multi-Phase Transport ({AMPT}) model calculations. The comparisons to viscous hydrodynamic models demonstrate that the different order harmonic correlations respond differently to the initial conditions and the temperature dependence of the ratio of shear viscosity to entropy density ($\eta/s$). A small average value of $\eta/s$ is favored independent of the specific choice of initial conditions in the models. The calculations with the AMPT initial conditions yield results closest to the measurements. Correlations between the magnitudes of $v_2$, $v_3$ and $v_4$ show moderate $p_{\rm T}$ dependence in mid-central collisions. Together with existing measurements of individual flow harmonics, the presented results provide further constraints on the initial conditions and the transport properties of the system produced in heavy-ion collisions.