I want to show that the Maxwell action $$S = -\frac{1}{4}\int d^4 x F_{\mu\nu} F^{\mu\nu}$$ is invariant under conformal transformations in $d=4$.
For this I considered the proof given in Zee's book on GR, section IX.9 page 621. However, I think he is making an error there. His proof goes as follows. For a conformal transformation $x'^\mu = x^\mu - \xi^\mu,$ the gauge boson transforms as the Lie derivative $\delta A_\mu = \mathcal{L}_\xi A_\mu$. From here he generalizes this to the field strength, i.e. $\delta F_{\mu\nu} = \mathcal{L}_\xi F_{\mu\nu}$. But then when computing the action on the Lagrangian density, he writes $$\delta(F_{\mu\nu}F^{\mu\nu}) = 2 F^{\mu\nu}\delta F_{\mu\nu}.\tag{IX.9.16}$$ This seems incorrect to me, because the Lie derivative transforms differently depending on whether the indices are upstairs or downstairs. I would expect to see $$\delta(F_{\mu\nu}F^{\mu\nu}) = F^{\mu\nu}\delta F_{\mu\nu} +F_{\mu\nu}\delta F^{\mu\nu}.$$
Now I want to use the fact that these variations are Lie derivatives found by
\begin{align} \delta F^{\mu\nu} &= -(\partial_\rho \xi^\mu) F^{\rho\nu} - (\partial_\rho \xi^\nu) F^{\mu\rho} + \xi^\alpha \partial_\alpha F^{\mu\nu}\\ \delta F_{\mu\nu} &= (\partial_\mu \xi^\rho) F_{\rho\nu} + (\partial_\nu \xi^\rho) F_{\mu\rho} + \xi^\alpha \partial_\alpha F_{\mu\nu}.\tag{IX.9.16} \end{align}
However, this results in
$$\delta(F_{\mu\nu}F^{\mu\nu}) = \partial_\alpha(\xi^\alpha F_{\mu\nu}F^{\mu\nu}) - (\partial \cdot \xi)F_{\mu\nu}F^{\mu\nu}.$$
This is decisively not what we should obtain.
How can I fix this proof?