Derivation of Euler-Lagrange equation from principle of least action When deriving the Euler-Lagrange equation for a field $\phi$ the term 
$$
\int\textrm{d}x^{\mu}~\partial_{\mu}\left( \dfrac{\partial \mathcal{L}
}{\partial(\partial_{\mu}\phi)}\right)\delta\phi
$$
is assumed to vanish as it can be turned into a surface integral and $\delta\phi$ is assumed to vanish on the surface. What are the exact steps to turn this into a surface integral?
 A: 123hoedjevan gives you a wrong answer. The principle of least action states that the physical configuration of the system of fields realizes a minimum of the action with respect to compactly supported variations of the fields which, by the very definition of compactly supportedness, must then vanish on the boundary of the support itself. This in turn means that, in the formula
$$
\delta S=\int_{\Omega}d^{d}x\ \bigg[\bigg(\frac{\partial\mathcal{L}}{\partial x^{\mu}}-\partial_{\mu}\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\bigg)\,\delta\phi\bigg]+\int_{\Omega}d^{d}x\ \partial_{\mu}\bigg(\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\,\delta\phi\bigg)
$$
($d=\dim \Omega$), $\Omega$ must be compact. Now, as for the last term, we get
$$
\int_{\Omega}d^{d}x\ \partial_{\mu}\bigg(\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\,\delta\phi\bigg)=\int_{\partial\Omega}d^{d-1}x_{\mu}\ \bigg(\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\,\delta\phi\bigg)
$$
where I denoted by $d^{d-1}x_{\mu}$ the (oriented) volume element of the boundary $\partial\Omega$. The identity follows from Stoke's theorem, which (in one of its many forms) states that if you have a function $f$ defined on a compact set $\Omega$, then
$$
\int_{\Omega} d^{d}x\ \big(\partial_{\mu}\,f\big)=\int_{\partial\Omega} d^{d-1}x\ \big(f\ n_{\mu}\big)
$$
where $\partial\Omega$ is the boundary of $\Omega$ and $n_{\mu}$ are the components of the vector field normal to $\partial \Omega$ (notice that $n_{\mu}\,d^{d-1}x$ and my definition of $d^{d-1}x_{\mu}$ are the same thing). The proof of the theorem can be easily found on standard textbooks or on the internet. Going back to our integral, as $\delta\phi$ is by definition (i.e. as part of the hypotheses of the theorem) zero on the boundary of $\Omega$,
$$
\int_{\Omega}d^{d}x\ \partial_{\mu}\bigg(\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\,\delta\phi\bigg)=\int_{\partial\Omega}d^{d-1}x_{\mu}\ \bigg(\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\,\delta\phi\bigg)=\int_{\partial\Omega}d^{d-1}x_{\mu}\ 0=0
$$
hence, as $\delta S=0$,
$$
\delta S=\int_{\Omega}d^{d}x\ \bigg[\bigg(\frac{\partial\mathcal{L}}{\partial x^{\mu}}-\partial_{\mu}\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}\bigg)\,\delta\phi\bigg]=0
$$
and since this must hold for any compact $\Omega$ and any compactly supported $\delta\phi$,
$$
\frac{\partial\mathcal{L}}{\partial x^{\mu}}-\partial_{\mu}\frac{\partial \mathcal{L}}{\partial\, \partial_{\mu}\phi}=0
$$
The compactness of $\Omega$ (and in turn the compactly supportedness of the variation of the fields) can indeed be removed from the hypotheses of the theorem, as long as the action integral is well-defined on $\Omega$. On the other hand, the fact that the variation must vanish on the boundary of the domain of integration cannot be removed from the hypotheses. Hence the former is still valid for non-compact $\Omega$'s.
When proving Noether's theorem (which is different from proving the equivalence between minimization and the Euler-Lagrange equations, and ultimately depends on this very proof), one allows for variations that do not vanish on the boundary of the domain of integration; moreover, in its usual formulation, Noether's theorem allows for the coordinates of the domain of integration to be varied. It is in this context that the Noether current arises as a divergenceless field, and the Noether current is defined as
$$
j^{\mu}=-T^{\mu}_{\nu}\ \delta x^{\nu}+\frac{\partial\mathcal{L}}{\partial\,\partial_{\mu}\phi}\ \delta\phi
$$
where
$$
T^{\mu}_{\nu}=\frac{\partial\mathcal{L}}{\partial\,\partial_{\mu}\phi}\ \partial_{\nu}\phi-\mathcal{L}\ \delta^{\mu}_{\nu}
$$
is the canonical energy-momentum tensor.
A: The real answer is: it doesn't really. Or rather: we can still extract some physics out of it! Lets derive
$$
S[\phi] = \int d^d x \mathcal{L}[\phi,\partial\phi]\\
\delta S[\phi] = \int d^dx \delta \mathcal{L} = \int d^dx \left( \frac{\partial\mathcal{L}}{\partial \phi} \delta \phi + \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta (\partial_\mu\phi) \right)
$$
we first use that $\delta$ and $\partial$ commute and then Leibniz rule to get
$$
= \int d^dx \left( \frac{\partial\mathcal{L}}{\partial \phi} \delta \phi + \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \partial_\mu(\delta \phi) \right)\\
= \int d^dx \left( \frac{\partial\mathcal{L}}{\partial \phi} \delta \phi + \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi \right)  - \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi}\right) \delta \phi  \right)\\
= \int d^dx \left( \underbrace{\left( \frac{\partial\mathcal{L}}{\partial \phi} - \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi}\right) \right)}_\text{EoM = 0} \delta \phi  + \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi \right)   \right) = 0\\
$$
So we are left with 
$$
\delta S [\phi]  = \int_{\mathbb{R}^n} d^d x \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi   \right) = 0 \\
$$
There are 2 ways this can be $0$. 


*

*First we see that by stokes we get the GLOBAL statement
$$
\delta S [\phi]  = \int d \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi   \right) = \left.\frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi \right|_{\partial \mathbb{R}^n} =0 \\
$$
where we usually assume all fields vanish at $\partial \mathbb{R}^n$, i.e. at infinity. 

*Another way to get  $\delta S = 0$ (and in a sense a more physical and mathematically rigorous way) is to demand the LOCAL statement $ \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi   \right) = 0$. This is exactly the conserved current demand: $\partial_\mu j^\mu = 0$. So we find that 
$$
j^\mu =  \frac{\partial\mathcal{L}}{\partial \partial_\mu\phi} \delta \phi
$$
is the conserved current associated with variation $\delta\phi$. 
