Results from a search for supersymmetry in events with four or more leptons including electrons, muons and taus are presented. The analysis uses a data sample corresponding to 20.3 $fb^{-1}$ of proton--proton collisions delivered by the Large Hadron Collider at $\sqrt{s}$ = 8 TeV and recorded by the ATLAS detector. Signal regions are designed to target supersymmetric scenarios that can be either enriched in or depleted of events involving the production of a $Z$ boson. No significant deviations are observed in data from Standard Model predictions and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits at the 95% confidence level on the masses of relevant supersymmetric particles are obtained. In R-parity-violating simplified models with decays of the lightest supersymmetric particle to electrons and muons, limits of 1350 GeV and 750 GeV are placed on gluino and chargino masses, respectively. In R-parity-conserving simplified models with heavy neutralinos decaying to a massless lightest supersymmetric particle, heavy neutralino masses up to 620 GeV are excluded. Limits are also placed on other supersymmetric scenarios.

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.

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.

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of $20.1 \rm{fb}^{-1}$ of proton-proton collision data at $\sqrt{s}=8$ TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via $\tilde{t} \rightarrow t \tilde{\chi}_{1}^{0}$ or $\tilde{t}\rightarrow b\tilde{\chi}_{1}^{\pm} \rightarrow b W^{\left(\ast\right)} \tilde{\chi}_{1}^{0}$, where $\tilde{\chi}_{1}^{0}$ ($\tilde{\chi}_{1}^{\pm}$) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of $\tilde{t} \rightarrow t \tilde{\chi}_{1}^{0}$. For a branching fraction of 100%, top squark masses in the range 270-645 GeV are excluded for $\tilde{\chi}_{1}^{0}$ masses below 30 GeV. For a branching fraction of 50% to either $\tilde{t} \rightarrow t \tilde{\chi}_{1}^{0}$ or $\tilde{t}\rightarrow b\tilde{\chi}_{1}^{\pm}$, and assuming the $\tilde{\chi}_{1}^{\pm}$ mass to be twice the $\tilde{\chi}_{1}^{0}$ mass, top squark masses in the range 250-550 GeV are excluded for $\tilde{\chi}_{1}^{0}$ masses below 60 GeV.

The results of a search for top squark (stop) pair production in final states with one isolated lepton, jets, and missing transverse momentum are reported. The analysis is performed with proton--proton collision data at $\sqrt{s} = 8$ TeV collected with the ATLAS detector at the LHC in 2012 corresponding to an integrated luminosity of $20$ fb$^{-1}$. The lightest supersymmetric particle (LSP) is taken to be the lightest neutralino which only interacts weakly and is assumed to be stable. The stop decay modes considered are those to a top quark and the LSP as well as to a bottom quark and the lightest chargino, where the chargino decays to the LSP by emitting a $W$ boson. A wide range of scenarios with different mass splittings between the stop, the lightest neutralino and the lightest chargino are considered, including cases where the $W$ bosons or the top quarks are off-shell. Decay modes involving the heavier charginos and neutralinos are addressed using a set of phenomenological models of supersymmetry. No significant excess over the Standard Model prediction is observed. A stop with a mass between $210$ and $640$ GeV decaying directly to a top quark and a massless LSP is excluded at $95$ % confidence level, and in models where the mass of the lightest chargino is twice that of the LSP, stops are excluded at $95$ % confidence level up to a mass of $500$ GeV for an LSP mass in the range of $100$ to $150$ GeV. Stringent exclusion limits are also derived for all other stop decay modes considered, and model-independent upper limits are set on the visible cross-section for processes beyond the Standard Model.

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.

Many extensions of the Standard Model predict the existence of charged heavy long-lived particles, such as $R$-hadrons or charginos. These particles, if produced at the Large Hadron Collider, should be moving non-relativistically and are therefore identifiable through the measurement of an anomalously large specific energy loss in the ATLAS pixel detector. Measuring heavy long-lived particles through their track parameters in the vicinity of the interaction vertex provides sensitivity to metastable particles with lifetimes from 0.6 ns to 30 ns. A search for such particles with the ATLAS detector at the Large Hadron Collider is presented, based on a data sample corresponding to an integrated luminosity of 18.4 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 8 TeV. No significant deviation from the Standard Model background expectation is observed, and lifetime-dependent upper limits on $R$-hadrons and chargino production are set. Gluino $R$-hadrons with 10 ns lifetime and masses up to 1185 GeV are excluded at 95$\%$ confidence level, and so are charginos with 15 ns lifetime and masses up to 482 GeV.

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.

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.

Results are reported of a search for new phenomena, such as supersymmetric particle production, that could be observed in high-energy proton--proton collisions. Events with large numbers of jets, together with missing transverse momentum from unobserved particles, are selected. The data analysed were recorded by the ATLAS experiment during 2015 using the 13 TeV centre-of-mass proton--proton collisions at the Large Hadron Collider, and correspond to an integrated luminosity of 3.2 fb$^{-1}$. The search selected events with various jet multiplicities from $\ge 7$ to $\ge 10$ jets, and with various $b$-jet multiplicity requirements to enhance sensitivity. No excess above Standard Model expectations is observed. The results are interpreted within two supersymmetry models, where gluino masses up to 1400 GeV are excluded at 95% confidence level, significantly extending previous limits.

A search for long-lived particles decaying into hadrons is presented. The analysis uses 139 fb$^{-1}$ of $pp$ collision data collected at $\sqrt{s} = 13$ TeV by the ATLAS detector at the LHC using events that contain multiple energetic jets and a displaced vertex. The search employs dedicated reconstruction techniques that significantly increase the sensitivity to long-lived particles decaying 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 used to set limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model, and on scenarios with pair-production of supersymmetric particles with long-lived electroweakinos that decay via a small $R$-parity-violating coupling. The pair-production of electroweakinos with masses below 1.5 TeV is excluded for mean proper lifetimes in the range from 0.03 ns to 1 ns. When produced in the decay of $m(\tilde{g})=2.4$ TeV gluinos, electroweakinos with $m(\tilde\chi^0_1)=1.5$ TeV are excluded with lifetimes in the range of 0.02 ns to 4 ns.

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.

The results of a search for vector-like top quarks using events with exactly one lepton, at least four jets, and large missing transverse momentum are reported. The search is optimised for pair production of vector-like top quarks in the $Z(\rightarrow \! \! \nu \nu) \, t + X$ decay channel. LHC pp collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded by the ATLAS detector in 2015 and 2016 are used, corresponding to an integrated luminosity of 36.1 $\mathrm{fb}^{-1}$. No significant excess over the Standard Model expectation is seen and upper limits on the production cross-section of a vector-like $T$ quark pair as a function of the $T$ quark mass are derived. The observed (expected) 95% CL lower limits on the $T$ mass are 870 GeV (890 GeV) for the weak-isospin singlet model, 1.05 TeV (1.06 TeV) for the weak-isospin doublet model and 1.16 TeV (1.17 TeV) for the pure $Zt$ decay mode. Limits are also set on the mass as a function of the decay branching ratios, excluding large parts of the parameter space for masses below 1 TeV.

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.

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.

We present a search for magnetic monopoles and high-electric-charge objects using LHC Run 2 $\sqrt{s} =$13 TeV proton$-$proton collisions recorded by the ATLAS detector. A total integrated luminosity of 138 fb$^{-1}$ was collected by a specialized trigger. No highly ionizing particle candidate was observed. Considering the Drell-Yan and photon-fusion pair production mechanisms as benchmark models, cross-section upper limits are presented for spin-0 and spin-$\frac{1}{2}$ magnetic monopoles of magnetic charge $1g_\textrm{D}$ and $2g_\textrm{D}$ and for high-electric-charge objects of electric charge $20 \leq |z| \leq 100$, for masses between 200 GeV and 4000 GeV. The search improves by approximately a factor of three the previous cross-section limits on the Drell-Yan production of magnetic monopoles and high-electric charge objects. Also, the first ATLAS limits on the photon-fusion pair production mechanism of magnetic monopoles and high-electric-charge objects have been obtained.

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 is reported for excited $\tau$-leptons and leptoquarks in events with two hadronically decaying $\tau$-leptons and two or more jets. The search uses proton-proton (pp) collision data at $\sqrt{s} = 13$ TeV recorded by the ATLAS experiment during the Run 2 of the Large Hadron Collider in 2015-2018. The total integrated luminosity is 139 fb$^{-1}$. The excited $\tau$-lepton is assumed to be produced and to decay via a four-fermion contact interaction into an ordinary $\tau$-lepton and a quark-antiquark pair. The leptoquarks are assumed to be produced in pairs via the strong interaction, and each leptoquark is assumed to couple to a charm or lighter quark and a $\tau$-lepton. No excess over the background prediction is observed. Excited $\tau$-leptons with masses below 2.8 TeV are excluded at 95% CL in scenarios with the contact interaction scale $\Lambda$ set to 10 TeV. At the extreme limit of model validity where $\Lambda$ is set equal to the excited $\tau$-lepton mass, excited $\tau$-leptons with masses below 4.6 TeV are excluded. Leptoquarks with masses below 1.3 TeV are excluded at 95% CL if their branching ratio to a charm quark and a $\tau$-lepton equals 1. The analysis does not exploit flavour-tagging in the signal region.

A search for physics beyond the Standard Model inducing periodic signals in the dielectron and diphoton invariant mass spectra is presented using 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data collected by the ATLAS experiment at the LHC. Novel search techniques based on continuous wavelet transforms are used to infer the frequency of periodic signals from the invariant mass spectra and neural network classifiers are used to enhance the sensitivity to periodic resonances. In the absence of a signal, exclusion limits are placed at the 95% confidence level in the two-dimensional parameter space of the clockwork gravity model. Model-independent searches for deviations from the background-only hypothesis are also performed.

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).

A search for magnetic monopoles and high-electric-charge objects is presented using 34.4 fb$^{-1}$ of 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC during 2015 and 2016. The considered signature is based upon high ionization in the transition radiation tracker of the inner detector associated with a pencil-shape energy deposit in the electromagnetic calorimeter. The data were collected by a dedicated trigger based on the tracker high-threshold hit capability. The results are interpreted in models of Drell-Yan pair production of stable particles with two spin hypotheses (0 and 1/2) and masses ranging from 200 GeV to 4000 GeV. The search improves by approximately a factor of five the constraints on the direct production of magnetic monopoles carrying one or two Dirac magnetic charges and stable objects with electric charge in the range $20\le|z|\le60$ and extends the charge range to $60<|z|\le100$.

A search for a new massive charged gauge boson, $W'$, is performed with the ATLAS detector at the LHC. The dataset used in this analysis was collected from proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =13$ TeV, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The reconstructed $tb$ invariant mass is used to search for a $W'$ boson decaying into a top quark and a bottom quark. The result is interpreted in terms of a $W'$ boson with purely right-handed or left-handed chirality in a mass range of 0.5-6 TeV. Different values for the coupling of the $W'$ boson to the top and bottom quarks are considered, taking into account interference with single-top-quark production in the $s$-channel. No significant deviation from the background prediction is observed. The results are expressed as upper limits on the $W' \rightarrow tb$ production cross-section times branching ratio as a function of the $W'$-boson mass and in the plane of the coupling vs the $W'$-boson mass.

Higgsinos with masses near the electroweak scale can solve the hierarchy problem and provide a dark matter candidate, while detecting them at the LHC remains challenging if their mass-splitting is $\mathcal{O}$(1 GeV). This Letter presents a novel search for nearly mass-degenerate higgsinos in events with an energetic jet, missing transverse momentum, and a low-momentum track with a significant transverse impact parameter using 140 fb$^{-1}$ of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment. For the first time since LEP, a range of mass-splittings between the lightest charged and neutral higgsinos from 0.3 GeV to 0.9 GeV is excluded at 95% confidence level, with a maximum reach of approximately 170 GeV in the higgsino mass.

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.

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.