Two-loop Feynman integrals for leading colour tt¯W production at hadron colliders

We compute a complete set of the two-loop Feynman integrals that are required for the next-to-next-to-leading order QCD corrections to on-shell top-pair production in association with a W boson at hadron colliders in the leading colour approximation. These Feynman integrals also contribute to Higgs or Z-boson production in association with a top pair. We employ the method of differential equations (DEs), facilitated by the use of finite field methods to handle the algebraic complexity stemming from the seven-scale kinematics. The presence of the top quark in the virtual propagators, in addition to the mass of the external W boson, gives rise to nested square roots and three elliptic curves. We obtain DEs that depend at most quadratically on the dimensional regulator epsilon for sectors where these analytic structures appear, and are epsilon-factorised otherwise. We express the DEs in terms of a minimal set of differential one-forms, separating the logarithmic ones. We solve the DEs numerically in the physical kinematic region, with the method of generalised power series expansions.