Cross-sections for the production of a $Z$ boson in association with two photons are measured in proton$-$proton collisions at a centre-of-mass energy of 13 TeV. The data used correspond to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS experiment during Run 2 of the LHC. The measurements use the electron and muon decay channels of the $Z$ boson, and a fiducial phase-space region where the photons are not radiated from the leptons. The integrated $Z(\rightarrow\ell\ell)\gamma\gamma$ cross-section is measured with a precision of 12% and differential cross-sections are measured as a function of six kinematic variables of the $Z\gamma\gamma$ system. The data are compared with predictions from MC event generators which are accurate to up to next-to-leading order in QCD. The cross-section measurements are used to set limits on the coupling strengths of dimension-8 operators in the framework of an effective field theory.

A search for physics beyond the standard model (SM) in the final state with a hadronically decaying tau lepton and a neutrino is presented. This analysis is based on data recorded by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC, corresponding to a total integrated luminosity of 138 fb$^{=1}$. The transverse mass spectrum is analyzed for the presence of new physics. No significant deviation from the SM prediction is observed. Limits are set on the production cross section of a W' boson decaying into a tau lepton and a neutrino. Lower limits are set on the mass of the sequential SM-like heavy charged vector boson and the mass of a quantum black hole. Upper limits are placed on the couplings of a new boson to the SM fermions. Constraints are put on a nonuniversal gauge interaction model and an effective field theory model. For the first time, upper limits on the cross section of $t$-channel leptoquark (LQ) exchange are presented. These limits are translated into exclusion limits on the LQ mass and on its coupling in the $t$-channel. The sensitivity of this analysis extends into the parameter space of LQ models that attempt to explain the anomalies observed in B meson decays. The limits presented for the various interpretations are the most stringent to date. Additionally, a model-independent limit is provided.

Proton-proton interactions resulting in final states with two photons are studied in a search for the signature of flavor-changing neutral current interactions of top quarks (t) and Higgs bosons (H). The analysis is based on data collected at a center-of-mass energy of 13 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. No significant excess above the background prediction is observed. Upper limits on the branching fractions ($\mathcal{B}$) of the top quark decaying to a Higgs boson and an up (u) or charm quark (c) are derived through a binned fit to the diphoton invariant mass spectrum. The observed (expected) 95% confidence level upper limits are found to be 0.019 (0.031)% for $\mathcal B$(t $\to$ Hu) and 0.073 (0.051)% for $\mathcal{B}$(t $\to$ Hc). These are the strictest upper limits yet determined.

A search is presented for displaced production of Higgs bosons or $Z$ bosons, originating from the decay of a neutral long-lived particle (LLP) and reconstructed in the decay modes $H\rightarrow \gamma\gamma$ and $Z\rightarrow ee$. The analysis uses the full Run 2 data set of proton$-$proton collisions delivered by the LHC at an energy of $\sqrt{s}=13$ TeV between 2015 and 2018 and recorded by the ATLAS detector, corresponding to an integrated luminosity of 139 fb$^{-1}$. Exploiting the capabilities of the ATLAS liquid argon calorimeter to precisely measure the arrival times and trajectories of electromagnetic objects, the analysis searches for the signature of pairs of photons or electrons which arise from a common displaced vertex and which arrive after some delay at the calorimeter. The results are interpreted in a gauge-mediated supersymmetry breaking model with pair-produced higgsinos that decay to LLPs, and each LLP subsequently decays into either a Higgs boson or a $Z$ boson. The final state includes at least two particles that escape direct detection, giving rise to missing transverse momentum. No significant excess is observed above the background expectation. The results are used to set upper limits on the cross section for higgsino pair production, up to a $\tilde\chi^0_1$ mass of 369 (704) GeV for decays with 100% branching ratio of $\tilde\chi^0_1$ to Higgs ($Z$) bosons for a $\tilde\chi^0_1$ lifetime of 2 ns. A model-independent limit is also set on the production of pairs of photons or electrons with a significant delay in arrival at the calorimeter.

A search for new phenomena has been performed in final states with at least one isolated high-momentum photon, jets and missing transverse momentum in proton--proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV. The data, collected by the ATLAS experiment at the CERN LHC, correspond to an integrated luminosity of 139 $fb^{-1}$. The experimental results are interpreted in a supersymmetric model in which pair-produced gluinos decay into neutralinos, which in turn decay into a gravitino, at least one photon, and jets. No significant deviations from the predictions of the Standard Model are observed. Upper limits are set on the visible cross section due to physics beyond the Standard Model, and lower limits are set on the masses of the gluinos and neutralinos, all at 95% confidence level. Visible cross sections greater than 0.022 fb are excluded and pair-produced gluinos with masses up to 2200 GeV are excluded for most of the NLSP masses investigated.

A search for diphoton resonances in the mass range between 10 and 70 GeV with the ATLAS experiment at the Large Hadron Collider (LHC) is presented. The analysis is based on $pp$ collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at a centre-of-mass energy of 13 TeV recorded from 2015 to 2018. Previous searches for diphoton resonances at the LHC have explored masses down to 65 GeV, finding no evidence of new particles. This search exploits the particular kinematics of events with pairs of closely spaced photons reconstructed in the detector, allowing examination of invariant masses down to 10 GeV. The presented strategy covers a region previously unexplored at hadron colliders because of the experimental challenges of recording low-energy photons and estimating the backgrounds. No significant excess is observed and the reported limits provide the strongest bound on promptly decaying axion-like particles coupling to gluons and photons for masses between 10 and 70 GeV.

A search for long-lived particles decaying into hadrons is presented. The analysis uses 139 fb$^{-1}$ of $pp$ collision data collected at $\sqrt{s} = 13$ TeV by the ATLAS detector at the LHC using events that contain multiple energetic jets and a displaced vertex. The search employs dedicated reconstruction techniques that significantly increase the sensitivity to long-lived particles decaying in the ATLAS inner detector. Background estimates for Standard Model processes and instrumental effects are extracted from data. The observed event yields are compatible with those expected from background processes. The results are used to set limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model, and on scenarios with pair-production of supersymmetric particles with long-lived electroweakinos that decay via a small $R$-parity-violating coupling. The pair-production of electroweakinos with masses below 1.5 TeV is excluded for mean proper lifetimes in the range from 0.03 ns to 1 ns. When produced in the decay of $m(\tilde{g})=2.4$ TeV gluinos, electroweakinos with $m(\tilde\chi^0_1)=1.5$ TeV are excluded with lifetimes in the range of 0.02 ns to 4 ns.

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.

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 chargino$-$neutralino pair production in three-lepton final states with missing transverse momentum is presented. The study is based on a dataset of $\sqrt{s} = 13$ TeV $pp$ collisions recorded with the ATLAS detector at the LHC, corresponding to an integrated luminosity of 139 fb$^{-1}$. No significant excess relative to the Standard Model predictions is found in data. The results are interpreted in simplified models of supersymmetry, and statistically combined with results from a previous ATLAS search for compressed spectra in two-lepton final states. Various scenarios for the production and decay of charginos ($\tilde\chi^\pm_1$) and neutralinos ($\tilde\chi^0_2$) are considered. For pure higgsino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair-production scenarios, exclusion limits at 95% confidence level are set on $\tilde\chi^0_2$ masses up to 210 GeV. Limits are also set for pure wino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair production, on $\tilde\chi^0_2$ masses up to 640 GeV for decays via on-shell $W$ and $Z$ bosons, up to 300 GeV for decays via off-shell $W$ and $Z$ bosons, and up to 190 GeV for decays via $W$ and Standard Model Higgs bosons.

Searches are performed for nonresonant and resonant di-Higgs boson production in the $b\bar{b}\gamma\gamma$ final state. The data set used corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. No excess above the expected background is found and upper limits on the di-Higgs boson production cross sections are set. A 95% confidence-level upper limit of 4.2 times the cross section predicted by the Standard Model is set on $pp \rightarrow HH$ nonresonant production, where the expected limit is 5.7 times the Standard Model predicted value. The expected constraints are obtained for a background hypothesis excluding $pp \rightarrow HH$ production. The observed (expected) constraints on the Higgs boson trilinear coupling modifier $\kappa_{\lambda}$ are determined to be $[-1.5, 6.7]$ $([-2.4, 7.7])$ at 95% confidence level, where the expected constraints on $\kappa_{\lambda}$ are obtained excluding $pp \rightarrow HH$ production from the background hypothesis. For resonant production of a new hypothetical scalar particle $X$ ($X \rightarrow HH \rightarrow b\bar{b}\gamma\gamma$), limits on the cross section for $pp \to X \to HH$ are presented in the narrow-width approximation as a function of $m_{X}$ in the range $251 \leq m_{X} \leq 1000$ GeV. The observed (expected) limits on the cross section for $pp \to X \to HH$ range from 640 fb to 44 fb (391 fb to 46 fb) over the considered mass range.

A novel search for exotic decays of the Higgs boson into pairs of long-lived neutral particles, each decaying into a bottom quark pair, is performed using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collision data collected with the ATLAS detector at the LHC. Events consistent with the production of a Higgs boson in association with a leptonically decaying $Z$ boson are analysed. Long-lived particle (LLP) decays are reconstructed from inner-detector tracks as displaced vertices with high mass and track multiplicity relative to Standard Model processes. The analysis selection requires the presence of at least two displaced vertices, effectively suppressing Standard Model backgrounds. The residual background contribution is estimated using a data-driven technique. No excess over Standard Model predictions is observed, and upper limits are set on the branching ratio of the Higgs boson to LLPs. Branching ratios above 10% are excluded at 95% confidence level for LLP mean proper lifetimes $c\tau$ as small as 4 mm and as large as 100 mm. For LLP masses below 40 GeV, these results represent the most stringent constraint in this lifetime regime.

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.

A search for long-lived particles decaying into jets is presented. Data were collected with the CMS detector at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search examines the distinctive topology of displaced tracks and secondary vertices. The selected events are found to be consistent with standard model predictions. For a simplified model in which long-lived neutral particles are pair produced and decay to two jets, pair production cross sections larger than 0.2 fb are excluded at 95% confidence level for a long-lived particle mass larger than 1000 GeV and proper decay lengths between 3 and 130 mm. Several supersymmetry models with gauge-mediated supersymmetry breaking or $R$-parity violation, where pair-produced long-lived gluinos or top squarks decay to several final-state topologies containing displaced jets, are also tested. For these models, in the mass ranges above 200 GeV, gluino masses up to 2300-2400 GeV and top squark masses up to 1350-1600 GeV are excluded for proper decay lengths approximately between 10 and 100 mm. These are the most restrictive limits to date on these models.

A search for the electroweak production of pairs of charged sleptons or charginos decaying into two-lepton final states with missing transverse momentum is presented. Two simplified models of $R$-parity-conserving supersymmetry are considered: direct pair-production of sleptons ($\tilde{\ell}\tilde{\ell}$), with each decaying into a charged lepton and a $\tilde{\chi}_1^0$ neutralino, and direct pair-production of the lightest charginos $(\tilde{\chi}_1^\pm\tilde{\chi}_1^\mp)$, with each decaying into a $W$-boson and a $\tilde{\chi}_1^0$. The lightest neutralino ($\tilde{\chi}_1^0$) is assumed to be the lightest supersymmetric particle (LSP). The analyses target the experimentally challenging mass regions where $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and $m(\tilde{\chi}_1^\pm)-m(\tilde{\chi}_1^0)$ are close to the $W$-boson mass (`moderately compressed' regions). The search uses 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider. No significant excesses over the expected background are observed. Exclusion limits on the simplified models under study are reported in the ($\tilde{\ell},\tilde{\chi}_1^0$) and ($\tilde{\chi}_1^\pm,\tilde{\chi}_1^0$) mass planes at 95% confidence level (CL). Sleptons with masses up to 150 GeV are excluded at 95% CL for the case of a mass-splitting between sleptons and the LSP of 50 GeV. Chargino masses up to 140 GeV are excluded at 95% CL for the case of a mass-splitting between the chargino and the LSP down to about 100 GeV.

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

A search for long-lived particles (LLPs) produced in decays of standard model (SM) Higgs bosons is presented. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, recorded at the LHC in 2016-2018. A novel technique is employed to reconstruct decays of LLPs in the endcap muon detectors. The search is sensitive to a broad range of LLP decay modes and to masses as low as a few GeV. No excess of events above the SM background is observed. The most stringent limits to date on the branching fraction of the Higgs boson to LLPs subsequently decaying to quarks and $\tau^+\tau^-$ are found for proper decay lengths greater than 6, 20, and 40 m, for LLP masses of 7, 15, and 40 GeV, respectively.

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

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.

Results are reported from a search for supersymmetric particles in the final state with multiple jets and large missing transverse momentum. The search uses a sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected with the CMS detector in 2016-2018, corresponding to an integrated luminosity of 137 fb$^{-1}$, representing essentially the full LHC Run 2 data sample. The analysis is performed in a four-dimensional search region defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No significant excess in the event yield is observed relative to the expected background contributions from standard model processes. Limits on the pair production of gluinos and squarks are obtained in the framework of simplified models for supersymmetric particle production and decay processes. Assuming the lightest supersymmetric particle to be a neutralino, lower limits on the gluino mass as large as 2000 to 2310 GeV are obtained at 95% confidence level, while lower limits on the squark mass as large as 1190 to 1630 GeV are obtained, depending on the production scenario.

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.

We present a search for magnetic monopoles and high-electric-charge objects using LHC Run 2 $\sqrt{s} =$13 TeV proton$-$proton collisions recorded by the ATLAS detector. A total integrated luminosity of 138 fb$^{-1}$ was collected by a specialized trigger. No highly ionizing particle candidate was observed. Considering the Drell-Yan and photon-fusion pair production mechanisms as benchmark models, cross-section upper limits are presented for spin-0 and spin-$\frac{1}{2}$ magnetic monopoles of magnetic charge $1g_\textrm{D}$ and $2g_\textrm{D}$ and for high-electric-charge objects of electric charge $20 \leq |z| \leq 100$, for masses between 200 GeV and 4000 GeV. The search improves by approximately a factor of three the previous cross-section limits on the Drell-Yan production of magnetic monopoles and high-electric charge objects. Also, the first ATLAS limits on the photon-fusion pair production mechanism of magnetic monopoles and high-electric-charge objects have been obtained.

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 for supersymmetry in events with four or more charged leptons (electrons, muons and $\tau$-leptons) is presented. The analysis uses a data sample corresponding to $139\,\mbox{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 $\tau$-leptons are designed to target several supersymmetric models, while a general five-lepton signal region targets any new physics phenomena leading to a final state with five charged leptons. Data yields are consistent with Standard Model expectations and results are used to set upper limits on contributions from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of general gauge-mediated supersymmetry, excluding higgsino masses up to $540$ GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of $1.6$ TeV, $1.2$ TeV, and $2.5$ TeV are placed on wino, slepton and gluino masses, respectively.

A search is presented for new particles produced at the LHC in proton-proton collisions at $\sqrt{s} =$ 13 TeV, using events with energetic jets and large missing transverse momentum. The analysis is based on a data sample corresponding to an integrated luminosity of 101 fb$^{-1}$, collected in 2017-2018 with the CMS detector. Machine learning techniques are used to define separate categories for events with narrow jets from initial-state radiation and events with large-radius jets consistent with a hadronic decay of a W or Z boson. A statistical combination is made with an earlier search based on a data sample of 36 fb$^{-1}$, collected in 2016. No significant excess of events is observed with respect to the standard model background expectation determined from control samples in data. The results are interpreted in terms of limits on the branching fraction of an invisible decay of the Higgs boson, as well as constraints on simplified models of dark matter, on first-generation scalar leptoquarks decaying to quarks and neutrinos, and on models with large extra dimensions. Several of the new limits, specifically for spin-1 dark matter mediators, pseudoscalar mediators, colored mediators, and leptoquarks, are the most restrictive to date.