Estimate of Γ(H to gg) at O(α_s^7) from Padé approximants with constraints from large-n_f

In this work we estimate the first unknown QCD correction to the decay of the Higgs boson into gluons, Hto gg, in the heavy-top limit using Padé approximants built to the perturbative series and to its Borel transform. In the process, we also obtain estimates for the first unknown coefficients of the QCD β function, of the QCD decoupling relation, and of the gluonium correlator. We also provide an updated prediction for the first unknown coefficient of Γ(Hto bb). We determine the polynomial dependence of the perturbative coefficients on the number of quark flavors, n_f. We point out that the large-n_f result for Hto gg can be obtained from results for the gluonium correlator already available in the literature. The coefficient of the highest power of n_f is then imposed as a constraint and is also used as a new test to discard unreliable approximants. The method is validated through the postdiction of the last known coefficient of the series. Our model-independent estimate for the yet unknown coefficient of order α_s^7 of Γ(Hto gg) is -312pm 15 (in the mathrm{MS} scheme with μ=m_H, where m_H is the Higgs mass). With our estimate the residual renormalization-scale dependence of Γ(Hto gg) is reduced by a factor of almost 1.5.