In July 2012, the ATLAS and CMS Collaborations at the CERN Large Hadron Collider announced the observation of a Higgs boson at a mass of around 125 GeV. Ten years later, and with the data corresponding to the production of 30 times larger number of Higgs bosons, we have learnt much more about the properties of the Higgs boson. The CMS experiment has observed the Higgs boson in numerous fermionic and bosonic decay channels, established its spin-parity quantum numbers, determined its mass and measured its production cross sections in various modes. Here the CMS Collaboration reports the most up-to-date combination of results on the properties of the Higgs boson, including the most stringent limit on the cross section for the production of a pair of Higgs bosons, on the basis of data from proton-proton collisions at a centre-of-mass energy of 13 TeV. Within the uncertainties, all these observations are compatible with the predictions of the standard model of elementary particle physics. Much evidence points to the fact that the standard model is a low-energy approximation of a more comprehensive theory. Several of the standard model issues originate in the sector of Higgs boson physics. An order of magnitude larger number of Higgs bosons, expected to be examined over the next fifteen years, will help deepen our understanding of this crucial sector.

A search for Higgs boson pair (HH) production with one Higgs boson decaying to two bottom quarks and the other to two W bosons are presented. The search is done using proton-proton collisions data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded by the CMS detector at the LHC from 2016 to 2018. The final states considered include at least one leptonically decaying W boson. No evidence for the presence of a signal is observed and corresponding upper limits on the HH production cross section are derived. The limit on the inclusive cross section of the nonresonant HH production, assuming that the distributions of kinematic observables are as expected in the standard model (SM), is observed (expected) to be 14 (18) times the value predicted by the SM, at 95% confidence level. The limits on the cross section are also presented as functions of various Higgs boson coupling modifiers, and anomalous Higgs boson coupling scenarios. In addition, limits are set on the resonant HH production via spin-0 and spin-2 resonances within the mass range 250-900 GeV.

A search for an exotic decay of the Higgs boson to a Z boson and a light pseudoscalar particle (a), decaying to a pair of leptons and a pair of photons, respectively, is presented. The search is based on proton-proton collision data at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, collected with the CMS detector and corresponding to an integrated luminosity of 138 fb$^{-1}$. The analysis probes pseudoscalar masses $m_\mathrm{a}$ between 1 and 30 GeV, leading to two pairs of well-isolated leptons and photons. Upper limits at 95% confidence level are set on the Higgs boson production cross section times its branching fraction to two leptons and two photons. The observed (expected) limits are in the range of 1.1-17.8 (1.7-17.9) fb within the probed $m_\mathrm{a}$ interval. An excess of data above the expected standard model background with a local (global) significance of 2.6 (1.3) standard deviations is observed for a mass hypothesis of $m_\mathrm{a}$ = 3 GeV. Limits on models involving axion-like particles, formulated as an effective field theory, are also reported.

Diboson production in association with jets is studied in the fully leptonic final states, pp $\to$ (Z$\gamma^*$)(Z/$\gamma^*$) + jets $\to$ 2$\ell$2$\ell'$ + jets, ($\ell,\ell'$ = e or $\mu$) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. Differential distributions and normalized differential cross sections are measured as a function of jet multiplicity, transverse momentum $p_\mathrm{T}$, pseudorapidity $\eta$, invariant mass and $\Delta\eta$ of the highest-$p_\mathrm{T}$ and second-highest-$p_\mathrm{T}$ jets, and as a function of invariant mass of the four-lepton system for events with various jet multiplicities. These differential cross sections are compared with theoretical predictions that mostly agree with the experimental data. However, in a few regions we observe discrepancies between the predicted and measured values. Further improvement of the predictions is required to describe the ZZ + jets production in the whole phase space.

A statistical combination of searches for heavy resonances decaying to pairs of bosons or leptons is presented. The data correspond to an integrated luminosity of 35.9 fb$^{-1}$ collected during 2016 by the CMS experiment at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV. The data are found to be consistent with expectations from the standard model background. Exclusion limits are set in the context of models of spin-1 heavy vector triplets and of spin-2 bulk gravitons. For mass-degenerate W' and Z' resonances that predominantly couple to the standard model gauge bosons, the mass exclusion at 95% confidence level of heavy vector bosons is extended to 4.5 TeV as compared to 3.8 TeV determined from the best individual channel. This excluded mass increases to 5.0 TeV if the resonances couple predominantly to fermions.

The first search for singly produced 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 and on 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 study of the production of prompt J/$\psi$ mesons contained in jets in proton-proton collisions at $\sqrt{s} =$ 8 TeV is presented. The analysis is based on data corresponding to an integrated luminosity of 19.1 fb$^{-1}$ collected with the CMS detector at the LHC. For events with at least one observed jet, the angular separation between the J/$\psi$ meson and the jet is used to test whether the J/$\psi$ meson is part of the jet. The analysis shows that most prompt J/$\psi$ mesons with energy above 15 GeV and rapidity $|y|<$ 1 are contained in jets with pseudorapidity $|\eta_{\text{jet}}|$ $<$ 1. The differential distributions of the probability to have a J/$\psi$ meson contained in a jet as a function of jet energy for a fixed J/$\psi$ energy fraction are compared to a theoretical model using the fragmenting jet function approach. The data agree best with fragmenting jet function calculations that use a long-distance matrix element parameter set in which prompt J/$\psi$ mesons are predicted to be unpolarized. This technique demonstrates a new way to test predictions for prompt J/$\psi$ production using nonrelativistic quantum chromodynamics.

A search for heavy neutral leptons (HNLs) decaying in the CMS muon system is presented. A data sample is used corresponding to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}$ = 13 TeV, recorded at the CERN LHC in 2016-2018. Decay products of long-lived HNLs could interact with the shielding materials in the CMS muon system and create hadronic and electromagnetic showers detected in the muon chambers. This distinctive signature provides a unique handle to search for HNLs with masses below 4 GeV and proper decay lengths of the order of meters. The signature is sensitive to HNL couplings to all three generations of leptons. Candidate events are required to contain a prompt electron or muon originating from a vertex on the beam axis and a displaced shower in the muon chambers. No significant deviations from the standard model background expectation are observed. In the electron (muon) channel, the most stringent limits to date are set for HNLs in the mass range of 2.1-3.0 (1.9-3.3) GeV, reaching mixing matrix element squared values as low as 8.6 (4.6) $\times$ 10$^{-6}$.

A search for high-mass dimuon resonance production in association with one or more b quark jets is presented. The study uses proton-proton collision data collected with the CMS detector at the LHC corresponding to an integrated luminosity of 138 fb$^{-1}$ at a center-of-mass energy of 13 TeV. Model-independent limits are derived on the number of signal events with exactly one or more than one b quark jet. Results are also interpreted in a lepton-flavor-universal model with Z$'$ boson couplings to a bb quark pair ($g_\mathrm{b}$), an sb quark pair ($g_\mathrm{b}\delta_\mathrm{bs}$), and any same-flavor charged lepton ($g_\ell$) or neutrino pair ($g_\nu$), with $\left|g_{\nu}\right| = \left|g_\ell\right|$. For a Z$'$ boson with a mass $m_{\mathrm{Z}'}$ = 350 GeV (2 TeV) and $\left|\delta_\mathrm{bs}\right|$$\lt$ 0.25, the majority of the parameter space with 0.0057 $\lt$$\left|g_\ell\right|$$\lt$ 0.35 (0.25 $\lt$$\left|g_\ell\right|$$\lt$ 0.43) and 0.0079 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.46 (0.34 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.57) is excluded at 95% confidence level. Finally, constraints are set on a Z$'$ model with parameters consistent with low-energy b $\to$ s$\ell\ell$ measurements. In this scenario, most of the allowed parameter space is excluded for a Z$'$ boson with 350 $\lt m_{\mathrm{Z}'}$ $\lt$ 500 GeV, while the constraints are less stringent for higher $m_{\mathrm{Z}'}$ hypotheses. This is the first dedicated search at the LHC for a high-mass dimuon resonance produced in association with multiple b quark jets, and the constraints obtained on models with this signature are the most stringent to date.

A measurement of the dijet production cross section is reported based on proton-proton collision data collected in 2016 at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of up to 36.3 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm for distance parameters of $R$ = 0.4 and 0.8. Cross sections are measured double-differentially (2D) as a function of the largest absolute rapidity $\lvert y_\text{max}\rvert$ of the two jets with the highest transverse momenta $p_\mathrm{T}$ and their invariant mass $m_{1,2}$, and triple-differentially (3D) as a function of the rapidity separation $y^*$, the total boost $y_\mathrm{b}$, and either $m_{1,2}$ or the average $p_\mathrm{T}$ of the two jets. The cross sections are unfolded to correct for detector effects and are compared with fixed-order calculations derived at next-to-next-to-leading order in perturbative quantum chromodynamics. The impact of the measurements on the parton distribution functions and the strong coupling constant at the mass of the Z boson is investigated, yielding a value of $\alpha_\mathrm{S}(m_\mathrm{Z})$ = 0.1179 $\pm$ 0.0019.