Noether's theorem

Emmy Noether was an influential German mathematician known for her groundbreaking contributions to abstract algebra and theoretical physics.

Noether's (first)^[1] theorem states that every differentiable symmetry of the action of a physical system has a corresponding conservation law. The theorem was proven by German mathematician Emmy Noether in 1915 and published in 1918.^[2] The action of a physical system is the integral over time of a Lagrangian function (which may or may not be an integral over space of a Lagrangian density function), from which the system's behavior can be determined by the principle of least action.

Noether's theorem is used in theoretical physics and the calculus of variations. A generalization of the formulations on constants of motion in Lagrangian and Hamiltonian mechanics (developed in 1788 and 1833, respectively), it does not apply to systems that cannot be modeled with a Lagrangian alone (e.g. systems with a Rayleigh dissipation function). In particular, dissipative systems with continuous symmetries need not have a corresponding conservation law.

Basic illustrations and background

As an illustration, if a physical system behaves the same regardless of how it is oriented in space, its Lagrangian is rotationally symmetric: from this symmetry, Noether's theorem dictates that the angular momentum of the system be conserved, as a consequence of its laws of motion. The physical system itself need not be symmetric; a jagged asteroid tumbling in space conserves angular momentum despite its asymmetry. It is the laws of its motion that are symmetric.

As another example, if a physical process exhibits the same outcomes regardless of place or time, then its Lagrangian is symmetric under continuous translations in space and time: by Noether's theorem, these symmetries account for the conservation laws of linear momentum and energy within this system, respectively.

Noether's theorem is important, both because of the insight it gives into conservation laws, and also as a practical calculational tool. It allows investigators to determine the conserved quantities (invariants) from the observed symmetries of a physical system. Conversely, it allows researchers to consider whole classes of hypothetical Lagrangians with given invariants, to describe a physical system. As an illustration, suppose that a physical theory is proposed which conserves a quantity X. A researcher can calculate the types of Lagrangians that conserve X through a continuous symmetry. Due to Noether's theorem, the properties of these Lagrangians provide further criteria to understand the implications and judge the fitness of the new theory.

There are numerous versions of Noether's theorem, with varying degrees of generality. The original version only applied to ordinary differential equations (particles) and not partial differential equations (fields). The original versions also assume that the Lagrangian only depends upon the first derivative, while later versions generalize the theorem to Lagrangians depending on the n^th derivative. There are natural quantum counterparts of this theorem, expressed in the Ward–Takahashi identities. Generalizations of Noether's theorem to superspaces are also available.

Informal statement of the theorem

All fine technical points aside, Noether's theorem can be stated informally

If a system has a continuous symmetry property, then there are corresponding quantities whose values are conserved in time.^[3]

A more sophisticated version of the theorem involving fields states that:

To every differentiable symmetry generated by local actions, there corresponds a conserved current.

The word "symmetry" in the above statement refers more precisely to the covariance of the form that a physical law takes with respect to a one-dimensional Lie group of transformations satisfying certain technical criteria. The conservation law of a physical quantity is usually expressed as a continuity equation.

The formal proof of the theorem utilizes the condition of invariance to derive an expression for a current associated with a conserved physical quantity. In modern (since ca. 1980^[4]) terminology, the conserved quantity is called the Noether charge, while the flow carrying that charge is called the Noether current. The Noether current is defined up to a solenoidal (divergenceless) vector field.

In the context of gravitation, Felix Klein's statement of Noether's theorem for action I stipulates for the invariants:^[5]

If an integral I is invariant under a continuous group G_ρ with ρ parameters, then ρ linearly independent combinations of the Lagrangian expressions are divergences.

Historical context

A conservation law states that some quantity X in the mathematical description of a system's evolution remains constant throughout its motion — it is an invariant. Mathematically, the rate of change of X (its derivative with respect to time) vanishes,

$\frac{dX}{dt} = 0 ~.$

Such quantities are said to be conserved; they are often called constants of motion (although motion per se need not be involved, just evolution in time). For example, if the energy of a system is conserved, its energy is invariant at all times, which imposes a constraint on the system's motion and may help in solving for it. Aside from insights that such constants of motion give into the nature of a system, they are a useful calculational tool; for example, an approximate solution can be corrected by finding the nearest state that satisfies the suitable conservation laws.

The earliest constants of motion discovered were momentum and energy, which were proposed in the 17th century by René Descartes and Gottfried Leibniz on the basis of collision experiments, and refined by subsequent researchers. Isaac Newton was the first to enunciate the conservation of momentum in its modern form, and showed that it was a consequence of Newton's third law. According to general relativity, the conservation laws of linear momentum, energy and angular momentum are only exactly true globally when expressed in terms of the sum of the stress–energy tensor (non-gravitational stress–energy) and the Landau–Lifshitz stress–energy–momentum pseudotensor (gravitational stress–energy). The local conservation of non-gravitational linear momentum and energy in a free-falling reference frame is expressed by the vanishing of the covariant divergence of the stress–energy tensor. Another important conserved quantity, discovered in studies of the celestial mechanics of astronomical bodies, is the Laplace–Runge–Lenz vector.

In the late 18th and early 19th centuries, physicists developed more systematic methods for discovering invariants. A major advance came in 1788 with the development of Lagrangian mechanics, which is related to the principle of least action. In this approach, the state of the system can be described by any type of generalized coordinates q; the laws of motion need not be expressed in a Cartesian coordinate system, as was customary in Newtonian mechanics. The action is defined as the time integral I of a function known as the Lagrangian L

$I = \int L(\mathbf{q}, \dot{\mathbf{q}}, t) \, dt ~,$

where the dot over q signifies the rate of change of the coordinates q,

$\dot{\mathbf{q}} = \frac{d\mathbf{q}}{dt} ~.$

Hamilton's principle states that the physical path q(t)—the one actually taken by the system—is a path for which infinitesimal variations in that path cause no change in I, at least up to first order. This principle results in the Euler–Lagrange equations,

$\frac{d}{dt} \left( \frac{\partial L}{\partial \dot{\mathbf{q}}} \right) = \frac{\partial L}{\partial \mathbf{q}} ~.$

Thus, if one of the coordinates, say q_k, does not appear in the Lagrangian, the right-hand side of the equation is zero, and the left-hand side requires that

$\frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}_k} \right) = \frac{dp_k}{dt} = 0~,$

where the momentum

$p_k = \frac{\partial L}{\partial \dot{q}_k}$

is conserved throughout the motion (on the physical path).

Thus, the absence of the ignorable coordinate q_k from the Lagrangian implies that the Lagrangian is unaffected by changes or transformations of q_k; the Lagrangian is invariant, and is said to exhibit a symmetry under such transformations. This is the seed idea generalized in Noether's theorem.

Several alternative methods for finding conserved quantities were developed in the 19th century, especially by William Rowan Hamilton. For example, he developed a theory of canonical transformations which allowed changing coordinates so that some coordinates disappeared from the Lagrangian, as above, resulting in conserved canonical momenta. Another approach, and perhaps the most efficient for finding conserved quantities, is the Hamilton–Jacobi equation.

Mathematical expression

Simple form using perturbations

The essence of Noether's theorem is generalizing the ignorable coordinates outlined.

Imagine that the action I defined above is invariant under small perturbations (warpings) of the time variable t and the generalized coordinates q; in a notation commonly used in physics,

$t \rightarrow t^{\prime} = t + \delta t$

$\mathbf{q} \rightarrow \mathbf{q}^{\prime} = \mathbf{q} + \delta \mathbf{q} ~,$

where the perturbations δt and δq are both small, but variable. For generality, assume there are (say) N such symmetry transformations of the action, i.e. transformations leaving the action unchanged; labelled by an index r = 1, 2, 3, …, N.

Then the resultant perturbation can be written as a linear sum of the individual types of perturbations,

$\delta t = \sum_r \varepsilon_r T_r \!$

$\delta \mathbf{q} = \sum_r \varepsilon_r \mathbf{Q}_r ~,$

where ε_r are infinitesimal parameter coefficients corresponding to each:

• generator T_r of time evolution, and
• generator Q_r of the generalized coordinates.

For translations, Q_r is a constant with units of length; for rotations, it is an expression linear in the components of q, and the parameters make up an angle.

Using these definitions, Noether showed that the N quantities

$\left(\frac{\partial L}{\partial \dot{\mathbf{q}}} \cdot \dot{\mathbf{q}} - L \right) T_r - \frac{\partial L}{\partial \dot{\mathbf{q}}} \cdot \mathbf{Q}_r$

(which have the dimensions of [energy]·[time] + [momentum]·[length] = [action]) are conserved (constants of motion).

Examples

Time invariance

For illustration, consider a Lagrangian that does not depend on time, i.e., that is invariant (symmetric) under changes t → t + δt, without any change in the coordinates q. In this case, N = 1, T = 1 and Q = 0; the corresponding conserved quantity is the total energy H^[6]

$H = \frac{\partial L}{\partial \dot{\mathbf{q}}} \cdot \dot{\mathbf{q}} - L.$

Translational invariance

Consider a Lagrangian which does not depend on an ("ignorable", as above) coordinate q_k; so it is invariant (symmetric) under changes q_k → q_k + δq_k. In that case, N = 1, T = 0, and Q_k = 1; the conserved quantity is the corresponding momentum p_k^[7]

$p_k = \frac{\partial L}{\partial \dot{q_k}}.$

In special and general relativity, these apparently separate conservation laws are aspects of a single conservation law, that of the stress–energy tensor,^[8] that is derived in the next section.

Rotational invariance

The conservation of the angular momentum L = r × p is analogous to its linear momentum counterpart.^[9] It is assumed that the symmetry of the Lagrangian is rotational, i.e., that the Lagrangian does not depend on the absolute orientation of the physical system in space. For concreteness, assume that the Lagrangian does not change under small rotations of an angle δθ about an axis n; such a rotation transforms the Cartesian coordinates by the equation

$\mathbf{r} \rightarrow \mathbf{r} + \delta\theta \mathbf{n} \times \mathbf{r}.$

Since time is not being transformed, T=0. Taking δθ as the ε parameter and the Cartesian coordinates r as the generalized coordinates q, the corresponding Q variables are given by

$\mathbf{Q} = \mathbf{n} \times \mathbf{r}.$

Then Noether's theorem states that the following quantity is conserved,

$\frac{\partial L}{\partial \dot{\mathbf{q}}} \cdot \mathbf{Q}_{r} = \mathbf{p} \cdot \left( \mathbf{n} \times \mathbf{r} \right) = \mathbf{n} \cdot \left( \mathbf{r} \times \mathbf{p} \right) = \mathbf{n} \cdot \mathbf{L}.$

In other words, the component of the angular momentum L along the n axis is conserved.

If n is arbitrary, i.e., if the system is insensitive to any rotation, then every component of L is conserved; in short, angular momentum is conserved.

Field theory version

Although useful in its own right, the version of Noether's theorem just given is a special case of the general version derived in 1915. To give the flavor of the general theorem, a version of the Noether theorem for continuous fields in four-dimensional space–time is now given. Since field theory problems are more common in modern physics than mechanics problems, this field theory version is the most commonly used version (or most often implemented) of Noether's theorem.

Let there be a set of differentiable fields $\phi$ defined over all space and time; for example, the temperature $T(\mathbf{x}, t)$ would be representative of such a field, being a number defined at every place and time. The principle of least action can be applied to such fields, but the action is now an integral over space and time

$I = \int \mathcal{L} \left(\phi, \partial_\mu \phi, x^\mu \right) \, d^4 x$

(the theorem can actually be further generalized to the case where the Lagrangian depends on up to the n^th derivative using jet bundles)

A continuous transformation of the fields $\phi$ can be written infinitesimally as

$\phi \mapsto \phi + \epsilon \Psi$,

where $\Psi$ is in general a function that may depend on both $x^\mu$ and $\phi$. The condition for $\Psi$ to generate a physical symmetry is that the action $I$ is left invariant. This will certainly be true if the Lagrangian density $\mathcal{L}$ is left invariant, but it will also be true if the Lagrangian changes by a divergence,

$\mathcal{L} \mapsto \mathcal{L} + \epsilon \partial_\mu \Lambda^\mu$,

since the integral of a divergence becomes a boundary term according to the divergence theorem. A system described by a given action might have multiple independent symmetries of this type, indexed by $r = 1, 2, \ldots N$, so the most general symmetry transformation would be written as

$\phi \mapsto \phi + \epsilon_r \Psi_r$

with the consequence

$\mathcal{L} \mapsto \mathcal{L} + \epsilon_r \partial_\mu \Lambda^\mu_r$

For such systems, Noether's theorem states that there are $N$ conserved current densities

$j^\nu_r = \Lambda^\nu_r - \frac{\partial \mathcal{L}}{\partial \phi_{,\nu}} \cdot \Psi_r$

(where the dot product is understood to contract the field indices, not the $\nu$ index or $r$ index)

In such cases, the conservation law is expressed in a four-dimensional way

$\partial_{\nu} j^\nu = 0$

which expresses the idea that the amount of a conserved quantity within a sphere cannot change unless some of it flows out of the sphere. For example, electric charge is conserved; the amount of charge within a sphere cannot change unless some of the charge leaves the sphere.

For illustration, consider a physical system of fields that behaves the same under translations in time and space, as considered above; in other words, $L \left(\boldsymbol\phi, \partial_\mu{\boldsymbol\phi}, x^\mu \right)$ is constant in its third argument. In that case, N = 4, one for each dimension of space and time. An infinitesimal translation in space, $x^\mu \mapsto x^\mu + \epsilon_r \delta^\mu_r$ (with $\delta$ denoting the Kronecker delta), affects the fields as $\phi(x^\mu) \mapsto \phi(x^\mu - \epsilon_r \delta^\mu_r)$: that is, relabelling the coordinates is equivalent to leaving the coordinates in place while translating the field itself, which in turn is equivalent to transforming the field by replacing its value at each point $x^\mu$ with the value at the point $x^\mu - \epsilon X^\mu$ "behind" it which would be mapped onto $x^\mu$ by the infinitesimal displacement under consideration. Since this is infinitesimal, we may write this transformation as

$\Psi_r = -\delta^\mu_r \partial_\mu \phi$

The Lagrangian density transforms in the same way, $\mathcal{L}(x^\mu) \mapsto \mathcal{L}(x^\mu - \epsilon_r \delta^\mu_r)$, so

$\Lambda^\mu_r = -\delta^\mu_r \mathcal{L}$

and thus Noether's theorem corresponds to the conservation law for the stress–energy tensor T_μ^ν,^[8] where we have used $\mu$ in place of $r$. To wit, by using the expression given earlier, and collecting the four conserved currents (one for each $\mu$) into a tensor $T$, Noether's theorem gives

$T_\mu{}^\nu = -\delta^\nu_\mu \mathcal{L} + \delta^\sigma_\mu \partial_\sigma \phi \frac{\partial \mathcal{L}}{\partial \phi_{,\nu}} = \left(\frac{\partial \mathcal{L}}{\partial \phi_{,\nu}}\right) \cdot \phi_{,\mu} - \delta^\nu_\mu \mathcal{L}$

with

$T_\mu{}^\nu{}_{,\nu} = 0$

(note that we relabelled $\mu$ as $\sigma$ at an intermediate step to avoid conflict). (However, note that the $T$ obtained in this way may differ from the symmetric tensor used as the source term in general relativity; see Canonical stress–energy tensor.)

The conservation of electric charge, by contrast, can be derived by considering Ψ linear in the fields φ rather than in the derivatives.^[10] In quantum mechanics, the probability amplitude ψ(x) of finding a particle at a point x is a complex field φ, because it ascribes a complex number to every point in space and time. The probability amplitude itself is physically unmeasurable; only the probability p = |ψ|² can be inferred from a set of measurements. Therefore, the system is invariant under transformations of the ψ field and its complex conjugate field ψ^* that leave |ψ|² unchanged, such as

$\psi \rightarrow e^{i\theta} \psi \ ,\ \psi^{*} \rightarrow e^{-i\theta} \psi^{*}~,$

a complex rotation. In the limit when the phase θ becomes infinitesimally small, δθ, it may be taken as the parameter ε, while the Ψ are equal to iψ and −iψ*, respectively. A specific example is the Klein–Gordon equation, the relativistically correct version of the Schrödinger equation for spinless particles, which has the Lagrangian density

$L = \psi_{,\nu} \psi^{*}_{,\mu} \eta^{\nu \mu} + m^2 \psi \psi^{*}.$

In this case, Noether's theorem states that the conserved (∂⋅j = 0) current equals

$j^{\nu} = i \left( \frac{\partial \psi}{\partial x^{\mu}} \psi^{*} - \frac{\partial \psi^{*}}{\partial x^{\mu}} \psi \right) \eta^{\nu \mu}~,$

which, when multiplied by the charge on that species of particle, equals the electric current density due to that type of particle. This "gauge invariance" was first noted by Hermann Weyl, and is one of the prototype gauge symmetries of physics.

Derivations

One independent variable

Consider the simplest case, a system with one independent variable, time. Suppose the dependent variables q are such that the action integral

$I = \int_{t_1}^{t_2} L [\mathbf{q} [t], \dot{\mathbf{q}} [t], t] \, dt$

is invariant under brief infinitesimal variations in the dependent variables. In other words, they satisfy the Euler–Lagrange equations

$\frac{d}{dt} \frac{\partial L}{\partial \dot{\mathbf{q}}} [t] = \frac{\partial L}{\partial \mathbf{q}} [t].$

And suppose that the integral is invariant under a continuous symmetry. Mathematically such a symmetry is represented as a flow, φ, which acts on the variables as follows

$t \rightarrow t' = t + \varepsilon T \!$

$\mathbf{q} [t] \rightarrow \mathbf{q}' [t'] = \phi [\mathbf{q} [t], \varepsilon] = \phi [\mathbf{q} [t' - \varepsilon T], \varepsilon]$

where ε is a real variable indicating the amount of flow, and T is a real constant (which could be zero) indicating how much the flow shifts time.

$\dot{\mathbf{q}} [t] \rightarrow \dot{\mathbf{q}}' [t'] = \frac{d}{dt} \phi [\mathbf{q} [t], \varepsilon] = \frac{\partial \phi}{\partial \mathbf{q}} [\mathbf{q} [t' - \varepsilon T], \varepsilon] \dot{\mathbf{q}} [t' - \varepsilon T] .$

The action integral flows to

$\begin{align} I' [\varepsilon] & = \int_{t_1 + \varepsilon T}^{t_2 + \varepsilon T} L [\mathbf{q}'[t'], \dot{\mathbf{q}}' [t'], t'] \, dt' \\[6pt] & = \int_{t_1 + \varepsilon T}^{t_2 + \varepsilon T} L [\phi [\mathbf{q} [t' - \varepsilon T], \varepsilon], \frac{\partial \phi}{\partial \mathbf{q}} [\mathbf{q} [t' - \varepsilon T], \varepsilon] \dot{\mathbf{q}} [t' - \varepsilon T], t'] \, dt' \end{align}$

which may be regarded as a function of ε. Calculating the derivative at ε = 0 and using the symmetry, we get

$\begin{align} 0 & = \frac{d I'}{d \varepsilon} [0] = L [\mathbf{q} [t_2], \dot{\mathbf{q}} [t_2], t_2] T - L [\mathbf{q} [t_1], \dot{\mathbf{q}} [t_1], t_1] T \\[6pt] & {} + \int_{t_1}^{t_2} \frac{\partial L}{\partial \mathbf{q}} \left( - \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} T + \frac{\partial \phi}{\partial \varepsilon} \right) + \frac{\partial L}{\partial \dot{\mathbf{q}}} \left( - \frac{\partial^2 \phi}{(\partial \mathbf{q})^2} {\dot{\mathbf{q}}}^2 T + \frac{\partial^2 \phi}{\partial \varepsilon \partial \mathbf{q}} \dot{\mathbf{q}} - \frac{\partial \phi}{\partial \mathbf{q}} \ddot{\mathbf{q}} T \right) \, dt. \end{align}$

Notice that the Euler–Lagrange equations imply

$\begin{align} \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} T \right) & = \left( \frac{d}{dt} \frac{\partial L}{\partial \dot{\mathbf{q}}} \right) \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \left( \frac{d}{dt} \frac{\partial \phi}{\partial \mathbf{q}} \right) \dot{\mathbf{q}} T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \ddot{\mathbf{q}} \, T \\[6pt] & = \frac{\partial L}{\partial \mathbf{q}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \left( \frac{\partial^2 \phi}{(\partial \mathbf{q})^2} \dot{\mathbf{q}} \right) \dot{\mathbf{q}} T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \ddot{\mathbf{q}} \, T. \end{align}$

Substituting this into the previous equation, one gets

$\begin{align} 0 & = \frac{d I'}{d \varepsilon} [0] = L [\mathbf{q} [t_2], \dot{\mathbf{q}} [t_2], t_2] T - L [\mathbf{q} [t_1], \dot{\mathbf{q}} [t_1], t_1] T - \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} [t_2] T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} [t_1] T \\[6pt] & {} + \int_{t_1}^{t_2} \frac{\partial L}{\partial \mathbf{q}} \frac{\partial \phi}{\partial \varepsilon} + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial^2 \phi}{\partial \varepsilon \partial \mathbf{q}} \dot{\mathbf{q}} \, dt. \end{align}$

Again using the Euler–Lagrange equations we get

$\frac{d}{d t} \left( \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \varepsilon} \right) = \left( \frac{d}{d t} \frac{\partial L}{\partial \dot{\mathbf{q}}} \right) \frac{\partial \phi}{\partial \varepsilon} + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial^2 \phi}{\partial \varepsilon \partial \mathbf{q}} \dot{\mathbf{q}} = \frac{\partial L}{\partial \mathbf{q}} \frac{\partial \phi}{\partial \varepsilon} + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial^2 \phi}{\partial \varepsilon \partial \mathbf{q}} \dot{\mathbf{q}}.$

Substituting this into the previous equation, one gets

$\begin{align} 0 & = L [\mathbf{q} [t_2], \dot{\mathbf{q}} [t_2], t_2] T - L [\mathbf{q} [t_1], \dot{\mathbf{q}} [t_1], t_1] T - \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} [t_2] T + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} [t_1] T \\[6pt] & {} + \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \varepsilon} [t_2] - \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \varepsilon} [t_1]. \end{align}$

From which one can see that

$\left( \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \mathbf{q}} \dot{\mathbf{q}} - L \right) T - \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \varepsilon}$

is a constant of the motion, i.e., it is a conserved quantity. Since φ[q, 0] = q, we get $\frac{\partial \phi}{\partial \mathbf{q}} = 1$ and so the conserved quantity simplifies to

$\left( \frac{\partial L}{\partial \dot{\mathbf{q}}} \dot{\mathbf{q}} - L \right) T - \frac{\partial L}{\partial \dot{\mathbf{q}}} \frac{\partial \phi}{\partial \varepsilon}.$

To avoid excessive complication of the formulas, this derivation assumed that the flow does not change as time passes. The same result can be obtained in the more general case.

Field-theoretic derivation

Noether's theorem may also be derived for tensor fields φ^A where the index A ranges over the various components of the various tensor fields. These field quantities are functions defined over a four-dimensional space whose points are labeled by coordinates x^μ where the index μ ranges over time (μ = 0) and three spatial dimensions (μ = 1, 2, 3). These four coordinates are the independent variables; and the values of the fields at each event are the dependent variables. Under an infinitesimal transformation, the variation in the coordinates is written

$x^\mu \rightarrow \xi^\mu = x^\mu + \delta x^\mu \!$

whereas the transformation of the field variables is expressed as

$\phi^A \rightarrow \alpha^A (\xi^\mu) = \phi^A (x^\mu) + \delta \phi^A (x^\mu)\,.$

By this definition, the field variations δφ^A result from two factors: intrinsic changes in the field themselves and changes in coordinates, since the transformed field α^A depends on the transformed coordinates ξ^μ. To isolate the intrinsic changes, the field variation at a single point x^μ may be defined

$\alpha^A (x^\mu) = \phi^A (x^\mu) + \bar{\delta} \phi^A (x^\mu)\,.$

If the coordinates are changed, the boundary of the region of space–time over which the Lagrangian is being integrated also changes; the original boundary and its transformed version are denoted as Ω and Ω’, respectively.

Noether's theorem begins with the assumption that a specific transformation of the coordinates and field variables does not change the action, which is defined as the integral of the Lagrangian density over the given region of spacetime. Expressed mathematically, this assumption may be written as

$\int_{\Omega^{\prime}} L \left( \alpha^A, {\alpha^A}_{,\nu}, \xi^{\mu} \right) d^{4}\xi - \int_{\Omega} L \left( \phi^A, {\phi^A}_{,\nu}, x^{\mu} \right) d^{4}x = 0$

where the comma subscript indicates a partial derivative with respect to the coordinate(s) that follows the comma, e.g.

${\phi^A}_{,\sigma} = \frac{\partial \phi^A}{\partial x^{\sigma}}\,.$

Since ξ is a dummy variable of integration, and since the change in the boundary Ω is infinitesimal by assumption, the two integrals may be combined using the four-dimensional version of the divergence theorem into the following form

$\int_{\Omega} \left\{ \left[ L \left( \alpha^A, {\alpha^A}_{,\nu}, x^\mu \right) - L \left( \phi^A, {\phi^A}_{,\nu}, x^\mu \right) \right] + \frac{\partial}{\partial x^\sigma} \left[ L \left( \phi^A, {\phi^A}_{,\nu}, x^\mu \right) \delta x^\sigma \right] \right\} d^4 x = 0 \,.$

The difference in Lagrangians can be written to first-order in the infinitesimal variations as

$\left[ L \left( \alpha^A, {\alpha^A}_{,\nu}, x^\mu \right) - L \left( \phi^A, {\phi^A}_{,\nu}, x^\mu \right) \right] = \frac{\partial L}{\partial \phi^A} \bar{\delta} \phi^A + \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \bar{\delta} {\phi^A}_{,\sigma} \,.$

However, because the variations are defined at the same point as described above, the variation and the derivative can be done in reverse order; they commute

$\bar{\delta} {\phi^A}_{,\sigma} = \bar{\delta} \frac{\partial \phi^A}{\partial x^\sigma} = \frac{\partial}{\partial x^\sigma} ( \bar{\delta} \phi^A) \,.$

Using the Euler–Lagrange field equations

$\frac{\partial}{\partial x^\sigma} \left( \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \right) = \frac{\partial L}{\partial \phi^A}$

the difference in Lagrangians can be written neatly as

$\left[ L \left( \alpha^A, {\alpha^A}_{,\nu}, x^\mu \right) - L \left( \phi^A, {\phi^A}_{,\nu}, x^{\mu} \right) \right] = \frac{\partial}{\partial x^\sigma} \left( \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \right) \bar{\delta} \phi^A + \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \bar{\delta} {\phi^A}_{,\sigma} = \frac{\partial}{\partial x^\sigma} \left( \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \bar{\delta} \phi^A \right) \,.$

Thus, the change in the action can be written as

$\int_\Omega \frac{\partial}{\partial x^\sigma} \left\{ \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \bar{\delta} \phi^A + L \left( \phi^A, {\phi^A}_{,\nu}, x^\mu \right) \delta x^\sigma \right\} d^{4}x = 0 \,.$

Since this holds for any region Ω, the integrand must be zero

$\frac{\partial}{\partial x^\sigma} \left\{ \frac{\partial L}{\partial {\phi^A}_{,\sigma}} \bar{\delta} \phi^A + L \left( \phi^A, {\phi^A}_{,\nu}, x^\mu \right) \delta x^\sigma \right\} = 0 \,.$

For any combination of the various symmetry transformations, the perturbation can be written

$\delta x^{\mu} = \varepsilon X^\mu\!$

$\delta \phi^A = \varepsilon \Psi^A = \bar{\delta} \phi^A + \varepsilon \mathcal{L}_X \phi^A$

where $\mathcal{L}_X \phi^A$ is the Lie derivative of φ^A in the X^μ direction. When φ^A is a scalar or ${X^\mu}_{,\nu} = 0 \,$,

$\mathcal{L}_X \phi^A = \frac{\partial \phi^A}{\partial x^{\mu}} X^\mu\,.$

These equations imply that the field variation taken at one point equals

$\bar{\delta} \phi^A = \varepsilon \Psi^A - \varepsilon \mathcal{L}_X \phi^A\,.$

Differentiating the above divergence with respect to ε at ε=0 and changing the sign yields the conservation law

$\frac{\partial }{\partial x^\sigma} j^\sigma = 0$

where the conserved current equals

$j^\sigma = \left[\frac{\partial L}{\partial {\phi^A}_{,\sigma}} \mathcal{L}_X \phi^A - L \, X^\sigma\right] - \left(\frac{\partial L}{\partial {\phi^A}_{,\sigma}} \right) \Psi^A\,.$

Manifold/fiber bundle derivation

Suppose we have an n-dimensional oriented Riemannian manifold, M and a target manifold T. Let $\mathcal{C}$ be the configuration space of smooth functions from M to T. (More generally, we can have smooth sections of a fiber bundle over M.)

Examples of this M in physics include:

• In classical mechanics, in the Hamiltonian formulation, M is the one-dimensional manifold R, representing time and the target space is the cotangent bundle of space of generalized positions.
• In field theory, M is the spacetime manifold and the target space is the set of values the fields can take at any given point. For example, if there are m real-valued scalar fields, $\phi_1,\ldots,\phi_m$, then the target manifold is R^m. If the field is a real vector field, then the target manifold is isomorphic to R³.

Now suppose there is a functional

$\mathcal{S}:\mathcal{C}\rightarrow \mathbf{R},$

called the action. (Note that it takes values into R, rather than C; this is for physical reasons, and doesn't really matter for this proof.)

To get to the usual version of Noether's theorem, we need additional restrictions on the action. We assume $\mathcal{S}[\phi]$ is the integral over M of a function

$\mathcal{L}(\phi,\partial_\mu\phi,x)$

called the Lagrangian density, depending on φ, its derivative and the position. In other words, for φ in $\mathcal{C}$

$\mathcal{S}[\phi]\,=\,\int_M \mathcal{L}[\phi(x),\partial_\mu\phi(x),x] \mathrm{d}^nx.$

Suppose we are given boundary conditions, i.e., a specification of the value of φ at the boundary if M is compact, or some limit on φ as x approaches ∞. Then the subspace of $\mathcal{C}$ consisting of functions φ such that all functional derivatives of $\mathcal{S}$ at φ are zero, that is:

$\frac{\delta \mathcal{S}[\phi]}{\delta \phi(x)}\approx 0$

and that φ satisfies the given boundary conditions, is the subspace of on shell solutions. (See principle of stationary action)

Now, suppose we have an infinitesimal transformation on $\mathcal{C}$, generated by a functional derivation, Q such that

$Q \left[ \int_N \mathcal{L} \, \mathrm{d}^n x \right] \approx \int_{\partial N} f^\mu [\phi(x),\partial\phi,\partial\partial\phi,\ldots] \, \mathrm{d}s_\mu$

for all compact submanifolds N or in other words,

$Q[\mathcal{L}(x)]\approx\partial_\mu f^\mu(x)$

for all x, where we set

$\mathcal{L}(x)=\mathcal{L}[\phi(x), \partial_\mu \phi(x),x].\$

If this holds on shell and off shell, we say Q generates an off-shell symmetry. If this only holds on shell, we say Q generates an on-shell symmetry. Then, we say Q is a generator of a one parameter symmetry Lie group.

Now, for any N, because of the Euler–Lagrange theorem, on shell (and only on-shell), we have

$Q\left[\int_N \mathcal{L} \, \mathrm{d}^nx \right]$	$=\int_N \left[\frac{\partial\mathcal{L}}{\partial\phi}- \partial_\mu\frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\right]Q[\phi] \, \mathrm{d}^nx + \int_{\partial N} \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}Q[\phi] \, \mathrm{d}s_\mu$
	$\approx\int_{\partial N} f^\mu \, \mathrm{d}s_\mu .$

Since this is true for any N, we have

$\partial_\mu\left[\frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}Q[\phi]-f^\mu\right]\approx 0.$

But this is the continuity equation for the current $J^\mu\,\!$ defined by:^[11]

$J^\mu\,=\,\frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}Q[\phi]-f^\mu,$

which is called the Noether current associated with the symmetry. The continuity equation tells us that if we integrate this current over a space-like slice, we get a conserved quantity called the Noether charge (provided, of course, if M is noncompact, the currents fall off sufficiently fast at infinity).

Comments

Noether's theorem is an on shell theorem: it relies on use of the equations of motion—the classical path. It reflects the relation between the boundary conditions and the variational principle. Assuming no boundary terms in the action, Noether's theorem implies that

$\int_{\partial N} J^\mu \mathrm{d}s_\mu \approx 0~.$

The quantum analogs of Noether's theorem involving expectation values, e.g. ⟨∫d⁴x ∂·J⟩ = 0, probing off shell quantities as well are the Ward–Takahashi identities.

Generalization to Lie algebras

Suppose say we have two symmetry derivations Q₁ and Q₂. Then, [Q₁, Q₂] is also a symmetry derivation. Let's see this explicitly. Let's say

$Q_1[\mathcal{L}]\approx\partial_\mu f_1^\mu$

and

$Q_2[\mathcal{L}]\approx\partial_\mu f_2^\mu$

Then,

$[Q_1,Q_2][\mathcal{L}]=Q_1[Q_2[\mathcal{L}]]-Q_2[Q_1[\mathcal{L}]]\approx\partial_\mu f_{12}^\mu$

where f₁₂ = Q₁[f₂^μ] − Q₂[f₁^μ]. So,

$j_{12}^\mu=\left(\frac{\partial}{\partial (\partial_\mu\phi)}\mathcal{L}\right)(Q_1[Q_2[\phi]]-Q_2[Q_1[\phi]])-f_{12}^\mu.$

This shows we can extend Noether's theorem to larger Lie algebras in a natural way.

Generalization of the proof

This applies to any local symmetry derivation Q satisfying QS ≈ 0, and also to more general local functional differentiable actions, including ones where the Lagrangian depends on higher derivatives of the fields. Let ε be any arbitrary smooth function of the spacetime (or time) manifold such that the closure of its support is disjoint from the boundary. ε is a test function. Then, because of the variational principle (which does not apply to the boundary, by the way), the derivation distribution q generated by q[ε][Φ(x)] = ε(x)Q[Φ(x)] satisfies q[ε][S] ≈ 0 for every ε, or more compactly, q(x)[S] ≈ 0 for all x not on the boundary (but remember that q(x) is a shorthand for a derivation distribution, not a derivation parametrized by x in general). This is the generalization of Noether's theorem.

To see how the generalization is related to the version given above, assume that the action is the spacetime integral of a Lagrangian that only depends on φ and its first derivatives. Also, assume

$Q[\mathcal{L}]\approx\partial_\mu f^\mu$

Then,

$\begin{align} q[\varepsilon][\mathcal{S}] & = \int q[\varepsilon][\mathcal{L}] \, \mathrm{d}^n x \\ & = \int \left\{ \left(\frac{\partial}{\partial \phi}\mathcal{L}\right) \varepsilon Q[\phi]+ \left[\frac{\partial}{\partial (\partial_\mu \phi)}\mathcal{L}\right]\partial_\mu(\varepsilon Q[\phi]) \right\} \, \mathrm{d}^n x \\ & = \int \left\{ \varepsilon Q[\mathcal{L}] + \partial_{\mu}\varepsilon \left[\frac{\partial}{\partial \left( \partial_\mu \phi\right)} \mathcal{L} \right] Q[\phi] \right\} \, \mathrm{d}^n x \\ & \approx \int \varepsilon \partial_\mu \left\{f^\mu-\left[\frac{\partial}{\partial (\partial_\mu\phi)}\mathcal{L}\right]Q[\phi]\right\} \, \mathrm{d}^n x \end{align}$

for all ε.

More generally, if the Lagrangian depends on higher derivatives, then

$\partial_\mu\left[f^\mu-\left[\frac{\partial}{\partial (\partial_\mu\phi)} \mathcal{L} \right] Q[\phi] - 2\left[\frac{\partial}{\partial (\partial_\mu \partial_\nu \phi)}\mathcal{L}\right]\partial_\nu Q[\phi]+\partial_\nu\left[\left[\frac{\partial}{\partial (\partial_\mu \partial_\nu \phi)}\mathcal{L}\right] Q[\phi]\right]-\,\cdots\right]\approx 0.$

Examples

Example 1: Conservation of energy

Looking at the specific case of a Newtonian particle of mass m, coordinate x, moving under the influence of a potential V, coordinatized by time t. The action, S, is:

$\begin{align} \mathcal{S}[x] & = \int L[x(t),\dot{x}(t)] \, dt \\ & = \int \left(\frac{m}{2}\sum_{i=1}^3\dot{x}_i^2-V(x(t))\right) \, dt. \end{align}$

The first term in the brackets is the kinetic energy of the particle, whilst the second is its potential energy. Consider the generator of time translations Q = ∂/∂t. In other words, $Q[x(t)]=\dot{x}(t)$. Note that x has an explicit dependence on time, whilst V does not; consequently:

$Q[L]=m \sum_i\dot{x}_i\ddot{x}_i-\sum_i\frac{\partial V(x)}{\partial x_i}\dot{x}_i = \frac{d}{dt}\left[\frac{m}{2}\sum_i\dot{x}_i^2-V(x)\right]$

so we can set

$f=\frac{m}{2} \sum_i\dot{x}_i^2-V(x).$

Then,

$\begin{align} j & = \sum_{i=1}^3\frac{\partial L}{\partial \dot{x}_i}Q[x_i]-f \\ & = m \sum_i\dot{x}_i^2 -\left[\frac{m}{2}\sum_i\dot{x}_i^2 -V(x)\right] \\ & = \frac{m}{2}\sum_i\dot{x}_i^2+V(x). \end{align}$

The right hand side is the energy, and Noether's theorem states that $\dot{j}=0$ (i.e. the principle of conservation of energy is a consequence of invariance under time translations).

More generally, if the Lagrangian does not depend explicitly on time, the quantity

$\sum_{i=1}^3 \frac{\partial L}{\partial \dot{x}_i}\dot{x_i}-L$

(called the Hamiltonian) is conserved.

Example 2: Conservation of center of momentum

Still considering 1-dimensional time, let

$\begin{align} \mathcal{S}[\vec{x}] & = \int \mathcal{L}[\vec{x}(t),\dot{\vec{x}}(t)] \, \mathrm{d}t \\ & = \int \left [\sum^N_{\alpha=1} \frac{m_\alpha}{2}(\dot{\vec{x}}_\alpha)^2 -\sum_{\alpha<\beta} V_{\alpha\beta}(\vec{x}_\beta-\vec{x}_\alpha)\right] \, \mathrm{d}t \end{align}$

i.e. N Newtonian particles where the potential only depends pairwise upon the relative displacement.

For $\vec{Q}$, let's consider the generator of Galilean transformations (i.e. a change in the frame of reference). In other words,

$Q_i[x^j_\alpha(t)]=t \delta^j_i. \,$

Note that

$\begin{align} Q_i[\mathcal{L}] & = \sum_\alpha m_\alpha \dot{x}_\alpha^i-\sum_{\alpha<\beta}\partial_i V_{\alpha\beta}(\vec{x}_\beta-\vec{x}_\alpha)(t-t) \\ & = \sum_\alpha m_\alpha \dot{x}_\alpha^i. \end{align}$

This has the form of $\frac{\mathrm{d}}{\mathrm{d}t}\sum_\alpha m_\alpha x^i_\alpha$ so we can set

$\vec{f}=\sum_\alpha m_\alpha \vec{x}_\alpha.$

Then,

$\vec{j}=\sum_\alpha \left(\frac{\partial}{\partial \dot{\vec{x}}_\alpha}\mathcal{L}\right)\cdot\vec{Q}[\vec{x}_\alpha]-\vec{f}$

$=\sum_\alpha (m_\alpha \dot{\vec{x}}_\alpha t-m_\alpha \vec{x}_\alpha)$

$=\vec{P}t-M\vec{x}_{CM}$

where $\vec{P}$ is the total momentum, M is the total mass and $\vec{x}_{CM}$ is the center of mass. Noether's theorem states:

$\dot{\vec{j}} = 0 \Rightarrow {\vec{P}}-M \dot{\vec{x}}_{CM} = 0.$

Example 3: Conformal transformation

Both examples 1 and 2 are over a 1-dimensional manifold (time). An example involving spacetime is a conformal transformation of a massless real scalar field with a quartic potential in (3 + 1)-Minkowski spacetime.

$\mathcal{S}[\phi]\,$	$=\int \mathcal{L}[\phi (x),\partial_\mu \phi (x)] \, \mathrm{d}^4x$
	$=\int \left( \frac{1}{2}\partial^\mu \phi \partial_\mu \phi -\lambda \phi^4\right ) \, \mathrm{d}^4x$

For Q, consider the generator of a spacetime rescaling. In other words,

$Q[\phi(x)]=x^\mu\partial_\mu \phi(x)+\phi(x). \!$

The second term on the right hand side is due to the "conformal weight" of φ. Note that

$Q[\mathcal{L}]=\partial^\mu\phi\left(\partial_\mu\phi+x^\nu\partial_\mu\partial_\nu\phi+\partial_\mu\phi\right)-4\lambda\phi^3\left(x^\mu\partial_\mu\phi+\phi\right).$

This has the form of

$\partial_\mu\left[\frac{1}{2}x^\mu\partial^\nu\phi\partial_\nu\phi-\lambda x^\mu \phi^4 \right] = \partial_\mu\left(x^\mu\mathcal{L}\right)$

(where we have performed a change of dummy indices) so set

$f^\mu=x^\mu\mathcal{L}.\,$

Then,

$j^\mu=\left[\frac{\partial}{\partial (\partial_\mu\phi)}\mathcal{L}\right]Q[\phi]-f^\mu$

$=\partial^\mu\phi\left(x^\nu\partial_\nu\phi+\phi\right)-x^\mu\left(\frac{1}{2}\partial^\nu\phi\partial_\nu\phi-\lambda\phi^4\right).$

Noether's theorem states that $\partial_\mu j^\mu = 0 \!$ (as one may explicitly check by substituting the Euler–Lagrange equations into the left hand side).

Note that if one tries to find the Ward–Takahashi analog of this equation, one runs into a problem because of anomalies.

Applications

Application of Noether's theorem allows physicists to gain powerful insights into any general theory in physics, by just analyzing the various transformations that would make the form of the laws involved invariant. For example:

• the invariance of physical systems with respect to spatial translation (in other words, that the laws of physics do not vary with locations in space) gives the law of conservation of linear momentum;
• invariance with respect to rotation gives the law of conservation of angular momentum;
• invariance with respect to time translation gives the well-known law of conservation of energy

In quantum field theory, the analog to Noether's theorem, the Ward–Takahashi identity, yields further conservation laws, such as the conservation of electric charge from the invariance with respect to a change in the phase factor of the complex field of the charged particle and the associated gauge of the electric potential and vector potential.

The Noether charge is also used in calculating the entropy of stationary black holes.^[12]

See also

• Charge (physics)
• Gauge symmetry
• Gauge symmetry (mathematics)
• Invariant (physics)
• Goldstone boson
• Symmetry in physics

Notes

1. See also Noether's second theorem.
2. Noether E (1918). "Invariante Variationsprobleme". Nachr. D. König. Gesellsch. D. Wiss. Zu Göttingen, Math-phys. Klasse 1918: 235–257.
3. Thompson, W.J. (1994). Angular Momentum: an illustrated guide to rotational symmetries for physical systems 1. Wiley. p. 5. ISBN 0-471-55264-X.
4. The term "Noether charge" occurs in Seligman, Group theory and its applications in physics, 1980: Latin American School of Physics, Mexico City, American Institute of Physics, 1981. It comes enters wider use during the 1980s, e.g. by G. Takeda in: Errol Gotsman, Gerald Tauber (eds.) From SU(3) to Gravity: Festschrift in Honor of Yuval Ne'eman, 1985, p. 196.
5. Nina Byers (1998) "E. Noether's Discovery of the Deep Connection Between Symmetries and Conservation Laws." in Proceedings of a Symposium on the Heritage of Emmy Noether, held on 2–4 December 1996, at the Bar-Ilan University, Israel, Appendix B.
6. Lanczos 1970, pp. 401–3
7. Lanczos 1970, pp. 403–4
8. Goldstein 1980, pp. 592–3
9. Lanczos 1970, pp. 404–5
10. Goldstein 1980, pp. 593–4
11. Michael E. Peskin, Daniel V. Schroeder (1995). An Introduction to Quantum Field Theory. Basic Books. p. 18. ISBN 0-201-50397-2.
12. Vivek Iyer; Wald (1995). "A comparison of Noether charge and Euclidean methods for Computing the Entropy of Stationary Black Holes". Physical Review D 52 (8): 4430–9. arXiv:gr-qc/9503052. Bibcode:1995PhRvD..52.4430I. doi:10.1103/PhysRevD.52.4430.

References

• Goldstein, Herbert (1980). Classical Mechanics (2nd ed.). Reading, MA: Addison-Wesley. pp. 588–596. ISBN 0-201-02918-9.
• Kosmann-Schwarzbach, Yvette (2010). The Noether theorems: Invariance and conservation laws in the twentieth century. Sources and Studies in the History of Mathematics and Physical Sciences. Springer-Verlag. ISBN 978-0-387-87867-6
• Lanczos, C. (1970). The Variational Principles of Mechanics (4th ed.). New York: Dover Publications. pp. 401–5. ISBN 0-486-65067-7.
• Olver, Peter (1993). Applications of Lie groups to differential equations. Graduate Texts in Mathematics 107 (2nd ed.). Springer-Verlag. ISBN 0-387-95000-1.
• Sardanashvily, G. (2016). Noether's Theorems. Applications in Mechanics and Field Theory. Springer-Verlag. ISBN 978-94-6239-171-0.

External links

• Emmy Noether; Mort Tavel (translator) (1971). "Invariant Variation Problems". Transport Theory and Statistical Physics 1 (3): 186–207. arXiv:physics/0503066. Bibcode:1971TTSP....1..186N. doi:10.1080/00411457108231446. (Original in Gott. Nachr. 1918:235–257)
• Emmy Noether (1918). "Invariante Variationsprobleme" (in German).
• Emmy Noether and The Fabric of Reality (video) on YouTube
• Byers, Nina (1998). "E. Noether's Discovery of the Deep Connection Between Symmetries and Conservation Laws". arXiv:physics/9807044 [physics.hist-ph].
• John Baez (2002) "Noether's Theorem in a Nutshell."
• Hanca, J.; Tulejab, S.; Hancova, M. (2004). "Symmetries and conservation laws: Consequences of Noether's theorem". American Journal of Physics 72 (4): 428–35. Bibcode:2004AmJPh..72..428H. doi:10.1119/1.1591764.
• Merced Montesinos; Ernesto Flores (2006). "Symmetric energy–momentum tensor in Maxwell, Yang–Mills, and Proca theories obtained using only Noether's theorem" (PDF). Revista Mexicana de Física 52: 29–36. arXiv:hep-th/0602190. Bibcode:2006RMxF...52...29M.
• Vladimir Cuesta; Merced Montesinos; José David Vergara (2007). "Gauge invariance of the action principle for gauge systems with noncanonical symplectic structures". Physical Review D 76: 025025. Bibcode:2007PhRvD..76b5025C. doi:10.1103/PhysRevD.76.025025.
• Sardanashvily (2009). "Gauge conservation laws in a general setting. Superpotential". International Journal of Geometric Methods in Modern Physics 6 (06): 1047. arXiv:0906.1732. Bibcode:2009arXiv0906.1732S. doi:10.1142/S0219887809003862.
• Neuenschwander, Dwight E. (2010). Emmy Noether's Wonderful Theorem. Johns Hopkins University Press. ISBN 978-0-8018-9694-1.
• Noether's Theorem at MathPages.