Two-loop bottom mass effects on the Higgs transverse momentum spectrum in top-induced gluon fusion

We compute bottom mass (m(b)) corrections to the transverse momentum (q(T)) spectrum of Higgs bosons produced by gluon fusion in the regime q(T) similar to m(b) MUCH LESS-THAN m(H) at leading power in m(b)/m(H) and q(T)/m(H), where the gluons couple to the Higgs via a top loop. To this end we calculate the quark mass dependence of the transverse momentum dependent gluon beam functions (aka gluon TMDPDFs) at two loops in the framework of SCET. These functions represent the collinear matrix elements in the factorized gluon-fusion cross section for small q(T). We discuss in detail technical subtleties regarding rapidity regulators and zero-bin subtractions in the calculation of the virtual corrections present for massive quarks. Combined with the known soft function for m(b) not equal 0 our results allow to determine the resummed Higgs q(T) distribution in the top-induced gluon fusion channel at NNLL ' (and eventually (NLL)-L-3) with full dependence on m(b)/q(T). We perform a first phenomenological analysis at fixed order, where the new corrections to the massless approximation lead to percent-level effects in the peak region of the Higgs q(T) spectrum. Upon resummation they may thus be relevant for state-of-the-art precision predictions for the LHC.