We present a fully differential computation of e + e − → γ γ at next-to-next-to-leading order in QED. The process has been implemented into McMule , completing its set of next-to-next-to-leading-order calculations for the most important 2 → 2 processes. The results allow for generic applications to electron-positron colliders with centre-of-mass energies up to a few GeV, particularly for luminosity measurements.
Engel et al. (Mon,) studied this question.