We study the low-momentum behaviour of Yang-Mills propagators obtained
from Landau-gauge Dyson-Schwinger equations (DSE) in the PT-BFM scheme. We compare
the ghost propagator numerical results with the analytical ones obtained by analyzing
the low-momentum behaviour of the ghost propagator DSE in Landau gauge, assuming for
the truncation a constant ghost-gluon vertex and a simple model for a massive gluon propagator.
The asymptotic expression obtained for the regular or decoupling ghost dressing
function up to the order O(q2) is proven to fit pretty well the numerical PT-BFM results.
Furthermore, when the size of the coupling renormalized at some scale approaches some
critical value, the numerical PT-BFM propagators tend to behave as the scaling ones. We
also show that the scaling solution, implying a diverging ghost dressing function, cannot
be a DSE solution in the PT-BFM scheme but an unattainable limiting case.The author is particularly indebted to Ph. Boucaud, J.P Leroy, A. Le Yaouanc, J. Micheli and O. Pene for very fruitful discussions at the initial stages of the work and to J. Papavassiliou and A. C. Aguilar also for very valuable discussions and comments, and specially for providing me with some unpublished results which were exploited in this paper. J. R-Q also acknowledges the Spanish MICINN for the support by the research project FPA2009-10773 and "Junta de Andalucia" by P07FQM02962