A search for Supersymmetry involving the pair production of gluinos decaying via third-generation squarks to the lightest neutralino is reported. It uses an LHC proton--proton dataset at a center-of-mass energy $\sqrt{s} = 13$ TeV with an integrated luminosity of 3.2 fb$^{-1}$ collected with the ATLAS detector in 2015. The signal is searched for in events containing several energetic jets, of which at least three must be identified as $b$-jets, large missing transverse momentum and, potentially, isolated electrons or muons. Large-radius jets with a high mass are also used to identify highly boosted top quarks. No excess is found above the predicted background. For neutralino masses below approximately 700 GeV, gluino masses of less than 1.78 TeV and 1.76 TeV are excluded at the 95% CL in simplified models of the pair production of gluinos decaying via sbottom and stop, respectively. These results significantly extend the exclusion limits obtained with the $\sqrt{s} = 8$ TeV dataset.

This Letter presents the measurement of the fiducial and differential cross-sections of the electroweak production of a $Z \gamma$ pair in association with two jets. The analysis uses 140 fb$^{-1}$ of LHC proton-proton collision data taken at $\sqrt{s}$=13 TeV recorded by the ATLAS detector during the years 2015-2018. Events with a $Z$ boson candidate decaying into either an $e^+e^-$ or $\mu^+ \mu^-$ pair, a photon and two jets are selected. The electroweak component is extracted by requiring a large dijet invariant mass and a large rapidity gap between the two jets and is measured with an observed and expected significance well above five standard deviations. The fiducial $pp \rightarrow Z \gamma jj$ cross-section for the electroweak production is measured to be 3.6 $\pm$ 0.5 fb. The total fiducial cross-section that also includes contributions where the jets arise from strong interactions is measured to be $16.8^{+2.0}_{-1.8}$ fb. The results are consistent with the Standard Model predictions. Differential cross-sections are also measured using the same events and are compared with parton-shower Monte Carlo simulations. Good agreement is observed between data and predictions.

A search for non-resonant Higgs boson pair ($HH$) production is presented, in which one of the Higgs bosons decays to a b-quark pair ($b\bar b$) and the other decays to $WW^*$, $ZZ^*$, or $\tau^+\tau^-$, with in each case a final state with $\ell^+\ell^- +$ neutrinos ($\ell = e, \mu$). The analysis targets separately the gluon-gluon fusion and vector boson fusion production modes. Data recorded by the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider, corresponding to an integrated luminosity of $140\mathrm{fb}^{-1}$, are used in this analysis. Events are selected to have exactly two $b$-tagged jets and two leptons with opposite electric charge and missing transverse momentum in the final state. These events are classified using multivariate analysis algorithms to separate the $HH$ events from other Standard Model processes. No evidence of the signal is found. The observed (expected) upper limit on the cross-section for non-resonant Higgs boson pair production is determined to be 9.7 (16.2) times the Standard Model prediction at 95% confidence level. The Higgs boson self-interaction coupling parameter $\kappa_\lambda$ and the quadrilinear coupling parameter $\kappa_{2V}$ are each separately constrained by this analysis to be within the ranges ${[-6.2, 13.3]}$ and ${[-0.17, 2.4]}$, respectively, at 95% confidence level, when all other parameters are fixed.

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 search for supersymmetry involving the pair production of gluinos decaying via off-shell third-generation squarks into the lightest neutralino ($\tilde\chi^0_1$) is reported. It exploits LHC proton$-$proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV with an integrated luminosity of 139 fb$^{-1}$ collected with the ATLAS detector from 2015 to 2018. The search uses events containing large missing transverse momentum, up to one electron or muon, and several energetic jets, at least three of which must be identified as containing $b$-hadrons. Both a simple kinematic event selection and an event selection based upon a deep neural-network are used. No significant excess above the predicted background is found. In simplified models involving the pair production of gluinos that decay via off-shell top (bottom) squarks, gluino masses less than 2.44 TeV (2.35 TeV) are excluded at 95% CL for a massless $\tilde\chi^0_1$. Limits are also set on the gluino mass in models with variable branching ratios for gluino decays to $b\bar{b}\tilde\chi^0_1$, $t\bar{t}\tilde\chi^0_1$ and $t\bar{b}\tilde\chi^-_1$ / $\bar{t}b\tilde\chi^+_1$.

Searches for the exclusive decays of the Higgs boson into $D^*\gamma$ and of the $Z$ boson into $D^0\gamma$ and $K^0_s\gamma$ can probe flavour-violating Higgs and $Z$ boson couplings to light quarks. Searches for these decays are performed with a $pp$ collision data sample corresponding to an integrated luminosity of $136.3$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV between 2016-2018 with the ATLAS detector at the CERN Large Hadron Collider. In the $D^*\gamma$ and $D^0\gamma$ channels, the observed (expected) 95$\%$ confidence-level upper limits on the respective branching fractions are ${\cal B}(H\rightarrow D^*\gamma)< 1.0 (1.2)\times 10^{-3}$, ${\cal B}(Z\rightarrow D^0\gamma)< 4.0 (3.4)\times 10^{-6}$, while the corresponding results in the $K^0_s\gamma$ channel are ${\cal B}(Z\rightarrow K^0_s\gamma)< 3.1 (3.0)\times 10^{-6}$.

Several extensions of the Standard Model predict associated production of dark-matter particles with a Higgs boson. Such processes are searched for in final states with missing transverse momentum and a Higgs boson decaying to a $b\bar b$ pair with the ATLAS detector using 36.1 fb$^{-1}$ of $pp$ collisions at a center-of-mass energy of 13 TeV at the LHC. The observed data are in agreement with the Standard Model predictions and limits are placed on the associated production of dark-matter particles and a Higgs boson.

The results of a search for new heavy $W^\prime$ bosons decaying to an electron or muon and a neutrino using proton-proton collision data at a centre-of-mass energy of $\sqrt{s} = 13$ TeV are presented. The dataset was collected in 2015 and 2016 by the ATLAS experiment at the Large Hadron Collider and corresponds to an integrated luminosity of 36.1 fb$^{-1}$. As no excess of events above the Standard Model prediction is observed, the results are used to set upper limits on the $W^\prime$ boson cross-section times branching ratio to an electron or muon and a neutrino as a function of the $W^\prime$ mass. Assuming a $W^\prime$ boson with the same couplings as the Standard Model $W$ boson, $W^\prime$ masses below 5.1 TeV are excluded at the 95% confidence level.

A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

A search for leptoquark pair production decaying into $te^- \bar{t}e^+$ or $t\mu^- \bar{t}\mu^+$ in final states with multiple leptons is presented. The search is based on a dataset of $pp$ collisions at $\sqrt{s}$=13 TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 139~fb$^{-1}$. Events are selected with two or more light leptons (electron or muon) and at least two jets out of which at least one jet is identified as coming from a $b$-hadron. Four signal regions, with the requirement of at least three light leptons, are considered based on the number of leptons of a given flavour. The main background processes are estimated using dedicated control regions in a simultaneous fit with the signal regions to data. No excess above the Standard Model background prediction is observed and 95% confidence level limits on the production cross section times branching ratio are derived as a function of the leptoquark mass. Under the assumption of exclusive decays into $te^{-}$ ($t\mu^{-}$), the corresponding lower limit on the scalar mixed-generation leptoquark mass $m_{\mathrm{LQ}_{\mathrm{mix}}^{\mathrm{d}}}$ is at 1.58 (1.59) TeV and on the vector leptoquark mass $m_{\tilde{U}_1}$ at 1.67 (1.67) TeV in the minimal coupling scenario and at 1.95 (1.95) TeV in the Yang-Mills scenario.