we find it necessary to include the full momentum dependence of the corresponding dressing functions $z_{\bar q A q}^{(i)}(p,q)$. It is only the combination of all these terms, including their gauge invariant partner structures in the quark-gluon vertex equation together with their momentum dependencies, that leads to the very good agreement with the lattice propagator. In particular this concerns the wave function renormalisation $Z_q(p)$ for small momenta, where an important contribution stems from mesonic fluctuations in the infrared. These fluctuations have been included with functional methods, e.g. in [@Jungnickel:1995fp; @Berges:2000ew; @Schaefer:2006ds; @Braun:2009gm; @Mitter:2013fxa; @Alkofer:1993gu; @Fischer:2007ze; @Fischer:2008sp; @Cloet:2008re]. Restricting the discussion only to the relative importance of quark-gluon vertex tensor structures, recent findings in Dyson-Schwinger studies [@Hopfer:2013np; @Williams:2014iea; @Aguilar:2014lha] agree with our, see also [@Alkofer:2008tt] for an earlier study. Moreover, our present findings suggest the inclusion of the STI-consistent higher quark-gluon interactions in future DSE-studies.

Finally we want to point out that one crucial contribution to the quark mass function comes from the addition to the flow of the Yukawa coupling, $\partial_t\Delta h$, due to dynamical hadronisation, see . As soon as one runs into the spontaneously broken phase of QCD, $\langle
\sigma\rangle\neq 0$, this term contributes to the quark mass 