Version 1:
An infinitesimal variation on the fields $\phi\mapsto\phi'$ is said to be a symmetry if $\delta \mathcal{L}:=\mathcal{L}(\phi',\partial\phi')-\mathcal{L}(\phi,\partial\phi)$ is a total derivative $\partial_\mu F^\mu$. If this is the case, let $\phi'(x)=\phi(x)+\delta\phi(x)$. Then $$0=\delta\mathcal{L}-\partial_\mu F^\mu=\frac{\partial\mathcal{L}}{\partial\phi}\delta\phi+\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\partial_\mu\delta\phi-\partial_\mu F^\mu=\delta\phi\left(\frac{\partial\mathcal{L}}{\partial\phi}-\partial_\mu\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\right)+\partial_\mu\left(\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\delta\phi-F^\mu\right).\tag{1}$$ Therefore, on-shell, the current $$j^\mu:=\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\delta\phi-F^\mu\tag{2}$$ is conserved.
Energy-momentum tensor: $\delta\phi(x)=\epsilon^\mu\partial_\mu\phi$ and $\delta\mathcal{L}=\epsilon^\mu\partial_\mu\mathcal{L}=\partial_\mu(\epsilon^\mu\mathcal{L})$, which yields $j^\mu=\epsilon^\nu T_\nu^\mu$, with $$T_\nu^\mu=\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\partial_\nu\phi-\delta^\mu_\nu\mathcal{L}.\tag{3}$$
Advantages:
- There is no mention to horizontal changes ($x\mapsto x'=x+\delta x$). In particular, the only data required is a vector field in the space of field configurations.
- It includes the possibility that the action is modified by a boundary term.
- The derivation is simple.
Disadvantages: It doesn't yield a method of computing $F^\mu$. I always get confused when I try to compute $F^\mu$ because what I end up doing is $\frac{\partial\mathcal{L}}{\partial\phi}\delta\phi+\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\partial_\mu\delta\phi.$
Version 2: An infinitesimal variation $x\mapsto x'=x+\delta x$ and $\phi\mapsto\phi'$ with $\phi'(x')=\phi(x)+\delta\phi(x)$ is a symmetry if $\delta S_\Omega(\phi):=S_{\Omega'}(\phi')-S_\Omega(\phi)=0$. After a computation one gets that for a general transformation (not necessarily a symmetry) $$\delta S_\Omega(\phi)=\int d^D x\left(-\partial_\mu\delta x^\nu T^\mu_\nu+\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\partial_\mu\delta\phi+\delta\phi\frac{\partial\mathcal{L}}{\partial\phi}\right).\tag{3}$$ Assuming $\delta x=0$ and that $\delta\phi$ vanishes on $\partial\Omega$, one gets to the Euler-Lagrange equations. On the other hand, if $\delta x^\mu=\omega^a X_a^\mu(x)$ and $\delta\phi(x)=\omega^a \mathcal{F}_a(\phi(x))$, one obtains $$\delta S_\Omega(\phi)=\int d^D x\omega^a\left(-\partial_\mu X_a^\nu T^\mu_\nu+\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\partial_\mu\mathcal{F}_a(\phi)+\mathcal{F}_a(\phi)\frac{\partial\mathcal{L}}{\partial\phi}\right)-\int d^D x\partial_\mu\omega^a j^\mu_a,\tag{4}$$ with $$j^\mu_a=T^\mu_\nu X^\nu_a-\frac{\partial\mathcal{L}}{\partial\partial_\mu\phi}\mathcal{F}_a(\phi).\tag{5}$$ Now, assume that the transformation is a symmetry when $\omega^a$ is constant. Then, by definition of symmetry the stuff in the first parenthesis vanishes (this is in fact what one needs to check to ensure this is a symmetry). Then, returning to the case $\omega^a$ may vary with $x$, we obtain $$\delta S_\omega(\phi)=-\int d^Dx\partial_\mu\omega^aj^\mu_a.\tag{6}$$ Then, by definition of the equations of motion (I have a comment on this below), this variation vanishes on shell. Choosing $\omega^a$ to vanish on $\partial\Omega$ one can integrate by parts and conclude that $j^\mu_a$ is conserved.
Energy-momentum tensor: $X^\mu_\nu=\delta^\mu_\nu$ and $\mathcal{F}_\nu(\phi(x))=0$.
Advantages:
- It shows explicitly how to check if a transformation is a symmetry.
- It yields an efficient method to compute the current (namely the coefficients of $\partial_\mu\omega^a$).
- It gives an explicit formula for the conserved current.
Disadvantages:
- The proof is long.
- It requires horizontal transformations
- At least in its present form it does not allow for variations of the action by boundary terms (I imagine this can be corrected easily).
Question: What is the relationship between these two formulations of Noether's theorem. I am particularly interested on why the first only requires the data of a vector field on the space of field configurations.
Side question: In version 2 there seems to be a loop hole. The vanishing of the variation of the action uses the on-shell condition. However, the Euler-Lagrange equations do not contemplate horizontal transformations. Then why can we guarantee that $\delta_\Omega S(\phi)=0$ on-shell?