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.

Two-particle Bose-Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data sample, containing 4.27 $\times$$10^{9}$ minimum bias events corresponding to an integrated luminosity of 0.607 nb$^{-1}$, was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter ($\alpha$), the Lévy scale parameter ($R$), and the correlation strength parameter ($\lambda$). The source shape, characterized by $\alpha$, is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy $R$ parameter. The $\lambda$ parameter is studied in terms of the core-halo model.

A search for long-lived particles decaying in the outer regions of the CMS silicon tracker or in the calorimeters is presented. The search is based on a data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018, corresponding to an integrated luminosity of 138 fb$^{-1}$. A novel technique, using trackless and out-of-time jet information combined in a deep neural network discriminator, is employed to identify decays of long-lived particles. The results are interpreted in a simplified model of chargino-neutralino production, where the neutralino is the next-to-lightest supersymmetric particle, is long-lived, and decays to a gravitino and either a Higgs or Z boson. This search is most sensitive to neutralino proper decay lengths of approximately 0.5 m, for which masses up to 1.18 TeV are excluded at 95% confidence level. The current search is the best result to date in the mass range from the kinematic limit imposed by the Higgs mass up to 1.8 TeV.

Pseudorapidity, transverse momentum, and multiplicity distributions are measured in the pseudorapidity range $|\eta| <$ 2.4 for charged particles with transverse momenta satisfying $p_\mathrm{T} >$ 0.5 GeV in proton-proton collisions at a center-of-mass energy of $\sqrt{s} =$ 13 TeV. Measurements are presented in three different event categories. The most inclusive of the categories corresponds to an inelastic pp data set, while the other two categories are exclusive subsets of the inelastic sample that are either enhanced or depleted in single diffractive dissociation events. The measurements are compared to predictions from Monte Carlo event generators used to describe high-energy hadronic interactions in collider and cosmic-ray physics.

This Letter reports measurements of differential cross sections for the production of two Z bosons in association with jets in proton-proton collisions at $\sqrt{s} =$ 8 and 13 TeV. The analysis is based on data samples collected at the LHC with the CMS detector, corresponding to integrated luminosities of 19.7 and 35.9 fb$^{-1}$ at 8 and 13 TeV, respectively. The measurements are performed in the leptonic decay modes ZZ $\to\ell^+ \ell^- \ell'^+ \ell'^-$, where $\ell,\ell' =$ e, $\mu$. The differential cross sections as a function of the jet multiplicity, the transverse momentum $p_\mathrm{T}$, and pseudorapidity of the $p_\mathrm{T}$-leading and subleading jets are presented. In addition, the differential cross sections as a function of variables sensitive to the vector boson scattering, such as the invariant mass of the two $p_\mathrm{T}$-leading jets and their pseudorapidity separation, are reported. The results are compared to theoretical predictions and found in good agreement within the theoretical and experimental uncertainties.

Measurements of inclusive isolated-photon and photon+jet production in proton-proton collisions at $\sqrt{s} =$ 13 TeV are presented. The analysis uses data collected by the CMS experiment in 2015, corresponding to an integrated luminosity of 2.26 fb$^{-1}$. The cross section for inclusive isolated photon production is measured as a function of the photon transverse energy in a fiducial region. The cross section for photon+jet production is measured as a function of the photon transverse energy in the same fiducial region with identical photon requirements and with the highest transverse momentum jet. All measurements are in agreement with predictions from next-to-leading-order perturbative QCD.

A measurement of the production cross section of top quark pairs in association with two b jets ($\mathrm{t\bar{t}}\mathrm{b\bar{b}}$) is presented using data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The cross section is measured in the all-jet decay channel of the top quark pair by selecting events containing at least eight jets, of which at least two are identified as originating from the hadronization of b quarks. A combination of multivariate analysis techniques is used to reduce the large background from multijet events not containing a top quark pair, and to help discriminate between jets originating from top quark decays and other additional jets. The cross section is determined for the total phase space to be 5.5 $\pm$ 0.3 (stat)${}^{+1.6}_{-1.3}$ (syst) pb and also measured for two fiducial $\mathrm{t\bar{t}}\mathrm{b\bar{b}}$ definitions. The measured cross sections are found to be larger than theoretical predictions by a factor of 1.5-2.4, corresponding to 1-2 standard deviations.

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.

A measurement of jet substructure observables is presented using \ttbar events in the lepton+jets channel from proton-proton collisions at $\sqrt{s}=$ 13 TeV recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Multiple jet substructure observables are measured for jets identified as bottom, light-quark, and gluon jets, as well as for inclusive jets (no flavor information). The results are unfolded to the particle level and compared to next-to-leading-order predictions from POWHEG interfaced with the parton shower generators PYTHIA 8 and HERWIG 7, as well as from SHERPA 2 and DIRE2. A value of the strong coupling at the Z boson mass, $\alpha_S(m_\mathrm{Z}) = $ 0.115$^{+0.015}_{-0.013}$, is extracted from the substructure data at leading-order plus leading-log accuracy.

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