  We compute the two-loop QCD helicity amplitudes for the production of a Higgs
boson in association with a bottom-quark pair at a hadron collider. We take the
approximations of leading colour and work in the five flavour scheme, where the
bottom quarks are massless while the Yukawa coupling is non-zero. We extract
analytic expressions from multiple numerical evaluations over finite fields and
present the results in terms of an independent set of special functions that
can be reliably evaluated over the full phase space.
