The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.

Results are presented from a search for CP violation in top quark pair production, using proton-proton collisions at a center-of-mass energy of 13 TeV. The data used for this analysis consist of final states with two charged leptons collected by the CMS experiment, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The search uses two observables, $\mathcal{O}_1$ and $\mathcal{O}_3$, which are Lorentz scalars. The observable $\mathcal{O}_1$ is constructed from the four-momenta of the charged leptons and the reconstructed top quarks, while $\mathcal{O}_3$ consists of the four-momenta of the charged leptons and the b quarks originating from the top quarks. Asymmetries in these observables are sensitive to CP violation, and their measurement is used to determine the chromoelectric dipole moment of the top quark. The results are consistent with the expectation from the standard model.

The path-length dependent parton energy loss within the dense partonic medium created in lead-lead collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV is studied by determining the azimuthal anisotropies for dijets with high transverse momentum. The data were collected by the CMS experiment in 2018 and correspond to an integrated luminosity of 1.69 nb$^{-1}$. For events containing back-to-back jets, correlations in relative azimuthal angle and pseudorapidity ($\eta$) between jets and hadrons, and between two hadrons, are constructed. The anisotropies are expressed as the Fourier expansion coefficients $v_n$, $n = $ 2-4 of these azimuthal distributions. The dijet $v_n$ values are extracted from long-range (1.5 $\lt$$\vert\Delta\eta\vert$$\lt$ 2.5) components of these correlations, which suppresses the background contributions from jet fragmentation processes. Positive dijet $v_2$ values are observed which increase from central to more peripheral events, while the $v_3$ and $v_4$ values are consistent with zero within experimental uncertainties.

The CMS experiment at the LHC has measured the differential cross sections of Z bosons decaying to pairs of leptons, as functions of transverse momentum and rapidity, in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The measured Z boson elliptic azimuthal anisotropy coefficient is compatible with zero, showing that Z bosons do not experience significant final-state interactions in the medium produced in the collision. Yields of Z bosons are compared to Glauber model predictions and are found to deviate from these expectations in peripheral collisions, indicating the presence of initial collision geometry and centrality selection effects. The precision of the measurement allows, for the first time, for a data-driven determination of the nucleon-nucleon integrated luminosity as a function of lead-lead centrality, thereby eliminating the need for its estimation based on a Glauber model.

The spectra measurements of charged hadrons produced in the shower of a parton originating in the same hard scattering with a leptonically decaying Z boson, are reported in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. Both PbPb and pp data sets are recorded by the CMS experiment at the LHC, and correspond to an integrated luminosity of 1.7 nb$^{-1}$ and 320 pb$^{-1}$, respectively. Hadronic collision data with one reconstructed Z boson candidate with the transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV/$c$ are analyzed. The Z boson constrains the initial energy and direction of the associated parton. In heavy ion events, azimuthal angular distributions of charged hadrons with respect to the direction of a Z boson are sensitive to modifications of the in-medium parton shower and medium response. Compared to reference data from pp interactions, the results for central PbPb collisions indicate a modification of the angular correlations. The measurements of the fragmentation functions and $p_\mathrm{T}$ spectra of charged particles in Z boson events, which are sensitive to medium modifications of the parton shower longitudinal structure, are also reported. Significant modifications in central PbPb events compared to pp reference data are also found for these observables.

A search for pairs of dijet resonances with the same mass is conducted in final states with at least four jets. Results are presented separately for the case where the four jet production proceeds via an intermediate resonant state and for nonresonant production. The search uses a data sample corresponding to an integrated luminosity of 138 fb$^{-1}$ collected by the CMS detector in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Model-independent limits, at 95% confidence level, are reported on the production cross section of four-jet and dijet resonances. These first LHC limits on resonant pair production of dijet resonances via high mass intermediate states are applied to a signal model of diquarks that decay into pairs of vector-like quarks, excluding diquark masses below 7.6 TeV for a particular model scenario. There are two events in the tails of the distributions, each with a four-jet mass of 8 TeV and an average dijet mass of 2 TeV, resulting in local and global significances of 3.9 and 1.6 standard deviations, respectively, if interpreted as a signal. The nonresonant search excludes pair production of top squarks with masses between 0.50 TeV to 0.77 TeV, with the exception of a small interval between 0.52 and 0.58 TeV, for supersymmetric $R$-parity-violating decays to quark pairs, significantly extending previous limits. Here, the most significant excess above the predicted background occurs at an average dijet mass of 0.95 TeV, for which the local and global significances are 3.6 and 2.5 standard deviations, respectively.

An inclusive search for long-lived exotic particles decaying to a pair of muons is presented. The search uses data collected by the CMS experiment at the CERN LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV in 2016 and 2018 and corresponding to an integrated luminosity of 97.6 fb$^{-1}$. The experimental signature is a pair of oppositely charged muons originating from a common secondary vertex spatially separated from the pp interaction point by distances ranging from several hundred $\mu$m to several meters. The results are interpreted in the frameworks of the hidden Abelian Higgs model, in which the Higgs boson decays to a pair of long-lived dark photons Z$_\mathrm{D}$, and of a simplified model, in which long-lived particles are produced in decays of an exotic heavy neutral scalar boson. For the hidden Abelian Higgs model with $m_\mathrm{Z_D}$ greater than 20 GeV and less than half the mass of the Higgs boson, they provide the best limits to date on the branching fraction of the Higgs boson to dark photons for $c\tau$(Z$_\mathrm{D}$) (varying with $m_\mathrm{Z_D}$) between 0.03 and ${\approx}$ 0.5 mm, and above ${\approx}$ 0.5 m. Our results also yield the best constraints on long-lived particles with masses larger than 10 GeV produced in decays of an exotic scalar boson heavier than the Higgs boson and decaying to a pair of muons.

A search is presented for long-lived particles produced in pairs in proton-proton collisions at the LHC operating at a center-of-mass energy of 13 TeV. The data were collected with the CMS detector during the period from 2015 through 2018, and correspond to a total integrated luminosity of 140 fb$^{-1}$. This search targets pairs of long-lived particles with mean proper decay lengths between 0.1 and 100 mm, each of which decays into at least two quarks that hadronize to jets, resulting in a final state with two displaced vertices. No significant excess of events with two displaced vertices is observed. In the context of $R$-parity violating supersymmetry models, the pair production of long-lived neutralinos, gluinos, and top squarks is excluded at 95% confidence level for cross sections larger than 0.08 fb, masses between 800 and 3000 GeV, and mean proper decay lengths between 1 and 25 mm.

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.

Pseudorapidity (eta) distributions of charged particles produced in proton-proton collisions at a centre-of-mass energy of 8 TeV are measured in the ranges abs(eta) < 2.2 and 5.3 < abs(eta) < 6.4 covered by the CMS and TOTEM detectors, respectively. The data correspond to an integrated luminosity of 45 inverse microbarns. Measurements are presented for three event categories. The most inclusive category is sensitive to 91-96% of the total inelastic proton-proton cross section. The other two categories are disjoint subsets of the inclusive sample that are either enhanced or depleted in single diffractive dissociation events. The data are compared to models used to describe high-energy hadronic interactions. None of the models considered provide a consistent description of the measured distributions.

A search is presented for vector-like T and B quark-antiquark pairs produced in proton-proton collisions at a center-of-mass energy of 13 TeV. Data were collected by the CMS experiment at the CERN LHC in 2016-2018, with an integrated luminosity of 138 fb$^{-1}$. Events are separated into single-lepton, same-sign charge dilepton, and multilepton channels. In the analysis of the single-lepton channel a multilayer neural network and jet identification techniques are employed to select signal events, while the same-sign dilepton and multilepton channels rely on the high-energy signature of the signal to distinguish it from standard model backgrounds. The data are consistent with standard model background predictions, and the production of vector-like quark pairs is excluded at 95% confidence level for T quark masses up to 1.54 TeV and B quark masses up to 1.56 TeV, depending on the branching fractions assumed, with maximal sensitivity to decay modes that include multiple top quarks. The limits obtained in this search are the strongest limits to date for $\mathrm{T\overline{T}}$ production, excluding masses below 1.48 TeV for all decays to third generation quarks, and are the strongest limits to date for $\mathrm{B\overline{B}}$ production with B quark decays to tW.

The charge-parity ($CP$) structure of the Yukawa interaction between the Higgs (H) boson and the top quark is measured in a data sample enriched in the $\mathrm{t\bar{t}}$ and tH associated production, using 138 fb$^{-1}$ of data collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC. The study targets events where the H boson decays via H $\to$ WW or H $\to$$\tau\tau$ and the top quarks decay via t $\to$ Wb: the W bosons decay either leptonically or hadronically, and final states characterized by the presence of at least two leptons are studied. Machine learning techniques are applied to these final states to enhance the separation of $CP$-even from $CP$-odd scenarios. Two-dimensional confidence regions are set on $\kappa_\mathrm{t}$ and $\tilde{\kappa}_\mathrm{t}$, which are respectively defined as the $CP$-even and $CP$-odd top-Higgs Yukawa coupling modifiers. No significant fractional $CP$-odd contributions, parameterized by the quantity $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ are observed; the parameter is determined to be $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ = 0.59 with an interval of (0.24, 0.81) at 68% confidence level. The results are combined with previous results covering the H $\to$ ZZ and H $\to$ $\gamma\gamma$ decay modes, yielding two- and one-dimensional confidence regions on $\kappa_\mathrm{t}$ and $\tilde{\kappa}_\mathrm{t}$, while $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ is determined to be $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ = 0.28 with an interval of $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ $\lt$ 0.55 at 68% confidence level, in agreement with the standard model $CP$-even prediction of $\lvert{f_{CP}^{\mathrm{Htt}}}\rvert$ = 0.

Results are presented from a search for the Higgs boson decay H $\to$ Z$\gamma$, where Z $\to$$\ell^+\ell^-$ with $\ell$ = e or $\mu$. The search is performed using a sample of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events are assigned to mutually exclusive categories, which exploit differences in both event topology and kinematics of distinct Higgs production mechanisms to enhance signal sensitivity. The signal strength $\mu$, defined as the product of the cross section and the branching fraction [$\sigma($pp $\to$ H$)\mathcal{B}($H $\to$ Z$\gamma)$] relative to the standard model prediction, is extracted from a simultaneous fit to the $\ell^+\ell^-\gamma$ invariant mass distributions in all categories and is found to be $\mu$ = 2.4 $\pm$ 0.9 for a Higgs boson mass of 125.38 GeV. The statistical significance of the observed excess of events is 2.7 standard deviations. This measurement corresponds to $\sigma($pp $\to$ H$)\mathcal{B}($H $\to$ Z$\gamma)$ = 0.21 $\pm$ 0.08 pb. The observed (expected) upper limit at 95% confidence level on $\mu$ is 4.1 (1.8). The ratio of branching fractions $\mathcal{B}($H $\to$ Z$\gamma) / \mathcal{B}($H $\to$ $\gamma\gamma)$ is measured to be 1.5 $^{+0.7}_{-0.6}$, which agrees with the standard model prediction of 0.69 $\pm$ 0.04 at the 1.5 standard deviation level.

A search is presented for long-lived charged particles that decay within the volume of the silicon tracker of the CMS experiment. Such particles can produce events with an isolated track that is missing hits in the outermost layers of the silicon tracker, and is also associated with little energy deposited in the calorimeters and no hits in the muon detectors. The search for events with this "disappearing track" signature is performed in a sample of proton-proton collisions recorded by the CMS experiment at the LHC with a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 101 fb$^{-1}$ recorded in 2017 and 2018. The observation of 48 events is consistent with the estimated background of 47.8 $_{-2.3}^{+2.7}$ (stat) $\pm$ 8.1 (syst) events. Upper limits are set on chargino production in the context of an anomaly-mediated supersymmetry breaking model for purely wino and higgsino neutralino scenarios. At 95% confidence level, the first constraint is placed on chargino masses in the higgsino case, excluding below 750 (175) GeV for a lifetime of 3 (0.05) ns. In the wino case, the results of this search are combined with a previous CMS search to produce a result representing the complete LHC data set recorded in 2015-2018, the most stringent constraints to date. At 95% confidence level, chargino masses in the wino case are excluded below 884 (474) GeV for a lifetime of 3 (0.2) ns.

A search for physics beyond the standard model (SM) in final states with an electron or muon and missing transverse momentum is presented. The analysis uses data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected with the CMS detector at the LHC in 2016–2018 and corresponding to an integrated luminosity of 138 fb−1. No significant deviation from the SM prediction is observed. Model-independent limits are set on the production cross section of W’ bosons decaying into lepton-plus-neutrino final states. Within the framework of the sequential standard model, with the combined results from the electron and muon decay channels a W’ boson with mass less than 5.7 TeV is excluded at 95% confidence level. Results on a SM precision test, the determination of the oblique electroweak W parameter, are presented using LHC data for the first time. These results together with those from the direct W’ resonance search are used to extend existing constraints on composite Higgs scenarios. This is the first experimental exclusion on compositeness parameters using results from LHC data other than Higgs boson measurements.

A search for physics beyond the standard model (SM) in the final state with a hadronically decaying tau lepton and a neutrino is presented. This analysis is based on data recorded by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC, corresponding to a total integrated luminosity of 138 fb$^{=1}$. The transverse mass spectrum is analyzed for the presence of new physics. No significant deviation from the SM prediction is observed. Limits are set on the production cross section of a W' boson decaying into a tau lepton and a neutrino. Lower limits are set on the mass of the sequential SM-like heavy charged vector boson and the mass of a quantum black hole. Upper limits are placed on the couplings of a new boson to the SM fermions. Constraints are put on a nonuniversal gauge interaction model and an effective field theory model. For the first time, upper limits on the cross section of $t$-channel leptoquark (LQ) exchange are presented. These limits are translated into exclusion limits on the LQ mass and on its coupling in the $t$-channel. The sensitivity of this analysis extends into the parameter space of LQ models that attempt to explain the anomalies observed in B meson decays. The limits presented for the various interpretations are the most stringent to date. Additionally, a model-independent limit is provided.

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.

Double-parton scattering is investigated using events with a Z boson and jets. The Z boson is reconstructed using only the dimuon channel. The measurements are performed with proton-proton collision data recorded by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected in the year 2016. Differential cross sections of Z + $\geq$ 1 jet and Z + $\geq$ 2 jets are measured with transverse momentum of the jets above 20 GeV and pseudorapidity $|\eta|$$\lt$ 2.4. Several distributions with sensitivity to double-parton scattering effects are measured as functions of the angle and the transverse momentum imbalance between the Z boson and the jets. The measured distributions are compared with predictions from several event generators with different hadronization models and different parameter settings for multiparton interactions. The measured distributions show a dependence on the hadronization and multiparton interaction simulation parameters, and are important input for future improvements of the simulations.

Three searches are presented for signatures of physics beyond the standard model (SM) in $\tau\tau$ final states in proton-proton collisions at the LHC, using a data sample collected with the CMS detector at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Upper limits at 95% confidence level (CL) are set on the products of the branching fraction for the decay into $\tau$ leptons and the cross sections for the production of a new boson $\phi$, in addition to the H(125) boson, via gluon fusion (gg$\phi$) or in association with b quarks, ranging from $\mathcal{O}$(10 pb) for a mass of 60 GeV to 0.3 fb for a mass of 3.5 TeV each. The data reveal two excesses for gg$\phi$ production with local $p$-values equivalent to about three standard deviations at $m_\phi$ = 0.1 and 1.2 TeV. In a search for $t$-channel exchange of a vector leptoquark U$_1$, 95% CL upper limits are set on the dimensionless U$_1$ leptoquark coupling to quarks and $\tau$ leptons ranging from 1 for a mass of 1 TeV to 6 for a mass of 5 TeV, depending on the scenario. In the interpretations of the $M_\mathrm{h}^{125}$ and $M_\mathrm{h, EFT}^{125}$ minimal supersymmetric SM benchmark scenarios, additional Higgs bosons with masses below 350 GeV are excluded at 95% CL.

Results are presented on a search for CP violation in the production and decay of top quark-antiquark pairs in the lepton+jets channel. The search is based on data from proton-proton collisions at 13 TeV, collected with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. Possible CP violation effects are evaluated by measuring uncorrected asymmetries in observables constructed from linearly independent four-momentum vectors of the final-state particles. The dimensionless chromoelectric dipole moment of the top quark obtained from the observed asymmetries is measured to be 0.04 $\pm$ 0.10 (stat) $\pm$ 0.07 (syst), and the asymmetries exhibit no evidence for CP-violating effects, consistent with expectations from the standard model.

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

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.

The first search for scalar leptoquarks produced in $\tau$-lepton-quark collisions is presented. It is based on a set of proton-proton collision data recorded with the CMS detector at the LHC at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The reconstructed final state consists of a jet, significant missing transverse momentum, and a $\tau$ lepton reconstructed through its hadronic or leptonic decays. Limits are set on the product of the leptoquark production cross section and branching fraction and interpreted as exclusions in the plane of the leptoquark mass and the leptoquark-$\tau$-quark coupling strength.

A search is presented for fractionally charged particles with charge below 1$e$, using their small energy loss in the tracking detector as a key variable to observe a signal. The analyzed data set corresponds to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions collected at $\sqrt{s}$ = 13 TeV in 2016-2018 at the CERN LHC. This is the first search at the LHC for new particles with charges between $e/$3 and $e$. Masses up to 640 GeV and charges as low as $e/$3 are excluded at 95% confidence level. These are the most stringent limits to date for the considered Drell-Yan-like production mode.

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 search for low-mass dilepton resonances in Higgs boson decays is conducted in the four-lepton final state. The decay is assumed to proceed via a pair of beyond the standard model particles, or one such particle and a Z boson. The search uses proton-proton collision data collected with the CMS detector at the CERN LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$, at a center-of-mass energy $\sqrt{s} =$ 13 TeV. No significant deviation from the standard model expectation is observed. Upper limits at 95% confidence level are set on model-independent Higgs boson decay branching fractions. Additionally, limits on dark photon and axion-like particle production, based on two specific models, are reported.

A search is reported for pairs of light Higgs bosons (H$_1$) produced in supersymmetric cascade decays in final states with small missing transverse momentum. A data set of LHC pp collisions collected with the CMS detector at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 138 fb$^{-1}$ is used. The search targets events where both H$_1$ bosons decay into $\mathrm{b\bar{b}}$ pairs that are reconstructed as large-radius jets using substructure techniques. No evidence is found for an excess of events beyond the background expectations of the standard model (SM). Results from the search are interpreted in the next-to-minimal supersymmetric extension of the SM, where a "singlino" of small mass leads to squark and gluino cascade decays that can predominantly end in a highly Lorentz-boosted singlet-like H$_1$ and a singlino-like neutralino of small transverse momentum. Upper limits are set on the product of the squark or gluino pair production cross section and the square of the $\mathrm{b\bar{b}}$ branching fraction of the H$_1$ in a benchmark model containing almost mass-degenerate gluinos and light-flavour squarks. Under the assumption of an SM-like H$_1$$\to$$\mathrm{b\bar{b}}$ branching fraction, H$_1$ bosons with masses in the range 40-120 GeV arising from the decays of squarks or gluinos with a mass of 1200 to 2500 GeV are excluded at 95% confidence level.

The first measurements of diboson production cross sections in proton-proton interactions at a center-of-mass energy of 5.02 TeV are reported. They are based on data collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 302 pb$^{-1}$. Events with two, three, or four charged light leptons (electrons or muons) in the final state are analyzed. The WW, WZ, and ZZ total cross sections are measured as $\sigma_\mathrm{WW} =$ 37.0 $^{+5.5}_{-5.2}$ (stat) $^{+2.7}_{-2.6}$ (syst) pb, $\sigma_\mathrm{WZ} =$ 6.4 $^{+2.5}_{-2.1}$ (stat) $^{+0.5}_{-0.3}$ (syst) pb, and $\sigma_\mathrm{ZZ} =$ 5.3 $^{+2.5}_{-2.1}$ (stat) $^{+0.5}_{-0.4}$ (syst) pb. All measurements are in good agreement with theoretical calculations at combined next-to-next-to-leading order quantum chromodynamics and next-to-leading order electroweak accuracy.

A search is performed for exclusive high-mass $\gamma\gamma$$\to$ WW and $\gamma\gamma$$\to$ ZZ production in proton-proton collisions using intact forward protons reconstructed in near-beam detectors, with both weak bosons decaying into boosted and merged jets. The analysis is based on a sample of proton-proton collisions collected by the CMS and TOTEM experiments at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 100 fb$^{-1}$. No excess above the standard model background prediction is observed, and upper limits are set on the pp $\to$ pWWp and pp $\to$ pZZp cross sections in a fiducial region defined by the diboson invariant mass $m$(VV) $\lt$ 1 TeV (with V = W, Z) and proton fractional momentum loss 0.04 $\lt$$\xi$$\lt$ 0.20. The results are interpreted as new limits on dimension-6 and dimension-8 anomalous quartic gauge couplings.

The measurements of the inclusive and differential fiducial cross sections of the Higgs boson decaying to a pair of photons are presented. The analysis is performed using proton-proton collisions data recorded with the CMS detector at the LHC at a centre-of-mass energy of 13 TeV and corresponding to an integrated luminosity of 137 fb$^{-1}$. The inclusive fiducial cross section is measured to be $\sigma_\mathrm{fid}$ = 73.4 $_{-5.3}^{+5.4}$ (stat) ${}_{-2.2}^{+2.4}$ (syst) fb, in agreement with the standard model expectation of 75.4 $\pm$ 4.1 fb. The measurements are also performed in fiducial regions targeting different production modes and as function of several observables describing the diphoton system, the number of additional jets present in the event, and other kinematic observables. Two double differential measurements are performed. No significant deviations from the standard model expectations are observed.

A measurement is presented for the electroweak production of a W boson, a photon ($\gamma$), and two jets (j) in proton-proton collisions. The leptonic decay of the W boson is selected by requiring one identified electron or muon and large missing transverse momentum. The two jets are required to have large invariant dijet mass and large separation in pseudorapidity. The measurement is performed with the data collected by the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross section for the electroweak W$\gamma$jj production is 23.5 $^{+4.9}_{-4.7}$ fb, whereas the total cross section for W$\gamma$jj production is 113 $\pm$ 13 fb. Differential cross sections are also measured with the distributions unfolded to the particle level. All results are in agreement with the standard model expectations. Constraints are placed on anomalous quartic gauge couplings (aQGCs) in terms of dimension-8 effective field theory operators. These are the most stringent limits to date on the aQGCs parameters $f_\mathrm{M,2-5}$$/$$\Lambda^4$ and $f_\mathrm{T,6-7}$$/$$\Lambda^4$.

Production of prompt J/$\psi$ meson pairs in proton-proton collisions at $\sqrt{s}$ = 7 TeV is measured with the CMS experiment at the LHC in a data sample corresponding to an integrated luminosity of about 4.7 inverse-femtobarns. The two J/$\psi$ mesons are fully reconstructed via their decays into $\mu^+\mu^-$ pairs. This observation provides for the first time access to the high-transverse-momentum region of J/$\psi$ pair production where model predictions are not yet established. The total and differential cross sections are measured in a phase space defined by the individual J/$\psi$ transverse momentum ($p_T^{J/\psi}$) and rapidity (|$y^{J/\psi}$|): |$y^{J/\psi}$| lower than 1.2 for $p_T^{J/\psi}$ greater than 6.5 GeV/c; |$y^{J/\psi}$| in [1.2,1.43] for a $p_T$ threshold that scales linearly with |$y^{J/\psi}$| from 6.5 to 4.5 GeV/c; and |$y^{J/\psi}$| in [1.43,2.2] for $p_T^{J/\psi}$ greater than 4.5 GeV/c. The total cross section, assuming unpolarized prompt J/$\psi$ pair production is 1.49 $\pm$ 0.07 (stat.) $\pm$ 0.13 (syst.) nb. Different assumptions about the J/$\psi$ polarization imply modifications to the cross section ranging from -31% to +27%.

A study of the associated production of a Z boson and a charm quark jet (Z + c), and a comparison to production with a b quark jet (Z + b), in pp collisions at a centre-of-mass energy of 8 TeV are presented. The analysis uses a data sample corresponding to an integrated luminosity of 19.7 fb$^{-1}$, collected with the CMS detector at the CERN LHC. The Z boson candidates are identified through their decays into pairs of electrons or muons. Jets originating from heavy flavour quarks are identified using semileptonic decays of c or b flavoured hadrons and hadronic decays of charm hadrons. The measurements are performed in the kinematic region with two leptons with $p_{\rm T}^{\ell} > $ 20 GeV, ${|\eta^{\ell}|} < $ 2.1, 71 $ < m_{\ell\ell} < $ 111 GeV, and heavy flavour jets with $p_{\rm T}^{{\rm jet}} > $ 25 GeV and ${|\eta^{{\rm jet}}|} < $ 2.5. The Z + c production cross section is measured to be $\sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{c} + \mathrm{X}) {\cal B}(\mathrm{Z} \rightarrow \ell^+\ell^-) = $ 8.8 $ \pm $ 0.5 (stat) $ \pm $ 0.6 (syst) pb. The ratio of the Z + c and Z + b production cross sections is measured to be $\sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{c} + \mathrm{X}) / \sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{b} + \mathrm{X}) = $ 2.0 $ \pm $ 0.2 (stat) $ \pm $ 0.2 (syst). The Z + c production cross section and the cross section ratio are also measured as a function of the transverse momentum of the Z boson and of the heavy flavour jet. The measurements are compared with theoretical predictions.

The observation of the production of four top quarks in proton-proton collisions is reported, based on a data sample collected by the CMS experiment at a center-of-mass energy of 13 TeV in 2016-2018 at the CERN LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with two same-sign, three, or four charged leptons (electrons and muons) and additional jets are analyzed. Compared to previous results in these channels, updated identification techniques for charged leptons and jets originating from the hadronization of b quarks, as well as a revised multivariate analysis strategy to distinguish the signal process from the main backgrounds, lead to an improved expected signal significance of 4.9 standard deviations above the background-only hypothesis. Four top quark production is observed with a significance of 5.6 standard deviations, and its cross section is measured to be 17.7 $^{+3.7}_{-3.5}$ (stat) $^{+2.3}_{-1.9}$ (syst) fb, in agreement with the available standard model predictions.

The transverse momentum ($p_\mathrm{T}$) distributions of $\Lambda$, $\Xi^-$, and $\Omega^-$ baryons, their antiparticles, and K$^0_\mathrm{S}$ mesons are measured in proton-proton (pp) and proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV over a broad rapidity range. The data, corresponding to integrated luminosities of 40.2 nb$^{-1}$ and 15.6 $\mu$b$^{-1}$ for pp and pPb collisions, respectively, were collected by the CMS experiment. The nuclear modification factor $R_\mathrm{pPb}$, defined as the ratio of the particle yield in pPb collisions and a scaled pp reference, is measured for each particle. A strong dependence on particle species is observed in the $p_\mathrm{T}$ range from 2 to 7 GeV, where $R_\mathrm{pPb}$ for K$^0_\mathrm{S}$ is consistent with unity, while an enhancement ordered by strangeness content and/or particle mass is observed for the three baryons. In pPb collisions, the strange hadron production is asymmetric about the nucleon-nucleon center-of-mass rapidity. Enhancements, which depend on the particle type, are observed in the direction of the Pb beam. The results are compared to predictions from EPOS LHC, which includes parametrized radial flow. The model is in qualitative agreement with the $R_\mathrm{pPb}$ data, but fails to describe the dependence on particle species in the yield asymmetries measured away from mid-rapidity in pPb collisions.

The production of four top quarks ($\mathrm{t\bar{t}t\bar{t}}$) is studied with LHC proton-proton collision data samples collected by the CMS experiment at a center-of-mass energy of 13 TeV, and corresponding to integrated luminosities of up to 138 fb$^{-1}$. Events that have no leptons (all-hadronic), one lepton, or two opposite-sign leptons (where lepton refers only to prompt electrons or prompt muons) are considered. This is the first $\mathrm{t\bar{t}t\bar{t}}$ measurement that includes the all-hadronic final state. The observed significance of the $\mathrm{t\bar{t}t\bar{t}}$ signal in these final states of 3.9 standard deviations (1.5 expected) provides evidence for $\mathrm{t\bar{t}t\bar{t}}$ production, with a measured cross section of 36 $^{+12}_{-11}$ fb. Combined with earlier CMS results in other final states, the signal significance is 4.0 standard deviations (3.2 expected). The combination returns an observed cross section of 17 $\pm$ 4 (stat) $\pm$ 3 (syst) fb, which is consistent with the standard model prediction.

Modifications to the distribution of charged particles with respect to high transverse momentum ($p_\mathrm{T}$) jets passing through a quark-gluon plasma are explored using the CMS detector. Back-to-back dijets are analyzed in lead-lead and proton-proton collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles in bins of relative pseudorapidity and angular distance from the leading and subleading jet axes. In comparing the lead-lead and proton-proton collision results, modifications to the charged-particle relative distance distribution and to the momentum distributions around the jet axis are found to depend on the dijet momentum balance $x_j$, which is the ratio between the subleading and leading jet $p_\mathrm{T}$. For events with $x_j$$\approx$ 1, these modifications are observed for both the leading and subleading jets. However, while subleading jets show significant modifications for events with a larger dijet momentum imbalance, much smaller modifications are found for the leading jets in these events.

For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV, almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of two below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC, QGSJET II, and HIJING predictions are all at least a factor of five lower than the data. The latter effect might be explained by a significant contribution of ultra-peripheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.

A search for long-lived particles (LLPs) produced in decays of standard model (SM) Higgs bosons is presented. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, recorded at the LHC in 2016-2018. A novel technique is employed to reconstruct decays of LLPs in the endcap muon detectors. The search is sensitive to a broad range of LLP decay modes and to masses as low as a few GeV. No excess of events above the SM background is observed. The most stringent limits to date on the branching fraction of the Higgs boson to LLPs subsequently decaying to quarks and $\tau^+\tau^-$ are found for proper decay lengths greater than 6, 20, and 40 m, for LLP masses of 7, 15, and 40 GeV, respectively.

A search for a charged Higgs boson H$^\pm$ decaying into a heavy neutral Higgs boson H and a W boson is presented. The analysis targets the H decay into a pair of tau leptons with at least one of them decaying hadronically and with an additional electron or muon present in the event. The search is based on proton-proton collision data recorded by the CMS experiment during 2016-2018 at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The data are consistent with standard model background expectations. Upper limits at 95% confidence level are set on the product of the cross section and branching fraction for an H$^\pm$ in the mass range of 300-700 GeV, assuming an H with a mass of 200 GeV. The observed limits range from 0.085 pb for an H$^\pm$ mass of 300 GeV to 0.019 pb for a mass of 700 GeV. These are the first limits on H$^\pm$ production in the H$^\pm$ $\to$ HW$^\pm$ decay channel at the LHC.

A search for the electroweak production of a vector-like quark T$'$, decaying to a top quark and a Higgs boson is presented. The search is based on a sample of proton-proton collision events recorded at the LHC at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. This is the first T$'$ search that exploits the Higgs boson decay to a pair of photons. For narrow isospin singlet T$'$ states with masses up to 1.1 TeV, the excellent diphoton invariant mass resolution of 1-2% results in an increased sensitivity compared to previous searches based on the same production mechanism. The electroweak production of a T$'$ quark with mass up to 960 GeV is excluded at 95% confidence level, assuming a coupling strength $\kappa_\mathrm{T}$ = 0.25 and a relative decay width $\Gamma/M_{\mathrm{T}'}$ $\lt$ 5%.

Measurements of the inclusive and normalised differential cross sections are presented for the production of single top quarks in association with a W boson in proton-proton collisions at a centre-of-mass energy of 13 TeV. The data used were recorded with the CMS detector at the LHC during 2016-2018, and correspond to an integrated luminosity of 138 fb$^{-1}$. Events containing one electron and one muon in the final state are analysed. For the inclusive measurement, a multivariate discriminant, exploiting the kinematic properties of the events is used to separate the signal from the dominant $\mathrm{t\bar{t}}$ background. A cross section of 79.2 $\pm$ 0.9 (stat) $^{+7.7}_{-8.0}$ (syst) $\pm$ 1.2 (lumi) pb is obtained, consistent with the predictions of the standard model. For the differential measurements, a fiducial region is defined according to the detector acceptance, and the requirement of exactly one jet coming from the fragmentation of a bottom quark. The resulting distributions are unfolded to particle level and agree with the predictions at next-to-leading order in perturbative quantum chromodynamics.

A measurement of the top quark mass is performed using a data sample enriched with single top quark events produced in the $t$ channel. The study is based on proton-proton collision data, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC in 2016. Candidate events are selected by requiring an isolated high-momentum lepton (muon or electron) and exactly two jets, of which one is identified as originating from a bottom quark. Multivariate discriminants are designed to separate the signal from the background. Optimized thresholds are placed on the discriminant outputs to obtain an event sample with high signal purity. The top quark mass is found to be 172.13 $^{+0.76}_{-0.77}$ GeV, where the uncertainty includes both the statistical and systematic components, reaching sub-GeV precision for the first time in this event topology. The masses of the top quark and antiquark are also determined separately using the lepton charge in the final state, from which the mass ratio and difference are determined to be 0.9952 $^{+0.0079}_{-0.0104}$ and 0.83 $^{+1.79}_{-1.35}$ GeV, respectively. The results are consistent with $CPT$ invariance.

A search is reported for heavy resonances and quantum black holes decaying into e$\mu$, e$\tau$, and $\mu\tau$ final states in proton-proton collision data recorded by the CMS experiment at the CERN LHC during 2016-2018 at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The e$\mu$, e$\tau$, and $\mu\tau$ invariant mass spectra are reconstructed, and no evidence is found for physics beyond the standard model. Upper limits are set at 95% confidence level on the product of the cross section and branching fraction for lepton flavor violating signals. Three benchmark signals are studied: resonant $\tau$ sneutrino production in $R$ parity violating supersymmetric models, heavy Z' gauge bosons with lepton flavor violating decays, and nonresonant quantum black hole production in models with extra spatial dimensions. Resonant $\tau$ sneutrinos are excluded for masses up to 4.2 TeV in the e$\mu$ channel, 3.7 TeV in the e$\tau$ channel, and 3.6 TeV in the $\mu\tau$ channel. A Z' boson with lepton flavor violating couplings is excluded up to a mass of 5.0 TeV in the e$\mu$ channel, up to 4.3 TeV in the e$\tau$ channel, and up to 4.1 TeV in the $\mu\tau$ channel. Quantum black holes in the benchmark model are excluded up to the threshold mass of 5.6 TeV in the e$\mu$ channel, 5.2 TeV in the e$\tau$ channel, and 5.0 TeV in the $\mu\tau$ channel. In addition, model-independent limits are extracted to allow comparisons with other models for the same final states and similar event selection requirements. The results of these searches provide the most stringent limits available from collider experiments for heavy particles that undergo lepton flavor violating decays.

Protons consist of three valence quarks, two up-quarks and one down-quark, held together by gluons and a sea of quark-antiquark pairs. Collectively, quarks and gluons are referred to as partons. In a proton-proton collision, typically only one parton of each proton undergoes a hard scattering - referred to as single-parton scattering - leaving the remainder of each proton only slightly disturbed. Here, we report the study of double- and triple-parton scatterings through the simultaneous production of three J/$\psi$ mesons, which consist of a charm quark-antiquark pair, in proton-proton collisions recorded with the CMS experiment at the Large Hadron Collider. We observed this process - reconstructed through the decays of J/$\psi$ mesons into pairs of oppositely charged muons - with a statistical significance above five standard deviations. We measured the inclusive fiducial cross section to be 272 $^{+141}_{-104}$ (stat) $\pm$ 17 (syst) fb, and compared it to theoretical expectations for triple-J/$\psi$ meson production in single-, double- and triple-parton scattering scenarios. Assuming factorization of multiple hard-scattering probabilities in terms of single-parton scattering cross sections, double- and triple-parton scattering are the dominant contributions for the measured process.

Measurements of jet substructure describing the composition of quark- and gluon-initiated jets are presented. Proton-proton (pp) collision data at $\sqrt{s}$ =13 TeV collected with the CMS detector are used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Generalized angularities are measured that characterize the jet substructure and distinguish quark- and gluon-initiated jets. These observables are sensitive to the distributions of transverse momenta and angular distances within a jet. The analysis is performed using a data sample of dijet events enriched in gluon-initiated jets, and, for the first time, a Z+jet event sample enriched in quark-initiated jets. The observables are measured in bins of jet transverse momentum, and as a function of the jet radius parameter. Each measurement is repeated applying a "soft drop" grooming procedure that removes soft and large angle radiation from the jet. Using these measurements, the ability of various models to describe jet substructure is assessed, showing a clear need for improvements in Monte Carlo generators.

A precise measurement of the invisible width of the Z boson produced in proton-proton collisions at a center-of-mass energy of 13 TeV is presented using data recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The result is obtained from a simultaneous fit to kinematic distributions for two data samples of Z boson plus jets: one dominated by Z boson decays to invisible particles and the other by Z boson decays to muon and electron pairs. The invisible width is measured to be 523 $\pm$ 3 (stat) $\pm$ 16 (syst) MeV. This result is the first precise measurement of the invisible width of the Z boson at a hadron collider, and is the single most precise direct measurement to date, competitive with the combined result of the direct measurements from the LEP experiments.

The first collider search for dark matter arising from a strongly coupled hidden sector is presented and uses a data sample corresponding to 138 fb$^{-1}$, collected with the CMS detector at the CERN LHC, at $\sqrt{s} =$ 13 TeV. The hidden sector is hypothesized to couple to the standard model (SM) via a heavy leptophobic Z' mediator produced as a resonance in proton-proton collisions. The mediator decay results in two "semivisible" jets, containing both visible matter and invisible dark matter. The final state therefore includes moderate missing energy aligned with one of the jets, a signature ignored by most dark matter searches. No structure in the dijet transverse mass spectra compatible with the signal is observed. Assuming the Z' has a universal coupling of 0.25 to the SM quarks, an inclusive search, relevant to any model that exhibits this kinematic behavior, excludes mediator masses of 1.5-4.0 TeV at 95% confidence level, depending on the other signal model parameters. To enhance the sensitivity of the search for this particular class of hidden sector models, a boosted decision tree (BDT) is trained using jet substructure variables to distinguish between semivisible jets and SM jets from background processes. When the BDT is employed to identify each jet in the dijet system as semivisible, the mediator mass exclusion increases to 5.1 TeV, for wider ranges of the other signal model parameters. These limits exclude a wide range of strongly coupled hidden sector models for the first time.

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.