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.

Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

The relative yields of $\Upsilon$ mesons produced in pp and PbPb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV and reconstructed via the dimuon decay channel are measured using data collected by the CMS experiment. Double ratios are formed by comparing the yields of the excited states, $\Upsilon$(2S) and $\Upsilon$(3S), to the ground state, $\Upsilon$(1S), in both PbPb and pp collisions at the same center-of-mass energy. The double ratios, [$\Upsilon$(nS)/$\Upsilon$(1S)]$_\mathrm{PbPb}$ / [$\Upsilon$(nS)/$\Upsilon$(1S)]$_\mathrm{pp}$, are measured to be 0.308 $\pm$ 0.055 (stat) $\pm$ 0.019 (syst) for the $\Upsilon$(2S) and less than 0.26 at 95% confidence level for the $\Upsilon$(3S). No significant $\Upsilon$(3S) signal is found in the PbPb data. The double ratios are studied as a function of collision centrality, as well as dimuon transverse momentum and rapidity. No significant dependencies are observed.

The transverse momentum ($p_\mathrm{t}$) spectrum of prompt D$^0$ mesons and their antiparticles has been measured via the hadronic decay channels D$^0 \to \mathrm{K}^- \pi^+$ and $\overline{\mathrm{D}}^0 \to \mathrm{K}^+ \pi^-$ in pp and PbPb collisions at a centre-of-mass energy of 5.02 TeV per nucleon pair with the CMS detector at the LHC. The measurement is performed in the D$^0$ meson $p_\mathrm{t}$ range of 2-100 GeV and in the rapidity range of $|y| < $1. The pp (PbPb) dataset used for this analysis corresponds to an integrated luminosity of 27.4 pb$^{-1}$ (530 $\mu$b$^{-1}$). The measured D$^0$ meson $p_\mathrm{t}$ spectrum in pp collisions is well described by perturbative QCD calculations. The nuclear modification factor, comparing D$^0$ meson yields in PbPb and pp collisions, was extracted for both minimum-bias and the 10% most central PbPb interactions. For central events, the D$^0$ meson yield in the PbPb collisions is suppressed by a factor of 5-6 compared to the pp reference in the $p_\mathrm{t}$ range of 6-10 GeV. For D$^0$ mesons in the high-$p_\mathrm{t}$ range of 60-100 GeV, a significantly smaller suppression is observed. The results are also compared to theoretical calculations.

The prompt D$^0$ meson azimuthal anisotropy coefficients, $v_2$ and $v_3$, are measured at midrapidity ($|y| < 1.0$) in PbPb collisions at a center-of-mass energy $\sqrt{s_\mathrm{NN}} = $5.02 TeV per nucleon pair with data collected by the CMS experiment. The measurement is performed in the transverse momentum ($p_\mathrm{T}$) range of 1 to 40 GeV/c, for central and midcentral collisions. The $v_2$ coefficient is found to be positive throughout the $p_\mathrm{T}$ range studied. The first measurement of the prompt D$^0$ meson $v_3$ coefficient is performed, and values up to 0.07 are observed for $p_\mathrm{T}$ around 4 GeV/c. Compared to measurements of charged particles, a similar $p_\mathrm{T}$ dependence, but smaller magnitude for $p_\mathrm{T} < $6 GeV/c, is found for prompt D$^0$ meson $v_2$ and $v_3$ coefficients. The results are consistent with the presence of collective motion of charm quarks at low $p_\mathrm{T}$ and a path length dependence of charm quark energy loss at high $p_\mathrm{T}$, thereby providing new constraints on the theoretical description of the interactions between charm quarks and the quark-gluon plasma.

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.

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 differential production cross sections of B+/- mesons are measured via the exclusive decay channels B+/- to J/psi K+/- to mu+ mu- K+/- as a function of transverse momentum in pp and PbPb collisions at a center-of-mass energy sqrt(s[NN]) = 5.02 TeV per nucleon pair with the CMS detector at the LHC. The pp (PbPb) dataset used for this analysis corresponds to an integrated luminosity of 28.0 inverse picobarns (351 inverse microbarns). The measurement is performed in the B+/- meson transverse momentum range of 7 to 50 GeV/c, in the rapidity interval abs(y) < 2.4. In this kinematic range, a strong suppression of the production cross section by about a factor of two is observed in the PbPb system in comparison to the expectation from pp reference data. These results are found to be roughly compatible with theoretical calculations incorporating beauty quark diffusion and energy loss in a quark-gluon plasma.

This paper reports the measurement of J/psi meson production in proton-proton (pp) and proton-lead (pPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV by the CMS experiment at the LHC. The data samples used in the analysis correspond to integrated luminosities of 28 inverse picobarns and 35 inverse nanobarns for pp and pPb collisions, respectively. Prompt and nonprompt J/psi mesons, the latter produced from the decay of B mesons, are measured in their dimuon decay channels. Differential cross sections are measured in the transverse momentum range of 2 < pt < 30 GeV/c, and center-of-mass rapidity ranges of abs(y[CM]) < 2.4 (pp) and -2.87 < y[CM] < 1.93 (pPb). The nuclear modification factor, R[pPb], is measured as functions of both pt and y[CM]. Small modifications of the J/psi cross sections are observed in pPb relative to pp collisions. The ratio of J/psi production cross sections in p-going and Pb-going directions, R[FB], studied as functions of pt and y[CM], shows a significant decrease for increasing transverse energy deposited at large pseudorapidities. These results, which cover a wide kinematic range, provide new insight on the role of cold nuclear matter effects on prompt and nonprompt J/psi production.

The relative modification of the prompt psi(2S) and J/psi yields from pp to PbPb collisions, at the center of mass energy of 5.02 TeV per nucleon pair, is presented. The analysis is based on pp and PbPb data samples collected by the CMS experiment at the LHC in 2015, corresponding to integrated luminosities of 28.0 inverse picobarns and 464 inverse microbarns, respectively. The double ratio of measured yields of prompt charmonia reconstructed through their decays into muon pairs, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/N[J/psi])[pp], is determined as a function of PbPb collision centrality and charmonium transverse momentum pt, in two kinematic intervals: abs(y) < 1.6 covering 6.5 < pt < 30 GeV/c and 1.6 < abs(y) < 2.4 covering 3 < pt < 30 GeV/c. The centrality-integrated double ratios are 0.36 +/- 0.08 (stat) +/-0.05 (syst) in the first interval and 0.24 +/- 0.22 (stat) +/- 0.09 (syst) in the second. The double ratio is lower than unity in all the measured bins, suggesting that the psi(2S) yield is more suppressed than the J/psi yield in the explored phase space.