A search for heavy resonances decaying into a pair of $Z$ bosons leading to $\ell^+\ell^-\ell'^+\ell'^-$ and $\ell^+\ell^-\nu\bar\nu$ final states, where $\ell$ stands for either an electron or a muon, is presented. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV collected from 2015 to 2018 that corresponds to the full integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Different mass ranges spanning 200 GeV to 2000 GeV for the hypothetical resonances are considered, depending on the final state and model. In the absence of a significant observed excess, the results are interpreted as upper limits on the production cross section of a spin-0 or spin-2 resonance. The upper limits for the spin-0 resonance are translated to exclusion contours in the context of Type-I and Type-II two-Higgs-doublet models, and the limits for the spin-2 resonance are used to constrain the Randall--Sundrum model with an extra dimension giving rise to spin-2 graviton excitations.

A search is presented for the pair production of heavy vector-like $T$ and $B$ quarks, primarily targeting the $T$ quark decays to a $W$ boson and a $b$-quark. The search is based on $36.1$ fb$^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Data are analysed in the lepton-plus-jets final state, including at least one $b$-tagged jet and a large-radius jet identified as originating from the hadronic decay of a high-momentum $W$ boson. No significant deviation from the Standard Model expectation is observed in the reconstructed $T$ mass distribution. The observed 95% confidence level lower limit on the $T$ mass are 1350 GeV assuming 100% branching ratio to $Wb$. In the SU(2) singlet scenario, the lower mass limit is 1170 GeV. This search is also sensitive to a heavy vector-like $B$ quark decaying to $Wt$ and other final states. The results are thus reinterpreted to provide a 95% confidence level lower limit on the $B$ quark mass at 1250 GeV assuming 100% branching ratio to $Wt$; in the SU(2) singlet scenario, the limit is 1080 GeV. Mass limits on both $T$ and $B$ production are also set as a function of the decay branching ratios. The 100% branching ratio limits are found to be applicable to heavy vector-like $Y$ and $X$ production that decay to $Wb$ and $Wt$, respectively.

A search for a heavy neutral Higgs boson, $A$, decaying into a $Z$ boson and another heavy Higgs boson, $H$, is performed using a data sample corresponding to an integrated luminosity of 139 fb$^{-1}$ from proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector at the LHC. The search considers the $Z$ boson decaying into electrons or muons and the $H$ boson into a pair of $b$-quarks or $W$ bosons. The mass range considered is 230-800 GeV for the $A$ boson and 130-700 GeV for the $H$ boson. The data are in good agreement with the background predicted by the Standard Model, and therefore 95% confidence-level upper limits for $\sigma \times B(A\rightarrow ZH) \times B(H\rightarrow bb$ or $H\rightarrow WW)$ are set. The upper limits are in the range 0.0062-0.380 pb for the $H\rightarrow bb$ channel and in the range 0.023-8.9 pb for the $H\rightarrow WW$ channel. An interpretation of the results in the context of two-Higgs-Doublet models is also given.

A search for flavour-changing neutral-current decays of a top quark into an up-type quark (either up or charm) and a light scalar particle $X$ decaying into a bottom anti-bottom quark pair is presented. The search focuses on top-quark pair production where one top quark decays to $qX$, with $X \rightarrow b\bar{b}$, and the other top quark decays according to the Standard Model, with the $W$ boson decaying leptonically. The final state is thus characterised by an isolated electron or muon and at least four jets. Events are categorised according to the multiplicity of jets and jets tagged as originating from $b$-quarks, and a neural network is used to discriminate between signal and background processes. The data analysed correspond to 139 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV, recorded with the ATLAS detector at the LHC. The 95% confidence-level upper limits between 0.019% and 0.062% are derived for the branching fraction $\mathcal{B}$($t \rightarrow uX$) and between 0.018% and 0.078% for the branching fraction $\mathcal{B}$($t \rightarrow cX$), for masses of the scalar particle $X$ between 20 and 160 GeV.

This article presents a search for new resonances decaying into a $Z$ or $W$ boson and a 125 GeV Higgs boson $h$, and it targets the $\nu\bar{\nu}b\bar{b}$, $\ell^+\ell^-b\bar{b}$, or $\ell^{\pm}{\nu}b\bar{b}$ final states, where $\ell=e$ or $\mu$, in proton-proton collisions at $\sqrt{s}=13$ TeV. The data used correspond to a total integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector during Run 2 of the LHC at CERN. The search is conducted by examining the reconstructed invariant or transverse mass distributions of $Zh$ or $Wh$ candidates for evidence of a localised excess in the mass range from 220 GeV to 5 TeV. No significant excess is observed and 95% confidence-level upper limits between 1.3 pb and 0.3 fb are placed on the production cross section times branching fraction of neutral and charged spin-1 resonances and CP-odd scalar bosons. These limits are converted into constraints on the parameter space of the Heavy Vector Triplet model and the two-Higgs-doublet model.

A search for a new $Z'$ gauge boson predicted by $L_{\mu}-L_{\tau}$ models, based on charged-current Drell-Yan production, $pp \rightarrow W^{\pm(*)} \rightarrow Z' \mu^{\pm} \nu \rightarrow \mu^{\pm}\mu^{\mp}\mu^{\pm}\nu$, is presented. The data sample used corresponds to an integrated luminosity of 140 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} = 13$ TeV recorded by the ATLAS detector at the Large Hadron Collider. The search examines a final state of $3\mu$ plus large missing transverse momentum. Upper limits are set on the $Z'$ production cross-section times branching ratio in the mass range of 5-81 GeV. After combining with the previous $Z'$ search using the neutral-current Drell-Yan production with a $4\mu$ final state, the most stringent exclusion limits to date are achieved in the parameter space of the $Z'$ coupling strength and mass.

A measurement of the $ZZ$ production in the $\ell^{-}\ell^{+}\ell^{\prime -}\ell^{\prime +}$ and $\ell^{-}\ell^{+}\nu\bar{\nu}$ channels $(\ell = e, \mu)$ in proton--proton collisions at $\sqrt{s} = 8$ TeV at the Large Hadron Collider at CERN, using data corresponding to an integrated luminosity of 20.3 fb$^{-1}$ collected by the ATLAS experiment in 2012 is presented. The fiducial cross sections for $ZZ\to\ell^{-}\ell^{+}\ell^{\prime -}\ell^{\prime +}$ and $ZZ\to \ell^{-}\ell^{+}\nu\bar{\nu}$ are measured in selected phase-space regions. The total cross section for $ZZ$ events produced with both $Z$ bosons in the mass range 66 to 116 GeV is measured from the combination of the two channels to be $7.3\pm0.4\textrm{(stat)}\pm0.3\textrm{(syst)}\pm0.2\textrm{(lumi)}$ pb, which is consistent with the Standard Model prediction of $6.6^{+0.7}_{-0.6}$ pb. The differential cross sections in bins of various kinematic variables are presented. The differential event yield as a function of the transverse momentum of the leading $Z$ boson is used to set limits on anomalous neutral triple gauge boson couplings in $ZZ$ production.

This paper presents a study of the production of $WW$ or $WZ$ boson pairs, with one $W$ boson decaying to $e\nu$ or $\mu\nu$ and one $W$ or $Z$ boson decaying hadronically. The analysis uses 20.2 fb$^{-1}$ of $\sqrt{s}=8$ TeV $pp$ collision data, collected by the ATLAS detector at the Large Hadron Collider. Cross-sections for $WW/WZ$ production are measured in high-$p_{T}$ fiducial regions defined close to the experimental event selection. The cross-section is measured for the case where the hadronically decaying boson is reconstructed as two resolved jets, and the case where it is reconstructed as a single jet. The transverse momentum distribution of the hadronically decaying boson is used to search for new physics. Observations are consistent with the Standard Model predictions, and $95\%$ confidence intervals are calculated for parameters describing anomalous triple gauge-boson couplings.

A search for heavy resonances decaying to a $W$ or $Z$ boson and a Higgs boson in the $q\bar{q}^{(\prime)}b\bar{b}$ final state is described. The search uses 36.1 fb$^{-1}$ of proton-proton collision data at $\sqrt{s} =$ 13 TeV collected by the ATLAS detector at the CERN Large Hadron Collider in 2015 and 2016. The data are in agreement with the Standard Model expectations, with the largest excess found at a resonance mass of 3.0 TeV with a local (global) significance of 3.3 (2.1) $\sigma$. The results are presented in terms of constraints on a simplified model with a heavy vector triplet. Upper limits are set on the production cross-section times branching ratio for resonances decaying to a $W$ ($Z$) boson and a Higgs boson, itself decaying to $b\bar{b}$, in the mass range between 1.1 and 3.8 TeV; the limits range between 83 and 1.6 fb (77 and 1.1 fb) at 95% confidence level.

A search is conducted for new resonances decaying into a $W$ or $Z$ boson and a 125 GeV Higgs boson in the $\nu\bar{\nu}b\bar{b}$, $\ell^{\pm}{\nu}b\bar{b}$, and $\ell^+\ell^-b\bar{b}$ final states, where $\ell ^{\pm}= e^{\pm}$ or $\mu^{\pm}$, in $pp$ collisions at $\sqrt s = 13$ TeV. The data used correspond to a total integrated luminosity of 36.1 fb$^{-1}$ collected with the ATLAS detector at the Large Hadron Collider during the 2015 and 2016 data-taking periods. The search is conducted by examining the reconstructed invariant or transverse mass distributions of $Wh$ and $Zh$ candidates for evidence of a localised excess in the mass range of 220 GeV up to 5 TeV. No significant excess is observed and the results are interpreted in terms of constraints on the production cross-section times branching fraction of heavy $W^\prime$ and $Z^\prime$ resonances in heavy-vector-triplet models and the CP-odd scalar boson $A$ in two-Higgs-doublet models. Upper limits are placed at the 95 % confidence level and range between $9.0\times 10^{-4}$ pb and $8.1\times 10^{-1}$ pb depending on the model and mass of the resonance.

A search for a charged Higgs boson, $H^{\pm}$, produced in top-quark decays, $t \rightarrow H^{\pm}b$, is presented. The search targets $H^{\pm}$ decays into a bottom and a charm quark, $H^{\pm} \rightarrow cb$. The analysis focuses on a selection enriched in top-quark pair production, where one top quark decays into a leptonically decaying $W$ boson and a bottom quark, and the other top quark decays into a charged Higgs boson and a bottom quark. This topology leads to a lepton-plus-jets final state, characterised by an isolated electron or muon and at least four jets. The search exploits the high multiplicity of jets containing $b$-hadrons, and deploys a neural network classifier that uses the kinematic differences between the signal and the background. The search uses a dataset of proton-proton collisions collected at a centre-of-mass energy $\sqrt{s}=13$ TeV between 2015 and 2018 with the ATLAS detector at CERN's Large Hadron Collider, amounting to an integrated luminosity of 139 fb$^{-1}$. Observed (expected) 95% confidence-level upper limits between 0.15% (0.09%) and 0.42% (0.25%) are derived for the product of branching fractions $\mathscr{B}(t\rightarrow H^{\pm}b) \times \mathscr{B}(H^{\pm}\rightarrow cb)$ for charged Higgs boson masses between 60 and 160 GeV, assuming the SM production of the top-quark pairs.

A search is presented for the pair production of heavy vector-like $B$ quarks, primarily targeting $B$ quark decays into a $W$ boson and a top quark. The search is based on $36.1$ $fb^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Data are analysed in the lepton-plus-jets final state, characterised by a high-transverse-momentum isolated electron or muon, large missing transverse momentum, and multiple jets, of which at least one is $b$-tagged. No significant deviation from the Standard Model expectation is observed. The 95% confidence level lower limit on the $B$ mass is 1350 GeV assuming a 100% branching ratio to $Wt$. In the SU(2) singlet scenario, the lower mass limit is 1170 GeV. This search is also sensitive to a heavy vector-like $B$ quark decaying into other final states ($Zb$ and $Hb$) and thus mass limits on $B$ production are set as a function of the decay branching ratios. The 100% branching ratio limits are found to be also applicable to heavy vector-like $X$ production, with charge $+$5/3, that decay into $Wt$.

Several extensions of the Standard Model predict the production of dark matter particles at the LHC. An uncharted signature of dark matter particles produced in association with $VV=W^\pm W^\mp$ or $ZZ$ pairs from a decay of a dark Higgs boson $s$ is searched for using 139 fb$^{-1}$ of $pp$ collisions recorded by the ATLAS detector at a center-of-mass energy of 13 TeV. The $s\to V(q\bar q)V(q\bar q)$ decays are reconstructed with a novel technique aimed at resolving the dense topology from boosted $VV$ pairs using jets in the calorimeter and tracking information. Dark Higgs scenarios with $m_s > 160$ GeV are excluded.

Results of a search for the electroweak associated production of charginos and next-to-lightest neutralinos, pairs of charginos or pairs of tau sleptons are presented. These processes are characterised by final states with at least two hadronically decaying tau leptons, missing transverse momentum and low jet activity. The analysis is based on an integrated luminosity of 20.3 fb$^{-1}$ of proton--proton collisions at $\sqrt{s}=8$ TeV recorded with the ATLAS experiment at the Large Hadron Collider. No significant excess is observed with respect to the predictions from Standard Model processes. Limits are set at 95% confidence level on the masses of the lighter chargino and next-to-lightest neutralino for various hypotheses for the lightest neutralino mass in simplified models. In the scenario of direct production of chargino pairs, with each chargino decaying into the lightest neutralino via an intermediate tau slepton, chargino masses up to 345 GeV are excluded for a massless lightest neutralino. For associated production of mass-degenerate charginos and next-to-lightest neutralinos, both decaying into the lightest neutralino via an intermediate tau slepton, masses up to 410 GeV are excluded for a massless lightest neutralino.

The results of a search for electroweakino pair production $pp \rightarrow \tilde\chi^\pm_1 \tilde\chi^0_2$ in which the chargino ($\tilde\chi^\pm_1$) decays into a $W$ boson and the lightest neutralino ($\tilde\chi^0_1$), while the heavier neutralino ($\tilde\chi^0_2$) decays into the Standard Model 125 GeV Higgs boson and a second $\tilde\chi^0_1$ are presented. The signal selection requires a pair of $b$-tagged jets consistent with those from a Higgs boson decay, and either an electron or a muon from the $W$ boson decay, together with missing transverse momentum from the corresponding neutrino and the stable neutralinos. The analysis is based on data corresponding to 139 $\mathrm{fb}^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collisions provided by the Large Hadron Collider and recorded by the ATLAS detector. No statistically significant evidence of an excess of events above the Standard Model expectation is found. Limits are set on the direct production of the electroweakinos in simplified models, assuming pure wino cross-sections. Masses of $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ up to 740 GeV are excluded at 95% confidence level for a massless $\tilde{\chi}^{0}_{1}$.

A search for a heavy neutral Higgs boson, $A$, decaying into a $Z$ boson and another heavy Higgs boson, $H$, is performed using a data sample corresponding to an integrated luminosity of 36.1 fb$^{-1}$ from proton-proton collisions at $\sqrt{s} = 13$ TeV recorded in 2015 and 2016 by the ATLAS detector at the Large Hadron Collider. The search considers the $Z$ boson decaying to electrons or muons and the $H$ boson into a pair of $b$-quarks. No evidence for the production of an $A$ boson is found. Considering each production process separately, the 95% confidence-level upper limits on the $pp\rightarrow A\rightarrow ZH$ production cross-section times the branching ratio $H\rightarrow bb$ are in the range of 14-830 fb for the gluon-gluon fusion process and 26-570 fb for the $b$-associated process for the mass ranges 130-700 GeV of the $H$ boson and process for the mass ranges 130-700 GeV of the $H$ boson and 230-800 GeV of the $A$ boson. The results are interpreted in the context of the two-Higgs-doublet model.

Measurements of jet activity in top-quark pair events produced in proton--proton collisions are presented, using 3.2 fb$^{-1}$ of $pp$ collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Events are chosen by requiring an opposite-charge $e\mu$ pair and two $b$-tagged jets in the final state. The normalised differential cross-sections of top-quark pair production are presented as functions of additional-jet multiplicity and transverse momentum, $p_{\mathrm T}$. The fraction of signal events that do not contain additional jet activity in a given rapidity region, the gap fraction, is measured as a function of the $p_{\mathrm T}$ threshold for additional jets, and is also presented for different invariant mass regions of the $e\mu b\bar{b}$ system. All measurements are corrected for detector effects and presented as particle-level distributions compared to predictions with different theoretical approaches for QCD radiation. While the kinematics of the jets from top-quark decays are described well, the generators show differing levels of agreement with the measurements of observables that depend on the production of additional jets.

A search for flavour-changing neutral current (FCNC) $tqH$ interactions involving a top quark, another up-type quark ($q=u$, $c$), and a Standard Model (SM) Higgs boson decaying into a $\tau$-lepton pair ($H\rightarrow \tau^+\tau^-$) is presented. The search is based on a dataset of $pp$ collisions at $\sqrt{s}=13$ TeV that corresponds to an integrated luminosity of 139 fb$^{-1}$ recorded with the ATLAS detector at the Large Hadron Collider. Two processes are considered: single top quark FCNC production in association with a Higgs boson ($pp\rightarrow tH$), and top quark pair production in which one of the top quarks decays into $Wb$ and the other decays into $qH$ through the FCNC interactions. The search selects events with two hadronically decaying $\tau$-lepton candidates ($\tau_{\text{had}}$) or at least one $\tau_{\text{had}}$ with an additional lepton ($e$, $\mu$), as well as multiple jets. Event kinematics is used to separate signal from the background through a multivariate discriminant. A slight excess of data is observed with a significance of 2.3$\sigma$ above the expected SM background, and 95% CL upper limits on the $t\to qH$ branching ratios are derived. The observed (expected) 95% CL upper limits set on the $t\to cH$ and $t\to uH$ branching ratios are $9.4 \times 10^{-4}$ $(4.8^{+2.2}_{-1.4}\times 10^{-4})$ and $6.9\times 10^{-4}$ $(3.5^{+1.5}_{-1.0}\times 10^{-4})$, respectively. The corresponding combined observed (expected) upper limits on the dimension-6 operator Wilson coefficients in the effective $tqH$ couplings are $C_{c\phi} <1.35$ $(0.97)$ and $C_{u\phi} <1.16$ $(0.82)$.

The production rates of prompt and non-prompt $J/\psi$ and $\psi(2\mathrm{S})$ mesons are measured using 2.1 $fb^{-1}$ and 11.4 $fb^{-1}$ of data collected with the ATLAS experiment at the LHC, in proton-proton collisions at $\sqrt{s}=7$ and 8 TeV respectively. Production cross-sections for both prompt and non-prompt production sources, ratios of $\psi(2\mathrm{S})$ to $J/\psi$ production, and fractions of non-prompt to inclusive production for $J/\psi$ and $\psi(2\mathrm{S})$ are measured double-differentially as a function of meson $p_{T}$ and rapidity. These measurements are made in a restricted fiducial volume and also corrected for geometrical acceptance after which they are compared to a variety of theoretical predictions.

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 is presented for direct top squark pair production in final states with one isolated electron or muon, jets, and missing transverse momentum in proton-proton collisions at sqrt(s) = 7 TeV. The measurement is based on 4.7 fb-1 of data collected with the ATLAS detector at the LHC. Each top squark is assumed to decay to a top quark and the lightest supersymmetric particle (LSP). The data are found to be consistent with Standard Model expectations. Top squark masses between 230 GeV and 440 GeV are excluded with 95% confidence for massless LSPs, and top squark masses around 400 GeV are excluded for LSP masses up to 125 GeV.

Measurements of fiducial cross sections for the electroweak production of two jets in association with a Z-boson are presented. The measurements are performed using 20.3 inverse femtobarns of proton-proton collision data collected at a centre-of-mass energy of sqrt(s)=8 TeV by the ATLAS experiment at the Large Hadron Collider. The electroweak component is extracted by a fit to the dijet invariant mass distribution in a fiducial region chosen to enhance the electroweak contribution over the dominant background in which the jets are produced via the strong interaction. The electroweak cross sections measured in two fiducial regions are in good agreement with the Standard Model expectations and the background-only hypothesis is rejected with significance above the 5 sigma level. The electroweak process includes the vector boson fusion production of a Z-boson and the data are used to place limits on anomalous triple gauge boson couplings. In addition, measurements of cross sections and differential distributions for inclusive Z-boson-plus-dijet production are performed in five fiducial regions, each with different sensitivity to the electroweak contribution. The results are corrected for detector effects and compared to predictions from the SHERPA and POWHEG event generators.

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.

High transverse momentum jets produced in pp collisions at a centre of mass energy of 7 TeV are used to measure the transverse energy-energy correlation function and its associated azimuthal asymmetry. The data were recorded with the ATLAS detector at the LHC in the year 2011 and correspond to an integrated luminosity of 158 $\mathrm{pb}^{-1}$. The selection criteria demand the average transverse momentum of the two leading jets in an event to be larger than 250 GeV. The data at detector level are well described by Monte Carlo event generators. They are unfolded to the particle level and compared with theoretical calculations at next-to-leading-order accuracy. The agreement between data and theory is good and provides a precision test of perturbative Quantum Chromodynamics at large momentum transfers. From this comparison, the strong coupling constant given at the $Z$ boson mass is determined to be $\alpha_{\mathrm{s}}(m_Z) = 0.1173 \pm 0.0010 \mbox{ (exp.) }^{+0.0065}_{-0.0026} \mbox{ (theo.)}$.

A search for new phenomena in final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair, jets, and large missing transverse momentum is presented. This analysis makes use of proton--proton collision data with an integrated luminosity of $36.1 \; \mathrm{fb}^{-1}$, collected during 2015 and 2016 at a centre-of-mass energy $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The search targets the pair production of supersymmetric coloured particles (squarks or gluinos) and their decays into final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair and the lightest neutralino ($\tilde{\chi}_1^0$) via one of two next-to-lightest neutralino ($\tilde{\chi}_2^0$) decay mechanisms: $\tilde{\chi}_2^0 \rightarrow Z \tilde{\chi}_1^0$, where the $Z$ boson decays leptonically leading to a peak in the dilepton invariant mass distribution around the $Z$ boson mass; and $\tilde{\chi}_2^0 \rightarrow \ell^+\ell^- \tilde{\chi}_1^0$ with no intermediate $\ell^+\ell^-$ resonance, yielding a kinematic endpoint in the dilepton invariant mass spectrum. The data are found to be consistent with the Standard Model expectation. Results are interpreted using simplified models, and exclude gluinos and squarks with masses as large as 1.85 TeV and 1.3 TeV at 95% confidence level, respectively.