A search for resonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126-139 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector at the Large Hadron Collider. The analysis is divided into two channels, targeting Higgs boson decays which are reconstructed as pairs of small-radius jets or as individual large-radius jets. Spin-0 and spin-2 benchmark signal models are considered, both of which correspond to resonant $HH$ production via gluon$-$gluon fusion. The data are consistent with Standard Model predictions. Upper limits are set on the production cross-section times branching ratio to Higgs boson pairs of a new resonance in the mass range from 251 GeV to 5 TeV.
A search is performed for resonant and non-resonant Higgs boson pair production in the $\gamma\gamma b\bar{b}$ final state. The data set used corresponds to an integrated luminosity of 36.1 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. No significant excess relative to the Standard Model expectation is observed. The observed limit on the non-resonant Higgs boson pair cross-section is 0.73 pb at 95% confidence level. This observed limit is equivalent to 22 times the predicted Standard Model cross-section. The Higgs boson self-coupling ($\kappa_\lambda = \lambda_{HHH} / \lambda_{HHH}^{\rm SM}$) is constrained at 95% confidence level to $-8.2 < \kappa_\lambda < 13.2$. For resonant Higgs boson pair production through X $\rightarrow$ HH $\rightarrow$ $\gamma\gamma b\bar{b}$, the limit is presented, using the narrow-width approximation, as a function of $m_X$ in the range 260 GeV $< m_X <$ 1000 GeV. The observed limits range from 1.1 pb to 0.12 pb over this mass range.
Searches for non-resonant and resonant Higgs boson pair production are performed in the $\gamma\gamma WW^{*}$ channel with the final state of $\gamma\gamma\ell\nu jj$ using 36.1 fb$^{-1}$ of proton-proton collision data recorded at a centre-of-mass energy of $\sqrt{s} = 13$ TeV by the ATLAS detector at the Large Hadron Collider. No significant deviation from the Standard Model prediction is observed. A 95% confidence-level observed upper limit of 7.7 pb is set on the cross section for non-resonant production, while the expected limit is 5.4 pb. A search for a narrow-width resonance $X$ decaying to a pair of Standard Model Higgs bosons $HH$ is performed with the same set of data, and the observed upper limits on $\sigma(pp \rightarrow X) \times B(X \rightarrow HH)$ range between 40.0 pb and 6.1 pb for masses of the resonance between 260 GeV and 500 GeV, while the expected limits range between 17.6 pb and 4.4 pb. When deriving the limits above, the Standard Model branching ratios of the $H \rightarrow \gamma\gamma$ and $H \rightarrow WW^{*}$ are assumed.
A search for the non-resonant production of Higgs boson pairs in the $HH\rightarrow b\bar{b}\tau^+\tau^-$ channel is performed using 140 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of $13$ TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The analysis strategy is optimised to probe anomalous values of the Higgs boson self-coupling modifier $\kappa_\lambda$ and of the quartic $HHVV$ ($V = W,Z$) coupling modifier $\kappa_{2V}$. No significant excess above the expected background from Standard Model processes is observed. An observed (expected) upper limit $\mu_{HH}<5.9$$(3.3)$ is set at 95% confidence-level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The coupling modifiers are constrained to an observed (expected) 95% confidence interval of $-3.1 < \kappa_\lambda < 9.0$ ($-2.5 < \kappa_\lambda < 9.3$) and $-0.5 < \kappa_{2V} < 2.7$ ($-0.2 < \kappa_{2V} < 2.4$), assuming all other Higgs boson couplings are fixed to the Standard Model prediction. The results are also interpreted in the context of effective field theories via constraints on anomalous Higgs boson couplings and Higgs boson pair production cross-sections assuming different kinematic benchmark scenarios.