Renormalization constants for Lattice QCD: new results from Numerical Stochastic Perturbation Theory

By making use of Numerical Stochastic Perturbation Theory (NSPT) we can compute renormalization constants for Lattice QCD to high orders, e.g. three or four loops for quark bilinears. We report on the status of our computations, which provide several results for Wilson quarks and in particular (values and/or ratios of) Z_V, Z_A, Z_S, Z_P. Results are given for various number of flavors (n_f = 0, 2, 3, 4). While we recall the care which is due for the computation of quantities for which an anomalous dimension is in place, we point out that our computational framework is well suited to a variety of other calculations and we briefly discuss the application of NSPT to other regularizations (in particular the Clover action).