The second- and third-order azimuthal anisotropy Fourier harmonics of charged particles produced in pPb collisions, at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV, are studied over a wide range of event multiplicities. Multiparticle correlations are used to isolate global properties stemming from the collision overlap geometry. The second-order "elliptic" harmonic moment is obtained with high precision through four-, six-, and eight-particle correlations and, for the first time, the third-order "triangular" harmonic moment is studied using four-particle correlations. A sample of peripheral PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV that covers a similar range of event multiplicities as the pPb results is also analyzed. Model calculations of initial-state fluctuations in pPb and PbPb collisions can be directly compared to the high precision experimental results. This work provides new insight into the fluctuation-driven origin of the $v_3$ coefficients in pPb and PbPb collisions, and into the dominating overall collision geometry in PbPb collisions at the earliest stages of heavy ion interactions.

Charge-dependent azimuthal correlations of same- and opposite-sign pairs with respect to the second- and third-order event planes have been measured in pPb collisions at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV and PbPb collisions at 5.02 TeV with the CMS experiment at the LHC. The measurement is motivated by the search for the charge separation phenomenon predicted by the chiral magnetic effect (CME) in heavy ion collisions. Three- and two-particle azimuthal correlators are extracted as functions of the pseudorapidity difference, the transverse momentum ($p_\mathrm{t}$) difference, and the $p_\mathrm{t}$ average of same- and opposite-charge pairs in various event multiplicity ranges. The data suggest that the charge-dependent three-particle correlators with respect to the second- and third-order event planes share a common origin, predominantly arising from charge-dependent two-particle azimuthal correlations coupled with an anisotropic flow. The CME is expected to lead to a $v_2$-independent three-particle correlation when the magnetic field is fixed. Using an event shape engineering technique, upper limits on the $v_2$-independent fraction of the three-particle correlator are estimated to be 13% for pPb and 7% for PbPb collisions at 95% confidence level. The results of this analysis, both the dominance of two-particle correlations as a source of the three-particle results and the similarities seen between PbPb and pPb, provide stringent constraints on the origin of charge-dependent three-particle azimuthal correlations and challenge their interpretation as arising from a chiral magnetic effect in heavy ion collisions.

The azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV pPb data are extracted via long-range two-particle correlations as a function of event multiplicity and compared to corresponding results in pp and PbPb collisions. Using a four-particle cumulant technique, $v_n$ correlations are measured for the first time in pp and pPb collisions. The $v_2$ and $v_4$ coefficients are found to be positively correlated in all collision systems. For high multiplicity pPb collisions an anticorrelation of $v_2$ and $v_3$ is observed, with a similar correlation strength as in PbPb data at the same multiplicity. The new correlation results strengthen the case for a common origin of the collectivity seen in pPb and PbPb collisions in the measured multiplicity range.

Despite the f$_0$(980) hadron having been discovered half a century ago, the question about its quark content has not been settled: it might be an ordinary quark-antiquark ($\mathrm{q\bar{q}}$) meson, a tetraquark ($\mathrm{q\bar{q}q\bar{q}}$) exotic state, a kaon-antikaon ($\mathrm{K\bar{K}}$) molecule, or a quark-antiquark-gluon ($\mathrm{q\bar{q}g}$) hybrid. This paper reports strong evidence that the f$_0$(980) state is an ordinary $\mathrm{q\bar{q}}$ meson, inferred from the scaling of elliptic anisotropies ($v_2$) with the number of constituent quarks ($n_\mathrm{q}$), as empirically established using conventional hadrons in relativistic heavy ion collisions. The f$_0$(980) state is reconstructed via its dominant decay channel f$_0$(980) $\to$$\pi^+\pi^-$, in proton-lead collisions recorded by the CMS experiment at the LHC, and its $v_2$ is measured as a function of transverse momentum ($p_\mathrm{T}$). It is found that the $n_q$ = 2 ($\mathrm{q\bar{q}}$ state) hypothesis is favored over $n_q$ = 4 ($\mathrm{q\bar{q}q\bar{q}}$ or $\mathrm{K\bar{K}}$ states) by 7.7, 6.3, or 3.1 standard deviations in the $p_\mathrm{T}$$\lt$ 10, 8, or 6 GeV/$c$ ranges, respectively, and over $n_\mathrm{q}$ = 3 ($\mathrm{q\bar{q}g}$ hybrid state) by 3.5 standard deviations in the $p_\mathrm{T}$$\lt$ 8 GeV/$c$ range. This result represents the first determination of the quark content of the f$_0$(980) state, made possible by using a novel approach, and paves the way for similar studies of other exotic hadron candidates.

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 azimuthal anisotropy of $\Upsilon$(1S) mesons in high-multiplicity proton-lead collisions is studied using data collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 8.16 TeV. The $\Upsilon$(1S) mesons are reconstructed using their dimuon decay channel. The anisotropy is characterized by the second Fourier harmonic coefficients, found using a two-particle correlation technique, in which the $\Upsilon$(1S) mesons are correlated with charged hadrons. A large pseudorapidity gap is used to suppress short-range correlations. Nonflow contamination from the dijet background is removed using a low-multiplicity subtraction method, and the results are presented as a function of $\Upsilon$(1S) transverse momentum. The azimuthal anisotropies are smaller than those found for charmonia in proton-lead collisions at the same collision energy, but are consistent with values found for $\Upsilon$(1S) mesons in lead-lead interactions at a nucleon-nucleon center-of-mass energy of 5.02 TeV.

A measurement of the elliptic flow ($v_2$) of prompt J/$\psi$ mesons in high-multiplicity pPb collisions is reported using data collected by the CMS experiment at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV. Prompt J/$\psi$ mesons decaying into two muons are reconstructed in the rapidity region in the nucleon-nucleon center-of-mass frame ($y_\mathrm{cm}$), corresponding to either $-$2.86 $\lt y_{\mathrm{cm}} <$ $-$1.86 or 0.94 $\lt y_{\mathrm{cm}} <$ 1.94. The average $v_2$ result from the two rapidity ranges is reported over the transverse momentum (p$_\mathrm{T}$) range from 0.2 to 10 GeV. Positive $v_2$ values are observed for the prompt J/$\psi$ meson, as extracted from long-range two-particle correlations with charged hadrons, for 2 $\lt p_{\mathrm{T}}<$ 8 GeV. The prompt J/$\psi$ results are compared with previous CMS measurements of elliptic flow for open charm mesons (D$^0$) and strange hadrons. From these measurements, constraints can be obtained on the collective dynamics of charm quarks produced in high-multiplicity events arising from small systems.

The elliptic azimuthal anisotropy coefficient ($v_2$) is measured for charm (D$^0$) and strange (K$_\mathrm{S}^0$, $\Lambda$, $\Xi^-$, and $\Omega^-$) hadrons, using a data sample of pPb collisions collected by the CMS experiment, at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_\mathrm{NN}}} =$ 8.16 TeV. A significant positive $v_2$ signal from long-range azimuthal correlations is observed for all particle species in high-multiplicity pPb collisions. The measurement represents the first observation of possible long-range collectivity for open heavy flavor hadrons in small systems. The results suggest that charm quarks have a smaller $v_2$ than the lighter quarks, probably reflecting a weaker collective behavior. This effect is not seen in the larger PbPb collision system at $\sqrt{s_{_\mathrm{NN}}} =$ 5.02 TeV, also presented.

Measurements of the second Fourier harmonic coefficient ($v_2$) of the azimuthal distributions of prompt and nonprompt D$^0$ mesons produced in pp and pPb collisions are presented. Nonprompt D$^0$ mesons come from beauty hadron decays. The data samples are collected by the CMS experiment at nucleon-nucleon center-of-mass energies of 13 and 8.16 TeV, respectively. In high multiplicity pp collisions, $v_2$ signals for prompt charm hadrons are reported for the first time, and are found to be comparable to those for light-flavor hadron species over a transverse momentum ($p_\mathrm{T}$) range of 2-6 GeV. Compared at similar event multiplicities, the prompt D$^0$ meson $v_2$ values in pp and pPb collisions are similar in magnitude. The $v_2$ values for open beauty hadrons are extracted for the first time via nonprompt D$^0$ mesons in pPb collisions. For $p_\mathrm{T}$ in the range of 2-5 GeV, the results suggest that $v_2$ for nonprompt D$^0$ mesons are smaller than those for prompt D$^0$ mesons. These new measurements indicate a positive charm hadron $v_2$ in pp collisions and suggest a mass dependence in $v_2$ between charm and beauty hadrons in the pPb system. These results provide insights into the origin of heavy-flavor quark collectivity in small systems.

Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum $p_\mathrm{T}$ and of a geometric variable $\phi^*$ are measured, where $\phi^*$ highly correlates with $p_\mathrm{T}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.