A bstract We present the electroweak corrections for the production of a photon pair through gluon fusion, focusing on the contribution from the first two generations of quarks. The two-loop amplitude is calculated using a series of projection operators which define scalar form factors. In order to evaluate the Master Integrals which appear in this process we employ both generalized polylogarithms and Chen-iterated integrals. In order to perform a phenomenological study we develop a semi-numerical evaluation of the Master Integrals employing a fitting procedure to speed up the evaluation of burdensome higher weight contributions. We present results for the LHC, finding corrections of around a couple of percent to the leading order gg → γγ process. Our results are implemented into the parton-level Monte Carlo code MCFM.
Fiore et al. (Thu,) studied this question.