An analytic result for the 0 → ggHHH amplitude

A bstract We present a fully analytic calculation of the leading-order one-loop amplitude for triple Higgs production via gluon fusion, gg → HHH , retaining full dependence on the mass of the heavy quark circulating in the loop. This amplitude provides a direct probe of the triple and quartic Higgs self-couplings, the measurement of which is a central goal of current and future colliders. The amplitude can be presented in compact form thanks to the use of analytic reconstruction techniques, based on finite-field and p -adic evaluations, multivariate partial fraction decompositions, and primary decompositions to identify common numerator factors. Although full analytic results are given in the text and in the supplementary material, the main thrust of this paper is to further test and illustrate these analytic reconstruction techniques in a concrete physical example. Our results provide a compact and efficient representation of the matrix element for this process, enabling evaluations that are more than an order of magnitude faster than existing numerical alternatives. Full analytic control of the leading-order, loop-induced amplitude is an important step towards handling more complex 2-loop or real-radiation corrections to this and related processes.