The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

The path-length dependent parton energy loss within the dense partonic medium created in lead-lead collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV is studied by determining the azimuthal anisotropies for dijets with high transverse momentum. The data were collected by the CMS experiment in 2018 and correspond to an integrated luminosity of 1.69 nb$^{-1}$. For events containing back-to-back jets, correlations in relative azimuthal angle and pseudorapidity ($\eta$) between jets and hadrons, and between two hadrons, are constructed. The anisotropies are expressed as the Fourier expansion coefficients $v_n$, $n = $ 2-4 of these azimuthal distributions. The dijet $v_n$ values are extracted from long-range (1.5 $\lt$$\vert\Delta\eta\vert$$\lt$ 2.5) components of these correlations, which suppresses the background contributions from jet fragmentation processes. Positive dijet $v_2$ values are observed which increase from central to more peripheral events, while the $v_3$ and $v_4$ values are consistent with zero within experimental uncertainties.

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 charged particle transverse momentum spectra in the midrapidity and forward pseudorapidity ranges are presented for pPb collisions at $\sqrt{s_{\rm NN}}=5.02$~TeV. The data sample corresponding to an integrated luminosity of 26~${\rm nb}^{-1}$ was collected with the CMS detector at the LHC. The nuclear modification factor is measured at midrapidity by normalizing the measured pPb spectrum to a pp reference spectrum constructed from previous measurements. In addition, the asymmetries in the charged particle yields between equivalent positive and negative pseudorapidity ranges in the center-of-mass frame are presented as a function of transverse momentum.

The transverse momentum balance of pairs of back-to-back b quark jets in PbPb and pp collisions recorded with the CMS detector at the LHC is reported. The center-of-mass energy in both collision systems is 5.02 TeV per nucleon pair. Compared to the pp collision baseline, b quark jets have a larger imbalance in the most central PbPb collisions, as expected from the jet quenching effect. The data are also compared to the corresponding measurement with inclusive dijets. In the most central collisions, the imbalance of b quark dijets is comparable to that of inclusive dijets.

The pseudorapidity distributions of dijets as a function of their average transverse momentum ($p_\mathrm{T}^\text{ave}$) are measured in proton-lead (pPb) and proton-proton (pp) collisions. The data samples were collected by the CMS experiment at the CERN LHC, at a nucleon-nucleon center-of-mass energy of 5.02 TeV. A significant modification of the pPb spectra with respect to the pp spectra is observed in all $p_\mathrm{T}^\text{ave}$ intervals investigated. The ratios of the pPb and pp distributions are compared to next-to-leading order perturbative quantum chromodynamics calculations with unbound nucleon and nuclear parton distribution functions (PDFs). These results give the first evidence that the gluon PDF at large Bjorken $x$ in lead ions is strongly suppressed with respect to the PDF in unbound nucleons.

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 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.

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.

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.

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 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.

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.

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.

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.

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.

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.

Modifications of the properties of jets in PbPb collisions, relative to those in pp collisions, are studied at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles with the jet axis in relative pseudorapidity ($\Delta \eta$), relative azimuth ($\Delta \phi$), and relative angular distance from the jet axis $\Delta \mathrm{r} = \sqrt{{(\Delta\eta)^{2}+(\Delta\phi)^{2}}}$. This analysis uses data collected with the CMS detector at the LHC, corresponding to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Charged particle number densities, jet fragmentation functions, and jet shapes are presented as a function of PbPb collision centrality and charged-particle track transverse momentum, providing a differential description of jet modifications due to interactions with the quark-gluon plasma.

The modification of jet shapes in PbPb collisions, relative to those in pp collisions, is studied for jets associated with an isolated photon. The data were collected with the CMS detector at the LHC at a nucleon-nucleon center-of-mass energy of 5.02 TeV. Jet shapes are constructed from charged particles with track transverse momenta ($p_\mathrm{T}$) above 1 GeV/c in annuli around the axes of jets with $p^\mathrm{jet}_\mathrm{T}$ $>$ 30 GeV/c associated with an isolated photon with $p^\gamma_\mathrm{T}$ $>$ 60 GeV/c. The jet shape distributions are consistent between peripheral PbPb and pp collisions, but are modified for more central PbPb collisions. In these central PbPb events, a larger fraction of the jet momentum is observed at larger distances from the jet axis compared to pp, reflecting the interaction between the partonic medium created in heavy ion collisions and the traversing partons.

Two-particle correlations are presented for K$^0_\mathrm{S}$, $\Lambda$, and $\overline\Lambda$ strange hadrons as a function of relative momentum in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The dataset corresponds to an integrated luminosity of 0.607 nb$^{-1}$ and was collected using the CMS detector at the CERN LHC. These correlations are sensitive to quantum statistics and to final-state interactions between the particles. The source size extracted from the K$^0_\mathrm{S}$K$^0_\mathrm{S}$ correlations is found to decrease from 4 to 1 fm in going from central to peripheral collisions. Strong interaction scattering parameters (i.e., scattering length and effective range) are determined from the $\Lambda$K$^0_\mathrm{S}$ and $\Lambda\Lambda$ (including their charge conjugates) correlations using the Lednick$\'y$-Lyuboshitz model and are compared to theoretical and other experimental results.

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 production cross sections of B$^0_\mathrm{s}$ mesons and charge conjugates are measured in proton-proton (pp) and PbPb collisions via the exclusive decay channel B$^0_\mathrm{s}$ $\to$ $\mu^+\mu^-$K$^+$K$^-$ at a centre-of-mass energy of 5.02 TeV per nucleon pair and within the rapidity range $|y|$ $\lt$ 2.4 using the CMS detector at the LHC. The pp measurement is performed as a function of transverse momentum (p$_\mathrm{T}$) of the B$^0_\mathrm{s}$ meson in the range of 7 to 50 GeV/$c$ and is compared to the predictions of perturbative QCD calculations. The B$^0_\mathrm{s}$ production yield in PbPb collisions is measured in two p$_\mathrm{T}$ intervals, 7 to 15 and 15 to 50 GeV/$c$, and compared to the yield in pp collisions in the same kinematic region. The nuclear modification factor ($R_\mathrm{AA}$) is found to be 1.5 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) for 7-15 GeV/$c$, and 0.87 $\pm$ 0.30 (stat) $\pm$ 0.17 (syst) for 15-50 GeV/$c$, respectively. Within current uncertainties, the results are consistent with models of strangeness enhancement and a suppression as observed for the B$^+$ mesons.

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.

Bose-Einstein correlations between identical particles are measured in samples of proton-proton collisions at 0.9 and 7 TeV centre-of-mass energies, recorded by the CMS experiment at the LHC. The signal is observed in the form of an enhancement of number of pairs of same-sign charged particles with small relative momentum. The dependence of this enhancement on kinematic and topological features of the event is studied.