Calculation of higher-twist evolution kernels for polarized deep inelastic scattering

Based on the non-local light-ray operator technique, we develop an algorithm for the computational calculation of evolution kernels for higher-twist operators in leading order of the perturbation theory. We compute the evolution kernel for the twist-3 operators in the flavour singlet channel. Our result confirms the local anomalous dimensions computed by Bukhvostov, Kuraev, and Lipatov as well as the non-local evolution kernels of Balitzky and Braun.