This Letter presents results from a combination of searches for Higgs boson pair production using 126$-$140 fb$^{-1}$ of proton-proton collision data at $\sqrt{s}=13$ TeV recorded with the ATLAS detector. At 95% confidence level (CL), the upper limit on the production rate is 2.9 times the standard model (SM) prediction, with an expected limit of 2.4 assuming no Higgs boson pair production. Constraints on the Higgs boson self-coupling modifier $\kappa_{\lambda}=\lambda_{HHH}/\lambda_{HHH}^\mathrm{SM}$, and the quartic $HHVV$ coupling modifier $\kappa_{2V}=g_{HHVV}/g_{HHVV}^\mathrm{SM}$, are derived individually, fixing the other parameter to its SM value. The observed 95% CL intervals are $-1.2 < \kappa_{\lambda} < 7.2$ and $0.6 < \kappa_{2V} < 1.5$, respectively, while the expected intervals are $-1.6 < \kappa_{\lambda} < 7.2$ and $0.4 < \kappa_{2V} < 1.6$ in the SM case. Constraints obtained for several interaction parameters within Higgs effective field theory are the strongest to date, offering insights into potential deviations from SM predictions.
A search is presented for non-resonant Higgs boson pair production, targeting the $bbZZ$, 4$V$ ($V$ = $W$ or $Z$), $VV\tau\tau$, 4$\tau$, $\gamma\gamma VV$ and $\gamma\gamma\tau\tau$ decay channels. Events are categorised based on the multiplicity of light charged leptons (electrons or muons), hadronically decaying tau leptons, and photons. The search is based on a data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 140 fb$^{-1}$. No evidence of the signal is found and the observed (expected) upper limit on the cross-section for non-resonant Higgs boson pair production is determined to be 17 (11) times the Standard Model predicted cross-section at 95% confidence level under the background-only hypothesis. The observed (expected) constraints on the $HHH$ coupling modifier, $\kappa_{\lambda}$, are determined to be $-6.2 < \kappa_{\lambda} < 11.6$ ($-4.5 < \kappa_{\lambda} < 9.6$) at 95% confidence level, assuming the Standard Model for the expected limits and that new physics would only affect $\kappa_{\lambda}$.
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.
Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.
A search for nonresonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}={13}$ TeV collected with the ATLAS detector at the Large Hadron Collider, and targets both the gluon-gluon fusion and vector-boson fusion production modes. No evidence of the signal is found and the observed (expected) upper limit on the cross-section for nonresonant Higgs boson pair production is determined to be 5.4 (8.1) times the Standard Model predicted cross-section at 95% confidence level. Constraints are placed on modifiers to the $HHH$ and $HHVV$ couplings. The observed (expected) $2\sigma$ constraints on the $HHH$ coupling modifier, $\kappa_\lambda$, are determined to be $[-3.5, 11.3]$ ($[-5.4, 11.4]$), while the corresponding constraints for the $HHVV$ coupling modifier, $\kappa_{2V}$, are $[-0.0, 2.1]$ ($[-0.1, 2.1]$). In addition, constraints on relevant coefficients are derived in the context of the Standard Model effective field theory and Higgs effective field theory, and upper limits on the $HH$ production cross-section are placed in seven Higgs effective field theory benchmark scenarios.