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

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

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.

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.

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 jets is presented. Data were collected with the CMS detector at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search examines the distinctive topology of displaced tracks and secondary vertices. The selected events are found to be consistent with standard model predictions. For a simplified model in which long-lived neutral particles are pair produced and decay to two jets, pair production cross sections larger than 0.2 fb are excluded at 95% confidence level for a long-lived particle mass larger than 1000 GeV and proper decay lengths between 3 and 130 mm. Several supersymmetry models with gauge-mediated supersymmetry breaking or $R$-parity violation, where pair-produced long-lived gluinos or top squarks decay to several final-state topologies containing displaced jets, are also tested. For these models, in the mass ranges above 200 GeV, gluino masses up to 2300-2400 GeV and top squark masses up to 1350-1600 GeV are excluded for proper decay lengths approximately between 10 and 100 mm. These are the most restrictive limits to date on these models.

A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

The production of single top quarks and top antiquarks via the $t$-channel exchange of a virtual $W$ boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the LHC using $140\,\mathrm{fb^{-1}}$ of ATLAS data. The total cross-sections are determined to be $\sigma(tq)=137^{+8}_{-8}\,\mathrm{pb}$ and $\sigma(\bar{t}q)=84^{+6}_{-5}\,\mathrm{pb}$ for top-quark and top-antiquark production, respectively. The combined cross-section is found to be $\sigma(tq+\bar{t}q)=221^{+13}_{-13}\,\mathrm{pb}$ and the cross-section ratio is $R_{t}=\sigma(tq)/\sigma(\bar{t}q)=1.636^{+0.036}_{-0.034}$. The predictions at next-to-next-to-leading-order in quantum chromodynamics are in good agreement with these measurements. The predicted value of $R_{t}$ using different sets of parton distribution functions is compared with the measured value, demonstrating the potential to further constrain the functions when using this result in global fits. The measured cross-sections are interpreted in an effective field theory approach, setting limits at the 95% confidence level on the strength of a four-quark operator and an operator coupling the third quark generation to the Higgs boson doublet: $-0.37 < C_{Qq}^{3,1}/\Lambda^2 < 0.06$ and $-0.87 < C_{\phi Q}^{3}/\Lambda^2 < 1.42$. The constraint $|V_{tb}|>0.95$ at the 95% confidence level is derived from the measured value of $\sigma(tq+\bar{t}q)$. In a more general approach, pairs of CKM matrix elements involving top quarks are simultaneously constrained, leading to confidence contours in the corresponding two-dimensional parameter spaces.

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

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb$^{-1}$, collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses.

A search for production of the supersymmetric partners of the top quark, top squarks, is presented. The search is based on proton-proton collision events containing multiple jets, no leptons, and large transverse momentum imbalance. The data were collected with the CMS detector at the CERN LHC at a center-of-mass energy of 13 TeV, and correspond to an integrated luminosity of 137 fb$^{-1}$. The targeted signal production scenarios are direct and gluino-mediated top squark production, including scenarios in which the top squark and neutralino masses are nearly degenerate. The search utilizes novel algorithms based on deep neural networks that identify hadronically decaying top quarks and W bosons, which are expected in many of the targeted signal models. No statistically significant excess of events is observed relative to the expectation from the standard model, and limits on the top squark production cross section are obtained in the context of simplified supersymmetric models for various production and decay modes. Exclusion limits as high as 1310 GeV are established at the 95% confidence level on the mass of the top squark for direct top squark production models, and as high as 2260 GeV on the mass of the gluino for gluino-mediated top squark production models. These results represent a significant improvement over the results of previous searches for supersymmetry by CMS in the same final state.

A search is reported for long-lived dark photons with masses between 0.1 GeV and 15 GeV, from exotic decays of Higgs bosons produced via vector-boson-fusion. Events that contain displaced collimated Standard Model fermions reconstructed in the calorimeter or muon spectrometer are probed. This search uses the full LHC Run 2 (2015-2018) data sample collected in proton-proton collisions at $\sqrt{s}=13$ TeV, corresponding to an integrated luminosity of 139 $fb^{-1}$. Dominant backgrounds from Standard Model processes and non-collision sources are estimated by using data-driven techniques. The observed event yields in the signal regions are consistent with the expected background. Upper limits on the Higgs boson to dark photon branching fraction are reported as a function of the dark-photon mean proper decay length or of the dark-photon mass and the coupling between the Standard Model and the potential dark sector. This search is combined with previous ATLAS searches obtained in the gluon-gluon fusion and \textit{WH} production modes. A branching fraction above 10% is excluded at 95% CL for a 125 GeV Higgs boson decaying into two dark photons for dark-photon mean proper decay lengths between 173 and 1296 mm and mass of 10 GeV.

A search for a new heavy boson produced via gluon-fusion in the four-lepton channel with missing transverse momentum or jets is performed. The search uses proton-proton collision data equivalent to an integrated luminosity of 139 fb$^{-1}$ at a centre-of-mass energy of 13 TeV collected by the ATLAS detector between 2015 and 2018 at the Large Hadron Collider. This study explores the decays of heavy bosons: $R\rightarrow SH$ and $A\rightarrow ZH$, where $R$ is a CP-even boson, $A$ is a CP-odd boson, $H$ is a CP-even boson, and $S$ is considered to decay into invisible particles that are candidates for dark matter. In these processes, $S\rightarrow \textrm{invisible}$ and $H\rightarrow ZZ$. The $Z$ boson associated with the heavy scalar boson $H$ decays into all decay channels of the $Z$ boson. The mass range under consideration is 390-1300 (320-1300) GeV for the $R$ ($A$) boson and 220-1000 GeV for the $H$ boson. No significant deviation from the Standard Model backgrounds is observed. The results are interpreted as upper limits at a 95% confidence level on the cross-section times the branching ratio of the heavy resonances.

Results are reported from a search for supersymmetric particles in the final state with multiple jets and large missing transverse momentum. The search uses a sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected with the CMS detector in 2016-2018, corresponding to an integrated luminosity of 137 fb$^{-1}$, representing essentially the full LHC Run 2 data sample. The analysis is performed in a four-dimensional search region defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No significant excess in the event yield is observed relative to the expected background contributions from standard model processes. Limits on the pair production of gluinos and squarks are obtained in the framework of simplified models for supersymmetric particle production and decay processes. Assuming the lightest supersymmetric particle to be a neutralino, lower limits on the gluino mass as large as 2000 to 2310 GeV are obtained at 95% confidence level, while lower limits on the squark mass as large as 1190 to 1630 GeV are obtained, depending on the production scenario.

A search for long-lived particles (LLPs) produced in decays of standard model (SM) Higgs bosons is presented. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, recorded at the LHC in 2016-2018. A novel technique is employed to reconstruct decays of LLPs in the endcap muon detectors. The search is sensitive to a broad range of LLP decay modes and to masses as low as a few GeV. No excess of events above the SM background is observed. The most stringent limits to date on the branching fraction of the Higgs boson to LLPs subsequently decaying to quarks and $\tau^+\tau^-$ are found for proper decay lengths greater than 6, 20, and 40 m, for LLP masses of 7, 15, and 40 GeV, respectively.

A search for long-lived particles (LLPs) produced in association with a Z boson is presented. The study is performed using data from proton-proton collisions with a center-of-mass energy of 13 TeV recorded by the CMS experiment during 2016-2018, corresponding to an integrated luminosity of 117 fb$^{-1}$. The LLPs are assumed to decay to a pair of standard model quarks that are identified as displaced jets within the CMS tracker system. Triggers and selections based on Z boson decays to electron or muon pairs improve the sensitivity to light LLPs (down to 15 GeV). This search provides sensitivity to beyond the standard model scenarios which predict LLPs produced in association with a Z boson. In particular, the results are interpreted in the context of exotic decays of the Higgs boson to a pair of scalar LLPs (H $\to$ SS). The Higgs boson decay branching fraction is constrained to values less than 6% for proper decay lengths of 10-100 mm and for LLP masses between 40 and 55 GeV. In the case of low-mass ($\approx$ 15 GeV) scalar particles that subsequently decay to a pair of b quarks, the search is sensitive to branching fractions $\mathcal{B}$(H $\to$ SS) $\lt$ 20% for proper decay lengths of 10-50 mm. The use of associated production with a Z boson increases the sensitivity to low-mass LLPs of this analysis with respect to gluon fusion searches. In the case of 15 GeV scalar LLPs, the improvement corresponds to a factor of 2 at a proper decay length of 30 mm.

The first collider search for dark matter arising from a strongly coupled hidden sector is presented and uses a data sample corresponding to 138 fb$^{-1}$, collected with the CMS detector at the CERN LHC, at $\sqrt{s} =$ 13 TeV. The hidden sector is hypothesized to couple to the standard model (SM) via a heavy leptophobic Z' mediator produced as a resonance in proton-proton collisions. The mediator decay results in two "semivisible" jets, containing both visible matter and invisible dark matter. The final state therefore includes moderate missing energy aligned with one of the jets, a signature ignored by most dark matter searches. No structure in the dijet transverse mass spectra compatible with the signal is observed. Assuming the Z' has a universal coupling of 0.25 to the SM quarks, an inclusive search, relevant to any model that exhibits this kinematic behavior, excludes mediator masses of 1.5-4.0 TeV at 95% confidence level, depending on the other signal model parameters. To enhance the sensitivity of the search for this particular class of hidden sector models, a boosted decision tree (BDT) is trained using jet substructure variables to distinguish between semivisible jets and SM jets from background processes. When the BDT is employed to identify each jet in the dijet system as semivisible, the mediator mass exclusion increases to 5.1 TeV, for wider ranges of the other signal model parameters. These limits exclude a wide range of strongly coupled hidden sector models for the first time.

A search for heavy neutral leptons (HNLs), the right-handed Dirac or Majorana neutrinos, is performed in final states with three charged leptons (electrons or muons) using proton-proton collision data collected by the CMS experiment at $\sqrt{s} =$ 13 TeV at the CERN LHC. The data correspond to an integrated luminosity of 138 fb$^{-1}$. The HNLs could be produced through mixing with standard model neutrinos $\nu$. For small values of the HNL mass ($\lt$ 20 GeV) and the square of the HNL-$\nu$ mixing parameter (10$^{-7}$-10$^{-2}$), the decay length of these particles can be large enough so that the secondary vertex of the HNL decay can be resolved with the CMS silicon tracker. The selected final state consists of one lepton emerging from the primary proton-proton collision vertex, and two leptons forming a displaced, secondary vertex. No significant deviations from the standard model expectations are observed, and constraints are obtained on the HNL mass and coupling strength parameters, excluding previously unexplored regions of parameter space in the mass range 1-20 GeV and squared mixing parameter values as low as 10$^{-7}$.

A search for exclusive two-photon production via photon exchange in proton-proton collisions, pp $\to$ p$\gamma\gamma$p with intact protons, is presented. The data correspond to an integrated luminosity of 9.4 fb$^{-1}$ collected in 2016 using the CMS and TOTEM detectors at a center-of-mass energy of 13 TeV at the LHC. Events are selected with a diphoton invariant mass above 350 GeV and with both protons intact in the final state, to reduce backgrounds from strong interactions. The events of interest are those where the invariant mass and rapidity calculated from the momentum losses of the forward-moving protons matches the mass and rapidity of the central, two-photon system. No events are found that satisfy this condition. Interpreting this result in an effective dimension-8 extension of the standard model, the first limits are set on the two anomalous four-photon coupling parameters. If the other parameter is constrained to its standard model value, the limits at 95% CL are $\lvert\zeta_1\rvert$ $\lt$ 2.9 $\times$ 10$^{-13}$ GeV$^{-4}$ and $\lvert\zeta_2\rvert$ $\lt$ 6.0 $\times$ 10$^{-13}$ GeV$^{-4}$.

A search for physics beyond the standard model (SM) in final states with an electron or muon and missing transverse momentum is presented. The analysis uses data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected with the CMS detector at the LHC in 2016–2018 and corresponding to an integrated luminosity of 138 fb−1. No significant deviation from the SM prediction is observed. Model-independent limits are set on the production cross section of W’ bosons decaying into lepton-plus-neutrino final states. Within the framework of the sequential standard model, with the combined results from the electron and muon decay channels a W’ boson with mass less than 5.7 TeV is excluded at 95% confidence level. Results on a SM precision test, the determination of the oblique electroweak W parameter, are presented using LHC data for the first time. These results together with those from the direct W’ resonance search are used to extend existing constraints on composite Higgs scenarios. This is the first experimental exclusion on compositeness parameters using results from LHC data other than Higgs boson measurements.

A search is reported for pairs of light Higgs bosons (H$_1$) produced in supersymmetric cascade decays in final states with small missing transverse momentum. A data set of LHC pp collisions collected with the CMS detector at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 138 fb$^{-1}$ is used. The search targets events where both H$_1$ bosons decay into $\mathrm{b\bar{b}}$ pairs that are reconstructed as large-radius jets using substructure techniques. No evidence is found for an excess of events beyond the background expectations of the standard model (SM). Results from the search are interpreted in the next-to-minimal supersymmetric extension of the SM, where a "singlino" of small mass leads to squark and gluino cascade decays that can predominantly end in a highly Lorentz-boosted singlet-like H$_1$ and a singlino-like neutralino of small transverse momentum. Upper limits are set on the product of the squark or gluino pair production cross section and the square of the $\mathrm{b\bar{b}}$ branching fraction of the H$_1$ in a benchmark model containing almost mass-degenerate gluinos and light-flavour squarks. Under the assumption of an SM-like H$_1$$\to$$\mathrm{b\bar{b}}$ branching fraction, H$_1$ bosons with masses in the range 40-120 GeV arising from the decays of squarks or gluinos with a mass of 1200 to 2500 GeV are excluded at 95% confidence level.

A search is presented for heavy bosons decaying to Z($\nu\bar{\nu}$)V(qq'), where V can be a W or a Z boson. A sample of proton-proton collision data at $\sqrt{s} =$ 13 TeV was collected by the CMS experiment during 2016-2018. The data correspond to an integrated luminosity of 137 fb$^{-1}$. The event categorization is based on the presence of high-momentum jets in the forward region to identify production through weak vector boson fusion. Additional categorization uses jet substructure techniques and the presence of large missing transverse momentum to identify W and Z bosons decaying to quarks and neutrinos, respectively. The dominant standard model backgrounds are estimated using data taken from control regions. The results are interpreted in terms of radion, W' boson, and graviton models, under the assumption that these bosons are produced via gluon-gluon fusion, Drell-Yan, or weak vector boson fusion processes. No evidence is found for physics beyond the standard model. Upper limits are set at 95% confidence level on various types of hypothetical new bosons. Observed (expected) exclusion limits on the masses of these bosons range from 1.2 to 4.0 (1.1 to 3.7) TeV.

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.