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

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 $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

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.

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.