Measurements of the suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

Jet shapes have been measured in inclusive jet production in proton-proton collisions at sqrt(s) = 7 TeV using 3 pb^{-1} of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed using the anti-kt algorithm with transverse momentum 30 GeV < pT < 600 GeV and rapidity in the region |y| < 2.8. The data are corrected for detector effects and compared to several leading-order QCD matrix elements plus parton shower Monte Carlo predictions, including different sets of parameters tuned to model fragmentation processes and underlying event contributions in the final state. The measured jets become narrower with increasing jet transverse momentum and the jet shapes present a moderate jet rapidity dependence. Within QCD, the data test a variety of perturbative and non-perturbative effects. In particular, the data show sensitivity to the details of the parton shower, fragmentation, and underlying event models in the Monte Carlo generators. For an appropriate choice of the parameters used in these models, the data are well described.

A measurement of the cross section for the inclusive production of isolated prompt photons in pp collisions at a centre-of-mass energy sqrt(s) = 7TeV is presented. The measurement covers the pseudorapidity ranges |eta|<1.37 and 1.52<|eta|<1.81 in the transverse energy range 15 < E_T <100 GeV. The results are based on an integrated luminosity of 880 nb-1, collected with the ATLAS detector at the Large Hadron Collider. Photon candidates are identified by combining information from the calorimeters and from the inner tracker. Residual background in the selected sample is estimated from data based on the observed distribution of the transverse isolation energy in a narrow cone around the photon candidate. The results are compared to predictions from next-to-leading order perturbative QCD calculations.

Dijet angular distributions from the first LHC pp collisions at center-of-mass energy sqrt(s) = 7 TeV have been measured with the ATLAS detector. The dataset used for this analysis represents an integrated luminosity of 3.1 pb-1. Dijet $\chi$ distributions and centrality ratios have been measured up to dijet masses of 2.8 TeV, and found to be in good agreement with Standard Model predictions. Analysis of the $\chi$ distributions excludes quark contact interactions with a compositeness scale $\Lambda$ below 3.4 TeV, at 95% confidence level, significantly exceeding previous limits.

A search for new heavy particles manifested as resonances in two-jet final states is presented. The data were produced in 7 TeV proton-proton collisions by the Large Hadron Collider (LHC) and correspond to an integrated luminosity of 315 nb^-1 collected by the ATLAS detector. No resonances were observed. Upper limits were set on the product of cross section and signal acceptance for excited-quark (q*) production as a function of q* mass. These exclude at the 95% CL the q* mass interval 0.30 < mq* < 1.26 TeV, extending the reach of previous experiments.

A search for flavor-changing neutral-current couplings between a top quark, an up or charm quark and a $Z$ boson is presented, using proton-proton collision data at $\sqrt{s} = 13$ TeV collected by the ATLAS detector at the Large Hadron Collider. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets both single-top-quark events produced as $gq\rightarrow tZ$ (with $q = u, c$) and top-quark-pair events, with one top quark decaying through the $t \rightarrow Zq$ channel. The analysis considers events with three leptons (electrons or muons), a $b$-tagged jet, possible additional jets, and missing transverse momentum. The data are found to be consistent with the background-only hypothesis and 95% confidence-level limits on the $t \rightarrow Zq$ branching ratios are set, assuming only tensor operators of the Standard Model effective field theory framework contribute to the $tZq$ vertices. These are $6.2 \times 10^{-5}$ ($13\times 10^{-5}$) for $t\rightarrow Zu$ ($t\rightarrow Zc$) for a left-handed $tZq$ coupling, and $6.6 \times 10^{-5}$ ($12\times 10^{-5}$) in the case of a right-handed coupling. These results are interpreted as 95% CL upper limits on the strength of corresponding couplings, yielding limits for $|C_{uW}^{(13)*}|$ and $|C_{uB}^{(13)*}|$ ($|C_{uW}^{(31)}|$ and $|C_{uB}^{(31)}|$) of 0.15 (0.16), and limits for $|C_{uW}^{(23)*}|$ and $|C_{uB}^{(23)*}|$ ($|C_{uW}^{(32)}|$ and $|C_{uB}^{(32)}|$) of 0.22 (0.21), assuming a new-physics energy scale $\Lambda_\text{NP}$ of 1 TeV.

Jet cross sections have been measured for the first time in proton-proton collisions at a centre-of-mass energy of 7 TeV using the ATLAS detector. The measurement uses an integrated luminosity of 17 nb-1 recorded at the Large Hadron Collider. The anti-kt algorithm is used to identify jets, with two jet resolution parameters, R = 0.4 and 0.6. The dominant uncertainty comes from the jet energy scale, which is determined to within 7% for central jets above 60 GeV transverse momentum. Inclusive single-jet differential cross sections are presented as functions of jet transverse momentum and rapidity. Dijet cross sections are presented as functions of dijet mass and the angular variable $\chi$. The results are compared to expectations based on next-to-leading-order QCD, which agree with the data, providing a validation of the theory in a new kinematic regime.

Jet quenching is the process of color-charged partons losing energy via interactions with quark-gluon plasma droplets created in heavy-ion collisions. The collective expansion of such droplets is well described by viscous hydrodynamics. Similar evidence of collectivity is consistently observed in smaller collision systems, including $pp$ and $p$+Pb collisions. In contrast, while jet quenching is observed in Pb+Pb collisions, no evidence has been found in these small systems to date, raising fundamental questions about the nature of the system created in these collisions. The ATLAS experiment at the Large Hadron Collider has measured the yield of charged hadrons correlated with reconstructed jets in 0.36 nb$^{-1}$ of $p$+Pb and 3.6 pb$^{-1}$ of $pp$ collisions at 5.02 TeV. The yields of charged hadrons with $p_\mathrm{T}^\mathrm{ch} >0.5$ GeV near and opposite in azimuth to jets with $p_\mathrm{T}^\mathrm{jet} > 30$ or $60$ GeV, and the ratios of these yields between $p$+Pb and $pp$ collisions, $I_{p\mathrm{Pb}}$, are reported. The collision centrality of $p$+Pb events is categorized by the energy deposited by forward neutrons from the struck nucleus. The $I_{p\mathrm{Pb}}$ values are consistent with unity within a few percent for hadrons with $p_\mathrm{T}^\mathrm{ch} >4$ GeV at all centralities. These data provide new, strong constraints which preclude almost any parton energy loss in central $p$+Pb collisions.

First measurements of the W -> lnu and Z/gamma* -> ll (l = e, mu) production cross sections in proton-proton collisions at sqrt(s) = 7 TeV are presented using data recorded by the ATLAS experiment at the LHC. The results are based on 2250 W -> lnu and 179 Z/gamma* -> ll candidate events selected from a data set corresponding to an integrated luminosity of approximately 320 nb-1. The measured total W and Z/gamma*-boson production cross sections times the respective leptonic branching ratios for the combined electron and muon channels are $\stotW$ * BR(W -> lnu) = 9.96 +- 0.23(stat) +- 0.50(syst) +- 1.10(lumi) nb and $\stotZg$ * BR(Z/gamma* -> ll) = 0.82 +- 0.06(stat) +- 0.05(syst) +- 0.09(lumi) nb (within the invariant mass window 66 < m_ll < 116 GeV). The W/Z cross-section ratio is measured to be 11.7 +- 0.9(stat) +- 0.4(syst). In addition, measurements of the W+ and W- production cross sections and of the lepton charge asymmetry are reported. Theoretical predictions based on NNLO QCD calculations are found to agree with the measurements.

This Letter reports the observation of $WZ\gamma$ production and a measurement of its cross-section using 140.1 $\pm$ 1.2 fb$^{-1}$ of proton-proton collision data recorded at a center-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The $WZ\gamma$ production cross-section, with both the $W$ and $Z$ bosons decaying leptonically, $pp \rightarrow WZ\gamma \rightarrow {\ell'}^{\pm}\nu\ell^{+}\ell^{-}\gamma$ ($\ell^{(')} = e, \mu$), is measured in a fiducial phase-space region defined such that the leptons and the photon have high transverse momentum and the photon is isolated. The cross-section is found to be 2.01 $\pm$ 0.30 (stat.) $\pm$ 0.16 (syst) fb. The corresponding Standard Model predicted cross-section calculated at next-to-leading order in perturbative quantum chromodynamics and at leading order in the electroweak coupling constant is 1.50 $\pm$ 0.06 fb. The observed significance of the $WZ\gamma$ signal is 6.3$\sigma$, compared with an expected significance of 5.0$\sigma$.

The exclusive production of pion pairs in the process $pp\to pp\pi^+\pi^-$ has been measured at $\sqrt{s}$ = 7 TeV with the ATLAS detector at the LHC, using 80 $\mu$b$^{-1}$ of low-luminosity data. The pion pairs were detected in the ATLAS central detector while outgoing protons were measured in the forward ATLAS ALFA detector system. This represents the first use of proton tagging to measure an exclusive hadronic final state at the LHC. A cross-section measurement is performed in two kinematic regions defined by the proton momenta, the pion rapidities and transverse momenta, and the pion-pion invariant mass. Cross section values of $4.8 \pm 1.0 \text{(stat.)} + {}^{+0.3}_{-0.2} \text{(syst.)}\mu$b and $9 \pm 6 \text{(stat.)} + {}^{+2}_{-2}\text{(syst.)}\mu$b are obtained in the two regions; they are compared with theoretical models and provide a demonstration of the feasibility of measurements of this type.

Studies of the correlations of the two highest transverse momentum (leading) jets in individual Pb+Pb collision events can provide information about the mechanism of jet quenching by the hot and dense matter created in such collisions. In Pb+Pb and pp collisions at $\sqrt{s_{_\text{NN}}}$ = 5.02 TeV, measurements of the leading dijet transverse momentum ($p_{\mathrm{T}}$) correlations are presented. Additionally, measurements in Pb+Pb collisions of the dijet pair nuclear modification factors projected along leading and subleading jet $p_{\mathrm{T}}$ are made. The measurements are performed using the ATLAS detector at the LHC with 260 pb$^{-1}$ of pp data collected in 2017 and 2.2 nb$^{-1}$ of Pb+Pb data collected in 2015 and 2018. An unfolding procedure is applied to the two-dimensional leading and subleading jet $p_{\mathrm{T}}$ distributions to account for experimental effects in the measurement of both jets. Results are provided for dijets with leading jet $p_{\mathrm{T}}$ greater than 100 GeV. Measurements of the dijet-yield-normalized $x_{\mathrm{J}}$ distributions in Pb+Pb collisions show an increased fraction of imbalanced jets compared to pp collisions; these measurements are in agreement with previous measurements of the same quantity at 2.76 TeV in the overlapping kinematic range. Measurements of the absolutely-normalized dijet rate in Pb+Pb and pp collisions are also presented, and show that balanced dijets are significantly more suppressed than imbalanced dijets in Pb+Pb collisions. It is observed in the measurements of the pair nuclear modification factors that the subleading jets are significantly suppressed relative to leading jets with $p_{\mathrm{T}}$ between 100 and 316 GeV for all centralities in Pb+Pb collisions.

The azimuthal variation of jet yields in heavy-ion collisions provides information about the path-length dependence of the energy loss experienced by partons passing through the hot, dense nuclear matter known as the quark-gluon plasma. This paper presents the azimuthal anisotropy coefficients $v_2$, $v_3$, and $v_4$ measured for jets in Pb+Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV using the ATLAS detector at the LHC. The measurement uses data collected in 2015 and 2018, corresponding to an integrated luminosity of 2.2 nb$^{-1}$. The $v_n$ values are measured as a function of the transverse momentum of the jets between 71 GeV and 398 GeV and the event centrality. A nonzero value of $v_2$ is observed in all but the most central collisions. The value of $v_2$ is largest for jets with lower transverse momentum, with values up to 0.05 in mid-central collisions. A smaller, nonzero value of $v_3$ of approximately 0.01 is measured with no significant dependence on jet $p_T$ or centrality, suggesting that fluctuations in the initial state play a small but distinct role in jet energy loss. No significant deviation of $v_4$ from zero is observed in the measured kinematic region.

A measurement of the top-quark mass ($m_t$) in the $t\bar{t}\rightarrow~\textrm{lepton}+\textrm{jets}$ channel is presented, with an experimental technique which exploits semileptonic decays of $b$-hadrons produced in the top-quark decay chain. The distribution of the invariant mass $m_{\ell\mu}$ of the lepton, $\ell$ (with $\ell=e,\mu$), from the $W$-boson decay and the muon, $\mu$, originating from the $b$-hadron decay is reconstructed, and a binned-template profile likelihood fit is performed to extract $m_t$. The measurement is based on data corresponding to an integrated luminosity of 36.1 fb$^{-1}$ of $\sqrt{s} = 13~\textrm{TeV}$$pp$ collisions provided by the Large Hadron Collider and recorded by the ATLAS detector. The measured value of the top-quark mass is $m_{t} = 174.41\pm0.39~(\textrm{stat.})\pm0.66~(\textrm{syst.})\pm0.25~(\textrm{recoil})~\textrm{GeV}$, where the third uncertainty arises from changing the PYTHIA8 parton shower gluon-recoil scheme, used in top-quark decays, to a recently developed setup.

A search for new phenomena in final states with hadronically decaying tau leptons, $b$-jets, and missing transverse momentum is presented. The analyzed dataset comprises $pp$~collision data at a center-of-mass energy of $\sqrt s = 13$ TeV with an integrated luminosity of 139/fb, delivered by the Large Hadron Collider and recorded with the ATLAS detector from 2015 to 2018. The observed data are compatible with the expected Standard Model background. The results are interpreted in simplified models for two different scenarios. The first model is based on supersymmetry and considers pair production of top squarks, each of which decays into a $b$-quark, a neutrino and a tau slepton. Each tau slepton in turn decays into a tau lepton and a nearly massless gravitino. Within this model, top-squark masses up to 1.4 TeV can be excluded at the 95% confidence level over a wide range of tau-slepton masses. The second model considers pair production of leptoquarks with decays into third-generation leptons and quarks. Depending on the branching fraction into charged leptons, leptoquarks with masses up to around 1.25 TeV can be excluded at the 95% confidence level for the case of scalar leptoquarks and up to 1.8 TeV (1.5 TeV) for vector leptoquarks in a Yang--Mills (minimal-coupling) scenario. In addition, model-independent upper limits are set on the cross section of processes beyond the Standard Model.

This paper presents a search for dark matter in the context of a two-Higgs-doublet model together with an additional pseudoscalar mediator, $a$, which decays into the dark-matter particles. Processes where the pseudoscalar mediator is produced in association with a single top quark in the 2HDM+$a$ model are explored for the first time at the LHC. Several final states which include either one or two charged leptons (electrons or muons) and a significant amount of missing transverse momentum are considered. The analysis is based on proton-proton collision data collected with the ATLAS experiment at $\sqrt{s} = 13$ TeV during LHC Run2 (2015-2018), corresponding to an integrated luminosity of 139 fb$^{-1}$. No significant excess above the Standard Model predictions is found. The results are expressed as 95% confidence-level limits on the parameters of the signal models considered.

This paper reports a search for Higgs boson pair ($hh$) production in association with a vector boson ($W$ or $Z$) using 139 $fb^{-1}$ of proton-proton collision data at $\sqrt{s}=$ 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed in final states in which the vector boson decays leptonically ($W\to\ell\nu, Z\to\ell\ell,\nu\nu$ with $\ell=e, \mu$) and the Higgs bosons each decay into a pair of $b$-quarks. It targets $Vhh$ signals from both non-resonant $hh$ production, present in the Standard Model (SM), and resonant $hh$ production, as predicted in some SM extensions. A 95% confidence-level upper limit of 183 (87) times the SM cross-section is observed (expected) for non-resonant $Vhh$ production when assuming the kinematics are as expected in the SM. Constraints are also placed on Higgs boson coupling modifiers. For the resonant search, upper limits on the production cross-sections are derived for two specific models: one is the production of a vector boson along with a neutral heavy scalar resonance $H$, in the mass range 260-1000 GeV, that decays into $hh$, and the other is the production of a heavier neutral pseudoscalar resonance $A$ that decays into a $Z$ boson and $H$ boson, where the $A$ boson mass is 360-800 GeV and the $H$ boson mass is 260-400 GeV. Constraints are also derived in the parameter space of two-Higgs-doublet models.

A search for a $WZ$ resonance, in the fully leptonic final state (electrons and muons), is performed using 139 fb$^{-1}$ of data collected at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The results are interpreted in terms of a singly charged Higgs boson of the Georgi$-$Machacek model, produced by $WZ$ fusion, and of a Heavy Vector Triplet, with the resonance produced by $WZ$ fusion or the Drell$-$Yan process. No significant excess over the Standard Model predictions is observed and limits are set on the production cross-section times branching ratio as a function of the resonance mass for these processes.

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 a long-lived, heavy neutral lepton ($\mathcal{N}$) in 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data collected by the ATLAS detector at the Large Hadron Collider is reported. The $\mathcal{N}$ is produced via $W \rightarrow \mathcal{N} \mu$ or $W \rightarrow \mathcal{N} e$ and decays into two charged leptons and a neutrino, forming a displaced vertex. The $\mathcal{N}$ mass is used to discriminate between signal and background. No signal is observed, and limits are set on the squared mixing parameters of the $\mathcal{N}$ with the left-handed neutrino states for the $\mathcal{N}$ mass range $3$ GeV $< m_{\mathcal{N}} < 15$ GeV. For the first time, limits are given for both single-flavor and multiflavor mixing scenarios motivated by neutrino flavor oscillation results for both the normal and inverted neutrino-mass hierarchies.

Several extensions of the Standard Model predict the production of dark matter particles at the LHC. A search for dark matter particles produced in association with a dark Higgs boson decaying into $W^{+}W^{-}$ in the $\ell^\pm\nu q \bar q'$ final states with $\ell=e,\mu$ is presented. This analysis uses 139 fb$^{-1}$ of $pp$ collisions recorded by the ATLAS detector at a centre-of-mass energy of 13 TeV. The $W^\pm \to q\bar q'$ decays are reconstructed from pairs of calorimeter-measured jets or from track-assisted reclustered jets, a technique aimed at resolving the dense topology from a pair of boosted quarks using jets in the calorimeter and tracking information. The observed data are found to agree with Standard Model predictions. Scenarios with dark Higgs boson masses ranging between 140 and 390 GeV are excluded.

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.

A search for Majorana neutrinos in same-sign $WW$ scattering events is presented. The analysis uses $\sqrt{s}= 13$ TeV proton-proton collision data with an integrated luminosity of 140 fb$^{-1}$ recorded during 2015-2018 by the ATLAS detector at the Large Hadron Collider. The analysis targets final states including exactly two same-sign muons and at least two hadronic jets well separated in rapidity. The modelling of the main backgrounds, from Standard Model same-sign $WW$ scattering and $WZ$ production, is constrained with data in dedicated signal-depleted control regions. The distribution of the transverse momentum of the second-hardest muon is used to search for signals originating from a heavy Majorana neutrino with a mass between 50 GeV and 20 TeV. No significant excess is observed over the background expectation. The results are interpreted in a benchmark scenario of the Phenomenological Type-I Seesaw model. In addition, the sensitivity to the Weinberg operator is investigated. Upper limits at the 95% confidence level are placed on the squared muon-neutrino-heavy-neutrino mass-mixing matrix element $\vert V_{\mu N} \vert^{2}$ as a function of the heavy Majorana neutrino's mass $m_N$, and on the effective $\mu\mu$ Majorana neutrino mass $|m_{\mu\mu}|$.

This paper presents a search for dark matter, $\chi$, using events with a single top quark and an energetic $W$ boson. The analysis is based on proton-proton collision data collected with the ATLAS experiment at $\sqrt{s}=$ 13 TeV during LHC Run 2 (2015-2018), corresponding to an integrated luminosity of 139 fb$^{-1}$. The search considers final states with zero or one charged lepton (electron or muon), at least one $b$-jet and large missing transverse momentum. In addition, a result from a previous search considering two-charged-lepton final states is included in the interpretation of the results. The data are found to be in good agreement with the Standard Model predictions and the results are interpreted in terms of 95% confidence-level exclusion limits in the context of a class of dark matter models involving an extended two-Higgs-doublet sector together with a pseudoscalar mediator particle. The search is particularly sensitive to on-shell production of the charged Higgs boson state, $H^{\pm}$, arising from the two-Higgs-doublet mixing, and its semi-invisible decays via the mediator particle, $a$: $H^{\pm} \rightarrow W^\pm a (\rightarrow \chi\chi)$. Signal models with $H^{\pm}$ masses up to 1.5 TeV and $a$ masses up to 350 GeV are excluded assuming a tan$\beta$ value of 1. For masses of $a$ of 150 (250) GeV, tan$\beta$ values up to 2 are excluded for $H^{\pm}$ masses between 200 (400) GeV and 1.5 TeV. Signals with tan$\beta$ values between 20 and 30 are excluded for $H^{\pm}$ masses between 500 and 800 GeV.

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.