number of flavors and $\Lambda_{\text{QCD}}$ denotes the corresponding QCD-$\Lambda$-scale. The data in Fig. \[peik\](b) show a slightly steeper slope at distances below $r\sqrt{\sigma}\simeq0.5$ compared to the pure gauge potential given in Ref. [@Necco:2001xg] indicating that the QCD coupling gets stronger in the entire distance range analyzed here when including dynamical quarks. This is in qualitative agreement with (\[runningcoupling\]). To include the effect of a stronger Coulombic part in the QCD potential we test the Cornell parameterization, $$\begin{aligned}
\frac{V(r)}{\sqrt{\sigma}} = -\frac{4}{3}\frac{\alpha}{r\sqrt{\sigma}} + r \sqrt{\sigma}
\label{t=0ansatz}\;,\end{aligned}$$ with a free parameter $\alpha$. From a best-fit analysis of Eq. (\[t=0ansatz\]) to the data ranging from $0.2\;\lsim\;r\sqrt{\sigma}\;\lsim\;2.6$ we find $$\begin{aligned}
\alpha&=&0.212 (3)\;.\label{res}\end{aligned}$$ This already may indicate that the logarithmic weakening of the coupling with decreasing distance will not too strongly influence the properties of the QCD potential at these distances, [*i.e.*]{} at $r\;\gsim\;0.1$ fm. However, the value of $\alpha$ is moderately larger than $\alpha_{\text{str}}\;\simeq\;0.196$ introduced above. To compare the relative size of $\alpha$ in full QCD to $\alpha$ in the quenched theory we again have performed a best-fit analysis of the quenched zero temperature potential given in [@Necco:2001xg] using the Ansatz given in Eq. (\[t=0ansatz\]) and a similar distance range. Here we find $\alpha_{\text{quenched}} = 