Next-to-next-to-leading order evolution of polarized parton densities in the Larin scheme

In many calculations involving polarized twist-2 parton densities to higher order in the strong coupling constant one uses the Larin scheme to describe chiral effects in dimensional regularization. Upon forming observables, the scheme dependence cancels. Still one needs a corresponding regularization scheme to compute the contributing building blocks, like massless and massive Wilson coefficients, as well as the massive three-loop operator matrix elements used in the variable flavor number scheme. These are matched to the evolved parton distribution functions in the Larin scheme. Starting with suitable input distributions we provide the solution of scale evolution of the different polarized parton distribution functions in Bjorken <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>x</a:mi></a:math> space for a wide range of virtualities <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:msup><c:mi>Q</c:mi><c:mn>2</c:mn></c:msup></c:math> in the Larin scheme, at next-to-leading, at next-to-next-to-leading order for the first time. We also illustrate the deviation between the parton distributions in the Larin and <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mover accent="true"><e:mi>MS</e:mi><e:mo stretchy="true">¯</e:mo></e:mover></e:math> schemes numerically.