Two-loop QCD corrections to $H\to b+\bar{b}+g$ at higher powers in the dimensional regulator

We compute the two-loop massless QCD corrections to the amplitude of Higgs boson decay to bottom quark pair and a gluon ($H\to b+\bar{b}+g$) in the higher powers of the dimensional regularization parameter $\epsilon$. The calculation is performed by projecting the amplitude onto the appropriate Lorentz structures related to the process. We also show the numerical behaviour of the form factors for a few sample phase-space points. These amplitudes are necessary ingredients for computing the three-loop virtual corrections to bottom-quark annihilation to Higgs plus jet production at the hadron collisions.