The production of $W^{\pm}Z$ events in proton--proton collisions at a centre-of-mass energy of 13 TeV is measured with the ATLAS detector at the LHC. The collected data correspond to an integrated luminosity of 3.2 fb$^{-1}$. The $W^{\pm}Z$ candidates are reconstructed using leptonic decays of the gauge bosons into electrons or muons. The measured inclusive cross section in the detector fiducial region for leptonic decay modes is $\sigma_{W^\pm Z \rightarrow \ell^{'} \nu \ell \ell}^{\textrm{fid.}} = 63.2 \pm 3.2$ (stat.) $\pm 2.6$ (sys.) $\pm 1.5$ (lumi.) fb. In comparison, the next-to-leading-order Standard Model prediction is $53.4^{+3.6}_{-2.8}$ fb. The extrapolation of the measurement from the fiducial to the total phase space yields $\sigma_{W^{\pm}Z}^{\textrm{tot.}} = 50.6 \pm 2.6$ (stat.) $\pm 2.0$ (sys.) $\pm 0.9$ (th.) $\pm 1.2$ (lumi.) pb, in agreement with a recent next-to-next-to-leading-order calculation of $48.2^{+1.1}_{-1.0}$ pb. The cross section as a function of jet multiplicity is also measured, together with the charge-dependent $W^+Z$ and $W^-Z$ cross sections and their ratio.

The results of a search for the stop, the supersymmetric partner of the top quark, in final states with one isolated electron or muon, jets, and missing transverse momentum are reported. The search uses the 2015 LHC $pp$ collision data at a center-of-mass energy of $\sqrt{s}=13$ TeV recorded by the ATLAS detector and corresponding to an integrated luminosity of 3.2 fb${}^{-1}$. The analysis targets two types of signal models: gluino-mediated pair production of stops with a nearly mass-degenerate stop and neutralino; and direct pair production of stops, decaying to the top quark and the lightest neutralino. The experimental signature in both signal scenarios is similar to that of a top quark pair produced in association with large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits on gluino and stop masses are set at 95% confidence level. The results extend the LHC Run-1 exclusion limit on the gluino mass up to 1460 GeV in the gluino-mediated scenario in the high gluino and low stop mass region, and add an excluded stop mass region from 745 to 780 GeV for the direct stop model with a massless lightest neutralino. The results are also reinterpreted to set exclusion limits in a model of vector-like top quarks.

A measurement of the total $pp$ cross section at the LHC at $\sqrt{s}=8$ TeV is presented. An integrated luminosity of $500$ $\mu$b$^{-1}$ was accumulated in a special run with high-$\beta^{\star}$ beam optics to measure the differential elastic cross section as a function of the Mandelstam momentum transfer variable $t$. The measurement is performed with the ALFA sub-detector of ATLAS. Using a fit to the differential elastic cross section in the $-t$ range from $0.014$ GeV$^2$ to $0.1$ GeV$^2$ to extrapolate $t\rightarrow 0$, the total cross section, $\sigma_{\mathrm{tot}}(pp\rightarrow X)$, is measured via the optical theorem to be: $\sigma_{\mathrm{tot}}(pp\rightarrow X) = {96.07} \; \pm 0.18 \; ({{stat.}}) \pm 0.85 \; ({{exp.}}) \pm 0.31 \; ({extr.}) \; {mb} \;,$ where the first error is statistical, the second accounts for all experimental systematic uncertainties and the last is related to uncertainties in the extrapolation $t\rightarrow 0$. In addition, the slope of the exponential function describing the elastic cross section at small $t$ is determined to be $B = 19.74 \pm 0.05 \; ({{stat.}}) \pm 0.23 \; ({{syst.}}) \; {GeV}^{-2}$.

The result of a search for pair production of the supersymmetric partner of the Standard Model bottom quark ($\tilde{b}_1$) is reported. The search uses 3.2 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=$13 TeV collected by the ATLAS experiment at the Large Hadron Collider in 2015. Bottom squarks are searched for in events containing large missing transverse momentum and exactly two jets identified as originating from $b$-quarks. No excess above the expected Standard Model background yield is observed. Exclusion limits at 95% confidence level on the mass of the bottom squark are derived in phenomenological supersymmetric $R$-parity-conserving models in which the $\tilde{b}_1$ is the lightest squark and is assumed to decay exclusively via $\tilde{b}_1 \rightarrow b \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ is the lightest neutralino. The limits significantly extend previous results; bottom squark masses up to 800 (840) GeV are excluded for the $\tilde{\chi}_1^0$ mass below 360 (100) GeV whilst differences in mass above 100 GeV between the $\tilde{b}_1$ and the $\tilde{\chi}_1^0$ are excluded up to a $\tilde{b}_1$ mass of 500 GeV.

A search for squarks and gluinos in final states containing hadronic jets, missing transverse momentum but no electrons or muons is presented. The data were recorded in 2015 by the ATLAS experiment in $\sqrt{s}=$ 13 TeV proton--proton collisions at the Large Hadron Collider. No excess above the Standard Model background expectation was observed in 3.2 fb$^{-1}$ of analyzed data. Results are interpreted within simplified models that assume R-parity is conserved and the neutralino is the lightest supersymmetric particle. An exclusion limit at the 95% confidence level on the mass of the gluino is set at 1.51 TeV for a simplified model incorporating only a gluino octet and the lightest neutralino, assuming the lightest neutralino is massless. For a simplified model involving the strong production of mass-degenerate first- and second-generation squarks, squark masses below 1.03 TeV are excluded for a massless lightest neutralino. These limits substantially extend the region of supersymmetric parameter space excluded by previous measurements with the ATLAS detector.

Measurements of the electroweak production of a $W$ boson in association with two jets at high dijet invariant mass are performed using $\sqrt{s} = 7$ and $8$ TeV proton-proton collision data produced by the Large Hadron Collider, corresponding respectively to 4.7 and 20.2 fb$^{-1}$ of integrated luminosity collected by the ATLAS detector. The measurements are sensitive to the production of a $W$ boson via a triple-gauge-boson vertex and include both the fiducial and differential cross sections of the electroweak process.

Measurements of the production cross section of a $Z$ boson in association with jets in proton-proton collisions at $\sqrt{s} = 13$ TeV are presented, using data corresponding to an integrated luminosity of 3.16 fb$^{-1}$ collected by the ATLAS experiment at the CERN Large Hadron Collider in 2015. Inclusive and differential cross sections are measured for events containing a $Z$ boson decaying to electrons or muons and produced in association with up to seven jets with $p_T > 30$ GeV and $|y| <2.5$. Predictions from different Monte Carlo generators based on leading-order and next-to-leading-order matrix elements for up to two additional partons interfaced with parton shower and fixed-order predictions at next-to-leading order and next-to-next-to-leading order are compared with the measured cross sections. Good agreement within the uncertainties is observed for most of the modelled quantities, in particular with the generators which use next-to-leading-order matrix elements and the more recent next-to-next-to-leading-order fixed-order predictions.

Two searches for new phenomena in final states containing a same-flavour opposite-lepton (electron or muon) pair, jets, and large missing transverse momentum are presented. These searches make use of proton--proton collision data, collected during 2015 and 2016 at a centre-of-mass energy $\sqrt{s}=13$ TeV by the ATLAS detector at the Large Hadron Collider, which correspond to an integrated luminosity of 14.7 fb$^{-1}$. Both searches target the pair production of supersymmetric particles, squarks or gluinos, which decay to final states containing a same-flavour opposite-sign lepton pair via one of two mechanisms: a leptonically decaying Z boson in the final state, leading to a peak in the dilepton invariant-mass distribution around the Z boson mass; and decays of neutralinos (e.g. $\tilde{\chi}_2^0 \rightarrow \ell^+\ell^- \tilde{\chi}_1^0$), 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 in simplified models of gluino-pair (squark-pair) production, and provide sensitivity to gluinos (squarks) with masses as large as 1.70 TeV (980 GeV).

We present charged-particle distributions sensitive to the underlying event, measured by the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV, in low-luminosity Large Hadron Collider fills corresponding to an integrated luminosity of 1.6 nb$^{-1}$. The distributions were constructed using charged particles with absolute pseudorapidity less than 2.5 and with transverse momentum greater than 500 MeV, in events with at least one such charged particle with transverse momentum above 1 GeV. These distributions characterise the angular distribution of energy and particle flows with respect to the charged particle with highest transverse momentum, as a function of both that momentum and of charged-particle multiplicity. The results have been corrected for detector effects and are compared to the predictions of various Monte Carlo event generators, experimentally establishing the level of underlying-event activity at LHC Run 2 energies and providing inputs for the development of event generator modelling. The current models in use for UE modelling typically describe this data to 5% accuracy, compared with data uncertainties of less than 1%.

A search for supersymmetry in events with large missing transverse momentum, jets, and at least one hadronically decaying tau lepton has been performed using 3.2 fb$^{-1}$ of proton-proton collision data at $\sqrt{s}=13$ TeV recorded by the ATLAS detector at the Large Hadron Collider in 2015. Two exclusive final states are considered, with either exactly one or at least two tau leptons. No excess over the Standard Model prediction is observed in the data. Results are interpreted in the context of gauge-mediated supersymmetry breaking and a simplified model of gluino pair production with tau-rich cascade decays, substantially improving on previous limits. In the GMSB model considered, supersymmetry-breaking scale ($\Lambda$) values below 92 TeV are excluded at the 95% confidence level, corresponding to gluino masses below 2000 GeV. For large values of $\tan\beta$, values of $\Lambda$ up to 107 TeV and gluino masses up to 2300 GeV are excluded. In the simplified model, gluino masses are excluded up to 1570 GeV for neutralino masses around 100 GeV. Neutralino masses up to 700 GeV are excluded for all gluino masses between 800 GeV and 1500 GeV, while the strongest exclusion of 750 GeV is achieved for gluino masses around 1400 GeV.

High-precision measurements by the ATLAS Collaboration are presented of inclusive $W^+\to\ell^+\nu$, $W^-\to\ell^-\bar{\nu}$ and $Z/\gamma^*\to\ell\ell$ ($\ell=e,\mu$) Drell-Yan production cross sections at the LHC. The data were collected in proton-proton collisions at $\sqrt{s} = 7$ TeV with an integrated luminosity of 4.6 fb$^{-1}$. Differential $W^+$ and $W^-$ cross sections are measured in a lepton pseudorapidity range $|\eta_{\ell}| = 2.5$. Differential $Z/\gamma^*$ cross sections are measured as a function of the absolute dilepton rapidity, for $|y_{\ell\ell}| < 3.6$, for three intervals of dilepton mass, $m_{\ell\ell}$, extending from 46 to 150 GeV. The integrated and differential electron- and muon-channel cross sections are combined and compared to theoretical predictions using recent sets of parton distribution functions. The data, together with the final inclusive $e^{\pm}p$ scattering cross-section data from H1 and ZEUS, are interpreted in a next-to-next-to-leading-order QCD analysis, and a new set of parton distribution functions, ATLAS-epWZ16, is obtained. The ratio of strange-to-light sea-quark densities in the proton is determined more accurately than in previous determinations based on collider data only, and is established to be close to unity in the sensitivity range of the data. A new measurement of the CKM matrix element $|V_{cs}|$ is also provided.

Measurements of normalized differential cross-sections of top quark pair ($t\bar t$) production are presented as a function of the mass, the transverse momentum and the rapidity of the $t\bar t$ system in proton-proton collisions at center-of-mass energies of $\sqrt{s}$ = 7 TeV and 8 TeV. The dataset corresponds to an integrated luminosity of 4.6 fb$^{-1}$ at 7 TeV and 20.2 fb$^{-1}$ at 8 TeV, recorded with the ATLAS detector at the Large Hadron Collider. Events with top quark pair signatures are selected in the dilepton final state, requiring exactly two charged leptons and at least two jets with at least one of the jets identified as likely to contain a $b$-hadron. The measured distributions are corrected for detector effects and selection efficiency to cross-sections at the parton level. The differential cross-sections are compared with different Monte Carlo generators and theoretical calculations of $t\bar t$ production. The results are consistent with the majority of predictions in a wide kinematic range.

A search for doubly charged Higgs bosons with pairs of prompt, isolated, highly energetic leptons with the same electric charge is presented. The search uses a proton-proton collision data sample at a centre-of-mass energy of 13 TeV corresponding to 36.1 $\mathrm{fb}^{-1}$ of integrated luminosity recorded in 2015 and 2016 by the ATLAS detector at the LHC. This analysis focuses on the decays $H^{\pm\pm}\rightarrow e^{\pm}e^{\pm}$, $H^{\pm\pm}\rightarrow e^{\pm}\mu^{\pm}$ and $H^{\pm\pm}\rightarrow \mu^{\pm}\mu^{\pm}$, fitting the dilepton mass spectra in several exclusive signal regions. No significant evidence of a signal is observed and corresponding limits on the production cross-section are derived at 95% confidence level. The observed lower limit on the mass of a doubly charged Higgs boson only coupling to left-handed leptons ($e$,$\mu$) varies from 770 GeV to 870 GeV (850 GeV expected) for $B(H^{\pm\pm}\rightarrow \ell^{\pm}\ell^{\pm})$ = 100% and both the expected and observed mass limits are above 450 GeV for $B(H^{\pm\pm}\rightarrow \ell^{\pm}\ell^{\pm})$ = 10% and any combination of partial branching ratios.

Results of a search for physics beyond the Standard Model in events containing an energetic photon and large missing transverse momentum with the ATLAS detector at the Large Hadron Collider are reported. As the number of events observed in data, corresponding to an integrated luminosity of 36.1 $\textrm fb^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV, is in agreement with the Standard Model expectations, exclusion limits in models where dark-matter candidates are pair-produced are determined. For dark-matter production via an axial-vector or a vector mediator in the s-channel, this search excludes mediator masses below 750-1200 GeV for dark-matter candidate masses below 230-480 GeV at 95% confidence level, depending on the couplings. In an effective theory of dark-matter production, the limits restrict the value of the suppression scale $M_{*}$ to be above 790 GeV at 95% confidence level. A limit is also reported on the production of a high-mass scalar resonance by processes beyond the Standard Model, in which the resonance decays to $Z\gamma$ and the Z boson subsequently decays into neutrinos.

The distributions of transverse momentum and longitudinal momentum fraction of charged particles in jets are measured in Pb+Pb and pp collisions with the ATLAS detector at the LHC. The distributions are measured as a function of jet transverse momentum and rapidity. The analysis utilises an integrated luminosity of 0.14 nb$^{-1}$ of Pb+Pb data and 4.0 pb$^{-1}$ of pp data collected in 2011 and 2013, respectively, at the same centre-of-mass energy of 2.76 TeV per colliding nucleon pair. The distributions measured in pp collisions are used as a reference for those measured in Pb+Pb collisions in order to evaluate the impact on the internal structure of jets from the jet energy loss of fast partons propagating through the hot, dense medium created in heavy-ion collisions. Modest but significant centrality-dependent modifications of fragmentation functions in Pb+Pb collisions with respect to those in pp collisions are seen. No significant dependence of modifications on jet $p_{\mathrm{T}}$ and rapidity selections is observed except for the fragments with the highest transverse momenta for which some reduction of yields is observed for more forward jets.

Measurements of differential cross sections of top quark pair production in association with jets by the ATLAS experiment at the LHC are presented. The measurements are performed as functions of the top quark transverse momentum, the transverse momentum of the top quark-antitop quark system and the out-of-plane transverse momentum using data from $pp$ collisions at $\sqrt{s}=13$ TeV collected by the ATLAS detector at the LHC in 2015 and corresponding to an integrated luminosity of 3.2 fb$^{-1}$. The top quark pair events are selected in the lepton (electron or muon) + jets channel. The measured cross sections, which are compared to several predictions, allow a detailed study of top quark production.

The results of a search for the direct pair production of top squarks, the supersymmetric partner of the top quark, in final states with one isolated electron or muon, several energetic jets, and missing transverse momentum are reported. The analysis also targets spin-0 mediator models, where the mediator decays into a pair of dark-matter particles and is produced in association with a pair of top quarks. The search uses data from proton-proton collisions delivered by the Large Hadron Collider in 2015 and 2016 at a centre-of-mass energy of $\sqrt{s}=13$ TeV and recorded by the ATLAS detector, corresponding to an integrated luminosity of 36 fb$^{-1}$. A wide range of signal scenarios with different mass-splittings between the top squark, the lightest neutralino and possible intermediate supersymmetric particles are considered, including cases where the W bosons or the top quarks produced in the decay chain are off-shell. No significant excess over the Standard Model prediction is observed. The null results are used to set exclusion limits at 95% confidence level in several supersymmetry benchmark models. For pair-produced top-squarks decaying into top quarks, top-squark masses up to 940 GeV are excluded. Stringent exclusion limits are also derived for all other considered top-squark decay scenarios. For the spin-0 mediator models, upper limits are set on the visible cross-section.

A detailed study of multi-particle azimuthal correlations is presented using $pp$ data at $\sqrt{s}=5.02$ and 13 TeV, and $p$+Pb data at $\sqrt{s_{\rm{NN}}}=5.02$ TeV, recorded with the ATLAS detector at the LHC. The azimuthal correlations are probed using four-particle cumulants $c_{n}\{4\}$ and flow coefficients $v_n\{4\}=(-c_{n}\{4\})^{1/4}$ for $n=2$ and 3, with the goal of extracting long-range multi-particle azimuthal correlation signals and suppressing the short-range correlations. The values of $c_{n}\{4\}$ are obtained as a function of the average number of charged particles per event, $\left\langle N_{\rm{ch}} \right\rangle$, using the recently proposed two-subevent and three-subevent cumulant methods, and compared with results obtained with the standard cumulant method. The three-subevent method is found to be least sensitive to short-range correlations, which originate mostly from jets with a positive contribution to $c_{n}\{4\}$. The three-subevent method gives a negative $c_{2}\{4\}$, and therefore a well-defined $v_2\{4\}$, nearly independent of $\left\langle N_{\rm{ch}} \right\rangle$, which provides direct evidence that the long-range multi-particle azimuthal correlations persist to events with low multiplicity. Furthermore, $v_2\{4\}$ is found to be smaller than the $v_2\{2\}$ measured using the two-particle correlation method, as expected for long-range collective behavior. Finally, the measured values of $v_2\{4\}$ and $v_2\{2\}$ are used to estimate the number of sources relevant for the initial eccentricity in the collision geometry.

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.

A search for the supersymmetric partners of the Standard Model bottom and top quarks is presented. The search uses 36.1 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13$ TeV collected by the ATLAS experiment at the Large Hadron Collider. Direct production of pairs of bottom and top squarks ($\tilde{b}_{1}$ and $\tilde{t}_{1}$) is searched for in final states with $b$-tagged jets and missing transverse momentum. Distinctive selections are defined with either no charged leptons (electrons or muons) in the final state, or one charged lepton. The zero-lepton selection targets models in which the $\tilde{b}_{1}$ is the lightest squark and decays via $\tilde{b}_{1} \rightarrow b \tilde{\chi}^{0}_{1}$, where $\tilde{\chi}^{0}_{1}$ is the lightest neutralino. The one-lepton final state targets models where bottom or top squarks are produced and can decay into multiple channels, $\tilde{b}_{1} \rightarrow b \tilde{\chi}^{0}_{1}$ and $\tilde{b}_{1} \rightarrow t \tilde{\chi}^{\pm}_{1}$, or $\tilde{t}_{1} \rightarrow t \tilde{\chi}^{0}_{1}$ and $\tilde{t}_{1} \rightarrow b \tilde{\chi}^{\pm}_{1}$, where $\tilde{\chi}^{\pm}_{1}$ is the lightest chargino and the mass difference $m_{\tilde{\chi}^{\pm}_{1}}- m_{\tilde{\chi}^{0}_{1}}$ is set to 1 GeV. No excess above the expected Standard Model background is observed. Exclusion limits at 95\% confidence level on the mass of third-generation squarks are derived in various supersymmetry-inspired simplified models.

A search for direct top squark pair production resulting in events with either a same-flavour opposite-sign dilepton pair with invariant mass compatible with a $Z$ boson or a pair of jets compatible with a Standard Model (SM) Higgs boson ($h$) is presented. Requirements on additional leptons, jets, jets identified as originating from $b$-quarks, and missing transverse momentum are imposed to target the other decay products of the top squark pair. The analysis is performed using proton-proton collision data at $\sqrt{s}=13$ TeV collected with the ATLAS detector at the LHC in 2015--2016, corresponding to an integrated luminosity of 36.1 fb$^{-1}$. No excess is observed in the data with respect to the SM predictions. The results are interpreted in two sets of models. In the first set, direct production of pairs of lighter top squarks ($\tilde{t}_1$) with long decay chains involving $Z$ or Higgs bosons is considered. The second set includes direct pair production of the heavier top squark pairs ($\tilde{t}_2$) decaying via $\tilde{t}_2 \rightarrow Z \tilde{t}_1 $ or $\tilde{t}_2 \rightarrow h \tilde{t}_1$. The results exclude at 95\% confidence level $\tilde{t}_2$ and $\tilde{t}_1$ masses up to about 800 GeV, extending the exclusion region of supersymmetric parameter space covered by previous LHC searches.

Results from a search for supersymmetry in events with four or more charged leptons (electrons, muons and taus) are presented. The analysis uses a data sample corresponding to 36.1 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 taus are designed to target a range of supersymmetric scenarios that can be either enriched in or depleted of events involving the production and decay of a $Z$ boson. Data yields are consistent with Standard Model expectations and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of General Gauge Mediated supersymmetry, where higgsino masses are excluded up to 295 GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of 1.46 TeV, 1.06 TeV, and 2.25 TeV are placed on wino, slepton and gluino masses, respectively.

A search for supersymmetric partners of top quarks decaying as $\tilde{t}_1\to c\tilde\chi^0_1$ and supersymmetric partners of charm quarks decaying as $\tilde{c}_1\to c\tilde\chi^0_1$, where $\tilde\chi^0_1$ is the lightest neutralino, is presented. The search uses 36.1 ${\rm fb}^{-1}$ $pp$ collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS experiment at the Large Hadron Collider and is performed in final states with jets identified as containing charm hadrons. Assuming a 100% branching ratio to $c\tilde\chi^0_1$, top and charm squarks with masses up to 850 GeV are excluded at 95% confidence level for a massless lightest neutralino. For $m_{\tilde{t}_1,\tilde{c}_1}-m_{\tilde\chi^0_1}

A combination of the searches for pair-produced vector-like partners of the top and bottom quarks in various decay channels ($T$$\rightarrow$$Zt/Wb/Ht$, $B$$\rightarrow$$Zb/Wt/Hb$) is performed using 36.1 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The observed data are found to be in good agreement with the Standard Model background prediction in all individual searches. Therefore, combined 95% confidence-level upper limits are set on the production cross-section for a range of vector-like quark scenarios, significantly improving upon the reach of the individual searches. Model-independent limits are set assuming the vector-like quarks decay to Standard Model particles. A singlet $T$ is excluded for masses below 1.31 TeV and a singlet $B$ is excluded for masses below 1.22 TeV. Assuming a weak isospin $(T,B)$ doublet and $|V_{Tb}| \ll |V_{tB}|$, $T$ and $B$ masses below 1.37 TeV are excluded.

A search for vectorlike quarks is presented, which targets their decay into a $Z$ boson and a third-generation Standard Model quark. In the case of a vectorlike quark $T$ ($B$) with charge $+2/3e$ ($-1/3e$), the decay searched for is $T \rightarrow Zt$ ($B \rightarrow Zb$). Data for this analysis were taken during 2015 and 2016 with the ATLAS detector at the Large Hadron Collider and correspond to an integrated luminosity of 36.1 fb$^{-1}$ of $pp$ collisions at $\sqrt{s} = 13$ TeV. The final state used is characterized by the presence of $b$-tagged jets, as well as a $Z$ boson with high transverse momentum, which is reconstructed from a pair of opposite-sign same-flavor leptons. Pair and single production of vectorlike quarks are both taken into account and are each searched for using optimized dileptonic exclusive and trileptonic inclusive event selections. In these selections, the high scalar sum of jet transverse momenta, the presence of high-transverse-momentum large-radius jets, as well as - in the case of the single-production selections - the presence of forward jets are used. No significant excess over the background-only hypothesis is found and exclusion limits at 95% confidence level allow masses of vectorlike quarks of $m_T > 1030$ GeV ($m_T > 1210$ GeV) and $m_B > 1010$ GeV ($m_B > 1140$ GeV) in the singlet (doublet) model. In the case of 100% branching ratio for $T\rightarrow Zt$ ($B\rightarrow Zb$), the limits are $m_T > 1340$ GeV ($m_B > 1220$ GeV). Limits at 95% confidence level are also set on the coupling to Standard Model quarks for given vectorlike quark masses.