A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018. The analyzed data correspond to an integrated luminosity of up to 137 fb$^{-1}$. The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage.

A search for the production of pairs of heavy Majorana neutrinos (N$_\ell$) from the decays of Z' bosons is performed using the CMS detector at the LHC. The data were collected in proton-proton collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, with an integrated luminosity of 138 fb$^{-1}$. The signature for the search is an excess in the invariant mass distribution of the final-state objects, two same-flavor leptons (e or $\mu$) and at least two jets. No significant excess of events beyond the expected background is observed. Upper limits at 95% confidence level are set on the product of the Z' production cross section and its branching fraction to a pair of N$_\ell$, as functions of N$_\ell$ and Z' boson masses ($m_{\mathrm{N}_\ell}$ and $m_\mathrm{Z'}$, respectively) for $m_\mathrm{Z'}$ from 0.4 to 4.6 TeV and $m_{\mathrm{N}_\ell}$ from 0.1 TeV to $m_\mathrm{Z'}$/2. In the theoretical framework of a left-right symmetric model, exclusion bounds in the $m_{\mathrm{N}_\ell}-m_\mathrm{Z'}$ plane are presented in both the electron and muon channels. The observed upper limit on $m_\mathrm{Z'}$ reaches up to 4.42 TeV. These are the most restrictive limits to date on the mass of N$_\ell$ as a function of the Z' boson mass.

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.

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.

The nuclear modification factors of J/$\psi$ and $\psi$(2S) mesons are measured in PbPb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC in 2015, corresponding to integrated luminosities of 464 $\mu$b$^{-1}$ and 28 pb$^{-1}$, respectively. The measurements are performed in the dimuon rapidity range of $|y| <$ 2.4 as a function of centrality, rapidity, and transverse momentum (p$_\mathrm{T}$) from p$_\mathrm{T}=$ 3 GeV/$c$ in the most forward region and up to 50 GeV/$c$. Both prompt and nonprompt (coming from b hadron decays) mesons are observed to be increasingly suppressed with centrality, with a magnitude similar to the one observed at $\sqrt{s_{\mathrm{NN}}}=$ 2.76 TeV for the two J/$\psi$ meson components. No dependence on rapidity is observed for either prompt or nonprompt J/$\psi$ mesons. An indication of a lower prompt J/$\psi$ meson suppression at p$_\mathrm{T} >$ 25 GeV/$c$ is seen with respect to that observed at intermediate p$_\mathrm{T}$. The prompt $\psi$(2S) meson yield is found to be more suppressed than that of the prompt J/$\psi$ mesons in the entire p$_\mathrm{T}$ range.

A search is presented for the decay of the 125 GeV Higgs boson (H) to a pair of new light pseudoscalar bosons (a), followed by the prompt decay of each a boson to a bottom quark-antiquark pair, H $\to$ aa $\to$$\mathrm{b\bar{b}b\bar{b}}$. The analysis is performed using a data sample of proton-proton collisions collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. To reduce the background from standard model processes, the search requires the Higgs boson to be produced in association with a leptonically decaying W or Z boson. The analysis probes the production of new light bosons in a 15 $\lt$$m_\mathrm{a}$$\lt$ 60 GeV mass range. Assuming the standard model predictions for the Higgs boson production cross sections for pp $\to$ WH and ZH, model independent upper limits at 95% confidence level are derived for the branching fraction $\mathcal{B}$(H $\to$ aa $\to$ $\mathrm{b\bar{b}b\bar{b}}$). The combined WH and ZH observed upper limit on the branching fraction ranges from 1.10 for $m_\mathrm{a} =$ 20 GeV to 0.36 for $m_\mathrm{a} =$ 60 GeV, complementing other measurements in the $\mu\mu\tau\tau$, $\tau\tau\tau\tau$ and bb$\ell\ell$ ($\ell=$ $\mu$,$\tau$) channels.

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 is presented for charged, long-lived supersymmetric particles in final states with one or more disappearing tracks. The search is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the CERN LHC between 2016 and 2018, corresponding to an integrated luminosity of 137 fb$^{-1}$. The search is performed over final states characterized by varying numbers of jets, b-tagged jets, electrons, and muons. The length of signal-candidate tracks in the plane perpendicular to the beam axis is used to characterize the lifetimes of wino- and higgsino-like charginos produced in the context of the minimal supersymmetric standard model. The d$E$/d$x$ energy loss of signal-candidate tracks is used to increase the sensitivity to charginos with a large mass and thus a small Lorentz boost. The observed results are found to be statistically consistent with the background-only hypothesis. Limits on the pair production cross section of gluinos and squarks are presented in the framework of simplified models of supersymmetric particle production and decay, and for electroweakino production based on models of wino and higgsino dark matter. The limits presented are the most stringent to date for scenarios with light third-generation squarks and a wino- or higgsino-like dark matter candidate capable of explaining the known dark matter relic density.

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.

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.