Two-loop leading-color QCD corrections for Higgs plus two-jet production in the heavy-top limit

We present the leading-color two-loop QCD corrections for Higgs-boson production in association with two jets through gluon fusion in the heavy-top effective theory. We provide analytic expressions for the finite remainders of the helicity amplitudes, written in terms of one-mass pentagon functions with spinor-helicity coefficients. These expressions are obtained by reconstructing the amplitudes from numerical finite-field samples computed within the numerical unitarity framework. The reconstruction is made possible by several advances in exploiting the analytic structure of the amplitudes, which both reduce the number of required samples and lead to compact representations. In particular, we introduce a new algorithm for multivariate partial fraction decomposition, based on a generic bivariate slice and a simplified treatment of ideal intersections. Using the resulting analytic expressions, we provide an efficient and stable implementation of their numerical evaluation, ready for phenomenological applications. Finally, we study the singularity structure of the remainders and confirm the existence of a threshold at non-degenerate physical momentum configurations, usually associated with massive virtual particle exchanges.