Study of Higgs boson pair production in the $HH\to b\bar{b}\gamma \gamma$ final state with 308 fb${}^{-1}$ of data collected at $\sqrt{s}=$ 13 TeV and 13.6 TeV by the ATLAS experiment

A search for Higgs boson pair production in the $b\bar{b}\gamma \gamma$ final state is performed. The proton-proton collision dataset corresponds to an integrated luminosity of 308 fb${}^{-1}$, consisting of two samples, 140 fb${}^{-1}$ at a centre-of-mass energy of $\sqrt{s}=13$ TeV and 168 fb${}^{-1}$ at $\sqrt{s}=13.6$ TeV, recorded between 2015 and 2024 by the ATLAS detector at the CERN Large Hadron Collider. In addition to a larger dataset, this analysis improves upon the previous search in the same final state through several methodological and technical developments. The Higgs boson pair production cross section divided by the Standard Model prediction is found to be ${\mu }_{HH}=0.9_{-1.1}^{+1.4}$ (${\mu }_{HH}=1_{-1.0}^{+1.3}$ expected), which translates into a 95% confidence-level upper limit of ${\mu }_{HH}<3.7$. At the same confidence level the Higgs self-coupling modifier is constrained to be in the range $-1.6<{\kappa }_{\lambda }<6.6$ ($-1.8<{\kappa }_{\lambda }<6.9$ expected).