Elusiveness of infrared critical exponents in Landau gauge Yang-Mills theories

We solve a truncated system of coupled Dyson-Schwinger equations for the gluon and ghost propagators in SU($N_c$) Yang-Mills theories in Faddeev-Popov quantization on a four-torus. This compact space-time manifold provides an efficient mean to solve the gluon and ghost Dyson-Schwinger equations without any angular approximations. We verify that analytically two power-like solutions in the very far infrared seem possible. However, only one of these solutions can be matched to a numerical solution for non-vanishing momenta. For a bare ghost-gluon vertex this implies that the gluon propagator is only weakly infrared vanishing, $D_{gl}(k^2) \propto (k^2)^{2\kappa -1}$, $\kappa \approx 0.595$, and the ghost propagator is infrared singular, $D_{gh}(k^2) \propto (k^2)^{-\kappa -1}$. For non-vanishing momenta our solutions are in agreement with the results of recent SU(2) Monte-Carlo lattice calculations. The running coupling possesses an infrared fixed point. We obtain $\alpha(0) = 8.92/N_c$ for all gauge groups SU($N_c$). Above one GeV the running coupling rapidly approaches its perturbative form.Comment: 30 pages, REVTEX4 style, 11 figure