Two-Loop Four-Gluon Amplitudes with the Numerical Unitarity Method

We present the first numerical computation of two-loop amplitudes based on the unitarity method. As a proof of principle, we compute the four-gluon process. We discuss the new method, analyze its numerical properties and apply it to reconstruct the analytic form of the amplitudes. The numerical method is universal, and can be automated to provide multi-scale two-loop computations for phenomenologically relevant signatures at hadron colliders.