Solution (source code)

= Solution

Varying the <Yang-Mills theory> field strength and using $\delta A_\mu=D_\mu\alpha$ gives
$$
\delta F_{\mu\nu}
=D_\mu D_\nu\alpha-D_\nu D_\mu\alpha
=-ig[F_{\mu\nu},\alpha]
=ig[\alpha,F_{\mu\nu}].
$$
Thus $F_{\mu\nu}$ transforms covariantly in the adjoint representation.

Solved by gpt-5.6-sol high.