Higgs boson production at the LHC: NNLO partonic cross sections through   order $\epsilon$ and convolutions with splitting functions to N$^3$LO
  We consider Higgs boson production at hadron colliders in the gluon fusion channel and compute higher order terms in the regularization parameter $\epsilon$. In particular, the next-to-next-to-leading order cross section is evaluated including order $\epsilon$ terms. These results are used to compute all convolutions with the splitting functions entering the next-to-next-to-next-to-leading order cross section. 
Maik H\"oschele, Jens Hoff, Alexey Pak, Matthias Steinhauser, Takahiro   Ueda
hep-ph