str=Directory[];
SetDirectory["/home/blumlein/OLDFILES/LADDER"];
<<Sigma.m;
<<HarmonicSums.m;
SetDirectory[str]

2-loop asymptotic massive polarized Wilson Coefficients in z space (single mass case)
In any use of these functions the following paper has to be qouted:

I.~Bierenbaum, J.~Bl \" umlein, A.~De Freitas, A.~Goedicke, S.~Klein and K.~Sch \" onwald,  $O (\alpha_s^2 $) Polarized Heavy Flavor Corrections to Deep - Inelastic Scattering at $Q^2 \gg  m^2 $,  [arXiv : 2211.15337[hep - ph]].

The results are expressed of harmonic  polylogarithms, with LL = Log[Q^2/m^2] and m the heavy quark mass.

Non-singlet Wilson coefficient

ΔLNS2=CF TF (-(4/3) (LM^2) (1+x)+(4/27) (-47-218 x+36 z2+36 x z2)+((265/9)+((2 LM)/3)+2 (LM^2)) (Subscript[δ, 1-x])-4 (1+x) ((H[0,x])^2)+LM (-(8/9) (-1+11 x)-(8/3) (1+x) H[0,x])-(8/9) (5+14 x) H[1,x]-(4/3) (1+x) ((H[1,x])^2)+H[0,x] (-(8/9) (13+28 x)-(8/3) (1+x) H[1,x])-(8/3) (1+x) H[0,1,x]+(1/xm)(((8 (LM^2))/3)-(2/27) (-359+144 z2)+8 ((H[0,x])^2)+LM ((80/9)+(16/3) H[0,x])+(116/9) H[1,x]+(8/3) ((H[1,x])^2)+H[0,x] ((268/9)+(16/3) H[1,x])+(16/3) H[0,1,x]))/.x->z/. 1/xm->1/(1-z)/.LM->LL/.CF->4/3/.TF->1/2;

Pure singlet Wilson coefficient

ΔHPS2=CF TF (-(592/3) (-1+z)-(32/3) (9-3 z+(z^2)) z2+16 (1+z) z3+(-(256/3) (-2+z)-32 (1+z) z2) H[0,z]-((32 (((1+z))^3) H[-1,z] H[0,z])/(3 z))+(8/3) (21+2 (z^2)) ((H[0,z])^2)+(16/3) (1+z) ((H[0,z])^3)+(LM^2) (20 (-1+z)-8 (1+z) H[0,z])+LM (-8 (-1+z)+8 (-1+3 z) H[0,z]-8 (1+z) ((H[0,z])^2))+(-88 (-1+z)-80 (-1+z) H[0,z]) H[1,z]-20 (-1+z) ((H[1,z])^2)+((32 (((1+z))^3) H[0,-1,z])/(3 z))+(16 (-1+3 z)+32 (1+z) H[0,z]) H[0,1,z]-32 (1+z) H[0,0,1,z]+16 (1+z) H[0,1,1,z])/.CF->4/3/.TF->1/2/.LM->LL;

Gluonic Wilson coefficient coupling to the heavy quark

ΔHQG2=(LL^2) (-(16/3) (TF^2) (-1+2 z)+CF TF (6+(-1+2 z) (-4 H[0,z]-8 H[1,z]))+CA TF (48 (-1+z)-16 (1+z) H[0,z]+8 (-1+2 z) H[1,z]))+LL ((TF^2) ((16/3) (-3+4 z)+(-1+2 z) ((16/3) H[0,z]+(16/3) H[1,z]))+CA TF (-8 (-12+11 z)-16 z2+8 (1+8 z) H[0,z]-32 (-1+z) H[1,z]-8 (-1+2 z) ((H[1,z])^2)+(1+2 z) (-16 H[-1,z] H[0,z]-8 ((H[0,z])^2)+16 H[0,-1,z]))+CF TF (4 (-17+13 z)-24 (-1+2 z) z2+16 (-3+2 z) H[0,z]+4 (-17+20 z) H[1,z]+(-1+2 z) (8 ((H[0,z])^2)+32 H[0,z] H[1,z]+16 ((H[1,z])^2)-8 H[0,1,z])))+CA TF (-(8/3) (-101+104 z)-8 (-1-10 z+4 (z^2)) z3+(4/3) (194-163 z+6 (z^2)) H[0,z]-((16 (2+3 z+9 (z^2)+11 (z^3)) H[-1,z] H[0,z])/(3 z))+16 (z^2) ((H[-1,z])^2) H[0,z]+((2/3) (126-48 z+41 (z^2))-4 (-3-6 z+2 (z^2)) H[-1,z]) ((H[0,z])^2)+(8/3) (3+4 z) ((H[0,z])^3)+4 (43-53 z+2 (z^2)) H[1,z]-4 (-53+56 z+(z^2)) H[0,z] H[1,z]+4 (3-6 z+2 (z^2)) ((H[0,z])^2) H[1,z]+2 (19-24 z+(z^2)) ((H[1,z])^2)+z2 (-(4/3) (114-84 z+47 (z^2))+16 (-1-2 z+(z^2)) H[-1,z]-32 (2+z) H[0,z]-16 (((-1+z))^2) H[1,z])+(((16 (2+3 z+9 (z^2)+11 (z^3)))/(3 z))-32 (z^2) H[-1,z]+8 (-3-6 z+2 (z^2)) H[0,z]) H[0,-1,z]+4 (-19+28 z+2 (z^2)) H[0,1,z]-8 (-7-10 z+2 (z^2)) H[0,z] H[0,1,z]+(-1+2 z) (-16 H[0,z] ((H[1,z])^2)+16 H[1,z] H[0,1,z])+32 (z^2) H[0,-1,-1,z]-8 (-3-6 z+2 (z^2)) H[0,0,-1,z]+8 (-9-10 z+2 (z^2)) H[0,0,1,z]+(1+2 z) (16 H[-1,z] H[0,1,z]-16 H[0,-1,1,z]-16 H[0,1,-1,z])+48 H[0,1,1,z])+CF TF (-(20/3) (-20+17 z)+8 (1+14 z+8 (z^2)) z3-(8/3) (-46+53 z+6 (z^2)) H[0,z]+((16 (4+12 (z^2)+13 (z^3)) H[-1,z] H[0,z])/(3 z))-32 (((1+z))^2) ((H[-1,z])^2) H[0,z]+(-(4/3) (-27+6 z+23 (z^2))+16 (((1+z))^2) H[-1,z]) ((H[0,z])^2)-4 (-47+41 z+4 (z^2)) H[1,z]+8 (8-14 z+(z^2)) H[0,z] H[1,z]-16 (z^2) ((H[0,z])^2) H[1,z]-2 (-33+40 z+2 (z^2)) ((H[1,z])^2)+z2 ((8/3) (-60+42 z+29 (z^2))-32 (((1+z))^2) H[-1,z]+32 (z^2) H[1,z]+(-1+2 z) (32 H[0,z]+16 H[1,z]))+(-((16 (4+12 (z^2)+13 (z^3)))/(3 z))+64 (((1+z))^2) H[-1,z]-32 (((-1+z))^2) H[0,z]) H[0,-1,z]-16 (-6+(z^2)) H[0,1,z]+32 (((-1+z))^2) H[0,z] H[0,1,z]-64 (((1+z))^2) H[0,-1,-1,z]+32 (1-6 z+(z^2)) H[0,0,-1,z]-32 (((-1+z))^2) H[0,0,1,z]+(-1+2 z) (-(8/3) ((H[0,z])^3)-16 H[0,z] ((H[1,z])^2)-8 ((H[1,z])^3)-16 H[1,z] H[0,1,z]+24 H[0,1,1,z]))/.CF->4/3/.TF->1/2/.CA->3;

Gluonic Wilson coefficient coupling to the light quarks

ΔLSG2=16/3*TF^2*(4*z-3+(2*z-1)*(H[0,z]+H[1,z]))*LL/.TF->1/2;