We have computed the even-N moments Nβ€20 of the gluon-to-quark
splitting function Pqgβ at the fourth order of perturbative QCD via the
renormalization of off-shell operator matrix elements. Our results, derived
analytically for a general gauge group, agree with all results obtained for
this function so far, in particular with the lowest five moments obtained via
physical cross sections. Using our new moments and all available endpoint
constraints, we construct approximations for the four-loop Pqgβ(x) that
should be sufficient for a wide range of collider-physics applications. The
N3LO corrections resulting from these and the corresponding quark-quark
splitting functions lead to a marked improvement of the perturbative accuracy
for the scale derivative of the singlet quark distribution, with effects of 1%
or less at xβ³10β4 at a standard reference scale with Ξ±sβ=0.2.Comment: 17 pages latex, 3 figures, 2 ancillary files (FORM file with results
and FORTRAN subroutine