Analytic two-loop amplitudes for qq¯→γγ and gg → γγ mediated by a heavy-quark loop

We address the analytic computation of the two-loop scattering amplitudes for the production of two photons in parton-parton scattering, mediated by loops of heavy quarks. Due to the presence of integrals of elliptic type, both partonic channels have been previously computed using semi-numerical methods. In this paper, leveraging new advances in the theory of differential equations for elliptic Feynman integrals, we derive a canonical basis for all integrals involved and compute them in terms of independent iterated integrals over elliptic and polylogarithmic differential forms. We use this representation to showcase interesting cancellations in the physical expressions for the scattering amplitudes. Furthermore, we address their numerical evaluation by producing series expansion representations for the whole amplitudes, which we demonstrate to be fast and numerically reliable across a large region of the phase space.