On the NLO QCD corrections to the production of the heaviest neutral   Higgs scalar in the MSSM
  We present a calculation of the two-loop top-stop-gluino contributions to Higgs production via gluon fusion in the MSSM. By means of an asymptotic expansion in the heavy particle masses, we obtain explicit and compact analytic formulae that are valid when the Higgs and the top quark are lighter than stops and gluino, without assuming a specific hierarchy between the Higgs mass and the top mass. Being applicable to the heaviest Higgs scalar in a significant region of the MSSM parameter space, our results complement earlier ones obtained with a Taylor expansion in the Higgs mass, and can be easily implemented in computer codes to provide an efficient and accurate determination of the Higgs production cross section. 
G. Degrassi, S. Di Vita, P. Slavich
hep-ph