I've got a question concerning non-Abelian supersymmetric gauge theories.
Consider supersymmetric non-Abelian theory realized on chiral superfields $\Phi_i$ in a representation $R$ with matrix generators $T_{i}^{aj}$. Let us define supergauge transformation as $$\Phi_i \rightarrow (e^{2\imath g_a \Omega^a T^a})_{i}{}^{j} \, \Phi_j.$$ The supergauge-invariant term in lagrangian is $$\mathcal{L} = \Bigl[\Phi^{*i}\,(e^V)_i{}^j \, \Phi_j\Bigr]_D.$$ For this to be gauge-invariant, the non-Abelian gauge transformation for the vector field must be $$e^V \rightarrow e^{\imath \Omega^\dagger}\,e^V\,e^{-\imath \Omega}.$$ Using Baker-Hausdorff formula, we obtain $$V^a \rightarrow V^a + \imath(\Omega^{a*}-\Omega^a)+g_a \, f^{abc}\,V^b(\Omega^{c*}+\Omega^c)+...$$ Usually at this moment they argue that since the second term on the right side does not depend on $V^a$, one can always do a supergauge transformation to Wess-Zumino gauge by choosing $\Omega^{a*}-\Omega^a$ appropriately.
This is the moment that I don't get. What does it mean? Strictly speaking, the latter expression is complicated non-linear equation on components of $V^a$ superfield.
I guess they mean, that since the second term on r.h.s. doesn't depend on $V^a$, it's possible to solve it within the framework of perturbation theory in the coupling constant(s) $g_a$. Is it correct? If so, how to prove it strictly in all orders?