Perturbative Wilson loops from unquenched Monte Carlo simulations at weak couplings

Perturbative expansions of several small Wilson loops are computed through next-to-next-to-leading order in unquenched lattice QCD, from Monte Carlo simulations at weak couplings. This approach provides a much simpler alternative to conventional diagrammatic perturbation theory, and is applied here for the first time to full QCD. Two different sets of lattice actions are considered: one set uses the unimproved plaquette gluon action together with the unimproved staggered-quark action; the other set uses the one-loop-improved Symanzik gauge-field action together with the so-called ``asqtad'' improved-staggered quark action. Simulations are also done with different numbers of dynamical fermions. An extensive study of the systematic uncertainties is presented, which demonstrates that the small third-order perturbative component of the observables can be reliably extracted from simulation data. We also investigate the use of the rational hybrid Monte Carlo algorithm for unquenched simulations with unimproved-staggered fermions. Our results are in excellent agreement with diagrammatic perturbation theory, and provide an important cross-check of the perturbation theory input to a recent determination of the strong coupling $\alpha_{\bar{\rm MS}}(M_Z)$ by the HPQCD collaboration.