peak position in the gluon dressing function $1/Z_A$, which translates to $$\begin{aligned}
 \alpha_s(20 \text{ GeV}) = 0.21\ .\end{aligned}$$ The dressing functions of the Yang-Mills three-point functions, $z_{\bar c A c},\, z_{A^3}$, are calculated momentum-dependently for a single momentum channel. The four-gluon vertex is approximated using the three-gluon vertex, see App. \[trunc:3g\]. This is a very good approximation down to semi-perturbative momenta [@Pelaez:2013cpa; @Blum:2014gna; @Eichmann:2014xya; @Binosi:2014kka; @Gracey:2014mpa; @Cyrol:2014kca], wheres deviations occur mostly for momenta where the glue gap implies already decoupling.\
*Matter sector*: We take into account the full momentum dependence of the quark propagator, parameterised by its wavefunction renormalisation $Z_q(p)$ and mass function $M_q(p)$, where we have for the current quark mass $$\begin{aligned}
M_q(20\text{ GeV}) = 1.3 \text{ MeV}, \end{aligned}$$ see App. \[trunc:qpr\] for details. The treatment of the quark-gluon vertex is of crucial importance for the whole system. Therefore we include the full momentum dependence of all eight linearly independent tensor structures ${\cal T}^{(i)}_{\bar q A q}$ of Landau gauge as described in App. \[trunc:qgl\]. Additionally we perform a gauge-invariant completion of any quark-gluon vertex tensor structure that contributes quantitatively, leading to the inclusion of two-quark–two-gluon and two-quark-three-gluon vertices that are approximated gauge-invariantly, see App. \[trunc:qgl\].

In the four-fermi 