Higgs Boson pair production at next-to-next-to-leading order in QCD

We compute the next-to-next-to-leading order QCD corrections for standard model Higgs boson pair production inclusive cross section at hadron colliders within the large top-mass approximation. We provide numerical results for the LHC, finding that the corrections are large, resulting in an increase of O(20%) with respect to the next-to-leading order result at c.m. energy sqrt[sH]=14 TeV. We observe a substantial reduction in the scale dependence, with overlap between the current and previous order prediction. All our results are normalized using the full top- and bottom-mass dependence at leading order. We also provide analytical expressions for the K factors as a function of sH.