Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

A search for supersymmetry is presented based on proton-proton collision events containing identified hadronically decaying top quarks, no leptons, and an imbalance $p_\mathrm{T}^\text{miss}$ in transverse momentum. The data were collected with the CMS detector at the CERN LHC at a center-of-mass energy of 13 TeV, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. Search regions are defined in terms of the multiplicity of bottom quark jet and top quark candidates, the $p_\mathrm{T}^\text{miss}$, the scalar sum of jet transverse momenta, and the $m_{\mathrm{T2}}$ mass variable. No statistically significant excess of events is observed relative to the expectation from the standard model. Lower limits on the masses of supersymmetric particles are determined at 95% confidence level in the context of simplified models with top quark production. For a model with direct top squark pair production followed by the decay of each top squark to a top quark and a neutralino, top squark masses up to 1020 GeV and neutralino masses up to 430 GeV are excluded. For a model with pair production of gluinos followed by the decay of each gluino to a top quark-antiquark pair and a neutralino, gluino masses up to 2040 GeV and neutralino masses up to 1150 GeV are excluded. These limits extend previous results.

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 WZ production cross section in proton-proton collisions at sqrt(s) = 13 TeV is measured with the CMS experiment at the LHC using a data sample corresponding to an integrated luminosity of 2.3 inverse femtobarns. The measurement is performed in the leptonic decay modes WZ to l nu l' l', where l, l'= e, mu. The measured cross section for the range 60 < m[l'l'] < 120 GeV is sigma(pp to WZ) = 39.9 +/- 3.2 (stat) +2.9/-3.1 (syst) +/- 0.4 (theo) +/- 1.3 (lumi) pb, consistent with the standard model prediction.

The azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV pPb data are extracted via long-range two-particle correlations as a function of event multiplicity and compared to corresponding results in pp and PbPb collisions. Using a four-particle cumulant technique, $v_n$ correlations are measured for the first time in pp and pPb collisions. The $v_2$ and $v_4$ coefficients are found to be positively correlated in all collision systems. For high multiplicity pPb collisions an anticorrelation of $v_2$ and $v_3$ is observed, with a similar correlation strength as in PbPb data at the same multiplicity. The new correlation results strengthen the case for a common origin of the collectivity seen in pPb and PbPb collisions in the measured multiplicity range.

A search for supersymmetry is presented based on multijet events with large missing transverse momentum produced in proton-proton collisions at a center-of-mass energy of sqrt(s) = 13 TeV. The data, corresponding to an integrated luminosity of 35.9 inverse femtobarns, were collected with the CMS detector at the CERN LHC in 2016. The analysis utilizes four-dimensional exclusive search regions defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No evidence for a significant excess of events is observed relative to the expectation from the standard model. Limits on the cross sections for the pair production of gluinos and squarks are derived in the context of simplified models. Assuming the lightest supersymmetric particle to be a weakly interacting neutralino, 95% confidence level lower limits on the gluino mass as large as 1800 to 1960 GeV are derived, and on the squark mass as large as 960 to 1390 GeV, depending on the production and decay scenario.

Transverse momentum spectra of charged pions, kaons, and protons are measured in proton-proton collisions at sqrt(s) = 13 TeV with the CMS detector at the LHC. The particles, identified via their energy loss in the silicon tracker, are measured in the transverse momentum range of pt ~ 0.1-1.7 GeV and rapidities abs(y) < 1. The pt spectra and integrated yields are compared to previous results at smaller sqrt(s) and to predictions of Monte Carlo event generators. The average pt increases with particle mass and charged particle multiplicity of the event. Comparisons with previous CMS results at sqrt(s) = 0.9, 2.76, and 7 TeV show that the average pt and the ratios of hadron yields feature very similar dependences on the particle multiplicity in the event, independently of the center-of-mass energy of the pp collision.

The top quark pair production cross section is measured for the first time in proton-proton collisions at sqrt(s) = 13 TeV by the CMS experiment at the CERN LHC, using data corresponding to an integrated luminosity of 43 inverse picobarns. The measurement is performed by analyzing events with at least one electron and one muon of opposite charge, and at least two jets. The measured cross section is 746 +/- 58 (stat) +/- 53 (syst) +/- 36 (lumi) pb, in agreement with the expectation from the standard model.

A measurement of the double-differential inclusive jet cross section as a function of jet transverse momentum pT and absolute jet rapidity |y| is presented. The analysis is based on proton-proton collisions collected by the CMS experiment at the LHC at a centre-of-mass energy of 13 TeV. The data samples correspond to integrated luminosities of 71 and 44 inverse picobarns for |y| < 3 and 3.2 < |y| < 4.7, respectively. Jets are reconstructed with the anti-kt clustering algorithm for two jet sizes, R, of 0.7 and 0.4, in a phase space region covering jet pT up to 2 TeV and jet rapidity up to |y| = 4.7. Predictions of perturbative quantum chromodynamics at next-to-leading order precision, complemented with electroweak and nonperturbative corrections, are used to compute the absolute scale and the shape of the inclusive jet cross section. The cross section difference in R, when going to a smaller jet size of 0.4, is best described by Monte Carlo event generators with next-to-leading order predictions matched to parton showering, hadronisation, and multiparton interactions. In the phase space accessible with the new data, this measurement provides a first indication that jet physics is as well understood at sqrt(s) = 13 TeV as at smaller centre-of-mass energies.

A search for top squark pair production in pp collisions at sqrt(s) = 13 TeV is performed using events with a single isolated electron or muon, jets, and a large transverse momentum imbalance. The results are based on data collected in 2016 with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 inverse femtobarns. No significant excess of events is observed above the expectation from standard model processes. Exclusion limits are set in the context of supersymmetric models of pair production of top squarks that decay either to a top quark and a neutralino or to a bottom quark and a chargino. Depending on the details of the model, we exclude top squarks with masses as high as 1120 GeV. Detailed information is also provided to facilitate theoretical interpretations in other scenarios of physics beyond the standard model.

Results on two-particle angular correlations for charged particles produced in pp collisions at a center-of-mass energy of 13 TeV are presented. The data were taken with the CMS detector at the LHC and correspond to an integrated luminosity of about 270 inverse nanobarns. The correlations are studied over a broad range of pseudorapidity (abs(eta) < 2.4) and over the full azimuth (phi) as a function of charged particle multiplicity and transverse momentum (pt). In high-multiplicity events, a long-range (abs(Delta eta) > 2.0), near-side (Delta phi approximately 0) structure emerges in the two-particle Delta eta-Delta phi correlation functions. The magnitude of the correlation exhibits a pronounced maximum in the range 1.0 < pt < 2.0 GeV/c and an approximately linear increase with the charged particle multiplicity, with an overall correlation strength similar to that found in earlier pp data at sqrt(s) = 7 TeV. The present measurement extends the study of near-side long-range correlations up to charged particle multiplicities of N[ch] approximately 180, a region so far unexplored in pp collisions. The observed long-range correlations are compared to those seen in pp, pPb, and PbPb collisions at lower collision energies.

Searches for invisible decays of the Higgs boson are presented. The data collected with the CMS detector at the LHC correspond to integrated luminosities of 5.1, 19.7, and 2.3 inverse femtobarns at centre-of-mass energies of 7, 8, and 13 TeV, respectively. The search channels target Higgs boson production via gluon fusion, vector boson fusion, and in association with a vector boson. Upper limits are placed on the branching fraction of the Higgs boson decay to invisible particles, as a function of the assumed production cross sections. The combination of all channels, assuming standard model production, yields an observed (expected) upper limit on the invisible branching fraction of 0.24 (0.23) at the 95% confidence level. The results are also interpreted in the context of Higgs-portal dark matter models.

The differential cross sections for inclusive production of B+ hadrons are measured as a function of the B+ transverse momentum pT[B] and rapidity y[B] in pp collisions at a centre-of-mass energy of 13 TeV, using data collected by the CMS experiment that correspond to an integrated luminosity of 49.4 inverse picobarns. The measurement uses the exclusive decay channel B+ to J/psi K+, with J/psi mesons that decay to a pair of muons. The results show a reasonable agreement with theoretical calculations within the uncertainties.

The differential cross section for inclusive particle production as a function of energy in proton-proton collisions at a center-of-mass energy of 13 TeV is measured in the very forward region of the CMS detector. The measurement is based on data collected with the CMS apparatus at the LHC, and corresponds to an integrated luminosity of 0.35 inverse microbarns. The energy is measured in the CASTOR calorimeter, which covers the pseudorapidity region -6.6 < eta < -5.2. The results are given as a function of the total energy deposited in CASTOR, as well as of its electromagnetic and hadronic components. The spectra are sensitive to the modeling of multiparton interactions in pp collisions, and provide new constraints for hadronic interaction models used in collider and in high energy cosmic ray physics.

This paper presents a measurement of the underlying event activity in proton-proton collisions at a center-of-mass energy of 13 TeV, performed using inclusive Z boson production events collected with the CMS experiment at the LHC. The analyzed data correspond to an integrated luminosity of 2.1 fb$^{-1}$. The underlying event activity is quantified in terms of the charged particle multiplicity, as well as of the scalar sum of the charged particles' transverse momenta in different topological regions defined with respect to the Z boson direction. The distributions are unfolded to the stable particle level and compared with predictions from various Monte Carlo event generators, as well as with similar CDF and CMS measurements at center-of-mass energies of 1.96 and 7 TeV respectively.

Normalized differential cross sections for top quark pair production are measured in the dilepton (e$^+$e$^-$, $\mu^+\mu^-$, and $\mu^\mp$e$^\pm$) decay channels in proton-proton collisions at a center-of-mass energy of 13 TeV. The measurements are performed with data corresponding to an integrated luminosity of 2.1 fb$^{-1}$ using the CMS detector at the LHC. The cross sections are measured differentially as a function of the kinematic properties of the leptons, jets from bottom quark hadronization, top quarks, and top quark pairs at the particle and parton levels. The results are compared to several Monte Carlo generators that implement calculations up to next-to-leading order in perturbative quantum chromodynamics interfaced with parton showering, and also to fixed-order theoretical calculations of top quark pair production up to next-to-next-to-leading order.

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 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.

A measurement is presented of the associated production of a single top quark and a Z boson. The study uses data from proton-proton collisions at $\sqrt{s} =$ 13 TeV recorded by the CMS experiment, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Using final states with three leptons (electrons or muons), the tZq production cross section is measured to be $\sigma$(pp$\to$tZq$\to$Wb$\ell^+\ell^-$q) = $123 ^{+33}_{-31}$ (stat) $^{+29}_{-23}$ (syst) fb, where $\ell$ stands for electrons, muons, or $\tau$ leptons, with observed and expected significances of 3.7 and 3.1 standard deviations, respectively.

Properties of the Higgs boson are measured in the H to ZZ to 4l (l= e, mu) decay channel. A data sample of proton-proton collisions at sqrt(s) = 13 TeV, collected with the CMS detector at the LHC and corresponding to an integrated luminosity of 35.9 inverse femtobarns is used. The signal strength modifier mu, defined as the ratio of the observed Higgs boson rate in the H to ZZ to 4l decay channel to the standard model expectation, is measured to be mu = 1.05 +0.19/-0.17 at m[H ]= 125.09 GeV, the combined ATLAS and CMS measurement of the Higgs boson mass. The signal strength modifiers for the individual Higgs boson production modes are also measured. The cross section in the fiducial phase space defined by the requirements on lepton kinematics and event topology is measured to be 2.92 +0.48/-0.44 (stat) +0.28/-0.24 (syst) fb, which is compatible with the standard model prediction of 2.76 +/- 0.14 fb. Differential cross sections are reported as a function of the transverse momentum of the Higgs boson, the number of associated jets, and the transverse momentum of the leading associated jet. The Higgs boson mass is measured to be m[H] = 125.26 +/- 0.21 GeV and the width is constrained using on-shell production to be Gamma[H] < 1.10 GeV, at 95% confidence level.

The cross section of top quark-antiquark pair production in proton-proton collisions at sqrt(s) = 13 TeV is measured by the CMS experiment at the LHC, using data corresponding to an integrated luminosity of 2.2 inverse femtobarns. The measurement is performed by analyzing events in which the final state includes one electron, one muon, and two or more jets, at least one of which is identified as originating from hadronization of a b quark. The measured cross section is 815 +/- 9 (stat) +/- 38 (syst) +/- 19 (lumi) pb, in agreement with the expectation from the standard model.

Differential and double-differential cross sections for the production of top quark pairs in proton-proton collisions at 13 TeV are measured as a function of jet multiplicity and of kinematic variables of the top quarks and the top quark-antiquark system. This analysis is based on data collected by the CMS experiment at the LHC corresponding to an integrated luminosity of 2.3 inverse femtobarns. The measurements are performed in the lepton+jets decay channels with a single muon or electron in the final state. The differential cross sections are presented at particle level, within a phase space close to the experimental acceptance, and at parton level in the full phase space. The results are compared to several standard model predictions.

Searches for resonant and nonresonant pair-produced Higgs bosons (HH) decaying respectively into ll nu nu, through either W or Z bosons, and bbbar are presented. The analyses are based on a sample of proton-proton collisions at sqrt(s) = 13 TeV, collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 inverse femtobarns. Data and predictions from the standard model are in agreement within uncertainties. For the standard model HH hypothesis, the data exclude at 95% confidence level a product of the production cross section and branching fraction larger than 72 fb, corresponding to 79 times the prediction, consistent with expectations. Constraints are placed on different scenarios considering anomalous couplings, which could affect the rate and kinematics of HH production. Upper limits at 95% confidence level are set on the production cross section of narrow-width spin-0 and spin-2 particles decaying to Higgs boson pairs, the latter produced with minimal gravity-like coupling.

Bose-Einstein correlations of charged hadrons are measured over a broad multiplicity range, from a few particles up to about 250 reconstructed charged hadrons in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The results are based on data collected using the CMS detector at the LHC during runs with a special low-pileup configuration. Three analysis techniques with different degrees of dependence on simulations are used to remove the non-Bose-Einstein background from the correlation functions. All three methods give consistent results. The measured lengths of homogeneity are studied as functions of particle multiplicity as well as average pair transverse momentum and mass. The results are compared with data from both CMS and ATLAS at $\sqrt{s} =$ 7 TeV, as well as with theoretical predictions.

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.

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 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.

Four-lepton production in proton-proton collisions, pp to (Z/gamma*)(Z/gamma*) to l+l-l'+l'-, where l, l' = e or mu, is studied at a center-of-mass energy of 13 TeV with the CMS detector at the LHC. The data sample corresponds to an integrated luminosity of 2.6 inverse femtobarns. The ZZ production cross section, sigma(pp to ZZ) = 14.6 +1.9/-1.8 (stat) +0.5/-0.3 (syst) +/- 0.2 (theo) +/- 0.4 (lumi) pb, is measured for events with two opposite-sign, same-flavor lepton pairs produced in the mass region 60 < m[l+l-], m[l'+l'-] < 120 GeV. The Z boson branching fraction to four leptons is measured to be B(Z to l+l-l'+l'-) = 4.9 +0.8/-0.7 (stat) +0.3/-0.2 (syst) +0.2/-0.1 (theo) +/- 0.1 (lumi) x E-6 for the four-lepton invariant mass in the range 80 < m[l+l-l'+l'-] < 100 GeV and dilepton mass m[l+l-] > 4 GeV for all opposite-sign, same-flavor lepton pairs. The results are in agreement with standard model predictions.

A measurement of the differential cross sections for a W boson produced in association with jets in the muon decay channel is presented. The measurement is based on 13 TeV proton-proton collision data corresponding to an integrated luminosity of 2.2 inverse femtobarns, recorded by the CMS detector at the LHC. The cross sections are reported as functions of jet multiplicity, jet transverse momentum pT, jet rapidity, the scalar pT sum of the jets, and angular correlations between the muon and the jet for different jet multiplicities. The measured cross sections are in agreement with predictions that include multileg leading-order (LO) and next-to-LO matrix element calculations interfaced with parton showers, as well as a next-to-next-to-LO calculation for the W boson and one jet production.

A search for lepton flavour violating decays of the Higgs boson in the $\mu\tau$ and e$\tau$ decay modes is presented. The search is based on a data set corresponding to an integrated luminosity of 35.9 fb$^{-1}$ of proton-proton collisions collected with the CMS detector in 2016, at a centre-of-mass energy of 13 TeV. No significant excess over the standard model expectation is observed. The observed (expected) upper limits on the lepton flavour violating branching fractions of the Higgs boson are $\mathcal{B}$(H$\to\mu\tau$) < 0.25% (0.25%) and $\mathcal{B}$(H$\to$e$\tau$) < 0.61% (0.37%), at 95% confidence level. These results are used to derive upper limits on the off-diagonal $\mu\tau$ and e$\tau$ Yukawa couplings $\sqrt{|{Y_{\mu\tau}}|^{2}+|{Y_{\tau\mu}}|^{2}} < 1.43 \times 10^{-3}$ and $\sqrt{|{Y_{\mathrm{e}\tau}}|^{2}+|{Y_{\tau\mathrm{e}}}|^{2}} < 2.26 \times 10^{-3}$ at 95% confidence level. The limits on the lepton flavour violating branching fractions of the Higgs boson and on the associated Yukawa couplings are the most stringent to date.

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.

Results are reported from a search for supersymmetric particles in proton-proton collisions in the final state with a single, high transverse momentum lepton; multiple jets, including at least one b-tagged jet; and large missing transverse momentum. The data sample corresponds to an integrated luminosity of 2.3 inverse femtobarns at sqrt(s) = 13 TeV, recorded by the CMS experiment at the LHC. The search focuses on processes leading to high jet multiplicities, such as gluino pair production with gluinos to t t-bar neutralino[1]. The quantity M[J], defined as the sum of the masses of the large-radius jets in the event, is used in conjunction with other kinematic variables to provide discrimination between signal and background and as a key part of the background estimation method. The observed event yields in the signal regions in data are consistent with those expected for standard model backgrounds, estimated from control regions in data. Exclusion limits are obtained for a simplified model corresponding to gluino pair production with three-body decays into top quarks and neutralinos. Gluinos with a mass below 1600 GeV are excluded at a 95% confidence level for scenarios with low neutralino[1] mass, and neutralinos with a mass below 800 GeV are excluded for a gluino mass of about 1300 GeV. For models with two-body gluino decays producing on-shell top squarks, the excluded region is only weakly sensitive to the top squark mass.

A search for physics beyond the standard model in final states with at least one photon, large transverse momentum imbalance, and large total transverse event activity is presented. Such topologies can be produced in gauge-mediated supersymmetry models in which pair-produced gluinos or squarks decay to photons and gravitinos via short-lived neutralinos. The data sample corresponds to an integrated luminosity of 35.9 inverse femtobarns of proton-proton collisions at sqrt(s) = 13 TeV recorded by the CMS experiment at the LHC in 2016. No significant excess of events above the expected standard model background is observed. The data are interpreted in simplified models of gluino and squark pair production, in which gluinos or squarks decay via neutralinos to photons. Gluino masses of up to 1.50-2.00 TeV and squark masses up to 1.30-1.65 TeV are excluded at 95% confidence level, depending on the neutralino mass and branching fraction.

A search for a heavy resonance decaying into a top quark and antiquark ($\mathrm{t\bar{t}}$) pair is performed using proton-proton collisions at $\sqrt{s} =$ 13 TeV. The search uses the data set collected with the CMS detector in 2016, which corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The analysis considers three exclusive final states and uses reconstruction techniques that are optimized for top quarks with high Lorentz boosts, which requires the use of nonisolated leptons and jet substructure techniques. No significant excess of events relative to the expected yield from standard model processes is observed. Upper limits on the production cross section of heavy resonances decaying to a $\mathrm{t\bar{t}}$ pair are calculated. Limits are derived for a leptophobic topcolor Z' resonance with widths of 1, 10, and 30%, relative to the mass of the resonance, and exclude masses up to 3.80, 5.25, and 6.65 TeV, respectively. Kaluza-Klein excitations of the gluon in the Randall-Sundrum model are excluded up to 4.55 TeV. To date, these are the most stringent limits on $\mathrm{t\bar{t}}$ resonances.

A measurement of the cross section of the associated production of a single top quark and a W boson in final states with a muon or electron and jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data correspond to an integrated luminosity of 36 fb$^{-1}$ collected with the CMS detector at the CERN LHC in 2016. A boosted decision tree is used to separate the tW signal from the dominant $\mathrm{t\bar{t}}$ background, whilst the subleading W+jets and multijet backgrounds are constrained using data-based estimates. This result is the first observation of the tW process in final states containing a muon or electron and jets, with a significance exceeding 5 standard deviations. The cross section is determined to be 89 $\pm$ 4 (stat) $\pm$ 12 (syst) pb, consistent with the standard model.

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.

This Letter presents the observation of the rare Z boson decay Z $\to\psi\ell^+\ell^-$. Here, $\psi$ represents contributions from direct J/$\psi$ and $\psi$(2S) $\to$ J/$\psi X$, $\ell^+\ell^-$ is a pair of electrons or muons, and the J/$\psi$ meson is detected via its decay to $\mu^+\mu^-$. The sample of proton-proton collision data, collected by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV, corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The signal is observed with a significance in excess of 5 standard deviations. After subtraction of the $\psi$(2S) $\to$ J/$\psi X$ contribution, the ratio of the branching fraction of the exclusive decay Z $\to\psi\ell^+\ell^-$ to the decay Z $\to\mu^+\mu^-\mu^+\mu^-$ within a fiducial phase space is measured to be $\mathcal{B}($Z $\to\psi\ell^+\ell^-) / \mathcal{B}($Z $\to\mu^+\mu^-\mu^+\mu^-) =$ 0.67 $\pm$ 0.18 (stat) $\pm$ 0.05 (syst).

Three of the most significant measured deviations from standard model predictions, the enhanced decay rate for B $\to$ D$^{(*)}\tau\nu$, hints of lepton universality violation in B $\to$ K$^{(*)}\ell\ell$ decays, and the anomalous magnetic moment of the muon, can be explained by the existence of leptoquarks (LQs) with large couplings to third-generation quarks and masses at the TeV scale. The existence of these states can be probed at the LHC in high energy proton-proton collisions. A novel search is presented for pair production of LQs coupled to a top quark and a muon using data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded by the CMS experiment. No deviation from the standard model prediction has been observed and scalar LQs decaying exclusively into t$\mu$ are excluded up to masses of 1420 GeV. The results of this search are combined with those from previous searches for LQ decays into t$\tau$ and b$\nu$, which excluded scalar LQs below masses of 900 and 1080 GeV. Vector LQs are excluded up to masses of 1190 GeV for all possible combinations of branching fractions to t$\mu$, t$\tau$ and b$\nu$. With this analysis, all relevant couplings of LQs with an electric charge of -1/3 to third-generation quarks are probed for the first time.

A search for resonances decaying into a W boson and a radion, where the radion decays into two W bosons, is presented. The data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$ recorded in proton-proton collisions with the CMS detector at $\sqrt{s} =$ 13 TeV. One isolated charged lepton is required, together with missing transverse momentum and one or two massive large-radius jets, containing the decay products of either two or one W bosons, respectively. No excess over the background estimation is observed. The results are combined with those from a complementary channel with an all-hadronic final state, described in an accompanying paper. Limits are set on parameters of an extended warped extra-dimensional model. These searches are the first of their kind at the LHC.

Multiparticle azimuthal correlations of prompt D$^0$ mesons are measured in PbPb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. For the first time, a four-particle cumulant method is used to extract the second Fourier coefficient of the azimuthal distribution ($v_2$) of D$^0$ mesons as a function of event centrality and the D$^0$ transverse momentum. The ratios of the four-particle $v_2$ values to previously measured two-particle cumulant results provide direct experimental access to event-by-event fluctuations of charm quark azimuthal anisotropies. These ratios are also found to be comparable to those of inclusive charged particles in the event. However, hints of deviations are seen in the most central and peripheral collisions. To investigate the origin of flow fluctuations in the charm sector, these measurements are compared with models implementing fluctuations of charm quark energy loss via collisional or radiative processes in the quark-gluon plasma. These models cannot quantitatively describe the data over the full transverse momentum and centrality ranges, although the calculations with collisional energy loss provide a better description of the data.

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}$ = 0.1179 $\pm$ 0.0019.

A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{-1}$ of data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum $k_\mathrm{T}$ and the splitting angle $\Delta R$ of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of $\ln(k_\mathrm{T}$/GeV) and $\ln(R/\Delta R)$ is measured for jets with distance parameters $R$ = 0.4 or 0.8, transverse momentum $p_\mathrm{T} \gt$ 700 GeV, and rapidity $\vert y\vert \lt $ 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.

A measurement of vector boson scattering and constraints on anomalous quartic gauge couplings from events with two Z bosons and two jets are presented. The analysis is based on a data sample of proton-proton collisions at sqrt(s) = 13 TeV collected with the CMS detector and corresponding to an integrated luminosity of 35.9 inverse femtobarns. The search is performed in the fully leptonic final state ZZ to lll'l', where l, l' = e, mu. The electroweak production of two Z bosons in association with two jets is measured with an observed (expected) significance of 2.7 (1.6) standard deviations. A fiducial cross section for the electroweak production is measured to be sigma[EW](pp -> ZZjj -> lll'l'jj) = 0.40 -0.16 +0.21 (stat) -0.09 +0.13 (syst) fb, which is consistent with the standard model prediction. Limits on anomalous quartic gauge couplings are determined in terms of the effective field theory operators T0, T1, T2, T8, and T9. This is the first measurement of vector boson scattering in the ZZ channel at the LHC.

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.

Results are reported from a search for supersymmetric particles in proton-proton collisions in the final state with a single lepton; multiple jets, including at least one b-tagged jet; and large missing transverse momentum. The search uses a sample of proton-proton collision data at sqrt(s) = 13 TeV recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 inverse femtobarns. The observed event yields in the signal regions are consistent with those expected from standard model backgrounds. The results are interpreted in the context of simplified models of supersymmetry involving gluino pair production, with gluino decay into either on- or off-mass-shell top squarks. Assuming that the top squarks decay into a top quark plus a stable, weakly interacting neutralino, scenarios with gluino masses up to about 1.9 TeV are excluded at 95% confidence level for neutralino masses up to about 1 TeV.

A search is performed for neutral non-standard-model Higgs bosons decaying to two muons in the context of the minimal supersymmetric standard model (MSSM). Proton-proton collision data recorded by the CMS experiment at the CERN Large Hadron Collider at a center-of-mass energy of 13 TeV were used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search is sensitive to neutral Higgs bosons produced via the gluon fusion process or in association with a $\mathrm{b\overline{b}}$ quark pair. No significant deviations from the standard model expectation are observed. Upper limits at 95% confidence level are set in the context of the $m_\mathrm{h}^{\text{mod+}}$ and phenomenological MSSM scenarios on the parameter $\tan\beta$ as a function of the mass of the pseudoscalar A boson, in the range from 130 to 600 GeV. The results are also used to set a model-independent limit on the product of the branching fraction for the decay into a muon pair and the cross section for the production of a scalar neutral boson, either via gluon fusion, or in association with b quarks, in the mass range from 130 to 1000 GeV.