Dirac-Operatoridentität, die Gammamatrizen trennt

Das versuche ich zu zeigen

$$(i\gamma^\mu\partial_\mu -e\gamma^\mu A_\mu)^2 = (i\partial_\mu -eA_\mu)^2 -\frac{e}{2} \sigma^{\mu\nu}F_{\mu\nu},$$ Wo $\sigma^{\mu\nu} = \frac{i}{2}[\gamma^\mu,\gamma^\nu]$.

Bisher habe ich

$$\begin{align} (i\gamma^\mu\partial_\mu -e\gamma^\mu A_\mu)^2 &= \gamma^\mu\gamma^\nu(i\partial_\mu-e A_\mu)(i\partial_\nu -e A_\nu) =(i\partial_\mu)^2+(eA_\mu)^2 - ie\gamma^\mu\gamma^\nu(\partial_\mu A_\nu + A_\mu\partial_\nu) \\ &=(i\partial_\mu-eA_\mu)^2+ie\partial_\mu A^\mu +ie A_\mu\partial^\mu - \frac{1}{2}ie(\gamma^\mu\gamma^\nu+2\eta^{\mu\nu}-\gamma^\nu\gamma^\mu)(\partial_\mu A_\nu + A_\mu\partial_\nu) \\ &= (i\partial_\mu-eA_\mu)^2- \frac{1}{2}ie(\gamma^\mu\gamma^\nu-\gamma^\nu\gamma^\mu)(\partial_\mu A_\nu + A_\mu\partial_\nu) \\ &= (i\partial_\mu-eA_\mu)^2 - e\sigma^{\mu\nu}(\partial_\mu A_\nu + A_\mu\partial_\nu), \end{align}$$ das ist nah, aber nicht ganz da. Das kann ich zeigen $\sigma^{\mu\nu}\partial_\mu A_\nu = \frac{1}{2} \sigma^{\mu\nu}F_{\mu\nu}$, also hätte ich das richtige Ergebnis, wenn die $A_\mu\partial_\nu$Begriff verschwindet. Aber ich habe kein gutes Argument dafür, warum das der Fall sein sollte. Warum sollte das sein?