On the next-to-next-to-leading order evolution of flavour-singlet fragmentation functions

We present the third-order contributions to the quark–gluon and gluon–quark timelike splitting functions for the evolution of fragmentation functions in perturbative QCD. These quantities have been derived by studying physical evolution kernels for photon- and Higgs-exchange structure functions in deep-inelastic scattering and their counterparts in semi-inclusive annihilation, together with constraints from the momentum sum rule and the supersymmetric limit. For this purpose we have also calculated the second-order coefficient functions for one-hadron inclusive Higgs decay in the heavy-top limit. A numerically tolerable uncertainty remains for the quark–gluon splitting function, which does not affect the endpoint logarithms for small and large momentum fractions. We briefly discuss these limits and illustrate the numerical impact of the third-order corrections. Compact and accurate parametrizations are provided for all third-order timelike splitting functions.