The J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 33.6 fb${-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1$^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst) ]$\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

A search is described for the production of a pair of bottom-type vector-like quarks (B VLQs) with mass greater than 1000 GeV. Each B VLQ decays into a b quark and a Higgs boson, a b quark and a Z boson, or a t quark and a W boson. This analysis considers both fully hadronic final states and those containing a charged lepton pair from a Z boson decay. The products of the H $to$ bb boson decay and of the hadronic Z or W boson decays can be resolved as two distinct jets or merged into a single jet, so the final states are classified by the number of reconstructed jets. The analysis uses data corresponding to an integrated luminosity of 138 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC from 2016 to 2018. No excess over the expected background is observed. Lower limits are set on the B VLQ mass at 95% confidence level. These depend on the B VLQ branching fractions and are 1570 and 1540 GeV for 100% B $\to$ bH and 100% B $\to$ bZ, respectively. In most cases, the mass limits obtained exceed previous limits by at least 100 GeV.

A search for exotic decays of the Higgs boson (H) with a mass of 125 GeV to a pair of light pseudoscalars $\mathrm{a}_1$ is performed in final states where one pseudoscalar decays to two b quarks and the other to a pair of muons or $\tau$ leptons. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded with the CMS detector is analyzed. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set at 95% confidence level (CL) on the Higgs boson branching fraction to $\mu\mu$bb and to $\tau\tau$bb, via a pair of $\mathrm{a}_1$s. The limits depend on the pseudoscalar mass $m_{\mathrm{a}_1}$ and are observed to be in the range (0.17-3.3) $\times$ 10$^{-4}$ and (1.7-7.7) $\times$ 10$^{2}$ in the $\mu\mu$bb and $\tau\tau$bb final states, respectively. In the framework of models with two Higgs doublets and a complex scalar singlet (2HDM+S), the results of the two final states are combined to determine model-independent upper limits on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$ $\to$ $\ell\ell$bb) at 95% CL, with $\ell$ being a muon or a $\tau$ lepton. For different types of 2HDM+S, upper bounds on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$) are extracted from the combination of the two channels. In most of the Type II 2HDM+S parameter space, $\mathcal{B}($H $\to$ $\mathrm{a}_1\mathrm{a}_1$) values above 0.23 are excluded at 95% CL for $m_{\mathrm{a}_1}$ values between 15 and 60 GeV.

A search for the central exclusive production of top quark-antiquark pairs ($\mathrm{t\bar{t}}$) is performed for the first time using proton-tagged events in proton-proton collisions at the LHC at a centre-of-mass energy of 13 TeV. The data correspond to an integrated luminosity of 29.4 fb$^{-1}$. The $\mathrm{t\bar{t}}$ decay products are reconstructed using the central CMS detector, while forward protons are measured in the CMS-TOTEM precision proton spectrometer. An observed (expected) upper bound on the production cross section of 0.59 (1.14) pb is set at 95% confidence level, for collisions of protons with fractional momentum losses between 2 and 20%.

A search is presented for a third-generation leptoquark (LQ) coupled exclusively to a $\tau$ lepton and a b quark. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV recorded with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with $\tau$ leptons and a varying number of jets originating from b quarks are considered, targeting the single and pair production of LQs, as well as nonresonant $t$-channel LQ exchange. An excess is observed in the data with respect to the background expectation in the combined analysis of all search regions. For a benchmark LQ mass of 2 TeV and an LQ-b-$\tau$ coupling strength of 2.5, the excess reaches a local significance of up to 2.8 standard deviations. Upper limits at the 95% confidence level are placed on the LQ production cross section in the LQ mass range 0.5-2.3 TeV, and up to 3 TeV for $t$-channel LQ exchange. Leptoquarks are excluded below masses of 1.22-1.88 TeV for different LQ models and varying coupling strengths up to 2.5. The study of nonresonant $\tau\tau$ production through $t$-channel LQ exchange allows lower limits on the LQ mass of up to 2.3 TeV to be obtained.

The strange quark content of the proton is probed through the measurement of the production cross section for a W boson and a charm (c) quark in proton-proton collisions at a center-of-mass energy of 13 TeV. The analysis uses a data sample corresponding to a total integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. The W bosons are identified through their leptonic decays to an electron or a muon, and a neutrino. Charm jets are tagged using the presence of a muon or a secondary vertex inside the jet. The W+c production cross section and the cross section ratio $R^\pm_\text{c}$ = $\sigma$(W$^+$+$\bar{\text{c}}$) / $\sigma$(W$^-$+$\text{c}$) are measured inclusively and differentially as functions of the transverse momentum and the pseudorapidity of the lepton originating from the W boson decay. The precision of the measurements is improved with respect to previous studies, reaching 1% in $R^\pm_\text{c}$. The precision of the measurements is improved with respect to previous studies, reaching 1% in $R^\pm_\text{c}$ = 0.950 $\pm$ 0.005 (stat) $\pm$ 0.010 (syst). The measurements are compared with theoretical predictions up to next-to-next-to-leading order in perturbative quantum chromodynamics.

The observation of the production of four top quarks in proton-proton collisions is reported, based on a data sample collected by the CMS experiment at a center-of-mass energy of 13 TeV in 2016-2018 at the CERN LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with two same-sign, three, or four charged leptons (electrons and muons) and additional jets are analyzed. Compared to previous results in these channels, updated identification techniques for charged leptons and jets originating from the hadronization of b quarks, as well as a revised multivariate analysis strategy to distinguish the signal process from the main backgrounds, lead to an improved expected signal significance of 4.9 standard deviations above the background-only hypothesis. Four top quark production is observed with a significance of 5.6 standard deviations, and its cross section is measured to be 17.7 $^{+3.7}_{-3.5}$ (stat) $^{+2.3}_{-1.9}$ (syst) fb, in agreement with the available standard model predictions.

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 for resonances in events with the $\gamma$+jet final state has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC. The total data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$. Models of excited quarks and quantum black holes are considered. Using a wide-jet reconstruction for the candidate jet, the $\gamma$+jet invariant mass spectrum measured in data is examined for the presence of resonances over the standard model continuum background. The background is estimated by fitting the mass distribution with a functional form. The data exhibit no statistically significant deviations from the expected standard model background. Exclusion limits at 95% confidence level on the resonance mass and other parameters are set. Excited light-flavor quarks (excited bottom quarks) are excluded up to a mass of 6.0 (3.8) TeV. Quantum black hole production is excluded for masses up to 7.5 (5.2) TeV in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model. These lower mass bounds are the most stringent to date among those obtained in the $\gamma$+jet final state.

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.