The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.

We report the first results of elliptic ($v_2$), triangular ($v_3$) and quadrangular flow ($v_4$) of charged particles in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE detector at the CERN Large Hadron Collider. The measurements are performed in the central pseudorapidity region $|\eta|<0.8$ and for the transverse momentum range $0.2<p_{\rm T}<5$ GeV/$c$. The anisotropic flow is measured using two-particle correlations with a pseudorapidity gap greater than one unit and with the multi-particle cumulant method. Compared to results from Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV, the anisotropic flow coefficients $v_{2}$, $v_{3}$ and $v_{4}$ are found to increase by ($3.0\pm0.6$)%, ($4.3\pm1.4$)% and ($10.2\pm3.8$)%, respectively, in the centrality range 0-50%. This increase can be attributed mostly to an increase of the average transverse momentum between the two energies. The measurements are found to be compatible with hydrodynamic model calculations. This comparison provides a unique opportunity to test the validity of the hydrodynamic picture and the power to further discriminate between various possibilities for the temperature dependence of shear viscosity to entropy density ratio of the produced matter in heavy-ion collisions at the highest energies.

The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

The inclusive transverse momentum ($p_{\rm T}$) distributions of primary charged particles are measured in the pseudo-rapidity range $|\eta|<0.8$ as a function of event centrality in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV with ALICE at the LHC. The data are presented in the $p_{\rm T}$ range $0.15<p_{\rm T}<50$ GeV/$c$ for nine centrality intervals from 70-80% to 0-5%. The Pb-Pb spectra are presented in terms of the nuclear modification factor $R_{\rm{AA}}$ using a pp reference spectrum measured at the same collision energy. We observe that the suppression of high-$p_{\rm T}$ particles strongly depends on event centrality. In central collisions (0-5%) the yield is most suppressed with $R_{\rm{AA}}\approx0.13$ at $p_{\rm T}=6$-7 GeV/$c$. Above $p_{\rm T}=7$ GeV/$c$, there is a significant rise in the nuclear modification factor, which reaches $R_{\rm{AA}} \approx0.4$ for $p_{\rm T}>30$ GeV/$c$. In peripheral collisions (70-80%), the suppression is weaker with $R_{\rm{AA}} \approx 0.7$ almost independently of $p_{\rm T}$. The measured nuclear modification factors are compared to other measurements and model calculations.

In this Letter, the ALICE Collaboration presents the first measurements of the charged-particle multiplicity density, $\rm{d}N_{\rm{ch}}/\rm{d}\eta$, and total charged-particle multiplicity, $N_{\rm{ch}}^{\rm{tot}}$, in Xe-Xe collisions at a centre-of-mass energy per nucleon--nucleon pair of $\sqrt{s_{\rm NN}}$ = 5.44 TeV. The measurements are performed as a function of collision centrality over a wide pseudorapidity range of $-3.5 < \eta < 5$. The values of $\rm{d}N_{\rm{ch}}/\rm{d}\eta$ at mid-rapidity and $N_{\rm{ch}}^{\rm{tot}}$ for central collisions, normalised to the number of nucleons participating in the collision ($N_{\rm{part}}$) as a function of $\sqrt{s_{\rm NN}}$, follow the trends established in previous heavy-ion measurements. The same quantities are also found to increase as a function of $N_{\rm{part}}$, and up to the 10% most central collisions the trends are the same as the ones observed in Pb-Pb at a similar energy. For more central collisions, the Xe-Xe scaled multiplicities exceed those in Pb-Pb for a similar $N_{\rm{part}}$. The results are compared to phenomenological models and theoretical calculations based on different mechanisms for particle production in nuclear collisions. All considered models describe the data reasonably well within 20%.

The pseudorapidity density of charged particles ($\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$) at mid-rapidity in Pb-Pb collisions has been measured at a center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 5.02 TeV. It increases with centrality and reaches a value of $1943 \pm 54$ in $|\eta|<0.5$ for the 5% most central collisions. A rise in $\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$ as a function of $\sqrt{s_{\rm NN}}$ for the most central collisions is observed, steeper than that observed in proton-proton collisions and following the trend established by measurements at lower energy. The centrality dependence of $\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$ as a function of the average number of participant nucleons, ${\langle N_\mathrm{part} \rangle}$, calculated in a Glauber model, is compared with the previous measurement at lower energy. A constant factor of about 1.2 describes the increase in $\frac{2}{\langle N_\mathrm{part} \rangle}\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$ from $\sqrt{s_{\rm NN}}$ = 2.76 TeV to $\sqrt{s_{\rm NN}}$ = 5.02 TeV for all centrality intervals, within the measured range of 0-80% centrality. The results are also compared to models based on different mechanisms for particle production in nuclear collisions.

The centrality dependence of the charged-particle multiplicity density at mid-rapidity in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is presented. The charged-particle density normalized per participating nucleon pair increases by about a factor 2 from peripheral (70-80%) to central (0-5%) collisions. The centrality dependence is found to be similar to that observed at lower collision energies. The data are compared with models based on different mechanisms for particle production in nuclear collisions.

We present the first wide-range measurement of the charged-particle pseudorapidity density distribution, for different centralities (the 0-5%, 5-10%, 10-20%, and 20-30% most central events) in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV at the LHC. The measurement is performed using the full coverage of the ALICE detectors, $-5.0 < \eta < 5.5$, and employing a special analysis technique based on collisions arising from LHC "satellite" bunches. We present the pseudorapidity density as a function of the number of participating nucleons as well as an extrapolation to the total number of produced charged particles ($N_{\rm ch} = 17165 \pm 772$ for the 0-5% most central collisions). From the measured ${\rm d}N_{\rm ch}/{\rm d}\eta$ distribution we derive the rapidity density distribution, ${\rm d}N_{\rm ch}/{\rm d}y$, under simple assumptions. The rapidity density distribution is found to be significantly wider than the predictions of the Landau model. We assess the validity of longitudinal scaling by comparing to lower energy results from RHIC. Finally the mechanisms of the underlying particle production are discussed based on a comparison with various theoretical models.

We present the charged-particle pseudorapidity density in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02\,\mathrm{Te\kern-.25exV}$ in centrality classes measured by ALICE. The measurement covers a wide pseudorapidity range from $-3.5$ to $5$, which is sufficient for reliable estimates of the total number of charged particles produced in the collisions. For the most central (0-5%) collisions we find $21\,400\pm 1\,300$ while for the most peripheral (80-90%) we find $230\pm 38$. This corresponds to an increase of $(27\pm4)\%$ over the results at $\sqrt{s_{\mathrm{NN}}}=2.76\,\mathrm{Te\kern-.25exV}$ previously reported by ALICE. The energy dependence of the total number of charged particles produced in heavy-ion collisions is found to obey a modified power-law like behaviour. The charged-particle pseudorapidity density of the most central collisions is compared to model calculations --- none of which fully describes the measured distribution. We also present an estimate of the rapidity density of charged particles. The width of that distribution is found to exhibit a remarkable proportionality to the beam rapidity, independent of the collision energy from the top SPS to LHC energies.

This publication describes the methods used to measure the centrality of inelastic Pb-Pb collisions at a center-of-mass energy of 2.76 TeV per colliding nucleon pair with ALICE. The centrality is a key parameter in the study of the properties of QCD matter at extreme temperature and energy density, because it is directly related to the initial overlap region of the colliding nuclei. Geometrical properties of the collision, such as the number of participating nucleons and number of binary nucleon-nucleon collisions, are deduced from a Glauber model with a sharp impact parameter selection, and shown to be consistent with those extracted from the data. The centrality determination provides a tool to compare ALICE measurements with those of other experiments and with theoretical calculations.

The centrality dependence of the charged-particle pseudorapidity density measured with ALICE in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ over a broad pseudorapidity range is presented. This Letter extends the previous results reported by ALICE to more peripheral collisions. No strong change of the charged-particle pseudorapidity density distributions with centrality is observed, and when normalised to the number of participating nucleons in the collisions, the evolution over pseudorapidity with centrality is likewise small. The broad pseudorapidity range allows precise estimates of the total number of produced charged particles which we find to range from $162\pm22$ (syst.) to $17170\pm770$ (syst.) in 80-90% and 0-5 central collisions, respectively. The total charged-particle multiplicity is seen to approximately scale with the number of participating nucleons in the collision. This suggests that hard contributions to the charged-particle multiplicity are limited. The results are compared to models which describe $\mbox{d}N_{\mbox{ch}}/\mbox{d}\eta$ at mid-rapidity in the most central Pb-Pb collisions and it is found that these models do not capture all features of the distributions.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

The spectra of charged particles produced within the pseudorapidity window abs(eta) < 1 at sqrt(s[NN]) = 5.02 TeV are measured using 404 inverse microbarns of PbPb and 27.4 inverse picobarns of pp data collected by the CMS detector at the LHC in 2015. The spectra are presented over the transverse momentum ranges spanning 0.5 < pt < 400 GeV in pp and 0.7 < pt < 400 GeV in PbPb collisions. The corresponding nuclear modification factor, R[AA], is measured in bins of collision centrality. The R[AA] in the 5% most central collisions shows a maximal suppression by a factor of 7-8 in the pt region of 6-9 GeV. This dip is followed by an increase, which continues up to the highest pt measured, and approaches unity in the vicinity of pt = 200 GeV. The R[AA] is compared to theoretical predictions and earlier experimental results at lower collision energies. The newly measured pp spectrum is combined with the pPb spectrum previously published by the CMS Collaboration to construct the pPb nuclear modification factor, R[pA], up to 120 GeV. For pt > 20 GeV, R[pA] exhibits weak momentum dependence and shows a moderate enhancement above unity.

We report the measurements of correlations between event-by-event fluctuations of amplitudes of anisotropic flow harmonics in nucleus-nucleus collisions, obtained for the first time using a new analysis method based on multiparticle cumulants in mixed harmonics. This novel method is robust against systematic biases originating from non-flow effects and by construction any dependence on symmetry planes is eliminated. We demonstrate that correlations of flow harmonics exhibit a better sensitivity to medium properties than the individual flow harmonics. The new measurements are performed in Pb-Pb collisions at the centre-of-mass energy per nucleon pair of $\sqrt{s_{_{\rm NN}}}=2.76$ TeV by the ALICE experiment at the Large Hadron Collider (LHC). The centrality dependence of correlation between event-by-event fluctuations of the elliptic, $v_2$, and quadrangular, $v_4$, flow harmonics, as well as of anti-correlation between $v_2$ and triangular, $v_3$, flow harmonics are presented. The results cover two different regimes of the initial state configurations: geometry-dominated (in mid-central collisions) and fluctuation-dominated (in the most central collisions). Comparisons are made to predictions from MC-Glauber, viscous hydrodynamics, AMPT and HIJING models. Together with the existing measurements of individual flow harmonics the presented results provide further constraints on initial conditions and the transport properties of the system produced in heavy-ion collisions.

A measurement is presented of the charged hadron multiplicity in hadronic PbPb collisions, as a function of pseudorapidity and centrality, at a collision energy of 2.76 TeV per nucleon pair. The data sample is collected using the CMS detector and a minimum-bias trigger, with the CMS solenoid off. The number of charged hadrons is measured both by counting the number of reconstructed particle hits and by forming hit doublets of pairs of layers in the pixel detector. The two methods give consistent results. The charged hadron multiplicity density dN(ch)/d eta, evaluated at eta=0 for head-on collisions, is found to be 1612 +/- 55, where the uncertainty is dominated by systematic effects. Comparisons of these results to previous measurements and to various models are also presented.

The directed flow of charged particles at mid-rapidity is measured in Pb-Pb collisions at $\sqrt{s_\rm{NN}}=2.76$ TeV relative to the collision symmetry plane defined by the spectator nucleons. A negative slope of the rapidity-odd directed flow component with approximately 3 times smaller magnitude than found at the highest RHIC energy is observed. This suggests a smaller longitudinal tilt of the initial system and disfavors the strong fireball rotation predicted for the LHC energies. The rapidity-even directed flow component is measured for the first time with spectators and found to be independent of pseudorapidity with a sign change at transverse momenta $p_{\rm T}$ between $1.2$ and $1.7$ GeV/$c$. Combined with the observation of a vanishing rapidity-even $p_{\rm T}$ shift along the spectator deflection this is strong evidence for dipole-like initial density fluctuations in the overlap zone of the nuclei. Similar trends in the rapidity-even directed flow and the estimate from two-particle correlations at mid-rapidity, which is larger by about a factor of 40, indicate a weak correlation between fluctuating participant and spectator symmetry planes. These observations open new possibilities for investigation of the initial conditions in heavy-ion collisions with spectator nucleons.

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.

We report the first measurement of charged particle elliptic flow in Pb-Pb collisions at 2.76 TeV with the ALICE detector at the CERN Large Hadron Collider. The measurement is performed in the central pseudorapidity region (|$\eta$|<0.8) and transverse momentum range 0.2< $p_{\rm T}$< 5.0 GeV/$c$. The elliptic flow signal v$_2$, measured using the 4-particle correlation method, averaged over transverse momentum and pseudorapidity is 0.087 $\pm$ 0.002 (stat) $\pm$ 0.004 (syst) in the 40-50% centrality class. The differential elliptic flow v$_2(p_{\rm T})$ reaches a maximum of 0.2 near $p_{\rm T}$ = 3 GeV/$c$. Compared to RHIC Au-Au collisions at 200 GeV, the elliptic flow increases by about 30%. Some hydrodynamic model predictions which include viscous corrections are in agreement with the observed increase.

We report on results obtained with the Event Shape Engineering technique applied to Pb-Pb collisions at $\sqrt{s_\rm{NN}}=2.76$ TeV. By selecting events in the same centrality interval, but with very different average flow, different initial state conditions can be studied. We find the effect of the event-shape selection on the elliptic flow coefficient $v_2$ to be almost independent of transverse momentum $p_\rm{T}$, as expected if this effect is due to fluctuations in the initial geometry of the system. Charged hadron, pion, kaon, and proton transverse momentum distributions are found to be harder in events with higher-than-average elliptic flow, indicating an interplay between radial and elliptic flow.

The event-by-event correlations between three flow amplitudes are measured for the first time in Pb--Pb collisions, using higher-order Symmetric Cumulants. We find that different three-harmonic correlations develop during the collective evolution of the medium, when compared with correlations that exist in the initial state. These new results cannot be interpreted in terms of previous lower-order flow measurements, since contributions from two-harmonic correlations are explicitly removed in the new observables. Comparison with Monte Carlo simulations provides new and independent constraints for the initial conditions and system properties of nuclear matter created in heavy-ion collisions.

The second-order azimuthal anisotropy Fourier harmonics, v2, are obtained in pPb and PbPb collisions over a wide pseudorapidity (eta) range based on correlations among six or more charged particles. The pPb data, corresponding to an integrated luminosity of 35 inverse nanobarns, were collected during the 2013 LHC pPb run at a nucleon-nucleon center-of-mass energy of 5.02 TeV by the CMS experiment. A sample of semi-peripheral PbPb collision data at sqrt(s[NN])= 2.76 TeV, corresponding to an integrated luminosity of 2.5 inverse microbarns and covering a similar range of particle multiplicities as the pPb data, is also analyzed for comparison. The six- and eight-particle cumulant and the Lee-Yang zeros methods are used to extract the v2 coefficients, extending previous studies of two- and four-particle correlations. For both the pPb and PbPb systems, the v2 values obtained with correlations among more than four particles are consistent with previously published four-particle results. These data support the interpretation of a collective origin for the previously observed long-range (large Delta[eta]) correlations in both systems. The ratios of v2 values corresponding to correlations including different numbers of particles are compared to theoretical predictions that assume a hydrodynamic behavior of a pPb system dominated by fluctuations in the positions of participant nucleons. These results provide new insights into the multi-particle dynamics of collision systems with a very small overlapping region.

A systematic study of the factorization of long-range azimuthal two-particle correlations into a product of single-particle anisotropies is presented as a function of pt and eta of both particles, and as a function of the particle multiplicity in PbPb and pPb collisions. The data were taken with the CMS detector for PbPb collisions at sqrt(s[NN]) = 2.76 TeV and pPb collisions at sqrt(s[NN]) = 5.02 TeV, covering a very wide range of multiplicity. Factorization is observed to be broken as a function of both particle pt and eta. When measured with particles of different pt, the magnitude of the factorization breakdown for the second Fourier harmonic reaches 20% for very central PbPb collisions but decreases rapidly as the multiplicity decreases. The data are consistent with viscous hydrodynamic predictions, which suggest that the effect of factorization breaking is mainly sensitive to the initial-state conditions rather than to the transport properties (e.g., shear viscosity) of the medium. The factorization breakdown is also computed with particles of different eta. The effect is found to be weakest for mid-central PbPb events but becomes larger for more central or peripheral PbPb collisions, and also for very high-multiplicity pPb collisions. The eta-dependent factorization data provide new insights to the longitudinal evolution of the medium formed in heavy ion collisions.

We report on the first measurement of the triangular $v_3$, quadrangular $v_4$, and pentagonal $v_5$ charged particle flow in Pb-Pb collisions at 2.76 TeV measured with the ALICE detector at the CERN Large Hadron Collider. We show that the triangular flow can be described in terms of the initial spatial anisotropy and its fluctuations, which provides strong constraints on its origin. In the most central events, where the elliptic flow $v_2$ and $v_3$ have similar magnitude, a double peaked structure in the two-particle azimuthal correlations is observed, which is often interpreted as a Mach cone response to fast partons. We show that this structure can be naturally explained from the measured anisotropic flow Fourier coefficients.

Anisotropic flow coefficients, $v_n$, non-linear flow mode coefficients, $\chi_{n,mk}$, and correlations among different symmetry planes, $\rho_{n,mk}$ are measured in Pb-Pb collisions at $\sqrt{s_\rm{NN}}=5.02$ TeV. Results obtained with multi-particle correlations are reported for the transverse momentum interval $0.2<p_\rm{T}<5.0$ GeV/$c$ within the pseudorapidity interval $0.4<|\eta|<0.8$ as a function of collision centrality. The $v_n$ coefficients and $\chi_{n,mk}$ and $\rho_{n,mk}$ are presented up to the ninth and seventh harmonic order, respectively. Calculations suggest that the correlations measured in different symmetry planes and the non-linear flow mode coefficients are dependent on the shear and bulk viscosity to entropy ratios of the medium created in heavy-ion collisions. The comparison between these measurements and those at lower energies and calculations from hydrodynamic models places strong constraints on the initial conditions and transport properties of the system.

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.