The discovery of the Higgs boson has led to new possible signatures for heavy resonance searches at the LHC. Since then, search channels including at least one Higgs boson plus another particle have formed an important part of the program of new physics searches. In this report, the status of these searches by the CMS Collaboration is reviewed. Searches are discussed for resonances decaying to two Higgs bosons, a Higgs and a vector boson, or a Higgs boson and another new resonance, with proton-proton collision data collected at $\sqrt{s}$ = 13 TeV in the years 2016-2018. A combination of the results of these searches is presented together with constraints on different beyond-the-standard model scenarios, including scenarios with extended Higgs sectors, heavy vector bosons and extra dimensions. Studies are shown for the first time by CMS on the validity of the narrow-width approximation in searches for the resonant production of a pair of Higgs bosons. The potential for a discovery at the High Luminosity LHC is also discussed.

The first search for narrow resonances decaying to three well-separated hadronic jets is presented. The search uses proton-proton collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected at the CERN LHC. No significant deviations from the background predictions are observed between 1.75-9.00 TeV. The results provide the first mass limits on a right-handed boson Z$_{\mathrm{R}}$ decaying to three gluons, an excited quark decaying via a vector boson to three quarks, as well as updated limits on a Kaluza-Klein gluon decaying via a radion to three gluons.

A test of lepton flavor universality in B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ and B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$ decays, as well as a measurement of differential and integrated branching fractions of a nonresonant B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ decay are presented. The analysis is made possible by a dedicated data set of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded in 2018, by the CMS experiment at the LHC, using a special high-rate data stream designed for collecting about 10 billion unbiased b hadron decays. The ratio of the branching fractions $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) to $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$) is determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the B$^\pm$$\to$ J/$\psi$K$^\pm$ with J/$\psi$$\to$$\mu^+\mu^-$ and e$^+$e$^-$ decays, which allow for significant cancellation of systematic uncertainties. The ratio $R$(K) is measured in the range 1.1 $\lt q^2 \lt$ 6.0 GeV$^2$, where $q$ is the invariant mass of the lepton pair, and is found to be $R$(K) = 0.78$^{+0.47}_{-0.23}$, in agreement with the standard model expectation $R$(K) $\approx$ 1. This measurement is limited by the statistical precision of the electron channel. The integrated branching fraction in the same $q^2$ range, $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) = (12.42 $\pm$ 0.68) $\times$ 10$^{-8}$, is consistent with the present world-average value and has a comparable precision.

Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$)/$\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).

The first evidence for the standard model production of a top quark in association with a W boson and a Z boson is reported. The measurement is performed in multilepton final states, where the Z boson is reconstructed via its decays to electron or muon pairs and the W boson decays either to leptons or hadrons. The analysed data were recorded by the CMS experiment at the CERN LHC in 2016-2018 in proton-proton collisions at $\sqrt{s}$ = 13 TeV, and correspond to an integrated luminosity of 138 fb$^{-1}$. The measured cross section is 354 $\pm$ 54 (stat) $\pm$ 95 (syst) fb, and corresponds to a statistical significance of 3.4 standard deviations.

A study of the anomalous couplings of the Higgs boson to vector bosons, including $CP$-violation effects, has been conducted using its production and decay in the WW channel. This analysis is performed on proton-proton collision data collected with the CMS detector at the CERN LHC during 2016-2018 at a center-of-mass energy of 13 TeV, and corresponds to an integrated luminosity of 138 fb$^{-1}$. The different-flavor dilepton (e$\mu$) final state is analyzed, with dedicated categories targeting gluon fusion, electroweak vector boson fusion, and associated production with a W or Z boson. Kinematic information from associated jets is combined using matrix element techniques to increase the sensitivity to anomalous effects at the production vertex. A simultaneous measurement of four Higgs boson couplings to electroweak vector bosons is performed in the framework of a standard model effective field theory. All measurements are consistent with the expectations for the standard model Higgs boson and constraints are set on the fractional contribution of the anomalous couplings to the Higgs boson production cross section.

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 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 long-lived particles (LLPs) decaying in the CMS muon detectors is presented. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded at the LHC in 2016-2018, is used. The decays of LLPs are reconstructed as high multiplicity clusters of hits in the muon detectors. In the context of twin Higgs models, the search is sensitive to LLP masses from 0.4 to 55 GeV and a broad range of LLP decay modes, including decays to hadrons, $\tau$ leptons, electrons, or photons. No excess of events above the standard model background is observed. The most stringent limits to date from LHC data are set on the branching fraction of the Higgs boson decay to a pair of LLPs with masses below 10 GeV. This search also provides the best limits for various intervals of LLP proper decay length and mass. Finally, this search sets the first limits at the LHC on a dark quantum chromodynamic sector whose particles couple to the Higgs boson through gluon, Higgs boson, photon, vector, and dark-photon portals, and is sensitive to branching fractions of the Higgs boson to dark quarks as low as 2$\times$10$^{-3}$.

A search for pair production of scalar and vector leptoquarks (LQs) each decaying to a muon and a bottom quark is performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV with the CMS detector at the CERN LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess above standard model expectation is observed. Scalar (vector) LQs with masses less than 1810 (2120) GeV are excluded at 95% confidence level, assuming a 100% branching fraction of the LQ decaying to a muon and a bottom quark. These limits represent the most stringent to date.