Toward the QCD $\beta$ Function with Dynamical Wilson Fermions
  We present data for the scaling behavior of lattice QCD with two flavors of light Wilson fermions. This is done by matching $\pi$ n and $\rho$ masses at the two lattice sizes of $16^3\times32 $ and $8^3\times16$. We find that at $\beta=6/g^2\equiv5.3$ the matching does not occur over a range extending down to $\beta=3.5$. For $\beta=5.5$ matching may be achieved at $\beta=4.8-4.9$, leading to a $\Delta \beta =0.6-0.7$ which is higher than the perturbative 2-loop value of 0.45. In both cases we conclude that the simulations are very far from the perturbative scaling region. 
K.M. Bitar, R.G. Edwards, U.M. Heller, A.D. Kennedy and P. Vranas