(* Mathematica script AdlerFunction4lMiniMOMexact.m *)
(* as AdlerFunction4lMiniMOM.m, but c2 and c3 of miniMOM are now exact to all digits *)
(* ************************************************************ *)

Print["AdlerFunction4lMiniMOM.m: a Mathematica script for analytic running"]
Print["It gives the coefficients d1, d2, d3 of the perturbation expansion"]
Print["for the nf=3 massless canonical Adler function d(Q^2)"]
Print["in the 4-loop Lambert-MiniMOM renormalization scheme"]
Print["as defined in paper arXiv:1703.01321"]
Print["Note that the perturbation expansion of the massless Adler function is:"]
Print["d(Q^2)_{pt} = a(Q^2) + d1 a(Q^2)^2 + d2 a(Q^2)^3 + d3 a(Q^2)^4"]
Print["where a(Q^2) = alpha_s(Q^2)/Pi in the 4-loop Lambert MiniMOM."]
Print["See Eqs. (29)-(30) of  arXiv:1703.01321"]
Print["In addition, it gives the parameters of the resummed expression"]
Print["for d(Q^2), Eq. (31) of arxiv:1703.01321:"]
Print["{tilde alpha1}=alt[2,1]; {tilde alpha2}=alt(2,2)=1-{tilde alpha1}"]
Print["kappa1= kapps(2,1); kappa2=kapps(2,2)"]
Print["be0, c1, c2MS, c3MS, c2, c3 are the RGE coefficients; the symbol explanation can be called by the command: ?symbolname"]

 (* ************************************************************ *)
(* ************************************************************ *)


(* When using the script, the symbol explanation can be called by the command: ?symbolname *)

 (* ::Section:: *)
(*Definitions*)
nf::usage= "nf is the number of active quarks in the QCD theory, here nf=3.";
be0::usage= "be0 = beta0 RGE coefficient, here be0=9/4 (nf=3)"
c1::usage= "c1 = beta1/beta0 RGE coefficient, here c1=16/9 (nf=3)"
c2MS::usage= "c2MS = beta2(MSbar)/beta0 RGE coefficient in MSbar scheme, here c2MS=3863/864=4.47106... (nf=3)"
c3MS::usage= "c3MS = beta3(MSbar)/beta0 RGE coefficient in MSbar scheme, here c3MS=20.990247..."
c2MM::usage= "c2 = beta2/beta0 RGE coefficient in 4-loop Lambert-MiniMOM scheme, here c2=9.297 (nf=3) "
c3MM::usage= "c3 = beta3/beta0 RGE coefficient in 4-loop Lambert-MiniMOM scheme, c3=c2^2/c1 = 48.650625... (nf=3)"
d1bar::usage= "the coefficient d1 in MSbar scheme"
d2bar::usage= "the coefficient d2 in MSbar scheme"
d3bar::usage= "the coefficient d3 in MSbar scheme"
td1bar::usage= "the coefficient {tilde d}1 in MSbar scheme"
td2bar::usage= "the coefficient {tilde d}2 in MSbar scheme"
td3bar::usage= "the coefficient {tilde d}3 in MSbar scheme"
d1::usage= "the coefficient d1 in 4-loop Lambert-MiniMOM scheme [cf.arXiv:1703.01321, Eqs.(29)-(30)]"
d2::usage= "the coefficient d2 in 4-loop Lambert-MiniMOM scheme [cf.arXiv:1703.01321, Eqs.(29)-(30)]"
d3::usage= "the coefficient d3 in 4-loop Lambert-MiniMOM scheme (cf.arXiv:1703.01321, Eqs.(29)-(30)]"
td1::usage= "the auxiliary coefficient {tilde d}1 in 4-loop Lambert-MiniMOM scheme"
td2::usage= "the auxiliary coefficient {tilde d}2 in 4-loop Lambert-MiniMOM scheme"
td3::usage= "the auxiliary coefficient {tilde d}3 in 4-loop Lambert-MiniMOM scheme"
alt::usage= "alt[2,1]=the coefficient {tilde alpha1}, in 4-loop Lambert-MiniMOM scheme; alt[2,1]={tilde alpha2}=1-{tilde alpha1}  [cf. arXiv:1703.01321, Eq.(31)]"
kapps::usage= "kapps[2,1]=the coefficient kappa1, in 4-loop Lambert-MiniMOM scheme; kapps[2,2]=kappa2  [cf.arXiv:1703.01321, Eq.(31)]"
(*--------------------------------------------------------------------------_*)       

nf = 3
 
be0 = 9/4
 
c1 = 16/9
 
c2MS = 3863/864
 
c3MS = 20.99023981042311
 
c2MM = 9.297028230447523
 
c3MM = 71.45379490569553
 
d1bar = 1/12 + (9*(11/2 - 4*Zeta[3]))/4
 
d2bar = -12521/288 + (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
     (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16
 
d3bar = 144939499/20736 - (5693495*Zeta[3])/864 + (5445*Zeta[3]^2)/8 + 
     9*(1045381/15552 - (40655*Zeta[3])/864 + (5*Zeta[3]^2)/6 - 
       (260*Zeta[5])/27) + 27*(-6131/5832 + (203*Zeta[3])/324 + 
       (5*Zeta[5])/18) + (65945*Zeta[5])/288 - (7315*Zeta[7])/48 + 
     3*(-13044007/10368 + (12205*Zeta[3])/12 - 55*Zeta[3]^2 + 
       (29675*Zeta[5])/432 + (665*Zeta[7])/72)
 
td1bar = 1/12 + (9*(11/2 - 4*Zeta[3]))/4
 
td2bar = -12521/288 - (16*(1/12 + (9*(11/2 - 4*Zeta[3]))/4))/9 + 
     (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
     (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16
 
td3bar = 144939499/20736 + (8891*(1/12 + (9*(11/2 - 4*Zeta[3]))/4))/2592 - 
     (5693495*Zeta[3])/864 + (5445*Zeta[3]^2)/8 + 
     9*(1045381/15552 - (40655*Zeta[3])/864 + (5*Zeta[3]^2)/6 - 
       (260*Zeta[5])/27) + 27*(-6131/5832 + (203*Zeta[3])/324 + 
       (5*Zeta[5])/18) + (65945*Zeta[5])/288 - 
     (40*(-12521/288 + (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
        (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16))/9 - (7315*Zeta[7])/48 + 
     3*(-13044007/10368 + (12205*Zeta[3])/12 - 55*Zeta[3]^2 + 
       (29675*Zeta[5])/432 + (665*Zeta[7])/72)
 
d1 = 1.6398212048969851
 
d2 = 1.545051067468254
 
d3 = 8.016491849146206
 
td1 = 1.6398212048969851
 
td2 = -1.370186630126386
 
td3 = -1.1392538295502659
 


  alt[1, 1] = 1.
 
alt[2, 1] = 0.5000000000000000000000000 - 0.9340290391494241481812213 I
   
alt[2, 2] =0.5000000000000000000000000 + 0.9340290391494241481812213 I
   
kapps[1, 1] = 0.482483080857287269382051734330
   
kapps[2, 1] = 0.96903758942788972533271187 - 0.43577379898357130601420551 I
   
kapps[2, 2] =0.96903758942788972533271187 + 0.43577379898357130601420551 I
  
