Perturbation theory vs. simulation for tadpole improvement factors in pure gauge theories

We calculate the mean link in Landau gauge for Wilson and improved SU(3) anisotropic gauge actions, using two loop perturbation theory and Monte Carlo simulation employing an accelerated Langevin algorithm. Twisted boundary conditions are employed, with a twist in all four lattice directions considerably improving the (Fourier accelerated) convergence to an improved lattice Landau gauge. Two loop perturbation theory is seen to predict the mean link extremely well even into the region of commonly simulated gauge couplings and so can be used remove the need for numerical tuning of self-consistent tadpole improvement factors. A three loop perturbative coefficient is inferred from the simulations and is found to be small. We show that finite size effects are small and argue likewise for (lattice) Gribov copies and double Dirac sheets.