Higgs pair production in gluon fusion to higher orders in Higgs Effective Field Theory

Higgs pair production offers the opportunity to probe correlations among the couplings of one or two Higgs bosons to fermions and gauge bosons. In this context, it serves as a powerful test of the underlying Effective Field Theory (EFT) framework. In particular, while such couplings remain correlated in the Standard Model Effective Field Theory (SMEFT) at dimension six, they can become fully de-correlated in Higgs Effective Field Theory (HEFT) already at leading order in the EFT expansion. In this work, we study Higgs pair production via gluon fusion within the HEFT framework. We demonstrate that adopting a consistent power counting in combination with next-to-leading order (NLO) diagrams necessitates the inclusion of higher-dimensional operators beyond the leading ones. We analyze their phenomenological impact and re-assess critically the kinematic benchmark scenarios commonly used in experimental non-resonant di-Higgs searches in light of these additional contributions.