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.

A search for nonresonant Higgs boson pair production in the $b\bar{b}\gamma\gamma$ final state is performed using 140 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. This analysis supersedes and expands upon the previous nonresonant ATLAS results in this final state based on the same data sample. The analysis strategy is optimised to probe anomalous values not only of the Higgs ($H$) boson self-coupling modifier $\kappa_\lambda$ but also of the quartic $HHVV$ ($V=W,Z$) coupling modifier $\kappa_{2V}$. No significant excess above the expected background from Standard Model processes is observed. An observed upper limit $\mu_{HH}<4.0$ is set at 95% confidence level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The 95% confidence intervals for the coupling modifiers are $-1.4<\kappa_\lambda<6.9$ and $-0.5<\kappa_{2V}<2.7$, assuming all other Higgs boson couplings except the one under study are fixed to the Standard Model predictions. The results are interpreted in the Standard Model effective field theory and Higgs effective field theory frameworks in terms of constraints on the couplings of anomalous Higgs boson (self-)interactions.

Statistical combinations of searches for charginos and neutralinos using various decay channels are performed using $139\,$fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13\,$TeV with the ATLAS detector at the Large Hadron Collider. Searches targeting pure-wino chargino pair production, pure-wino chargino-neutralino production, or higgsino production decaying via Standard Model $W$, $Z$, or $h$ bosons are combined to extend the mass reach to the produced SUSY particles by 30-100 GeV. The depth of the sensitivity of the original searches is also improved by the combinations, lowering the 95% CL cross-section upper limits by 15%-40%.

A search for physics beyond the Standard Model inducing periodic signals in the dielectron and diphoton invariant mass spectra is presented using 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data collected by the ATLAS experiment at the LHC. Novel search techniques based on continuous wavelet transforms are used to infer the frequency of periodic signals from the invariant mass spectra and neural network classifiers are used to enhance the sensitivity to periodic resonances. In the absence of a signal, exclusion limits are placed at the 95% confidence level in the two-dimensional parameter space of the clockwork gravity model. Model-independent searches for deviations from the background-only hypothesis are also performed.

This paper presents a study of $Z \to ll\gamma~$decays with the ATLAS detector at the Large Hadron Collider. The analysis uses a proton-proton data sample corresponding to an integrated luminosity of 20.2 fb$^{-1}$ collected at a centre-of-mass energy $\sqrt{s}$ = 8 TeV. Integrated fiducial cross-sections together with normalised differential fiducial cross-sections, sensitive to the kinematics of final-state QED radiation, are obtained. The results are found to be in agreement with state-of-the-art predictions for final-state QED radiation. First measurements of $Z \to ll\gamma\gamma$ decays are also reported.

A search for high-mass charged and neutral bosons decaying to $W\gamma$ and $Z\gamma$ final states is presented in this paper. The analysis uses a data sample of $\sqrt{s} = 13$ TeV proton-proton collisions with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector during LHC Run 2 operation. The sensitivity of the search is determined using models of the production and decay of spin-1 charged bosons and spin-0/2 neutral bosons. The range of resonance masses explored extends from 1.0 TeV to 6.8 TeV. At these high resonance masses, it is beneficial to target the hadronic decays of the $W$ and $Z$ bosons because of their large branching fractions. The decay products of the high-momentum $W/Z$ bosons are strongly collimated and boosted-boson tagging techniques are employed to improve the sensitivity. No evidence of a signal above the Standard Model backgrounds is observed, and upper limits on the production cross-sections of these bosons times their branching fractions to $W\gamma$ and $Z\gamma$ are derived for various boson production models.

Searches for new resonances are performed using an unsupervised anomaly-detection technique. Events with at least one electron or muon are selected from 140 fb$^{-1}$ of $pp$ collisions at $\sqrt{s} = 13$ TeV recorded by ATLAS at the Large Hadron Collider. The approach involves training an autoencoder on data, and subsequently defining anomalous regions based on the reconstruction loss of the decoder. Studies focus on nine invariant mass spectra that contain pairs of objects consisting of one light jet or $b$-jet and either one lepton ($e$, $\mu$), photon, or second light jet or $b$-jet in the anomalous regions. No significant deviations from the background hypotheses are observed.

A search for light long-lived neutral particles with masses in the $O$(MeV-GeV) range is presented. The analysis targets the production of long-lived dark photons in the decay of a Higgs boson produced via gluon-gluon fusion or in association with a $W$ boson. Events that contain displaced collimated Standard Model fermions reconstructed in the calorimeter or muon spectrometer are selected in 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV $pp$ collision data collected by the ATLAS detector at the LHC. Background estimates for contributions from Standard Model processes and instrumental effects are extracted from data. The observed event yields are consistent with the expected background. Exclusion limits are reported on the production cross-section times branching fraction as a function of the mean proper decay length $c\tau$ of the dark photon, or as a function of the dark-photon mass and kinetic mixing parameter that quantifies the coupling between the Standard Model and potential hidden (dark) sectors. A Higgs boson branching fraction above 1% is excluded at 95% CL for a Higgs boson decaying into two dark photons for dark-photon mean proper decay lengths between 10 mm and 250 mm and dark photons with masses between 0.4 GeV and 2 GeV.

A search for forward proton scattering in association with light-by-light scattering mediated by an axion-like particle is presented, using the ATLAS Forward Proton spectrometer to detect scattered protons and the central ATLAS detector to detect pairs of outgoing photons. Proton-proton collision data recorded in 2017 at a centre-of-mass energy of $\sqrt{s} = 13$ TeV were analysed, corresponding to an integrated luminosity of 14.6 fb$^{-1}$. A total of 441 candidate signal events were selected. A search was made for a narrow resonance in the diphoton mass distribution, corresponding to an axion-like particle (ALP) with mass in the range 150-1600 GeV. No excess is observed above a smooth background. Upper limits on the production cross section of a narrow resonance are set as a function of the mass, and are interpreted as upper limits on the ALP production coupling constant, assuming 100% decay branching ratio into a photon pair. The inferred upper limit on the coupling constant is in the range 0.04-0.09 TeV$^{-1}$ at 95%confidence level.

A search for flavor-changing neutral-current couplings between a top quark, an up or charm quark and a $Z$ boson is presented, using proton-proton collision data at $\sqrt{s} = 13$ TeV collected by the ATLAS detector at the Large Hadron Collider. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets both single-top-quark events produced as $gq\rightarrow tZ$ (with $q = u, c$) and top-quark-pair events, with one top quark decaying through the $t \rightarrow Zq$ channel. The analysis considers events with three leptons (electrons or muons), a $b$-tagged jet, possible additional jets, and missing transverse momentum. The data are found to be consistent with the background-only hypothesis and 95% confidence-level limits on the $t \rightarrow Zq$ branching ratios are set, assuming only tensor operators of the Standard Model effective field theory framework contribute to the $tZq$ vertices. These are $6.2 \times 10^{-5}$ ($13\times 10^{-5}$) for $t\rightarrow Zu$ ($t\rightarrow Zc$) for a left-handed $tZq$ coupling, and $6.6 \times 10^{-5}$ ($12\times 10^{-5}$) in the case of a right-handed coupling. These results are interpreted as 95% CL upper limits on the strength of corresponding couplings, yielding limits for $|C_{uW}^{(13)*}|$ and $|C_{uB}^{(13)*}|$ ($|C_{uW}^{(31)}|$ and $|C_{uB}^{(31)}|$) of 0.15 (0.16), and limits for $|C_{uW}^{(23)*}|$ and $|C_{uB}^{(23)*}|$ ($|C_{uW}^{(32)}|$ and $|C_{uB}^{(32)}|$) of 0.22 (0.21), assuming a new-physics energy scale $\Lambda_\text{NP}$ of 1 TeV.

Differential cross-section measurements of $Z\gamma$ production in association with hadronic jets are presented, using the full 139 fb$^{-1}$ dataset of $\sqrt{s}=13$ TeV proton-proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO$_\text{PS}$, as well as next-to-leading-order predictions from MadGraph5_aMC@NLO and Sherpa.

Differential and double-differential distributions of kinematic variables of leptons from decays of top-quark pairs ($t\bar{t}$) are measured using the full LHC Run 2 data sample collected with the ATLAS detector. The data were collected at a $pp$ collision energy of $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 140 fb$^{-1}$. The measurements use events containing an oppositely charged $e\mu$ pair and $b$-tagged jets. The results are compared with predictions from several Monte Carlo generators. While no prediction is found to be consistent with all distributions, a better agreement with measurements of the lepton $p_{\text{T}}$ distributions is obtained by reweighting the $t\bar{t}$ sample so as to reproduce the top-quark $p_{\text{T}}$ distribution from an NNLO calculation. The inclusive top-quark pair production cross-section is measured as well, both in a fiducial region and in the full phase-space. The total inclusive cross-section is found to be \[ \sigma_{t\bar{t}} = 829 \pm 1\;(\textrm{stat}) \pm 13\;(\textrm{syst}) \pm 8\;(\textrm{lumi}) \pm 2\; (\textrm{beam})\ \textrm{pb}, \] where the uncertainties are due to statistics, systematic effects, the integrated luminosity and the beam energy. This is in excellent agreement with the theoretical expectation.

New particles with large masses that decay into hadronically interacting particles are predicted by many models of physics beyond the Standard Model. A search for a massive resonance that decays into pairs of dijet resonances is performed using 140 fb$^{-1}$ of proton$-$proton collisions at $\sqrt{s}=13$ TeV recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Resonances are searched for in the invariant mass of the tetrajet system, and in the average invariant mass of the pair of dijet systems. A data-driven background estimate is obtained by fitting the tetrajet and dijet invariant mass distributions with a four-parameter dijet function and a search for local excesses from resonant production of dijet pairs is performed. No significant excess of events beyond the Standard Model expectation is observed, and upper limits are set on the production cross-sections of new physics scenarios.

A search for heavy right-handed Majorana or Dirac neutrinos $N_{\mathrm{R}}$ and heavy right-handed gauge bosons $W_{\mathrm{R}}$ is performed in events with energetic electrons or muons, with the same or opposite electric charge, and energetic jets. The search is carried out separately for topologies of clearly separated final-state products (``resolved'' channel) and topologies with boosted final states with hadronic and/or leptonic products partially overlapping and reconstructed as a large-radius jet (``boosted'' channel). The events are selected from $pp$ collision data at the LHC with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector at $\sqrt{s}$ = 13 TeV. No significant deviations from the Standard Model predictions are observed. The results are interpreted within the theoretical framework of a left-right symmetric model, and lower limits are set on masses in the heavy right-handed $W_{\mathrm{R}}$ boson and $N_{\mathrm{R}}$ plane. The excluded region extends to about $m(W_{\mathrm{R}}) = 6.4$ TeV for both Majorana and Dirac $N_{\mathrm{R}}$ neutrinos at $m(N_{\mathrm{R}})<1$ TeV. $N_{\mathrm{R}}$ with masses of less than 3.5 (3.6) TeV are excluded in the electron (muon) channel at $m(W_{\mathrm{R}})=4.8$ TeV for the Majorana neutrinos, and limits of $m(N_{\mathrm{R}})$ up to 3.6 TeV for $m(W_{\mathrm{R}}) = 5.2$ (5.0) TeV in the electron (muon) channel are set for the Dirac neutrinos. These constitute the most stringent exclusion limits to date for the model considered.

A measurement of single top-quark production in the s-channel is performed in proton$-$proton collisions at a centre-of-mass energy of 13 TeV with the ATLAS detector at the CERN Large Hadron Collider. The dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The analysis is performed on events with an electron or muon, missing transverse momentum and exactly two $b$-tagged jets in the final state. A discriminant based on matrix element calculations is used to separate single-top-quark s-channel events from the main background contributions, which are top-quark pair production and $W$-boson production in association with jets. The observed (expected) signal significance over the background-only hypothesis is 3.3 (3.9) standard deviations, and the measured cross-section is $\sigma=8.2^{+3.5}_{-2.9}$ pb, consistent with the Standard Model prediction of $\sigma^{\mathrm{SM}}=10.32^{+0.40}_{-0.36}$ pb.

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.

Two related searches for phenomena beyond the standard model (BSM) are performed using events with hadronic jets and significant transverse momentum imbalance. The results are based on a sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2016-2018 and corresponding to an integrated luminosity of 137 fb$^{-1}$. The first search is inclusive, based on signal regions defined by the hadronic energy in the event, the jet multiplicity, the number of jets identified as originating from bottom quarks, and the value of the kinematic variable $M_\mathrm{T2}$ for events with at least two jets. For events with exactly one jet, the transverse momentum of the jet is used instead. The second search looks in addition for disappearing tracks produced by BSM long-lived charged particles that decay within the volume of the tracking detector. No excess event yield is observed above the predicted standard model background. This is used to constrain a range of BSM models that predict the following: the pair production of gluinos and squarks in the context of supersymmetry models conserving $R$-parity, with or without intermediate long-lived charginos produced in the decay chain; the resonant production of a colored scalar state decaying to a massive Dirac fermion and a quark; or the pair production of scalar and vector leptoquarks each decaying to a neutrino and a top, bottom, or light-flavor quark. In most of the cases, the results obtained are the most stringent constraints to date.

Cross-section measurements of top-quark pair production where the hadronically decaying top quark has transverse momentum greater than $355$ GeV and the other top quark decays into $\ell \nu b$ are presented using 139 fb$^{-1}$ of data collected by the ATLAS experiment during proton-proton collisions at the LHC. The fiducial cross-section at $\sqrt{s}=13$ TeV is measured to be $\sigma = 1.267 \pm 0.005 \pm 0.053$ pb, where the uncertainties reflect the limited number of data events and the systematic uncertainties, giving a total uncertainty of $4.2\%$. The cross-section is measured differentially as a function of variables characterising the $t\bar{t}$ system and additional radiation in the events. The results are compared with various Monte Carlo generators, including comparisons where the generators are reweighted to match a parton-level calculation at next-to-next-to-leading order. The reweighting improves the agreement between data and theory. The measured distribution of the top-quark transverse momentum is used to set limits on the Wilson coefficients of the dimension-six operators $O_{tG}$ and $O_{tq}^{(8)}$ in the effective field theory framework.

A search for leptoquarks decaying into the $b\tau$ final state is performed using Run 2 proton-proton collision data from the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb$^{-1}$ at $\sqrt{s} = 13$ TeV recorded by the ATLAS detector. The benchmark models considered in this search are vector leptoquarks with electric charge of 2/3e and scalar leptoquarks with an electric charge of 4/3e. No significant excess above the Standard Model prediction is observed, and 95% confidence level upper limits are set on the cross-section times branching fraction of leptoquarks decaying into $b\tau$. For the vector leptoquark production two models are considered: the Yang-Mills and Minimal coupling models. In the Yang-Mills (Minimal coupling) scenario, vector leptoquarks with a mass below 1.58 (1.35) TeV are excluded for a gauge coupling of 1.0 and below 2.05 (1.99) TeV for a gauge coupling of 2.5. In the case of scalar leptoquarks, masses below 1.28 TeV (1.53 TeV) are excluded for a Yukawa coupling of 1.0 (2.5). Finally, an interpretation of the results with minimal model dependence is performed for each of the signal region categories, and limits on the visible cross-section for beyond the Standard Model processes are provided.

Measurements of differential cross sections are presented for inclusive isolated-photon production in $pp$ collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 fb$^{-1}$ of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region.

This Letter reports the observation of single top quarks produced together with a photon, which directly probes the electroweak coupling of the top quark. The analysis uses 139 fb$^{-1}$ of 13 TeV proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. Requiring a photon with transverse momentum larger than 20 GeV and within the detector acceptance, the fiducial cross section is measured to be 688 $\pm$ 23 (stat.) $^{+75}_{-71}$ (syst.) fb, to be compared with the standard model prediction of 515 $^{+36}_{-42}$ fb at next-to-leading order in QCD.

Measurements of the suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

A search for the leptonic charge asymmetry ($A_\text{c}^{\ell}$) of top-quark$-$antiquark pair production in association with a $W$ boson ($t\bar{t}W$) is presented. The search is performed using final states with exactly three charged light leptons (electrons or muons) and is based on $\sqrt{s} = 13$ TeV proton$-$proton collision data collected with the ATLAS detector at the Large Hadron Collider at CERN during the years 2015$-$2018, corresponding to an integrated luminosity of 139 fb$^{-1}$. A profile-likelihood fit to the event yields in multiple regions corresponding to positive and negative differences between the pseudorapidities of the charged leptons from top-quark and top-antiquark decays is used to extract the charge asymmetry. At reconstruction level, the asymmetry is found to be $-0.123 \pm 0.136$ (stat.) $\pm \, 0.051$ (syst.). An unfolding procedure is applied to convert the result at reconstruction level into a charge-asymmetry value in a fiducial volume at particle level with the result of $-0.112 \pm 0.170$ (stat.) $\pm \, 0.054$ (syst.). The Standard Model expectations for these two observables are calculated using Monte Carlo simulations with next-to-leading-order plus parton shower precision in quantum chromodynamics and including next-to-leading-order electroweak corrections. They are $-0.084 \, ^{+0.005}_{-0.003}$ (scale) $\pm\, 0.006$ (MC stat.) and $-0.063 \, ^{+0.007}_{-0.004}$ (scale) $\pm\, 0.004$ (MC stat.) respectively, and in agreement with the measurements.

This paper presents a statistical combination of searches targeting final states with two top quarks and invisible particles, characterised by the presence of zero, one or two leptons, at least one jet originating from a $b$-quark and missing transverse momentum. The analyses are searches for phenomena beyond the Standard Model consistent with the direct production of dark matter in $pp$ collisions at the LHC, using 139 fb$^{-\text{1}}$ of data collected with the ATLAS detector at a centre-of-mass energy of 13 TeV. The results are interpreted in terms of simplified dark matter models with a spin-0 scalar or pseudoscalar mediator particle. In addition, the results are interpreted in terms of upper limits on the Higgs boson invisible branching ratio, where the Higgs boson is produced according to the Standard Model in association with a pair of top quarks. For scalar (pseudoscalar) dark matter models, with all couplings set to unity, the statistical combination extends the mass range excluded by the best of the individual channels by 50 (25) GeV, excluding mediator masses up to 370 GeV. In addition, the statistical combination improves the expected coupling exclusion reach by 14% (24%), assuming a scalar (pseudoscalar) mediator mass of 10 GeV. An upper limit on the Higgs boson invisible branching ratio of 0.38 (0.30$^{+\text{0.13}}_{-\text{0.09}}$) is observed (expected) at 95% confidence level.

A search for events with a dark photon produced in association with a dark Higgs boson via rare decays of the Standard Model $Z$ boson is presented, using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collision data recorded by the ATLAS detector at the Large Hadron Collider. The dark Higgs boson decays into a pair of dark photons, and at least two of the three dark photons must each decay into a pair of electrons or muons, resulting in at least two same-flavor opposite-charge lepton pairs in the final state. The data are found to be consistent with the background prediction, and upper limits are set on the dark photon's coupling to the dark Higgs boson times the kinetic mixing between the Standard Model photon and the dark photon, $\alpha_{D}\varepsilon^2$, in the dark photon mass range of $[5, 40]$ GeV except for the $\Upsilon$ mass window $[8.8, 11.1]$ GeV. This search explores new parameter space not previously excluded by other experiments.