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.

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

A search for the production of three massive vector bosons in proton--proton collisions is performed using data at $\sqrt{s}=13\,TeV$ recorded with the ATLAS detector at the Large Hadron Collider in the years 2015--2017, corresponding to an integrated luminosity of $79.8\,\text{fb}^{-1}$. Events with two same-sign leptons $\ell$ (electrons or muons) and at least two reconstructed jets are selected to search for $WWW\to\ell\nu\ell\nu qq$. Events with three leptons without any same-flavour opposite-sign lepton pairs are used to search for $WWW\to\ell\nu\ell\nu\ell\nu$, while events with three leptons and at least one same-flavour opposite-sign lepton pair and one or more reconstructed jets are used to search for $WWZ\to\ell\nu qq \ell\ell$. Finally, events with four leptons are analysed to search for $WWZ\to\ell\nu\ell\nu\ell\ell$ and $WZZ\to qq \ell\ell\ell\ell$. Evidence for the joint production of three massive vector bosons is observed with a significance of 4.0 standard deviations, where the expectation is 3.1 standard deviations.

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.

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 search is presented for photonic signatures, motivated by generalized models of gauge-mediated supersymmetry breaking. This search makes use of proton-proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded by the ATLAS detector at the LHC, and it explores models dominated by both strong and electroweak production of supersymmetric partner states. Experimental signatures incorporating an isolated photon and significant missing transverse momentum are explored. These signatures include events with an additional photon or additional jet activity not associated with any specific underlying quark flavor. No significant excess of events is observed above the Standard Model prediction, and 95% confidence-level upper limits of between 0.083 fb and 0.32 fb are set on the visible cross section of contributions from physics beyond the Standard Model. These results are interpreted in terms of lower limits on the masses of gluinos, squarks, and gauginos in the context of generalized models of gauge-mediated supersymmetry, which reach as high as 2.3 TeV for strongly produced and 1.3 TeV for weakly produced supersymmetric partner pairs.

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.

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 the four-lepton invariant mass spectrum is made with the ATLAS detector, using an integrated luminosity of 36.1 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}$ = 13 TeV delivered by the Large Hadron Collider. The differential cross-section is measured for events containing two same-flavour opposite-sign lepton pairs. It exhibits a rich structure, with different mass regions dominated in the Standard Model by single $Z$ boson production, Higgs boson production, and $Z$ boson pair production, and non-negligible interference effects at high invariant masses. The measurement is compared with state-of-the-art Standard Model calculations, which are found to be consistent with the data. These calculations are used to interpret the data in terms of $gg\rightarrow ZZ \rightarrow 4\ell$ and $Z \rightarrow 4\ell$ subprocesses, and to place constraints on a possible contribution from physics beyond the Standard Model.

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 result of a search for pair production of the supersymmetric partner of the Standard Model bottom quark ($\tilde{b}_1$) is reported. The search uses 3.2 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=$13 TeV collected by the ATLAS experiment at the Large Hadron Collider in 2015. Bottom squarks are searched for in events containing large missing transverse momentum and exactly two jets identified as originating from $b$-quarks. No excess above the expected Standard Model background yield is observed. Exclusion limits at 95% confidence level on the mass of the bottom squark are derived in phenomenological supersymmetric $R$-parity-conserving models in which the $\tilde{b}_1$ is the lightest squark and is assumed to decay exclusively via $\tilde{b}_1 \rightarrow b \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ is the lightest neutralino. The limits significantly extend previous results; bottom squark masses up to 800 (840) GeV are excluded for the $\tilde{\chi}_1^0$ mass below 360 (100) GeV whilst differences in mass above 100 GeV between the $\tilde{b}_1$ and the $\tilde{\chi}_1^0$ are excluded up to a $\tilde{b}_1$ mass of 500 GeV.

Results from a search for supersymmetry in events with four or more charged leptons (electrons, muons and taus) are presented. The analysis uses a data sample corresponding to 36.1 fb$^{-1}$ of proton-proton collisions delivered by the Large Hadron Collider at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying taus are designed to target a range of supersymmetric scenarios that can be either enriched in or depleted of events involving the production and decay of a $Z$ boson. Data yields are consistent with Standard Model expectations and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of General Gauge Mediated supersymmetry, where higgsino masses are excluded up to 295 GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of 1.46 TeV, 1.06 TeV, and 2.25 TeV are placed on wino, slepton and gluino masses, respectively.

The prevalence of hadronic jets at the LHC requires that a deep understanding of jet formation and structure is achieved in order to reach the highest levels of experimental and theoretical precision. There have been many measurements of jet substructure at the LHC and previous colliders, but the targeted observables mix physical effects from various origins. Based on a recent proposal to factorize physical effects, this Letter presents a double-differential cross-section measurement of the Lund jet plane using 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collision data collected with the ATLAS detector using jets with transverse momentum above 675 GeV. The measurement uses charged particles to achieve a fine angular resolution and is corrected for acceptance and detector effects. Several parton shower Monte Carlo models are compared with the data. No single model is found to be in agreement with the measured data across the entire plane.

This letter presents a combination of searches for Higgs boson pair production using up to 36.1 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV recorded with the ATLAS detector at the LHC. The combination is performed using six analyses searching for Higgs boson pairs decaying into the bbbb, bbWW, bb$\tau\tau$, WWWW, bb$\gamma \gamma$ and WW$\gamma\gamma$ final states. Results are presented for non-resonant and resonant Higgs boson pair production modes. No statistically significant excess in data above the Standard Model predictions is found. The combined observed (expected) limit at 95% confidence level on the non-resonant Higgs boson pair production cross-section is 6.9 (10) times the predicted Standard Model cross-section. Limits are also set on the ratio ($ \kappa_{\lambda} $) of the Higgs boson self-coupling to its Standard Model value. This ratio is constrained at 95% confidence level in observation (expectation) to $ -5.0 < \kappa_{\lambda} <12.0 $ ($ -5.8 < \kappa_{\lambda} <12.0 $). In addition, limits are set on the production of narrow scalar resonances and spin-2 Kaluza-Klein Randall-Sundrum gravitons. Exclusion regions are also provided in the parameter space of the habemus Minimal Supersymmetric Standard Model and the Electroweak Singlet Model.

Jet substructure quantities are measured using jets groomed with the soft-drop grooming procedure in dijet events from 32.9 fb$^{-1}$ of $pp$ collisions collected with the ATLAS detector at $\sqrt{s} = 13$ TeV. These observables are sensitive to a wide range of QCD phenomena. Some observables, such as the jet mass and opening angle between the two subjets which pass the soft-drop condition, can be described by a high-order (resummed) series in the strong coupling constant $\alpha_S$. Other observables, such as the momentum sharing between the two subjets, are nearly independent of $\alpha_S$. These observables can be constructed using all interacting particles or using only charged particles reconstructed in the inner tracking detectors. Track-based versions of these observables are not collinear safe, but are measured more precisely, and universal non-perturbative functions can absorb the collinear singularities. The unfolded data are directly compared with QCD calculations and hadron-level Monte Carlo simulations. The measurements are performed in different pseudorapidity regions, which are then used to extract quark and gluon jet shapes using the predicted quark and gluon fractions in each region. All of the parton shower and analytical calculations provide an excellent description of the data in most regions of phase space.

A search for production of the supersymmetric partners of the top quark, top squarks, is presented. The search is based on proton-proton collision events containing multiple jets, no leptons, and large transverse momentum imbalance. 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 137 fb$^{-1}$. The targeted signal production scenarios are direct and gluino-mediated top squark production, including scenarios in which the top squark and neutralino masses are nearly degenerate. The search utilizes novel algorithms based on deep neural networks that identify hadronically decaying top quarks and W bosons, which are expected in many of the targeted signal models. No statistically significant excess of events is observed relative to the expectation from the standard model, and limits on the top squark production cross section are obtained in the context of simplified supersymmetric models for various production and decay modes. Exclusion limits as high as 1310 GeV are established at the 95% confidence level on the mass of the top squark for direct top squark production models, and as high as 2260 GeV on the mass of the gluino for gluino-mediated top squark production models. These results represent a significant improvement over the results of previous searches for supersymmetry by CMS in the same final state.

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.

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

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.

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.

A search for direct pair production of top squarks in final states with two tau leptons, $b$-jets, and missing transverse momentum is presented. The analysis is based on proton-proton collision data at $\sqrt{s} = 13$ TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. Two exclusive channels with either two hadronically decaying tau leptons or one hadronically and one leptonically decaying tau lepton are considered. No significant deviation from the Standard Model predictions is observed in the data. The analysis results are interpreted in terms of model-independent limits and used to derive exclusion limits on the masses of the top squark $\tilde t_1$ and the tau slepton $\tilde \tau_1$ in a simplified model of supersymmetry with a nearly massless gravitino. In this model, masses up to $m(\tilde t_1) = 1.16$ TeV and $m(\tilde \tau_1) = 1.00$ TeV are excluded at 95% confidence level.

A search for the direct production of charginos and neutralinos in final states with at least two hadronically decaying tau leptons is presented. The analysis uses a dataset of $pp$ collisions corresponding to an integrated luminosity of 36.1 fb$^{-1}$, recorded with the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. No significant deviation from the Standard Model background expectation is observed. Limits are derived in scenarios of $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ pair production and of $\tilde{\chi}_{1}^{\pm}\tilde{\chi}_{2}^{0}$ and $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ production in simplified models where the neutralinos and charginos decay solely via intermediate left-handed staus and tau sneutrinos, and the mass of the $\tilde{\tau}_{\mathrm L}$ state is set to be halfway between the masses of the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$. Chargino masses up to 630 GeV are excluded at 95% confidence level in the scenario of direct production of $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ for a massless $\tilde{\chi}_{1}^{0}$. Common $\tilde{\chi}_{1}^{\pm}$, $\tilde{\chi}_{2}^{0}$ masses up to 760 GeV are excluded in the case of production of $\tilde{\chi}_{1}^{\pm}\tilde{\chi}_{2}^{0}$ and $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ assuming a massless $\tilde{\chi}_{1}^{0}$. Exclusion limits for additional benchmark scenarios with large and small mass-splitting between the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$ are also studied by varying the $\tilde{\tau}_{\mathrm L}$ mass between the masses of the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$.

Narrow resonances decaying into $WW$, $WZ$ or $ZZ$ boson pairs are searched for in 36.7 fb $^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. The diboson system is reconstructed using pairs of large-radius jets with high transverse momentum and tagged as compatible with the hadronic decay of high-momentum $W$ or $Z$ bosons, using jet mass and substructure properties. The search is sensitive to diboson resonances with masses in the range 1.2-5.0 TeV. No significant excess is observed in any signal region. Exclusion limits are set at the 95% confidence level on the production cross section times branching ratio to dibosons for a range of theories beyond the Standard Model. Model-dependent lower limits on the mass of new gauge bosons are set, with the highest limit set at 3.5 TeV in the context of mass-degenerate resonances that couple predominantly to bosons.

Multi-particle cumulants and corresponding Fourier harmonics are measured for azimuthal angle distributions of charged particles in $pp$ collisions at $\sqrt{s}$ = 5.02 and 13 TeV and in $p$+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, and compared to the results obtained for low-multiplicity Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These measurements aim to assess the collective nature of particle production. The measurements of multi-particle cumulants confirm the evidence for collective phenomena in $p$+Pb and low-multiplicity Pb+Pb collisions. On the other hand, the $pp$ results for four-particle cumulants do not demonstrate collective behaviour, indicating that they may be biased by contributions from non-flow correlations. A comparison of multi-particle cumulants and derived Fourier harmonics across different collision systems is presented as a function of the charged-particle multiplicity. For a given multiplicity, the measured Fourier harmonics are largest in Pb+Pb, smaller in $p$+Pb and smallest in $pp$ collisions. The $pp$ results show no dependence on the collision energy, nor on the multiplicity.

A search for magnetic monopoles and high-electric-charge objects is presented using 34.4 fb$^{-1}$ of 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC during 2015 and 2016. The considered signature is based upon high ionization in the transition radiation tracker of the inner detector associated with a pencil-shape energy deposit in the electromagnetic calorimeter. The data were collected by a dedicated trigger based on the tracker high-threshold hit capability. The results are interpreted in models of Drell-Yan pair production of stable particles with two spin hypotheses (0 and 1/2) and masses ranging from 200 GeV to 4000 GeV. The search improves by approximately a factor of five the constraints on the direct production of magnetic monopoles carrying one or two Dirac magnetic charges and stable objects with electric charge in the range $20\le|z|\le60$ and extends the charge range to $60<|z|\le100$.

The dynamics of isolated-photon plus two-jet production in $pp$ collisions at a centre-of-mass energy of 13 TeV are studied with the ATLAS detector at the LHC using a dataset corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Cross sections are measured as functions of a variety of observables, including angular correlations and invariant masses of the objects in the final state, $\gamma+jet+jet$. Measurements are also performed in phase-space regions enriched in each of the two underlying physical mechanisms, namely direct and fragmentation processes. The measurements cover the range of photon (jet) transverse momenta from 150 GeV (100 GeV) to 2 TeV. The tree-level plus parton-shower predictions from SHERPA and PYTHIA as well as the next-to-leading-order QCD predictions from SHERPA are compared with the measurements. The next-to-leading-order QCD predictions describe the data adequately in shape and normalisation except for regions of phase space such as those with high values of the invariant mass or rapidity separation of the two jets, where the predictions overestimate the data.

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.

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.

Results of a search for the pair production of photon-jets$-$collimated groupings of photons$-$in the ATLAS detector at the Large Hadron Collider are reported. Highly collimated photon-jets can arise from the decay of new, highly boosted particles that can decay to multiple photons collimated enough to be identified in the electromagnetic calorimeter as a single, photonlike energy cluster. Data from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 36.7 fb$^{-1}$, were collected in 2015 and 2016. Candidate photon-jet pair production events are selected from those containing two reconstructed photons using a set of identification criteria much less stringent than that typically used for the selection of photons, with additional criteria applied to provide improved sensitivity to photon-jets. Narrow excesses in the reconstructed diphoton mass spectra are searched for. The observed mass spectra are consistent with the Standard Model background expectation. The results are interpreted in the context of a model containing a new, high-mass scalar particle with narrow width, $X$, that decays into pairs of photon-jets via new, light particles, $a$. Upper limits are placed on the cross section times the product of branching ratios $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow \gamma \gamma)^{2}$ for 200 GeV $< m_{X} <$ 2 TeV and for ranges of $ m_a $ from a lower mass of 100 MeV up to between 2 and 10 GeV, depending upon $ m_X $. Upper limits are also placed on $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow 3\pi^{0})^{2}$ for the same range of $ m_X $ and for ranges of $ m_a $ from a lower mass of 500 MeV up to between 2 and 10 GeV.

Jet substructure observables have significantly extended the search program for physics beyond the Standard Model at the Large Hadron Collider. The state-of-the-art tools have been motivated by theoretical calculations, but there has never been a direct comparison between data and calculations of jet substructure observables that are accurate beyond leading-logarithm approximation. Such observables are significant not only for probing the collinear regime of QCD that is largely unexplored at a hadron collider, but also for improving the understanding of jet substructure properties that are used in many studies at the Large Hadron Collider. This Letter documents a measurement of the first jet substructure quantity at a hadron collider to be calculated at next-to-next-to-leading-logarithm accuracy. The normalized, differential cross-section is measured as a function of log$_{10}\rho^2$, where $\rho$ is the ratio of the soft-drop mass to the ungroomed jet transverse momentum. This quantity is measured in dijet events from 32.9 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collisions recorded by the ATLAS detector. The data are unfolded to correct for detector effects and compared to precise QCD calculations and leading-logarithm particle-level Monte Carlo simulations.

This paper presents a measurement of quantities related to the formation of jets from high-energy quarks and gluons (fragmentation). Jets with transverse momentum 100 GeV $<p_T<$ 2.5 TeV and pseudorapidity $|\eta| < 2.1$ from an integrated luminosity of 33 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions are reconstructed with the ATLAS detector at the Large Hadron Collider. Charged-particle tracks with $p_T > 500$ MeV and $|\eta| < 2.5$ are used to probe the detailed structure of the jet. The fragmentation properties of the more forward and the more central of the two leading jets from each event are studied. The data are unfolded to correct for detector resolution and acceptance effects. Comparisons with parton shower Monte Carlo generators indicate that existing models provide a reasonable description of the data across a wide range of phase space, but there are also significant differences. Furthermore, the data are interpreted in the context of quark- and gluon-initiated jets by exploiting the rapidity dependence of the jet flavor fraction. A first measurement of the charged-particle multiplicity using model-independent jet labels (topic modeling) provides a promising alternative to traditional quark and gluon extractions using input from simulation. The simulations provide a reasonable description of the quark-like data across the jet $p_T$ range presented in this measurement, but the gluon-like data have systematically fewer charged particles than the simulations.

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.

A search is presented for the Higgs boson decay to a pair of electrons (e$^+$e$^-$) in proton-proton collisions at $\sqrt{s}$ = 13 TeV. The data set was collected with the CMS experiment at the LHC between 2016 and 2018, corresponding to an integrated luminosity of 138 fb$^{-1}$. The analysis uses event categories targeting Higgs boson production via gluon fusion and vector boson fusion. The observed upper limit on the Higgs boson branching fraction to an electron pair is 3.0 $\times$ 10$^{-4}$ (3.0 $\times$ 10$^{-4}$ expected) at the 95% confidence level, which is the most stringent limit on this branching fraction to date.

A search for resonances in events with the $\gamma$+jet final state has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC. The total data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$. Models of excited quarks and quantum black holes are considered. Using a wide-jet reconstruction for the candidate jet, the $\gamma$+jet invariant mass spectrum measured in data is examined for the presence of resonances over the standard model continuum background. The background is estimated by fitting the mass distribution with a functional form. The data exhibit no statistically significant deviations from the expected standard model background. Exclusion limits at 95% confidence level on the resonance mass and other parameters are set. Excited light-flavor quarks (excited bottom quarks) are excluded up to a mass of 6.0 (3.8) TeV. Quantum black hole production is excluded for masses up to 7.5 (5.2) TeV in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model. These lower mass bounds are the most stringent to date among those obtained in the $\gamma$+jet final state.

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

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb$^{-1}$. Candidate events are selected by requiring two oppositely-charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (u or c), and a down-type quark (d, s, or b). The results improve the previous bounds by three to six orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.