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.

The first search for scalar leptoquarks produced in $\tau$-lepton-quark collisions is presented. It is based on a set of proton-proton collision data recorded with the CMS detector at the LHC at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The reconstructed final state consists of a jet, significant missing transverse momentum, and a $\tau$ lepton reconstructed through its hadronic or leptonic decays. Limits are set on the product of the leptoquark production cross section and branching fraction and interpreted as exclusions in the plane of the leptoquark mass and the leptoquark-$\tau$-quark coupling strength.

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 mass of the top quark is measured using a sample of $\mathrm{t\overline{t}}$ events collected by the CMS detector using proton-proton collisions at $\sqrt{s} =$ 13 TeV at the CERN LHC. Events are selected with one isolated muon or electron and at least four jets from data corresponding to an integrated luminosity of 35.9 fb$^{-1}$. For each event the mass is reconstructed from a kinematic fit of the decay products to a $\mathrm{t\overline{t}}$ hypothesis. Using the ideogram method, the top quark mass is determined simultaneously with an overall jet energy scale factor (JSF), constrained by the mass of the W boson in $\mathrm{q\overline{q}'}$ decays. The measurement is calibrated on samples simulated at next-to-leading order matched to a leading-order parton shower. The top quark mass is found to be 172.25 $\pm$ 0.08 (stat+JSF) $\pm$ 0.62 (syst) GeV. The dependence of this result on the kinematic properties of the event is studied and compared to predictions of different models of $\mathrm{t\overline{t}}$ production, and no indications of a bias in the measurements are observed.

A search for production of the supersymmetric partners of the top quark, top squarks, is presented. The search is based on proton-proton collision events containing multiple jets, no leptons, and large transverse momentum imbalance. The data were collected with the CMS detector at the CERN LHC at a center-of-mass energy of 13 TeV, and correspond to an integrated luminosity of 137 fb$^{-1}$. The targeted signal production scenarios are direct and gluino-mediated top squark production, including scenarios in which the top squark and neutralino masses are nearly degenerate. The search utilizes novel algorithms based on deep neural networks that identify hadronically decaying top quarks and W bosons, which are expected in many of the targeted signal models. No statistically significant excess of events is observed relative to the expectation from the standard model, and limits on the top squark production cross section are obtained in the context of simplified supersymmetric models for various production and decay modes. Exclusion limits as high as 1310 GeV are established at the 95% confidence level on the mass of the top squark for direct top squark production models, and as high as 2260 GeV on the mass of the gluino for gluino-mediated top squark production models. These results represent a significant improvement over the results of previous searches for supersymmetry by CMS in the same final state.

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.

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.

A search for the standard model Higgs boson decaying to a charm quark-antiquark pair, H $\to$$\mathrm{c\bar{c}}$, produced in association with a leptonically decaying V (W or Z) boson is presented. The search is performed with proton-proton collisions at $\sqrt{s}$ = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb$^{-1}$. Novel charm jet identification and analysis methods using machine learning techniques are employed. The analysis is validated by searching for Z $\to$$\mathrm{c\bar{c}}$ in VZ events, leading to its first observation at a hadron collider with a significance of 5.7 standard deviations. The observed (expected) upper limit on $\sigma$(VH)$ \mathcal{B}$(H $\to$$\mathrm{c\bar{c}}$) is 0.94 (0.50 $^{+0.22}_{-0.15}$) pb at 95% confidence level (CL), corresponding to 14 (7.6 $^{+3.4}_{-2.3}$) times the standard model prediction. For the Higgs-charm Yukawa coupling modifier, $\kappa_\mathrm{c}$, the observed (expected) 95% CL interval is 1.1 $\lt$ $\vert\kappa_\mathrm{C}\vert$ $\lt$ 5.5 ($\vert\kappa_\mathrm{c}\vert$ $\lt$ 3.4), the most stringent constraint to date.

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 observation of WW$\gamma$ production in proton-proton collisions at a center-of-mass energy of 13 TeV with an integrated luminosity of 138 fb$^{-1}$ is presented. The observed (expected) significance is 5.6 (5.1) standard deviations. Events are selected by requiring exactly two leptons (one electron and one muon) of opposite charge, moderate missing transverse momentum, and a photon. The measured fiducial cross section for WW$\gamma$ is 5.9 $\pm$ 0.8 (stat) $\pm$ 0.8 (syst) $\pm$ 0.7 (modeling) fb, in agreement with the next-to-leading order quantum chromodynamics prediction. The analysis is extended with a search for the associated production of the Higgs boson and a photon, which is generated by a coupling of the Higgs boson to light quarks. The result is used to constrain the Higgs boson couplings to light quarks.

A search for exclusive two-photon production via photon exchange in proton-proton collisions, pp $\to$ p$\gamma\gamma$p with intact protons, is presented. The data correspond to an integrated luminosity of 9.4 fb$^{-1}$ collected in 2016 using the CMS and TOTEM detectors at a center-of-mass energy of 13 TeV at the LHC. Events are selected with a diphoton invariant mass above 350 GeV and with both protons intact in the final state, to reduce backgrounds from strong interactions. The events of interest are those where the invariant mass and rapidity calculated from the momentum losses of the forward-moving protons matches the mass and rapidity of the central, two-photon system. No events are found that satisfy this condition. Interpreting this result in an effective dimension-8 extension of the standard model, the first limits are set on the two anomalous four-photon coupling parameters. If the other parameter is constrained to its standard model value, the limits at 95% CL are $\lvert\zeta_1\rvert$ $\lt$ 2.9 $\times$ 10$^{-13}$ GeV$^{-4}$ and $\lvert\zeta_2\rvert$ $\lt$ 6.0 $\times$ 10$^{-13}$ GeV$^{-4}$.

A search is reported for near-threshold structures in the J/$\psi$J/$\psi$ invariant mass spectrum produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV from data collected by the CMS experiment, corresponding to an integrated luminosity of 135 fb$^{-1}$. Three structures are found, and a model with quantum interference among these structures provides a good description of the data. A new structure is observed with a significance above 5 standard deviations at a mass of 6638 $^{+43}_{-38}$ (stat) $^{+16}_{-31}$ (syst) MeV. Another structure with even higher significance is found at a mass of 6847 $^{+44}_{-28}$ (stat) $^{+48}_{-20}$ (syst) MeV, which is consistent with the X(6900) resonance reported by the LHCb experiment and confirmed by the ATLAS experiment. Evidence for another new structure, with a local significance of 4.7 standard deviations, is found at a mass of 7134 $^{+48}_{-25}$ (stat) $^{+41}_{-15}$ (syst) MeV. Results are also reported for a model without interference, which does not fit the data as well and shows mass shifts up to 150 MeV relative to the model with interference.

A search for Kaluza-Klein excited vector boson resonances, $W_\mathrm{KK}$, decaying in cascade to three W bosons via a scalar radion $R, W_\mathrm{KK}\to WR \to WWW$, with two or three massive jets is presented. The search is performed with proton-proton collision data recorded at $\sqrt{s} =$ 13 TeV, collected by the CMS experiment at the CERN LHC, during 2016-2018, corresponding to an integrated luminosity of 138 fb$^{-1}$. Two final states are simultaneously probed, one where the two W bosons produced by the R decay are reconstructed as separate, large-radius, massive jets, and one where they are merged in a single large-radius jet. The observed data are in agreement with the standard model expectations. Limits are set on the product of the $W_\mathrm{KK}$ resonance cross section and branching fraction to three W bosons in an extended warped extra-dimensional model and are the first of their kind at the LHC.

A search for invisible decays of the Higgs boson produced via vector boson fusion (VBF) has been performed with 101 fb$^{-1}$ of proton-proton collisions delivered by the LHC at $\sqrt{s} =$ 13 TeV and collected by the CMS detector in 2017 and 2018. The sensitivity to the VBF production mechanism is enhanced by constructing two analysis categories, one based on missing transverse momentum, and a second based on the properties of jets. In addition to control regions with Z and W boson candidate events, a highly populated control region, based on the production of a photon in association with jets, is used to constrain the dominant irreducible background from the invisible decay of a Z boson produced in association with jets. The results of this search are combined with all previous measurements in the VBF topology, based on data collected in 2012 (at $\sqrt{s} =$ 8 TeV), 2015, and 2016, corresponding to integrated luminosities of 19.7, 2.3, and 36.3 fb$^{-1}$, respectively. The observed (expected) upper limit on the invisible branching fraction of the Higgs boson is found to be 0.18 (0.10) at the 95% confidence level, assuming the standard model production cross section. The results are also interpreted in the context of Higgs-portal models.

A fiducial cross section for W$\gamma$ production in proton-proton collisions is measured at a center-of-mass energy of 13 TeV in 137 fb$^{-1}$ of data collected using the CMS detector at the LHC. The W $\to$ e$\nu$ and $\mu\nu$ decay modes are used in a maximum-likelihood fit to the lepton-photon invariant mass distribution to extract the combined cross section. The measured cross section is compared with theoretical expectations at next-to-leading order in quantum chromodynamics. In addition, 95% confidence level intervals are reported for anomalous triple-gauge couplings within the framework of effective field theory.

A search for dark matter in events with a displaced nonresonant muon pair and missing transverse momentum is presented. The analysis is performed using an integrated luminosity of 138 fb$^{-1}$ of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV produced by the LHC in 2016-2018. No significant excess over the predicted backgrounds is observed. Upper limits are set on the product of the inelastic dark matter production cross section $\sigma$(pp $\to$ A' $\to$$\chi_1$$\chi_2$) and the decay branching fraction $\mathcal{B}$($\chi_2$$\to$$\chi_1 \mu^+ \mu^-$), where A' is a dark photon and $\chi_1$ and $\chi_2$ are states in the dark sector with near mass degeneracy. This is the first dedicated collider search for inelastic dark matter.

A search is presented for physics beyond the standard model, based on measurements of dijet angular distributions in proton-proton collisions at $\sqrt{s}=$ 13 TeV. The data collected with the CMS detector at the LHC correspond to an integrated luminosity of 35.9 fb$^{-1}$. The observed distributions are found to be in agreement with predictions from perturbative quantum chromodynamics that include electroweak corrections. Constraints are placed on models containing quark contact interactions, extra spatial dimensions, quantum black holes, or dark matter using the detector-level distributions. In a benchmark model where only left-handed quarks participate, contact interactions are excluded at the 95% confidence level up to a scale of 12.8 or 17.5 TeV, for destructive or constructive interference, respectively. The most stringent lower limits to date are set on the ultraviolet cutoff in the Arkani-Hamed-Dimopoulos-Dvali model of extra dimensions. In the Giudice-Rattazzi-Wells convention, the cutoff scale is excluded up to 10.1 TeV. The production of quantum black holes is excluded for masses below 5.9 and 8.2 TeV, depending on the model. For the first time, lower limits between 2.0 and 4.6 TeV are set on the mass of a dark matter mediator for (axial-)vector mediators, for the universal quark coupling $g_\mathrm{q}\geq$ 1.

Many new physics models, including versions of supersymmetry characterized by $R$-parity violation (RPV), compressed mass spectra, long decay chains, or additional hidden sectors, predict the production of events with top quarks, low missing transverse momentum, and many additional quarks or gluons. The results of a search for new physics in events with two top quarks and additional jets are reported. The search is performed using events with at least seven jets and exactly one electron or muon. No requirement on missing transverse momentum is imposed. The study is based on a sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV corresponding to 137 fb$^{-1}$ of integrated luminosity collected with the CMS detector at the LHC in 2016-2018. The data are used to determine best fit values and upper limits on the cross section for pair production of top squarks in scenarios of RPV and stealth supersymmetry. Top squark masses up to 670 (870) GeV are excluded at 95% confidence level for the RPV (stealth) scenario, and the maximum observed local signal significance is 2.8 standard deviations for the RPV scenario with top squark mass of 400 GeV.

A search is presented for lepton-flavor violating decays of the Higgs boson to $\mu\tau$ and e$\tau$. The data set corresponds to an integrated luminosity of 137 fb$^{-1}$ collected at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV. No significant excess has been found, and the results are interpreted in terms of upper limits on lepton-flavor violating branching fractions of the Higgs boson. The observed (expected) upper limits on the branching fractions are, respectively, $\mathcal{B}($H $\to\mu\tau)$$\lt$ 0.15 (0.15)% and $\mathcal{B}($H$\to$e$\tau)$ $\lt$ 0.22 (0.16)% at 95% confidence level.

A search for the production of a top quark in association with a photon and additional jets via flavor changing neutral current interactions is presented. The analysis uses proton-proton collision data recorded by the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The search is performed by looking for processes where a single top quark is produced in association with a photon, or a pair of top quarks where one of the top quarks decays into a photon and an up or charm quark. Events with an electron or a muon, a photon, one or more jets, and missing transverse momentum are selected. Multivariate analysis techniques are used to discriminate signal and standard model background processes. No significant deviation is observed over the predicted background. Observed (expected) upper limits are set on the branching fractions of top quark decays: $\mathcal{B}$(t $\to$ u$\gamma$) $\lt$ 0.95 $\times$ 10$^{-5}$ (1.20 $\times$ 10$^{-5}$) and $\mathcal{B}$(t $\to$ c$\gamma$) $\lt$ 1.51 $\times$ 10$^{-5}$ (1.54 $\times$ 10$^{-5}$) at 95% confidence level, assuming a single nonzero coupling at a time. The obtained limit for $\mathcal{B}$(t $\to$ u$\gamma$) is similar to the current best limit, while the limit for $\mathcal{B}$(t $\to$ c$\gamma$) is significantly tighter than previous results.

A search for the lepton-flavor violating decay of the Higgs boson and potential additional Higgs bosons with a mass in the range 110-160 GeV to an e$^{\pm}\mu^{\mp}$ pair is presented. The search is performed with a proton-proton collision dataset at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess is observed for the Higgs boson. The observed (expected) upper limit on the e$^{\pm}\mu^{\mp}$ branching fraction for it is determined to be 4.4 (4.7) $\times$ 10$^{-5}$ at 95% confidence level, the most stringent limit set thus far from direct searches. The largest excess of events over the expected background in the full mass range of the search is observed at an e$^{\pm}\mu^{\mp}$ invariant mass of approximately 146 GeV with a local (global) significance of 3.8 (2.8) standard deviations.

Protons consist of three valence quarks, two up-quarks and one down-quark, held together by gluons and a sea of quark-antiquark pairs. Collectively, quarks and gluons are referred to as partons. In a proton-proton collision, typically only one parton of each proton undergoes a hard scattering - referred to as single-parton scattering - leaving the remainder of each proton only slightly disturbed. Here, we report the study of double- and triple-parton scatterings through the simultaneous production of three J/$\psi$ mesons, which consist of a charm quark-antiquark pair, in proton-proton collisions recorded with the CMS experiment at the Large Hadron Collider. We observed this process - reconstructed through the decays of J/$\psi$ mesons into pairs of oppositely charged muons - with a statistical significance above five standard deviations. We measured the inclusive fiducial cross section to be 272 $^{+141}_{-104}$ (stat) $\pm$ 17 (syst) fb, and compared it to theoretical expectations for triple-J/$\psi$ meson production in single-, double- and triple-parton scattering scenarios. Assuming factorization of multiple hard-scattering probabilities in terms of single-parton scattering cross sections, double- and triple-parton scattering are the dominant contributions for the measured process.

A search for resonances decaying into a W boson and a radion, where the radion decays into two W bosons, is presented. The data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$ recorded in proton-proton collisions with the CMS detector at $\sqrt{s} =$ 13 TeV. One isolated charged lepton is required, together with missing transverse momentum and one or two massive large-radius jets, containing the decay products of either two or one W bosons, respectively. No excess over the background estimation is observed. The results are combined with those from a complementary channel with an all-hadronic final state, described in an accompanying paper. Limits are set on parameters of an extended warped extra-dimensional model. These searches are the first of their kind at the LHC.

Multiparticle azimuthal correlations of prompt D$^0$ mesons are measured in PbPb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. For the first time, a four-particle cumulant method is used to extract the second Fourier coefficient of the azimuthal distribution ($v_2$) of D$^0$ mesons as a function of event centrality and the D$^0$ transverse momentum. The ratios of the four-particle $v_2$ values to previously measured two-particle cumulant results provide direct experimental access to event-by-event fluctuations of charm quark azimuthal anisotropies. These ratios are also found to be comparable to those of inclusive charged particles in the event. However, hints of deviations are seen in the most central and peripheral collisions. To investigate the origin of flow fluctuations in the charm sector, these measurements are compared with models implementing fluctuations of charm quark energy loss via collisional or radiative processes in the quark-gluon plasma. These models cannot quantitatively describe the data over the full transverse momentum and centrality ranges, although the calculations with collisional energy loss provide a better description of the data.

The leptonic and inclusive hadronic decay branching fractions of the W boson are measured using proton-proton collision data collected at $\sqrt{s} =$ 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Events characterized by the production of one or two W bosons are selected and categorized based on the multiplicity and flavor of reconstructed leptons, the number of jets, and the number of jets identified as originating from the hadronization of b quarks. A binned maximum likelihood estimate of the W boson branching fractions is performed simultaneously in each event category. The measured branching fractions of the W boson decaying into electron, muon, and tau lepton final states are (10.83 $\pm$ 0.10)%, (10.94 $\pm$ 0.08)%, and (10.77 $\pm$ 0.21)%, respectively, consistent with lepton flavor universality for the weak interaction. The average leptonic and inclusive hadronic decay branching fractions are estimated to be (10.89 $\pm$ 0.08)% and (67.32 $\pm$ 0.23)%, respectively. Based on the hadronic branching fraction, three standard model quantities are subsequently derived: the sum of squared elements in the first two rows of the Cabibbo-Kobayashi-Maskawa (CKM) matrix $\sum_{ij}\vert V_{ij}\vert^{2}$ = 1.984 $\pm$ 0.021, the CKM element $\vert V_\mathrm{cs}\vert$ = 0.967 $\pm$ 0.011, and the strong coupling constant at the W boson mass scale, $\alpha_\mathrm{S}(m^2_\mathrm{W})$ = 0.095 $\pm$ 0.033.