We report on the first measurement of inclusive J/$\psi$ elliptic flow, $v_2$, in heavy-ion collisions at the LHC. The measurement is performed with the ALICE detector in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV in the rapidity range $2.5 < y < 4.0$. The dependence of the J/$\psi$ $v_2$ on the collision centrality and on the J/$\psi$ transverse momentum is studied in the range $0 < p_{\rm T} < 10$ GeV/$c$. For semi-central Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV, an indication of non-zero $v_2$ is observed with a maximum value of $v_2 = 0.116 \pm 0.046 {\rm (stat.)} \pm 0.029 {\rm (syst.)}$ for J/$\psi$ in the transverse momentum range $2 < p_{\rm T} < 4$ GeV/$c$. The elliptic flow measurement complements the previously reported ALICE results on the inclusive J/$\psi$ nuclear modification factor and favors the scenario of a significant fraction of J/$\psi$ production from charm quarks in a deconfined partonic phase.

We present the first wide-range measurement of the charged-particle pseudorapidity density distribution, for different centralities (the 0-5%, 5-10%, 10-20%, and 20-30% most central events) in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV at the LHC. The measurement is performed using the full coverage of the ALICE detectors, $-5.0 < \eta < 5.5$, and employing a special analysis technique based on collisions arising from LHC "satellite" bunches. We present the pseudorapidity density as a function of the number of participating nucleons as well as an extrapolation to the total number of produced charged particles ($N_{\rm ch} = 17165 \pm 772$ for the 0-5% most central collisions). From the measured ${\rm d}N_{\rm ch}/{\rm d}\eta$ distribution we derive the rapidity density distribution, ${\rm d}N_{\rm ch}/{\rm d}y$, under simple assumptions. The rapidity density distribution is found to be significantly wider than the predictions of the Landau model. We assess the validity of longitudinal scaling by comparing to lower energy results from RHIC. Finally the mechanisms of the underlying particle production are discussed based on a comparison with various theoretical models.

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

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.

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.

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.

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.

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.

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.

The first evidence for the standard model production of a top quark in association with a W boson and a Z boson is reported. The measurement is performed in multilepton final states, where the Z boson is reconstructed via its decays to electron or muon pairs and the W boson decays either to leptons or hadrons. The analysed data were recorded by the CMS experiment at the CERN LHC in 2016-2018 in proton-proton collisions at $\sqrt{s}$ = 13 TeV, and correspond to an integrated luminosity of 138 fb$^{-1}$. The measured cross section is 354 $\pm$ 54 (stat) $\pm$ 95 (syst) fb, and corresponds to a statistical significance of 3.4 standard deviations.

Measurements of the charge-dependent two-particle angular correlation function in proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at$\sqrt{s_\mathrm{NN}}$ = 5.02 TeV are reported. The pPb and PbPb datasets correspond to integrated luminosities of 186\nbinv and 0.607 nb$^{-1}$, respectively, and were collected using the CMS detector at the CERN LHC. The charge-dependent correlations are characterized by balance functions of same- and opposite-sign particle pairs. The balance functions, which contain information about the creation time of charged particle pairs and the development of collectivity, are studied as functions of relative pseudorapidity ($\Delta \eta$) and relative azimuthal angle ($\Delta \phi$), for various multiplicity and transverse momentum ($p_\mathrm{T}$) intervals. A multiplicity dependence of the balance function is observed in $\Delta \eta$ and $\Delta \phi$ for both systems. The width of the balance functions decreases towards high-multiplicity collisions in the momentum region $\lt$2 GeV, for pPb and PbPb results. No multiplicity dependence is observed at higher transverse momentum. The data are compared with HYDJET, HIJING and AMPT generator predictions, none of which capture completely the multiplicity dependence seen in the data.

The central exclusive production of charged-hadron pairs in pp collisions at a centre-of-mass energy of 13\TeV is examined, based on data collected in a special high-$\beta^*$ run of the LHC. The nonresonant continuum processes are studied with the invariant mass of the centrally produced two-pion system in the resonance-free region, $m_{\pi^+\pi^-}$$\lt$ 0.7 GeV or $m_{\pi^+\pi^-}$$\gt$ 1.8 GeV. Differential cross sections as functions of the azimuthal angle between the surviving protons, squared exchanged four-momenta, and $m_{\pi^+\pi^-}$ are measured in a wide region of scattered proton transverse momenta, between 0.2 and 0.8 GeV, and for pion rapidities $\lvert y\rvert$$\lt$ 2. A rich structure of interactions related to double-pomeron exchange is observed. A parabolic minimum in the distribution of the two-proton azimuthal angle is observed for the first time. It can be interpreted as an effect of additional pomeron exchanges between the protons from the interference between the bare and the rescattered amplitudes. After model tuning, various physical quantities are determined that are related to the pomeron cross section, proton-pomeron and meson-pomeron form factors, pomeron trajectory and intercept, and coefficients of diffractive eigenstates of the proton.

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.

The leptonic and inclusive hadronic decay branching fractions of the W boson are measured using proton-proton collision data collected at $\sqrt{s} =$ 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Events characterized by the production of one or two W bosons are selected and categorized based on the multiplicity and flavor of reconstructed leptons, the number of jets, and the number of jets identified as originating from the hadronization of b quarks. A binned maximum likelihood estimate of the W boson branching fractions is performed simultaneously in each event category. The measured branching fractions of the W boson decaying into electron, muon, and tau lepton final states are (10.83 $\pm$ 0.10)%, (10.94 $\pm$ 0.08)%, and (10.77 $\pm$ 0.21)%, respectively, consistent with lepton flavor universality for the weak interaction. The average leptonic and inclusive hadronic decay branching fractions are estimated to be (10.89 $\pm$ 0.08)% and (67.32 $\pm$ 0.23)%, respectively. Based on the hadronic branching fraction, three standard model quantities are subsequently derived: the sum of squared elements in the first two rows of the Cabibbo-Kobayashi-Maskawa (CKM) matrix $\sum_{ij}\vert V_{ij}\vert^{2}$ = 1.984 $\pm$ 0.021, the CKM element $\vert V_\mathrm{cs}\vert$ = 0.967 $\pm$ 0.011, and the strong coupling constant at the W boson mass scale, $\alpha_\mathrm{S}(m^2_\mathrm{W})$ = 0.095 $\pm$ 0.033.

The inclusive jet cross section is measured as a function of jet transverse momentum $p_\mathrm{T}$ and rapidity $y$. The measurement is performed using proton-proton collision data at $\sqrt{s}$ = 5.02 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 27.4 pb$^{-1}$. The jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm using a distance parameter of $R$ = 0.4, within the rapidity interval $\lvert y\rvert$$\lt$ 2, and across the kinematic range 0.06 $\lt$$p_\mathrm{T}$$\lt$ 1 TeV. The jet cross section is unfolded from detector to particle level using the determined jet response and resolution. The results are compared to predictions of perturbative quantum chromodynamics, calculated at both next-to-leading order and next-to-next-to-leading order. The predictions are corrected for nonperturbative effects, and presented for a variety of parton distribution functions and choices of the renormalization/factorization scales and the strong coupling $\alpha_\mathrm{S}$.

The hydrodynamic flow-like behavior of charged hadrons in high-energy lead-lead collisions is studied through multiparticle correlations. The elliptic anisotropy values based on different orders of multiparticle cumulants, $v_{2}\{2k\}$, are measured up to the tenth order ($k$ = 5) as functions of the collision centrality at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. A hierarchy is observed between the coefficients, with $v_{2}\{2\} > v_{2}\{4\} \gtrsim v_{2}\{6\} \gtrsim v_{2}\{8\} \gtrsim v_{2}\{10\}$. Based on these results, centrality-dependent moments for the fluctuation-driven event-by-event $v_{2}$ distribution are determined, including the skewness, kurtosis and, for the first time, superskewness. Assuming a hydrodynamic expansion of the produced medium, these moments directly probe the initial-state geometry in high-energy nucleus-nucleus collisions.

We present an observation of photon-photon production of $\tau$ lepton pairs in ultraperipheral lead-lead collisions. The measurement is based on a data sample with an integrated luminosity of 404 $\mu$b$^{-1}$ collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The $\gamma\gamma$$\to$$\tau^+\tau^-$ process is observed for $\tau\tau$ events with a muon and three charged hadrons in the final state. The measured fiducial cross section is $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$ = 4.8 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) $\mu$b, in agreement with leading-order QED predictions. Using $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$, we estimate a model-dependent value of the anomalous magnetic moment of the $\tau$ lepton of $a_\tau$ = 0.001 $^{+0.055}_{-0.089}$.

A search is presented for charged, long-lived supersymmetric particles in final states with one or more disappearing tracks. The search is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the CERN LHC between 2016 and 2018, corresponding to an integrated luminosity of 137 fb$^{-1}$. The search is performed over final states characterized by varying numbers of jets, b-tagged jets, electrons, and muons. The length of signal-candidate tracks in the plane perpendicular to the beam axis is used to characterize the lifetimes of wino- and higgsino-like charginos produced in the context of the minimal supersymmetric standard model. The d$E$/d$x$ energy loss of signal-candidate tracks is used to increase the sensitivity to charginos with a large mass and thus a small Lorentz boost. The observed results are found to be statistically consistent with the background-only hypothesis. Limits on the pair production cross section of gluinos and squarks are presented in the framework of simplified models of supersymmetric particle production and decay, and for electroweakino production based on models of wino and higgsino dark matter. The limits presented are the most stringent to date for scenarios with light third-generation squarks and a wino- or higgsino-like dark matter candidate capable of explaining the known dark matter relic density.

The azimuthal anisotropy of $\Upsilon$(1S) mesons in high-multiplicity proton-lead collisions is studied using data collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 8.16 TeV. The $\Upsilon$(1S) mesons are reconstructed using their dimuon decay channel. The anisotropy is characterized by the second Fourier harmonic coefficients, found using a two-particle correlation technique, in which the $\Upsilon$(1S) mesons are correlated with charged hadrons. A large pseudorapidity gap is used to suppress short-range correlations. Nonflow contamination from the dijet background is removed using a low-multiplicity subtraction method, and the results are presented as a function of $\Upsilon$(1S) transverse momentum. The azimuthal anisotropies are smaller than those found for charmonia in proton-lead collisions at the same collision energy, but are consistent with values found for $\Upsilon$(1S) mesons in lead-lead interactions at a nucleon-nucleon center-of-mass energy of 5.02 TeV.

A search for heavy neutral leptons (HNLs), the right-handed Dirac or Majorana neutrinos, is performed in final states with three charged leptons (electrons or muons) using proton-proton collision data collected by the CMS experiment at $\sqrt{s} =$ 13 TeV at the CERN LHC. The data correspond to an integrated luminosity of 138 fb$^{-1}$. The HNLs could be produced through mixing with standard model neutrinos $\nu$. For small values of the HNL mass ($\lt$ 20 GeV) and the square of the HNL-$\nu$ mixing parameter (10$^{-7}$-10$^{-2}$), the decay length of these particles can be large enough so that the secondary vertex of the HNL decay can be resolved with the CMS silicon tracker. The selected final state consists of one lepton emerging from the primary proton-proton collision vertex, and two leptons forming a displaced, secondary vertex. No significant deviations from the standard model expectations are observed, and constraints are obtained on the HNL mass and coupling strength parameters, excluding previously unexplored regions of parameter space in the mass range 1-20 GeV and squared mixing parameter values as low as 10$^{-7}$.

The measurement of the cross section for the production of a Z boson, decaying to dielectrons or dimuons, in association with at least one bottom quark jet are performed with proton-proton collision data at $\sqrt{s} =$ 13 TeV. The data sample corresponds to an integrated luminosity of 137 fb$^{-1}$, collected by the CMS experiment at the LHC during 2016-2018. The integrated cross sections for Z + $\ge$ 1 b jet and Z + $\ge$ 2 b jets are reported for the electron, muon, and combined channels. The fiducial cross sections in the combined channel are 6.52 $\pm$ 0.04 (stat) $\pm$ 0.40 (syst) $\pm$ 0.14 (theo) pb for Z + $\ge$ 1 b jet and 0.65 $\pm$ 0.03 (stat) $\pm$ 0.07 (syst) $\pm$ 0.02 (theo) pb for Z + $\ge$ 2 b jets. The differential cross section distributions are measured as functions of various kinematic observables that are useful for precision tests of perturbative quantum chromodynamics predictions. The ratios of integrated and differential cross sections for Z + $\ge$ 2 b jets and Z + $\ge$ 1 b jet processes are also determined. The value of the integrated cross section ratio measured in the combined channel is 0.100 $\pm$ 0.005 (stat) $\pm$ 0.007 (syst) $\pm$ 0.003 (theo). All measurements are compared with predictions from various event generators.

The first observation of the production of W$^\pm$W$^\pm$ bosons from double parton scattering processes using same-sign electron-muon and dimuon events in proton-proton collisions is reported. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ recorded at a center-of-mass energy of 13 TeV using the CMS detector at the CERN LHC. Multivariate discriminants are used to distinguish the signal process from the main backgrounds. A binned maximum likelihood fit is performed to extract the signal cross section. The measured cross section for production of same-sign W bosons decaying leptonically is 80.7 $\pm$ 11.2 (stat) $^{+9.5}_{-8.6}$ (syst) $\pm$ 12.1 (model) fb, whereas the measured fiducial cross section is 6.28 $\pm$ 0.81 (stat) $\pm$ 0.69 (syst) $\pm$ 0.37 (model) fb. The observed significance of the signal is 6.2 standard deviations above the background-only hypothesis.

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.

A search is presented for high-mass exclusive diphoton production via photon-photon fusion in proton-proton collisions at $\sqrt{s}$ = 13 TeV in events where both protons survive the interaction. The analysis utilizes data corresponding to an integrated luminosity of 103 fb$^{-1}$ collected in 2016-2018 with the central CMS detector and the CMS and TOTEM precision proton spectrometer (PPS). Events that have two photons with high transverse momenta ($p_\mathrm{T}^\gamma > $ 100 GeV), back-to-back in azimuth, and with a large diphoton invariant mass ($m_{\gamma\gamma} \gt$ 350 GeV) are selected. To remove the dominant inclusive diphoton backgrounds, the kinematic properties of the protons detected in PPS are required to match those of the central diphoton system. Only events having opposite-side forward protons detected with a fractional momentum loss between 0.035 and 0.15 (0.18) for the detectors on the negative (positive) side of CMS are considered. One exclusive diphoton candidate is observed for an expected background of 1.1 events. Limits at 95% confidence level are derived for the four-photon anomalous coupling parameters $\lvert\zeta_1\rvert \lt$ 0.073 TeV$^{-4}$ and $\lvert\zeta_2\rvert \lt$ 0.15 TeV$^{-4}$, using an effective field theory. Additionally, upper limits are placed on the production of axion-like particles with coupling strength to photons $f^{-1}$ that varies from 0.03 TeV$^{-1}$ to 1 TeV$^{-1}$ over the mass range from 500 to 2000 GeV.

The first measurement of the top quark pair ($\mathrm{t\bar{t}}$) production cross section in proton-proton collisions at $\sqrt{s}$ = 13.6 TeV is presented. Data recorded with the CMS detector at the CERN LHC in Summer 2022, corresponding to an integrated luminosity of 1.21 fb$^{-1}$, are analyzed. Events are selected with one or two charged leptons (electrons or muons) and additional jets. A maximum likelihood fit is performed in event categories defined by the number and flavors of the leptons, the number of jets, and the number of jets identified as originating from b quarks. An inclusive $\mathrm{t\bar{t}}$ production cross section of 881 $\pm$ 23 (stat+syst) $\pm$ 20 (lumi) pb is measured, in agreement with the standard model prediction of 924 $^{+32}_{-40}$ pb.

The polarization of $\tau$ leptons is measured using leptonic and hadronic $\tau$ lepton decays in Z $\to$$\tau^+\tau^-$ events in proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by CMS at the CERN LHC with an integrated luminosity of 36.3 fb$^{-1}$. The measured $\tau^-$ polarization at the Z boson mass pole is $\mathcal{P}_{\tau}$(Z) = $-$0.144 $\pm$ 0.006 (stat) $\pm$ 0.014 (syst) = $-$0.144 $\pm$ 0.015, in good agreement with the measurement of the $\tau$ lepton asymmetry parameter of $A_{\tau}$ = 0.1439 $\pm$ 0.0043 = $-\mathcal{P}_{\tau}$(Z) at LEP. The $\tau$ polarization depends on the ratio of the vector to axial-vector couplings of the $\tau$ leptons in the neutral current expression, and thus on the effective weak mixing angle $\sin^{2}\theta_\mathrm{W}^{\text{eff}}$, independently of the Z boson production mechanism. The obtained value $\sin^{2}\theta_\mathrm{W}^{\text{eff}}$ = 0.2319 $\pm$ 0.0008 (stat) $\pm$ 0.0018 (syst) = 0.2319 $\pm$ 0.0019 is in good agreement with measurements at e$^+$e$^-$ colliders.

The first study of the shapes of jets arising from bottom (b) quarks in heavy ion collisions is presented. Jet shapes are studied using charged hadron constituents as a function of their radial distance from the jet axis. Lead-lead (PbPb) collision data at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV were recorded by the CMS detector at the LHC, with an integrated luminosity of 1.69 nb$^{-1}$. Compared to proton-proton collisions, a redistribution of the energy in b jets to larger distances from the jet axis is observed in PbPb collisions. This medium-induced redistribution is found to be substantially larger for b jets than for inclusive jets.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

A search for partonic collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged constituents using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity $\lvert\eta\rvert$$\lt$ 1.6. Two-particle correlations among the charged constituents within the jets are studied as functions of the particles' azimuthal angle and pseudorapidity separations ($\Delta\phi^*$ and $\Delta\eta^*$) in a jet coordinate basis, where constituents' $\eta^*$, $\phi^*$ are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to $N_\text{ch}^\mathrm{j}$$\approx$ 100. Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for $\lvert\Delta\eta^*\rvert$$\gt$ 2. For low-multiplicity jets, the long-range elliptic anisotropic harmonic, $v^*_2$, is observed to decrease with $N_\text{ch}^\mathrm{j}$. This trend is well described by Monte Carlo event generators. However, a rising trend for $v^*_2$ emerges at $N_\text{ch}^\mathrm{j}$$\gtrsim$ 80, hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of parton fragmentation processes in the vacuum.

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.

Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark-gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we extracted the speed of sound in this medium created using lead-lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 $\pm$ 0.002 (stat) $\pm$ 0.016 (syst) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 $\pm$ 8 (syst) MeV. The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

Measurements of the elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients are presented for D$^0$ mesons produced in b hadron decays (nonprompt D$^0$ mesons) in lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The results are compared with previously published charm meson anisotropies measured using prompt D$^0$ mesons. The data were collected with the CMS detector in 2018 with an integrated luminosity of 0.58 nb$^{-1}$. Azimuthal anisotropy is sensitive to the interactions of quarks with the hot and dense medium created in heavy ion collisions. Comparing results for prompt and nonprompt D$^0$ mesons can assist in understanding the mass dependence of these interactions. The nonprompt results show lower magnitudes of $v_2$ and $v_3$ and weaker dependences on the meson transverse momentum and collision centrality than those found for prompt D$^0$ mesons. The results are in agreement with theoretical predictions that include a mass dependence in the interactions of quarks with the medium.

The top quark pair production cross section is measured in proton-proton collisions at a center-of-mass energy of 5.02 TeV. The data were collected in a special LHC low-energy and low-intensity run in 2017, and correspond to an integrated luminosity of 302 pb$^{-1}$. The measurement is performed using events with one electron and one muon of opposite charge, and at least two jets. The measured cross section is 60.7 $\pm$ 5.0 (stat) $\pm$ 2.8 (syst) $\pm$ 1.1 (lumi) pb. To reduce the statistical uncertainty, a combination with the result in the single lepton + jets channel, based on data collected in 2015 at the same center-of-mass energy and corresponding to an integrated luminosity of 27.4 pb$^{-1}$, is then performed. The resulting measured value is 63.0 $\pm$ 4.1 (stat) $\pm$ 3.0 (syst+lumi) pb, in agreement with the standard model prediction of 66.8 $^{+2.9}_{-3.1}$ pb.

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q $\to$ e$\mu$t) and decay (t $\to$ e$\mu$q) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an e$\mu$ pair and an up (charm) quark of 0.13 $\times$ 10$^{-6}$ (1.31 $\times$ 10$^{-6}$), 0.07 $\times$ 10$^{-6}$ (0.89 $\times$ 10$^{-6}$), and 0.25 $\times$ 10$^{-6}$ (2.59 $\times$ 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.

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

A search is presented for a right-handed W boson (W$_\mathrm{R}$) and a heavy neutrino (N), in a final state consisting of two same-flavor leptons (ee or $\mu\mu$) and two quarks. The search is performed with the CMS experiment at the CERN LHC using a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The search covers two regions of phase space, one where the decay products of the heavy neutrino are merged into a single large-area jet, and one where the decay products are well separated. The expected signal is characterized by an excess in the invariant mass distribution of the final-state objects. No significant excess over the standard model background expectations is observed. The observations are interpreted as upper limits on the product of W$_\mathrm{R}$ production cross sections and branching fractions assuming that couplings are identical to those of the standard model W boson. For N masses $m_\mathrm{N}$ equal to half the W$_\mathrm{R}$ mass $m_\mathrm{W_R}$ ($m_\mathrm{N}$ = 0.2 TeV), $m_\mathrm{W_R}$ is excluded at 95% confidence level up to 4.7 (4.8) and 5.0 (5.4) TeV for the electron and muon channels, respectively. This analysis provides the most stringent limits on the W$_\mathrm{R}$ mass to date.

Multijet events at large transverse momentum ($p_\mathrm{T}$) are measured at $\sqrt{s}$ = 13 TeV using data recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with $p_\mathrm{T}$$>$ 50 GeV that are produced in association with a high-$p_\mathrm{T}$ dijet system is measured in various ranges of the $p_\mathrm{T}$ of the jet with the highest transverse momentum and as a function of the azimuthal angle difference $\Delta\phi_{1,2}$ between the two highest $p_\mathrm{T}$ jets in the dijet system. The differential production cross sections are measured as a function of the transverse momenta of the four highest $p_\mathrm{T}$ jets. The measurements are compared with leading and next-to-leading order matrix element calculations supplemented with simulations of parton shower, hadronization, and multiparton interactions. In addition, the measurements are compared with next-to-leading order matrix element calculations combined with transverse-momentum dependent parton densities and transverse-momentum dependent parton shower.

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.