Given the complete supersymmetric lagrangian of a free abelian gauge multiplet $$ \mathcal{L} = -\frac{1}{4} F_{\mu\nu} F^{\mu\nu} + i \bar{\lambda} \bar{\sigma}^\mu \partial_\mu \lambda + \frac{1}{2} D^2 + \xi D $$ With the following variations $$ \delta_\zeta A_\mu = \bar{\zeta} \bar{\sigma}^\mu \lambda + \bar{\lambda} \bar{\sigma}^\mu \zeta $$ $$ \delta_\zeta \lambda_\alpha = \frac{i}{2} F_{\mu\nu} (\sigma^\mu \bar{\sigma}^\nu \zeta)_\alpha + D \zeta_\alpha $$ $$ \delta_\zeta D = -i \bar{\zeta} \bar{\sigma}^\mu \partial_\mu \lambda + i \partial_\mu \bar{\lambda} \bar{\sigma}^\mu \zeta $$
Where $\zeta_\alpha$ is an infinitesimal Grassmann number.
Now to show that the algebra closes I took the commutator of the variations, acting on each of the fields.
Showing that the algebra closes for $\lambda$ and $D$ is relatively start forward algebra, resulting in $$ [\delta_\psi, \delta_\zeta] \lambda_\alpha = -2i(\bar{\zeta} \bar{\sigma}^\mu \psi - \bar{\psi} \bar{\sigma}^\mu \zeta)\partial_\mu \lambda_\alpha $$
$$ [\delta_\psi, \delta_\zeta] D = -2i(\bar{\zeta} \bar{\sigma}^\mu \psi - \bar{\psi} \bar{\sigma}^\mu \zeta)\partial_\mu D $$
However when I do it for the gauge fields I find, $$ [\delta_\psi, \delta_\zeta] A_\mu = -2i(\bar{\zeta} \bar{\sigma}^\alpha \psi - \bar{\psi} \bar{\sigma}^\alpha \zeta)(\partial_\alpha A_\mu - \partial_\mu A_\alpha) $$
Which almosts looks like the translation that one would expect. However it has that extra piece there. Is there something that makes that $\partial_\mu A_\alpha$ term zero that I am missing?
In http://arxiv.org/abs/hep-th/0109172 Jose suggests (pg 12) that the auxiliary fields are added to close the algebra off-shell, which in the case of the gauginos it did. He also mentions on (pg 21) that in the case of having gauge fields the algebra will only close up to a gauge transformation. Now he did that before adding auxiliary fields, so maybe it isn't true...
However if this is true, it seems odd to me that the algebra would only close up to a gauge transformation. To me that would imply that gauge fixing is required and that doesn't really seem kosher. For we would need to choose a gauge for the theory to be valid, otherwise we don't have a true gauge symmetry. How do we get around that problem?
EDIT:
@Qmechanic I believe a good reference is http://arxiv.org/abs/hep-ph/9709356, just a note in advanced he uses the mostly plus convention where I use the mostly minus, essentially the translations and sigma matrices just differ by an overall minus sign, but the spirit is there. I would look at sections 2 and 3 in particular section 3.3 (for the abelian case set a=1 and the covariant derivatives just become standard derivatives). Also note that I absorbed the factor $\frac{1}{2}$ into the infinitesimal parameter and he's choice of sign is different on the variation of the gauge field.
Notes of Convention: My conventions are nearly identical to the SUSY Primer reference except for the metric. $$ \eta^{\mu\nu} = diag(+ - - -) $$
Clarification of the Question Provided that my arithmetic is correct, how can one see that the $\partial_\mu A_\alpha$ term dies?
After that, even without checking my work (thank you if you do!), why is it for the non-abelian case equation (3.3.9), do we see a covariant derivative? If the algebra were to truly be closed I would expect that to be a partial derivative which are the spacetime translations that we expect to result (see eq 3.1.29). If we should see a covariant derivative, then why?
Additional Comment
It turns out that the turn doesn't die and is indeed correct. It gives the form $$ [\delta_\psi, \delta_\zeta] A_\mu^a = -2i(\bar{\zeta} \bar{\sigma}^\alpha \psi - \bar{\psi} \bar{\sigma}^\alpha \zeta)F_{\alpha\mu}^a $$
Which holds in the case of the Non-Abelian case (the other terms go from being of the form in (3.3.9) of SUSY Primer)
In Freedman's and De Wit's paper journals.aps.org/prd/abstract/10.1103/PhysRevD.12.2286 they discuss the implications of what having a gauge symmetry does to the commutator of two covariant SUSY transformations. It now generates a new transformation which can be described as a translation plus a gauge transformation, or what they call a gauge-covariant translation.
How is this considered closing the graded lie algebra, which should consist of Translations, Boosts/Rotations, and a SUSY transformation?