Double insertions of SMEFT operators in gluon fusion Higgs boson production

Deviations from the Standard Model (SM) can be parameterized in terms of the SM effective field theory (SMEFT), which is typically truncated at dimension-6. Including higher dimension operators -- as well as considering simultaneous insertions of multiple dimension-6 operators -- may be necessary in some processes, in order to correctly capture the properties of the underlying UV theory. As a step towards clarifying this in the Higgs boson production in gluon fusion process, we study double insertions of dimension-6 operators in the 1-loop virtual amplitude. We present needed Feynman rules up to $\mathcal{O}(1/\Lambda^4)$ and we numerically study the impact of various approximations to the $\mathcal{O}(1/\Lambda^4)$ expansion.