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 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 is described for the production of a pair of bottom-type vector-like quarks (VLQs), each decaying into a b or $\mathrm{\bar{b}}$ quark and either a Higgs or a Z boson, with a mass greater than 1000 GeV. The analysis is based on data from proton-proton collisions at a 13 TeV center-of-mass energy recorded at the CERN LHC, corresponding to a total integrated luminosity of 137 fb$^{-1}$. As the predominant decay modes of the Higgs and Z bosons are to a pair of quarks, the analysis focuses on final states consisting of jets resulting from the six quarks produced in the events. Since the two jets produced in the decay of a highly Lorentz-boosted Higgs or Z boson can merge to form a single jet, nine independent analyses are performed, categorized by the number of observed jets and the reconstructed event mode. No signal in excess of the expected background is observed. Lower limits are set on the VLQ mass at 95% confidence level equal to 1570 GeV in the case where the VLQ decays exclusively to a b quark and a Higgs boson, 1390 GeV for when it decays exclusively to a b quark and a Z boson, and 1450 GeV for when it decays equally in these two modes. These limits represent significant improvements over the previously published VLQ limits.

A generic search is presented for the associated production of a Z boson or a photon with an additional unspecified massive particle X, pp $\to$ pp + Z/$\gamma$ + X, in proton-tagged events from proton-proton collisions at $\sqrt{s}$ = 13 TeV, recorded in 2017 with the CMS detector and the CMS-TOTEM precision proton spectrometer. The missing mass spectrum is analysed in the 600-1600 GeV range and a fit is performed to search for possible deviations from the background expectation. No significant excess in data with respect to the background predictions has been observed. Model-independent upper limits on the visible production cross section of pp $\to$ pp + Z/$\gamma$ + X are set.

The first measurement of the CP structure of the Yukawa coupling between the Higgs boson and $\tau$ leptons is presented. The measurement is based on data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. The analysis uses the angular correlation between the decay planes of $\tau$ leptons produced in Higgs boson decays. The effective mixing angle between CP-even and CP-odd $\tau$ Yukawa couplings is found to be $-$1 $\pm$ 19$^\circ$, compared to an expected value of 0 $\pm$ 21$^\circ$ at the 68.3% confidence level. The data disfavour the pure CP-odd scenario at 3.0 standard deviations. The results are compatible with predictions for the standard model Higgs boson.

Azimuthal correlations between the two jets with the largest transverse momenta $ {p_{\mathrm{T}}} $ in inclusive 2-, 3-, and 4-jet events are presented for several regions of the leading jet $ {p_{\mathrm{T}}} $ up to 4 TeV. For 3- and 4-jet scenarios, measurements of the minimum azimuthal angles between any two of the three or four leading $ {p_{\mathrm{T}}} $ jets are also presented. The analysis is based on data from proton-proton collisions collected by the CMS Collaboration at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Calculations based on leading-order matrix elements supplemented with parton showering and hadronization do not fully describe the data, so next-to-leading-order calculations matched with parton shower and hadronization models are needed to better describe the measured distributions. Furthermore, we show that azimuthal jet correlations are sensitive to details of the parton showering, hadronization, and multiparton interactions. A next-to-leading-order calculation matched with parton showers in the MC@NLO method, as implemented in HERWIG 7, gives a better overall description of the measurements than the POWHEG method.

The production of Z bosons associated with jets is measured in pp collisions at $\sqrt{s}$ = 13 TeV with data recorded with the CMS experiment at the LHC corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV is measured for different regions of the Z boson's $p_\mathrm{T}$(Z), from lower than 10 GeV to higher than 100 GeV. The azimuthal correlation $\Delta \phi$ between the Z boson and the leading jet, as well as the correlations between the two leading jets are measured in three regions of $p_\mathrm{T}$(Z). The measurements are compared with several predictions at leading and next-to-leading orders, interfaced with parton showers. Predictions based on transverse-momentum dependent parton distributions and corresponding parton showers give a good description of the measurement in the regions where multiple parton interactions and higher jet multiplicities are not important. The effects of multiple parton interactions are shown to be important to correctly describe the measured spectra in the low $p_\mathrm{T}$(Z) regions.

A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018. The analyzed data correspond to an integrated luminosity of up to 137 fb$^{-1}$. The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage.

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.

The dependence of inclusive jet production in proton-proton collisions with a center-of-mass energy of 13 TeV on the distance parameter $R$ of the anti-$k_\mathrm{T}$ algorithm is studied using data corresponding to integrated luminosities up to 35.9 fb$^{-1}$ collected by the CMS experiment in 2016. The ratios of the inclusive cross sections as functions of transverse momentum $p_\mathrm{T}$ and rapidity $y$, for $R$ in the range 0.1 to 1.2 to those using $R = $ 0.4 are presented in the region 84 $\lt$ $p_\mathrm{T}$ $\lt$ 1588 GeV and $|y|$ $\lt$ 2.0. The results are compared to calculations at leading and next-to-leading order in the strong coupling constant using different parton shower models. The variation of the ratio of cross sections with $R$ is well described by calculations including a parton shower model, but not by a leading-order quantum chromodynamics calculation including nonperturbative effects. The agreement between the data and the theoretical predictions for the ratios of cross sections is significantly improved when next-to-leading order calculations with nonperturbative effects are used.