Real-virtual corrections to Higgs boson pair production at NNLO: three closed top quark loops

We compute the real-radiation corrections, i.e. the virtual corrections to the single real emission of a parton, to Higgs boson pair production at next-to-next-to-leading order in QCD, in an expansion for large top quark mass. We concentrate on the radiative corrections to the interference contribution from the next-to-leading order one-particle reducible and the leading order amplitudes. This is a well defined and gauge invariant subset of the full real-virtual corrections to the inclusive cross section. We obtain analytic results for all phase-space master integrals both as an expansion around the threshold and in an exact manner in terms of Goncharov polylogarithms. We demonstrate that for many applications it is sufficient to use the expanded expressions.