Two-loop corrections to the triple Higgs boson production cross section
  In this paper we compute the QCD corrections for the triple Higgs boson production cross section via gluon fusion, within the heavy-top approximation. We present, for the first time, analytical results for the next-to-leading order corrections, and also compute the soft and virtual contributions of the next-to-next-to-leading order cross section. We provide predictions for the total cross section and the triple Higgs invariant mass distribution. We find that the QCD corrections are large at both perturbative orders, and that the scale uncertainty is substantially reduced when the second order perturbative corrections are included. 
Daniel de Florian and Javier Mazzitelli
hep-ph