A search for the pair production of the lightest supersymmetric partner of the top quark, the top squark ($\tilde{\mathrm{t}}_1$), is presented. The search targets the four-body decay of the $\tilde{\mathrm{t}}_1$, which is preferred when the mass difference between the top squark and the lightest supersymmetric particle is smaller than the mass of the W boson. This decay mode consists of a bottom quark, two other fermions, and the lightest neutralino ($\tilde{\chi}^0_1$), which is assumed to be the lightest supersymmetric particle. The data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 TeV collected by the CMS experiment at the CERN LHC. Events are selected using the presence of a high-momentum jet, an electron or muon with low transverse momentum, and a significant missing transverse momentum. The signal is selected based on a multivariate approach that is optimized for the difference between $m(\tilde{\mathrm{t}}_1)$ and $m(\tilde{\chi}^0_1)$. The contribution from leading background processes is estimated from data. No significant excess is observed above the expectation from standard model processes. The results of this search exclude top squarks at 95% confidence level for masses up to 480 and 700 GeV for $m(\tilde{\mathrm{t}}_1) - m(\tilde{\chi}^0_1$) = 10 and 80 GeV, respectively.

Measurements of light-by-light scattering (LbL, $γγ$$\to$$γγ$) and the Breit-Wheeler process (BW, $γγ$$\to$$\mathrm{e^+e^-}$) are reported in ultraperipheral PbPb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV. The data sample, corresponding to an integrated luminosity of 1.7 nb$^{-1}$, was collected by the CMS experiment at the CERN LHC in 2018. Events with an exclusively produced $γγ$ or $\mathrm{e^+e^-}$ pair with invariant masses $m^{γγ,\mathrm{ee}}$$>$ 5 GeV, along with other fiducial criteria, are selected. The measured BW fiducial production cross section, $σ_\text{fid}$($γγ$$\to$$\mathrm{e^+e^-}$) = 263.5 $\pm$ 1.8 (stat) $\pm$ 17.8 (syst) $μ$b, as well as the differential distributions for various kinematic observables, are in agreement with leading-order quantum electrodynamics predictions complemented with final-state photon radiation. The measured differential BW cross sections allow discriminating between different theoretical descriptions of the photon flux of the lead ion. In the LbL final state, 26 exclusive diphoton candidate events are observed compared with 12.0 $\pm$ 2.9 expected for the background. Combined with previous results, the observed significance of the LbL signal with respect to the background-only hypothesis is above five standard deviations. The measured fiducial LbL scattering cross section, $σ_\text{fid} (γγ$$\to$$γγ)$ = 107 $\pm$ 24 (stat) $\pm$ 13 (syst) nb, is in agreement with next-to-leading-order predictions. Limits on the production of axion-like particles coupled to photons are set over the mass range 5-100 GeV, including the most stringent limits to date in the range of 5-10 GeV.

The production cross sections of B$^0_\mathrm{s}$ and B$^+$ mesons are reported in proton-proton (pp) collisions recorded by the CMS experiment at the CERN LHC with a center-of-mass energy of 5.02 TeV. The data sample corresponds to an integrated luminosity of 302 pb$^{-1}$. The cross sections are based on measurements of the B$^0_\mathrm{s}$$\to$ J/$ψ(μ^+μ^-)ϕ$(1020) (K$^+$K$^-$) and B$^+$$\to$ J/$ψ(μ^+μ^-)$K$^+$ decay channels. Results are presented in the transverse momentum ($p_\mathrm{T}$) range 7-50 GeV/$c$ and the rapidity interval $\lvert y \rvert$$\lt$ 2.4 for the B mesons. The measured $p_\mathrm{T}$-differential cross sections of B$^+$ and B$^0_\mathrm{s}$ in pp collisions are well described by fixed-order plus next-to-leading logarithm perturbative quantum chromodynamics calculations. Using previous PbPb collision measurements at the same nucleon-nucleon center-of-mass energy, the nuclear modification factors, $R_\mathrm{AA}$, of the B mesons are determined. For $p_\mathrm{T}$$\lt$ 10 GeV/$c$, both mesons are found to be suppressed in PbPb collisions (with $R_\mathrm{AA}$ values significantly below unity), with less suppression observed for the B$^0_\mathrm{s}$ mesons. In this $p_\mathrm{T}$ range, the $R_\mathrm{AA}$ values for the B$^+$ mesons are consistent with those for inclusive charged hadrons and D$^0$ mesons. Below 10 GeV/$c$, both B$^+$ and B$^0_\mathrm{s}$ are found to be less suppressed than either inclusive charged hadrons or D$^0$ mesons, with the B$^0_\mathrm{s}$$R_\mathrm{AA}$ value consistent with unity. The $R_\mathrm{AA}$ values found for the B$^+$ and B$^0_\mathrm{s}$ are compared to theoretical calculations, providing constraints on the mechanism of bottom quark energy loss and hadronization in the quark-gluon plasma, the hot and dense matter created in ultrarelativistic heavy ion collisions.

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 new algorithm is presented to discriminate reconstructed hadronic decays of tau leptons ($\tau_\mathrm{h}$) that originate from genuine tau leptons in the CMS detector against $\tau_\mathrm{h}$ candidates that originate from quark or gluon jets, electrons, or muons. The algorithm inputs information from all reconstructed particles in the vicinity of a $\tau_\mathrm{h}$ candidate and employs a deep neural network with convolutional layers to efficiently process the inputs. This algorithm leads to a significantly improved performance compared with the previously used one. For example, the efficiency for a genuine $\tau_\mathrm{h}$ to pass the discriminator against jets increases by 10-30% for a given efficiency for quark and gluon jets. Furthermore, a more efficient $\tau_\mathrm{h}$ reconstruction is introduced that incorporates additional hadronic decay modes. The superior performance of the new algorithm to discriminate against jets, electrons, and muons and the improved $\tau_\mathrm{h}$ reconstruction method are validated with LHC proton-proton collision data at $\sqrt{s} =$ 13 TeV.

We present the first direct search for exotic Higgs boson decays H $\to$$\mathcal{A}\mathcal{A}$, $\mathcal{A}$$\to$$\gamma\gamma$ in events with two photonlike objects. The hypothetical particle $\mathcal{A}$ is a low-mass spin-0 particle decaying promptly to a merged diphoton reconstructed as a single photonlike object. We analyze the data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 136 fb$^{-1}$. No excess above the estimated background is found. We set upper limits on the branching fraction $\mathcal{B}$(H $\to$$\mathcal{A}\mathcal{A}$$\to$ 4$\gamma$) of (0.9-3.3) $\times$ 10$^{-3}$ at 95% confidence level for masses of $\mathcal{A}$ in the range 0.1-1.2 GeV.

A search for the nonresonant production of Higgs boson pairs (HH) via gluon-gluon and vector boson fusion processes in final states with two bottom quarks and two tau leptons is presented. The search uses data from proton-proton collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events in which at least one tau lepton decays hadronically are considered and multiple machine learning techniques are used to identify and extract the signal. The data are found to be consistent, within uncertainties, with the standard model (SM) predictions. Upper limits on the HH production cross section are set to constrain the parameter space for anomalous Higgs boson couplings. The observed (expected) upper limit at 95% confidence level corresponds to 3.3 (5.2) times the SM prediction for the inclusive HH cross section and to 124 (154) times the SM prediction for the vector boson fusion HH cross section. At 95% confidence level, the Higgs field self-coupling is constrained to be within -1.7 and 8.7 times the SM expectation, and the coupling of two Higgs bosons to two vector bosons is constrained to be within -0.4 and 2.6 times the SM expectation.

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

A measurement is presented for the electroweak production of a W boson, a photon ($\gamma$), and two jets (j) in proton-proton collisions. The leptonic decay of the W boson is selected by requiring one identified electron or muon and large missing transverse momentum. The two jets are required to have large invariant dijet mass and large separation in pseudorapidity. The measurement is performed with the data collected by the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross section for the electroweak W$\gamma$jj production is 23.5 $^{+4.9}_{-4.7}$ fb, whereas the total cross section for W$\gamma$jj production is 113 $\pm$ 13 fb. Differential cross sections are also measured with the distributions unfolded to the particle level. All results are in agreement with the standard model expectations. Constraints are placed on anomalous quartic gauge couplings (aQGCs) in terms of dimension-8 effective field theory operators. These are the most stringent limits to date on the aQGCs parameters $f_\mathrm{M,2-5}$$/$$\Lambda^4$ and $f_\mathrm{T,6-7}$$/$$\Lambda^4$.