Measurements of the charge-dependent two-particle angular correlation function in proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at$\sqrt{s_\mathrm{NN}}$ = 5.02 TeV are reported. The pPb and PbPb datasets correspond to integrated luminosities of 186\nbinv and 0.607 nb$^{-1}$, respectively, and were collected using the CMS detector at the CERN LHC. The charge-dependent correlations are characterized by balance functions of same- and opposite-sign particle pairs. The balance functions, which contain information about the creation time of charged particle pairs and the development of collectivity, are studied as functions of relative pseudorapidity ($\Delta \eta$) and relative azimuthal angle ($\Delta \phi$), for various multiplicity and transverse momentum ($p_\mathrm{T}$) intervals. A multiplicity dependence of the balance function is observed in $\Delta \eta$ and $\Delta \phi$ for both systems. The width of the balance functions decreases towards high-multiplicity collisions in the momentum region $\lt$2 GeV, for pPb and PbPb results. No multiplicity dependence is observed at higher transverse momentum. The data are compared with HYDJET, HIJING and AMPT generator predictions, none of which capture completely the multiplicity dependence seen in the data.

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.

The hydrodynamic flow-like behavior of charged hadrons in high-energy lead-lead collisions is studied through multiparticle correlations. The elliptic anisotropy values based on different orders of multiparticle cumulants, $v_{2}\{2k\}$, are measured up to the tenth order ($k$ = 5) as functions of the collision centrality at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. A hierarchy is observed between the coefficients, with $v_{2}\{2\} > v_{2}\{4\} \gtrsim v_{2}\{6\} \gtrsim v_{2}\{8\} \gtrsim v_{2}\{10\}$. Based on these results, centrality-dependent moments for the fluctuation-driven event-by-event $v_{2}$ distribution are determined, including the skewness, kurtosis and, for the first time, superskewness. Assuming a hydrodynamic expansion of the produced medium, these moments directly probe the initial-state geometry in high-energy nucleus-nucleus collisions.

Two-particle Bose-Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data sample, containing 4.27 $\times$$10^{9}$ minimum bias events corresponding to an integrated luminosity of 0.607 nb$^{-1}$, was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter ($\alpha$), the Lévy scale parameter ($R$), and the correlation strength parameter ($\lambda$). The source shape, characterized by $\alpha$, is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy $R$ parameter. The $\lambda$ parameter is studied in terms of the core-halo model.

Differential measurements of charged particle azimuthal anisotropy are presented for lead-lead collisions at sqrt(s_NN) = 2.76 TeV with the ATLAS detector at the LHC, based on an integrated luminosity of approximately 8 mb^-1. This anisotropy is characterized via a Fourier expansion of the distribution of charged particles in azimuthal angle (phi), with the coefficients v_n denoting the magnitude of the anisotropy. Significant v_2-v_6 values are obtained as a function of transverse momentum (0.5<pT<20 GeV), pseudorapidity (|eta|<2.5) and centrality using an event plane method. The v_n values for n>=3 are found to vary weakly with both eta and centrality, and their pT dependencies are found to follow an approximate scaling relation, v_n^{1/n}(pT) \propto v_2^{1/2}(pT). A Fourier analysis of the charged particle pair distribution in relative azimuthal angle (Dphi=phi_a-phi_b) is performed to extract the coefficients v_{n,n}=<cos (n Dphi)>. For pairs of charged particles with a large pseudorapidity gap (|Deta=eta_a-eta_b|>2) and one particle with pT<3 GeV, the v_{2,2}-v_{6,6} values are found to factorize as v_{n,n}(pT^a,pT^b) ~ v_n(pT^a)v_n(pT^b) in central and mid-central events. Such factorization suggests that these values of v_{2,2}-v_{6,6} are primarily due to the response of the created matter to the fluctuations in the geometry of the initial state. A detailed study shows that the v_{1,1}(pT^a,pT^b) data are consistent with the combined contributions from a rapidity-even v_1 and global momentum conservation. A two-component fit is used to extract the v_1 contribution. The extracted v_1 is observed to cross zero at pT\sim1.0 GeV, reaches a maximum at 4-5 GeV with a value comparable to that for v_3, and decreases at higher pT.

The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

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.

The distributions of event-by-event harmonic flow coefficients v_n for n=2-4 are measured in sqrt(s_NN)=2.76 TeV Pb+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using charged particles with transverse momentum pT> 0.5 GeV and in the pseudorapidity range |eta|<2.5 in a dataset of approximately 7 ub^-1 recorded in 2010. The shapes of the v_n distributions are described by a two-dimensional Gaussian function for the underlying flow vector in central collisions for v_2 and over most of the measured centrality range for v_3 and v_4. Significant deviations from this function are observed for v_2 in mid-central and peripheral collisions, and a small deviation is observed for v_3 in mid-central collisions. It is shown that the commonly used multi-particle cumulants are insensitive to the deviations for v_2. The v_n distributions are also measured independently for charged particles with 0.5<pT<1 GeV and pT>1 GeV. When these distributions are rescaled to the same mean values, the adjusted shapes are found to be nearly the same for these two pT ranges. The v_n distributions are compared with the eccentricity distributions from two models for the initial collision geometry: a Glauber model and a model that includes corrections to the initial geometry due to gluon saturation effects. Both models fail to describe the experimental data consistently over most of the measured centrality range.

Measurements of two-particle angular correlations between an identified strange hadron (K0S or Lambda/anti-Lambda) and a charged particle, emitted in pPb collisions, are presented over a wide range in pseudorapidity and full azimuth. The data, corresponding to an integrated luminosity of approximately 35 inverse nanobarns, were collected at a nucleon-nucleon center-of-mass energy (sqrt(s[NN])) of 5.02 TeV with the CMS detector at the LHC. The results are compared to semi-peripheral PbPb collision data at sqrt(s[NN]) = 2.76 TeV, covering similar charged-particle multiplicities in the events. The observed azimuthal correlations at large relative pseudorapidity are used to extract the second-order (v[2]) and third-order (v[3]) anisotropy harmonics of K0S and Lambda/anti-Lambda particles. These quantities are studied as a function of the charged-particle multiplicity in the event and the transverse momentum of the particles. For high-multiplicity pPb events, a clear particle species dependence of v[2] and v[3] is observed. For pt < 2 GeV, the v[2] and v[3] values of K0S particles are larger than those of Lambda/anti-Lambda particles at the same pt. This splitting effect between two particle species is found to be stronger in pPb than in PbPb collisions in the same multiplicity range. When divided by the number of constituent quarks and compared at the same transverse kinetic energy per quark, both v[2] and v[3] for K0S particles are observed to be consistent with those for Lambda/anti-Lambda particles at the 10% level in pPb collisions. This consistency extends over a wide range of particle transverse kinetic energy and event multiplicities.

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.