Measurements of the suppression of the yield per nucleon of J/Psi and Psi' production for 800 GeV/c protons incident on heavy relative to light nuclear targets have been made with very broad coverage in xF and pT. The observed suppression is smallest at xF values of 0.25 and below and increases at larger values of xF. It is also strongest at small pT. Substantial differences between the Psi' and J/Psi are observed for the first time in p-A collisions. The suppression for the Psi' is stronger than that for the J/Psi for xF near zero, but becomes comparable to that for the J/Psi for xF > 0.6.

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 differential cross sections of $\Lambda_\mathrm{c}^+$ baryon production are measured via the exclusive decay channel $\Lambda_\mathrm{c}^+ \to $pK$^-\pi^+$, as a function of transverse momentum ($p_\mathrm{T}$) in proton-proton (pp) and lead-lead (PbPb) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV with the CMS detector at the LHC. The measurement is performed within the $\Lambda_\mathrm{c}^+$ rapidity interval $|y|<$1.0 in the $p_\mathrm{T}$ range of 5-20 GeV/$c$ in pp and 10-20 GeV/$c$ in PbPb collisions. The observed yields of $\Lambda_\mathrm{c}^+$ for $p_\mathrm{T}$ of 10-20 GeV/$c$ suggest a possible suppression in central PbPb collisions compared to pp collisions. The $\Lambda_\mathrm{c}^+/$D$^0$ production ratio in pp collisions is compared to theoretical models. In PbPb collisions, this ratio is consistent with the result from pp collisions in their common $p_\mathrm{T}$ range.

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 momentum-weighted sum of the electric charges of particles inside a jet, known as jet charge, is sensitive to the electric charge of the particle initiating the parton shower. This paper presents jet charge distributions in $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV lead-lead (PbPb) and proton-proton (pp) collisions recorded with the CMS detector at the LHC. These data correspond to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Leveraging the sensitivity of the jet charge to fundamental differences in the electric charges of quarks and gluons, the jet charge distributions from simulated events are used as templates to extract the quark- and gluon-like jet fractions from data. The modification of these jet fractions is examined by comparing pp and PbPb data as a function of the overlap of the colliding Pb nuclei (centrality). This measurement tests the color charge dependence of jet energy loss due to interactions with the quark-gluon plasma. No significant modification between different centrality classes and with respect to pp results is observed in the extracted fractions of quark- and gluon-like jet fractions.

The first measurements of the Fourier coefficients ($V_{n\Delta}$) of the azimuthal distributions of charged hadrons emitted from photon-proton ($\gamma$p) interactions at the LHC are presented. The data are extracted from 68.8 nb$^{-1}$ of ultra-peripheral proton-lead (pPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV using the CMS detector. The high energy lead ions produce a flux of photons that can interact with the oncoming proton. This $\gamma$p system provides a set of unique initial conditions with multiplicity lower than in photon-lead collisions but comparable to recent electron-positron and electron-proton data. The $V_{n\Delta}$ coefficients are presented in ranges of event multiplicity and transverse momentum ($p_\mathrm{T}$) and are compared to corresponding hadronic minimum bias pPb results. For a given multiplicity range, the mean $p_\mathrm{T}$ of charged particles is smaller in $\gamma$p than in pPb collisions. For both the $\gamma$p and pPb samples, $V_{1\Delta}$ is negative, $V_{2\Delta}$ is positive, and $V_{3\Delta}$ consistent with 0. For each multiplicity and $p_\mathrm{T}$ range, $V_{2\Delta}$ is larger for $\gamma$p events. The $\gamma$p data are consistent with model predictions that have no collective effects.

The transverse energy ($E_\mathrm{T}^{\gamma}$) spectra of photons isolated from other particles are measured using proton-proton (pp) and lead-lead (PbPb) collisions at the LHC at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV with integrated luminosities of 27.4 pb$^{-1}$and 404 $\mu$b$^{-1}$ for pp and PbPb data, respectively. The results are presented for photons with 25 $<$ $E_\mathrm{T}^{\gamma}$ $<$ 200 GeV in the pseudorapidity range $|\eta|$ $<$ 1.44, and for different centrality intervals for PbPb collisions. Photon production in PbPb collisions is consistent with that in pp collisions scaled by the number of binary nucleon-nucleon collisions, demonstrating that photons do not interact with the quark-gluon plasma. Therefore, isolated photons can provide information about the initial energy of the associated parton in photon+jet measurements. The results are compared with predictions from the next-to-leading-order JETPHOX generator for different parton distribution functions (PDFs) and nuclear PDFs (nPDFs). The comparisons can help to constrain the nPDFs global fits.

The first evidence for X(3872) production in relativistic heavy ion collisions is reported. The X(3872) production is studied in lead-lead (PbPb) collisions at a center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV per nucleon pair, using the decay chain X(3872) $\to$ J$/\psi\, \pi^+\pi^- \to$ $\mu^+\mu^-\pi^+\pi^-$. The data were recorded with the CMS detector in 2018 and correspond to an integrated luminosity of 1.7 nb$^{-1}$. The measurement is performed in the rapidity and transverse momentum ranges $|y|$ $\lt$ 1.6 and 15 $\lt$ $p_\mathrm{T}$ $\lt$ 50 GeV$/c$. The significance of the inclusive X(3872) signal is 4.2 standard deviations. The prompt X(3872) to $\psi$(2S) yield ratio is found to be $\rho^\mathrm{PbPb} = $ 1.08 $\pm$ 0.49 (stat) $\pm$ 0.52 (syst), to be compared with typical values of 0.1 for pp collisions. This result provides a unique experimental input to theoretical models of the X(3872) production mechanism, and of the nature of this exotic state.

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.

The second-order Fourier coefficients ($v_2$) characterizing the azimuthal distribution of $\Upsilon$(1S) and $\Upsilon$(2S) mesons arising from PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV are studied. The $\Upsilon$ mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The data set corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The scalar product method is used to extract the $v_2$ coefficients of the azimuthal distribution. Results are reported for the rapidity range $|y|$$\lt$ 2.4, with the transverse momentum 0 $\lt$$p_\mathrm{T}$$\lt$ 50 GeV/$c$, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/$\psi$ mesons, the measured $v_2$ values for the $\Upsilon$ mesons are found to be consistent with zero.

Jet production in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV is studied with the CMS detector at the LHC, using PbPb and pp data samples corresponding to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$, respectively. Jets with different areas are reconstructed using the anti-$k_\mathrm{T}$ algorithm by varying the distance parameter $R$. The measurements are performed using jets with transverse momenta ($p_\mathrm{T}$) greater than 200 GeV and in a pseudorapidity range of $|\eta|$$\lt$ 2. To reveal the medium modification of the jet spectra in PbPb collisions, the properly normalized ratio of spectra from PbPb and pp data is used to extract jet nuclear modification factors as functions of the PbPb collision centrality, $p_\mathrm{T}$ and, for the first time, as a function of $R$ up to 1.0. For the most central collisions, a strong suppression is observed for high-$p_\mathrm{T}$ jets reconstructed with all distance parameters, implying that a significant amount of jet energy is scattered to large angles. The dependence of jet suppression on $R$ is expected to be sensitive to both the jet energy loss mechanism and the medium response, and so the data are compared to several modern event generators and analytic calculations. The models considered do not fully reproduce the data.

We report on the production characteristics and total cross section for 9 beauty hadron pairs produced by a 600 GeV/ c π − beam, the first such information in this energy region. The events were detected in the hybrid emulsion spectrometer of Fermilab Experiment E653. The measured pair cross section for all χ F , assuming linear A dependence, is 33±11 (stat.)±6(syst.) nb/nucleon. Fits of the inclusive single-hadron production distribution to the forms d σ d χ F ∝ (1−|χ F −χ 0 |) n and d σ d p T 2 ∝ exp (−bp t 2 ) give n=5.0 −2.1−1.7 +2.7+1.7 , χ 0 =0.06 −0.07−0.03 +0.06+0.02 , and b=0.13 −0.04−0.02 +0.05+0.02 ( GeV /c −2 . .The pairs tend to be produced back-to-back.

We present total and differential cross sections for charm mesons produced in 600 GeV/ c π - emulsion interactions. Fits to d 2 σ / dx F dp T 2 ∞ (1−| x F |) n exp (- bp T 2 ) for 676 electronically reconstructed D mesons with x F >0 give n =4.25±0.24 ( stat .)±0.23 ( syst .) and b =0.76±0.03±0.03 ( GeV / c ) -2 . The total inclusive D + and D 0 cross sections are σ ( π - N → D ± ; x F >0) = 8.66±0.46±1.96 μb nucleon and σ(π - N→D 0 D 0 ; x F >0)=22.05±1.37±4.82μb nucleonk, where a linear dependence on the mean atomic weight of the target is assumed. These results are compared to next-to-leading order QCD predictions.

We report results on D 0 and D + production in proton-emulsion interactions at s =38.7 GeV. A fit to the form (1−| x F |) n exp (−bp 2 T ) yields n=6.9 +1.9 −1.8 and b=0.84 +0.10 −0.08 (GeV/ c ) −2 . The total inclusive cross section, is assuming linear A dependence, is measured to be 38±3(stat.) ±13 (sys.) μ b for the D 0 and 38±9±14 μ b for the D + . A comparison of these results with previous measurements indicates that nuclear effects do not strongly influence charm production. The predictions of QCD are in good agreement with our data.

We present results on charm pair correlations measured in proton-emulsion interactions at s =38.7 GeV. The predictions of leading order QCD for the distributions in invariant mass, rapidity gap, x F , and polar angle in the charm pair CMS are qualitatively consistent with our measurements. The mean p T of the pairs is equal within errors to that measured in dilepton production at the same energy and mass range.

Ultrarelativistic heavy ion collisions recreate in the laboratory the thermodynamical conditions prevailing in the early universe up to 10$^{-6}$ seconds, thereby allowing the study of the quark-gluon plasma (QGP), a state of quantum chromodynamics (QCD) matter with deconfined partons. The top quark, the heaviest elementary particle known, is accessible in nucleus-nucleus collisions at the CERN LHC, and constitutes a novel probe of the QGP. Here, we report the first-ever evidence for the production of top quarks in nucleus-nucleus collisions, using lead-lead collision data at a nucleon-nucleon centre-of-mass energy of 5.02 TeV recorded by the CMS experiment. Two methods are used to measure the cross section for top quark pair production ($\sigma_\mathrm{t\bar{t}}$) via the decay into charged leptons (electrons or muons) and bottom quarks. One method relies on the leptonic information alone, and the second one exploits, in addition, the presence of bottom quarks. The measured cross sections, $\sigma_\mathrm{t\bar{t}} = $ 2.54 $^{+0.84}_{-0.74}$ and 2.03 $^{+0.71}_{-0.64}$ $\mu$b, respectively, are compatible with expectations from scaled proton-proton data and QCD predictions.

A search for charm production in the coherent diffractive dissociation reaction pSi→XSi was carried out for the modes D 0 → K − π + , D 0 → K − π + π + π − , and D + → K − π + π + . No charm signals were observed, and the 90% confidence level upper limit for coherent charm pair production was determined to be 26 μ b per silicon nucleus. The results are interpreted as an upper limit of 0.2% on the amount of intrinsic charm in the proton.

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.

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$$\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$$\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

The first measurement of the dependence of $\gamma\gamma$$\to$$\mu^{+}\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\lt$$m_{\mu\mu}$$\lt$ 60 GeV are extracted for events including 0, 1, or at least 2 neutrons detected in the forward pseudorapidity range $|\eta|$$\gt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.

The spectra measurements of charged hadrons produced in the shower of a parton originating in the same hard scattering with a leptonically decaying Z boson, are reported in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. Both PbPb and pp data sets are recorded by the CMS experiment at the LHC, and correspond to an integrated luminosity of 1.7 nb$^{-1}$ and 320 pb$^{-1}$, respectively. Hadronic collision data with one reconstructed Z boson candidate with the transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV/$c$ are analyzed. The Z boson constrains the initial energy and direction of the associated parton. In heavy ion events, azimuthal angular distributions of charged hadrons with respect to the direction of a Z boson are sensitive to modifications of the in-medium parton shower and medium response. Compared to reference data from pp interactions, the results for central PbPb collisions indicate a modification of the angular correlations. The measurements of the fragmentation functions and $p_\mathrm{T}$ spectra of charged particles in Z boson events, which are sensitive to medium modifications of the parton shower longitudinal structure, are also reported. Significant modifications in central PbPb events compared to pp reference data are also found for these observables.

Jets containing a prompt J$/\psi$ meson are studied in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV, using the CMS detector at the LHC. Jets are selected to be in the transverse momentum range of 30 $\lt$$p_\mathrm{T}$$\lt$ 40 GeV. The J$/\psi$ yield in these jets is evaluated as a function of the jet fragmentation variable $z$, the ratio of the J$/\psi$$p_\mathrm{T} $ to the jet $p_\mathrm{T}$. The nuclear modification factor, $R_\mathrm{AA}$, is then derived by comparing the yield in lead-lead collisions to the corresponding expectation based on proton-proton data, at the same nucleon-nucleon center-of-mass energy. The suppression of the J$/\psi$ yield shows a dependence on $z$, indicating that the interaction of the J$/\psi$ with the quark-gluon plasma formed in heavy ion collisions depends on the fragmentation that gives rise to the J$/\psi$ meson.

The CMS experiment at the LHC has measured the differential cross sections of Z bosons decaying to pairs of leptons, as functions of transverse momentum and rapidity, in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The measured Z boson elliptic azimuthal anisotropy coefficient is compatible with zero, showing that Z bosons do not experience significant final-state interactions in the medium produced in the collision. Yields of Z bosons are compared to Glauber model predictions and are found to deviate from these expectations in peripheral collisions, indicating the presence of initial collision geometry and centrality selection effects. The precision of the measurement allows, for the first time, for a data-driven determination of the nucleon-nucleon integrated luminosity as a function of lead-lead centrality, thereby eliminating the need for its estimation based on a Glauber model.

Modifications to the distribution of charged particles with respect to high transverse momentum ($p_\mathrm{T}$) jets passing through a quark-gluon plasma are explored using the CMS detector. Back-to-back dijets are analyzed in lead-lead and proton-proton collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles in bins of relative pseudorapidity and angular distance from the leading and subleading jet axes. In comparing the lead-lead and proton-proton collision results, modifications to the charged-particle relative distance distribution and to the momentum distributions around the jet axis are found to depend on the dijet momentum balance $x_j$, which is the ratio between the subleading and leading jet $p_\mathrm{T}$. For events with $x_j$$\approx$ 1, these modifications are observed for both the leading and subleading jets. However, while subleading jets show significant modifications for events with a larger dijet momentum imbalance, much smaller modifications are found for the leading jets in these events.

For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV, almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of two below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC, QGSJET II, and HIJING predictions are all at least a factor of five lower than the data. The latter effect might be explained by a significant contribution of ultra-peripheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.

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.

Data from heavy ion collisions suggest that the evolution of a parton shower is modified by interactions with the color charges in the dense partonic medium created in these collisions, but it is not known where in the shower evolution the modifications occur. The momentum ratio of the two leading partons, resolved as subjets, provides information about the parton shower evolution. This substructure observable, known as the splitting function, reflects the process of a parton splitting into two other partons and has been measured for jets with transverse momentum between 140 and 500 GeV, in pp and PbPb collisions at a center-of-mass energy of 5.02 GeV per nucleon pair. In central PbPb collisions, the splitting function indicates a more unbalanced momentum ratio, compared to peripheral PbPb and pp collisions. The measurements are compared to various predictions from event generators and analytical calculations.

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.

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.

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.

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.

A search for partonic collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged constituents using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity $\lvert\eta\rvert$$\lt$ 1.6. Two-particle correlations among the charged constituents within the jets are studied as functions of the particles' azimuthal angle and pseudorapidity separations ($\Delta\phi^*$ and $\Delta\eta^*$) in a jet coordinate basis, where constituents' $\eta^*$, $\phi^*$ are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to $N_\text{ch}^\mathrm{j}$$\approx$ 100. Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for $\lvert\Delta\eta^*\rvert$$\gt$ 2. For low-multiplicity jets, the long-range elliptic anisotropic harmonic, $v^*_2$, is observed to decrease with $N_\text{ch}^\mathrm{j}$. This trend is well described by Monte Carlo event generators. However, a rising trend for $v^*_2$ emerges at $N_\text{ch}^\mathrm{j}$$\gtrsim$ 80, hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of parton fragmentation processes in the vacuum.

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.

Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark-gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we extracted the speed of sound in this medium created using lead-lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 $\pm$ 0.002 (stat) $\pm$ 0.016 (syst) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 $\pm$ 8 (syst) MeV. The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

Quasireal photons exchanged in relativistic heavy ion interactions are powerful probes of the gluonic structure of nuclei. The coherent J/$\psi$ photoproduction cross section in ultraperipheral lead-lead collisions is measured as a function of photon-nucleus center-of-mass energies per nucleon (W$^\text{Pb}_{\gamma\text{N}}$), over a wide range of 40 $\lt$ W$^\text{Pb}_{\gamma\text{N}}$$\lt$ 400 GeV. Results are obtained using data at the nucleon-nucleon center-of-mass energy of 5.02 TeV collected by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 1.52 nb$^{-1}$. The cross section is observed to rise rapidly at low W$^\text{Pb}_{\gamma\text{N}}$, and plateau above W$^\text{Pb}_{\gamma\text{N}}$$\approx$ 40 GeV, up to 400 GeV, a new regime of small Bjorken-$x$ ($\approx$ 6 $\times$ 10$^{-5}$) gluons being probed in a heavy nucleus. The observed energy dependence is not predicted by current quantum chromodynamic models.

The production of prompt $\Lambda^+_\mathrm{c}$ baryons is measured via the exclusive decay channel $\Lambda^+_\mathrm{c}$$\to$ pK$^-\pi^+$ at a center-of-mass energy per nucleon pair of 5.02 TeV, using proton-proton (pp) and lead-lead (PbPb) collision data collected by the CMS experiment at the CERN LHC. The pp and PbPb data were obtained in 2017 and 2018 with integrated luminosities of 252 and 0.607 nb$^{-1}$, respectively. The measurements are performed within the $\Lambda^+_\mathrm{c}$ rapidity interval $\vert y \vert$$\lt$ 1 with transverse momentum ($p_\mathrm{T}$) ranges of 3-30 and 6-40 GeV/$c$ for pp and PbPb collisions, respectively. Compared to the yields in pp collisions scaled by the expected number of nucleon-nucleon interactions, the observed yields of $\Lambda^+_\mathrm{c}$ with $p_\mathrm{T}$$\gt$ 10 GeV/$c$ are strongly suppressed in PbPb collisions. The level of suppression depends significantly on the collision centrality. The $\Lambda^+_\mathrm{c}$ / D$^0$ production ratio is similar in PbPb and pp collisions at $p_\mathrm{T}$$\gt$ 10 GeV/$c$, suggesting that the coalescence process does not play a dominant role in prompt $\Lambda^+_\mathrm{c}$ baryon production at higher $p_\mathrm{T}$.

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.

Evidence for the light-by-light scattering process, $\gamma\gamma$ $\to$ $\gamma\gamma$, in ultraperipheral PbPb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV is reported. The analysis is conducted using a data sample corresponding to an integrated luminosity of 390 $\mu$b$^{-1}$ recorded by the CMS experiment at the LHC. Light-by-light scattering processes are selected in events with two photons exclusively produced, each with transverse energy E$_\mathrm{T}^{\gamma}$ $>$ 2 GeV, pseudorapidity $|\eta^{\gamma}|$ $\lt$ 2.4, diphoton invariant mass $m^{\gamma\gamma}$ $\gt$ 5 GeV, diphoton transverse momentum $p_\mathrm{T}^{\gamma\gamma}$ $\lt$ 1 GeV, and diphoton acoplanarity below 0.01. After all selection criteria are applied, 14 events are observed, compared to expectations of 9.0 $\pm$ 0.9 (theo) events for the signal and 4.0 $\pm$ 1.2 (stat) for the background processes. The excess observed in data relative to the background-only expectation corresponds to a significance of 3.7 standard deviations, and has properties consistent with those expected for the light-by-light scattering signal. The measured fiducial light-by-light scattering cross section, $\sigma_\mathrm{fid} (\gamma\gamma$ $\to$ $\gamma\gamma) =$ 120 $\pm$ 46 (stat) $\pm$ 28 (syst) $\pm$ 12 (theo) nb, is consistent with the standard model prediction. The $m^{\gamma\gamma}$ distribution is used to set new exclusion limits on the production of pseudoscalar axion-like particles, via the $\gamma\gamma$ $\to$ a $\to$ $\gamma\gamma$ process, in the mass range $m_{\mathrm{a}} =$ 5-90 GeV.

The first study of the shapes of jets arising from bottom (b) quarks in heavy ion collisions is presented. Jet shapes are studied using charged hadron constituents as a function of their radial distance from the jet axis. Lead-lead (PbPb) collision data at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV were recorded by the CMS detector at the LHC, with an integrated luminosity of 1.69 nb$^{-1}$. Compared to proton-proton collisions, a redistribution of the energy in b jets to larger distances from the jet axis is observed in PbPb collisions. This medium-induced redistribution is found to be substantially larger for b jets than for inclusive jets.

We present an observation of photon-photon production of $\tau$ lepton pairs in ultraperipheral lead-lead collisions. The measurement is based on a data sample with an integrated luminosity of 404 $\mu$b$^{-1}$ collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The $\gamma\gamma$$\to$$\tau^+\tau^-$ process is observed for $\tau\tau$ events with a muon and three charged hadrons in the final state. The measured fiducial cross section is $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$ = 4.8 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) $\mu$b, in agreement with leading-order QED predictions. Using $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$, we estimate a model-dependent value of the anomalous magnetic moment of the $\tau$ lepton of $a_\tau$ = 0.001 $^{+0.055}_{-0.089}$.

The production of $\Upsilon$(2S) and $\Upsilon$(3S) mesons in lead-lead (PbPb) and proton-proton (pp) collisions is studied in their dimuon decay channel using the CMS detector at the LHC. The $\Upsilon$(3S) meson is observed for the first time in PbPb collisions, with a significance above five standard deviations. The ratios of yields measured in PbPb and pp collisions are reported for both the $\Upsilon$(2S) and $\Upsilon$(3S) mesons, as functions of transverse momentum and PbPb collision centrality. These ratios, when appropriately scaled, are significantly less than unity, indicating a suppression of $\Upsilon$ yields in PbPb collisions. This suppression increases from peripheral to central PbPb collisions. Furthermore, the suppression is stronger for $\Upsilon$(3S) mesons compared to $\Upsilon$(2S) mesons, extending the pattern of sequential suppression of quarkonium states in nuclear collisions previously seen for the $\psi$/J, $\psi$(2S), $\Upsilon$(1S), and $\Upsilon$(2S) mesons.

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.

Measurements of the elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients are presented for D$^0$ mesons produced in b hadron decays (nonprompt D$^0$ mesons) in lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The results are compared with previously published charm meson anisotropies measured using prompt D$^0$ mesons. The data were collected with the CMS detector in 2018 with an integrated luminosity of 0.58 nb$^{-1}$. Azimuthal anisotropy is sensitive to the interactions of quarks with the hot and dense medium created in heavy ion collisions. Comparing results for prompt and nonprompt D$^0$ mesons can assist in understanding the mass dependence of these interactions. The nonprompt results show lower magnitudes of $v_2$ and $v_3$ and weaker dependences on the meson transverse momentum and collision centrality than those found for prompt D$^0$ mesons. The results are in agreement with theoretical predictions that include a mass dependence in the interactions of quarks with the medium.

Azimuthal correlations of charged particles in xenon-xenon collisions at a center-of-mass energy per nucleon pair of $ \sqrt{s_{_\mathrm{NN}}} =$ 5.44 TeV are studied. The data were collected by the CMS experiment at the LHC with a total integrated luminosity of 3.42 $\mu$b$^{-1}$. The collective motion of the system formed in the collision is parameterized by a Fourier expansion of the azimuthal particle density distribution. The azimuthal anisotropy coefficients $v_{2}$, $v_{3}$, and $v_{4}$ are obtained by the scalar-product, two-particle correlation, and multiparticle correlation methods. Within a hydrodynamic picture, these methods have different sensitivities to non-collective and fluctuation effects. The dependence of the Fourier coefficients on the size of the colliding system is explored by comparing the xenon-xenon results with equivalent lead-lead data. Model calculations that include initial-state fluctuation effects are also compared to the experimental results. The observed angular correlations provide new constraints on the hydrodynamic description of heavy ion collisions.

The transverse momentum spectra of D$^0$ mesons from b hadron decays are measured at midrapidity ($|y|$ $<$ $1$) in pp and PbPb collisions at a nucleon-nucleon center of mass energy of 5.02 TeV with the CMS detector at the LHC. The D$^0$ mesons from b hadron decays are distinguished from prompt D$^0$ mesons by their decay topologies. In PbPb collisions, the B $\to$ D$^0$ yield is found to be suppressed in the measured $p_\mathrm{T}$ range from 2 to 100 GeV$/c$ as compared to pp collisions. The suppression is weaker than that of prompt D$^0$ mesons and charged hadrons for $p_\mathrm{T}$ around 10 GeV$/c$. While theoretical calculations incorporating partonic energy loss in the quark-gluon plasma can successfully describe the measured B $\to$ D$^0$ suppression at higher $p_\mathrm{T}$, the data show an indication of larger suppression than the model predictions in the range of 2 $\lt$ $p_\mathrm{T}$ $\lt$ 5 GeV$/c$.

Measurements of the pseudorapidity distributions of charged hadrons produced in xenon-xenon collisions at a nucleon-nucleon centre-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.44 TeV are presented. The measurements are based on data collected by the CMS experiment at the LHC. The yield of primary charged hadrons produced in xenon-xenon collisions in the pseudorapidity range $|\eta|$ $<$ 3.2 is determined using the silicon pixel detector in the CMS tracking system. For the 5% most central collisions, the charged-hadron pseudorapidity density in the midrapidity region $|\eta|$ $<$ 0.5 is found to be 1187 $\pm$ 36 (syst), with a negligible statistical uncertainty. The rapidity distribution of charged hadrons is also presented in the range $|y|$ $<$ 3.2 and is found to be independent of rapidity around $y =$ 0. Existing Monte-Carlo event generators are unable to simultaneously describe both results. Comparisons of charged-hadron multiplicities between xenon-xenon and lead-lead collisions at similar collision energies show that particle production at midrapidity is strongly dependent on the collision geometry in addition to the system size and collision energy.

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.