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

A search for a heavy resonance decaying into a top quark and a W boson in proton-proton collisions at $\sqrt{s} =$ 13 TeV is presented. The data analyzed were recorded with the CMS detector at the LHC and correspond to an integrated luminosity of 138 fb$^{-1}$. The top quark is reconstructed as a single jet and the W boson, from its decay into an electron or muon and the corresponding neutrino. A top quark tagging technique based on jet clustering with a variable distance parameter and simultaneous jet grooming is used to identify jets from the collimated top quark decay. The results are interpreted in the context of two benchmark models, where the heavy resonance is either an excited bottom quark b$^*$ or a vector-like quark B. A statistical combination with an earlier search by the CMS Collaboration in the all-hadronic final state is performed to place upper cross section limits on these two models. The new analysis extends the lower range of resonance mass probed from 1.4 down to 0.7 TeV. For left-handed, right-handed, and vector-like couplings, b$^*$ masses up to 3.0, 3.0, and 3.2 TeV are excluded at 95% confidence level, respectively. The observed upper limits represent the most stringent constraints on the b$^*$ model to date.

Using a data sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment in 2016-2018, the B$^0_\mathrm{s}$ $\to$ X(3872)$\phi$ decay is observed. Decays into J/$\psi\,\pi^+\pi^-$ and K$^+$K$^-$ are used to reconstruct, respectively, the X(3872) and $\phi$. The ratio of the product of branching fractions $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi)\,\mathcal{B}($X(3872) $\to$ J$/\psi\,\pi^+\pi^-)$ to the product $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ $\psi$(2S)$\phi)\,\mathcal{B}(\psi$(2S) $\to$ J/$\psi\,\pi^+\pi^-)$ is measured to be (2.21 $\pm$ 0.29 (stat) $\pm$ 0.17 (syst))%. The ratio $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi) / \mathcal{B}($B$^{0}$ $\to$ X(3872)K$^0)$ is found to be consistent with one, while the ratio $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi) / \mathcal{B}($B$^+$ $\to$ X(3872)K$^+)$ is two times smaller. This suggests a difference in the production dynamics of the X(3872) in B$^0$ and B$^0_\mathrm{s}$ meson decays compared to B$^+$. The reported observation may shed new light on the nature of the X(3872) particle.

A search is conducted for a low-mass charged Higgs boson produced in a top quark decay and subsequently decaying into a charm and a strange quark. The data sample was recorded in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS experiment at the LHC and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The search is performed in the process of top quark pair production, where one top quark decays to a bottom quark and a charged Higgs boson, and the other to a bottom quark and a W boson. With the W boson decaying to a charged lepton (electron or muon) and a neutrino, the final state comprises an isolated lepton, missing transverse momentum, and at least four jets, of which two are tagged as b jets. To enhance the search sensitivity, one of the jets originating from the charged Higgs boson is required to satisfy a charm tagging selection. No significant excess beyond standard model predictions is found in the dijet invariant mass distribution. An upper limit in the range 1.68-0.25% is set on the branching fraction of the top quark decay to the charged Higgs boson and bottom quark for a charged Higgs boson mass between 80 and 160 GeV.

The production of Z boson pairs in proton-proton (pp) collisions, pp $\to$ (Z/$\gamma^*$)(Z/$\gamma^*$) $\to$ 2$\ell$ 2$\ell'$, where $\ell,\ell'$ = e or $\mu$, is studied at a center-of-mass energy of 13 TeV with the CMS detector at the CERN LHC. The data sample corresponds to an integrated luminosity of 137 fb$^{-1}$, collected during 2016-2018. The ZZ production cross section, $\sigma_{\text{tot}}$ (pp $\to$ ZZ) = 17.2 $\pm$ 0.3 (stat) $\pm$ 0.5 (syst) $\pm$ 0.4 (theo) $\pm$ 0.3 (lumi) pb, measured for events with two pairs of opposite-sign, same-flavor leptons produced in the mass region 60 $\lt$ $m_{\ell^+\ell^-}$ $\lt $ 120 GeV is consistent with standard model predictions. Differential cross sections are also measured and agree with theoretical predictions. The invariant mass distribution of the four-lepton system is used to set limits on anomalous ZZZ and ZZ$\gamma$ couplings.

The production cross sections for prompt open-charm mesons in proton-proton collisions at a center-of-mass energy of 13 TeV are reported. The measurement is performed using a data sample collected by the CMS experiment corresponding to an integrated luminosity of 29 nb$^{-1}$. The differential production cross sections of the D$^{*\pm}$, D$^\pm$, and D$^0$ ($\overline{\mathrm{D}}^{0}$) mesons are presented in ranges of transverse momentum and pseudorapidity 4 $\lt$$p_\mathrm{T}$$\lt$ 100 GeV and $\lvert\eta\rvert$$\lt$ 2.1, respectively. The results are compared to several theoretical calculations and to previous measurements.

Measurements of the associated production of a W boson and a charm (c) quark in proton-proton collisions at a centre-of-mass energy of 8 TeV are reported. The analysis uses a data sample corresponding to a total integrated luminosity of 19.7 fb$^{-1}$ collected by the CMS detector at the LHC. The W bosons are identified through their leptonic decays to an electron or a muon, and a neutrino. Charm quark jets are selected using distinctive signatures of charm hadron decays. The product of the cross section and branching fraction $\sigma$(pp $\to$ W + c + X) $\mathcal{B}$(W $\to$$\ell\nu$), where $\ell$ = e or $\mu$, and the cross section ratio $\sigma$(pp $\to$ W$^+$ + c + X) / $\sigma$(pp $\to$ W$^-$ + $\mathrm{\bar{c}}$ + X) are measured inclusively and differentially as functions of the pseudorapidity and of the transverse momentum of the lepton from the W boson decay. The results are compared with theoretical predictions. The impact of these measurements on the determination of the strange quark distribution is assessed.

A first observation is presented for the electroweak production of a W boson, a photon, and two jets in proton-proton collisions. The W boson decays are selected by requiring one identified electron or muon and an imbalance in transverse momentum. The two jets are required to have a high dijet mass and a large separation in pseudorapidity. The measurement is based on data collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The observed (expected) significance for this process is 4.9 (4.6) standard deviations. After combining with previously reported CMS results at 8 TeV, the observed (expected) significance is 5.3 (4.8) standard deviations. The cross section for the electroweak W$\gamma$jj production in a restricted fiducial region is measured as 20.4 $\pm$ 4.5 fb and the total cross section for W$\gamma$ production in association with 2 jets in the same fiducial region is 108 $\pm$ 16 fb. All results are in good agreement with recent theoretical predictions. Constraints are placed on anomalous quartic gauge couplings in terms of dimension-8 effective field theory operators.

A search for new physics in events with two highly Lorentz-boosted Z bosons and large missing transverse momentum is presented. The analyzed proton-proton collision data, corresponding to an integrated luminosity of 137 fb$^{-1}$, were recorded at $\sqrt{s} =$ 13 TeV by the CMS experiment at the CERN LHC. The search utilizes the substructure of jets with large radius to identify quark pairs from Z boson decays. Backgrounds from standard model processes are suppressed by requirements on the jet mass and the missing transverse momentum. No significant excess in the event yield is observed beyond the number of background events expected from the standard model. For a simplified supersymmetric model in which the Z bosons arise from the decay of gluinos, an exclusion limit of 1920 GeV on the gluino mass is set at 95% confidence level. This is the first search for beyond-standard-model production of pairs of boosted Z bosons plus large missing transverse momentum.

Evidence is presented for the electroweak (EW) production of two jets (jj) in association with two Z bosons and constraints on anomalous quartic gauge couplings are set. The analysis is based on a data sample of proton-proton collisions at $\sqrt{s} = $ 13 TeV collected with the CMS detector in 2016-2018, and corresponding to an integrated luminosity of 137 fb$^{-1}$. The search is performed in the fully leptonic final state ZZ $\to$ $\ell\ell\ell'\ell'$, where $\ell,\ell' = $ e, $\mu$. The EW production of two jets in association with two Z bosons is measured with an observed (expected) significance of 4.0 (3.5) standard deviations. The cross sections for the EW production are measured in three fiducial volumes and the result is $\sigma_{\mathrm{EW}}$(pp $\to$ ZZjj $\to$ $\ell\ell\ell'\ell'$jj) = 0.33 $^{+0.11}_{-0.10}$ (stat) $^{+0.04}_{-0.03}$ (syst) fb in the most inclusive volume, in agreement with the standard model prediction of 0.275 $\pm$ 0.021 fb. Measurements of total cross sections for jj production in association with two Z bosons are also reported. Limits on anomalous quartic gauge couplings are derived in terms of the effective field theory operators T0, T1, T2, T8, and T9.

Collinear (small-angle) and large-angle, as well as soft and hard radiations are investigated in three-jet and Z + two-jet events collected in proton-proton collisions at the LHC. The normalized production cross sections are measured as a function of the ratio of transverse momenta of two jets and their angular separation. The measurements in the three-jet and Z + two-jet events are based on data collected at a center-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.8 fb$^{-1}$. The Z + two-jet events are reconstructed in the dimuon decay channel of the Z boson. The three-jet measurement is extended to include $\sqrt{s} =$ 13 TeV data corresponding to an integrated luminosity of 2.3 fb$^{-1}$. The results are compared to predictions from event generators that include parton showers, multiple parton interactions, and hadronization. The collinear and soft regions are in general well described by parton showers, whereas the regions of large angular separation are often best described by calculations using higher-order matrix elements.

Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum $p_\mathrm{T}$ and of a geometric variable $\phi^*$ are measured, where $\phi^*$ highly correlates with $p_\mathrm{T}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$$\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$$\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

The first measurement of the dependence of $\gamma\gamma$$\to$$\mu^{+}\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\lt$$m_{\mu\mu}$$\lt$ 60 GeV are extracted for events including 0, 1, or at least 2 neutrons detected in the forward pseudorapidity range $|\eta|$$\gt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.

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

Jets containing a prompt J$/\psi$ meson are studied in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV, using the CMS detector at the LHC. Jets are selected to be in the transverse momentum range of 30 $\lt$$p_\mathrm{T}$$\lt$ 40 GeV. The J$/\psi$ yield in these jets is evaluated as a function of the jet fragmentation variable $z$, the ratio of the J$/\psi$$p_\mathrm{T} $ to the jet $p_\mathrm{T}$. The nuclear modification factor, $R_\mathrm{AA}$, is then derived by comparing the yield in lead-lead collisions to the corresponding expectation based on proton-proton data, at the same nucleon-nucleon center-of-mass energy. The suppression of the J$/\psi$ yield shows a dependence on $z$, indicating that the interaction of the J$/\psi$ with the quark-gluon plasma formed in heavy ion collisions depends on the fragmentation that gives rise to the J$/\psi$ meson.

A search is presented for a heavy vector resonance decaying into a Z boson and the standard model Higgs boson, where the Z boson is identified through its leptonic decays to electrons, muons, or neutrinos, and the Higgs boson is identified through its hadronic decays. The search is performed in a Lorentz-boosted regime and is based on data collected from 2016 to 2018 at the CERN LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. Upper limits are derived on the production of a narrow heavy resonance Z', and a mass below 3.5 and 3.7 TeV is excluded at 95% confidence level in models where the heavy vector boson couples exclusively to fermions and to bosons, respectively. These are the most stringent limits placed on the Heavy Vector Triplet Z' model to date. If the heavy vector boson couples exclusively to standard model bosons, upper limits on the product of the cross section and branching fraction are set between 23 and 0.3 fb for a Z' mass between 0.8 and 4.6 TeV, respectively. This is the first limit set on a heavy vector boson coupling exclusively to standard model bosons in its production and decay.

The first measurements of production cross sections of polarized same-sign W$^\pm$W$^\pm$ boson pairs in proton-proton collisions are reported. The measurements are based on a data sample collected with the CMS detector at the LHC at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 137 fb$^{-1}$. Events are selected by requiring exactly two same-sign leptons, electrons or muons, moderate missing transverse momentum, and two jets with a large rapidity separation and a large dijet mass to enhance the contribution of same-sign W$^\pm$W$^\pm$ scattering events. An observed (expected) 95% confidence level upper limit of 1.17 (0.88) fb is set on the production cross section for longitudinally polarized same-sign W$^\pm$W$^\pm$ boson pairs. The electroweak production of same-sign W$^\pm$W$^\pm$ boson pairs with at least one of the W bosons longitudinally polarized is measured with an observed (expected) significance of 2.3 (3.1) standard deviations.

A search is performed for the rare decay W$^\pm\to\pi^\pm\gamma$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV. Data corresponding to an on W integrated luminosity of 137 fb$^{-1}$ were collected during 2016 to 2018 with the CMS detector. This analysis exploits a novel search strategy based on W boson production in top quark pair events. An inclusive search for the W$^\pm\to\pi^\pm\gamma$ decay is not optimal at the LHC because of the high trigger thresholds. Instead, a trigger selection is exploited in which the W boson originating from one of the top quarks is used to tag the event in a leptonic decay. The W boson emerging from the other top quark is used to search for the W$^\pm\to\pi^\pm\gamma$ signature. Such decays are characterized by an isolated track pointing to a large energy deposit, and by an isolated photon of large transverse momentum. The presence of b quark jets reduces the background from the hadronization of light-flavor quarks and gluons. The W$^\pm\to\pi^\pm\gamma$ decay is not observed. An upper exclusion limit is set to this branching fraction, corresponding to 1.50 $\times$ 10$^{-5}$ at 95% confidence level, whereas the expected upper limit exclusion limit is 0.85 $^{+0.52}_{-0.29}$ $\times$ 10$^{-5}$.

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb$^{-1}$, collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses.

The rate for Higgs (H) bosons production in association with either one (tH) or two ($\mathrm{t\bar{t}}$H) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton-proton collisions recorded at a center-of-mass energy of 13 TeV by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137 fb$^{-1}$. The analysis is aimed at events that contain H $\to$ WW, H $\to$$\tau\tau$, or H $\to$ ZZ decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among the tH, the $\mathrm{t\bar{t}}$H, and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the $\mathrm{t\bar{t}}$H and tH signals correspond to 0.92 $\pm$ 0.19 (stat) $^{+0.17}_{-0.13}$ (syst) and 5.7 $\pm$ 2.7 (stat) $\pm$ 3.0 (syst) of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for $\mathrm{t\bar{t}}$H, and to 1.4 (0.3) for tH production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $y_{\mathrm{t}}$ of the Higgs boson to the top quark divided by its SM expectation, $\kappa_\mathrm{t}$ = $y_\mathrm{t} / y_\mathrm{t}^\mathrm{SM}$, is constrained to be within $-$0.9 $\lt$$\kappa_\mathrm{t}$$\lt$$-$0.7 or 0.7 $\lt$$\kappa_\mathrm{t}$$\lt$ 1.1, at 95% confidence level. This result is the most sensitive measurement of the $\mathrm{t\bar{t}}$H production rate to date.

A search for exclusive two-photon production via photon exchange in proton-proton collisions, pp $\to$ p$\gamma\gamma$p with intact protons, is presented. The data correspond to an integrated luminosity of 9.4 fb$^{-1}$ collected in 2016 using the CMS and TOTEM detectors at a center-of-mass energy of 13 TeV at the LHC. Events are selected with a diphoton invariant mass above 350 GeV and with both protons intact in the final state, to reduce backgrounds from strong interactions. The events of interest are those where the invariant mass and rapidity calculated from the momentum losses of the forward-moving protons matches the mass and rapidity of the central, two-photon system. No events are found that satisfy this condition. Interpreting this result in an effective dimension-8 extension of the standard model, the first limits are set on the two anomalous four-photon coupling parameters. If the other parameter is constrained to its standard model value, the limits at 95% CL are $\lvert\zeta_1\rvert$ $\lt$ 2.9 $\times$ 10$^{-13}$ GeV$^{-4}$ and $\lvert\zeta_2\rvert$ $\lt$ 6.0 $\times$ 10$^{-13}$ GeV$^{-4}$.

Evidence for Higgs boson decay to a pair of muons is presented. This result combines searches in four exclusive categories targeting the production of the Higgs boson via gluon fusion, via vector boson fusion, in association with a vector boson, and in association with a top quark-antiquark pair. The analysis is performed using proton-proton collision data at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 137 fb$^{-1}$, recorded by the CMS experiment at the CERN LHC. An excess of events over the background expectation is observed in data with a significance of 3.0 standard deviations, where the expectation for the standard model (SM) Higgs boson with mass of 125.38 GeV is 2.5. The combination of this result with that from data recorded at $\sqrt{s} =$ 7 and 8 TeV, corresponding to integrated luminosities of 5.1 and 19.7 fb$^{-1}$, respectively, increases both the expected and observed significances by 1%. The measured signal strength, relative to the SM prediction, is 1.19 $^{+0.40}_{-0.39}$ (stat) $^{+0.15}_{-0.14}$ (syst). This result constitutes the first evidence for the decay of the Higgs boson to second generation fermions and is the most precise measurement of the Higgs boson coupling to muons reported to date.

A search is described for the production of a pair of bottom-type vector-like quarks (VLQs), each decaying into a b or $\mathrm{\bar{b}}$ quark and either a Higgs or a Z boson, with a mass greater than 1000 GeV. The analysis is based on data from proton-proton collisions at a 13 TeV center-of-mass energy recorded at the CERN LHC, corresponding to a total integrated luminosity of 137 fb$^{-1}$. As the predominant decay modes of the Higgs and Z bosons are to a pair of quarks, the analysis focuses on final states consisting of jets resulting from the six quarks produced in the events. Since the two jets produced in the decay of a highly Lorentz-boosted Higgs or Z boson can merge to form a single jet, nine independent analyses are performed, categorized by the number of observed jets and the reconstructed event mode. No signal in excess of the expected background is observed. Lower limits are set on the VLQ mass at 95% confidence level equal to 1570 GeV in the case where the VLQ decays exclusively to a b quark and a Higgs boson, 1390 GeV for when it decays exclusively to a b quark and a Z boson, and 1450 GeV for when it decays equally in these two modes. These limits represent significant improvements over the previously published VLQ limits.

A search for leptoquarks produced singly and in pairs in proton-proton collisions is presented. We consider the leptoquark (LQ) to be a scalar particle of charge -1/3$e$ coupling to a top quark plus a tau lepton ($\mathrm{t}\tau$) or a bottom quark plus a neutrino ($\mathrm{b}\nu$), or a vector particle of charge +2/3$e$, coupling to $\mathrm{t}\nu$ or $\mathrm{b}\tau$. These choices are motivated by models that can explain a series of anomalies observed in the measurement of B meson decays. In this analysis the signatures $\mathrm{t}\tau\nu\mathrm{b}$ and $\mathrm{t}\tau\nu$ are probed, using data recorded by the CMS experiment at the CERN LHC at $\sqrt{s} =$ 13 TeV and that correspond to an integrated luminosity of 137 fb$^{-1}$. These signatures have not been previously explored in a dedicated search. The data are found to be in agreement with the standard model prediction. Lower limits at 95% confidence level are set on the LQ mass in the range 0.98-1.73 TeV, depending on the LQ spin and its coupling $\lambda$ to a lepton and a quark, and assuming equal branching fractions for the two LQ decay modes considered. These are the most stringent constraints to date on the existence of leptoquarks in this scenario.