A bstract In this paper, we present analytical results for the two-loop QCD corrections to the production of two partons or a photon and a parton in hadronic collisions, mediated by loops of massive quarks. These amplitudes involve Feynman integrals defined on an elliptic curve. We compute them by generalizing our recent results for the production of two photons to include additional crossings of the corresponding master integrals, which we compute in terms of the same basis of independent iterated integrals. We discuss the analytical properties of the amplitudes, highlighting the cancellations of a large number of elliptic differential forms in their finite remainders. Finally, we elaborate on a strategy for their numerical evaluation based on generalized series expansions at singular points of the physical amplitude, through the introduction of suitable sets of variables that allow us to resolve all singularities.
Coro et al. (2026) studied this question.