A measurement of beauty hadron production at mid-rapidity in proton-lead collisions at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}}=5.02$ TeV is presented. The semi-inclusive decay channel of beauty hadrons into J/$\psi$ is considered, where the J/$\psi$ mesons are reconstructed in the dielectron decay channel at mid-rapidity down to transverse momenta of 1.3 GeV/$c$. The ${\rm {b\overline{b}}}$ production cross section at mid-rapidity, ${\rm d}\sigma_{\rm {b\overline{b}}}/{\rm d} y$, and the total cross section extrapolated over full phase space, $\sigma_{\rm {b\overline{b}}}$, are obtained. This measurement is combined with results on inclusive J/$\psi$ production to determine the prompt J/$\psi$ cross sections. The results in p-Pb collisions are then scaled to expectations from pp collisions at the same centre-of-mass energy to derive the nuclear modification factor $R_{\rm pPb}$, and compared to models to study possible nuclear modifications of the production induced by cold nuclear matter effects. $R_{\rm pPb}$ is found to be smaller than unity at low $p_{\rm T}$ for both J/$\psi$ coming from beauty hadron decays and prompt J/$\psi$.

A search for new phenomena has been performed in final states with at least one isolated high-momentum photon, jets and missing transverse momentum in proton--proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV. The data, collected by the ATLAS experiment at the CERN LHC, correspond to an integrated luminosity of 139 $fb^{-1}$. The experimental results are interpreted in a supersymmetric model in which pair-produced gluinos decay into neutralinos, which in turn decay into a gravitino, at least one photon, and jets. No significant deviations from the predictions of the Standard Model are observed. Upper limits are set on the visible cross section due to physics beyond the Standard Model, and lower limits are set on the masses of the gluinos and neutralinos, all at 95% confidence level. Visible cross sections greater than 0.022 fb are excluded and pair-produced gluinos with masses up to 2200 GeV are excluded for most of the NLSP masses investigated.

This search, a type not previously performed at ATLAS, uses a comparison of the production cross sections for $e^+ \mu^-$ and $e^- \mu^+$ pairs to constrain physics processes beyond the Standard Model. It uses $139 \text{fb}^{-1}$ of proton$-$proton collision data recorded at $\sqrt{s} = 13$ TeV at the LHC. Targeting sources of new physics which prefer final states containing $e^{+}\mu^{-}$ to $e^{-}\mu^{+}$, the search contains two broad signal regions which are used to provide model-independent constraints on the ratio of cross sections at the 2% level. The search also has two special selections targeting supersymmetric models and leptoquark signatures. Observations using one of these selections are able to exclude, at 95% confidence level, singly produced smuons with masses up to 640 GeV in a model in which the only other light sparticle is a neutralino when the $R$-parity-violating coupling $\lambda'_{231}$ is close to unity. Observations using the other selection exclude scalar leptoquarks with masses below 1880 GeV when $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=1$, at 95% confidence level. The limit on the coupling reduces to $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=0.46$ for a mass of 1420 GeV.

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.

A search is presented for pair-production of long-lived neutral particles using 33 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collision data, collected during 2016 by the ATLAS detector at the LHC. This search focuses on a topology in which one long-lived particle decays in the ATLAS inner detector and the other decays in the muon spectrometer. Special techniques are employed to reconstruct the displaced tracks and vertices in the inner detector and in the muon spectrometer. One event is observed that passes the full event selection, which is consistent with the estimated background. Limits are placed on scalar boson propagators with masses from 125 GeV to 1000 GeV decaying into pairs of long-lived hidden-sector scalars with masses from 8 GeV to 400 GeV. The limits placed on several low-mass scalars extend previous exclusion limits in the range of proper lifetimes $c \tau$ from 5 cm to 1 m.

A search for supersymmetry through the pair production of electroweakinos with mass splittings near the electroweak scale and decaying via on-shell $W$ and $Z$ bosons is presented for a three-lepton final state. The analyzed proton-proton collision data taken at a center-of-mass energy of $\sqrt{s}$ = 13 TeV were collected between 2015 and 2018 by the ATLAS experiment at the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb$^{-1}$. A search, emulating the recursive jigsaw reconstruction technique with easily reproducible laboratory-frame variables, is performed. The two excesses observed in the 2015-2016 data recursive jigsaw analysis in the low-mass three-lepton phase space are reproduced. Results with the full dataset are in agreement with the Standard Model expectations. They are interpreted to set exclusion limits at 95% confidence level on simplified models of chargino-neutralino pair production for masses up to 345 GeV.

A search for the direct production of the supersymmetric partners of $\tau$-leptons (staus) in final states with two hadronically decaying $\tau$-leptons is presented. The analysis uses a dataset of $pp$ collisions corresponding to an integrated luminosity of $139$ fb$^{-1}$, recorded with the ATLAS detector at the Large Hadron Collider at a center-of-mass energy of 13 TeV. No significant deviation from the expected Standard Model background is observed. Limits are derived in scenarios of direct production of stau pairs with each stau decaying into the stable lightest neutralino and one $\tau$-lepton in simplified models where the two stau mass eigenstates are degenerate. Stau masses from 120 GeV to 390 GeV are excluded at 95% confidence level for a massless lightest neutralino.

The results of a search for direct pair production of top squarks and for dark matter in events with two opposite-charge leptons (electrons or muons), jets and missing transverse momentum are reported, using 139 fb$^{-1}$ of integrated luminosity from proton-proton collisions at $\sqrt{s} = 13$ TeV, collected by the ATLAS detector at the Large Hadron Collider during Run 2 (2015-2018). This search considers the pair production of top squarks and is sensitive across a wide range of mass differences between the top squark and the lightest neutralino. Additionally, spin-0 mediator dark-matter models are considered, in which the mediator is produced in association with a pair of top quarks. The mediator subsequently decays to a pair of dark-matter particles. No significant excess of events is observed above the Standard Model background, and limits are set at 95% confidence level. The results exclude top squark masses up to about 1 TeV, and masses of the lightest neutralino up to about 500 GeV. Limits on dark-matter production are set for scalar (pseudoscalar) mediator masses up to about 250 (300) GeV.

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.

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.

The associated production of a Higgs boson and a top-quark pair is measured in events characterised by the presence of one or two electrons or muons. The Higgs boson decay into a $b$-quark pair is used. The analysed data, corresponding to an integrated luminosity of 139 fb$^{-1}$, were collected in proton-proton collisions at the Large Hadron Collider between 2015 and 2018 at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The measured signal strength, defined as the ratio of the measured signal yield to that predicted by the Standard Model, is $0.35^{+0.36}_{-0.34}$. This result is compatible with the Standard Model prediction and corresponds to an observed (expected) significance of 1.0 (2.7) standard deviations. The signal strength is also measured differentially in bins of the Higgs boson transverse momentum in the simplified template cross-section framework, including a bin for specially selected boosted Higgs bosons with transverse momentum above 300 GeV.

Cross-sections for the production of a $Z$ boson in association with two photons are measured in proton$-$proton collisions at a centre-of-mass energy of 13 TeV. The data used correspond to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS experiment during Run 2 of the LHC. The measurements use the electron and muon decay channels of the $Z$ boson, and a fiducial phase-space region where the photons are not radiated from the leptons. The integrated $Z(\rightarrow\ell\ell)\gamma\gamma$ cross-section is measured with a precision of 12% and differential cross-sections are measured as a function of six kinematic variables of the $Z\gamma\gamma$ system. The data are compared with predictions from MC event generators which are accurate to up to next-to-leading order in QCD. The cross-section measurements are used to set limits on the coupling strengths of dimension-8 operators in the framework of an effective field theory.

Searches are performed for nonresonant and resonant di-Higgs boson production in the $b\bar{b}\gamma\gamma$ final state. The data set used corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. No excess above the expected background is found and upper limits on the di-Higgs boson production cross sections are set. A 95% confidence-level upper limit of 4.2 times the cross section predicted by the Standard Model is set on $pp \rightarrow HH$ nonresonant production, where the expected limit is 5.7 times the Standard Model predicted value. The expected constraints are obtained for a background hypothesis excluding $pp \rightarrow HH$ production. The observed (expected) constraints on the Higgs boson trilinear coupling modifier $\kappa_{\lambda}$ are determined to be $[-1.5, 6.7]$ $([-2.4, 7.7])$ at 95% confidence level, where the expected constraints on $\kappa_{\lambda}$ are obtained excluding $pp \rightarrow HH$ production from the background hypothesis. For resonant production of a new hypothetical scalar particle $X$ ($X \rightarrow HH \rightarrow b\bar{b}\gamma\gamma$), limits on the cross section for $pp \to X \to HH$ are presented in the narrow-width approximation as a function of $m_{X}$ in the range $251 \leq m_{X} \leq 1000$ GeV. The observed (expected) limits on the cross section for $pp \to X \to HH$ range from 640 fb to 44 fb (391 fb to 46 fb) over the considered mass range.

A search is presented for displaced production of Higgs bosons or $Z$ bosons, originating from the decay of a neutral long-lived particle (LLP) and reconstructed in the decay modes $H\rightarrow \gamma\gamma$ and $Z\rightarrow ee$. The analysis uses the full Run 2 data set of proton$-$proton collisions delivered by the LHC at an energy of $\sqrt{s}=13$ TeV between 2015 and 2018 and recorded by the ATLAS detector, corresponding to an integrated luminosity of 139 fb$^{-1}$. Exploiting the capabilities of the ATLAS liquid argon calorimeter to precisely measure the arrival times and trajectories of electromagnetic objects, the analysis searches for the signature of pairs of photons or electrons which arise from a common displaced vertex and which arrive after some delay at the calorimeter. The results are interpreted in a gauge-mediated supersymmetry breaking model with pair-produced higgsinos that decay to LLPs, and each LLP subsequently decays into either a Higgs boson or a $Z$ boson. The final state includes at least two particles that escape direct detection, giving rise to missing transverse momentum. No significant excess is observed above the background expectation. The results are used to set upper limits on the cross section for higgsino pair production, up to a $\tilde\chi^0_1$ mass of 369 (704) GeV for decays with 100% branching ratio of $\tilde\chi^0_1$ to Higgs ($Z$) bosons for a $\tilde\chi^0_1$ lifetime of 2 ns. A model-independent limit is also set on the production of pairs of photons or electrons with a significant delay in arrival at the calorimeter.

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 novel search for exotic decays of the Higgs boson into pairs of long-lived neutral particles, each decaying into a bottom quark pair, is performed using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collision data collected with the ATLAS detector at the LHC. Events consistent with the production of a Higgs boson in association with a leptonically decaying $Z$ boson are analysed. Long-lived particle (LLP) decays are reconstructed from inner-detector tracks as displaced vertices with high mass and track multiplicity relative to Standard Model processes. The analysis selection requires the presence of at least two displaced vertices, effectively suppressing Standard Model backgrounds. The residual background contribution is estimated using a data-driven technique. No excess over Standard Model predictions is observed, and upper limits are set on the branching ratio of the Higgs boson to LLPs. Branching ratios above 10% are excluded at 95% confidence level for LLP mean proper lifetimes $c\tau$ as small as 4 mm and as large as 100 mm. For LLP masses below 40 GeV, these results represent the most stringent constraint in this lifetime regime.

A search for diphoton resonances in the mass range between 10 and 70 GeV with the ATLAS experiment at the Large Hadron Collider (LHC) is presented. The analysis is based on $pp$ collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at a centre-of-mass energy of 13 TeV recorded from 2015 to 2018. Previous searches for diphoton resonances at the LHC have explored masses down to 65 GeV, finding no evidence of new particles. This search exploits the particular kinematics of events with pairs of closely spaced photons reconstructed in the detector, allowing examination of invariant masses down to 10 GeV. The presented strategy covers a region previously unexplored at hadron colliders because of the experimental challenges of recording low-energy photons and estimating the backgrounds. No significant excess is observed and the reported limits provide the strongest bound on promptly decaying axion-like particles coupling to gluons and photons for masses between 10 and 70 GeV.

A search for the electroweak production of charginos and sleptons decaying into final states with two electrons or muons is presented. The analysis is based on 139 fb$^{-1}$ of proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider at $\sqrt{s}=13$ TeV. Three $R$-parity-conserving scenarios where the lightest neutralino is the lightest supersymmetric particle are considered: the production of chargino pairs with decays via either $W$ bosons or sleptons, and the direct production of slepton pairs. The analysis is optimised for the first of these scenarios, but the results are also interpreted in the others. No significant deviations from the Standard Model expectations are observed and limits at 95 % confidence level are set on the masses of relevant supersymmetric particles in each of the scenarios. For a massless lightest neutralino, masses up to 420 GeV are excluded for the production of the lightest-chargino pairs assuming $W$-boson-mediated decays and up to 1 TeV for slepton-mediated decays, whereas for slepton-pair production masses up to 700 GeV are excluded assuming three generations of mass-degenerate sleptons.

A search for chargino$-$neutralino pair production in three-lepton final states with missing transverse momentum is presented. The study is based on a dataset of $\sqrt{s} = 13$ TeV $pp$ collisions recorded with the ATLAS detector at the LHC, corresponding to an integrated luminosity of 139 fb$^{-1}$. No significant excess relative to the Standard Model predictions is found in data. The results are interpreted in simplified models of supersymmetry, and statistically combined with results from a previous ATLAS search for compressed spectra in two-lepton final states. Various scenarios for the production and decay of charginos ($\tilde\chi^\pm_1$) and neutralinos ($\tilde\chi^0_2$) are considered. For pure higgsino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair-production scenarios, exclusion limits at 95% confidence level are set on $\tilde\chi^0_2$ masses up to 210 GeV. Limits are also set for pure wino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair production, on $\tilde\chi^0_2$ masses up to 640 GeV for decays via on-shell $W$ and $Z$ bosons, up to 300 GeV for decays via off-shell $W$ and $Z$ bosons, and up to 190 GeV for decays via $W$ and Standard Model Higgs bosons.

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 is performed for the electroweak pair production of charginos and associated production of a chargino and neutralino, each of which decays through an $R$-parity-violating coupling into a lepton and a $W$, $Z$, or Higgs boson. The trilepton invariant-mass spectrum is constructed from events with three or more leptons, targeting chargino decays that include an electron or muon and a leptonically decaying $Z$ boson. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data produced by the Large Hadron Collider at a center-of-mass energy of $\sqrt{s}$ = 13 TeV and collected by the ATLAS experiment between 2015 and 2018. The data are found to be consistent with predictions from the Standard Model. The results are interpreted as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model. Limits are also set on the production of charginos and neutralinos for a Minimal Supersymmetric Standard Model with an approximate $B$-$L$ symmetry. Charginos and neutralinos with masses between 100 GeV and 1100 GeV are excluded depending on the assumed decay branching fractions into a lepton (electron, muon, or $\tau$-lepton) plus a boson ($W$, $Z$, or Higgs).

Narrow resonances decaying into $WW$, $WZ$ or $ZZ$ boson pairs are searched for in 139 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the Large Hadron Collider from 2015 to 2018. The diboson system is reconstructed using pairs of high transverse momentum, large-radius jets. These jets are built from a combination of calorimeter- and tracker-inputs compatible with the hadronic decay of a boosted $W$ or $Z$ boson, using jet mass and substructure properties. The search is performed for diboson resonances with masses greater than 1.3 TeV. No significant deviations from the background expectations are observed. Exclusion limits at the 95% confidence level are set on the production cross-section times branching ratio into dibosons for resonances in a range of theories beyond the Standard Model, with the highest excluded mass of a new gauge boson at 3.8 TeV in the context of mass-degenerate resonances that couple predominantly to gauge bosons.

A search for supersymmetry in events with four or more charged leptons (electrons, muons and $\tau$-leptons) is presented. The analysis uses a data sample corresponding to $139\,\mbox{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 $\tau$-leptons are designed to target several supersymmetric models, while a general five-lepton signal region targets any new physics phenomena leading to a final state with five charged leptons. Data yields are consistent with Standard Model expectations and results are used to set upper limits on contributions from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of general gauge-mediated supersymmetry, excluding higgsino masses up to $540$ GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of $1.6$ TeV, $1.2$ TeV, and $2.5$ TeV are placed on wino, slepton and gluino masses, respectively.

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 non-resonant Higgs boson pair production, as predicted by the Standard Model, is presented, where one of the Higgs bosons decays via the $H\rightarrow bb$ channel and the other via one of the $H \rightarrow WW^*/ZZ^*/\tau\tau$ channels. The analysis selection requires events to have at least two $b$-tagged jets and exactly two leptons (electrons or muons) with opposite electric charge in the final state. Candidate events consistent with Higgs boson pair production are selected using a multi-class neural network discriminant. The analysis uses 139 fb$^{-1}$ of $pp$ collision data recorded at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. An observed (expected) upper limit of 1.2 ($0.9^{+0.4}_{-0.3}$) pb is set on the non-resonant Higgs boson pair production cross-section at 95% confidence level, which is equivalent to 40 ($29^{+14}_{-9}$) times the value predicted in the Standard Model.