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.

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.

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.

The production of $Z$ bosons in association with a high-energy photon ($Z\gamma$ production) is studied in the neutrino decay channel of the $Z$ boson using $pp$ collisions at $\sqrt{s}$ = 13 TeV. The analysis uses a data sample with an integrated luminosity of 36.1 fb$^{-1}$ collected by the ATLAS detector at the LHC in 2015 and 2016. Candidate $Z\gamma$ events with invisible decays of the $Z$ boson are selected by requiring significant transverse momentum ($p_{T}$) of the dineutrino system in conjunction with a single isolated photon with large transverse energy ($E_{T}$). The rate of $Z\gamma$ production is measured as a function of photon $E_{T}$, dineutrino system $p_{T}$ and jet multiplicity. Evidence of anomalous triple gauge-boson couplings is sought in $Z\gamma$ production with photon $E_{T}$ greater than 600 GeV. No excess is observed relative to the Standard Model expectation, and upper limits are set on the strength of $ZZ\gamma$ and $Z\gamma\gamma$ couplings.

A search for heavy resonances decaying into a $W$ or $Z$ boson and a Higgs boson produced in proton$-$proton collisions at the Large Hadron Collider at $\sqrt{s} = 13$ TeV is presented. The analysis utilizes the dominant $W \to q \bar{q}^\prime$ or $Z \to q \bar{q}$ and $H \to b \bar{b}$ decays with substructure techniques applied to large-radius jets. A sample corresponding to an integrated luminosity of 139 fb$^{-1}$ collected with the ATLAS detector is analyzed and no significant excess of data is observed over the background prediction. The results are interpreted in the context of the Heavy Vector Triplet model with spin-1 $W^\prime$ and $Z^\prime$ bosons. Upper limits on the cross section are set for resonances with mass between 1.5 and 5.0 TeV, ranging from 6.8 to 0.53 fb for $W^\prime \to WH$ and from 8.7 to 0.53 fb for $Z^\prime \to ZH$ at the 95 % confidence level.

A search for new resonances decaying into jets containing b-hadrons in $pp$ collisions with the ATLAS detector at the LHC is presented in the dijet mass range from 0.57 TeV to 7 TeV. The dataset corresponds to an integrated luminosity of up to 36.1 fb$^{-1}$ collected in 2015 and 2016 at $\sqrt{s} = 13$ TeV. No evidence of a significant excess of events above the smooth background shape is found. Upper cross-section limits and lower limits on the corresponding signal mass parameters for several types of signal hypotheses are provided at 95% CL. In addition, 95% CL upper limits are set on the cross-sections for new processes that would produce Gaussian-shaped signals in the di-b-jet mass distributions.

A search for resonant and non-resonant pair production of Higgs bosons in the $b\bar{b}\tau^+\tau^-$ final state is presented. The search uses 36.1 fb$^{-1}$ of $pp$ collision data with $\sqrt{s}= 13$ TeV recorded by the ATLAS experiment at the LHC in 2015 and 2016. The semileptonic and fully hadronic decays of the $\tau$-lepton pair are considered. No significant excess above the expected background is observed in the data. The cross-section times branching ratio for non-resonant Higgs boson pair production is constrained to be less than 30.9 fb, 12.7 times the Standard Model expectation, at 95% confidence level. The data are also analyzed to probe resonant Higgs boson pair production, constraining a model with an extended Higgs sector based on two doublets and a Randall-Sundrum bulk graviton model. Upper limits are placed on the resonant Higgs boson pair production cross-section times branching ratio, excluding resonances $X$ in the mass range $305~{\rm GeV} < m_X < 402~{\rm GeV}$ in the simplified hMSSM minimal supersymmetric model for $\tan\beta=2$ and excluding bulk Randall-Sundrum gravitons $G_{\mathrm{KK}}$ in the mass range $325~{\rm GeV} < m_{G_{\mathrm{KK}}} < 885~{\rm GeV}$ for $k/\overline{M}_{\mathrm{Pl}} = 1$.

Searches for scalar leptoquarks pair-produced in proton-proton collisions at $\sqrt{s}=13$ TeV at the Large Hadron Collider are performed by the ATLAS experiment. A data set corresponding to an integrated luminosity of 36.1 fb$^{-1}$ is used. Final states containing two electrons or two muons and two or more jets are studied, as are states with one electron or muon, missing transverse momentum and two or more jets. No statistically significant excess above the Standard Model expectation is observed. The observed and expected lower limits on the leptoquark mass at 95% confidence level extend up to 1.29 TeV and 1.23 TeV for first- and second-generation leptoquarks, respectively, as postulated in the minimal Buchm\"uller-R\"uckl-Wyler model, assuming a branching ratio into a charged lepton and a quark of 50%. In addition, measurements of particle-level fiducial and differential cross sections are presented for the $Z\rightarrow ee$, $Z\rightarrow\mu\mu$ and $t\bar{t}$ processes in several regions related to the search control regions. Predictions from a range of generators are compared with the measurements, and good agreement is seen for many of the observables. However, the predictions for the $Z\rightarrow\ell\ell$ measurements in observables sensitive to jet energies disagree with the data.

A search for new charged massive gauge bosons, $W^\prime$, is performed with the ATLAS detector at the LHC. Data were collected in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV and correspond to an integrated luminosity of 36.1 $\textrm{fb}^{-1}$. This analysis searches for $W^\prime$ bosons in the $W^\prime \rightarrow t\bar{b}$ decay channel in final states with an electron or muon plus jets. The search covers resonance masses between 0.5 and 5.0 TeV and considers right-handed $W^\prime$ bosons. No significant deviation from the Standard Model (SM) expectation is observed and upper limits are set on the $W^\prime \rightarrow t\bar{b}$ cross section times branching ratio and the $W^\prime$ boson effective couplings as a function of the $W^\prime$ boson mass. For right-handed $W^\prime$ bosons with coupling to the SM particles equal to the SM weak coupling constant, masses below 3.15 TeV are excluded at the 95% confidence level. This search is also combined with a previously published ATLAS result for $W^\prime \rightarrow t\bar{b}$ in the fully hadronic final state. Using the combined searches, right-handed $W^\prime$ bosons with masses below 3.25 TeV are excluded at the 95% confidence level.

A search for a chargino$-$neutralino pair decaying via the 125 GeV Higgs boson into photons is presented. The study is based on the data collected between 2015 and 2018 with the ATLAS detector at the LHC, corresponding to an integrated luminosity of 139 fb$^{-1}$ of $pp$ collisions at a centre-of-mass energy of 13 TeV. No significant excess over the expected background is observed. Upper limits at 95% confidence level for a massless $\tilde{\chi}^{0}_{1}$ are set on several electroweakino production cross-sections and the visible cross-section for beyond the Standard Model processes. In the context of simplified supersymmetric models, 95% confidence-level limits of up to 310 GeV in $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$, where $m(\tilde{\chi}^{0}_{1})=0.5$ GeV, are set. Limits at 95% confidence level are also set on the $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ cross-section in the mass plane of $m(\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2})$ and $m(\tilde{\chi}^{0}_{1})$, and on scenarios with gravitino as the lightest supersymmetric particle. Upper limits at the 95% confidence-level are set on the higgsino production cross-section. Higgsino masses below 380 GeV are excluded for the case of the higgsino fully decaying into a Higgs boson and a gravitino.

Results of a search for new particles decaying into eight or more jets and moderate missing transverse momentum are presented. The analysis uses 139 fb$^{-1}$ of proton$-$proton collision data at $\sqrt{s} = 13$ TeV collected by the ATLAS experiment at the Large Hadron Collider between 2015 and 2018. The selection rejects events containing isolated electrons or muons, and makes requirements according to the number of $b$-tagged jets and the scalar sum of masses of large-radius jets. The search extends previous analyses both in using a larger dataset and by employing improved jet and missing transverse momentum reconstruction methods which more cleanly separate signal from background processes. No evidence for physics beyond the Standard Model is found. The results are interpreted in the context of supersymmetry-inspired simplified models, significantly extending the limits on the gluino mass in those models. In particular, limits on the gluino mass are set at 2 TeV when the lightest neutralino is nearly massless in a model assuming a two-step cascade decay via the lightest chargino and second-lightest neutralino.

A search for electroweak production of supersymmetric particles is performed in two-lepton and three-lepton final states using recursive jigsaw reconstruction. The search uses data collected in 2015 and 2016 by the ATLAS experiment in $\sqrt{s}$ = 13 TeV proton--proton collisions at the CERN Large Hadron Collider corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Chargino-neutralino pair production, with decays via W/Z bosons, is studied in final states involving leptons and jets and missing transverse momentum for scenarios with large and intermediate mass-splittings between the parent particle and lightest supersymmetric particle, as well as for the scenario where this mass splitting is close to the mass of the Z boson. The latter case is challenging since the vector bosons are produced with kinematic properties that are similar to those in Standard Model processes. Results are found to be compatible with the Standard Model expectations in the signal regions targeting large and intermediate mass-splittings, and chargino-neutralino masses up to 600 GeV are excluded at 95% confidence level for a massless lightest supersymmetric particle. Excesses of data above the expected background are found in the signal regions targeting low mass-splittings, and the largest local excess amounts to 3.0 standard deviations.

A search for flavour-changing neutral current (FCNC) events via the coupling of a top quark, a photon, and an up or charm quark is presented using 81 fb$^{-1}$ of proton-proton collision data taken at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC. Events with a photon, an electron or muon, a $b$-tagged jet, and missing transverse momentum are selected. A neural network based on kinematic variables differentiates between events from signal and background processes. The data are consistent with the background-only hypothesis, and limits are set on the strength of the $tq\gamma$ coupling in an effective field theory. These are also interpreted as 95% CL upper limits on the cross section for FCNC $t\gamma$ production via a left-handed (right-handed) $tu\gamma$ coupling of 36 fb (78 fb) and on the branching ratio for $t\rightarrow \gamma u$ of $2.8\times 10^{-5}$ ($6.1\times 10^{-5}$). In addition, they are interpreted as 95% CL upper limits on the cross section for FCNC $t\gamma$ production via a left-handed (right-handed) $tc\gamma$ coupling of 40 fb (33 fb) and on the branching ratio for $t\rightarrow \gamma c$ of $22\times 10^{-5}$ ($18\times 10^{-5}$).

A search for long-lived particles decaying into hadrons and at least one muon is presented. The analysis selects events that pass a muon or missing-transverse-momentum trigger and contain a displaced muon track and a displaced vertex. The analyzed dataset of proton-proton collisions at $\sqrt{s} = 13$ TeV was collected with the ATLAS detector and corresponds to 136 fb$^{-1}$. The search employs dedicated reconstruction techniques that significantly increase the sensitivity to long-lived particle decays that occur 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 presented as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model, and interpreted as exclusion limits in scenarios with pair-production of long-lived top squarks that decay via a small $R$-parity-violating coupling into a quark and a muon. Top squarks with masses up to 1.7 TeV are excluded for a lifetime of 0.1 ns, and masses below 1.3 TeV are excluded for lifetimes between 0.01 ns and 30 ns.

The inclusive top quark pair ($t\bar{t}$) production cross-section $\sigma_{t\bar{t}}$ has been measured in proton$-$proton collisions at $\sqrt{s}=13$ TeV, using $36.1$ fb$^{-1}$ of data collected in 2015$-$16 by the ATLAS experiment at the LHC. Using events with an opposite-charge $e\mu$ pair and $b$-tagged jets, the cross-section is measured to be: \begin{equation}\nonumber \sigma_{t\bar{t}} = 826.4 \pm 3.6\,\mathrm{(stat)}\ \pm 11.5\,\mathrm{(syst)}\ \pm 15.7\,\mathrm{(lumi)}\ \pm 1.9\,\mathrm{(beam)}\,\mathrm{pb}, \end{equation} where the uncertainties reflect the limited size of the data sample, experimental and theoretical systematic effects, the integrated luminosity, and the LHC beam energy, giving a total uncertainty of 2.4%. The result is consistent with theoretical QCD calculations at next-to-next-to-leading order. It is used to determine the top quark pole mass via the dependence of the predicted cross-section on $m_t^{\mathrm{pole}}$, giving $m_t^{\mathrm{pole}}=173.1^{+2.0}_{-2.1}$ GeV. It is also combined with measurements at $\sqrt{s}=7$ TeV and $\sqrt{s}=8$ TeV to derive ratios and double ratios of $t\bar{t}$ and $Z$ cross-sections at different energies. The same event sample is used to measure absolute and normalised differential cross-sections as functions of single-lepton and dilepton kinematic variables, and the results compared with predictions from various Monte Carlo event generators.

A search for Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is carried out with up to 36.1 $\mathrm{fb}^{-1}$ of LHC proton--proton collision data collected at $\sqrt{s}$ = 13 TeV with the ATLAS detector in 2015 and 2016. Three benchmark signals are studied: a spin-2 graviton decaying into a Higgs boson pair, a scalar resonance decaying into a Higgs boson pair, and Standard Model non-resonant Higgs boson pair production. Two analyses are carried out, each implementing a particular technique for the event reconstruction that targets Higgs bosons reconstructed as pairs of jets or single boosted jets. The resonance mass range covered is 260--3000 GeV. The analyses are statistically combined and upper limits on the production cross section of Higgs boson pairs times branching ratio to $b\bar{b}b\bar{b}$ are set in each model. No significant excess is observed; the largest deviation of data over prediction is found at a mass of 280 GeV, corresponding to 2.3 standard deviations globally. The observed 95% confidence level upper limit on the non-resonant production is 13 times the Standard Model prediction.

This Letter presents a search for the production of a long-lived neutral particle ($Z_d$) decaying within the ATLAS hadronic calorimeter, in association with a Standard Model (SM) $Z$ boson produced via an intermediate scalar boson, where $Z\to l^+l^-$ ($l=e,\mu$). The data used were collected by the ATLAS detector during 2015 and 2016 $pp$ collisions with a center-of-mass energy of $\sqrt{s} = 13$ TeV at the Large Hadron Collider and corresponds to an integrated luminosity of $36.1\pm0.8$ fb$^{-1}$. No significant excess of events is observed above the expected background. Limits on the production cross section of the scalar boson times its decay branching fraction into the long-lived neutral particle are derived as a function of the mass of the intermediate scalar boson, the mass of the long-lived neutral particle, and its $c\tau$ from a few centimeters to one hundred meters. In the case that the intermediate scalar boson is the SM Higgs boson, its decay branching fraction to a long-lived neutral particle with a $c\tau$ approximately between 0.1 m and 7 m is excluded with a 95% confidence level up to 10% for $m_{Z_d}$ between 5 and 15 GeV.

A search for supersymmetry in events with large missing transverse momentum, jets, and at least one hadronically decaying $\tau$-lepton is presented. Two exclusive final states with either exactly one or at least two $\tau$-leptons are considered. The analysis is based on proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ delivered by the Large Hadron Collider and recorded by the ATLAS detector in 2015 and 2016. No significant excess is observed over the Standard Model expectation. At 95% confidence level, model-independent upper limits on the cross section are set and exclusion limits are provided for two signal scenarios: a simplified model of gluino pair production with $\tau$-rich cascade decays, and a model with gauge-mediated supersymmetry breaking (GMSB). In the simplified model, gluino masses up to 2000 GeV are excluded for low values of the mass of the lightest supersymmetric particle (LSP), while LSP masses up to 1000 GeV are excluded for gluino masses around 1400 GeV. In the GMSB model, values of the supersymmetry-breaking scale are excluded below 110 TeV for all values of $\tan\beta$ in the range $2 \leq \tan\beta \leq 60$, and below 120 TeV for $\tan\beta>30$.

A search for a narrow scalar resonance decaying into an opposite-sign muon pair produced in events with and without $b$-tagged jets is presented in this paper. The search uses 36.1 fb$^{-1}$ of $\sqrt{s}= 13$ TeV proton-proton collision data recorded by the ATLAS experiment at the LHC. No significant excess of events above the expected Standard Model background is observed in the investigated mass range of 0.2 to 1.0 TeV. The observed upper limits at 95$\%$ confidence level on the cross section times branching ratio for $b$-quark associated production and gluon-gluon fusion are between 1.9 and 41 fb and 1.6 and 44 fb respectively, which is consistent with expectations.

A search for a right-handed gauge boson $W_{\mathrm{R}}$, decaying into a boosted right-handed heavy neutrino $N_{\mathrm{R}}$, in the framework of Left-Right Symmetric Models is presented. It is based on data from proton-proton collisions with a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider during the years 2015, 2016 and 2017, corresponding to an integrated luminosity of 80 fb$^{-1}$. The search is performed separately for electrons and muons in the final state. A distinguishing feature of the search is the use of large-radius jets containing electrons. Selections based on the signal topology result in smaller background compared with to expected signal. No significant deviation from the Standard Model prediction is observed and lower limits are set in the $W_{\mathrm{R}}$ and $N_{\mathrm{R}}$ mass plane. Mass values of the $W_{\mathrm{R}}$ smaller than 3.8-5 TeV are excluded for $N_{\mathrm{R}}$ in the mass range 0.1-1.8 TeV.

This paper presents measurements of $t\bar{t}$ production in association with additional $b$-jets in $pp$ collisions at the LHC at a centre-of-mass energy of 13 TeV. The data were recorded with the ATLAS detector and correspond to an integrated luminosity of 36.1 fb$^{-1}$. Fiducial cross-section measurements are performed in the dilepton and lepton-plus-jets $t\bar{t}$ decay channels. Results are presented at particle level in the form of inclusive cross-sections of $t\bar{t}$ final states with three and four $b$-jets as well as differential cross-sections as a function of global event properties and properties of $b$-jet pairs. The measured inclusive fiducial cross-sections generally exceed the $t\bar{t}b\bar{b}$ predictions from various next-to-leading-order matrix element calculations matched to a parton shower but are compatible within the total uncertainties. The experimental uncertainties are smaller than the uncertainties in the predictions. Comparisons of state-of-the-art theoretical predictions with the differential measurements are shown and good agreement with data is found for most of them.

This paper reports on a search for the electroweak diboson ($WW/WZ/ZZ$) production in association with a high-mass dijet system, using data from proton-proton collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV. The data, corresponding to an integrated luminosity of 35.5 fb$^{-1}$, were recorded with the ATLAS detector in 2015 and 2016 at the Large Hadron Collider. The search is performed in final states in which one boson decays leptonically, and the other boson decays hadronically. The hadronically decaying $W/Z$ boson is reconstructed as either two small-radius jets or one large-radius jet using jet substructure techniques. The electroweak production of $WW/WZ/ZZ$ in association with two jets is measured with an observed (expected) significance of 2.7 (2.5) standard deviations, and the fiducial cross section is measured to be $45.1 \pm 8.6(\mathrm{stat.}) ^{+15.9} _{-14.6} (\mathrm{syst.})$ fb.

A search for Higgs boson pair production via vector-boson fusion (VBF) in the $b\bar{b}b\bar{b}$ final state is carried out with the ATLAS experiment using 126 fb$^{-1}$ of proton-proton collision data delivered at $\sqrt{s} = 13$ TeV by the Large Hadron Collider. This search is sensitive to VBF production of additional heavy bosons that may decay into Higgs boson pairs, and in a non-resonant topology it can constrain the quartic coupling between the Higgs bosons and vector bosons. No significant excess relative to the Standard Model expectation is observed, and limits on the production cross-section are set at the 95 % confidence level for a heavy scalar resonance in the context of an extended Higgs sector, and for non-resonant Higgs boson pair production. Interpretation in terms of the coupling between a Higgs boson pair and two vector bosons is also provided: coupling values normalised to the Standard Model expectation of $\kappa_{2V} < -0.43$ and $\kappa_{2V} > 2.56$ are excluded at the 95 % confidence level in data.