NLO QCD corrections to pseudoscalar Higgs production in the MSSM
  We present a calculation of the two-loop quark-squark-gluino contributions to pseudoscalar Higgs boson production via gluon fusion in the MSSM. We regularize the loop integrals using the Pauli-Villars method, and obtain explicit and compact analytic results based on an expansion in the heavy particle masses. Our results - valid when the pseudoscalar Higgs boson is lighter than squarks and gluinos - can be easily implemented in computer codes for an efficient and accurate determination of the pseudoscalar production cross section. 
G. Degrassi, S. Di Vita and P. Slavich
hep-ph