D D P T: Los Alamos Electronic Archives: Physics/9909035
D D P T: Los Alamos Electronic Archives: Physics/9909035
r
X
i
v
:
p
h
y
s
i
c
s
/
9
9
0
9
0
3
5
v
1
1
9
S
e
p
1
9
9
9
Los Alamos Electronic Archives: physics/9909035
CLASSICAL MECHANICS
HARET C. ROSU
rosu@ifug3.ugto.mx
a
d
d
p
t
( ) r
a
p
V
graduate course
Copyright c 1999 H.C. Rosu
Le on, Guanajuato, Mexico
v1: September 1999.
1
CONTENTS
1. THE MINIMUM PRINCIPLES ... 3.
2. MOTION IN CENTRAL FORCES ... 19.
3. RIGID BODY ... 32.
4. SMALL OSCILLATIONS ... 52.
5. CANONICAL TRANSFORMATIONS ... 70.
6. POISSON PARENTHESES... 79.
7. HAMILTON-JACOBI EQUATIONS... 82.
8. ACTION-ANGLE VARIABLES ... 90.
9. PERTURBATION THEORY ... 96.
10. ADIABATIC INVARIANTS ... 111.
11. MECHANICS OF CONTINUOUS SYSTEMS ... 116.
2
1. THE MINIMUM PRINCIPLES
Forward: The history of minimum principles in physics is long and inter-
esting. The study of such principles is based on the idea that the nature acts
always in such a way that the important physical quantities are minimized
whenever a real physical process takes place. The mathematical background
for these principles is the variational calculus.
CONTENTS
1. Introduction
2. The principle of minimum action
3. The principle of DAlembert
4. Phase space
5. The space of congurations
6. Constraints
7. Hamiltons equations of motion
8. Conservation laws
9. Applications of the action principle
3
1. Introduction
The empirical evidence has shown that the motion of a particle in an inertial
system is correctly described by Newtons second law
F = d p/dt, whenever
possible to neglect the relativistic eects. When the particle happens not to
be forced to a complicated motion, the Cartesian coordinates are sucient
to describe the movement. If none of these conditions are fullled, rather
complicated equations of motion are to be expected.
In addition, when the particle moves on a given surface, certain forces called
constraint forces must exist to maintain the particle in contact with the
surface. Such forces are not so obvious from the phenomenological point
of view; they require a separate postulate in Newtonian mechanics, the one
of action and reaction. Moreover, other formalisms that may look more
general have been developed. These formalisms are equivalent to Newtons
laws when applied to simple practical problems, but they provide a general
approach for more complicated problems. The Hamilton principle is one
of these methods and its corresponding equations of motion are called the
Euler-Lagrange equations.
If the Euler-Lagrange equations are to be a consistent and correct descrip-
tion of the dynamics of particles, they should be equivalent to Newtons
equations. However, Hamiltons principle can also be applied to phenom-
ena generally not related to Newtons equations. Thus, although HP does
not give a new theory, it unies many dierent theories which appear as
consequences of a simple fundamental postulate.
The rst minimum principle was developed in the eld of optics by Heron
of Alexandria about 2,000 years ago. He determined that the law of the
reection of light on a plane mirror was such that the path taken by a light
ray to go from a given initial point to a given nal point is always the shortest
one. However, Herons minimum path principle does not give the right law
of reection. In 1657, Fermat gave another formulation of the principle by
stating that the light ray travels on paths that require the shortest time.
Fermats principle of minimal time led to the right laws of reection and
refraction. The investigations of the minimum principles went on, and in
the last half of the XVII century, Newton, Leibniz and Bernoulli brothers
initiated the development of the variational calculus. In the following years,
Lagrange (1760) was able to give a solid mathematical base to this principle.
4
In 1828, Gauss developed a method of studying Mechanics by means of his
principle of minimum constraint. Finally, in a sequence of works published
during 1834-1835, Hamilton presented the dynamical principle of minimum
action.This principle has always been the base of all Mechanics and also of
a big part of Physics.
Action is a quantity of dimensions of length multiplied by the momentum
or energy multiplied by time.
2. The action principle
The most general formulation of the law of motion of mechanical systems
is the action or Hamilton principle. According to this principle every me-
chanical system is characterized by a function dened as:
L
_
q
1
, q
2
, ..., q
s
,
q
1
,
q
2
,
q
s
, t
_
,
or shortly L
_
q,
q
, t
_
, and the motion of the system satises the following
condition: assume that at the moments t
1
and t
2
the system is in the posi-
tions given by the set of coordinates q
(1)
y q
(2)
; the system moves between
these positions in such a way that the integral
S =
_
t
2
t
1
L
_
q,
q
, t
_
dt (1)
takes the minimum possible value. The function L is called the Lagrangian
of the system, and the integral (1) is known as the action of the system.
The Lagrange function contains only q and
q
, and no other higher-order
derivatives. This is because the mechanical state is completely dened by
its coordinates and velocities.
Let us establish now the difererential equations that determine the minimum
of the integral (1). For simplicity we begin by assuming that the system has
only one degree of freedom, therefore we are looking for only one function
q (t). Let q = q (t) be the function for which S is a minimum. This means
that S grows when one q (t) is replaced by an arbitrary function
q (t) +q (t) , (2)
where q (t) is a small function through the interval from t
1
to t
2
[it is called
the variation of the function q (t)]. Since at t
1
and t
2
all the functions (2)
should take the same values q
(1)
and q
(2)
, one gets:
5
q (t
1
) = q (t
2
) = 0. (3)
What makes S change when q is replaced by q + q is given by:
_
t
2
t
1
L
_
q + q,
q
+
q
, t
_
dt
_
t
2
t
1
L
_
q,
q
, t
_
dt.
An expansion in series of this dierence in powers of q and
q
begins by
terms of rst order. The necessary condition of minimum (or, in general,
extremum) for S is that the sum of all terms turns to zero; Thus, the action
principle can be written down as follows:
S =
_
t
2
t
1
L
_
q,
q
, t
_
dt = 0, (4)
or by doing the variation:
_
t
1
t
2
_
L
q
q +
L
q
_
dt = 0 .
Taking into account that
q
= d/dt (q), we make an integration by parts
to get:
S =
_
L
q
q
_
t
2
t
1
+
_
t
1
t
2
_
L
q
d
dt
L
q
_
qdt = 0 . (5)
Considering the conditions (3), the rst term of this expresion disappears.
Only the integral remains that should be zero for all values of q. This is
possible only if the integrand is zero, which leads to the equation:
L
q
d
dt
L
q
= 0 .
For more degrees of freedom, the s dierent functions q
i
(t) should vary
independently. Thus, it is obvious that one gets s equations of the form:
d
dt
_
L
q
i
_
L
q
i
= 0 (i = 1, 2, ..., s) (6)
These are the equations we were looking for; in Mechanics they are called
Euler-Lagrange equations. If the Lagrangian of a given mechanical system
6
is known, then the equations (6) form the relationship between the accel-
erations, the velocities and the coordinates; in other words, they are the
equations of motion of the system. From the mathematical point of view,
the equations (6) form a system of s dierential equations of second order
for s unknown functions q
i
(t). The general solution of the system contains
2s arbitrary constants. To determine them means to completely dene the
movement of the mechanical system. In order to achieve this, it is neces-
sary to know the initial conditions that characterize the state of the system
at a given moment (for example, the initial values of the coordinates and
velocities.
3. DAlembert principle
The virtual displacement of a system is the change in its congurational
space under an arbitrary innitesimal variation of the coordinates r
i
,which
is compatible with the forces and constraints imposed on the system at the
given instant t. It is called virtual in order to distinguish it from the real
one, which takes place in a time interval dt, during which the forces and the
constraints can vary.
The constraints introduce two types of diculties in solving mechanics prob-
lems:
(1) Not all the coordinates are independent.
(2) In general, the constraint forces are not known a priori; they are some
unknowns of the problem and they should be obtained from the solution
looked for.
In the case of holonomic constraints the diculty (1) is avoided by introduc-
ing a set of independent coordinates (q
1,
q
2,...,
q
m
, where m is the number of
degrees of freedom involved). This means that if there are mconstraint equa-
tions and 3N coordenates (x
1
, ..., x
3N
), we can eliminate these n equations
by introducing the independent variables (q
1
, q
2
, .., , q
n
). A transformation
of the following form is used
x
1
= f
1
(q
1
, ..., q
m
, t)
.
.
.
x
3N
= f
3N
(q
1
, ..., q
n
, t) ,
where n = 3N m.
To avoid the diculty (2) Mechanics needs to be formulated in such a way
that the forces of constraint do not occur in the solution of the problem.
This is the essence of the principle of virtual work.
7
Virtual work: We assume that a system of N particles is described by
3N coordenates (x
1
, x
2
, ..., x
3N
) and let F
1,
F
2,...,
F
3N
be the components of
the forces acting on each particle. If the particles of the system display
innitesimal and instantaneous displacements x
1
, x
2
, ..., x
3N
under the
action of the 3N forces, then the performed work is:
W =
3N
j=1
F
j
x
j
. (7)
Such displacements are known as virtual displacements and W is called
virtual work; (7) can be also written as:
W =
N
=1
F
r . (8)
Forces of constraint: besides the applied forces F
(e)
.
The principle of virtual work: Let F
= F
(e)
+ R
. (9)
and if the system is in equilibrium, then
F
= F
(e)
+R
= 0 . (10)
Thus, the virtual work due to all possible forces F
is:
W =
N
=1
F
=
N
=1
_
F
(e)
+R
_
r
= 0 . (11)
If the system is such that the constraint forces do not make virtual work,
then from (11) we obtain:
N
=1
F
(e)
= 0 . (12)
Taking into account the previous denition, we are now ready to introduce
the DAlembert principle. According to Newton, the equation of motion is:
8
F
= 0 ,
which tells that the particles of the system would be in equilibrium under the
action of a force equal to the real one plus an inverted force
p
i
. Instead
of (12) we can write
N
=1
_
F
_
r
= 0 (13)
and by doing the same decomposition in applied and constraint forces (f
),
we obtain:
N
=1
_
F
(e)
_
r
+
N
=1
f
= 0 .
Again, let us limit ourselves to systems for which the virtual work due to
the forces of constraint is zero leading to
N
=1
_
F
(e)
_
r
= 0 , (14)
which is the DAlemberts principle. However, this equation does not have
a useful form yet for getting the equations of motion of the system. There-
fore, we should change the principle to an expression entailing the virtual
displacements of the generalized coordinates, which being independent from
each other, imply zero coecients for q
=
dr
dt
=
k
r
q
k
q
k
+
r
t
where r
= r
(q
1
, q
2
, ..., q
n
, t) .
Similarly, the arbitrary virtual displacement r
=
j
r
q
j
q
j
.
9
Then, the virtual work F
=1
F
=
j,
F
q
j
q
j
=
j
Q
j
q
j
, (15)
where the Q
j
are the so-called components of the generalized force, dened
in the form
Q
j
=
q
j
.
Now if we see eq. (14) as:
p
r
(16)
and by substituting in the previous results we can see that (16) can be
written:
_
d
dt
_
m
q
j
_
m
q
j
_
=
j
__
d
dt
_
T
q
j
_
T
q
j
_
Q
j
_
q
j
= 0 .
(17)
The variables q
j
can be an arbitrary system of coordinates describing the
motion of the system. However, if the constraints are holonomic, it is possi-
ble to nd systems of independent coordinates q
j
containing implicitly the
constraint conditions already in the equations of transformation x
i
= f
i
if
one nullies the coecients by separate:
d
dt
_
T
q
_
T
q
= Q
j
. (18)
There are m equations. The equations (18) are sometimes called the La-
grange equations, although this terminology is usually applied to the form
they take when the forces are conservative (derived from a scalar potential
V)
F
=
i
V.
Then Q
j
can be written as:
10
Q
j
=
V
q
j
.
The equations (18) can also be written in the form:
d
dt
_
T
q
j
_
(T V )
q
j
= 0 (19)
and dening the Lagrangian L in the form L = T V one gets
d
dt
_
L
q
j
_
L
q
j
= 0 . (20)
These are the Lagrange equations of motion.
4. - Phase space
In the geometrical interpretation of mechanical phenomena the concept of
phase space is very much in use. It is a space of 2s dimensions whose axes
of coordinates are the s generalized coordinates and the s momenta of the
given system. Each point of this space corresponds to a denite mechanical
state of the system. When the system is in motion, the representative point
in phase space performs a curve called phase trajectory.
5. - Space of congurations
The state of a system composed of n particles under the action of m con-
straints connecting some of the 3n cartesian coordinates is completely de-
termined by s = 3n m generalized coordinates. Thus, it is possible to
descibe the state of such a system by a point in the s dimensional space
usually called the conguration space, for which each of its dimensions cor-
responds to one q
j
. The time evolution of the system will be represented by a
curve in the conguration space made of points describing the instantaneous
conguration of the system.
6. - Constraints
One should take into account the constraints that act on the motion of the
system. The constraints can be classied in various ways. In the general
case in which the constraint equations can be written in the form:
i
c
i
q
i
= 0 ,
11
where the c
i
are functions of only the coordinates (the index counts
the constraint equations). If the rst members of these equations are not
total derivatives with respect to the time they cannot be integrated. In other
words, they cannot be reduced to relationships between only the coordinates,
that might be used to express the position by less coordinates, corresponding
to the real number of degrees of freedom. Such constraints are called non
holonomic (in contrast to the previous ones which are holonomic and which
connect only the coordinates of the system).
7. Hamiltons equations of motion
The formulation of the laws of Mechanics by means of the Lagrangian as-
sumes that the mechanical state of the system is determined by its gener-
alized coordinates and velocities. However, this is not the unique possible
method; the equivalent description in terms of its generalized coordinates
and momenta has a number of advantages.
Turning fromone set of independent variables to another one can be achieved
by what in mathematics is called Legendre transformation. In this case the
transformation takes the following form where the total dierential of the
Lagrangian as a function of coordinates and velocities is:
dL =
i
L
q
i
dq
i
+
i
L
q
i
d
q
i
,
that can be written as:
dL =
i
p
i
dq
i
+
i
p
i
d
q
i
, (21)
where we already know that the derivatives L/
q
i
, are by denition the
generalized momenta and moreover L / q
i
=
p
i
by Lagrange equations.
The second term in eq. (21)can be written as follows
i
p
i
d
q
i
= d
_
p
i
q
i
_
q
i
dq
i
.
By attaching the total diferential d
_
p
i
q
i
_
to the rst term and changing
the signs one gets from (21):
d
_
p
i
q
i
L
_
=
p
i
dq
i
+
p
i
q
i
. (22)
12
The quantity under the diferential is the energy of the system as a func-
tion of the coordinates and momenta and is called Hamiltonian function or
Hamiltonian of the system:
H (p, q, t) =
i
p
i
q
i
L . (23)
Then from ec. (22)
dH =
p
i
dq
i
+
p
i
q
i
where the independent variables are the coordinates and the momenta, one
gets the equations
q
i
=
H
p
i
p
i
=
H
q
i
. (24)
These are the equations of motion in the variables q y p and they are called
Hamiltons equations.
8. Conservation laws
8.1 Energy
Consider rst the conservation theorem resulting from the homogeneity of
time. Because of this homogeneity, the Lagrangian of a closed system does
not depend explicitly on time. Then, the total time diferential of the La-
grangian (not depending explicitly on time) can be written:
dL
dt
=
i
L
q
i
q
i
+
i
L
q
i
q
i
and according to the Lagrange equations we can rewrite the previous equa-
tion as follows:
dL
dt
=
q
i
d
dt
_
L
q
i
_
+
i
L
q
i
q
i
=
i
d
dt
_
q
i
L
q
i
_
,
or
i
d
dt
_
q
i
L
q
i
L
_
= 0 .
From this one concludes that the quantity
E
q
i
L
q
i
L (25)
13
remains constant during the movement of the closed system, that is it is
an integral of motion. This constant quantity is called the energy E of the
system.
8.2 Momentum
The homogeneity of space implies another conservation theorem. Because
of this homogeneity, the mechanical properties of a cosed system do not
vary under a parallel displacement of the system as a whole through space.
We consider an innitesimal displacement (i.e., the position vectors r
turned into r
a
+) and look for the condition for which the Lagrangian does
not change. The variation of the function L resulting from the innitesi-
mal change of the coordinates (maintaining constant the velocities of the
particles) is given by:
L =
a
L
r
a
r
a
=
a
L
r
a
,
extending the sum over all the particles of the system. Since is arbitrary,
the condition L = 0 is equivalent to
a
L
r
a
= 0 (26)
and taking into account the already mentioned equations of Lagrange
a
d
dt
_
L
v
a
_
=
d
dt
a
L
v
a
= 0 .
Thus, we reach the conclusion that for a closed mechanical system the vec-
torial quantity called impetus/momentum
P
a
L
v
a
remains constant during the motion.
8.3 Angular momentum
Let us study now the conservation theorem coming out from the isotropy of
space. For this we consider an innitesimal rotation of the system and look
for the condition under which the Lagrangian does not change.
We shall call an innitesimal rotation vector a vector of modulus equal
to the angle of rotation and whoose direction coincide with that of the
rotation axis. We shall look rst to the increment of the position vector
14
of a particle in the system, by taking the origin of coordinates on the axis
of rotation. The lineal displacement of the position vector as a function of
angle is
[r[ = r sin ,
(see the gure). The direction of the vector r is perpendicular to the plane
dened by r and , and therefore,
r = r . (27)
r
r
O
The rotation of the system changes not only the directions of the position
vectors but also the velocities of the particles that are modied by the same
rule for all the vectors. The velocity increment with respect to a xed frame
system will be:
v = v .
We apply now to these expressions the condition that the Lagrangian does
not vary under rotation:
L =
a
_
L
r
a
r
a
+
L
v
a
v
a
_
= 0
and substituting the denitions of the derivatives L/v
a
por p
a
and L/r
a
from the Lagrange equations by
p
a
; we get
a
_
p
a
r
a
+p
a
v
a
_
= 0 ,
15
or by circular permutation of the factors and getting out of the sum:
a
_
r
a
p
a
+v
a
p
a
_
=
d
dt
a
r
a
p
a
= 0 ,
because is arbitrary, one gets
d
dt
a
r
a
p
a
= 0
Thus, the conclusion is that during the motion of a closed system the vec-
torial quantity called the angular (or kinetic) momentum is conserved.
M
a
r
a
p
a
.
9.- Applications of the action principle
a) Equations of motion
Find the eqs of motion for a pendular mass sustained by a resort, by directly
applying Hamiltons principle.
k
g
m
For the pendulum in the gure the Lagrangian function is
16
L =
1
2
m(
r
2
+r
2
2
) +mgr cos
1
2
k(r r
o
)
2
,
therefore
_
t
2
t
1
Ldt =
_
t
2
t
1
_
m
_
r
r +r
2
+r
2
_
+mgr cos mgr sin k(r r
o
)r
_
dt
m
r
r dt = m
r d(r) = d
_
m
r r
_
mr
r dt .
In the same way
mr
2
dt = d
_
mr
2
_
d
_
mr
2
_
dt
dt
= d
_
mr
2
_
mr
2
+2mr
r
_
dt .
Therefore, the previous integral can be written
_
t
2
t
1
__
m
r mr
2
mg cos + k (r r
o
)
_
+
_
mr
2
+2mr
r
+mgr sin
_
_
dt
_
t
2
t
1
_
d
_
m
r r
_
+d
_
mr
2
2
__
= 0 .
Assuming that both r and are equal zero at t
1
and t
2
, the second integral
is obviously nought. Since r and are completely independent of each
other, the rst integral can be zero only if
m
r mr
2
mg cos +k(r r
o
) = 0
and
mr
2
+2mr
r
+mgr sin = 0 ,
These are the equations of motion of the system.
17
b) Exemple of calculating a minimum value
Prove that the shortest line between two given points p
1
and p
2
on a cilinder
is a helix.
The length S of an arbitrary line on the cilinder between p
1
and p
2
is given
by
S =
_
p
2
p
1
_
1 +r
2
_
d
dz
_
2
_
1/2
dz ,
where r, and z are the usual cilindrical coordinates for r = const. A
relationship between and z can be determined for which the last integral
has an extremal value by means of
d
dz
_
= 0 ,
where =
_
1 + r
2
1/2
y
=
d
dz
, but since / = 0 we have
=
_
1 + r
2
2
_
1/2
r
2
= c
1
= const. ,
therefore r
= c
2
. Thus, r = c
2
z +c
3
, which is the parametric equation of
a helix. Assuming that in p
1
we have = 0 and z = 0, then c
3
= 0. In p
2
,
make =
2
and z = z
2
, therefore c
2
= r
2
/z
2
, and r = (r
2
/z
2
) z is the
nal equation.
References
L. D. Landau and E. M Lifshitz, Mechanics, Theoretical Physics, vol I,
(Pergammon, 1976)
H. Goldstein, Classical Mechanics, (Addison-Wesley, 1992)
18
2. MOTION IN CENTRAL FORCES
Forward: Because of astronomical reasons, the motion under the action of
central forces has been the physical problem on which pioneer researchers
focused more, either from the observational standpoint or by trying to disen-
tangle the governing laws of motion. This movement is a basic example for
many mathematical formalisms. In its relativistic version, Keplers problem
is yet an area of much interest.
CONTENTS:
2.1 The two-body problem: reduction to the one-body problem
2.2 Equations of motion
2.3 Dierential equation of the orbit
2.4 Keplers problem
2.5 Dispertion by a center of forces (with example)
19
2.1 Two-body problem: Reduction to the one-body problem
Consider a system of two material points of masses m
1
and m
2
, in which
there are forces due only to an interaction potential V . We suppose that V
is a function of any position vector between m
1
and m
2
, r
2
r
1
, or of their
relative velocities
r
2
r
1
, or of the higher-order derivatives of r
2
r
1
. Such
a system has 6 degrees of freedom and therefore 6 independent generalized
coordinates.
We suppose that these are the vector coordinates of the center-of-mass R,
plus the three components of the relative dierence vector r = r
2
r
1
. The
Lagrangian of the system can be written in these coordinates as follows:
L = T(
R, r) V (r, r, r, .....). (1)
The kinetic energy T is the sum of the kinetic energy of the center-of-mass
plus the kinetic energy of the motion around it, T:
T =
1
2
(m
1
+ m
2
)
R
2
+T,
being
T=
1
2
m
1
r
2
1
+
1
2
m
2
r
2
2
.
Here, r
1
and r
2
are the position vectors of the two particles with respect to
the center-of-mass, and they are related to r by means of
r
1
=
m
2
m
1
+ m
2
r, r
2
=
m
1
m
1
+ m
2
r . (2)
Then, Ttakes the form
T=
1
2
m
1
m
2
m
1
+ m
2
r
2
and the total Lagrangian as given by equation (1) is:
L =
1
2
(m
1
+m
2
)
R
2
+
1
2
m
1
m
2
m
1
+ m
2
r
2
V (r, r, r, .....) , (3)
where from the reduced mass is dened as
=
m
1
m
2
m
1
+m
2
o
1
=
1
m
1
+
1
m
2
.
Then, the equation (3) can be written as follows
L =
1
2
(m
1
+m
2
)
R
2
+
1
2
r
2
V (r, r, r, .....).
20
From this equation we see that the coordinates
R are cyclic implying that
the center-of-mass is either xed or in uniform motion.
Now, none of the equations of motion for r will contain a term where R or
R will occur. This term is exactly what we will have if a center of force
would have been located in the center of mass with an additional particle
at a distance r away of mass (the reduced mass).
Thus, the motion of two particles around their center of mass, which is due
to a central force can be always reduced to an equivalent problem of a single
body.
2.2 Equations of motion
Now we limit ourselves to conservative central forces for which the potential
is a function of only r, V (r), so that the force is directed along r. Since
in order to solve the problem we only need to tackle a particle of mass
m moving around the xed center of force, we can put the origin of the
reference frame there. As the potential depends only on r, the problem has
spherical symmetry, that is any arbitrary rotation around a xed axis has no
eect on the solution. Therefore, an angular coordinate representing that
rotation should be cyclic providing another considerable simplication to
the problem. Due to the spherical symmetry, the total angular momentum
L = r p
is conserved. Thus, it can be inferred that r is perpendicular to the xed
axis of L. Now, if L = 0 the motion should be along a line passing through
the center of force, since for L = 0 r and r are parallel. This happens only
in the rectilinear motion, and therefore central force motions proceed in one
plane.
By taking the z axis as the direction of L, the motion will take place in the
(x, y) plane. The spherical angular coordinate will have the constant value
/2 and we can go on as foollows. The conservation of the angular momen-
tum provides three independent constants of motion. As a matter of fact,
two of them, expressing the constant direction of the angular momentum,
are used to reduce the problem of three degrees of freedom to only two. The
third coordinate corresponds to the conservation of the modulus of L.
In polar coordinates the Lagrangian is
L =
1
2
m( r
2
+ r
2
2
) V (r) . (4)
21
As we have seen, is a cyclic coordinate whose canonically conjugate mo-
mentum is the angular momentum
p
=
L
= mr
2
,
then, one of the equations of motion will be
p
=
d
dt
(mr
2
) = 0 . (5)
This leads us to
mr
2
= l = cte , (6)
where l is constant modulus of the angular momentum. From equation (5)
one also gets
d
dt
_
r
2
2
_
= 0. (7)
The factor 1/2 is introduced because (r
2
)/2 is the areolar velocity (the area
covered by the position vector per unit of time).
The conservation of the angular momentum is equivalent to saying that the
areolar velocity is constant. This is nothing else than a proof of Keplers
second law of planetary motion: the position vector of a planet covers equal
areas in equal time intervals. However, we stress that the constancy of the
areolar velocity is a property valid for any central force not only for inverse
square ones.
The other Lagrange equation for the r coordinates reads
d
dt
(m r) mr
2
+
V
r
= 0 . (8)
Denoting the force by f(r), we can write this equation as follows
m r mr
2
= f(r) . (9)
Using the equation (6), the last equation can be rewritten as
m r
l
2
mr
3
= f(r). (10)
Recalling now the conservation of the total energy
E = T +V =
1
2
m( r
2
+r
2
2
) +V (r) . (11)
22
we say that E is a constant of motion. This can be derived from the equa-
tions of motion. The equation (10) can be written as follows
m r =
d
dr
_
V (r) +
1
2
l
2
mr
2
_
, (12)
and by multiplying by r both sides, we get
m r r =
d
dt
(
1
2
m r) =
d
dt
_
V (r) +
1
2
l
2
mr
2
_
,
or
d
dt
_
1
2
m r
2
+ V (r) +
1
2
l
2
mr
2
_
= 0 .
Thus
1
2
m r
2
+ V (r) +
1
2
l
2
mr
2
= cte (13)
and since (l
2
/2mr
2
) = (mr
2
/2), the equation (13) is reduced to (11).
Now, let us solve the equations of motion for r and . Taking r from equation
(13), we have
r =
2
2
m
(E V
l
2
2mr
2
) , (14)
or
dt =
dr
2
_
2
m
(E V
l
2
2mr
2
)
. (15)
Let r
0
be the value of r at t = 0. The integral of the two terms of the
equation reads
t =
_
r
r
0
dr
2
_
2
m
(E V
l
2
2mr
2
)
. (16)
This equation gives t as a function of r and of the constants of integration
E, l and r
0
. It can be inverted, at least in a formal way, to give r as a
function of t and of the constants. Once we have r, there is no problem to
get starting from equation (6), that can be written as follows
d =
ldt
mr
2
. (17)
If
0
is the initial value of , then (17) will be
= l
_
t
0
dt
mr
2
(t)
+
0
. (18)
23
Thus, we have already get the equations of motion for the variables r and
.
2.3 The dierential equation of the orbit
A change of our standpoint regarding the approach of real central force prob-
lems prove to be convenient. Till now, solving the problem meant seeking r
and as functions of time and some constants of integration such as E, l,
etc. However, quite often, what we are really looking for is the equation of
the orbit, that is the direct dependence between r and , by eliminating the
time parameter t. In the case of central force problems, this elimination is
particularly simple because t is to be found in the equations of motion only
in the form of a variable with respect to which the derivatives are performed.
Indeed, the equation of motion (6) gives us a denite relationship between
dt and d
ldt = mr
2
d. (19)
The corresponding relationship between its derivatives with respect to t and
is
d
dt
=
l
mr
2
d
d
. (20)
This relationship can be used to convert (10) in a dierential equation for
the orbit. At the same time, one can solve for the equations of motion and
go on to get the orbit equation. For the time being, we follow up the rst
route.
From equation (20) we can write the second derivative with respect to t
d
2
dt
2
=
d
d
l
mr
2
_
d
d
l
mr
2
_
and the Lagrange equation for r, (10), will be
l
r
2
d
d
_
l
mr
2
dr
d
_
l
mr
3
= f(r) . (21)
But
1
r
2
dr
d
=
d(1/r)
d
.
Employing the change of variable u = 1/r, we have
l
2
u
2
m
_
d
2
u
d
2
+u
_
= f
_
1
u
_
. (22)
24
Since
d
du
=
dr
d
d
dr
=
1
u
2
d
dr
,
equation (22) can be written as follows
d
2
u
d
2
+u =
m
l
2
d
du
V
_
1
u
_
. (23)
Any of the equations (22) or (23) is the dierential equation of the orbit if
we know the force f or the potential V . Vice versa, if we know the orbit
equation we can get f or V .
For an arbitrary particular force law, the orbit equation can be obtained
by integrating the equation (22). Since a great deal of work has been done
when solving (10), we are left with the task of eliminating t in the solution
(15) by means of (19),
d =
ldr
mr
2
2
_
2
m
_
E V (r)
l
2
2mr
2
_
, (24)
or
=
_
r
r
0
dr
r
2
2
_
2mE
l
2
2mU
l
2
1
r
2
+
0
. (25)
By the change of variable u = 1/r,
=
0
_
u
u
0
du
2
_
2mE
l
2
2mU
l
2
u
2
, (26)
which is the formal solution for the orbit equation.
2.4 Keplers problem: the case of inverse square force
The inverse square central force lawis the most important of all and therefore
we shall pay more attention to this case. The force and the potential are:
f =
k
r
2
y V =
k
r
. (27)
To integrate the orbit equation we put (23) in (22),
d
2
u
d
2
+u =
mf(1/u)
l
2
u
2
=
mk
l
2
. (28)
25
Now, we perform the change of variable y = u
mk
l
2
, in order that the
dierential equation be written as follows
d
2
y
d
2
+y = 0 ,
possessing the solution
y = Bcos( ) ,
where B and are the corresponding integration constants. The solution in
terms of r is
1
r
=
mk
l
2
[1 + e cos( )] , (29)
where
e = B
l
2
mk
.
We can get the orbit equation from the formal solution (26). Although the
procedure is longer than solving the equation (28), it is nevertheless to do it
since the integration constant e is directly obtained as a function of E and
l.
We write equation (26) as follows
=
_
du
2
_
2mE
l
2
2mU
l
2
u
2
, (30)
where now one deals with a denite integral. Then of (30) is an integration
constant determined through the initial conditions and is not necessarily the
initial angle
0
at t = 0. The solution for this type of integrals is
_
dx
2
_
+x + x
2
=
1
2
arccos
_
+ 2x
2
q
_
, (31)
where
q =
2
4.
In order to apply this type of solutions to the equation (30) we should make
=
2mE
l
2
, =
2mk
l
2
, = 1,
and the discriminant q will be
q =
_
2mk
l
2
_
2
_
1 +
2El
2
mk
2
_
.
26
With these substitutions, (30) is
= arccos
_
_
l
2
u
mk
1
2
_
1 +
2El
2
mk
2
_
_
.
For u 1/r, the resulting orbit equation is
1
r
=
mk
l
2
_
_
1 +
2
1 +
2El
2
mk
2
cos( )
_
_
. (32)
Comparing (32) with the equation (29) we notice that the value of e is:
e =
2
1 +
2El
2
mk
2
. (33)
The type of orbit depends on the value of e according to the following table:
e > 1, E > 0 : hyperbola,
e = 1, E = 0 : parabola,
e < 1, E < 0 : elipse,
e = 0 E =
mk
2
2l
2
: circumference.
2.5 Dispersion by a center of force
From a historical point of view, the interest on central forces was related to
the astronomical problem of planetary motions. However, there is no reason
to consider them only under these circumstances. Another important issue
that one can study within Classical Mechanics is the dispersion of particles
by central forces. Of course, if the particles are of atomic size, we should keep
in mind that the classical formalism may not give the right result because
of quantum eects that begin to be important at those scales. Despite this,
there are classical predictions that continue to be correct. Moreover, the
main concepts of the dispersion phenomena are the same in both Classical
Mechanics and Quantum Mechanics; thus, one can learn this scientic idiom
in the classical picture, usually considered more convenient.
In its one-body formulation, the dispersion problem refers to the action
of the center of force on the trajectories of the coming particles. Let us
consider a uniform beam of particles, (say electrons, protons, or planets and
comets), but all of the same mass and energy impinging on a center of force.
27
We can assume that the force diminishes to zero at large distances. The
incident beam is characterized by its intensity I (also called ux density),
which is the number of particles that pass through per units of time and
normal surface. When one particle comes closer and closer to the center of
force will be attracted or repelled, and its orbit will deviate from the initial
rectilinear path. Once it passed the center of force, the perturbative eects
will diminish such that the orbit tends again to a streight line. In general,
the nal direction of the motion does not coincide with the incident one.
One says that the particle has been dispersed. By denition, the dierential
cross section () is
()d =
dN
I
, (34)
where dN is the number of particles dispersed per unit of time in the element
of solid angle d around the direction. In the case of central forces there
is a high degree of symmetry around the incident beam axis. Therefore, the
element of solid angle can be written
d = 2 sind, (35)
where is the angle between two incident dispersed directions, and is called
the dispersion angle.
For a given arbitrary particle the constants of the orbit and therefore the
degree of dispersion are determined by its energy and angular momentum.
It is convenient to express the latter in terms of a function of energy and
the so-called impact parameter s, which by denition is the distance from
the center of force to the straight suport line of the incident velocity. If u
0
is the incident velocity of the particle, we have
l = mu
0
s = s
2
2mE. (36)
Once E and s are xed, the angle of dispersion is uniquely determined.
For the time being, we suppose that dierent values of s cannot lead to
the same dispersion angle. Therefore, the number of dispersed particles in
the element of solid angle d between and + d should be equal to
the number of incident particles whose impact parameter ranges within the
corresponding s and s +ds:
2Is [ds[ = 2()I sin[d[ . (37)
In the equation (37) we have introduced absolute values because while the
number of particles is always positive, s and can vary in opposite direc-
tions. If we consider s as a function of the energy and the corresponding
28
dispersion angle,
s = s(, E),
the dependence of the cross section of will be given by
() =
s
sin
ds
d
. (38)
From the orbit equation (25), one can obtain directly a formal expression
for the dispersion angle. In addition, for the sake of simplicity, we tackle
the case of a pure repulsive dispersion. Since the orbit should be symmetric
with respect to the direction of the periapsis, the dispersion angle is
= 2 , (39)
where is the angle between the direction of the incident asymptote and
the direction of the periapsis. In turn, can be obtained from the equation
(25) by making r
0
= when
0
= (incident direction). Thus, =
when r = r
m
, the closest distance of the particle to the center of force.
Then, one can easily obtain
=
_
rm
dr
r
2
2
_
2mE
l
2
2mV
l
2
1
r
2
. (40)
Expressing l as a function of the impact parameter s (eq. (36)), the result
is
= 2
_
rm
sdr
r
2
_
r
2
_
1
V (r)
E
_
s
2
, (41)
or
= 2
_
um
0
sdu
2
_
1
v(u)
E
s
2
u
2
. (42)
The equations (41) and (42) are used rarely, as they do not enter in a direct
way in the numerical calculation of the dispersion angle. However, when an
analytic expression for the orbits is available, one can often get, merely by
simple inspection, a relationship between and s.
EXAMPLE:
This example is very important from the historical point of view. It refers
to the repulsive dispersion of charged particles in a Coulomb eld. The eld
29
is produced by a xed charge Ze and acts on incident particles of charge
Ze; therefore, the force can be written as follows
f =
ZZe
2
r
2
,
that is, one deals with a repulsive inverse square force. The constant is
k = ZZe
2
. (43)
The energy E is positive implying a hyperbolic orbit of eccentricity
=
2
1 +
2El
2
m(ZZe
2
)
2
=
2
1 +
_
2Es
ZZe
2
_
2
, (44)
where we have taken into account the equation (36). If the angle is taken
to be , then from the equation (29) we come to the conclusion that the
periapse corresponds to = 0 and the orbit equation reads
1
r
=
mZZe
2
l
2
[ cos 1] . (45)
The direction of the incident asymptote is thus determined by the condi-
tion r :
cos =
1
,
that is, according to equation (39),
sin
2
=
1
.
Then,
cot
2
2
=
2
1,
and by means of equation (44)
cot
2
=
2Es
ZZe
2
.
The functional relationship between the impact parameter and the disper-
sion angle will be
s =
ZZe
2
2E
cot
2
, (46)
30
and by eecting the transformation required by the equation (38) we nd
that () is given by
() =
1
4
_
ZZe
2
2E
_
2
csc
4
2
. (47)
The equation (47) gives the famous Rutherford scattering cross section de-
rived by him for the dispersion of particles on atomic nuclei. In the non-
relativistic limit, the same result is provided by the quantum mechanical
calculations.
The concept of total cross section
T
is very important in atomic physics.
Its denition is
T
=
_
4
()d = 2
_
0
()d .
However, if we calculate the total cross section for the Coulombian dispersion
by substituting the equation (47) in the denition above we get an innite
result. The physical reason is easy to see. According to the denition, the
total cross section is the number of particles per unit of incident intensity
that are dispersed in all directions. The Coulombian eld is an example of
long-range force; its eects are still present in the innite distance limit. The
small deviation limit is valid only for particles of large impact parameter.
Therefore, for an incident beam of innite lateral extension all the particles
will be dispersed and should be included in the total cross section. It is clear
that the innite value of
T
is not a special property of the Coulombian eld
and occurs for any type of long-range eld.
Further reading
L.S. Brown, Forces giving no orbit precession, Am. J. Phys. 46, 930 (1978)
H. Goldstein, More on the prehistory of the Laplace-Runge-Lenz vector, Am.
J. Phys. 44, 1123 (1976)
31
3. THE RIGID BODY
Forward: Due to its particular features, the study of the motion of the
rigid body has generated several interesting mathematical techniques and
methods. In this chapter, we briey present the basic rigid body concepts.
CONTENTS:
3.1 Denition
3.2 Degrees of freedom
3.3 Tensor of inertia (with example)
3.4 Angular momentum
3.5 Principal axes of inertia (with example)
3.6 The theorem of parallel axes (with 2 examples)
3.7 Dynamics of the rigid body (with example)
3.8 Symmetrical top free of torques
3.9 Euler angles
3.10 Symmetrical top with a xed point
32
3.1 Denition
A rigid body (RB) is dened as a system of particles whose relative distances
are forced to stay constant during the motion.
3.2 Degrees of freedom
In order to describe the general motion of a RB in the three-dimensional
space one needs six variables, for example the three coordinates of the center
of mass measured with respect to an inertial frame and three angles for
labeling the orientation of the body in space (or of a xed system within
the body with the origin in the center of mass). In other words, in the
three-dimensional space the RB can be described by at most six degrees of
freedom.
The number of degrees of freedom may be less when the rigid body is sub-
jected to various conditions as follows:
If the RB rotates around a single axis there is only one degree of
freedom (one angle).
If the RB moves in a plane, its motion can be described by ve degrees
of freedom (two coordinates and three angles).
3.3 Tensor of inertia.
We consider a body made of N particles of masses m
, = 1, 2, 3..., N.
If the body rotates at angular velocity around a xed point in the body
and this point, in turn, moves at velocity v with respect to a xed inertial
system, then the velocity of the th particle w.r.t. the inertial system is
given by
v
= v + r
. (1)
The kinetic energy of the th particle is
T
=
1
2
m
v
2
, (2)
where
v
2
= v
= (v + r
) (v + r
) (3)
= v v + 2v ( r
) + ( r
) ( r
)
= v
2
+ 2v( r
) + ( r
)
2
. (4)
33
Then the total energy is
T =
1
2
m
v
2
+
[v ( r
)] +
+
1
2
( r
)
2
;
T =
1
2
Mv
2
+v [
] +
1
2
( r
)
2
.
If the origin is xed to the solid body, we can take it in the center of mass.
Thus,
R =
M
= 0,
and therefore we get
T =
1
2
Mv
2
+
1
2
( r
)
2
(5)
T = T
trans
+ T
rot
(6)
where
T
trans
=
1
2
v
2
=
1
2
Mv
2
(7)
T
rot
=
1
2
( r
)
2
. (8)
In Eq. (8) we use the vectorial identity
(AB)
2
= A
2
B
2
(A B)
2
(9)
to get the following form of the equation
T
rot
=
1
2
2
r
2
( r
)
2
_
,
which in terms of the components of and r
= (
1
,
2
,
3
) and r
= (x
1
, x
2
, x
3
)
can be written as follows
T
rot
=
1
2
_
_
2
i
__
k
x
2
k
_
i
x
i
_
_
j
x
j
__
.
Now, we introduce
i
=
j
ij
j
34
T
rot
=
1
2
ij
m
ij
_
k
x
2
k
_
j
x
i
x
j
_
(10)
T
rot
=
1
2
ij
ij
k
x
2
k
x
i
x
j
_
. (11)
We can write T
rot
as follows
T
rot
=
1
2
ij
I
ij
j
(12)
where
I
ij
=
ij
k
x
2
k
x
i
x
j
_
. (13)
The nine quantities I
ij
are the components of a new mathematical entity, de-
noted by I
ij
and called tensor of inertia. It can be written in a convenient
way as a (3 3) matrix
I
ij
=
_
_
_
I
11
I
12
I
13
I
21
I
22
I
23
I
31
I
32
I
33
_
_
_ =
=
_
_
_
(x
2
2
+ x
2
3
)
x
1
x
2
x
1
x
3
x
2
x
1
(x
2
1
+ x
2
3
)
x
2
x
3
x
3
x
1
x
3
x
2
(x
2
1
+ x
2
2
)
_
_
_.
(14)
We note that I
ij
= I
ji
, and therefore I
ij
is a symmetric tensor, implying
that only six of the components are independent. The diagonal elements of
I
ij
are called moments of inertia with respect to the axes of coordinates,
whereas the negatives of the nondiagonal elements are called the products
of inertia. For a continuous distribution of mass, of density (r), I
ij
is
written in the following way
I
ij
=
_
V
(r)
_
ij
k
x
2
k
x
i
x
j
_
dV. (15)
EXAMPLE:
Find the elements I
ij
of the tensor of inertia I
ij
for a cube of uniform
density of side b, mass M, with one corner placed at the origin.
I
11
=
_
V
_
x
2
1
+x
2
2
+x
2
3
x
1
x
1
_
dx
1
dx
2
dx
3
=
b
_
0
b
_
0
b
_
0
(x
2
2
+x
2
3
)dx
1
dx
2
dx
3
.
35
The result of the three-dimensional integral is I
11
=
2
3
(b
3
)
2
=
2
3
Mb
2
.
I
12
=
_
V
(x
1
x
2
)dV =
b
_
0
b
_
0
b
_
0
(x
1
x
2
)dx
1
dx
2
dx
3
=
1
4
b
5
=
1
4
Mb
2
.
We see that all the other integrals are equal, so that
I
11
= I
22
= I
33
=
2
3
Mb
2
I
ij
i=j
=
1
4
Mb
2
,
leading to the following form of the matrix
I
ij
=
_
_
_
2
3
Mb
2
1
4
Mb
2
1
4
Mb
2
1
4
Mb
2 2
3
Mb
2
1
4
Mb
2
1
4
Mb
2
1
4
Mb
2 2
3
Mb
2
_
_
_.
3.4 Angular Momentum
The angular momentum of a RB made of N particles of masses m
is given
by
L =
, (16)
where
p
= m
= m
( r
) . (17)
Substituting (17) in (16), we get
L =
( r
) .
Employing the vectorial identity
A(B A) = (A A)B(A B)A = A
2
B(A B)A ,
leads to
L =
(r
2
( r
).
Considering the ith component of the vector L
L
i
=
k
x
2
k
_
x
i
_
j
x
j
j
_
,
36
and introducing the equation
i
=
j
ij
,
we get
L
i
=
j
ij
k
x
2
k
_
j
x
j
x
j
j
_
(18)
=
j
j
ij
_
k
x
2
k
x
i
x
j
_
(19)
=
j
j
ij
k
x
2
k
x
i
x
j
_
. (20)
Comparing with the equation (13) leads to
L
i
=
j
I
ij
j
. (21)
This equation can also be written in the form
L =I
ij
, (22)
or
_
_
_
L
1
L
2
L
3
_
_
_ =
_
_
_
I
11
I
12
I
13
I
21
I
22
I
23
I
31
I
32
I
33
_
_
_
_
_
_
3
_
_
_. (23)
The rotational kinetic energy, T
rot
, can be related to the angular momentum
as follows: rst, multiply the equation ( 21) by
1
2
i
1
2
L
i
=
1
2
j
I
ij
j
, (24)
and next summing over all the i indices, gives
i
1
2
L
i
i
=
1
2
ij
I
ij
j
.
Compararing this equation with (12), we see that the second term is T
rot
.
Therefore
T
rot
=
I
1
2
L
i
i
=
1
2
L . (25)
Now, we substitute (22) in the equation (25), getting the relationship be-
tween T
rot
and the tensor of inertia
T
rot
=
1
2
I
ij
. (26)
37
3.5 Principal axes of inertia
Taking the tensor of inertia I
ij
diagonal, that is I
ij
= I
i
ij
, the rotational
kinetic energy and the angular momentum are expressed as follows
T
rot
=
1
2
ij
I
ij
j
=
1
2
ij
ij
I
i
j
T
rot
=
1
2
i
I
i
2
i
(27)
and
L
i
=
j
I
ij
j
=
j
ij
I
i
j
= I
i
i
L = I. (28)
To seek a diagonal form of I
ij
is equivalent to nding a new system of
three axes for which the kinetic energy and the angular momentum take the
form given by (27) and (28). In this case the axes are called principal axes
of inertia. That means that given an inertial reference system within the
body, we can pass from it to the principal axes by a particular orthogonal
transformation, which is called transformation to the principal axes.
Making equal the components of (22) and (28), we have
L
1
= I
1
= I
11
1
+I
12
2
+I
13
3
(29)
L
2
= I
2
= I
21
1
+I
22
2
+I
23
3
(30)
L
3
= I
3
= I
31
1
+I
32
2
+I
33
3
. (31)
This is a system of equations that can be rewritten as
(I
11
I)
1
+I
12
2
+I
13
3
= 0 (32)
I
21
1
+ (I
22
I)
2
+I
23
3
= 0
I
31
1
+ I
32
2
+ (I
33
I)
3
= 0 .
38
To get nontrivial solutions, the determinant of the system should be zero
(I
11
I)
1
I
12
2
I
13
3
I
21
1
(I
22
I)
2
I
23
3
I
31
1
I
32
2
(I
33
I)
3
= 0 . (33)
This determinant leads to a polynomial of third order in I, known as the
characteristic polynomial. The equation (33) is called the secular equation or
characteristic equation. In practice, the principal moments of inertia, being
the eigenvalues of I, are obtained as solutions of the secular equation.
EXAMPLE:
Determine the principal axes of inertia for the cube of the previous example.
Substituting the values obtained in the previous example in the equation
(33) we get:
_
_
_
(
2
3
I)
1
4
1
4
1
4
(
2
3
I)
1
4
1
4
1
4
(
2
3
I)
_
_
_
= 0 ,
where = Mb
2
. Thus, the characteristic equation will be
_
11
12
I
__
11
12
I
__
1
6
I
_
= 0 .
The solutions, i.e., the principal moments of inertia are:
I
1
=
1
6
, I
2
= I
3
=
11
12
,
whose corresponding eigenvalues are given by
I =
1
6
1
2
3
_
_
_
1
1
1
_
_
_, I
2
, I
3
=
11
12
1
2
2
_
_
_
_
_
1
1
0
_
_
_,
_
_
_
1
0
1
_
_
_
_
_
.
The matrix that diagonalizes I
ij
is:
=
2
_
1
3
_
_
_
_
_
1
2
_
3
2
2
_
3
2
1
2
_
3
2
0
1 0
2
_
3
2
_
_
_
_
_
.
39
The diagonalized I
ij
will be
I
ij
diag
= ()
I
ij
=
_
_
_
1
6
0 0
0
11
12
0
0 0
11
12
_
_
_.
3.6 The theorem of parallel axes
We suppose that the system x
1
, x
2
, x
3
has the origin in the center of mass
of the RB. A second system X
1
, X
2
, X
3
, has the origin in another position
w.r.t. the rst system. The only imposed condition on them is to be parallel.
We dene the vectors r = (x
1
, x
2
, x
3
), R = (X
1
, X
2
, X
3
) y a = (a
1
, a
2
, a
3
)
in such a way that R = r +a, or in component form
X
i
= x
i
+a
i
. (34)
Let J
ij
be the components of the tensor of inertia w.r.t. the system X
1
X
2
X
3
,
J
ij
=
ij
k
X
2
k
X
i
X
j
_
. (35)
We substitute (34) in (35),
J
ij
=
ij
k
(x
k
+ a
k
)
2
(x
i
+a
i
)(x
j
+ a
j
)
_
=
_
ij
k
(x
k
)
2
x
i
x
j
_
+
ij
k
a
2
k
a
i
a
j
_
(36)
+[
k
2a
k
ij
(
x
k
) a
j
(
x
j
) a
i
(
x
i
)] .
Since the center of mass coordinate is dened as
x =
M
we take into account that we have already set the origin in the center of
mass, i.e.,
( x
1
, x
2
, x
3
) = (0, 0, 0) .
Now, if we also compare the rst term in (36) with the equation (13), we
have
J
ij
= I
ij
+M(a
2
ij
a
i
a
j
) (37)
40
and therefore the elements of the tensor of inertia I
ij
for the center of mass
system will be given by:
I
ij
= J
ij
M(
ij
a
2
a
i
a
j
) . (38)
This is known as the theorem of the parallel axes.
EXAMPLE:
Find I
ij
for the previous cube w.r.t. a reference system parallel to the system
in the rst example and with the origin in the center of mass.
We already know from the previous example that:
J
ij
=
_
_
_
2
3
1
4
1
4
1
4
2
3
1
4
1
4
1
4
2
3
_
_
_ .
Now, since the vector a = (
b
2
,
b
2
,
b
2
) and a
2
=
3
4
b
2
, we can use the equation
(38) and the fact that = Mb
2
to get,
I
11
= J
11
M(a
2
a
2
1
) =
1
6
Mb
2
(39)
I
22
= J
22
M(a
2
a
2
2
) =
1
6
Mb
2
(40)
I
33
= J
33
M(a
2
a
2
3
) =
1
6
Mb
2
(41)
I
12
= J
12
M(a
1
a
2
) = 0 (42)
I
12
= I
13
= I
23
= 0 . (43)
Therefore
I =
_
_
_
1
6
Mb
2
0 0
0
1
6
Mb
2
0
0 0
1
6
Mb
2
_
_
_.
EXAMPLE:
We consider the case for which the vector a = (0,
b
2
,
b
2
) and a
2
=
b
2
2
. Then,
the new tensor of inertia will be:
I
11
= J
11
M(a
2
a
2
1
) =
_
2
3
Mb
2
_
M
_
b
2
2
0
_
=
1
6
Mb
2
(44)
I
22
= J
22
M(a
2
a
2
2
) =
_
2
3
Mb
2
_
M
_
b
2
2
b
2
4
_
=
5
12
Mb
2
(45)
I
33
= J
33
M(a
2
a
2
3
) =
_
2
3
Mb
2
_
M
_
b
2
2
b
2
4
_
=
5
12
Mb
2
(46)
I
12
= J
12
M(a
1
a
2
) =
_
1
4
Mb
2
_
M(0) =
1
4
Mb
2
(47)
41
I
13
= J
13
M(a
1
a
3
) =
_
1
4
Mb
2
_
M(0) =
1
4
Mb
2
(48)
I
23
= J
23
M(a
2
a
3
) =
_
1
4
Mb
2
_
M(
1
4
Mb
2
) = 0 . (49)
It follows that I
ij
is equal to:
I
ij
=
_
_
_
1
6
Mb
2
1
4
Mb
2
1
4
Mb
2
1
4
Mb
2 5
12
Mb
2
0
1
4
Mb
2
0
5
12
Mb
2
_
_
_ .
3.7 The dynamics of the rigid body
The rate of change in time of the angular momentum L is given by:
_
dL
dt
_
inertial
= N
(e)
. (50)
For the description w.r.t. the system xed to the body we have to use the
operator identity
_
d
dt
_
inertial
=
_
d
dt
_
body
+ . (51)
Applying this operator to the equation (50)
_
dL
dt
_
inertial
=
_
dL
dt
_
body
+ L. (52)
Then, instead of (50) we shall have
_
dL
dt
_
body
+ L = N. (53)
Now we project the equation (53) onto the principal axes of inertia, that we
call (x
1
, x
2
, x
3
); then T
rot
and L take by far simpler forms, e.g.,
L
i
= I
i
i
. (54)
The ith component of (53) is
dL
i
dt
+
ijk
j
L
k
= N
i
. (55)
42
Now projecting onto the principal axes of inertia and using the equation
(54), one can put (55) in the form:
I
i
d
i
dt
+
ijk
k
I
k
= N
i
(56)
since the principal elements of inertia are independent of time. Thus, we
obtain the following system of equations known as Eulers equations
I
1
1
+
2
3
(I
2
I
3
) = N
1
(57)
I
2
2
+
3
1
(I
3
I
1
) = N
2
I
3
3
+
1
2
(I
1
I
2
) = N
3
.
EXAMPLE:
For the rolling and sliding of a billiard ball, prove that after a horizontal
kick the ball slips a distance
x
1
=
12u
2
0
49g
,
where u
0
is the initial velocity. Then, it starts rolling without gliding at the
time
t
1
=
2u
0
7g
.
SOLUTION: When the impulsive force stops, the initial conditions are
x
0
= 0, x
0
= u
0
= 0,
= 0 .
The friction force is given by
F
f
= ge
1
,
and the equation of motion reads
x = gM. (58)
The equation for L is
dL
3
dt
= I
3
= N
3
(59)
43
where I
3
is
I
3
=
_
(r)
_
x
2
1
x
2
2
_
dx
1
dx
2
dx
3
=
2
5
Ma
2
and
N
3
= F
f
a = Mga .
Substituting I
3
and N
3
in (59), one gets
a
=
5
2
g. (60)
Integrating once both (58) and (60), we get
x = gt + C
1
(61)
a
=
5
2
gt + C
2
. (62)
To these equations we apply the initial conditions to put them in the form
x(t) = gt + u
0
(63)
a
(t) =
5
2
gt. (64)
The condition of pure rolling (no friction) is
x(t) = a
(t). (65)
From (64) and (65) evaluated at t
1
we get
5
2
gt
1
= gt
1
+u
0
t
1
=
2u
0
7g
. (66)
Now we integrate (63) once again and applying the initial conditions we get
x(t) = g
t
2
2
+u
0
t . (67)
By evaluating (67) and (63) at time t
1
we are led to
x =
12u
2
49g
44
x =
5
7
u
0
.
3.8 Symmetrical top free of torques
A symmetric top is any solid of revolution. If the moments of inertia are
I
1
= I
2
= I
3
spherical top
I
1
= I
2
,= I
3
symmetric top
I
1
,= I
2
,= I
3
asymmetric top.
Let us consider the symmetric top I
1
= I
2
,= I
3
. In this case the axis X
3
is the axis of symmetry. The Euler equations projected onto the principal
axes of inertia read
I
1
1
+
2
3
(I
2
I
3
) = N
1
(68)
I
2
2
+
3
1
(I
3
I
1
) = N
2
(69)
I
3
3
+
1
2
(I
1
I
2
) = N
3
. (70)
Since the system we consider here is free of torques
N
1
= N
2
= N
3
= 0 , (71)
we use I
1
= I
2
in (71) to get
I
1
1
+
2
3
(I
2
I
3
) = 0 (72)
I
2
2
+
3
1
(I
3
I
1
) = 0 (73)
I
3
3
= 0. (74)
The equation (74) implies that
3
= const.
The equations (72) and (73) are rewritten as follows:
1
=
2
where =
3
_
I
3
I
1
I
1
_
(75)
2
=
1
. (76)
45
Multiplying (76) by i and summing it to (75), we have
(
1
+ i
2
) = (
2
i
1
)
(
1
+ i
2
) = i(
1
+i
2
).
If we write (t) =
1
(t) + i
2
(t), then
(t) i(t) = 0 .
The solution is
(t) = Aexp(it) .
This implies
(
1
+i
2
) = Acos(t) + i sin(t) .
Thus,
1
= Acos(t) (77)
2
= Asin(t). (78)
The modulus of the vector does not change in time
= [[[[ =
2
1
+
2
+
3
=
2
_
A
2
+
2
3
= const .
This vector performs a precessional motion of precession frequency given
by
=
3
_
I
3
I
1
I
1
_
.
Moreover, we notice that is constant.
If we denote by the angle between and X
3
the equations (77) and (78)
take the form
1
= sin cos(t)
2
= sin sin(t)
1
= cos ,
where A = sin.
For a attened body of revolution I
1
= I
2
= I
12
and I
3
> I
1
. For example,
for the case of the earth
=
3
_
I
3
I
12
I
12
_
3
305
.
46
The observations point to a mean value of fourteen months 450 days.
(This is due to the fact that the earth is not strictly a RB; there is also an
internal liquid structure).
3.9 Euler angles
As we already know, a rotation can be described by a rotation matrix by
means of the equation
x = x. (79)
x represents the set of axes of the system rotated w.r.t. the system whose
axes are represented by x. The rotation can be accomplished through
a set of partial rotations =
1
2
...
n
. There are many possibilities to
choose these s. One of them is the set of angles , and called Euler
angles. The partial rotations are in this case the following:
A rotation around the X
3
axis of angle (in the positive trigonometric
sense). The corresponding matrix is:
=
_
_
_
cos sin 0
sin cos 0
0 0 1
_
_
_ .
A rotation of angle around the X
1
axis (positive sense).The associ-
ated matrix is:
=
_
_
_
1 0 0
0 cos sin
0 sin cos
_
_
_.
A rotation of angle around the X
3
axis (positive sense); the assoiated
matrix is:
=
_
_
_
cos sin 0
sin cos 0
0 0 1
_
_
_.
The full transformation of the system of axes X
1
, X
2
, X
3
to the system
of axes X
1
, X
2
, X
3
is given by (79), where
=
.
Doing the product of matrices, we get
11
= cos cos cos sinsin
47
21
= sincos cos sin cos
31
= sin sin
12
= cos sin + cos cos sin
22
= sin sin + cos cos sin
32
= sincos
13
= sin cos
23
= cos sin
33
= cos
where
=
_
_
_
11
12
13
21
22
23
31
32
33
_
_
_.
Now, we take into account that:
is along the X
3
(xed) axis.
is along the so-called line of nodes.
is along the X
3
axis (of the body).
This allows to write the three components of each of the three vectors in the
system X
1
, X
2
, X
3
as follows:
1
=
sin sin,
1
=
cos
1
= 0
2
=
sin cos ,
2
=
sin
2
= 0
3
=
cos ,
3
= 0
3
= .
Then,
=
+
+
=
__
1
+
1
+
1
_
,
_
2
+
2
+
2
_
,
_
3
+
3
+
3
__
.
Thus, we are led to the following components of :
1
=
sin sin +
cos
2
=
sin cos
sin
3
=
cos + .
48
3.10 Symmetrical top with a xed point
As a more complicated example of the methods used to describe the dynam-
ics of the rigid body, we shall consider the motion of a symmetric body in a
unifom gravitational eld when a point of the axis of symmetry is xed in
the space.
The axis of symmetry is of course one of the principal axes, and we shall
take it as the z axis of the body-xed reference system. Since there is a
xed point, the conguration of the top will be determined by the three
Euler angles: measuring the deviation of z from the vertical, , giving the
azimuth of the top w.r.t. the vertical, and , which is the rotation angle of
the top w.r.t. its proper z. The distance from the center of gravity to the
xed point will be denoted by l. To get a solution to the motion of the top
we shall use the method of Lagrange instead of the Euler equations.
The kinetic energy is:
T =
1
2
I
1
(
2
1
+
2
2
) +
1
2
I
3
2
3
,
or, in terms of the Euler angles:
T =
1
2
I
1
(
2
sin
2
+
2
) +
1
2
I
3
(
cos + )
2
.
According to an elementary theorem, in a constant gravitational eld the
potential energy of a body is the same with that of a material point of equal
mass concentrated in its center of mass. A formal proof is as follows. The
potential energy of the body is the sum of the potential energies of all its
particles:
V = m
i
r
i
g , (80)
where g is the constant acceleration of gravity. According to the denition
of the center of mass, this is equivalent to
V = MR
i
g, (81)
thus proving the theorem. The potential energy is a function of the Euler
angles:
V = Mgl cos , (82)
and the Lagrangian will be
L =
1
2
I
1
(
2
sin
2
+
2
) +
1
2
I
3
(
cos + )
2
Mgl cos . (83)
49
We note that and are cyclic coordinates, and therefore p
and p
are
constants of motion.
p
=
L
= I
3
( +
cos ) = const (84)
and
p
=
L
= I
1
sin
2
+I
3
(
cos
2
+ cos ) = const. (85)
From the equation (84) we get
=
p
I
3
cos
I
3
, (86)
that we substitute in (85)
p
=
L
= I
1
sin
2
+ I
3
(
cos
2
+
p
I
3
cos
I
3
cos ) = const.
p
= I
1
sin
2
+ p
cos ,
where from we get
=
p
cos
I
1
sin
2
. (87)
Substituting it in (86) one gets
=
p
I
3
cos
I
1
sin
2
cos . (88)
Now, since the system is conservative, another integral of motion is the
energy
E = T + V =
1
2
I
1
(
2
sin
2
+
2
) +
1
2
I
3
(
cos + )
2
+Mgl cos .
The quantity I
3
3
= p
3
we get
I
3
p
2
3
= p
2
3
I
2
3
3
3
= p
2
3
1
2
I
3
2
3
=
1
2
p
2
I
3
.
50
The quantity
1
2
I
3
2
3
is a constant. Therefore, we can dene the quantity
E = E
1
2
I
3
2
3
= const.
=
1
2
I
1
2
+
1
2
I
1
2
sin
2
+Mgl cos ,
wherefrom we can identify
V () =
1
2
2
sin
2
+Mgl cos
V () =
1
2
I
1
_
p
cos
I
1
sin
2
_
2
sin
2
+ Mgl cos . (89)
Thus, Eis:
E=
1
2
I
1
2
+ V () .
From this equation we get
d
dt
=
_
2
I
1
(EV ())
_
1/2
, which leads to
t() =
_
d
2
_
_
2
I
1
_
(EV ())
. (90)
Performing the integral in (90) one gets t = f(), and therefore, in principle,
one can get (t). Then, (t) is replaced by
and (in eqs. (87) and (88))
and integrating them we can obtain the complete solution of the problem.
References
H. Goldstein, Classical Mechanics, (Addison-Wesley, 1992).
L. D. Landau & E. M. Lifshitz, Mechanics, (Pergammon, 1976).
J. B. Marion & S.T. Thornton, Classical Dynamics of Particles and
Systems, (Harcourt Brace, 1995).
W. Wrigley & W.M. Hollister, The Gyroscope: Theory and application,
Science 149, 713 (Aug. 13, 1965).
51
4. SMALL OSCILLATIONS
Forward: A familiar type of motion in mechanical and many other systems
are the small oscillations (vibrations). They can be met as atomic and
molecular vibrations, electric circuits, acoustics, and so on. In general, any
motion in the neighborhood of stable equilibria is vibrational.
CONTENTS:
4.1 THE SIMPLE HARMONIC OSCILLATOR
4.2 FORCED HARMONIC OSCILLATOR
4.3 DAMPED HARMONIC OSCILLATORS
4.4 NORMAL MODES
4.5 PARAMETRIC RESONANCE
52
4.1 THE SIMPLE HARMONIC OSCILLATOR
A system is at stable equilibrium when its potential energy U(q) is at min-
imum; when the system is slightly displaced from the equilibrium position,
a force dU/dq occurs which acts to restore the equilibrium. Let q
0
be the
value of the generalized coordinate corresponding to the equilibrium posi-
tion. Expanding U(q) U(q
0
) in a Taylor series of qq
0
for small deviations
from the equilibrium
U(q) U(q
0
)
=
1
2
k(q q
0
)
2
,
donde:
U
q
= 0
U(q) = 0 ,
This means that there are no external forces acting on the system and the
zero has been chosen at the equilibrium position; moreover, higher-order
terms have been neglected. The coecient k represents the value of the
second derivative of U(q) for q=q
0
. For simplicity reasons we denote
x = q q
0
for which the potential energy can be written as:
U(x) =
1
2
kx
2
. (1)
For simplicity reasons we denote
x = q q
0
for which the potential energy can be written as:
U(x) =
1
2
kx
2
. (2)
The kinetic energy of a system is
T =
1
2
m
x
2
, (3)
and using (2) and (3) we get the Lagrangian of a system performing linear
oscillations (so-called linear oscillator):
53
L =
1
2
m
x
2
1
2
kx
2
. (4)
The equation of motion corresponding to this L is:
m
x +kx = 0 ,
or
x +w
2
x = 0 , (5)
where w
2
=
_
k/m. This dierential equation has two independent solu-
tions: cos wt and sinwt, from which one can form the general solution:
x = c
1
cos wt +c
2
sin wt , (6)
or, we can also write the solution in the form:
x = a cos(wt +) . (7)
Since cos(wt+) = cos wt cos sinwtsin, by comparing with (6), one can
see that the arbitrary constants a and are related to c
1
and c
2
as follows:
a =
_
(c
2
1
+c
2
2
), y tan = c
1
/c
2
.
Thus, any system in the neighborhood of the stable equilibrium position
performs harmonic oscillatory motion. The a coecient in (7) is the am-
plitude of the oscillations, whereas the argument of the cosine function is
the phase of the harmonic oscillation; is the initial value of the phase,
which depends on the chosen origin of time. The quantity w is the angular
frequency of the oscillations, which does not depend on the initial conditions
of the system, being a proper characteristic of the harmonic oscillations.
Quite often the solution is expressed as the real part of a complex quantity
x = Re [Aexp(iwt)]
where A is the complex amplitude, whose modulus gives the ordinary am-
plitude:
A = a exp(i) .
The energy of a system in small oscillatory motion is:
54
E =
1
2
m
x
2
+
1
2
kx
2
,
or by substituting (7)
E =
1
2
mw
2
a
2
.
Now, we consider the case of n degrees of freedom. In this case, taking the
sum of exterior forces as zero, the generalized force will be given by
Q
i
=
U
q
i
= 0 . (8)
Repeating the procedure for the case of a single degree of freedom, we expand
the potential energy in Taylor series taking the minimum of the potential
energy at q
i
= q
i0
. Introducing small oscillation coordinates
x
i
= q
i
q
i0
,
we can write the series as follows
U(q
1
, q
2
, ..., q
n
) = U(q
10
, q
20
, ..., q
n0
)+
_
U
q
i
_
0
x
i
+
1
2!
_
2
U
q
i
q
j
_
0
x
i
x
j
+....
(9)
Under the same considerations as given for (2), we obtain:
U(q
1
, q
2
, ..., q
n
) = U =
1
2
i,j
k
ij
x
i
x
j
. (10)
From (9) one notes that k
ij
= k
ji
, i.e., they are symmetric w.r.t. their
subindices.
Let us look now to the kinetic energy, which, in general, is of the form
1
2
a
ij
(q)
x
i
x
j
,
where the a
ij
are functions of the coordinates only. Denoting them by
a
ij
= m
ij
the kinetic energy will be
T =
1
2
i,j
m
ij
x
i
x
j
. (11)
55
We can pass now to the Lagrangian for the system of n degrees of freedom
L = T U =
1
2
i,j
(m
ij
x
i
x
j
k
ij
x
i
x
j
) . (12)
This Lagrangian leads to the following set of simultaneous dierential equa-
tions of motion
d
dt
L
x
i
L
x
i
= 0 (13)
or
(m
ij
x
j
+k
ij
x
j
) = 0 . (14)
This is a linear system of homogeneous equations, which can be considered
as the n components of the matricial equation
(M)(
X) + (K)(X) = 0 , (15)
where the matrices are dened by:
(M) =
_
_
_
_
_
_
m
11
m
12
... m
1n
m
21
m
22
... m
2n
.
.
.
.
.
.
m
n1
m
n2
... m
nn
_
_
_
_
_
_
(16)
(K) =
_
_
_
_
_
_
k
11
k
12
... k
1n
k
21
k
22
... k
2n
.
.
.
.
.
.
k
n1
k
n2
... k
nn
_
_
_
_
_
_
(17)
(
X) =
d
2
dt
2
_
_
_
_
_
_
x
1
x
2
.
.
.
x
n
_
_
_
_
_
_
(18)
(X) =
_
_
_
_
_
_
x
1
x
2
.
.
.
x
n
_
_
_
_
_
_
. (19)
56
Similarly to the one dof system, we look for n unknown functions x
j
(t) of
the form
x
j
= A
j
exp(iwt) , (20)
where A
j
are constants to be determined. Substituting (20) in (14) and
dividing by exp(iwt), one gets a linear system of algebraic and homogeneous
equations, which should be fullled by A
j
.
j
(w
2
m
ik
+k
ik
)A
k
= 0 . (21)
This system has nponzero solutions if the determinant of its coecients is
zero.
k
ij
w
2
m
ij
2
= 0 . (22)
This is the characteristic equation of order n w.r.t. w
2
. In general, it has n
dierent real and positive solutions w
are called
proper frequencies of the system. Multiplying by A
i
and summing over i
one gets
j
(w
2
m
ij
+k
ij
)A
i
A
j
= 0 ,
where from
w
2
=
k
ij
A
i
A
i
/
m
ij
A
i
A
i
.
Since the coecients k
ij
and m
ij
are real and symmetric, the quadratic forms
of the numerator and denominator are real, and being essentially positive
one concludes that w
2
are equally positive.
EXAMPLE
As an example we model the equations of motion of a double pendulum.
The potential energy of this system with two degrees of freedom is
U = m
1
gl
1
(1 cos
1
) +m
2
gl
1
(1 cos
1
) +m
2
gl
2
(1 cos
2
) .
Applying (9), one gets
U =
1
2
(m
1
+m
2
)gl
1
2
1
+
1
2
m
2
gl
2
2
2
.
57
Comparing with (10), we identify
k
11
= (m
1
+ m
2
)l
2
1
k
12
= k
21
= 0
k
22
= m
2
gl
2
.
For the kinetic energy one gets
T =
1
2
(m
1
+m
2
)l
2
1
.
2
1
+
1
2
m
2
l
2
2
.
2
2
+m
2
l
1
l
2
.
1
.
2
.
Identifying terms from the comparison with (11) we nd
m
11
= (m
1
+m
2
)l
2
1
m
12
= m
21
= m
2
l
1
l
2
m
22
= m
2
l
2
2
.
Substituting the energies in (12) one obtains the Lagrangian for the double
pendulum oscillator and as the nal result the equations of motion for this
case:
_
m
11
m
12
m
21
m
22
__
..
1
..
2
_
+
_
k
11
0
0 k
22
__
1
2
_
= 0 .
4.2 FORCED HARMONIC OSCILLATOR
If an external weak force acts on an oscillator system the oscillations of the
system are known as forced oscillations.
Besides its proper potential energy the system gets a supplementary po-
tential energy U
e
(x, t) due to the external eld. Expanding the latter in a
Taylor series of small amplitudes x we get:
U
e
(x, t)
= U
e
(0, t) +x
_
U
e
x
_
x=0
.
The second term is the external force acting on the system at its equilibrium
position, that we denote by F(t). Then, the Lagrangian reads
L =
1
2
m
x
2
1
2
kx
2
+ xF(t) . (23)
The corresponding equation of motion is
58
m
x +kx = F(t) ,
or
x +w
2
x = F(t)/m , (24)
where w is the frequency of the proper oscillations. The general solution of
this equation is the sum of the solution of the homogeneous equation and a
particular solution of the nonhomogeneous equation
x = x
h
+ x
p
.
We shall study the case in which the external force is periodic in time of
frequency forma
F(t) = f cos(t +) .
The particular solution of (24) is sought in the form x
1
= b cos(t + ) and
by substituting it one nds that the relationship b = f/m(w
2
2
) should
be fullled. Adding up both solutions, one gets the general solution
x = a cos(wt + ) +
_
f/m(w
2
2
)
_
cos(t + ) . (25)
This result is a sum of two oscillations: one due to the proper frequency and
another at the frequency of the external force.
The equation (24) can in general be integrated for an arbitrary external
force. Writing it in the form
d
dt
(
x +iwx) iw(
x +iwx) =
1
m
F(t) ,
and making =
x +iwx, we have
d
dt
iw = F(t)/m .
The solution to this equation is = A(t) exp(iwt); for A(t) one gets
A= F(t) exp(iwt)/m .
Integrating it leads to the solution
59
= exp(iwt)
_
t
0
1
m
F(t) exp(iwt)dt +
o
. (26)
This is the general solution we look for; the function x(t) is given by the
imaginary part of the general solution divided by w.
EXAMPLE
We give here an example of employing the previous equation.
Determine the nal amplitude of oscillations of a system acted by an extenal
force F
0
= const. during a limited time T. For this time interval we have
=
F
0
m
exp(iwt)
_
T
0
exp(iwt)dt ,
=
F
0
iwm
[1 exp(iwt)] exp(iwt) .
Using [[
2
= a
2
w
2
we obtain
a =
2F
0
mw
2
sin(
1
2
wT) .
4.3 DAMPED HARMONIC OSCILLATOR
Until now we have studied oscillatory motions in free space (the vacuum),
or when the eects of the medium through which the oscillator moves are
negligeable. However, when a system moves through a medium its motion
is retarded by the reaction of the latter. There is a dissipation of the energy
in heat or other forms of energy. We are interested in a simple description
of the dissipation phenomena.
The reaction of the medium can be imagined in terms of friction forces.
When they are small we can expand them in powers of the velocity. The
zero-order term is zero because there is no friction force acting on a body
at rest. Thus, the lowest order nonzero term is proportional to the velocity,
and moreover we shall neglect all higher-order terms
f
r
=
x ,
where x is the generalized coordinate and is a positive coecient; the
minus sign shows the oposite direction to that of the moving system. Adding
this force to the equation of motion we get
m
..
x= kx
x ,
60
or
..
x= kx/m
x /m . (27)
Writing k/m = w
2
o
and /m = 2; where w
o
is the frequency of free oscilla-
tions of the system and is the damping coecient. Therefore
..
x +2
x +w
2
o
x = 0 .
The solution of this equation is sought of the type x = exp(rt), which we
substitute back in the equation to get the characteristic equation for r. Thus
r
2
+ 2 +w
2
o
= 0 ,
where from
r
1,2
=
_
(
2
w
2
o
) .
We are thus led to the following general solution of the equation of motion
x = c
1
exp(r
1
t) + c
2
exp(r
2
t) .
Among the roots r we shall look at the following particular cases:
(i) < w
o
. One gets complex conjugate solutions. The solution is
x = Re
_
Aexp
_
t +i
_
(w
2
o
2
)
__
,
where A is an arbitrary complex constant. The solution can be written of
the form
x = a exp(t) cos(wt +), where w =
_
(w
2
o
2
) , (28)
where a and are real constants. Thus, one can say that a damped oscil-
lation is a harmonic oscillation with an exponentially decreasing amplitude.
The rate of decreasing of the amplitude is determined by the exponent .
Moreover, the frequency w is smaller than that of free oscillations.
(ii) > w
o
. Then, both r are real and negative. The general form of the
solution is:
x = c
1
exp
_
_
_
(
2
w
2
o
)
_
t
_
+ c
2
exp
_
_
+
_
(
2
w
2
o
)
_
t
_
.
61
If the friction is large, the motion is just a monotone decaying amplitude
asymptotically (t ) tending to the equilibrium position (without oscil-
lations). This type of motion is called aperiodic.
(iii) = w
o
. Then r = , whose general solution is
x = (c
1
+c
2
t) exp(t) .
If we generalize to systems of n degrees of freedom, the generalized fric-
tion forces corresponding to the coordinates x
i
are linear functions of the
velocities
f
r,i
=
j
ij
.
x
i
. (29)
Using
ik
=
ki
, one can also write
f
r,i
=
F
.
x
i
,
where F =
1
2
i,j
ij
.
x
i
.
x
j
is called the dissipative function. The dierential
equation is obtained by adding up all these forces to (14)
(m
ij
x
j
+k
ij
x
j
) =
ij
.
x
i
. (30)
Employing
x
k
= A
k
exp(rt)
in (30) and deviding by exp(rt), one can obtain the following system of linear
algebraic equations for the constants A
j
j
(m
ij
r
2
+
ij
r + k
ij
)A
j
= 0 .
Making equal to zero the determinant of this system, one gets the corre-
sponding characteristic equation
m
ij
r
2
+
ij
r + k
ij
= 0 . (31)
This is an equation for r of degree 2n.
62
4.4 NORMAL MODES
Before dening the normal modes, we rewrite (15) as follows
M
..
X
_
+ K[X) = 0 ,
where [X) is the n-dimensional vector whose matrix representation is (19);
M and K are two operators having the matrix representation given by (16)
and (17), respectively. We have thus an operatorial equation. Since M is
a nonsingular and symmetric operator, the inverse operator M
1
and the
operators M
1/2
and M
1/2
are well dened. In this case, we can express the
operatorial equation in the form
d
2
dt
2
M
1/2
[X) = M
1/2
KM
1/2
M
1/2
[X) ,
or more compactly
d
2
dt
2
X
_
=
X
_
, (32)
where
X
_
= M
1/2
[X)
and
= M
1/2
KM
1/2
.
Since M
1/2
and K are symmetric operators, then is also symmetric. If
we use orthogonal eigenvectors as a vectorial base (for example, the three-
dimensional Euclidean space), the matrix representation of the operator can
be diagonal, e.g.,
ij
=
i
ij
.
Let us consider the following eigenvalue problem
[
i
) =
i
[
i
) , (33)
where [
i
) is an orthogonal set of eigenvectors. Or
M
1/2
KM
1/2
[
i
) =
i
[
i
) .
63
The eigenvalues are obtained by multiplying both sides by
i
[, leading to
i
=
i
[ M
1/2
KM
1/2
[
i
)
i
[
i
)
.
Since the potential and kinetic energies are considered positive quantities,
one should take
i
[ M
1/2
KM
1/2
[
i
))0
and therefore
i
> 0 .
This leads to the set
i
= w
2
i
.
If we express the vector
X
_
in terms of these eigenvectors of ,
X
_
=
i
y
i
X
_
,
where
y
i
=
i
X
_
. (34)
Inserting this result in the equation of motion (32), we obtain
d
2
dt
2
i
y
i
[
i
) =
X
_
=
i
y
i
[
i
) .
The scalar product of this equation with the constant eigenvector
j
[ leads
to the equation of motion for the generalized coordinate y
j
d
2
dt
2
y
j
= w
2
j
y
j
.
The solution of this equation reads
y
j
= A
j
cos(w
j
t +
j
) . (35)
64
Use of these new generalized harmonic coordinates lead to a set of indepen-
dent equations of motion. The relationship between y
j
and x
i
is given by
(34)
y
j
=
j1
x
1
+
j2
x
2
+... +
jn
x
n
.
The components
jl
(l = 1, 2, .., n) are determined by solving the eigenvalue
problem given by (33). The new coordinates are called normal coordinates
and the w
j
are known as the normal frequencies. The equivalent matrix
form (35) is
_
_
_
_
_
_
_
x
(j)
1
x
(j)
2
.
.
.
x
(j)
n
_
_
_
_
_
_
_
= A
j
cos(w
j
t +
j
)
_
_
_
_
_
_
j1
j2
.
.
.
jn
_
_
_
_
_
_
. (36)
These are the normal vibrational modes of the system. One reason for
introducing the coordinates y
j
is found from the expression for the kinetic
energy, which is seen to be invariant under the rotation to the new axes.
T =
1
2
n
j=1
M
j
.
y
2
j
.
EXAMPLE
Apply the matricial procedure as already shown, given the following equa-
tions of motion
d
2
dt
2
_
_
_
x
1
x
2
x
3
_
_
_ =
_
_
_
5 0 1
0 2 0
1 0 5
_
_
_
_
_
_
x
1
x
2
x
3
_
_
_ .
Comparing with (32), we identify the operator . To nd the eigenvectores,
we use (33) getting
_
_
_
5 0 1
0 2 0
1 0 5
_
_
_
_
_
_
3
_
_
_ =
i
_
_
_
3
_
_
_ .
The characteristic equation for
i
is
det(
i
I) = 0 ,
65
and by substituting the values
5 0 1
0 2 0
1 0 5
= 0 .
Solving the equation one gets
i
= 2, 4, 6. For = 4
_
_
_
5 0 1
0 2 0
1 0 5
_
_
_
_
_
_
3
_
_
_ = 4
_
_
_
3
_
_
_
we have the following set of equations
(5 4)
1
+
3
= 0
2
2
4
2
= 0
1
+ (5 4)
3
= 0 .
Taking into account the normalization condition, one is led to the following
values
1
=
3
=
1
2
= 0 .
Therefore
[
=4
) =
1
2
_
_
_
1
0
1
_
_
_ .
By the same means one gets
[
=6
) =
1
2
_
_
_
1
0
1
_
_
_
[
=2
) =
_
_
_
0
1
0
_
_
_ .
66
Thus, the new vectorial space is determined by
[
i
) =
_
_
_
1
2
1
2
0
0 0 1
2
1
2
0
_
_
_ ,
where from
i
[ =
_
_
_
1
2
0
1
2
1
2
0
1
2
0 1 0
_
_
_ .
Thus, the normal coordinates are given by (34)
_
_
_
y
1
y
2
y
3
_
_
_ =
_
_
_
1
2
0
1
2
1
2
0
1
2
0 1 0
_
_
_
_
_
_
x
1
x
2
x
3
_
_
_ .
4.5 PARAMETRIC RESONANCE
The important phenomenon of parametric resonance shows up for systems
initially at rest in un unstable equilibrium position, say x = 0; thus, the
slightest deviation fromthis position produces a displacement growing rapidly
(exponentially) in time. This is dierent from the ordinary resonances,
where the displacement grows only linearly in time.
The parameters of a linear system are the coecients m and k of the La-
grangian (4); if they are functions of time, the equation of motion is:
d
dt
(m
x) +kx = 0 . (37)
If we take a constant mass, the previous equation can be written in the form
d
2
x
dt
2
+ w
2
(t)x = 0 . (38)
The function w(t) is given by the problem at hand. Assuming it a periodic
function of frequency (of period T = 2/), i.e.,
w(t + T) = w(t) ,
any equation of the type (38) is invariant w.r.t. the transformation t t+T.
Thus, if x(t) is one of its solutions, x(t +T) is also a solution. Let x
1
(t) and
67
x
2
(t) be two independent solutions of 4.6.2). They should change to itselves
in a linear combination when t t +T. Thus, one gets
x
1
(t +T) =
1
x(t) (39)
x
2
(t +T) =
2
x(t) ,
or, in general
x
1
(t) =
t/T
1
F(t)
x
2
(t) =
t/T
2
G(t) ,
where F(t) and G(t) are periodical functions in time of period T. The
relationship between these constants can be obtained by manipulating the
following equations
..
x
1
+w
2
(t)x
1
= 0
..
x
2
+w
2
(t)x
2
= 0 .
Multiplying by x
2
and x
1
, respectively, and substracting term by term, we
get
..
x
1
x
2
..
x
2
x
1
=
d
dt
(
.
x
1
x
2
.
x
2
x
1
) = 0 ,
or
.
x
1
x
2
.
x
2
x
1
= const. .
Substituting t by t + T in the previous equation, the right hand side is
multiplied by
1
2
(see eqs. (39)); thus, it is obvious that the following
condition holds
2
= 1 , (40)
where one should take into account (38) and the fact that the coecients
are real. If x(t) is one integral of this equation, then x
1
,
2
. This leads to either
1
=
2
or
1
and
2
both real. In the rst case, based on (40) one gets
1
= 1/
1
,
68
that is equivalent to [
1
[
2
= [
2
[
2
= 1. In the second case, the two solutions
are of the form
x
1
(t) =
t/T
F(t)
x
2
(t) =
t/T
G(t) .
One of these functions grows exponentially in time, which is the character-
istic feature of the parametric resonance.
REFERENCES AND FURTHER READING
* H. Goldstein, Classical mechanics, Second ed. (Addison-Wesley, 1981).
* L. D. Landau & E. M. Lifshitz, Mechanics, (Pergammon, 1976).
* W. Hauser, Introduction to the principles of mechanics, (Wesley, 1965).
* E.I. Butikov, Parametric Resonance, Computing in Science &
Engineering, May/June 1999, pp. 76-83 (http://computer.org).
69
5. CANONICAL TRANSFORMATIONS
Forward: The main idea of canonical transformations is to nd all those
coordinate systems in the phase space for which the form of the Hamilton
eqs is invariant for whatever Hamiltonian. In applications one chooses the
coordinate system that allows a simple solution of the problem at hand.
CONTENTS:
5.1 Denitions, Hamiltonians and Kamiltonians
5.2 Necessary and sucient conditions for canonicity
5.3 Example of application of a canonical transformation
70
5.1 Denitions, Hamiltonians and Kamiltonians
For the time-independent and time-dependent cases, respectively, one de-
nes a canonical transformation as follows
Denition 1: A time-independent transformation Q = Q(q, p), and P =
P(q, p) is called canonical if and only if there is a function F(q, p) such that
dF(q, p) =
i
p
i
dq
i
i
P
i
(q, p)dQ
i
(q, p) .
Denition 2: A time-dependent transformation Q = Q(q, p, t), and P =
P(q, p, t) is called canonical if and only if there is a function F(q, p, t) such
that for an arbitrary xed time t = t
0
dF(p, q, t
0
) =
i
p
i
dq
i
i
P
i
(q, p, t
0
)dQ
i
(p, q, t
0
) ,
where
dF(p, q, t
0
) =
i
F(p, q, t
0
)
q
i
dq
i
+
i
F(p, q, t
0
)
p
i
dp
i
and
dQ(p, q, t
0
) =
i
Q(p, q, t
0
)
q
i
dq
i
+
i
Q(p, q, t
0
)
p
i
dp
i
Example: Prove that the following transformation is canonical
P =
1
2
(p
2
+ q
2
)
Q = Tan
1
_
q
p
_
.
Solution: According to the rst denition we have to check that pdqPdQ
is an exact dierential. Substituting P and Q in the denition we get
pdq PdQ = pdq
1
2
(p
2
+ q
2
)
pdq qdq
p
2
+ q
2
= d
_
pq
2
_
.
We can see that indeed the given transformation is canonical. We know
that a dynamical system is usually characterized by a Hamiltonian H =
H(q, p, t), where q = q(q
1
, q
2
, ..., q
n
), and p = p(p
1
, p
2
, ..., p
n
). Therefore,
the dynamics of the system fullls a set of 2n rst-order dierential eqs
71
(Hamiltons eqs.)
q
i
=
H
p
i
(1)
p
i
=
H
q
i
. (2)
Let us denote the coordinate transformations in the phase space by
Q
j
= Q
j
(q, p, t) (3)
P
j
= P
j
(q, p, t) . (4)
According to the aforementioned principle for the set of canonical trans-
formations denoted by (3) and (4), analogously to (1) and (2), there is a
function K = K(Q, P, t) such that we can write
Q
i
=
K
P
i
(5)
P
i
=
K
Q
i
. (6)
The relationship between the Hamiltonian H and the Kamiltonian K
1
can
be obtained arguing as follows
2
.
According to Hamiltons principle, the real trajectory of a classical system
can be obtained from the variation of the action integral
_
(
i
p
i
dq
i
Hdt) = 0 . (7)
If the transformation is canonical, the Kamiltonian K should fulll a rela-
tionship similar to (7). In other words, for the new set of variables Q and
P we still have
_
(
i
P
i
dQ
i
Kdt) = 0 . (8)
1
Here we follow the terminology of Goldstein by referring to K = K(Q, P, t), which
is dierent from the Hamiltonian H = H(p, q, t) by an additive time derivative, as the
Kamiltonian.
2
An alternative derivation has been given by G. S. S. Ludford and D. W. Yannitell,
Am. J. Phys. 36, 231 (1968).
72
Moreover, according to the Legendre transformation,
i
p
i
dq
i
Hdt =
L(q, q, t)dt, (7) - like (8)- is equivalent to
_
t
2
t
1
L(q, q, t)dt = 0 . (9)
In addition, (9) does not change if L is replaced by L = L+
dF(q,t)
dt
because
in this case
_
t
2
t
1
Ldt =
_
t
2
t
1
(L +
dF(q, t)
dt
)dt , (10)
or similarly
_
t
2
t
1
Ldt =
_
t
2
t
1
L(q, q, t)dt + F(q
(2)
, t
2
) F(q
(1)
, t
1
) . (11)
Thus, (10) and (11) dier only by constant terms whose variation is zero in
Hamiltons principle.
It follows that the Hamiltonian H and the Kamiltonian K are related by
the equation
3
p
i
q
i
H = P
i
Q
i
K +
dF
dt
. (12)
The function F is the so-called generating function. It can be expressed
as a function of any arbitrary set of independent variables. However, some
very convenient results are obtained if F is expressed as a function of the n
old variables and the n new ones, plus the time. The results are especially
convenient if the n old variables are exactly the n q
i
- or the n p
i
-, and if
the new variables are all of them the n Q
i
- or the n P
i
.
Using these coordinates, the possible combinations of n old variables and n
new variables -including t- in the generating function are
4
F
1
= F
1
(Q, q, t) (13)
F
2
= F
2
(P, q, t)
F
3
= F
3
(Q, p, t)
F
4
= F
4
(P, p, t) .
3
Some authors add to the right hand side of this equation a constant multiplicative
factor A that does not change (9). Here we use A = 1, that is, we decided to work
with the so-called reduced canonical transformations, since this simpler case is sucient
to illustrate the structure of the canonical transformations.
4
We shall use the convention of Goldstein to denote each of the dierent combinations
of the new and old variables in the generating function.
73
On the other hand, if we multiply (12) by dt we get:
p
i
dq
i
Hdt = PdQ
i
Kdt +dF . (14)
Making the change F F
1
above, and recalling that dQ
i
, dq
i
, and dt are
independent variables, we get:
P
i
=
F
1
Q
i
p
i
=
F
1
q
i
K = H +
F
1
t
.
Using now some algebraic manipulation it is possible to obtain analogous
expressions to the previous one for the rest of the generating functions. The
results are the following:
F
2
: Q
i
=
F
2
P
i
p
i
=
F
2
q
i
K = H +
F
2
t
F
3
: P
i
=
F
3
Q
i
q
i
=
F
3
p
i
K = H +
F
3
t
F
4
: Q
i
=
F
4
P
i
q
i
=
F
4
p
i
K = H +
F
4
t
.
In practice, one usually applies a useful theorem (see below) that allows,
together with the denitions we gave in the introduction for canonical trans-
formations, to solve any mechanical problem of interest
5
.
Theorem 5.1 We consider a system acted by a given external force. We
also suppose that the dynamical state of the system is determined by a set
of variables q, p = q
1
, q
2
, ..., q
n
, p
1
, p
2
, ..., p
n
and that the Hamiltonian of the
system is H = H(q, p, t). The time evolution of the variables q and p is
given by Hamiltons eqs.
q
i
=
H(q, p, t)
p
i
p
i
=
H(q, p, t)
q
i
.
If we now perform a transformation to the new variables
Q = Q(q, p, t) ; P = P(q, p, t)
5
For an example, see the nal section of this chapter.
74
and if the transformation is canonical, i.e., there exists a function F(q, p, t)
such that for a xed arbitrary time t = t
0
we have
dF(q, p, t
0
) =
i
y
i
dx
i
i
Y
i
dX
i
,
where x
i
, y
i
= q
i
, p
i
or p
i
, q
i
y X
i
, Y
i
= Q
i
, P
i
, or P
i
, Q
i
, then the equa-
tions of motion in terms of the variables Q and P are
Q
i
=
K(Q, P, t)
P
i
P
i
=
K(Q, P, t)
Q
i
,
where
K H +
F(q, p, t)
t
+
i
Y
i
X
i
(q, p, t)
t
.
Moreover, if the determinant of the matrix [
X
i
y
j
] is dierent of zero, then
the latter equation takes the form
K H +
F(x, X, t)
t
.
5.2 Necessary and sucient conditions for a trans-
formation to be canonical
We have already mentioned that by a canonical transformation we mean a
transformation, which, independently of the form of the Hamiltonian, keeps
unchanged the form of Hamiltons equations. However, one should be very
careful with this issue because some transformations fulll this requirement
only for a particular Hamiltonian
6
. Some authors call such transformations
canonical transformations w.r.t. H
7
.
To illustrate this point we use the following example, given in the paper
of J. Hurley: Let us consider a particular physical system whose Hamilto-
nian is
H =
p
2
2m
(15)
6
See, for example, J. Hurley Am. J. Phys. 40, 533 (1972).
7
See, for example, R. A. Matzner and L. C. Shepley, Classical Mechanics (Prentice
Hall, 1991).
75
and the following transformations
P = p
2
Q = q .
(16)
It is easy to show that the Kamiltonian K given by
K =
2P
3/2
3m
leads to
P = 2p p = 0 =
K
Q
and
Q = q =
p
m
=
P
1/2
m
=
K
P
. (17)
On the other hand, if we choose the Hamiltonian
H =
p
2
2m
+q
2
,
then it is possible to nd a Kamiltonian K for which the usage of the trans-
formation equations (16) maintains unchanged the form of Hamiltons equa-
tions. Thus, the equations (16) keeps unchanged the form of Hamiltons
equations only for a particular Hamiltonian.
It can be shown that the necessary and sucient conditions for the canonic-
ity of transformations of the form (3) and (4), that is, to keep unchanged the
form of Hamiltons equations whatever the Hamiltonian, are the following
[Q
i
, P
j
] = (18)
[P
i
, P
j
] = 0 (19)
[Q
i
, Q
j
] = 0 , (20)
where is an arbitrary constant related to scale changes. Finally, we would
like to make a few important comments before closing this section. First, we
should keep in mind that Q and P are not variables dening the conguration
of the system, i.e., they are not in general a set of generalized coordinates
8
. To distinguish Q and P from the generalized coordinates q and p, one
8
Except for the trivial case in which the canonical transformation is Q = q and P = p.
76
calls them canonical variables. In addition, the equations of motion -similar
in form to the Hamiltonian ones- for the generalized coordinates q and p-
that one gets for Q and P are called canonical Hamilton equations. Second,
although we did not check here, if the transformation Q = Q(q, p, t) and
P = P(q, p, t) is canonical, then its inverse q = q(Q, P, t) and p = p(Q, P, t)
is also canonical
9
.
5.3 Example of application of TC
As we already mentioned in the introduction, the main idea in performing
a canonical transformation is to nd a phase space coordinate system for
which the form of the Hamilton eqs is maintaind whatever the Hamiltonian
and to choose the one that makes easy the solution of the problem. We
illustrate this important fact with the following example.
EXAMPLE:
The Hamiltonian of a physical system is given by H =
2
p(q +t)
2
, where
is a constant. Determine q as a function of time.
Solution:
1. Solving the Hamilton equations for the variables q and p. Applying (1)
and (2) to the given Hamiltonian we get
2
(q + t)
2
= q, 2
2
p(q + t) = p .
This system is not easy to solve. However, we can get the solution by means
of an appropriate canonical transformation as we show in the following.
2. Using Q = q +t, P = p. According to the theorem given in section (5.1),
since
Q
p
= 0
P
(q)
= 0 ,
then the Kamiltonian K of the system is given by
K = H +
F(q, p, t)
t
+ P
Q
t
Q
P
t
. (21)
9
For a proof see, for example, E. A. Desloge, Classical Mechanics, Volume 2 (John
Wiley & Sons, 1982).
77
The form of the function F(q, p, t) can be obtained from the canonical trans-
formation given in section 5.1 (the case corresponds to a time-dependent
canonical transformation). Therefore, we substitute Q = q +t, P = p in
dF(q, p, t) = pdq PdQ ,
to get without diculty
F(q, p, t) = c, c= constant .
On the other hand,
P
t
= 0
Q
t
= 1 .
Finally, substituting these results in (21) (and also Q = q + t, P = p in H)
we get
K = P(
2
Q
2
+ 1) .
Moreover, from (5) we nd
Q =
2
Q
2
+ 1 .
This dierential equation is now easy to solve, leading to
q =
1
tan(t + ) t ,
where is an arbitrary phase.
78
6. POISSON BRACKETS
Forward: The Poisson brackets are very useful analytical tools for the
study of any dynamical system. Here, we dene them, give some of their
properties, and nally present several applications.
CONTENTS:
1. Denition and properties
2. Poisson formulation of the equations of motion
3. The constants of motion in Poisson formulation
79
1. Denition and properties of Poisson brackets
If u and v are any two quantities that depend on the dynamical state of a
system, i.e., on p and q) and possibly on time, the Poisson bracket of u and
v w.r.t. a set of canonical variables q and p
10
is dened as follows
[u, v]
i
_
u(q, p, t)
q
i
v(q, p, t)
p
i
u(q, p, t)
p
i
v(q, p, t)
q
i
_
. (1)
The Poisson brackets have the following properties (for u, v, and w arbitrary
functions of q, p, and t; a is an arbitrary constant, and r is any of q
i
, p
i
or
t)
11
:
1. [u, v] [v, u]
2. [u, u] 0
3. [u, v +w] [u, v] + [u, w]
4. [u, vw] v[u, w] + [u, v]w
5. a[u, v] [au, v] [u, av]
6.
[u,v]
r
[
u
r
, v] + [u,
v
r
]
7. The Jacobi identity, [u, [v, w]] + [v, [w, u]] + [w, [u, v]] 0 .
Another very important property of PBs is the content of the following
theorem
Theorem 6.1 If the transformation Q = Q(q, p, t), P = P(q, p, t) is a
canonical transformation, the PB of u and v w.r.t. the variables q, p is
equal to the PB of u and v w.r.t. the set of variables Q, P, i.e.,
i
_
u(q, p, t)
q
i
v(q, p, t)
p
i
u(q, p, t)
p
i
v(q, p, t)
q
i
_
=
i
_
u(q, p, t)
Q
i
v(q, p, t)
P
i
u(q, p, t)
P
i
v(q, p, t)
Q
i
_
.
10
As in the previous chapter, by q and p we mean q = q1, q2, ..., qn y p = p1, p2, ..., pn.
11
The proof of these properties can be obtained by using the denition of the PBs in
order to express each term of these identities by means of partial derivatives of u, v, and
w, and noticing by inspection that the resulting equations do hold.
80
2. Poisson formulation of the equations of motion
In the following, we outline as theorem-like statements the most important
results on the PB formulation of the eqs of motion of the dynamical systems
12
:
Theorem 6.2 Consider a system whose dynamical state is dened by the
canonical variables q, p and whose dynamical behaviour is dened by the
Hamiltonian H = H(q, p, t). Let F be an arbitrary quantity depending on
the dynamical state of the system, i.e., on q, p, and possibly on t. The rate
of change in time of F is given by
F = [F, H] +
F(q, p, t)
t
,
where [F, H] is the PB of F and H.
Theorem 6.3 (Poisson formulation of the eqs. of motion). Con-
sider a system described in terms of the canonical variables q, p, and whose
Hamiltonian is H = H(q, p, t). The motion of the system is governed in this
case by the equations
q
i
= [q
i
, H]
p
i
= [p
I
, H] .
3. Constants of motion in Poissons formulation
We shall use again a theorem-like sketch of the basic results on the constants
of motion in Poissons formulation. These results are the following.
Theorem 6.4 If one dynamical quantity F is not an explicit function of time
and if the PB of F and H is zero, [F, H] = 0, then F is a constant/integral
of motion as one can see from the theorem 6.2.
Corollary 6.4. If the Hamiltonian is not an explicit function of time, then
it is a constant of motion.
12
The proofs have been omitted as being well known. See, for example, E. A. Desloge,
Classical Mechanics, Volume 2 (John Wiley & Sons, 1982).
81
7. HAMILTON-JACOBI EQUATIONS
Forward: It is known from the previous chapters that in principle it is
possible to reduce the complexity of many dynamical problems by choosing
an appropriate canonical transformation. In particular, we can try to look
for those canonical transformations for which the Kamiltonian K is zero, a
situation leading to the Hamilton-Jacobi equations.
CONTENTS:
7.1 Introduction
7.2 Time-dependent Hamilton-Jacobi equations
7.3 Time-independent Hamilton-Jacobi equations
7.4 Generalization of the Hamilton-Jacobi equations
7.5 Example of application of the Hamilton-Jacobi equations
82
7.1 Introduction
In order to reach the goals of this chapter we need to make use of the
following result allowing us to nd the set of canonical variables for which
the Kamiltonian takes a particular form.
Theorem 7.1 Consider a system whose dynamical state is dened by p, q
and whose behaviour under the action of a given force is governed by the
Hamiltonian H = H(q, p, t). Let K = K(Q, P, t) be a known function of
the canonical variables Q, P, and time. Then, any function F(q, Q, t) that
satises the partial dierential equation
K
_
Q,
F(q, Q, t)
Q
, t
_
= H
_
q,
F(q, Q, t)
q
, t
_
+
F(q, Q, t)
t
and also the condition
2
F(q, Q, t)
q
j
Q
j
,= 0
is a generating function for a canonical transformation of q, p to Q, P, and
the corresponding Kamiltonian is K = K(Q, P, t).
In the following sections we shall use this theorem to nd those canonical
transformations whose Kamiltonian is zero
13
, that leads us to the Hamilton-
Jacobi equations.
7.2 Time-dependent HJ equations.
As a consequence of Theorem 7.1 and of requiring a zero Kamiltonian we
get the following theorem
Theorem 7.2 Consider a system of f degrees of freedom dened by the set
of variables q, p and of Hamiltonian H = H(q, p, t). If we build the partial
dierential equation
H
_
q,
S(q, t)
q
, t
_
+
S(q, t)
t
= 0 (1)
and if we are able to nd a solution of the form
S = S(q, , t) ,
13
More exactly, K
_
Q,
F(q,Q,t)
Q
, t
= 0.
83
where =
1
,
2
, ...,
f
is a set of constants and if in addition the solution
satises the condition
2
S(q, , t)
q
i
,= 0 ,
then q(t) can be obtained from the equations
S(q, , t)
i
=
i
, (2)
where =
1
,
2
, ...,
f
is a set of constants. The set of equations (2) provide
us with f algebraic equations in the f unknown variables q
1
, q
2
, ..., q
f
. The
values of the constants and are determined by the boundary conditions.
Moreover, if q(t) is given it is possible to nd p(t) starting from
p
i
=
S(q, , t)
q
i
. (3)
The partial dierential equation (1) is called the time-dependent Hamilton-
Jacobi equation. The function S(q, , t) is known as Hamiltons principal
function.
To achieve a better meaning of the theorem, as well as of the constants
and , we proceed with its proof.
Proof of the Theorem 7.2. According to the Theorem 7.1, any function
F(q, Q, t) satisfying the partial dierential equation
H
_
q,
F(q, Q, t)
q
, t
_
+
F(q, Q, t)
t
= 0
and also the condition
2
F(q, Q, t)
q
j
Q
j
,= 0
should be a generating function of a set of canonical variables Q, P for which
the Kamiltonian K is zero, i.e., K(Q, P, t) = 0. The function
F(q, Q, t) = [S(q, , t)]
=Q
S(q, Q, t)
belongs to this class. Then S(q, Q, t) is the generating function for a canon-
ical transformation leading to the new set of canonical variables Q, P, for
84
which the Kamiltonian K is identically zero. The transformation equations
associated to S(q, Q, t) are
p
i
=
S(q, Q, t)
q
i
(4)
P
i
=
S(q, Q, t)
Q
i
(5)
and because K(Q, P, t) 0, the equations of motion are
Q
i
=
K(Q,P,t)
P
i
= 0
P
i
=
K(Q,P,t)
Q
i
= 0
.
From these equations we infer that
Q
i
=
i
(6)
P
i
= , (7)
where
i
and
i
are constants. The choice of the negative sign for in (7)
is only a convention. If now we substitute the equations (6) and (7) in (5)
we get
i
=
_
S(q, Q, t)
Q
i
_
Q=
=
S(q, , t)
i
,
which reduces to (2). If, in addition, we substitute (6) in (4) we get (2).
This complets the proof.
7.3 Time-independent HJ equations
If the Hamiltonian does not depend explicitly on time, we can partially solve
the time-dependent Hamilton-Jacobi equation. This result can be spelled
out as the following theorem
Theorem 7.3 Consider a system of f degrees of freedom dened in terms
of q, p, and whose behaviour under a given force is governed by the time-
independent Hamiltonian H(q, p).
If we build the partial dierential equation
H
_
q,
W(q)
q
_
= E , (8)
85
where E is a constant whose value for a particular set of conditions is equal to
the value of the integral of motion H(q, p) for the given boundary conditions,
and if we can nd a solution to this equation of the form
W = W(q, ) ,
where
1
,
2
, ...,
f
is a set of constants that explicitly or implicitly
include the constant E, i.e., E = E(), and if the solution satises the
condition
2
W(q, )
q
i
,= 0 ,
then the equations of motion are given by
S(q, , t)
i
=
i
(9)
where
S(q, , t) W(q, ) E()t
and =
1
,
2
, ...,
f
is a set of constants. The set of equations (9) pro-
vide f algebraic equations in the f unknown variables q
1
, q
2
, ..., q
f
. The
values of the constants and are determined by the boundary conditions.
The partial dierential equation (8) is the time-independent Hamilton-Jacobi
equation, and the function W(q, ) is known as the characteristic Hamilton
function.
7.4 Generalization of the HJ equations
The Hamilton-Jacobi equation can be generalized according to the follow-
ing theorem allowing sometimes the simplication of some Hamilton-Jacobi
problems.
Theorem 7.4 Consider a system of f degrees of freedom whose dynamics
is dened by x, y, where x
i
, y
i
= q
i
, p
i
or p
i
, q
i
, and whose behaviour under
the action of a given force is governed by the Hamiltonian H(x, y, t). If we
write the partial dierential equation
H
_
x,
S(x, t)
x
, t
_
+
S(x, t)
t
= 0
and if we can nd a solution of this equation of the form
S = S(x, , t)
86
where
1
,
2
, ...,
f
is a set of constants and in addition the solution
satises the condition
2
S(x, , t)
x
j
,= 0 ,
then the laws of motion of the system can be obtained from the equations
S(x, , t)
x
i
= y
i
(10)
S(x, , t)
i
=
i
, (11)
where
1
,
2
, ...,
f
is a set of constants.
7.5 Example of application of the HJ equations
We shall solve the problem of the one-dimensional harmonic oscillator of
mass m, using the Hamilton-Jacobi method.
We know that the Hamiltonian of the system is
H =
p
2
2m
+
kx
2
2
(12)
According to Theorem 7.2 the Hamilton-Jacobi equation for the system is
1
2m
_
F
q
_
+
kq
2
2
+
F
t
= 0 (13)
We assume a solution of (13) of the form F = F
1
(q) +F
2
(t). Therefore, (13)
converts to
1
2m
_
dF
1
dq
_
2
+
kq
2
2
=
dF
2
dt
(14)
Making each side of the previous equation equal to , we nd
1
2m
_
dF
1
dq
_
2
+
kq
2
2
= (15)
dF
2
dt
= (16)
For zero constants of integration, the solutions are
F
1
=
_
2m(
kq
2
2
)dq (17)
F
2
= t (18)
87
Thus, the generating function F is
F =
_
2m(
kq
2
2
)dq t . (19)
According to (2), q(t) is obtained starting from
=
_
_
_
_
2m(
kq
2
2
)dq t
_
_
_
(20)
=
2m
2
_
dq
_
kq
2
2
t (21)
and eecting the integral we get
_
m
k
sin
1
(q
_
k/2) = t + , (22)
from which q is nally obtained in the form
q =
_
2
k
sin
_
k/m(t +) . (23)
In addition, we can give a physical interpretation to the constant according
to the following argument.
The factor
_
2
k
should correspond to the amplitude A of the oscillator. On
the other hand, the total energy E of a one-dimensional harmonic oscillator
of amplitude A is given by
E =
1
2
kA
2
=
1
2
k
_
_
2
k
_
2
= .
In other words, is physically the total energy E of the one-dimensional
harmonic oscillator.
FURTHER READING
C.C. Yan, Simplied derivation of the HJ eq., Am. J. Phys. 52, 555 (1984)
88
N. Anderson & A.M. Arthurs, Note on a HJ approach to the rocket pb., Eur.
J. Phys. 18, 404 (1997)
M.A. Peterson, Analogy between thermodynamics and mechanics, Am. J.
Phys. 47, 488 (1979)
Y. Hosotani & R. Nakayama, The HJ eqs for strings and p-branes, hep-
th/9903193 (1999)
89
8. ACTION-ANGLE VARIABLES
Forward: The Hamilton-Jacobi equation provides a link to going from a
set of canonical variables q, p to a second set Q, P, where both are constants
of motion.
In this chapter, we briey present another important procedure by which
one goes from an initial pair of canonical variables to a nal one, where not
both variables are simultaneously constants of motion.
CONTENS:
8.1 Separable systems
8.2 Cyclic systems
8.3 Action-angle variables
8.4 Motion in action-angle variables
8.5 Importance of action-angle variables
8.6 Example: the harmonic oscillator
90
8.1 Separable systems
Separable systems are those ones for which the Hamiltonian is not an ex-
plicit function of time, i.e.,
H = H(q, p) ,
allowing, in addition, to nd a solution of the time-independent Hamilton-
Jacobi of the form
W(q, ) =
i
W
i
(q
i
, ) .
8.2 Cyclic systems
We know that the dynamical state of a system is characterized by a set of
generalized coordinates q q
1
, q
2
, ..., q
f
and momenta p p
1
, p
2
, ..., p
f
. A
system in motion will describe an orbit in the phase space q, p At the same
time, there is an orbit in each of the subspaces q
i
, p
i
. In every plane q
i
, p
i
,
the orbit can be represented by an equation of the form p
i
= p
i
(q
i
) or a pair
of equations p
i
= p
i
(t), q
i
= q
i
(t). If for each value of i, the orbit p
i
= p
i
(q
i
)
is a closed curve in the plane q
i
p
i
, then we say that the system is cyclic.
In the gure 8 we show the two possibilities for a system to be cyclic. In
8.1a, the system is cyclic because q
i
oscillates between the limits dened by
q
i
= a and q
i
= b, whereas in gure 8.1b, the system is cyclic because q
i
moves from q
i
= a to q
i
= b, and repeats the same motion afterwards.
b a
p
i
q
i
a
b
p
i
q
i
Fig.8.1
At this point, it is worthwhile to make two helpful remarks.
Remark 1: The cyclic term has been introduced only for simplifying
the notation in the next sections. One should not interpret this term as if
91
the system is cyclic in each subspace q
i
, pi. The system should come back
to its initial state only in the global space q, p.
Remark 2: If the cyclic system has only one degree of freedom, the time
required by the system to accomplish the cycle in q p is constant; therefore
the motion in the space q p will be periodic in time. If the system has
more degrees of freedom, then, in general, the time required for a particular
cycle in one of the subspaces q
i
, p
i
will not be a constant, but will depend on
the motion of the other coordinates. As a result, the motion in the subspace
q
i
, p
i
will not be periodic in time. One should be careful with this point,
since not all the motions in the subspaces q
i
, p
i
are periodic.
8.3 Action-angle variables
We consider now a cyclic system of f degrees of freedom, whose dynamical
state is characterized by the canonical set q, p. Let H(q, p) be the Hamilto-
nian of the system and let
W(q, )
i
W
i
(q
i
, ) ,
(where =
1
,
2
, ...,
f
are constants) be a solution of the time-independent
Hamilton-Jacobi equation
H(q,
W
q
) = E .
Let J J
1
, J
2
, ..., J
f
be the set of constants dened by the equations
J
i
() =
_
W
i
(q
i
, )
q
i
dq
i
, (1)
where the integral is along a complete cycle in the variable q
i
. If we use the
function
W(q, ) W[q, (J)]
i
W
i
[q
i
, (J)]
i
W
i
(q
i
, J)
as the generating function of a canonical transformation of q, p to a new set
of coordinates w and momenta J, i.e., if we dene the variables w and J by
92
the transformation equations
p
i
=
W(q, )
q
i
=
W
i
(q
i
, J)
q
i
(2)
w
i
=
W(q, J)
J
i
, (3)
then the new coordinates w
1
, w
2
, ..., w
f
are called angle variables, and the
new momenta J
1
, J
2
, ..., J
f
are called action variables.
From (2) we get
p
i
(q
i
, ) =
W
i
[q
i
, J()]
q
i
=
W
i
(q
i
, )
q
i
. (4)
Substituting (4) in (1) one gets
J
i
() =
_
p
i
(q
i
, )dq
i
. (5)
The equation p
i
= p
i
(q
i
, ) gives the projected orbit p = p(q) on the subspace
p
i
, q
i
. The integral in the right hand side of the equation (5) is thus the area
bordered by the closed orbit, or beneath the orbit, as shown in gure 8.1.
Thus, the function J
i
() has a geometric interpretation as the area covered
in the subspace q
i
, p
i
during a complete cycle in the subspace. This area
depends on the constants or equivalently on the initial conditions and can
be arbitrary
14
.
8.4 Motion in terms of action-angle variables
We give the following theorem-like statement for the motion of a system in
terms of action-angle variables.
Theorem 8.4
Consider a separable cyclic system of f degrees of freedom whose motion
is described by the variables q, p q
1
, q
2
, ..., q
f
, p
1
, p
2
, ..., p
f
, together with
the Hamiltonian H(q, p). If we transform the motion to the action-angle
variables J, w, then the Hamiltonian H is a function of J alone, i.e.,
H = H(J)
14
Historically, the rst intents to pass from the classical mechanics to quantum mechan-
ics was related to the assumption that the value of Ji could be only a multiple of h/2,
where h is Plancks constant.
93
and the equations of motion will be
J
i
=
i
w
i
=
i
t +
i
,
where
i
y
i
are constants determined by the initial conditions, whereas
the
i
are constants known as the frequencies of the system being dened as
follows
i
=
_
H(J)
J
i
_
J=
i
.
8.5 Importance of the action-angle variables
The importance of the action-angle variables resides in providing a powerful
technique for directly getting the frequencies of periodic motions without
solving for the equations of motion of the system.
This important conclusion can be derived through the following argument.
Consider the change of w when q describes a complete cycle
w =
_
w
q
dq
On the other hand, we know that
w =
W
J
,
and therefore
w =
_
2
W
qJ
dq
=
d
dJ
_
W
q
dq
=
d
dJ
_
pdq
= 1 .
This result shows that w changes by a unity when q varies within a complete
period.
From the relationship
w = t + ,
94
we infer that in a period
w = 1
= .
This means that we can identify the constant with the inverse of the period
=
1
.
8.6 Example: The simple harmonic oscillator
Using the action-angle formalism prove that the frequency of the simple
one-dimensional harmonic oscillator is given by =
_
k/m/2.
Since H is a constant of motion, the orbit in the space q p is given
p
2
2m
+
kq
2
2
= E ,
where E is the energy. This is the equation of an ellipse of semiaxes
2mE
and
_
2E/k. The area enclosed by the ellipse is equal to the action J.
Therefore,
J =
2mE
2E
k
= 2
_
m
k
E .
It follows that
H(J) = E =
_
k/m
2
J .
Thus, the frequency will be
=
H(J)
J
=
_
k/m
2
.
95
9. CANONICAL PERTURBATION THEORY
Forward: The great majority of problems that we have to solve in Physics
are not exactly solvable. Because of this and taking into account that we
live in the epoch of computers the last decades have seen a lot of progress in
developing techniques leading to approximate solutions. The perturbation
method is used for not exactly solvable Hamiltonian problems when the
Hamiltonian diers slightly from an exactly solvable one. The dierence
between the two Hamiltonians is known as the perturbation Hamiltonian.
All perturbation methods are based on the smallness of the latter with
respect to both Hamiltonians.
CONTENTS:
9.1 Time-dependent perturbation theory (with two examples)
9.2 Time-independent perturbation theory (with an example)
96
9.1 Time-dependent perturbation theory
The most appropriate formulation of classical mechanics for the develop-
ment of perturbation methods is the Hamilton-Jacobi approach. Thus, we
take H
0
(p, q, t) as the Hamiltonian corresponding to the solvable (unper-
turbed) problem and consider the Hamilton principal function S(q,
0
, t), as
the generating function of a canonical transformation of (p, q) to the new
canonical pair (
0
,
0
) for which the new Hamiltonian (or Kamiltonian) K
0
of the unperturbed system is zero
S
t
+ H
0
(
S
q
, q, t) = K
0
= 0 . (1)
This is the Hamilton-Jacobi equation, where we used p = S/q. The whole
set of new canonical coordinates (
0
,
0
) are constant in the unperturbed
case because K
0
= 0 and:
0
=
K
0
0
,
0
=
K
0
0
. (2)
We now consider the Hamiltonian of the perturbed system written as follows:
H(q, p, t) = H
0
(q, p, t) + H(q, p, t); (H H
0
). (3)
Although (
0
,
0
) are still canonical coordinates (since the transformation
generated by S is independent of the particular form of the Hamiltonian),
they will not be constant and the Kamiltonian K of the perturbed system
will not be zero. In order not to forget that in the perturbed system the
transformed coordinates are not constant, we denote them by and in-
stead of
0
and
0
, which are the corresponding constants in the unperturbed
system. Thus, for the perturbed system we have:
K(, , t) = H +
S
t
= (H
0
+
S
t
) + H = H(, , t). (4)
The equations of motion for the transformed variables in the perturbed
system will be:
97
i
=
H(, , t)
i
=
H(, , t)
i
, (5)
where i = 1, 2, ..., n and n is the number of degrees of freedom of the system.
These are rigorous equations. If the system of 2n equations could be solved
for
i
and
i
as functions of time, the transformation equation (p, q) (, )
would give p
i
and q
i
as functions of time and the problem would be solved.
However, the exact solution of the equations (5) is not easier to get w.r.t.
the original equations. From (5) we see that when and are not constant
their variation in time is slow if we assume that H changes innitesimally
w.r.t. y . A rst approximation for the temporal variations of (, )
can be obtained by substituting y in the second terms of (5) by their
constant unperturbed values
i1
=
H(, , t)
i1
=
H(, , t)
0
, (6)
where
i1
and
i1
are the rst-order solutions, i.e., in the rst power of the
perturbation to
i
and
i
, and the vertical bars with zero subindices show
that after the derivation one should substitute and by their constant
unperturbed values. Once this is done, the equations (6) can be integrated
leading to
i
and
i
as functions of time (in the rst order). Next, using
the equations of transformation one can get p and q as functions of time in
the same rst-order approximation. The second-order approximation can
be obtained now by substituting in the second terms of (6) the rst-order
approximation of and w.r.t. time. In general, the perturbation solution
of order N is obtained by integrating the following equations
iN
=
H(, , t)
N1
iN
=
H(, , t)
N1
. (7)
Example 1
Let us consider the simple case of a free particle that next is subject to
a harmonic perturbation. Although this example is trivial can be used to
illustrate the aforementioned procedure. The unperturbed Hamiltonian is
98
H
0
=
p
2
2m
. (8)
Since H
0
,= H
0
(x), x is a cyclic variable and p =
0
is a constant of motion
in the unperturbed system. Recalling that p = S/x, we substitute in (1):
1
2m
_
S
x
_
2
+
S
t
= 0. (9)
Since the system is conservative, it is convenient to consider the principal
function of the form
S = S(x) +F(t). (10)
This type of separation of variables is quite useful when the Hamiltonian
does not depend explicitly on time. Then, one writes F(t) = Et, where E
is the total energy of the system
15
. Putting (10) in (9) we obtain
1
2m
_
dS
dx
_
2
= E, S =
2mEx =
0
x. (11)
Substituting (11) in (10), together with the fact that in this case the Hamil-
tonian is equal to the energy, we can write the principal function of Hamilton
as follows
S =
0
x
2
0
t
2m
. (12)
If the transformed momentum is
0
, the transformed coordinate (which is
also constant in the unperturbed system) is
0
=
S
0
= x
0
t
m
.
Therefore, the transformation generated by S is given by the following
eqs
15
See, M.R. Spiegel, Theoretical Mechanics, pp. 315, 316.
99
p =
0
,
x =
0
t
m
+
0
. (13)
They represent the solution for the motion of free particle. What we have
given up to now is only the procedure to obtain the equations of motion using
the Hamilton-Jacobi formulation. Now, we introduce the perturbation
H =
kx
2
2
=
m
2
x
2
2
, (14)
or, in terms of the transformed coordinates, using (13)
H =
m
2
2
_
t
m
+
_
2
. (15)
Notice that we have renounced at the subindices 0 for the transformed co-
ordinates since we already study the perturbed system.
Substituting (15) in (5) we get
= m
2
_
t
m
+
_
,
=
2
t
_
t
m
+
_
. (16)
As one may expect, these equations have an exact solution of the harmonic
type. To be sure of that we perform the time derivative of the rst equa-
tion allowing us to conclude that has a simple harmonic variation. The
same holds for x as a consequence of the transformations (13), which are
invariant in form in the perturbed system (up to the subindices of the trans-
formed coordinates). However, we are interested to illustrate the pertur-
bation method, so that we consider that k (the elastic constant) is a small
parameter. We seek approximate solutions in dierent perturbative orders,
without missing the fact that the transformed variables (, ) in the per-
turbed system are not constants of motion. In other words, even though
(, ) contain information on the unperturbed system, the eect of the per-
turbation is to make these parameters varying in time.
The rst-order perturbation is obtained in general as given by (6). Thus, we
100
have to substitute and by their unperturbed values in the second terms
of (16). To simplify, we take x(t = 0) = 0 and therefore
0
= 0, leading to
1
=
2
0
t,
1
=
0
2
t
2
m
, (17)
which integrated reads
1
=
0
0
t
2
2
,
1
=
0
2
t
3
3m
. (18)
The rst-order solutions for x and p are obtained by putting
1
and
1
in
the transformation eqs. (13), where from
x =
0
m
_
t
3
t
3
6
_
,
p =
0
_
1
2
t
2
2
_
. (19)
To obtain the approximate solution in the second perturbative order we have
to nd
2
and
2
, as was indicated in (7), by substituting in the second terms
of (16), and by
1
and
1
as given in (18). Integrating
2
and
2
and
using again the transformation eqs. (13), we get the second-order solutions
for x and p:
x =
0
m
_
t
3
t
3
3!
+
5
t
5
5!
_
,
p =
0
_
1
2
t
2
2!
+
4
t
4
4!
_
. (20)
In the limit in which the perturbation order N tends to innity, we obtain
the expected solutions compatible with the initial conditions
x
0
m
sint, p
0
cos t. (21)
The transformed variables (, ) contain information about the unperturbed
101
orbit parameters. For example, if we consider as unperturbed system that
corresponding to the Kepler problem, a convenient coordinate pair (, )
could be the (J, ) variables, which are the action and the phase angle of the
angle w, respectively (remember that w = t +, where is the frequency).
These variables are related to the set of orbital parameters such as semimajor
axis, eccentricity, inclination, and so on.
The eect of the perturbation is to produce a time variation of all these
parameters. If the perturbation is small, the variation of the parameters
during a period of the unperturbed motion will also be small. In such a
case, for small time intervals, the system moves along the so-called osculat-
ing orbit, having the same functional form as the orbit of the unperturbed
system; the dierence is that the parameters of the osculating curve vary in
time.
The osculating parameters can vary in two ways
Periodic variation: if the parameter comes back to its initial value
after a time interval that in rst approximation is usually the unper-
turbed period. These periodic eects of the perturbation do not alter
the mean values of the parameters. Therefore, the trajectory is quite
similar to the unperturbed one. These eects can be eliminated by
taking the average of the perturbations in a period of the unperturbed
motion.
Secular variation: At the end of every succesive orbital period there
is a net increment of the value of the parameter. Therefore, at the
end of many periods, the orbital parameters can be very dierent of
their unperturbed values. The instantaneous value of the variation of
a parameter, for example the frequency, is seldom of interest, because
its variation is very small in almost all cases in which the perturbation
formalism works. (This variation is so small that it is practically im-
possible to detect it in a single orbital period. This is why the secular
variation is measured after at least several periods.)
Example 2
From the theory of the Kepler two-body problem it is known that if we add
a potential 1/r
2
, the orbit of the motion of negative energy is a rotating
ellipse whose periapsis is precessing. In this example, we nd the precession
velocity for a more general perturbation
V =
k
r
h
r
n
, (22)
102
where n 2 is an entire number, and h is such that the second term of
the potential is a small perturbation of the rst one. The perturbative
Hamiltonian is
H =
h
r
n
. (23)
In the unperturbed problem, the angular position of the periapsis in the
orbit plane is given by the constant = 2w
2
. In the perturbed case, we
have
= 2
H
J
2
=
H
l
, (24)
where we have used J
2
= 2l. Moreover, J
2
and w
2
are two of the ve
integrals of motion that can be obtained when one uses the action-angle
variables to solve the Kepler problem.
We need to know the mean of in a period of the unperturbed orbit
)
1
_
0
H
l
dt =
l
_
1
_
0
H dt
_
=
H)
l
. (25)
The temporal mean of the unperturbed Hamiltonian is
H) = h
1
r
n
) =
h
_
0
dt
r
n
. (26)
On the other hand, since l = mr
2
(d/dt), we can get dt and plunge it in
(27). This leads to
H) =
mh
l
_
2
0
d
r
n2
=
mh
l
_
mk
l
2
_
n2
_
2
0
[1 +e cos ( )]
n2
d . (27)
where is a constant phase, e is the eccentricity, and where we expressed r
as a function of using the general equation of the orbit with the origin in
one of the focal points of the corresponding conic
103
1
r
=
mk
l
2
_
1 + e cos ( )
_
(28)
.
For n = 2:
H) =
2mh
l
,
) =
2mh
l
2
. (29)
For n = 3:
H) =
2m
2
hk
l
3
,
) =
6m
2
hk
l
4
. (30)
The latter case, n = 3, is of particular importance because the theory of
General Relativity predicts a correction of the Newtonian motion precisely of
r
3
order. This prediction is related to the famous problem of the precession
of the orbit of Mercury. Substituting the appropriate values of the period,
mass, great semiaxis, that is included in h, and so on, (30) predicts a mean
precession velocity of
) = 42.98 arcsec./century .
The mean value is by far larger than the aforementioned one (by a factor
larger than one hundred). But before making any comparison one should
eliminate from the mean value the contributions due to the following factors:
a) the eect known as the precession of the equinoxes (the motion of the ref-
erence point of longitudes w.r.t. the Milky Galaxy), b) the perturbations of
the Mercury orbit due to the interaction with the other planets. Once these
are eliminated, of which the rst is the most signicant, one may hopefully
obtain the contribution of the relativistic eect. In 1973, this contribution
has been estimated as 41.40.9 arcsec./century. This is consistent with the
prediction given by (30).
9.2 Time-independent perturbation theory
While in the time-dependent perturbation theory one seeks the time depen-
dence of the parameters of the unperturbed system initially considered as
104
constant, the aim of the time-independent approach is to nd the constant
quantities of the perturbed system. This theory can be applied only to con-
servative periodic systems (both in the perturbed and unperturbed state).
For example, it can be applied to planetary motion when one introduces any
type of conservative perturbation to the Kepler problem, in which case it is
known as the von Zeipel or Poincare method.
Here, we consider here the case of a periodic system of one degree of freedom
and time-independent Hamiltonian of the form
H = H(p, q, ), (31)
where is a small constant specifying the strength of the perturbation. We
assume that
H
0
(p, q) = H(p, q, 0) (32)
corresponds to a system that has an exact (closed-form) unperturbed solu-
tion in the action-angle variables (J
0
, w
0
), i.e.,
H
0
(p, q) = K
0
(J
0
)
0
= w
0
=
K
0
J
0
; (w
0
=
0
t +
0
). (33)
Since the canonical transformation from (p, q) to (J
0
, w
0
) is independent of
the particular formof the Hamiltonian, the perturbed Hamiltonian H(p, q, )
can be written as H(J
0
, w
0
, ). Due to the fact that the perturbed Hamil-
tonian depends on w
0
, J
0
, it is not constant any more. On the other hand,
in principle, one can get new action-angle variables (J, w) that may be more
appropriate for the perturbed system, such as
H(p, q, ) = E(J, )
= w =
E
J
(34)
J =
E
w
= 0 ; (J = constant).
Since the transformation connecting (p, q) to (J
0
, w
0
) is known, we have
to nd the canonical transformation S connecting (J
0
, w
0
) to (J, w). If we
assume that is small the transformation we look for should not dier much
from the identity transformation. Thus, we write the following expansion
105
S = S(w
0
, J, ) = S
0
(w
0
, J) + S
1
(w
0
, J) +
2
S
2
(w
0
, J) + ... . (35)
For = 0 we ask S to provide an identity transformation that leads to
S
0
= w
0
J . (36)
The canonical transformations generated by S read
w =
S
J
= w
0
+
S
1
J
(w
0
, J) +
2
S
2
J
(w
0
, J) +... .
J
0
=
S
w
0
= J +
S
1
w
0
(w
0
, J) +
2
S
2
w
0
(w
0
, J) + ... . (37)
Due to the fact that w
0
is an angle variable of the unperturbed system we
know that w
0
= 1 over a cycle. On the other hand, we know that the
canonical transformations have the property to conserve the phase space
volume. Therefore, we can write:
J =
_
pdq =
_
J
0
dw
0
. (38)
Integrating the second equation of (37) along an orbit of the perturbed
system, we get
_
J
0
dw
0
=
_
Jdw
0
+
n=1
n
_
S
n
w
0
dw
0
, (39)
and substituting (39) in (38) leads to
J = Jw
0
+
n=1
n
_
S
n
w
0
dw
0
. (40)
Since w
0
= 1, one gets
n=1
n
_
S
n
w
0
dw
0
= 0, (41)
or
106
_
S
n
w
0
dw
0
= 0. (42)
Moreover, the Hamiltonian can be expanded in as a function of w
0
and
J
0
:
H(w
0
, J
0
, ) = K
0
(J
0
) +K
1
(w
0
, J
0
) +
2
K
2
(w
0
, J
0
) +... , (43)
where the K
i
are known because H is a known function of w
0
and J
0
for a
given . On the other hand, one can write
H(p, q, ) = H(w
0
, J
0
, )
= E(J, ) , (44)
which is the expression for the energy in the new action-angle coordinates
(where J will be constant and w will be a linear function of time).
E can also be expanded in powers of :
E(J, ) = E
0
(J) + E
1
(J) +
2
E
2
(J) + ... . (45)
Taking into account (44) we can obtain equalities for the coecients of the
same powers of in (43) and (45). However, these expressions for the energy
depend on two dierent sets of variables. To solve this issue we express H
0
in terms of J by writing a Taylor expansion of H(w
0
, J
0
, ) w.r.t. J
0
in the
innitesimal neighborhood of J:
H(w
0
, J
0
, ) = H(w
0
, J, ) + (J
0
J)
H
J
+
(J
0
J)
2
2
2
H
J
2
+ ... , (46)
The derivatives of this Taylor expansion, which in fact are the derivatives
w.r.t. J
0
calculated for J
0
= J, can also be written as derivatives w.r.t. J,
once we substitute J
0
by J in H
0
(J
0
). In the previous equation, all the terms
containing J
0
should be rewritten in terms of J by using the transformation
dened by (37) connecting the coordinates (J
0
, w
0
) and (J, w). Thus, from
the second equation in (37) we obtain (J
0
J), which introduced in (46)
gives:
107
H(w
0
, J
0
, ) = H(w
0
, J, ) +
H
J
_
S
1
w
0
+
2
S
2
w
0
+...
_
+
1
2
2
H
J
2
2
_
S
1
w
0
_
2
+O(
3
) .(47)
Next, we can make use of (43) to write H(w
0
, J, ) = H(w
0
, J
0
, ) [
J
0
=J
,
that we put in (47) in order to get
H(w
0
, J
0
, ) = K
0
(J) + K
1
(w
0
, J) +
2
K
2
(w
0
, J) + ...
+
S
1
w
0
_
K
0
(J)
J
+
K
1
(w
0
, J)
J
+...
_
+
2
_
K
0
(J)
J
S
2
w
0
+
1
2
2
K
0
J
2
_
S
1
w
0
_
2
+...
_
E(J, ) (48)
= E
0
(J) + E
1
(J) +
2
E
2
(J) + ... .
Now, we can solve the problem in terms of the coecients E
i
(J) allowing us
to calculate the frequency of the perturbed motion in various perturbation
orders. Since the expansion of the terms of E
i
does not imply a dependence
on w
0
, then the occurrence of w
0
in (48) is articial. The K
i
(w
0
, J) of
(48) are known functions, whereas the S
i
(w
0
, J) and E
i
(J) are the unknown
quantities.
Making equal the corresponding powers of , we get
E
0
(J) = K
0
(J)
E
1
(J) = K
1
(w
0
, J) +
S
1
w
0
K
0
(J)
J
E
2
(J) = K
2
(w
0
, J) +
S
1
w
0
K
1
(w
0
, J)
J
+
1
2
_
S
1
w
0
_
2
2
K
0
(J)
J
2
+
S
2
w
0
K
0
(J)
J
. (49)
We can see that to obtain E
1
we need to know not only K
1
but also S
1
.
Moreover, we should keep in mind that the E
i
are constant, being functions
of J only, which is a constant of motion. We also have to notice that K
0
/J
does not depend on w
0
(since K
0
= K
0
(J) = K
0
(J
0
) [
J
0
=J
). Averaging over
w
0
on both sides of the second equation in (49), one gets
108
E
1
= E
1
)
= K
1
) +
K
0
J
S
1
w
0
) . (50)
But we have already seen that S
i
/w
0
) =
_
(S
i
/w
0
)dw
0
= 0. Therefore,
E
1
= E
1
) = K
1
) . (51)
Introducing (51) in the left hand side of the second equation of (49) we get
(S
1
/w
0
) as follows
S
1
w
0
=
K
1
) K
1
0
(J)
, (52)
where we used
0
= K
0
/J.
The solution for S
1
can now be found by a direct integration. In general,
once we assume that we already have E
n1
), the procedure to obtain E
n
is
the following
Perform the average on both sides of the nth equation of (49).
Introduce the obtained mean value of E
n
, in the complete equation for
E
n
given by (49) (the one before averaging).
The only remaining unknown S
n
can be obtained by integrating the
following relationship
S
n
w
0
= known function of w
0
and J .
Substitute S
n
in the complete equation for E
n
.
Once all this has been done, the procedure can be repeated for n + 1.
As one can see, getting the energy at a particular order n is possible if and
only if one has obtained the explicit form of S
n1
. On the other hand, S
n
can be obtained only when E
n
has been already found.
The time-independent perturbation theory is very similar to the Rayleigh-
Schroedinger perturbation scheme in wave mechanics, where one can get
109
E
n
only if the wavefunction is known at the n 1 order. Moreover, the
wavefunction of order n can be found only if E
n
has been calculated.
Bibliography
H. Goldstein, Classical Mechanics, Second ed., Spanish version, Edi-
torial Reverte S.A., 1992.
R.A. Matzner and L.C. Shepley, Classical Mechanics, Prentice-Hall
Inc., U.S.A., 1991.
R. Murray Spiegel, Theoretical Mechanics, Spanish version, McGraw-
Hill S.A. de C.V. Mexico, 1976.
L.D. Landau and E.M. Lifshitz, Mechanics, 3rd. ed. Course of Theo-
retical Physics, Volume 1. Pergammon Press, Ltd. 1976.
110
10. ADIABATIC INVARIANTS
Forward: An adiabatic invariant is a function of the parameters and con-
stants of motion of a system, which remains almost constant in the limit in
which the parameters change innitesimally in time, even though in the end
they may change by large amounts.
CONTENTS:
10.1 BRIEF HISTORY
10.1 GENERALITIES
111
10.1 BRIEF HISTORY
The notion of adiabatic invariance goes back to the early years of quantum
mechanics. Around 1910 people studying the emission and absorption of
radiation noticed that the atoms could live in quasi-stable states in which
their energy was alsmost constant. At the 1911 Solvay Congress, the problem
of adiabatic invariance became widely known due to Einstein that draw the
attention of the physics community to the adiabatic invariant E/ of the one-
dimensional pendulum of slowly varying length. He suggested that similar
invariants could occur for atomic systems setting their stability limits. Later,
Ehrenfest was able to nd such adiabatic invariants and their employment
led to the rst quantum approach of Bohr and Sommerfeld.
The method of adiabatic invariants resurged after several decades in the area
of magnetospheric physics. The Scandinavian scientists were especially in-
terested in boreal aurora phenomena, i.e., the motion of electrons and ions
in the terrestrial magnetosphere. One of them, H. Alfven, showed in his
book Cosmic Electrodynamics that under appropriate conditions a certain
combination of dynamical parameters of the charged particles remains con-
stant in the rst order . Apparently, Alfven did not realized that he found
an adiabatic invariant. This was pointed out by L. Landau and E. Lifshitz
who discussed in detail this issue.
10.2 GENERALITIES
We rst show in a simple way which is the adiabatic invariant in the case of
the one-dimensional harmonic oscillator. The employed method will be to
prove that in the limit of the innitesimally slow variation of the parameters,
the adiabatic invariant of the one-dimensional system goes to a quantity that
is exactly conserved in a corresponding two-dimensional system.
Consider a particle at the end of a rope of negligible mass that rotates uni-
formly on a table. Let a be the radius of the circle, m the mass of the
particle, and the angular frequency. If the origin is in the center of the
circle the projection of the motion on the x axis corresponds to the motion
of a simple harmonic oscillator. The motion continues to be harmonic even
when the rope is slowly made shorter by pulling it through a small hole
drilled at the origin. The problem is to nd the invariant quantity in this
case. It is by far more dicult to get the answer for the circular motion
in comparison with the corresponding projected oscillator. Since there are
only central forces even when the rope is slowly shortened the conserved
quantity is the angular momentum l. Therefore, in the slow limit one can
write
112
l = ma
2
. (1)
From this equation one can see that although l is an exact invariant ma
2
is only a slow limit, i.e., adiabatic invariant. More exactly, it is only in the
slow limit that we can claim that the two quantities are equal and we are
allowed to think of the right hand side as an invariant. In the slow limit,
the projected one-dimensional oscillator has a slowly varying amplitude, a
frequency = /2, and a total energy E =
1
2
m
2
a
2
. Therefore, the adia-
batic invariance of the right hand side of (1) implies the adiabatic invariance
of the quotient E/ for the projected oscillator.
A short analytic proof of this fact is the following. One writes the two-
dimensional Hamiltonian
H = (2m)
1
(p
2
x
+ p
2
y
) +
1
2
m
2
(x
2
+ y
2
) . (2)
Hamiltons equations lead to
x +
2
x = o, y +
2
y = 0 . (3)
Changing to polar coordinates we note that l = m
r
2
. The previous Hamil-
ton equations can be written as a single complex equation z +
2
z = 0,
where z = x + iy. We consider a solution of the form z = a expi(t + b),
where a, b, and are real functions that are slowly varying in time. The real
and imaginary parts are x = a cos(t+b) and y = a sin(t+b), respectively.
Since , a, and b vary slowly, l can be written in a very good approximation
as l = ma
2
. The energy of the projected oscillator is E =
1
2
m x
2
+
1
2
m
2
x
2
.
In the slow varying limit of the parameters, we can substitute x in E by the
cosine function to get E =
1
2
m
2
a
2
. Thus, if ma
2
is an adiabatic invariant
then E/ is an adiabatic invariant for the one-dimensional oscillator.
We now clarify the meaning of the term almost constant in the denition of
the adiabatic invariant in the forward to this chapter.
Consider a system of one degree of freedom, which initially is consevative
and periodic and depends on an initially constant parameter a. The slow
variation of this parameter due, for example, to a low amplitude perturba-
tion does not alter the periodic nature of the motion. By a slow variation
we mean the one for which a varies slowly during a period of the motion:
(da/dt) a . (4)
However, even when the variations of a are small during a given period, after
a sucently long time the motion may display large changes.
113
When the parameter a is constant, the system will be described by the
action-angle variables (w
0
, J
0
) and the Hamiltonian H = H(J
0
, a). Assume
now that the generating function of the transformation (q, p) (w
0
, J
0
) is
of the form W
(q, w
0
, a).
When a is allowed to vary in time, (w
0
, J
0
) will still be valid canonical
variables, but W
t
= H(J
0
, a) + a
W
a
. (5)
The second term in (5) can be seen as a perturbation. Then, the temporal
dependence of J
0
is given by
J
0
=
K
w
0
= a
w
0
_
W
a
_
. (6)
Proceeding similarly to the time-dependent perturbation theory, we seek
the variation to rst-order of the mean value of
J
0
during the period of the
unperturbed motion. Since a varies slowly, we can think of it as constant
during this interval. Thus, we can write
J
0
) =
1
a
w
0
_
W
a
_
dt =
a
w
0
_
W
a
_
dt +O( a
2
, a). (7)
One can prove that W
is a periodic function of w
0
, and therefore, can be
written, together with its derivative w.r.t. a, as a Fourier series
W
a
=
k
A
k
e
2ikw
0
. (8)
Substituting (8) in (7) we get
J
0
) =
a
k
2ikA
k
e
2ikw
0
dt +O( a
2
, a). (9)
Since the integrant does not contain any constant term the integral is zero.
This leads to
114
J
0
) = 0 +O( a
2
, a). (10)
Thus,
J
0
) will not display secular variations in the rst order, i.e., in a,
which is one of the basic properties of the adiabatic invariance. In this way,
the term almost constant in the denition of the adiabatic invariant should
be interpreted as constant in the rst order.
Further reading
L. Parker, Adiabatic invariance in simple harmonic motion, Am. J.
Phys. 39 (1971) pp. 24-27.
A.E. Mayo, Evidence for the adiabatic invariance of the black hole
horizon area, Phys. Rev. D58 (1998) 104007 [gr-qc/9805047].
115
11. MECHANICS OF CONTINUOUS SYSTEMS
Forward: All the formulations of mechanics up to now have dealt with
systems having a nite number of degrees of freedom or innitely countable.
However, many physical systems are continuous; for example, an elastic
solid in vibrational motion. Every point of such a solid participates in the
oscillation(s) and the total motion can be described only by specifying the
coordinates of all points. It is not so dicult to modify the previous formu-
lations in order to get a formalism that works for continuous media. The
most direct method consists in approximating the continuum by a collection
of discrete subunits (particles) and then study how the equations of motion
change when one goes to the continuous limit.
CONTENTS:
11.1 Lagrangian formulation: from the discrete to the continuous case
11.2 Lagrangian formulation for the continuous systems
11.3 Hamiltonian formulation and Poisson brackets
11.4 Noethers theorem
116
11.1 Lagrangian formulation: from discrete to the
continuous case
As one of the most simple case for which one can go easily from a discrete
system to a continuous counterpart we consider an innitely long elastic
rod doing longitudinal vibrations, i.e., oscillations of its points along the
rod axis. A system made of discrete particles that may be considered as
a discrete approximation of the rod is an innite chain of material points
separated by the same distance a and connected through identical massless
resorts of elastic constant k (see the gure).
i i i 1 +1
en
quilibrio
desplazado
del
e
equilibrio
We suppose that the material points can move only along the chain. Thus,
the discrete system is an extension of the polyatomic lineal molecule pre-
sented in chapter 6 of Goldsteins textbook. Thus, the equations of motion
of the chain can be obtained by applying the common techniques used in
the study of small oscillations. Denoting by
i
the displacement of the ith
particle with respect to its equilibrium position, the kinetic energy is
T =
1
2
i
m
.
2
i
, (1)
where m is the mass of each particle. The corresponding potential energy
is the sum of the potential energies of each of the resorts:
V =
1
2
i
k (
i+1
i
)
2
. (2)
From the equations (1) and (2) we get the Lagrangian of the system
L = T V =
1
2
i
_
m
.
2
i
k (
i+1
i
)
2
_
, (3)
which can also be written in the form
L =
1
2
i
a
_
m
a
.
2
i
ka
_
i+1
i
a
_
2
_
=
i
aL
i
, (4)
where a is the equilibrium distance between the points. The Euler-Lagrange
equations of motion for the
i
coordinates read
m
a
i
ka
_
i+1
i
a
2
_
+ka
_
i
i1
a
2
_
= 0. (5)
117
The particular form of L in equation (4) and the corresponding equations
of motion have been chosen as being the appropriate ones for passing to the
continuous limit by means of a 0. it is clear that m/a turns into mass per
lenght unit of the continuous system, but the limiting value of ka is not
so obvious. We recall that in the case of an elastic rod for which Hookes
law holds, the enlargement of the rod per unit of length is proportional to
the force (tension) according to
F = Y ,
where is the enlargement per length unit and Y is the Young modulus.
On the other hand, the relative enlargement of the length a of a discrete
system is given by = (
i+1
i
) /a. The corresponding force required to
act on the resort will be
F = k (
i+1
i
) = ka
_
i+1
i
a
_
,
and therefore a should correspond to the Young modulus de Young of the
continuous rod. When one passes from the discrete to the continuous case
the integer index i for the particular material point of the system turns into
the continuous position coordinate x; the discrete variable
i
is replaced by
(x). Moreover, the following quantity
i+1
i
a
=
(x + a) (x)
a
entering L
i
goes obviously to the limit
d
dx
,
when a 0. Finally, the sum over the number of discrete particles turns
into an integral over x, the length of the rod, and the Lagrangian (4) takes
the form
L =
1
2
_
_
2
Y
_
d
dx
_
2
_
dx. (6)
In the limit a 0, the last two terms in the equation of motion (5) are
Lim
a0
Y
a
_
_
d
dx
_
x
_
d
dx
_
xa
_
,
and taking again the limit a 0 the last expression denes the second
derivative of . Thus, the equation of motion for the elastic rod will be
d
2
dt
2
Y
d
2
dx
2
= 0, (7)
which is the known one-dimensional wave equation of propagational velocity
=
. (8)
118
The equation (8) is the well-known formula for the sound velocity (the prop-
agational velocity of the elastic longitudinal waves).
This simple example is sucient to illustrate the main features of the tran-
sition from a discrete to a continuous system. The most important fact that
we have to understand is the role of the coordinate x. It is not a gener-
alized coordinate; it is merely a continuous index substituting the discrete
one i. To any value of x corresponds a generalized coordinate (x). Since
also depends on the continuous variable t, maybe we have to write it in
a more precise way than (x.t) by showing that x, equally to t, can be
viewed as a parameter of the Lagrangian. If the system would have been
3-dimensional and not one-dimensional, the generalized coordinates would
have to be distinguished by three continuous indices x, y, z and write in the
form (x, y, z, t). We note that the quantities x, y, z, t are totally indepen-
dent and show up only in as explicit variables. The derivatives of with
respect to any of them could thus be written as total derivatives without
any ambiguity. The equation (6) also shows the Lagrangian is an integral
over the continuous index x; in the three-dimensional case, the Lagrangian
would write
L =
_ _ _
Ldxdydz, (9)
where L is called the Lagrangian density. In the case of the longitudinal
vibrations of the continuous rod, the Lagrangian density is
L =
1
2
_
_
d
dt
_
2
Y
_
d
dx
_
2
_
, (10)
and corresponds to the continuous limit of L
i
appearing in equation (4).
Thus, it is more the Lagrangian density than the Lagrangian that we used
to describe the motion in this case.
11.2 Lagrangian formulation for continuous systems
We note that in eq. (9) the L for the elastic rod dependes on
= /t, the
spatial derivative of , /x; x and t play a role similar to its parameters.
If besides the interactions between rst neighbours there would have been
some local forces, L would have been a function of . In general, for any
continuous system L can be an explicit function of x and t. Therefore, the
Lagrangian density for any one-dimensional continuous system should be of
the form
L = L
_
,
d
dx
,
d
dt
, x, t
_
. (11)
The total Lagrangian, following eq.(9), will be
L =
_
Ldx,
119
and Hamiltons principle in the limit of the continuous system takes the
form
I =
_
2
1
_
Ldxdt = 0. (12)
From Hamiltons principle for the continuous system one would expect to
get the continuous limit of the equations of motion. For this (as in section
2-2 of Goldstein) we can use a varied path for a convenient integration, by
choosing from a family of functions of depending on a parameter as
follows
(x, t; ) = (x, t; 0) + (x, t) . (13)
where (x, t; 0) is the correct function satisfying Hamiltons principle and
is an arbitrary function of good behaviour that is zero in the extreme
points of t and x. If we consider I as a function of , then in order to be an
extremal solution for the derivative of I with respect to it should become
zero in = 0. By directly deriving I we get
dI
da
=
_
t
2
t
1
_
x
2
x
1
dxdt
_
L
+
L
d
dt
_
d
dt
_
+
_
L
d
dx
_
_
d
dx
_
dt
_
.
(14)
Since the variation of , , should be zero at the end points, by integrating
by parts on x and t we obtain the relationships
_
t
2
t
1
L
d
dt
_
d
dt
_
dt =
_
t
2
t
1
d
dt
_
L
d
dt
_
d
d
dt,
and
_
x
2
x
1
L
d
dx
_
d
dx
_
dx =
_
x
2
x
1
d
dx
_
L
d
dx
_
d
d
dx.
From this, Hamiltons principle could be written as follows
_
t
2
t
1
_
x
2
x
1
dxdt
_
L
d
dt
_
L
d
dt
_
d
dx
_
L
d
dx
__
_
_
0
= 0 . (15)
Since the varied path is arbitrary the expression in the curly brackets is zero:
d
dt
_
L
d
dt
_
+
d
dx
_
L
d
dx
_
= 0. (16)
This equation is precisely the right equation of motion as given by Hamiltons
principle.
In the particular case of longitudinal vibrations along an elastic rod, the
form of the Lagrangian density given by equation (10) shows that
L
d
dt
=
d
dt
,
L
d
dx
= Y
d
dx
,
L
= 0.
120
Thus, as we would have liked, the Euler-Lagrange equation (16) can be
reduced to the equation of motion (7).
The Lagrange formulation that we presented up to now is valid for contin-
uous systems. It can be easily generalized to two, three, and more dimen-
sions. It is convenient to think of a four-dimensional space of coordinates
x
o
= t, x
1
= x, x
2
= y, x
3
= z.
In addition, we introduce the following notation
,
d
dx
;
,j
d
dx
j
;
i,
d
2
i
dx
dx
. (17)
Employing this notation and the four x coordinates the Lagrangian density
(11) takes the form:
L = L(
,
,
,
) . (18)
Thus, the total Lagrangian is an integral extended to all three-dimensional
space:
L =
_
L(dx
i
) . (19)
In the case of Hamiltons principle the integral is extended to a region of
the four-dimensional space
I =
_
L(dx
) = 0, (20)
where the variation of the
(x
; ) =
(x
) + (x
) .
They depend on a single parameter and reduce to
(x
+
L
_
(dx
) = 0. (21)
Integrating by parts equation (21), we get
dI
d
=
_
_
L
d
dx
_
L
,
__
(dx
)+
_
(dx
)
d
dx
_
L
_
= 0,
and taking the limit 0 the previous expression turns into:
_
dI
d
_
0
=
_
(dx
)
_
L
d
dx
_
L
,
__
_
_
0
= 0. (22)
121
Since the variations of each
_
L
,
_
= 0. (23)
The equations (23) are a system of partial dierential equations for the eld
quantities, with as many equations as s are.
Example: Given the Lagrangian density of an acoustical eld
L =
1
2
_
0
.
2
+2P
0
P
0
( )
2
_
.
0
is the equilibrium mass density and P
0
is the equilibrium pressure of the
gas. The rst term of L is the kinetic energy density, while the rest of the
terms represent the change in the potential energy of the gas per volume unit
due to the work done on the gas o por el curso de las contracciones y expan-
siones que son la marca de las vibraciones ac usticas, es el cociente entre
los calores molares a presi on y a volumen constante obtener las ecuaciones
de movimiento.
Solution:
In the four-dimensional notation, the form of the Lagrangian density is
L =
1
2
(
0
i,0
i,0
+ 2P
0
i,i
P
0
i,i
j,j
) . (24)
From the equation (23) the following equations of motion are obtained
j,00
P
0
i,ij
= 0, j = 1, 2, 3. (25)
Coming back to the vectorial notation the equations (25) can be written as
follows
0
d
2
dt
2
P
0
= 0. (26)
Using the fact that the vibrations are of small amplitude the relative varia-
tion of the gas density is given by the relationship
= .
Applying the divergence and using the previous equation we obtain
2
0
P
0
d
2
dt
2
= 0
which is a three-dimensional wave equation, where
=
P
0
0
is the sound velocity in gases.
122
11.3 Hamiltonian formulation and Poisson brackets
11.3.1 Hamiltonian formulation
Hamiltons formulation for continuous systems is similar to that for dis-
crete systems. To show the procedure we go back to the chain of material
points we considered at the beginning of the chapter, where for each
i
one
introduces a canonical momentum
p
i
=
L
i
= a
L
i
i
. (27)
The Hamiltonian of the system will be
H p
i
.
i
L = a
L
i
i
.
i
L, (28)
or
H = a
_
L
i
i
.
i
L
i
_
. (29)
Recalling that in the limit a 0, L L and the sum in the equation (29)
turns into an integral the Hamiltonian takes the form:
H =
_
dx
_
L
L
_
. (30)
The individual canonical momenta p
i
, given by equation (27), go to zero in
the continuous limit, nevertheless we can dene a momentum density that
remains nite:
Lim
a0
p
i
a
=
L
. (31)
The equation (30) has the form of a spatial integral of the Hamiltonian
density H dened as
H =
.
L . (32)
Even when one can introduce a Hamiltonian formulation in this direct way
for classical elds, we should keep in mind the the procedure has to give a
special treatment to the time variable. This is dierent of the Lagrangian
formulation where all the independent variables were considered on the same
foot. This is why the Hamilton method will be treated in a distinct manner.
The obvious way to do a three-dimensional generalizeation of the eld
is
the following.
We dene a canonical momentum
_
x
_
=
L
. (33)
where
(x
i
, t) ,
(x
i
, t) together, dene the phase space of innite dimen-
sions describing the classical eld and its time evolution.
123
Similarly to a discrete system we can seek a conservation theorem for that
look like the corresponding canonical momentum of the discrete systems. If a
given eld
,
= 0
that is
d
dt
d
dx
i
L
,i
= 0 . (34)
Thus, if
=
_
dV
(x
i
, t) .
The generalization for the density (eq. (32)) in the case of the Hamiltonian
density is
H
_
,
,i
,
, x
_
=
L, (35)
where it is assumed that the functional dependence of
.
can be eliminated
by inverting the eqs. (33). From this denition, one gets
H
=
.
=
.
(36)
as a consequence of eq. (33). In a similar way, we obtain
H
=
L
. (37)
Now, using Lagrange equations, eq. (37) turns into
H
=
d
dx
_
L
,
_
=
.
d
dx
i
_
L
,i
_
. (38)
Due to the occurrence of L, we still do not have a useful form. However, by
a similar procedure as used for getting the terms
H
and
H
for
H
,i
we
have
H
,i
=
,i
,i
,i
=
L
,i
. (39)
Thus, by substituting (39) in (38) we get
H
d
dx
i
_
H
,i
_
=
.
. (40)
124
The equations (36) y (40) can be rewritten using a notation closer to the
Hamilton ones for a discrete system. This is possible by employing the
concept of functional derivative
d
dx
i
,i
. (41)
Since H is not a function of
,i
the equations ( 36) and (40) can be written
.
=
H
,
.
=
H
. (42)
Now, by employing the same notation, the Lagrange eqs. (23) can be written
as follows
d
dt
_
L
= 0. (43)
It is fair to say that the almost unique advantage of the functional derivative
is the similarity of the formulas with the discrete case. Moreover, one can
see the parallel treatment of space and time variables.
11.3.2 Poisson brackets
We can get other properties of H by developing the total time derivative of
eq. (35), remembering that
.
es funci on de
,
,j
,
. Thus, we have
that
dH
dt
=
.
d
.
dt
L
d
.
dt
L
,i
d
,i
dt
L
t
.
In this expression, the second and the forth terms nullie each other because
of the denition (33). The derivative simplies to
dH
dt
=
.
,i
d
,i
dt
L
t
. (44)
On the other hand, considering H as a function of
,
,j
,
and
, the
total time derivative is
dH
dt
=
.
+
H
+
H
,i
d
,i
dt
+
H
t
, (45)
where we wrote the expression in such a manner to get an easy comparison
with the second terms of eq. (44), and where using the eqs. (36), (37) and
(39) we obtain
H
t
=
L
t
, (46)
which is an analog to the corresponding one for discrete systems.
125
On the other hand, the equality of the total and partial time derivatives does
not hold. Using Hamiltons equations of motion (eq. (36) yand(40)) and
interchanging the order of derivation, the eq.(45) can be written as follows
dH
dt
=
H
d
dx
i
_
H
,i
_
+
H
,i
d
.
dx
i
+
H
t
.
Now, employing eq. (46) and combining the terms we nally have
dH
dt
=
d
dx
i
_
.
,i
_
+
H
t
, (47)
which is the closest we can approximate to the corresponding equations for
discrete systems.
When L does not depend explicitly on time t, it will not be in H as well.
This implies the existence of a conservative current and consequently the
conservation of an integral quantity, which in this case is
H =
_
HdV . (48)
Thus, if H is not an explicit function of time, the conserved quantity is not
H, but the integral H.
The Hamiltonian is just an example of functions that are volume integrals
of densities. A general formalism can be provided for the time derivatives of
such integral quantities. Consider a given density | and let it be a function
of the coordinates of the phase space (
:
| = |
_
,
,i
,
,i
, x
_
. (49)
The corresponding integral quantity is
U (t) =
_
|dV (50)
where the volume integral is extended to all the space limited by the contour
surface on which
+
|
,i
.
,i
+
|
+
|
,i
.
,i
+
|
t
_
dV . (51)
Let us consider the term
_
dV
|
,i
.
,i
=
_
dV
|
,i
d
.
dx
i
.
Integrating by parts and taking into account the nullity of
,i
.
,i
=
_
dV
.
d
dx
i
_
|
,i
_
.
126
For the term in
.
,i
one uses a similar technique. Substituting the obtained
expressions and grouping appropriately the coecients of
.
and
.
, respec-
tively, and using the functional derivative notation equation (51) is reduced
to
dU
dt
=
_
dV
_
|
+
|
+
|
t
_
. (52)
Finally, introducing the canonical equations of motion (42), we have
dU
dt
=
_
dV
_
|
_
+
_
dV
|
t
. (53)
The rst integral in the rhs corresponds clearly to the Poisson brackets. If
| and J are two functions of density, these considerations allows us to take
as the denition of the Poisson bracket for integral quantities as
[U, W] =
_
dV
_
|
_
. (54)
We dene the partial time derivative of U through the following expression
U
t
=
_
dV
|
t
. (55)
Thus, the eq. (53) podr a could be written
dU
dt
= [U, H] +
U
t
, (56)
which exactly corresponds, in this notation, to the equation for the discrete
systems. Since by denition the Poisson bracket of H with itself is zero, the
eq. (46) turns into
dH
dt
=
H
t
, (57)
which is the integral form of eq. (47). Thus, the Poisson bracket formalism
ocurr as a consequence of the Hamiltonian formulation. However, one cannot
perform a description in terms of Poisson brackets for eld theories by a step
by step correspondence with the discrete case.
Nevertheless there is one way to work out the classical elds which includes
almost all the ingredients of Hamiltons formulation and Poisson brackets of
the discrete case. The basic idea is to replace the continuous spatial variable
or the continuous index by a countable discrete index.
The requirement that be zero at the end points is a contour condition that
might be achieved physically only by xing the rod between two rigid walls.
Then, the amplitud of oscillation can be represented by means of a Fourier
series:
(x) =
n=0
q
n
sin
2n(x x
1
)
2L
. (58)
127
Instead of the continuous index x we have the discrete . We could use this
representation of x only when (x) is a regular function, which happens for
many eld quantities.
We assume that there is only one real eld that can be developed in a
three-dimensional Fourier series
(
r , t) =
1
V
1/2
k=0
q
k
(t) exp
_
i
r
_
(59)
Here,
k is a wave vector that can take only discrete modulus and directions,
so that in only one lineal dimension there will be an integer (or sometimes,
half-integer) wavelengths. We say that
_
.r
dV =
k,k
. (60)
As a matter of fact, the allowed values of k are those for which the condition
(60) is satised, and the coecients q
k
(t) are given by
q
k
(t) =
1
V
1/2
_
e
i
r
(
r , t) dV . (61)
Similarly, for the density of the canonical momentum we have
(
r , t) =
1
V
1/2
k
p
k
(t) e
i
r
(62)
with p
k
(t) dened as
p
k
(t) =
1
V
1/2
_
e
i
r
(
r , t) dV . (63)
Both q
k
and p
k
are integral quantities. Thus, we can look for their Poisson
brackets. Since the exponentials do not contain the elds we have, according
to (54)
_
q
k
, p
k
=
1
V
_
dV e
i
r
_
_
=
1
V
_
dV e
i
r
that is, by equation (60),
_
q
k
, p
k
=
k,k
. (64)
128
From the denition of the Poisson brackets it is obvious that
_
q
k
, q
k
=
_
p
k
, p
k
= 0 . (65)
The time dependence of q
k
is sought starting from
.
q
k
(t) = [q
k
, H] =
1
V
1/2
_
dV e
i
r
_
_
that is
.
q
k
(t) =
1
V
1/2
_
e
i
r
H
. (66)
On the other hand, we have
H
p
k
=
_
dV
H
p
k
(67)
and therefore we get
p
k
=
1
V
1/2
e
i
r
. (68)
Comparando las ecuaciones (67) y (66) tenemos
.
q
k
(t) =
H
p
k
. (69)
In a similar way we can obtain the equation of motion for p
k
.
p
k
=
H
q
k
. (70)
Thus, p
k
y q
k
, obey the Hamilton equations of motion.
11.4 Noethers theorem
We have already seen that the properties of the Lagrangian (or of the Hamil-
tonian) imply the existence of conservative quantities. Thus, if the La-
grangian does not contain explicitly a particular displacement coordinate,
the corresponding caonical momentum is conserved. The absence of an ex-
plicit dependence on a coordinate means that the Lagrangian is not changed
by a transformation that alter the value of that coordinate; we say that it
is invariant or symmetric for that transformation.
The symmetry under a coordinate transformation refers to the eects of an
innitesimal transformation as follows
x
= x
+ x
, (71)
where the variation x
. Noethers
theorem deals with the eect of the transformation of the eld quantities
itselves. Such a transformation can be written
(x
_
x
_
=
_
x
_
+
_
x
_
. (72)
129
Here
_
x
_
is a measure of the eect of the variations of x
and of
. It
can be a function of all the other elds
is a dierent quantity
_
x
_
=
_
x
_
+
_
x
_
. (73)
The characterization of the transformations by means of innitesimal vari-
ations, starting from untransformed quantities means that we consider only
continuous transformations. Therefore, the symmetry under inversion of the
three-dimensional space is not a symmetry of the continuous type to which
Noethers theorem can be applied. As a consequence of the transforma-
tions both in the coordinates and the elds, the Lagrangian will, in general,
appear as a dierent function of the eld and spacetime coordinates:
L
_
_
x
_
,
,
_
x
_
, x
_
L
_
x
_
,
,
_
x
_
, x
_
. (74)
The version of Noethers theorem that we shall present is not of the most
general form possible, but makes easier the proof, without loosing too much
of its generality and the usefulness of the conclusions. We shall suppose the
following three conditions:
1. The spacetime is Euclidean, meaning that the relativity is reduced to
the Minkowski space, which is complex but at.
2. The Lagrangian density is of the same functional form for the trans-
formed quantities as for the original ones, that is
L
_
x
_
,
,
_
x
_
, x
_
= L
_
_
x
_
,
,
_
x
_
, x
_
. (75)
3. The value of the action integral is invariant under the transformation
I
_
dx
_
L
_
x
_
,
,
_
x
_
, x
_
=
_
L
_
_
x
_
,
,
_
x
_
, x
_
.
(76)
Combining the equations (75) and (76) we get the condition
_
_
dx
_
L
_
_
x
_
,
,
_
x
_
, x
L
_
_
x
_
,
,
_
x
_
, x
_
= 0 .
(77)
From the invariance condition, the equation (77) becomes
_
dx
L
_
, x
dx
L
_
, x
_
(78)
=
_
dx
_
L
_
, x
_
L
_
, x
_
_
+
_
s
L() x
dS
= 0 .
Here, L
_
, x
_
is a shorthand notation for the total functional dependence,
S is the three-dimensional surface of the region and x
is the dierence
130
vector between the points of S and the corresponding points of the trans-
formed surface S
dx
_
_
L
_
, x
_
L
_
, x
_
_
+
d
dx
_
L
_
, x
_
x
_
_
. (79)
Now, using the equation (73), the term in the brackets can be written in the
rst-order of approximation as follows
L
_
_
x
_
,
,
_
x
_
, x
_
L
_
_
x
_
,
,
_
x
_
, x
_
=
L
+
L
,
.
Using the Lagrange eld equations
L
_
, x
_
L
_
, x
_
=
d
dx
_
L
_
.
Then, the invariance condition (79) ocurr as
_
_
dx
_ d
dx
_
L
Lx
_
= 0, (80)
which already has the form of an equation for the conservation of a current.
It is useful to develop more the condition giving the form of the innites-
imal transformation as a function of the R innitesimal parameters
r,
r =
1, 2, ..., R, such that the variations of x
be lineal in
r
:
x
=
r
X
r
,
=
r
r
. (81)
By substituting these conditions in eq. (80) we get
_
r
d
dx
__
L
,
L
_
X
r
r
_
_
dx
_
= 0 .
Since the
r
parameters are arbitrary, there are rconservative currents as
solutions of the dierential conservation theorems:
d
dx
__
L
,
L
_
X
r
r
_
= 0 . (82)
The equations (82) are the main conclusion of Noethers theorem, telling
that if the system has symmetry properties fullling the conditions (1) and
(2) for transformations of the type given by the equations (81), then there
exist r conserved quantities.
Further reading
R.D. Kamien, Poisson bracket formulation of nematic polymer dynamics,
cond-mat/9906339 (1999)
131