Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

The anisotropy of the azimuthal distributions of charged particles produced in PbPb collisions with a nucleon-nucleon center-of-mass energy of 2.76 TeV is studied with the CMS experiment at the LHC. The elliptic anisotropy parameter defined as the second coefficient in a Fourier expansion of the particle invariant yields, is extracted using the event-plane method, two- and four-particle cumulants, and Lee--Yang zeros. The anisotropy is presented as a function of transverse momentum (pt), pseudorapidity (eta) over a broad kinematic range: 0.3 < pt < 20 GeV, abs(eta) < 2.4, and in 12 classes of collision centrality from 0 to 80%. The results are compared to those obtained at lower center-of-mass energies, and various scaling behaviors are examined. When scaled by the geometric eccentricity of the collision zone, the elliptic anisotropy is found to obey a universal scaling with the transverse particle density for different collision systems and center-of-mass energies.

Spectra of identified charged hadrons are measured in pPb collisions with the CMS detector at the LHC at sqrt(sNN) = 5.02 TeV. Charged pions, kaons, and protons in the transverse-momentum range pt approximately 0.1-1.7 GeV and laboratory rapidity abs(y) < 1 are identified via their energy loss in the silicon tracker. The average pt increases with particle mass and the charged multiplicity of the event. The increase of the average pt with charged multiplicity is greater for heavier hadrons. Comparisons to Monte Carlo event generators reveal that EPOS LHC, which incorporates additional hydrodynamic evolution of the created system, is able to reproduce most of the data features, unlike HIJING and AMPT. The pt spectra and integrated yields are also compared to those measured in pp and PbPb collisions at various energies. The average transverse momentum and particle ratio measurements indicate that particle production at LHC energies is strongly correlated with event particle multiplicity.

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.

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 quark-gluon plasma is studied via medium-induced changes to correlations between jets and charged particles in PbPb collisions compared to pp reference data. This analysis uses data sets from PbPb and pp collisions with integrated luminosities of 166 inverse microbarns and 5.3 inverse picobarns, respectively, collected at sqrt(s[NN]) = 2.76 TeV. The angular distributions of charged particles are studied as a function of relative pseudorapidity (Delta eta) and relative azimuthal angle (Delta phi) with respect to reconstructed jet directions. Charged particles are correlated with all jets with transverse momentum (pt) above 120 GeV, and with the leading and subleading jets (the highest and second-highest in pt, respectively) in a selection of back-to-back dijet events. Modifications in PbPb data relative to pp reference data are characterized as a function of PbPb collision centrality and charged particle pt. A centrality-dependent excess of low-pt particles is present for all jets studied, and is most pronounced in the most central events. This excess of low-pt particles follows a Gaussian-like distribution around the jet axis, and extends to large relative angles of Delta eta approximately 1 and Delta phi approximately 1.

Results on two-particle angular correlations for charged particles produced in pp collisions at a center-of-mass energy of 13 TeV are presented. The data were taken with the CMS detector at the LHC and correspond to an integrated luminosity of about 270 inverse nanobarns. The correlations are studied over a broad range of pseudorapidity (abs(eta) < 2.4) and over the full azimuth (phi) as a function of charged particle multiplicity and transverse momentum (pt). In high-multiplicity events, a long-range (abs(Delta eta) > 2.0), near-side (Delta phi approximately 0) structure emerges in the two-particle Delta eta-Delta phi correlation functions. The magnitude of the correlation exhibits a pronounced maximum in the range 1.0 < pt < 2.0 GeV/c and an approximately linear increase with the charged particle multiplicity, with an overall correlation strength similar to that found in earlier pp data at sqrt(s) = 7 TeV. The present measurement extends the study of near-side long-range correlations up to charged particle multiplicities of N[ch] approximately 180, a region so far unexplored in pp collisions. The observed long-range correlations are compared to those seen in pp, pPb, and PbPb collisions at lower collision energies.

Measurements of strange hadron (K0s, Lambda + anti-Lambda, and Xi+ + Xi-) transverse momentum spectra in pp, pPb, and PbPb collisions are presented over a wide range of rapidity and event charged-particle multiplicity. The data were collected with the CMS detector at the CERN LHC in pp collisions at sqrt(s) = 7 TeV, pPb collisions at sqrt(s[NN]) = 5.02 TeV, and PbPb collisions at sqrt(s[NN]) = 2.76 TeV. The average transverse kinetic energy is found to increase with multiplicity, at a faster rate for heavier strange particle species in all systems. At similar multiplicities, the difference in average transverse kinetic energy between different particle species is observed to be larger for pp and pPb events than for PbPb events. In pPb collisions, the average transverse kinetic energy is found to be slightly larger in the Pb-going direction than in the p-going direction for events with large multiplicity. The spectra are compared to models motivated by hydrodynamics.

Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

For the first time a principle-component analysis is used to separate out different orthogonal modes of the two-particle correlation matrix from heavy ion collisions. The analysis uses data from sqrt(s[NN]) = 2.76 TeV PbPb and sqrt(s[NN]) = 5.02 TeV pPb collisions collected by the CMS experiment at the LHC. Two-particle azimuthal correlations have been extensively used to study hydrodynamic flow in heavy ion collisions. Recently it has been shown that the expected factorization of two-particle results into a product of the constituent single-particle anisotropies is broken. The new information provided by these modes may shed light on the breakdown of flow factorization in heavy ion collisions. The first two modes ("leading" and "subleading") of two-particle correlations are presented for elliptical and triangular anisotropies in PbPb and pPb collisions as a function of pt over a wide range of event activity. The leading mode is found to be essentially equivalent to the anisotropy harmonic previously extracted from two-particle correlation methods. The subleading mode represents a new experimental observable and is shown to account for a large fraction of the factorization breaking recently observed at high transverse momentum. The principle-component analysis technique has also been applied to multiplicity fluctuations. These also show a subleading mode. The connection of these new results to previous studies of factorization is discussed.