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. All analyses use 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.

A search is presented for an extended Higgs sector with two new particles, X and $\phi$, in the process $X \to\phi\phi\to(\gamma\gamma)(\gamma\gamma)$. Novel neural networks classify events with diphotons that are merged and determine the diphoton masses. The search uses LHC proton-proton collision data at $\sqrt{s}$ = 13 TeV collected with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. No evidence of such resonances is seen. Upper limits are set on the production cross section for $m_X$ between 300 and 3000 GeV and $m_\phi / m_X$ between 0.5 and 2.5%, representing the most sensitive search in this channel.

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.

A search for the rare $\eta$$\to$$\mu^+\mu^-\mu^+\mu^-$ double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb$^{-1}$. A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the $\eta$$\to$$\mu^+ \mu^-$ decay as normalization, the branching fraction $\mathcal{B}(\eta$$\to$$\mu^+\mu^-\mu^+\mu^-)$ = [5.0 $\pm$ 0.8 (stat) $\pm$ 0.7 (syst) $\pm$ 0.7 ($\mathcal{B}_{2\mu}$)] $\times$ 10$^{-9}$ is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.

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.

A search for the production of long-lived particles in proton-proton collisions at a center-of-mass energy of 13 TeV at the CERN LHC is presented. The search is based on data collected by the CMS experiment in 2016-2018, corresponding to a total integrated luminosity of 137 fb$^{-1}$. This search is designed to be sensitive to long-lived particles with mean proper decay lengths between 0.1 and 1000 mm, whose decay products produce a final state with at least one displaced vertex and missing transverse momentum. A machine learning algorithm, which improves the background rejection power by more than an order of magnitude, is applied to improve the sensitivity. The observation is consistent with the standard model background prediction, and the results are used to constrain split supersymmetry (SUSY) and gauge-mediated SUSY breaking models with different gluino mean proper decay lengths and masses. This search is the first CMS search that shows sensitivity to hadronically decaying long-lived particles from signals with mass differences between the gluino and neutralino below 100 GeV. It sets the most stringent limits to date for split-SUSY models and gauge-mediated SUSY breaking models with gluino proper decay length less than 6 mm.

Measurements of fiducial and total inclusive cross sections for W and Z boson production are presented in proton-proton collisions at $\sqrt{s}$ = 5.02 and 13 TeV. Electron and muon decay modes ($\ell$ = e or $\mu$) are studied in the data collected with the CMS detector in 2017, in dedicated runs with reduced instantaneous luminosity. The data sets correspond to integrated luminosities of 298 $\pm$ 6 pb$^{-1}$ at 5.02 TeV and 206 $\pm$ 5 pb$^{-1}$ at 13 TeV. Measured values of the products of the total inclusive cross sections and the branching fractions at 5.02 TeV are $\sigma$(pp $\to$ W+X) $\mathcal{B}$(W $\to$$\ell\nu$) = 7300 $\pm$ 10 (stat) $\pm$ 60 (syst) $\pm$ 140 (lumi) pb, and $\sigma$(pp $\to$ Z+X) $\mathcal{B}$(Z $\to$$\ell^+\ell^-$) = 669 $\pm$ 2 (stat) $\pm$ 6 (syst) $\pm$ 13 (lumi) pb for the dilepton invariant mass in the range of 60-120 GeV. The corresponding results at 13 TeV are 20480 $\pm$ 10 (stat) $\pm$ 170 (syst) $\pm$ 470 (lumi) pb and 1952 $\pm$ 4 (stat) $\pm$ 18 (syst) $\pm$ 45 (lumi) pb. The measured values agree with cross section calculations at next-to-next-to-leading-order in perturbative quantum chromodynamics. Fiducial and total inclusive cross sections, ratios of cross sections of W$^+$ and W$^-$ production as well as inclusive W and Z boson production, and ratios of these measurements at 5.02 and 13 TeV are reported.

Nuclear medium effects on B$^+$ meson production are studied using the binary-collision scaled cross section ratio between events of different charged-particle multiplicities from proton-lead collisions. Data, collected by the CMS experiment in 2016 at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV, corresponding to an integrated luminosity of 175 nb$^{-1}$, were used. The scaling factors in the ratio are determined using a novel approach based on the Z $\to$$\mu^+\mu^-$ cross sections measured in the same events. The scaled ratio for B$^+$ is consistent with unity for all event multiplicities, putting stringent constraints on nuclear modification for heavy flavor.

Incoherent J/$\psi$ photoproduction in heavy ion ultraperipheral collisions (UPCs), in which the photon interacts with localized, fluctuating gluonic hotspots rather than the entire nucleus, provides a unique probe of those fluctuations. This study presents the first measurement of the dependence of this photoproduction yield on the photon-nucleon center-of-mass energy ($W_\gamma^\mathrm{N}$), using PbPb UPCs at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The data corresponds to an integrated luminosity of 1.52 nb$^{-1}$, recorded by the CMS experiment. The measurement covers a wide $W_\gamma^\mathrm{N}$ range of $\approx$40-400 GeV, probing gluons carrying a fraction of nucleon momentum $x$ in the range between 5.9 $\times$ 10$^{-3}$ and 6.5 $\times$ 10$^{-5}$. The measured incoherent J/$\psi$ photoproduction cross section is suppressed relative to theoretical predictions without nuclear effects. However, the ratio of incoherent to coherent photoproduction remains constant across the probed $W_\gamma^\mathrm{N}$ and $x$ range. Together, these results pose significant challenges to current theoretical models which include gluon saturation or nuclear shadowing effects.

A measurement is performed of Higgs bosons produced with high transverse momentum ($p_\mathrm{T}$) via vector boson or gluon fusion in proton-proton collisions. The result is based on a data set with a center-of-mass energy of 13 TeV collected in 2016-2018 with the CMS detector at the LHC and corresponds to an integrated luminosity of 138 fb$^{-1}$. The decay of a high-$p_\mathrm{T}$ Higgs boson to a boosted bottom quark-antiquark pair is selected using large-radius jets and employing jet substructure and heavy-flavor taggers based on machine learning techniques. Independent regions targeting the vector boson and gluon fusion mechanisms are defined based on the topology of two quark-initiated jets with large pseudorapidity separation. The signal strengths for both processes are extracted simultaneously by performing a maximum likelihood fit to data in the large-radius jet mass distribution. The observed signal strengths relative to the standard model expectation are 4.9 $^{+1.9}_{-1.6}$ and 1.6 $^{+1.7}_{-1.5}$ for the vector boson and gluon fusion mechanisms, respectively. A differential cross section measurement is also reported in the simplified template cross section framework.