M{\o}ller scattering at NNLO

We present a calculation of the full set of next-to-next-to-leading-order QED corrections to unpolarised M{\o}ller scattering. This encompasses photonic, leptonic, and non-perturbative hadronic corrections and includes electron mass effects as well as hard photon radiation. The corresponding matrix elements are implemented in the Monte Carlo framework McMule allowing for the computation of fully-differential observables. As a first application we show results tailored to the kinematics and detector design of the PRad II experiment where a high-precision theory prediction for M{\o}ller scattering is required to achieve the targeted precision. We observe that the corrections become essential to reliably calculate the corresponding differential distributions especially in regions where the leading-order contribution is absent.