delZnm["scalar"] =
aspi^2*((cf*nm*tf)/(8*ep^2) - (5*cf*nm*tf)/(48*ep)) + 
 aspi^3*(((121*ca*cf*nm*tf)/432 - (cf^2*nm*tf)/192 + 
     cf*((-5*nh*nm)/108 - (5*nl*nm)/108 - (5*nm^2)/216)*tf^2)/ep^2 + 
   ((-11*ca*cf*nm*tf)/72 + (3*cf^2*nm*tf)/32 + 
     cf*((nh*nm)/18 + (nl*nm)/18 + nm^2/36)*tf^2)/ep^3 + 
   (cf*((-35*nh*nm)/648 - (35*nl*nm)/648 - (35*nm^2)/1296)*tf^2 + 
     ca*cf*tf*((-139*nm)/1296 - (nm*Zeta[3])/4) + 
     cf^2*tf*((-23*nm)/96 + (nm*Zeta[3])/4))/ep);
