Thursday, July 16, 2026

Hamiltonian - General Approach

Notes to a video lecture on UNIZOR.COM

General Approach
to Hamiltonian


In the previous introductory lecture we have examined the way to simplify the system of n Euler-Lagrange differential equations of second order into a system of 2n equations of the first order using coordinates q={qi} and momenta p={pi}.

For simple classical systems the Lagrangian Mechanics uses Euler-Lagrange equations
∂L/∂qi = d/dt ∂L/∂q̇i
where Lagrangian L=T−U,
T is the full kinetic energy of a system (independent of coordinates),
U is its potential energy (independent of velocities).

Instead, we constructed mathematically equivalent system of Hamiltonian equations
i = ∂(T+U)/∂pi
−ṗi = ∂(T+U)/∂qi
and introduced a new function called Hamiltonian
H(q,p) = T(p) + U(q)
where q=(q1,...,qn) are coordinates
and p=(p1,...,pn) are momenta.

This transformation increased the number of variables from n to 2n, but reduced the complexity of differential equations to the first order and made them symmetric and esthetically more appealing.

However, this construction relied on special properties of the system:
kinetic energy T depended only on momenta,
potential energy U depended only on coordinates.

In general mechanical systems, the Lagrangian L(q,q̇,t) may depend on coordinates and velocities in a more complicated way, kinetic and potential energies might not be separable and the Lagrangian might not be written as a sum of a function of p and a function of q.
Therefore, the construction H=T+U is not applicable to the general case.

The Hamiltonian approach is attractive for these simple mechanical systems, so it is natural to ask whether it can be generalized.
More specifically, given a Lagrangian L(q,q̇,t), we'll construct a function H(q,p,t) whose partial derivatives generate the differential equations similar to above and equivalent to Euler-Lagrange equations of motion.

To accomplish this, we will reformulate the expression H=T+U to avoid explicitly using energies and, instead, express H in terms of coordinates, momenta and a Lagrangian.

Recall the calculations of the Hamiltonian for simple mechanical system in the previous introductory lecture:
T = Σi½mi·q̇i² = Σi½pi²/mi

Partial derivative of T by pi produces
∂T/∂pi = pi/mi = q̇i
Therefore,
i = ∂T/∂pi
Our expression for H=T+U in that case can be represented as
H = 2T − (T−U) = 2T − L =
=
Σimi·q̇i² − L =
=
Σi(mi·q̇i)·q̇i − L =
=
Σipi·q̇i − L

The above expression
H(q,p,t) = Σipi·q̇i − L(q,q̇,t)
formally seems to correspond to our task of constructing a function H(q,p,t) from coordinates, momenta and a given Lagrangian L(q,q̇,t).
It does not explicitly rely on kinetic or potential energy. Since it is written only in terms of the Lagrangian, generalized coordinates, generalized velocities, and generalized momenta, it is a natural candidate for the Hamiltonian in the general case.

It's very important to understand the change of a viewpoint on relationship between arguments, which is a key idea of the Hamiltonian Mechanics.
While in the Lagrangian formulation the independent variables are q, and t, in the Hamiltonian approach the independent variables are q, p and t with generalized velocities become dependent on them:
i = q̇i(q,p,t)

Let's check if this expression corresponds to differential equations involving the Hamiltonian.
∂H(q,p,t)/∂pk =
[take into consideration that now qi and pi are independent variables, while generalized velocities i are functions of these variables: i=q̇i(q,p,t)]
= ∂/∂pk
[Σipi·q̇i − L] =
[use chain rule]
= q̇k +
Σipi·(∂q̇i/∂pk) −
Σi(∂L/∂q̇i)·(∂q̇i/∂pk) =
= q̇k+
Σi[pi−∂L/∂q̇i]·(∂q̇i/∂pk) =
[recall that pi=∂L/∂q̇i by definition]
= q̇k

as required.

∂H(q,p,t)/∂qk =
=
Σipi·(∂q̇i/∂qk) − ∂L/∂qk
Σi(∂L/∂q̇i)·(∂q̇i/∂qk) =
=
Σi[pi−∂L/∂q̇i]·(∂q̇i/∂qk) −
− ∂L/∂qk =
[recall that pi=∂L/∂q̇i by definition]
= − ∂L/∂qk =
[using Euler-Lagrange equation]
= − d/dt ∂L/∂q̇k =
[using the definition of generalized momentum]
= − d/dt pk = −ṗk

as required.

So, from the definitions of momenta as
pi = ∂L/∂q̇i
the Hamiltonian as
H(q,p,t) = Σipi·q̇i − L(q,q̇,t)
and Euler-Lagrange equation for L
∂L/∂qi = d/dt ∂L/∂q̇i
follows that both differential equations for H are satisfied:
∂H/∂pi = q̇i
∂H/∂qi = −ṗ

In reverse, assuming Hamiltonian equations are satisfied, we can derive Euler-Lagrang equations.
∂H(q,p,t)/∂qk =
= ∂/∂qk
[Σipi·q̇i − L] =
=
Σipi·(∂q̇i/∂qk) − ∂L/∂qk
Σi(∂L/∂q̇i)·(∂q̇i/∂qk) =
=
Σi[pi−∂L/∂q̇i]·(∂q̇i/∂qk) −
− ∂L/∂qk =
[recall that pi=∂L/∂q̇i by definition]
= − ∂L/∂qk

Since we assumed that Hamiltonian differential equations are true,
∂H(q,p,t)/∂qk = −ṗk
Therefore,
− ∂L/∂qk = −ṗk

From the definition
k = dpk/dt = d/dt ∂L/∂q̇k
follows:
∂L/∂qk = d/dt ∂L/∂q̇k
which is the Euler-Lagrange equation.
That completes the proof that a system of n second order Euler-Lagrange differential equations is equivalent to 2n first order Hamiltonian equations.

The construction of the Hamiltonian
H = Σipi·q̇i − L
from the Lagrangian L is called the Legendre transformation.

Hamiltonian Introduction

Notes to a video lecture on UNIZOR.COM

Introduction to Hamiltonian

We assume you have a pretty good understanding of Lagrangian Mechanics. If not, the previous chapters of this course Lagrangian and Noether Theorem provide a description of its basic principles.

The Hamiltonian Mechanics is built upon Lagrangian Mechanics. The differences can be summarized as follows.

In Lagrangian Mechanics our main dynamic variables were generalized coordinates q=(q1,...,qn) and their time derivatives - generalized velocities q̇=(q̇1,...,q̇n) (here we use Newtonian 'dot notation' to indicate a time derivative).
In Hamiltonian Mechanics it's the same generalized coordinates {qi} and, separately from coordinates, generalized momenta p=(p1,...,pn) instead of velocities.

Generalized momenta of a mechanical system with the Lagrangian L(q,q̇,t) are defined as a set of components
pi=∂L/∂q̇i
This definition was already introduced in the lecture Noether p=m·v const of a previous chapter of this course.

Using this definition of generalized momentum, the Euler-Lagrange differential equation of the second degree
d/dt ∂L/∂q̇i = ∂L/∂qi (i∈[1,n])
would look simpler
d/dt pi = ∂L/∂qi or, shorter,
i = ∂L/∂qi

In Lagrangian Mechanics we had n functions of time {qi(t)} (generalized coordinates) and a system of n Euler-Lagrange differential equations of the second order.
The number of unknowns was equal to the number of equations.

Instead, in this momentum-based approach, we have 2n functions of time {qi(t)} (generalized coordinates) and {pi(t)} (generalized momenta) with only n differential equations
(A) i = ∂L/∂qi
We need n more equations to obtain a system of 2n first-order differential equations equivalent to n Euler-Lagrange equations of the second-order.

Consider a simple case of a conservative system of one point mass m in three-dimensional Euclidean space with Cartesian coordinates (q1,q2,q3), velocities i and time-independent Lagrangian L equaled to a difference between kinetic T and potential U energies
L = T − U

Since kinetic energy T is independent of position qi, the same equation (A) above can be written as
i = ∂(−U)/∂qi
or
−ṗi = ∂U/∂qi

In addition to these n differential equations, we can construct n more using the classical definition of the vector of momentum pi=m·q̇i and the kinetic energy expressed in terms of momenta p as
T = Σi½mi·q̇i² = Σi½pi²/mi

Partial derivative of T by pi produces
∂T/∂pi = pi/mi = q̇i
Therefore,
(B) i = ∂T/∂pi
We have constructed n more differential equations to complete the system.

The differential equation for a time derivative of the generalized momentum, as we stated above, is
i = ∂L/∂qi = ∂(T−U)/∂qi
Since kinetic energy is independent of coordinates, we can exclude it
(C) −ṗi = ∂U/∂qi

Equations (B) and (C) constitute 2n differential equations of the first order with 2n unknowns - coordinates and momenta
i = ∂T/∂pi
−ṗi = ∂U/∂qi
The problem is, the first n equations depend on kinetic energy T, while the second group of n equations depends on potential energy U.

We would like a formulation in which both sets of equations are generated by a single function of the same variables (q,p), whose partial derivatives with respect to p give time derivatives of coordinates, that is velocities , and with respect to q give time derivatives of the generalized momenta .

Recall that in our simple case the kinetic energy T is independent of positions qi (∂T/∂qi=0) and potential energy U is independent of momenta pi (∂U/∂pi=0).
Based on this property, we can add to partial differentiation of the first equation the potential energy U and add to partial differentiation of the second equation the kinetic energy T.
i = ∂(T+U)/∂pi
−ṗi = ∂(T+U)/∂qi

Let's introduce a new function called Hamiltonian
H(q,p) = T(p) + U(q)
where q=(q1,...,qn) are coordinates
and p=(p1,...,pn) are momenta.
With this notation our system of 2n differential equations of the first degree looks quite symmetrical (some might say "beautiful")
i = ∂H/∂pi
−ṗi = ∂H/∂qi

In this case of a simple mechanical system the Hamiltonian H=T+U represents a total (kinetic and potential) energy of our mechanical system, which makes our system of equation more related to real physical characteristics of a system.

In more general systems, however, the definition of the Hamiltonian is broader and not necessarily coincides with a total energy.

Looking ahead, let us state that the symmetric form of Hamilton's equations makes Hamiltonian Mechanics especially suitable for advanced topics such as canonical transformations, statistical mechanics and quantum mechanics.

Friday, July 3, 2026

Noether's Theorem Conservation: UNIZOR.COM -> Physics+ 4 All -> Lagrangian -> Noether's Theorem -> Conservation

Notes to a video lecture on UNIZOR.COM

Noether Theorem and
Conservation Laws



Background

Motion of a mechanical system is represented by a curve in extended configuration space with coordinates y={yi(x)} - a set of time-space coordinates describing a motion curve parameterized by independent variable x∈[a,b] for i∈[0,n].

These coordinates have physical meaning:
y0 is time t;
{yi} are a set of generalized coordinates {qi} for i∈[1,n];

The action functional
Φ[t,q] =
t(b)
t(a)
L(t,q,q')dt

was re-parameterized as an integral by independent parameter x∈[a,b]
Φ[t,q] =
b
a
L(t,q,qx/tx)·tx·dx

and expressed as
Φ[y] =
b
a
𝓛(y,yx)·dx

where
y(x)={yi(x)} (i∈[0,n]) signifies a set of all time-space coordinates parameterized by x∈[a,b], that is a trajectory in extended configuration space, with y0(x)=t(x), and yi(x)=qi(x) for i0 and
yx(x)={yix(x)} (i∈[0,n]) signifies a set of all derivatives of time-space coordinates by parameter x with y0x(x)=tx(x), and yix(x)=qix(x) for i0
and a new function 𝓛() is defined for i∈[0,n] as
𝓛(y,yx) = 𝓛({yi},{yix}) =
= L(t,q,qx/tx)·tx =
= L(t,
{qi},{qix/tx})·tx

The conclusion of the previous lecture:
d/dx Σi𝓛yix·ζi = 0 for i∈[0,n] where
ζ={ζi}={dyi(ε)/dε|ε=0} is a set of generators for each time-space coordinate
and where
x in a subscript indicates a derivative of a corresponding function by parameter x:
tx=dt/dx and
{qix}={dqi/dx} for i∈[1,n]

Linear Momentum Conservation

Let's choose a single kth space coordinate (k∈[1,n]) and consider the following ε-transformation of coordinates:
t(ε) = t,
qk(ε) = qk + ε,
qi(ε) = qi for all i≠k.
This represents a uniform movement along the kth space coordinate qk.

The corresponding generators of this transformation are
ζk = dyk/dε|ε=0 = 1 and
ζi = 0 for i≠k.

The main result of the Noether's theorem was:
as long as the action functional is invariant under the transformation,
an expression
J = Σi𝓛yix·ζi
is a constant of motion along a trajectory and its x-derivative equals to zero for all x∈[a,b]:
dJ/dx = 0

Considering all ζi=0 for i∈[0,n] except ζk=1,
J = 𝓛ykx
and, therefore,
dJ/dx =
= d/dx
𝓛ykx({yi},{yix}) = 0
or
d/dx ∂/∂ykx 𝓛({yi},{yix}) = 0

To see what follows from this equation, let's return to the original Lagrangian
L(t,{qi},{qi'}),
where apostrophe at q indicates a time-derivative, using yk=qk equivalence for k≠0 and taking into account our definition of function 𝓛:
𝓛(y,yx) = L(t,{qi},{qix/tx})·tx
Notice that we represented qi'=dqi/dt (a generalized velocity) as (dqi/dx)/(dt/dx)=qix/tx.

Now
∂/∂ykx 𝓛({yi},{yix}) =
= ∂/∂qkx
[L(t,{qi},{qix/tx})·tx]

As we indicated above, an expression qix/tx is a generalized velocity along ith space coordinate because
qix/tx = (dqi/dx)/(dt/dx) =
= dqi/dt = qi' = vi


Therefore, using the chain rule, we can write the expression above as
∂/∂qkx [L(t,{qi},{qix/tx})·tx] =
= ∂/∂qkx
[L(t,{qi},{vi})·tx] =
=
[∂/∂vk L(t,q,v)]·[∂vk/∂qkx]·tx =
[recall, vk=qkx/tx]
=
[∂/∂vk L(t,q,v)]·(1/tx)·tx =
= ∂/∂vk L(t,q,v) =
= ∂/∂qk' L(t,q,q')


Since pk=∂/∂qk' L(t,q,q') is a definition of generalized momentum, we conclude that under the ε-transformation of a single space coordinate qk described above that leaves the action functional invariant or, in other words, possesses translational symmetry the generalized momentum is conserved.

The conservation of momentum pk could be derived directly from the Euler–Lagrange equations once we know that the Lagrangian is independent of a coordinate qk. Indeed, from the Euler-Lagrange equation
d/dt ∂L/∂qk' = ∂L/qk
follows that, if the right-hand side is zero (independence of Lagrangian L of coordinate qk), the left-hand side is zero as well, which means that generalized momentum
pk = ∂L/∂qk'
is constant (conserved).

Noether's theorem is important because it reveals that momentum conservation is a consequence of a continuous symmetry of the action and extends this principle to every continuous symmetry. Momentum conservation is therefore not an isolated fact but one example of a universal connection between symmetry and conservation laws.

CONCLUSION

Assuming the translation
qk ⟶ qk + ε
leaves action functional invariant,
the Noether conserved quantity would be
J = ∂𝓛/∂ykx = ∂L/∂qk' = pk
which is a generalized momentum along kth generalized coordinate. Therefore,
dpk/dt = 0
which is the law of conservation of linear momentum.


Angular Momentum Conservation

Consider a rigid body rotating about a fixed axis with only two extended generalized coordinates describing its motion
y0 = t is time,
y1 = θ is an angle of rotation.

The ε-transformation (rotation) we would like to consider is the uniform rotation that can be expressed as
t(ε) = t,
θ(ε) = θ + ε.

Since the angle θ is simply a generalized coordinate, this situation is the same as in the previous one
t(ε) = t,
qk(ε) = qk + ε
with n=k=1,
generalized space coordinate being the angle of rotation q1=y1
and the derivation of Noether's conserved quantity is identical to the derivation for linear momentum.

Therefore, everything we derived for linear momentum in the above case is valid for angular momentum

Assuming the rotation
θ(=y1) ⟶ θ + ε (=y1 + ε)
leaves action functional invariant,
the Noether conserved quantity would be
J = ∂𝓛/∂θx = ∂L/∂θ' = ℒ
which is a generalized angular momentum.
Therefore, the angular momentum is conserved
dℒ/dt = 0
which is the law of conservation of angular momentum.


Energy Conservation

The third important application of Noether's theorem is the law of conservation of energy, which follows from the invariance of the action under translations of time.

Let's choose a uniform translation of the time coordinate that does not affect any space coordinates:
y0(ε) = y0 + ε,
which is time transformation
t(ε) = t + ε
and
yi(ε) = yi for all i∈[1,n],
which means that all generalized coordinates remain unchanged
qi(ε) = qi for all i∈[1,n].

For this kind of transformation the corresponding generators are
ζ0 = dy0/dε|ε=0 = 1 and
ζi = 0 for i≠0.

According to Noether's theorem, the conserved quantity is
J = Σi𝓛yix·ζi = 𝓛y0x·ζ0 =
= ∂
𝓛/∂tx = ∂/∂tx[L(t,q,qx/tx)·tx]
which is a constant of motion along a trajectory.

Let's perform all the required computations.
J = ∂/∂tx[L(t,q,qx/tx)·tx] =
[recall, vi=qix/tx=qi' - time derivative of a generalized coordinate]
= ∂/∂tx
[L(t,{qi},{qix/tx})·tx] =
[using a formula of a derivative of a product of two functions]
=
[∂/∂txL(t,{qi},{qix/tx})]·tx + L=
[apply the chain rule, taking into account that the dependence on tx enters only through the generalized velocities {qix/tx} = {vi} and using vi as a placeholder for qix/tx to shorter the notation]
=
[Σi(∂L/∂vi)·(d(qix/tx)/dtx)]·tx + L =
=
[Σi(∂L/∂vi)·(−qix/t²x)]·tx + L =
[substitute qix/tx²=vi·tx/tx²=vi/tx]
= −
Σi(∂L/∂vi)·vi + L =
[recall, ∂L/∂vi=∂L/∂qi' is s generalized momentum pi]
= −
Σi(pi·vi) + L

The final formula for a conserved quantity J is:
J = −Σi(pi·vi) + L

In all conservative mechanical systems considered in this course, the Lagrangian L has the form
L = T − U
where the kinetic energy T is a quadratic homogeneous function of generalized velocities and U is potential energy of a system.

Recall that
pi = ∂L/∂vi = ∂(T−U)/∂vi
Since potential energy U does not depend on generalized velocities,
pi = ∂T/∂vi

In classical mechanics the kinetic energy is a homogeneous quadratic (that is, of degree 2) function of the generalized velocities
T = Σi,jAijvivj

For any given quadratic homogeneous function
T = Σi,jAijvivj
the sum Σi[∂T/∂vi]·vi is equal to 2T as follows from the Euler theorem about homogeneous functions.
Here is a simple and elegant proof.
In case of quadratic homogeneous function
T(v1,...,vn) = Σi,jAijvivj
T(λ·v1,...,λ·vn) = λ²·T(v1,...,vn)
where λ - any real number.
Let's differentiate both sides by λ applying the chain rule for the left side
Σi[∂T/∂(λ·vi)]·vi = 2λ·T
Set λ=1 that results in λ·vi=vi, and the result is
Σi[∂T/∂vi]·vi = 2·T

Since generalized momentum is defined by
pi = ∂T/∂vi
and for ordinary mechanical systems L=T−U while U does not depend on the velocities, ∂L/∂vi = ∂T/∂vi
therefore,
Σi(pi·vi) = 2T
J = −2T + (T−U) = −(T+U)
which is a negative total energy of the system, whose conservation is equivalent to conservation of the total system's energy itself.

Therefore, the total energy is conserved under a time transformation that preserves the action functional, as described above.

Sunday, June 21, 2026

Noether's Theorem Derivation: UNIZOR.COM -> Physics+ 4 All -> Lagrangian -> Noether's Theorem -> Derivation

Notes to a video lecture on UNIZOR.COM

Derivation of Noether Theorem


Background

The previous lectures of the Noether Theorem part of the course Physics+ 4 All have introduced the concepts of symmetry, parameterized group of continuous transformations and a concept of an extended configuration space that combines time and generalized coordinates into one set of coordinates.

We suggested that the symmetries relevant to the laws of motion are transformations of extended time-space coordinates that leave the action functional invariant.

This lecture is about mathematical derivation of certain conservation laws as logical consequences from the symmetries of transformations.


Summary of Assumptions

(1) Let us consider an extended configuration space of a mechanical system with coordinates {t,q}, where t is time and q is a set of generalized coordinates qi(t) (i∈[1,n]).

(2) This system is described by its Lagrangian L(t,q(t),q'(t)) where q'(t) is a set of generalized velocities {qi'(t)} (time derivatives of generalized coordinates).

(3) The trajectory of the movement of this system, a curve in the extended configuration space, is described by parameterized functions
t(x) and
q(x)={qi(x)}, i∈[1,n]
where x is an abstract parameter changing from real value x=a to x=b with {t(a),q(a)} being the start and {t(b),q(b)} being the finish point of the movement.

(4) Given a group of continuous transformations of the extended configuration space parameterized by ε
tt(ε)=T(ε,t,q)
qq(ε)=Q(ε,t,q)
with ε=0 causing a transformation to be the identity transformation, that is t(0)=t and q(0)=q.
We assume, the transformation functions T(ε,t,q) and Q(ε,t,q) are sufficiently differentiable.

(5) This transformation of points (t,q)⟶(t(ε),q(ε))
induces the transformation of every trajectory
{t(x),q(x)}⟶{t(ε)(x),q(ε)(x)}
where x∈[a,b].
We further assume that the transformed trajectory
{t(ε)(x),q(ε)(x)}
belongs to the same class of physically admissible trajectories of a mechanical system defined by its properties and the laws of physics.

(6) Let's assume that the action functional of the movement of this mechanical system
Φ[t,q] =
t(b)
t(a)
L(t,q,q')dt

is invariant under this induced transformation of trajectories as parameter of transformation ε is infinitesimally changing from zero.
This assumption can be expressed as
(d/dε)Φ[t(ε),q(ε)]|ε=0 = 0
and we assume sufficient differentiability of the action functional by parameter ε.


Derivation of Noether Theorem

The problem with the above representation of the action functional is that not only an expression under an integral is transformed, but the limits of integration t(a) and t(b) change as well, which significantly complicates the analysis of the behavior of the action functional under ε-transformation of time and generalized coordinates.
The road to simplification is the parameterized representation of a trajectory {t(x),q(x)}.
Using this, we can replace
(a) dt = (dt/dx)·dx
(b) q' = dq/dt = (dq/dx)/(dt/dx)
(c) integration by t on [t(a),t(b)] can be now replaced with an integration by x - the parameter changing on a fixed segment [a,b].

Let's rewrite the action functional as the integral by x using abbreviations
dt/dx=tx (so, dt=tx·dx) and
dq/dx={dqi/dx}={qix}=qx
for brevity
Φ[t,q] =
b
a
L(t,q,qx/tx)·tx·dx

At this point we would like to bring some time-space uniformity.
Since both time t and generalized space coordinates q={qi} (i∈[1,n]) are all functions of one parameter x (x∈[a,b]) and all have equal standing as coordinates in an extended time-space configuration space, it makes sense to use a single letter
y={yi} (i∈[0,n]) with
y0=t and
yi=qi for all i∈[1,n]).

Also, we replace derivatives
dt/dx=tx with dy0/dx=y0x
and
dq/dx={dqi/dx}={qix}=qx
for i≠0 with
dy/dx={dyi/dx}={yix}=yx.

So, a set of all derivatives {tx,qix} can be written as
dy/dx={dyi/dx}={yix}=yx
where i∈[0,n].

Now we can simplify the formula for action functional by replacing the Lagrangian under the integration with a function that treats all time-space coordinates equally:
Φ[y] =
b
a
𝓛(y,yx)·dx

where
y(x)={yi(x)} (i∈[0,n]) signifies a set of all time-space coordinates parameterized by x∈[a,b], that is a trajectory in extended configuration space, with y0(x)=t(x), and yi(x)=qi(x) for i0 and
yx(x)={yix(x)} (i∈[0,n]) signifies a set of all derivatives of time-space coordinates by parameter x with y0x(x)=tx(x), and yix(x)=qix(x) for i0
and a new function 𝓛() is defined for i∈[0,n] as
𝓛(y,yx) = 𝓛({yi},{yix}) =
= L(t,q,qx/tx)·tx =
= L(t,
{qi},{qix/tx})·tx

This representation of the same action functional is simpler because a new function 𝓛() under the integral symmetrically depends on n+1 functions {yi(x)} (i∈[0,n]) that encompass t(x) and all {qi(x)} (i∈[1,n]) functions of parameter x and n+1 derivatives of these functions by x, and the parameter x is not a subject of ε-transformation.

In addition, the limits of integration by x are from a to b which are constant and not affected by the ε-transformation.

The latter form of function 𝓛() under an integral allows to express the assumption about the invariance of the action functional under ε-transformations
y⟶y(ε)
which means t⟶t(ε), q⟶q(ε)),
as
(d/dε)Φ[y(ε)]|ε=0 = 0
in a symmetrical way relative to all time-space coordinates in extended configuration space and use the known apparatus of Calculus to perform all the required operations.

Since our transformations are continuous, ε-transformations with infinitesimal ε are infinitesimal, that is the increments in coordinates
Δy(ε)=y(ε)−y
are infinitesimal as well.
At the same time, the ε-derivatives of the changing coordinates characterize the speed of their change by a transformation.
Expressions
ζ={ζi}={dyi(ε)/dε|ε=0}
are called generators of the transformation.
We will use them below.

Also,
d/dε[dy(ε)/dx]|ε=0 =
= d/dx
[dy(ε)/dε]|ε=0 = dζ/dx

Let's apply our assumption about invariance of the action functional under ε-transformation and equate the ε-derivative of the action functional at ε=0 to zero and do the calculations.
We'll abbreviate derivatives with subscriptors for brevity (like ζx for dζ/dx) and omit the |ε=0 to shorten the formulas:
0 = (d/dε)Φ[y(ε)] =
= (d/dε)
b
a
𝓛(y(ε),yx(ε))·dx =
where y and yx are group variables representing all time-space coordinates in extended configuration space {yi} and {yix} with i∈[0,n].
We can change the order of differentiation by ε and integration by x because we assumed that the function under the integral is sufficiently smooth, so the convergence theorem holds.

=
b
a
(d/dε) 𝓛(y(ε),yx(ε))·dx =

use the rules for differentiation of multi-variable functions, subscriptions to indicate the corresponding derivative, definition ζ=dy(ε)/dε|ε=0 and the rule for interchanging the differentiation by two different variables
dyx(ε)/dε = d/dε[dy(ε)/dx] =
= d/dx
[dy(ε)/dε] = dζ/dx = ζx
=
b
a
(𝓛y·dy(ε)/dε+𝓛yx·dyx(ε)/dε)·dx
=
b
a
(𝓛y·ζ+𝓛yx·ζx)·dx

where group item 𝓛y·ζ represents
a sum Σi𝓛yi·ζi with i∈[0,n]
which in expanded form is
Σi𝓛/∂yi·[dyi(ε)/dε|ε=0]
and group item 𝓛yx·ζx represents
a sum Σi(𝓛yix·ζix) with i∈[0,n]
which expands analogously with a subscript x indicating a derivative by x

We have derived with a fundamental identity


b
a
[Σ(i𝓛yi·ζi)+Σi(𝓛yix·ζix)]·dx = 0
where summation by i is for i∈[0,n].

Recall the following rules for integration by parts.

(d/dx)[𝓛yx·ζ] =
= (d
𝓛yx/dx)·ζ + 𝓛yx·ζx

b
a
(d/dx)[𝓛yx·ζ]·dx =

=
b
a
(d𝓛yx/dx)·ζ·dx +
b
a
𝓛yx·ζx·dx


[𝓛yx·ζ]|[a,b] =
=
b
a
(d𝓛yx/dx)·ζ·dx +
b
a
𝓛yx·ζx·dx


[𝓛yx·ζ]|[a,b]
b
a
(d𝓛yx/dx)·ζ·dx =

=
b
a
𝓛yx·ζx·dx


The above expression for
b
a
𝓛yx·ζx·dx

can be substituted into the fundamental identity presented above
b
a
(𝓛y·ζ+𝓛yx·ζx)·dx = 0

where group item
𝓛y·ζ represents Σi𝓛yi·ζi

getting
0 = [𝓛yx·ζ]|[a,b] +
+
b
a
[𝓛y − d/dx(𝓛yx)]·ζ·dx


Since the ε-transformation leaves the action functional invariant, it preserves the set of stationary trajectories of the action mapping one stationary trajectory onto itself (shift along a trajectory) or to another one (jump to another trajectory).
Since all physical trajectories are stationary for the action functional and, therefore, are the solutions of the Euler–Lagrange equations, the transformed trajectory, being stationary as well, also satisfies the Euler–Lagrange equation:
𝓛y = 𝓛/∂y = d/dx(∂𝓛/∂yx) =
= d
𝓛yx/dx
which nullifies an integral in the last identity.

Therefore,
[𝓛yx·ζ]|[a,b] = 0 or
𝓛yx(b)·ζ(b) − 𝓛yx(a)·ζ(a) = 0
where group parameter
𝓛yx·ζ represents Σi𝓛yix·ζi with the sum by i∈[0,n].

This means that the value of 𝓛yx·ζ is the same at x=a and x=b.
This holds for any subinterval [a,b] of any physically admissable trajectory.
Since the endpoints can be chosen arbitrarily along the trajectory, the above expression is constant along the trajectory and
d/dx[𝓛yx·ζ] = 0

Noether Theorem

Every continuous symmetry of the action functional (every ε-transformation of the extended configuration space that leaves the action functional invariant) is associated with a conserved quantity along physical trajectories.
In extended configuration space, the conserved quantity is
J = 𝓛yx·ζ
or, in expanded by coordinates format,
J = Σi𝓛yix·ζi with i∈[0,n]
Recall that
y={yi} is a set of time-space coordinates with i∈[0,n];
y0 is time t;
{yi} are a set of generalized coordinates {qi} with i∈[1,n];
ζ={ζi}={dyi(ε)/dε|ε=0} is a set of generators for each time-space coordinate;
𝓛(y,yx) = L(t,tx,q,qx/tx)·tx
where x in a subscript indicates a derivative of a corresponding function by parameter x:
tx=dt/dx and
{qix}={dqi/dx}

We have just proven that
J is a constant of motion and, therefore,
dJ/dx = 0

Sunday, April 12, 2026

Noether's Theorem Math 1: UNIZOR.COM -> Physics+ 4 All -> Lagrangian -> Noether's Theorem Math 1: ε−Transformation

Notes to a video lecture on UNIZOR.COM

Noether Math 1:
ε−Transformation


In 1918 Emmy Noether rigorously proved that the conservation laws are consequences of certain continuous symmetries of time, space and other symmetries appearing in classical mechanics, field theory, quantum mechanics, general relativity and other parts of theoretical physics.

In a very simplified form she proved that, if the laws of physics are the same today and tomorrow, or here and there, then some physical characteristics of a mechanical system must remain invariant, that is unchanged.

Symmetry

Imagine a physicist performs an experiment in a box.
We then rotate the box. If the physicist repeats this experiment with exactly the same initial configurations and finds no difference in the results, we say the laws governing the experiment are rotationally symmetric.

Similarly, if there are no differences in experiments under spatial or temporal translation, then the laws describing the experiment are spatially or temporally invariant.
Spatial translation is when the physicist repeats the experiment at a different location.
Temporal translation is when the physicist repeats the same experiment after a certain waiting period.

The mechanical laws are represented by the equations of motion.
So, "the same results of an experiment" means that the laws (that is, the equations of motion) are the same (invariant) relative to a corresponding transformation.

In Lagrangian Mechanics, this idea of symmetry takes a precise formulation.

A symmetry is a continuous transformation of time t, space coordinates q={qi}, or other dynamical variables under which the Lagrangian L describing the mechanical system remains unchanged or changes only by a total time derivative dF(t,q(t),q'(t))/dt of some smooth function F(t,q(t),q'(t)), where q'(t) signifies a time-derivative of a coordinate q, that is velocity.

A nuance related to a total time derivative of some smooth function F(t,q(t),q'(t)) might need an explanation.
If the Lagrangian will change by dF(t,q(t),q'(t))/dt, the action functional
Φ[L(t)] = ab[L(t,q,q')+dF/dt]dt
defined on a trajectory with endpoints q(a)=A and q(b)=B fixed will change by the boundary term F(b)−F(a).
Considering the end points of a trajectory are fixed as initial conditions, the extremals will be preserved, leading to the same Euler–Lagrange equations and the same physical trajectories.

Let's rigorously define the symmetry as a continuous transformation we are talking about.
Consider a configuration space of n degrees of freedom with generalized coordinates q={q1,...,qn} and time parameter t.

Consider a family of transformation functions parameterized by variable ε
t(ε)=T(ε,t,q1,...,qn)=T(ε,t,q)
and (for each i∈[1,n])
qi(ε)=Qi(ε,t,q1,...,qn)=Qi(ε,t,q)
define a transformations from original space-time coordinates {t,q} to the new ones {t(ε),q(ε)} if t(0)=t and qi(0)=qi for each i∈[1,n].

Sometimes for brevity we will use
q(ε)=Q(ε,t,q)
implying qi(ε) on the left and Qi(ε,t,q) on the right side for all i from 1 to n.

Functions T(ε,t,q) and Qi(ε,t,q) are assumed to be continuous with respect to the first argument ε for any t and q.

Usually, we will require even a stronger requirements for transformational functions T(ε,t,q) and Qi(ε,t,q) to uniformly converge to T(0,t,q) and Qi(0,t,q) as ε→0.

In addition, to define it more rigorously, we require that all transformations make up a group (more precisely, the Lie group) in a sense that transformations can be combined:
{t,q}→{t(α),q(α)}→{t(β),q(β)}
should be equivalent to
{t,q}→{t(α+β),q(α+β)}

In particular, any (ε)-transformation is reversible using (−ε)-transformation because of the additivity mentioned above:
(ε)+(−ε)=(0).
and, as we stated above,
t = t(0) = T(0,t,q)
q = q(0) = Q(0,t,q)

Functions T(ε,t,q) and Q(ε,t,q) define the parameterized by ε transformation of time and position.
Since, ultimately, we need to check the Lagrangian for invariance to transformations, we also need to know how velocities are changing with transformation because Lagrangian depends on them.

The derivation of a formula that defines this change is rather complex and the rest of this lecture is dedicated precisely to dealing with transformation of velocity.

If transformation does not involve time then the transformed velocity is just a time-derivative of the transformed position
qi'(ε) = dqi(ε)/dt = qi(ε)'
Notice the order of ε-transformation and apostrophe that indicates a time derivative.
Thus, in this particular case
qi'(ε) = dQ(ε,t,q)/dt

The complication in our case is, not only coordinates are changing, but time as well.
So, the correct definition of velocity transformation is
qi'(ε) = dqi(ε)/dt(ε)

Let's start with some simple manipulation.
dqi(ε)/dt(ε) =
=
[dqi(ε)/dt] / [dt(ε)/dt] =
=
[dQi(ε,t,q)/dt] / [dT(ε,t,q)/dt]

The numerator, by chain rule, equals to
∂Qi/∂t+Σk(∂Qi/∂qk)·qk'
The denominator, by chain rule, equals to
∂T/∂t+Σk(∂T/∂qk)·qk'

Therefore,
qi'(ε) =
∂Qi /∂t+Σk(∂Qi /∂qk)·qk'
∂T/∂t+Σk(∂T/∂qk)·qk'
where
Qi = Qi(ε,t,q) = qi(ε)
T = T(ε,t,q) = t(ε)

As ε0, we can represent ε-transformed values {t(ε),qi(ε)} in terms of initial values {t(0)=t,qi(0)=qi} and infinitesimal increments like in the Taylor's formula:
t(ε)= t + ε·[∂t(ε)/∂ε|ε=0] +o(ε)
qi(ε)= qi + ε·[∂qi(ε)/∂ε|ε=0] +o(ε)

An expression
τ(t,q) = ∂t(ε)/∂ε|ε=0 =
= ∂T(ε,t,q)/∂ε|ε=0

is a derivative of the time transformation function by parameter of transformation ε at the initial point of transformation ε=0 and is called the infinitesimal generator of time transformation.
So, in terms of this generator,
t(ε)= t + ε·τ(t,q) +o(ε)

An expression
ξi(t,q) = ∂qi(ε)/∂ε|ε=0 =
= ∂Qi(ε,t,q)/∂ε|ε=0

is a derivative of the space transformation function by parameter of transformation ε at the initial point of transformation ε=0 and is called the infinitesimal generator of space transformation.
So, in terms of this generator,
qi(ε)= qi + ε·ξi(t,q) +o(ε)

Using the above representations, we can evaluate the terms of the formula for qi'(ε).

∂T(ε,t,q)/∂t =
= ∂t(ε)/∂t =
= ∂
[t+ε·τ(t,q)+o(ε)]/∂t
derivative of t by t is 1;
ε is an independent of t parameter

= 1 + ε·∂τ(t,q)/∂t


∂T(ε,t,q)/∂qk =
= ∂t(ε)/∂qk =
= ∂
[t+ε·τ(t,q)+o(ε)]/∂qk
derivative of t by qk is 0;
ε is an independent of qk parameter

= ε·∂τ(t,q)/∂qk


∂Qi (ε,t,q)/∂t =
= ∂qi(ε)/∂t =
= ∂
[qi+ε·ξi(t,q)+o(ε)]/∂t
qi is the initial position in space (before transformation is applied) and, when used in partial derivatives of function that depends on it and other parameters, like time t, is treated as independent of other parameters, thus derivative ∂qi/dt=0;
ε is an independent of t parameter

= ε·∂ξi(t,q)/∂t


∂Qi (ε,t,q)/∂qk =
= ∂qi(ε)/∂qk =
= ∂
[qi+ε·ξi(t,q)+o(ε)]/∂qk
derivative of qi by qk is 1 for i=k and 0 otherwise ⇒ (using math symbol δ) derivative is equal to δik;
ε is an independent of t parameter

= δik + ε·∂ξi(t,q)/∂qk


We are ready to evaluate numerator and denominator of the above expression for qi'(ε).

Numerator =
= ε·∂ξi(t,q)/∂t +
+
Σk[δik+ε·∂ξi(t,q)/∂qk]·qk' =
= ε·∂ξi(t,q)/∂t +
+ qi' + ε·
Σk[∂ξi(t,q)/∂qk]·qk' =
= qi' +
+ε·
{∂ξi(t,q)/∂t+ Σk[∂ξi(t,q)/∂qk]·qk'}

Denominator =
= 1 + ε·∂τ(t,q)/∂t +
+
Σk[ε·∂τ(t,q)/∂qk]·qk' =
= 1 +
+ε·
{∂τ(t,q)/∂t+ Σk[∂τ(t,q)/∂qk]·qk'}

Now the formula for a transformation of velocity looks like
qi'(ε) =
qi' + ε·A
1 + ε·B
where
A = ∂ξi(t,q)/∂t +
+
Σk[∂ξi(t,q)/∂qk]·qk'
and
B = ∂τ(t,q)/∂t +
+
Σk[∂τ(t,q)/∂qk]·qk'

Let's not forget our goal - to establish how the velocity is transformed by ε-transformation, which implies differentiation of the transformed velocity by a transformation parameter ε at point ε=0 to find instantaneous velocity rate of change at the initial position {t,q}.

If the above expression would be of the form
qi'(ε) − qi' = ε·X + o(ε)
then we could say that
limε→0 (qi'(ε) − qi')/ε =
= d(qi'(ε))/dε|ε=0 = X


To achieve this goal, we will convert our formula
qi'(ε) =
qi' + ε·A
1 + ε·B
into more desirable form using this easily verifiable identity
(1+ε·B)·[qi' + ε·(A−B·qi')] =
= qi' + ε·A + o(ε)


Dividing both parts by (1+ε·B) will result in
qi' + ε·(A−B·qi') =
=
[qi' + ε·A + o(ε)] / (1+ε·B) =
=
[qi' + ε·A] / (1+ε·B) + o(ε) =
= qi'(ε) + o(ε)

Therefore,
qi'(ε) = qi' + ε·(A−B·qi') + o(ε)

Hence, the instantaneous rate of change of velocity at point {t,q} is
d(qi'(ε))/dε|ε=0 = A−B·qi' =
=
{∂ξi(t,q)/∂t +
+
Σk[∂ξi(t,q)/∂qk]·qk'}
− qi
{∂τ(t,q)/∂t +
+
Σk[∂τ(t,q)/∂qk]·qk'}

In the above expression we can recognize two full time derivatives
i(t,q)/dt = ∂ξi(t,q)/∂t +
+
Σk[∂ξi(t,q)/∂qk]·qk'
and
dτ(t,q)/dt = ∂τ(t,q)/∂t +
+
Σk[∂τ(t,q)/∂qk]·qk'
which makes our formula for a rate of change of velocity much simpler
d(qi'(ε))/dε|ε=0 =
= dξi(t,q)/dt − qi'·dτ(t,q)/dt

where
ξi(t,q)=∂qi(ε)/∂ε|ε=0
represents the rate of change of position by infinitesimal ε-transformation, whose time derivative (the first term in the above formula) is the rate of change of velocity without taking into consideration the transformation of time
and
τ(t,q)=∂t(ε)/∂ε|ε=0
whose time derivative multiplied by qi' (the second term in the above formula) represents an adjustment caused by ε-transformation of time.

Tuesday, March 17, 2026

Noether's Theorem and Angular Momentum Conservation:UNIZOR.COM -> Physics+ 4 All -> Lagrangian - Noether l = m·r⨯v const

Notes to a video lecture on UNIZOR.COM

Noether's Theorem
Angular Momentum Conservation

History

In 1918, German mathematician Emmy Noether published the proof of a very important theorem named after her.
Albert Einstein and other physicists considered her one of the most significant mathematicians of her time.
Noether's theorem has been called one of the most important mathematical theorems guiding the development of modern physics.
This lecture is about a particular case of Noether's theorem as it applies to classical mechanics using the Lagrangian approach.

Background

Recall these important items from the prior lectures of this course Physics+ 4 All on UNIZOR.COM.

1. Configuration space with generalized coordinates q(t)={qi(t)}, where i∈[1,n] and n is the number of degrees of freedom.
Generalized velocities (time derivatives of generalized coordinates) q'(t)={qi'(t)}.

2. Lagrangian
L(q,q') = T(q,q') − U(q)
- the difference between the total kinetic T and total potential U energies of a system.
Optionally, the Lagrangian might be explicitly dependent on time, but in this and subsequent lectures we assume that energies and, therefore, the Lagrangian do not explicitly depend on time, but only on positions and velocities. This makes the Lagrangian depend on time implicitly through the motion q(t).
Systems with the Lagrangian explicitly depending on time (like in case of variable gravitational field, driven oscillations or time-dependent electromagnetic field) are not considered here.

3. Action functional
Φ[L(q(t),q'(t))t∈[t1,t2]] =
= [t1,t2] L(q(t),q'(t))dt

4. Euler-Lagrange equations
These equations are established for each generalized coordinate
d/dt ∂L/∂qi' = ∂L/∂qi (i∈[1,n])
These are differential equations with unknown position functions q(t).
Their solutions are the extremals of the action functional above and, at the same time, the only candidates for real physical trajectories of a mechanical system in its configuration space with generalized coordinates.

5. Generalized Momentum
By definition, the generalized momentum is a set of partial derivatives of the Lagrangian by generalized velocities:
pi = ∂L/∂qi'
Previous lecture was dedicated to the proof of the Momentum Conservation law for each component pi, for which the Lagrangian is invariant (symmetrical) under a translation (shift) of the corresponding generalized cyclic coordinate qiqi:
∂L/∂qi=0 pi=const
This momentum conservation law is essential for this lecture about angular momentum.

Angular Momentum Conservation

Our purpose is to prove a particular case of Noether's theorem about the law of conservation of angular momentum of a mechanical system.

But instead of directly proving this for a particular case of rotational symmetry, we will use the already proven in the previous lecture result for a transformation of generalized coordinates.

Recall the general theorem proven in the previous lecture:
If the Lagrangian
L(q1,...,qn,q1',...,qn',t)
is invariant under translation of qk by infinitesimal value ε
qkqk + ε
then the kth coordinate of the generalized momentum
pk = ∂L/∂qk'
is conserved.

Let's use this theorem for a special case of a conservative planar system (two-dimensional system on the Euclidean plane) moving in a central field (the one with a potential depending only on a distance from a central point, like gravitational field) with polar coordinates:
radius r that we will interpret as a generalized coordinate q1 and
angle θ that we will interpret as a generalized coordinate q2.

A rotation of a physical system in the plane corresponds to a translation of the angular coordinate
θθ + ε

Let's express the Lagrangian L(r,θ,r',θ') in terms of these variables.
L = T − U
is the difference between kinetic and potential energies.

Kinetic energy
T(r,θ,r',θ') = ½m·v²
depends on the mass m and the magnitude of its speed v that in polar coordinates can be calculated using Cartesian coordinates (x,y) as follows:
x = r·cos(θ)
y = r·sin(θ)
x' = r'·cos(θ)−r·sin(θ)·θ'
y' = r'·sin(θ)+r·cos(θ)·θ'
v² = (x')² + (y')²
(x')² = (r')²·cos²(θ) −
− 2r'·cos(θ)·r·sin(θ)·θ' +
+ r²·sin²(θ)·(θ')²

(y')² = (r')²·sin²(θ) +
+ 2r'·sin(θ)·r·cos(θ)·θ' +
+ r²·cos²(θ)·(θ')²


After cancelling plus and minus of the middle term in a sum of two last expressions above and taking into consideration that
sin²(θ) + cos²(θ) = 1
we obtain
v² = (x')² + (y')² =
= (r')² + r²·(θ')²


Therefore, kinetic energy T is
T = ½m·[(r')² + r²·(θ')²]
which is independent of angle θ and invariant to its translation θθ+ε.

Potential energy U of a central field is, as we mentioned above, depends only on the distance from the central point. If our polar system of coordinates has an origin in that point, we can express the potential energy as U(r) - also independent of angle θ and, therefore, invariant to its translation.

A rotation of a physical system in the plane corresponds to a translation of the angular coordinate θθ+ε, but the Lagrangian of a system in a central field is invariant under the translation θθ+ε, that is angle θ is cyclic and
∂L/∂θ = 0

From this and the generalized Momentum Conservation law proven in the previous lecture follows that the corresponding momentum
pθ = ∂L/∂θ'
is conserved.

Let's express the momentum pθ in term of polar coordinates, taking into account that potential energy U does not depend on velocities.
pθ = ∂L/∂θ' = ∂(T−U)/∂θ' =
= ∂T/∂θ' =
= ∂/∂θ'
{½m·[(r')² + r²·(θ')²]} =
= m·r²·θ'


Thus,
pθ = m·r²·θ'
is a constant of motion of a conservative planar system in a central field.

The expression m·r²·θ' is exactly the magnitude of the angular momentum vector of a point-mass moving about the origin of polar coordinates (see Note 1 below), which in classical mechanics is defined as the mass multiplied by a vector product of radius and velocity vectors and directed perpendicular to a plane of rotation
𝓁 = m·rv

Thus the conserved generalized momentum corresponding to rotational symmetry is the angular momentum.

As a conclusion for this and the previous lecture, the translational symmetry leads to conservation of linear momentum, while rotational symmetry leads to conservation of angular momentum
_________
Note 1
The magnitude of the angular momentum vector can be calculated from its definition
𝓁 = m·rv
by expressing vectors in their Cartesian coordinate form, directing the X-axis along vector r, Y-axis to be perpendicular to it within a plane of motion and Z-axis to be perpendicular to other two axes and, therefore, perpendicular to an entire plane of motion.
Then
r = (r,0,0)
v = (r',r·θ',0)
Their vector product is the determinant of a matrix with i, j and k being unit vectors along the corresponding axes
ijk
r00
r'r·θ'0
This determinat equals to
i·0 + j·0 + k·r·r·θ'
Which means that the angular momentum vector is perpendicular to the plane of motion and its magnitude equals to r²·θ'.