  The Bethe-Salpeter equation for a pseudoscalar bound-system, with i) a ladder
kernel with massive gluons, ii) dynamically-dressed quark mass function and
iii) an extended quark-gluon vertex, is solved in Minkowski space by using the
Nakanishi integral representation of the Bethe-Salpeter amplitude. The quark
dressing is implemented through a phenomenological ansatz, which was tuned by
lattice QCD calculations of the quark running mass. The latter were also used
for assigning the range of the gluon mass and the parameter featuring the
extended color density. This framework allows to investigate the gluon dynamics
that manifest itself in the quark dressing, quark-gluon vertex and the binding,
directly in the physical space. We present the first results for low-density
pseudoscalar systems in order to elucidate the onset of the interplay between
the above mentioned gluonic phenomena, and we discuss both static and dynamical
quantities, like valence longitudinal and transverse distributions.
