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 new phenomena is performed in final states containing one or more jets and an imbalance in transverse momentum in pp collisions at a centre-of-mass energy of 13 TeV. The analysed data sample, recorded with the CMS detector at the CERN LHC, corresponds to an integrated luminosity of 2.3 inverse femtobarns. Several kinematic variables are employed to suppress the dominant background, multijet production, as well as to discriminate between other standard model and new physics processes. The search provides sensitivity to a broad range of new-physics models that yield a stable weakly interacting massive particle. The number of observed candidate events is found to agree with the expected contributions from standard model processes, and the result is interpreted in the mass parameter space of fourteen simplified supersymmetric models that assume the pair production of gluinos or squarks and a range of decay modes. For models that assume gluino pair production, masses up to 1575 and 975 GeV are excluded for gluinos and neutralinos, respectively. For models involving the pair production of top squarks and compressed mass spectra, top squark masses up to 400 GeV are excluded.

The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

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.

A measurement for inclusive 2- and 3-jet events of the azimuthal correlation between the two jets with the largest transverse momenta, $\Delta\phi_{12}$ , is presented. The measurement considers events where the two leading jets are nearly collinear ("back-to-back") in the transverse plane and is performed for several ranges of the leading jet transverse momentum. Proton-proton collision data collected with the CMS experiment at a center-of-mass energy of 13 TeV and corresponding to an integrated luminosity of 35.9 fb$^{-1}$ are used. Predictions based on calculations using matrix elements at leading-order and next-to-leading-order accuracy in perturbative quantum chromodynamics supplemented with leading-log parton showers and hadronization are generally in agreement with the measurements. Discrepancies between the measurement and theoretical predictions are as large as 15%, mainly in the region 177$^\circ$ $<$ $\Delta\phi_{12}$ $<$ 180$^\circ$. The 2- and 3-jet measurements are not simultaneously described by any of models.

Measurements of primary charged hadron multiplicity distributions are presented for non-single-diffractive events in proton-proton collisions at centre-of-mass energies of sqrt(s) = 0.9, 2.36, and 7 TeV, in five pseudorapidity ranges from |eta|<0.5 to |eta|<2.4. The data were collected with the minimum-bias trigger of the CMS experiment during the LHC commissioning runs in 2009 and the 7 TeV run in 2010. The multiplicity distribution at sqrt(s) = 0.9 TeV is in agreement with previous measurements. At higher energies the increase of the mean multiplicity with sqrt(s) is underestimated by most event generators. The average transverse momentum as a function of the multiplicity is also presented. The measurement of higher-order moments of the multiplicity distribution confirms the violation of Koba-Nielsen-Olesen scaling that has been observed at lower energies.

Azimuthal correlations of charged particles in xenon-xenon collisions at a center-of-mass energy per nucleon pair of $ \sqrt{s_{_\mathrm{NN}}} =$ 5.44 TeV are studied. The data were collected by the CMS experiment at the LHC with a total integrated luminosity of 3.42 $\mu$b$^{-1}$. The collective motion of the system formed in the collision is parameterized by a Fourier expansion of the azimuthal particle density distribution. The azimuthal anisotropy coefficients $v_{2}$, $v_{3}$, and $v_{4}$ are obtained by the scalar-product, two-particle correlation, and multiparticle correlation methods. Within a hydrodynamic picture, these methods have different sensitivities to non-collective and fluctuation effects. The dependence of the Fourier coefficients on the size of the colliding system is explored by comparing the xenon-xenon results with equivalent lead-lead data. Model calculations that include initial-state fluctuation effects are also compared to the experimental results. The observed angular correlations provide new constraints on the hydrodynamic description of heavy ion collisions.

A combination of searches for top squark pair production using proton-proton collision data at a center-of-mass energy of 13 TeV at the CERN LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment, is presented. Signatures with at least 2 jets and large missing transverse momentum are categorized into events with 0, 1, or 2 leptons. New results for regions of parameter space where the kinematical properties of top squark pair production and top quark pair production are very similar are presented. Depending on the model, the combined result excludes a top squark mass up to 1325 GeV for a massless neutralino, and a neutralino mass up to 700 GeV for a top squark mass of 1150 GeV. Top squarks with masses from 145 to 295 GeV, for neutralino masses from 0 to 100 GeV, with a mass difference between the top squark and the neutralino in a window of 30 GeV around the mass of the top quark, are excluded for the first time with CMS data. The results of theses searches are also interpreted in an alternative signal model of dark matter production via a spin-0 mediator in association with a top quark pair. Upper limits are set on the cross section for mediator particle masses of up to 420 GeV.

The pseudorapidity distributions of dijets as a function of their average transverse momentum ($p_\mathrm{T}^\text{ave}$) are measured in proton-lead (pPb) and proton-proton (pp) collisions. The data samples were collected by the CMS experiment at the CERN LHC, at a nucleon-nucleon center-of-mass energy of 5.02 TeV. A significant modification of the pPb spectra with respect to the pp spectra is observed in all $p_\mathrm{T}^\text{ave}$ intervals investigated. The ratios of the pPb and pp distributions are compared to next-to-leading order perturbative quantum chromodynamics calculations with unbound nucleon and nuclear parton distribution functions (PDFs). These results give the first evidence that the gluon PDF at large Bjorken $x$ in lead ions is strongly suppressed with respect to the PDF in unbound nucleons.

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.