QCD corrections to tri-boson production

We present a computation of the next-to-leading order QCD corrections to the production of three Z bosons at the LHC. We calculate these corrections using a completely numerical method that combines sector decomposition to extract infrared singularities with contour deformation of the Feynman parameter integrals to avoid internal loop thresholds. The NLO QCD corrections to pp -> ZZZ are approximately 50%, and are badly underestimated by the leading order scale dependence. However, the kinematic dependence of the corrections is minimal in phase space regions accessible at leading order.