
Citation 
 Permanent Link:
 https://ufdc.ufl.edu/UF00097926/00001
Material Information
 Title:
 A Study of accidental degeneracy in Hamiltonian mechanics
 Creator:
 Dulock, Victor August, 1939 ( Dissertant )
Broyles, Arthur A. ( Thesis advisor )
Thomas, B. S. ( Reviewer )
McIntosh, Harold V. ( Reviewer )
 Place of Publication:
 Gainesville, Fla.
 Publisher:
 University of Florida
 Publication Date:
 1964
 Copyright Date:
 1964
 Language:
 English
 Physical Description:
 vi, 108 leaves : illus. ; 28 cm.
Subjects
 Subjects / Keywords:
 Coordinate systems ( jstor )
Eigenfunctions ( jstor ) Eigenvalues ( jstor ) Harmonic oscillators ( jstor ) Kinetics ( jstor ) Magnetic fields ( jstor ) Oscillators ( jstor ) Quantum mechanics ( jstor ) Sine function ( jstor ) Symmetry ( jstor ) Dissertations, Academic  Physics  UF Dynamics ( lcsh ) Physics thesis Ph. D
 Genre:
 bibliography ( marcgt )
nonfiction ( marcgt )
Notes
 Abstract:
 One ordinarily defines the accidental degeneracy of a system as
that which does not follow from an analysis of the obvious geometrical
symmetry of the system in configuration space. Frequently, nevertheless
it can be shown that there is a higher order symmetry group in
the phase space of the system, formed with the help of certain hidden
symmetries, which is adequate to account for this accidental degeneracy.
A well known example of a system possessing this property,
and one of the first to be investigated historically, both classically
and quantum mechanically, is the nonrelat i vist ic Kepler problem, whose
quantum mechanical analogue is the problem of a oneelectron atom.
In the quantum mechanical case, the invariance of the Hamiltonian
with respect to the three dimensional rotation group accounts for the
degeneracy in the "m" quantum number, but not for that in the "1"
quantum number.
The purpose of this work is to review the more important results
obtained on the accidental degeneracy of harmonic oscillators and the
Kepler problem and to find some general principles with which the
symmetry groups of Hamiltonian systems can be found. In Chapter 1 a
basic mathematical foundation is laid. This foundation is based on
the fact that all quadratic functions over the phase space of a system
form a Lie algebra. On this basis a method is developed which enables one to find, rather easily, the constants of the motion for systems
with quadratic Hamiltonians or for those which can be made quadratic
through a canonical transformation.
This method is then applied to various systems in the succeeding
two chapters. The harmonic oscillator is treated in Chapter 2, both
in the two and in the n dimensional case. Chapter 3 deals with a
general Hamiltonian which may be interpreted as that of a charged
mass point in a plane harmonic oscillator potential and uniform magnetic
field. Finally, in Chapter 4 the Kepler problem is discussed
and found to have an SU3 symmetry besides the R4generally associated
with it.
All calculations are carried out classically. The results obtained
seem to be valid in the quantum mechanical case except as noted.
The transition to the quantum mechanical version is carried out in the
usual manner by replacing Poisson bracket relations with commutator
brackets and replacing functions by their corresponding quantum mechanical
operators.
 Thesis:
 ThesisUniversity of Florida, 1964.
 Bibliography:
 Bibliography: leaves 105107.
 General Note:
 Manuscript copy.
 General Note:
 Vita.
Record Information
 Source Institution:
 University of Florida
 Holding Location:
 University of Florida
 Rights Management:
 Copyright [name of dissertation author]. Permission granted to the University of Florida to digitize, archive and distribute this item for nonprofit research and educational purposes. Any reuse of this item in excess of fair use or other copyright exemptions requires permission of the copyright holder.
 Resource Identifier:
 022361215 ( AlephBibNum )
13645229 ( OCLC ) ACZ5051 ( NOTIS )

Downloads 
This item has the following downloads:

Full Text 
A STUDY OF ACCIDENTAL DEGENERACY
IN HAMILTONIAN MECHANICS
By
VICTOR AUGUST DULOCK, JR.
A DISSERTATION PRESENTED TO THE GRADUATE COUNCIL OF
THE UNIVERSITY OF FLORIDA
IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE
DEGREE OF DOCTOR OF PHILOSOPHY
UNIVERSITY OF FLORIDA
June, 1964
In memory of John F. Kennedy
ACKNOWLEDGMENTS
The author wishes to express his sincere appreciation to his
supervisory committee and in particular to Professor H.V. Mclntosh,
who directed this research. He would also like to acknowledge the
counsel and guidance given him by Dr. B.S. Thomas throughout his
graduate program.
The author also wishes to thank the Department of Physics and
the Graduate School of the University of Florida for their financial
support during his graduate career.
Finally, the author gratefully acknowledges the help and en
couragement given him by his wife, Maryann.
TABLE OF CONTENTS
Page
ACKNOWLEDGMENTS . . . . . . . . ... .... . ii
LIST OF FIGURES . . . . . . . . ... . . .. i.
INTRODUCTION . . . . . . . . . . . . 1
Purpose . . . . . . . . . 1
Historical Sketch . . . . . . . . . 2
CHAPTER
1. MATHEMATICAL FOUNDATIONS . . . . . . . 7
Introduction . . . . . . . .. . . 7
Lie Algebra Defined by Poisson Brackets . . . 8
2. THE HARMONIC OSCILLATOR . . . . . . .... 20
Introduction .................. . 20
The Plane Isotropic Oscillator . . . . .. 20
The Hopf Mapping . . . . . ..... 24
Canonical Coordinates ............... 29
The Anisotropic Harmonic Oscillator in
Two Dimensions . . . . . . . . ... 31
Commensurability and Constants of the Motion .... 34
The Isotropic Oscillator in n Dimensions . . .. 36
The Hopf Mapping . . . . . . . ... 38
Other Canonical Coordinates . . . . .... 42
The Anisotropic Oscillator in n Dimensions .... 46
3. THE HARMONIC OSCILLATOR IN A UNIFORM MAGNETIC FIELD . 50
Introduction . . . . . . . . . . 50
Classical Zeeman Effect for the Harmonic
Oscillator ..... . . . . . ... . 52
Limiting Cases . . . . . . ..... 55
Canonical Coordinates .. . . . . . . 59
Rotating Coordinates . . . . . . ... 61
Discussion of the Constants of the Motion . . .. 70
Gauge Transformations in Uniform Magnetic Fields . 77
TABLE OF CONTENTSContinued
Page
CHAPTER
4. THE KEPLER PROBLEM ............... ... .. .82
The Two Dimensional Case in Polar Coordinates . 82
The Kepler Problem In Three Dimensions . . . .. 87
The Two Dimensional Kepler Problem In Parabolic
Coordinates ......... ... ........ .. .94
5. SUMMARY .... ................... 98
LIST OF REFERENCES ..................... .. 105
BIOGRAPHY . .. . .. . .. . .. . . .. 108
LIST OF FIGURES
Figure Page
2.1 Gnomonic Projection . . . . . . . ... 26
2.2 Inverse Stereographic Projection . . . . .. 27
3.1 Particle orbit In a strong uniform magnetic field
and harmonic oscillator potential . . . . .. 64
3.2 Particle orbit in a weak uniform magnetic field
and harmonic oscillator potential . . . . .. 65
3.3 Particle orbit in a fixed uniform magnetic field
and harmonic oscillator potential with a large
initial tangential velocity . . . . . ... 66
3.4 Particle orbit in a fixed uniform magnetic field
and harmonic oscillator potential with a small
initial tangential velocity . . . . . ... 67
3.5 Particle orbit shown in Fig. 3.1 as viewed from a
rotating coordinate system . . . . . .... . 71
3.6 Particle orbit shown in Fig. 3.2 as viewed from a
rotating coordinate system . . . . . .... 72
4.1 Three dimensional coordinate system for the Kepler
problem . . . . . . . . ... .. ... .90
INTRODUCTION
Purpose
One ordinarily defines the accidental degeneracy of a system as
that which does not follow from an analysis of the obvious geometrical
symmetry of the system in configuration space. Frequently, neverthe
less it can be shown that there is a higher order symmetry group in
the phase space of the system, formed with the help of certain hidden
symmetries, which is adequate to account for this accidental degen
eracy. A well known example of a system possessing this property,
and one of the first to be investigated historically, both classically
and quantum mechanically, is the nonrelativistic Kepler problem, whose
quantum mechanical analogue is the problem of a oneelectron atom.
In the quantum mechanical case, the invariance of the Hamiltonian
with respect to the three dimensional rotation group accounts for the
degeneracy in the "m" quantum number, but not for that in the "1"
quantum number.
The purpose of this work is to review the more important results
obtained on the accidental degeneracy of harmonic oscillators and the
Kepler problem and to find some general principles with which the
symmetry groups of Hamiltonian systems can be found. In Chapter 1 a
basic mathematical foundation is laid. This foundation is based on
the fact that all quadratic functions over the phase space of a system
form a Lie algebra. On this basis a method is developed which enables
1
one to find, rather easily, the constants of the motion for systems
with quadratic Hamiltonians or for those which can be made quadratic
through a canonical transformation.
This method is then applied to various systems in the succeeding
two chapters. The harmonic oscillator is treated in Chapter 2, both
in the two and in the n dimensional case. Chapter 3 deals with a
general Hamiltonian which may be interpreted as that of a charged
mass point in a plane harmonic oscillator potential and uniform mag
netic field. Finally, in Chapter 4 the Kepler problem is discussed
and found to have an SU3 symmetry besides the R4 generally associated
with it.
All calculations are carried out classically. The results ob
tained seem to be valid in the quantum mechanical case except as noted.
The transition to the quantum mechanical version is carried out in the
usual manner by replacing Poisson bracket relations with commutator
brackets and replacing functions by their corresponding quantum mechani
cal operators.
Historical Sketch
W. Pauli [1] was perhaps the first to associate the accidental de
generacy of the Kepler problem with the less well known second vector
constant of the motion, the Runge vector R, which had been discussed
earlier by W. Lenz [2]. Pauli defines this vector classically as lying
in the plane of the orbit with a direction from the force center to the
aphelion and a length equal to the eccentricity of the classical orbit
and he also gives the quantum mechanical operator definition. The
commutation rules between the components of this vector and the angular
momentum are also stated explicitly. As a consequence of the the com
mutation rules he shows that the system is degenerate in "1" also.
However, L. Hulthen [3] observed the fact that these commutation rules
are the same as those of the generators of the four dimensional rotation
group. 0. Klein was given credit for first recognizing this fact.
V. Fock [4], in 1935, using a stereographic projection in momentum
space, solved the Schridinger equation in momentum space and showed
that the solutions are spherical functions of the four dimensional
sphere and hence that the symmetry group is R4. Later V. Bargmann [5],
in a discussion of Fock's results, showed that Fock's group was gener
ated by the components of the two vector constants R and L, the angular
momentum.
It was about this time that 0. Laporte [6] obtained a differential
operator identical to the Runge vector and presented a method for ob
taining all the eigenfunctions of the hydrogen atom by differentiation,
similar to the method used to find the eigenfunctions of L2 by using
the raising and lowering operators L+ and L_. He made extensive use
of the stereographic parameters, developed by him and G.Y. Rainich [7]
about the same time, in momentum space. The wave functions were found
to be proportional to hyperspherical surface harmonics.
In 1939 and 1940 E.L. Hill and his student, J.M. Jauch [8, 9, 10]
published two articles and a Ph.D. thesis dealing with accidental de
generacy. These authors considered both the two and three dimensional
Kepler problem and the harmonic oscillator in two and three dimensions.
The Kepler problem was found to have R3 and R4 symmetry groups in two
and three dimensions respectively, while the harmonic oscillators had
the symmetry of SU2 and SU3, the groups being generated by the constants
of the motion in all cases. The transformations induced by the gener
ators were found to have an analogy in classical mechanics, namely that
they represented transformations of one orbit of phase space into
another of the same energy. They also showed the existence of a corre
spondence between transformations in classical mechanics and quantum
mechanics.
In 1949, A.W. Saenz [11] wrote a Ph.D. thesis under Laporte, in
which the symmetry groups explaining the accidental degeneracy of vari
ous problems were found by reducing the problems to force free motion
on hyperspherical surfaces of various dimensions. In this manner he
considered the Kepler problem, harmonic oscillator and rigid rotor.
More recently, the n dimensional isotropic oscillator has been
discussed by G.A. Baker [12]. He found the synmetry group to be SUn
and demonstrated that this group could account for all the degeneracy
found. This had been demonstrated earlier by Y.N. Demkov [13] from a
slightly different point of view. Whereas Baker attacked the problem
by finding the most general group which left the Hamiltonian invariant,
Demkov used the generators of infinitesimal transformations which com
muted with the Hamiltonian and then showed that this defined an SUn
group. In a later article Demkov [14] demonstrated the equivalence of
his and Baker's works.
The connection between accidental degeneracy and hidden symmetry
was confirmed by A.P. Alliluev [15] in 1958. He considered two ex
amples, the two dimensional oscillator and the Kepler problem in n
dimensions, which reduced to Fock's result in the case where n = 3.
In the general case, the hidden symmetry of the Kepler problem in n
dimensions was found to be that of the n + 1 dimensional hypersphere.
The one dimensional Kepler problem has been treated by R. Loudon [16]
and was found to be doubly degenerate in all levels except the ground
state, an apparent violation of the theorem that the energy levels of
one dimensional systems are nondegenerate. Loudon discusses his results
in the light of this theorem and shows that it holds only when the po
tential is free of singularities. In 1959, H.V. Mclntosh [17] wrote a
qualitative review of the work up to that time and stated the results.
The most recent article to appear is one by Demkov [18] in which the
concept of excessive symmetry groups is introduced. He gives a pre
scription for finding the minimal symmetry group and as an example
applies it to the anisotropic two dimensional oscillator with incom
mensurable frequencies.
In addition to the works cited above dealing with thestudy of
accidental degeneracy itself, there has been a considerable amount of
work in using the symmetry groups for various problems and applying
them to the many body problem, particularly to the shell model of
the nucleus. One of the first of these was by J.P. Elliott [19, 20]
who undertook to study collective motion in the shell model utilizing
a coupling scheme associated with the degeneracy of the three dimensional
harmonic oscillator. The direct result was a classification of the
states according to the group SU3. Following this treatment was a
series of papers by V. Bargmann and M. Moshinsky [21, 22] on the
group theory of harmonic oscillators, in which they developed a
classification scheme for states of n particles in a common harmonic
oscillator potential. They were primarily concerned with finding a
scheme which would exhibit explicitly the collective nature of the
states. Similar work was also done by S. Goshen and H.J. Lipkin
[23, 24] with the exception that the particles were in a one dimen
sional potential well. This system could also be described from a
collective viewpoint and these states could be classified using the
group of transformations.
In recent times, L.C. Biedenharn has extensively studied the
Kepler problem and its symmetries, both in the nonrelativistic and
relativistic cases. In the nonrelativistic [251 case he has obtained
relations for the representations of R4 exploiting the fact that R4 is
locally isomorphic to R3 x R3. The nonrelativistic problem has also
been considered by Moshinsky [26], where use has been made of the
accidental degeneracy in the coulomb problem to obtain correlated wave
functions for a system of n particles in this potential. In the rela
tivistic case, Biedenharn [27] and Biedenharn and Swamy [28] have ob
tained operators which are analogous to the angular momentum and Runge
vector in the nonrelativistic problem.
CHAPTER 1
MATHEMATICAL FOUNDATIONS
Introduction
In classical and quantum mechanics "quadratic" Hamiltonians play
an exceptional role, a fact which is traceable in the end to the fact
that the Poisson bracket operation satisfies axioms whose expression
is tantamount to saying that the Poisson bracket is an alternating bi
linear functional which is a derivative in each of its arguments.
If H(n) designates the vector space of homogeneous polynomials
of degree n, whose basis may be chosen to be composed of homogeneous
monomials, the Poisson bracket operation maps the cartesian product
H(m) x H(n) into H(m+n2), as will be shown later. The Poisson bracket
operation is a bilinear functional, which means that it defines a linear
transformation when either argument is held fixed, and so it always has
a matrix representation. For the case m = 2, one has a mapping from
H(n) into itself, representable by a square matrix. In this latter case,
the extensive lore of matrix theory is available to help discuss the
transformations. In particular, the mapping may be described in terms
of its eigenvalues and eigenvectors.
Thus the mapping
Tq(p) = (q, p), (1.1)
for a fixed q E H(2), and argument p H(n), can be expected to have
7
specialized properties which become all the more significant when
n = 2. In this case the elements of H(2) form the basis of a Lie
algebra under the Poisson bracket operation. In any event, the gen
eral mapping Tq of Eq. (1.1) yields representations of this algebra
as it acts on polynomials, homogeneous of various degrees as shall
be shown later for the case n = 2. Thus the allusion to the excep
tional role of quadratic Hamiltonians refers to the expectation that
the algebraic structure of the Lie algebra will have its reflection
in the dynamical properties of the system which it represents.
Lie Algebra Defined by Poisson Brackets
The set of 2n linearly independent coordinates and moment of an
arbitrary physical system is said to form its phase space (. These
variables may be paired in such a way that they can be indexed by fi,
i = n, ...1, 1, ...n, where the positive indices indicate coordinates
and the negative ones the corresponding conjugate moment. The Poisson
bracket of any two functions of these variables is defined by
n
(g,h) = 6g 9 (1.2)
1 1 I I
It follows from Eq. (1.2) that the Poisson bracket operation satisfies
the axioms of a Lie algebra [29, 37]:
(f,g) = (g,f) (alternating rule),
(1.3)
(f, ag + Bh) = a{f, g) + B(f, h) (linear), (1.4)
(f, (g, h)) + (g, (h, f)) + (h, (f, g)) = 0 (Jacobi identity).
(1.5)
The Poisson bracket also acts like a derivative in that it satisfies
the relation
if, gh) = g(f, h) + (f, g)h (1.6)
Although all second degree polynomials form a Lie algebra, the
linear monomial functions of the coordinates in phase space itself
also form a linear vector space 0, where the basis may be chosen to
be the 2n fi. Denoting the Poisson bracket operation ( ) as the
symplectic inner product in 4, it is noted that this product does
not satisfy the usual axioms for an inner product since it is anti
symmetric under interchange of the arguments. This difficulty may be
remedied in the following manner. Using the definition of the Poisson
bracket it is found that
1 i = j >0
(fi f = 0 i f j
1 i = j < 0 (1.7)
From this it is observed that if i > 0, the dual vector fi to f may
be taken to be f1. Similarly for i < 0 the dual vector is taken to
be f .. Thus the dual basis to
I
(1.8a)
T = (fn, ... fli fi .... fn)
is then
r = (f ... f1 .fl' fn). (1.8b)
With these definitions the inner product may be defined as
(f, g) = (f, g), (1.9)
where f and g are linear combinations of the fi and
f= = afi (1.10)
i
It is evident then that
(fi' fj) = (fi' fj) = ij (1.11)
and
(f, f) = f, f) = ai2 (1.12)
The space D also has the peculiar property of being selfdual in a
special way, namely the mapping Tf is a linear mapping and is in a one
to one correspondence to the vectors in the dual space D+.
Considering tensor products of @, it is observed that (f, g) also
defines a mapping on these products. If one defines
H(n) = 0 (D (&) .... ( = (e
and if f E H(m) and g E H(n), then (f, g) E H(m+n2). In particular,
when m = 2 i.e., f E Q = H(2), a mapping is obtained which is represent
able by a square matrix.
Such a mapping may be expected to have eigenfunctions and eigen
values. Let h be a fixed element of Q (a quadratic Hamiltonian) and
let gi and Xi be the eigenfunctions and eigenvalues such that
(h, gi) = Th(gi) = Xigi, (1.14)
where the gi are taken to be linear i.e. gi E H(1). If h is represent
able as a normal matrix, the gi form a basis for the space H(1), and
are simply linear combinations of the fi. In particular, since the
gi E D = H(1), they form a basis for all the tensor spaces as well
where all possible products of the gi are taken i.e.,
((gi)) = H(1)
((gig ))= H(2)
((g .... gm)) = H(n). (1.15)
On account of the derivative rule Eq. (1.6) the following relation is
satisfied
(h, gij = gi(h, gj) + (h, gi) gj
= (Xi + Xj) gigj, (1.16)
and hence the monomials in gi are eigenfunctions of h in every H(n)
and form a basis even in the case of degeneracy, for all possible eigen
functions.
In cases where the operator Tf is normal, a complete set of or
thogonal eigenvectors exists which will be denoted gi. Hence, if the
eigenfunctions are indexed such that gi = gi then the eigenvectors of
any normal operator satisfy the rule
(gi, 9j) = i,j. (1.17)
The eigenvalues Xi corresponding to the eigenvectors gi of the normal
operator h occur in negative pairs. This follows from a consideration
of the Jacobi identity
(h,(g ,i gi )) + (g1,' (gi h)) + (gi,(h, gi} = 0,
thus
(h, 1) + (Xi) (g ,i' gi) + i (gi, g_ ) = 0,
and
(Xi + Xi) fgi, gi) = 0, (1.18)
and hence, i. = i since (g., gi) i 0.
Since h H(2), h may be written as
n
h = C ijgigj, (1.19)
i,j=n
i j
where the gi are eigenfunctions of h belonging to eigenvalues Xi and
are indexed so that
Xi = X i (1.20)
Assuming that the eigenvalues of h are distinct and using the identity
(h, h) = 0 one has, with the aid of Eq. (1.19)
0Z
i,j
i=j
z
i,j
i'j
Cij h, gigj)
Cij (xi + .j) gig..
Ii I J J'3 3
(1.21)
Hence each coefficient must vanish and either Cij = 0 or (Xi + Xj) = 0.
However Xi = Xj only if j = i. Therefore
(1.22)
To further evaluate these coefficients recall that
(h, gk) = .kgk,
thus
(1.23)
n
k'k = Cii gigi' 9gk
i=l
n
= ) C [g (9g 9, + 9i i k]9 9 "
i,ii E. 'i 9k + g g k
i=l
Using the fact that the coefficients C.. are symmetric in the
the coefficient of gk on the right side is
xk = Ck,k gk' gk)
(1.24)
indices
(1.25)
h =X Ci,i gigi.
and hence
C (1.26)
cII = gi gi)
The final conclusion is that
h = g i9i (1.27)
(g g;)
i=l
Consider now the case where h may have degenerate eigenvalues.
Using the general form for h, Eq. (1.19), and calculating (h, h) as
before Cij = 0 unless Xi = JX. However in the present instance it is
possible to have Xk = 'm, k # m and hence Xk m= '. Hence
n n
h = C ii gigi Ck,m gkgm, (1.28)
i=l k,m=l
=x
k m
k m
Calculating once again (h, gj) results in the following equations:
n n
Xj9j C g Ck,m gk C m gj} (1.29)
i=l k,m=i
Xk=Xm
k m
k m
n
= i i [gi gi' gj} + gi gi' gj)]
i=l
n
+ Ck,m [gk9m' j)+ mk' gj ]
k,m=1
Xk=
k m
Consider the first term in the second sum. If nmj then kfm implies
k#j and hence Ck, must be zero. Therefore, as before
n
h = ' gigi (1.30)
(9i, 9i)
i=l
It should also be noted that the transformation to eigenvector
coordinates for h is a canonical transformation since the eigenvectors
satisfy (gi, g) = 5i,j, which is a necessary and sufficient condition
for g. and gi to be canonical variables.
Consider now the canonical form of a nonnormal quadratic operator
h, of the Jordan canonical form. Again
n n
h = C Cjgj, (1.31)
i=n j=n
i>_j
where the gi are the principal vectors of h and indexed such that the
same orthogonality conditions hold as before. Using the property that
(h, gi} = gi+1, (1.32)
it is found that
n n
gk+l = Ci
i=n j=l
i>j
z C ij
i=n j=n
ij
Ck+l,k Sk+l
Ck,k+l gk+l
9{gigj 90k
[gigj gki + gj gi' i k
k > 0
k < 0 .
Hence all Cij = 8i,j+_ and
n
h= ggi+l 1
i=l
The sympletic norm also produces certain symmetries and anti
symmetries. For example, the dual mapping R(fi), which maps the
coordinates into moment and the moment into the negatives of the
coordinates,
R(fi) = f
i > 0
R(f i) = fi
preserves the inner product; i.e.
(R(fi). R(fj)) = (fi fj)
On the other hand the mapping P(fi) which simply exchanges coordinates
and moment
(1.33)
(1.34)
(1.35)
(1.36)
P(fi) = fi, (1.37)
or the time reversal mapping T(fi)
T(fi) fi
i > 0 (1.38)
T(fi) = fi
reverses the sign of the symplectic norm
(P(fi)' P(fj)) = (f, fj)
(T(fi), T(fj)) = (fi fj} (1.39)
These operators may be expressed in terms of the eigenvectors, gi of
the normal operator h, previously defined. Recall also that the gi
form a set of canonical coordinates and moment. Using the convention
that positive indices represent coordinates and negative indices the
corresponding canonical moment, these operators have the following
canonical forms:
n
R = ggi +gi9i (1.40)
2(9_ gi
n
i=1
P gig; g ;g ;
=1 2(gi, gi)
n
T = gg' (1.42)
i=1 ^i' gi
These operators have the following effects on
gkTT(
TT 9k) = (T, gk) =
gk
Tp(gk) = Pgk) = gk
TR(g = (R, gk = k
gk
The composite
equivalent to
eigenvectors of h:
k > 0
(1.43)
k < 0,
all k,
k > 0
k < 0
(1.44)
(1.45)
of any two of the operations on an eigenfunction is
operating with the third operator:
T p(TT (g9)) k= = TR(gk)
TR(TTk gk k
TR(TT(gk)) = gk = Tp(gk)
k > 0 TT
k < O = T(g).
k
It is also interesting to note that
following Poisson bracket relations
(P, T)
(R, T)
(P, R)
the three operators satisfy the
among themselves:
= 2R
= 2P
= 2T
(1.49)
(1.50)
(1.51)
The composite of the normal operator h with each of these operators
on the eigenvectors is also worth noting. It is relatively easy to show
gk
_gk
(1.46)
(1.47)
(1.48)
Tp(TR(gk)) =
that
Th(TT(k)) = TT(Th(gk)), (1.52)
and hence the matrices representing h and T commute. However one finds
that P and R anticommute with h, i.e.
T (Th(gk)) = Th(T (gk)) (1.53)
and
T (T (g)) Th (T (g )). (1.54)
Rhk hRk
The existence of either one of these last equations is sufficient to
insure that the eigenvalues appear in negative pairs [30].
In the next two chapters the method described above will be applied
to harmonic oscillator Hamiltonians in two and n dimensions and to vari
ations on the oscillator Hamiltonian.
CHAPTER 2
THE HARMONIC OSCILLATOR
Introduction
The accidental degeneracy and symmetries of the harmonic oscillators
has been discussed by several authors. A particularly thorough account
has been given in the thesis by J.M. Jauch [9], in which he considers the
isotropic oscillator in two, three and n dimensions, and finds the sym
metry group to be SUn, i.e., the Hamiltonian is invariant under transfor
mations induced by the generators of the group. A short discussion of
the anisotropic oscillator in two dimensions with commensurable frequen
cies is also given. The main results of this work are also contained in
the notes of E.L. Hill [33]. More recent discussions of the n dimensional
isotropic oscillator have been given by G.A. Baker [12] and Y.N. Demkov
[13, 14, 18].
The Plane Isotropic Oscillator
In units such that the mass and spring constant are unity, the
Hamiltonian of the plane isotropic harmonic oscillator is
H = (P2 + 2 + + y2)/2 (2.1)
x y
With respect to the basis (x,YPx ,Py), the matrix representation of H
as an operator under the Poisson bracket relation is,
0 0
S= 0 0
1 0
0 1
0
0
0
(2.2)
This representation is found by calculating the effect of H on each
member of the basis under Poisson bracket. With the convention that
x is a column matrix with a 1 in the first row, y with a 1 in the
second row and so on, the matrix representation of H is then easily
constructed. The eigenvalues of the matrix H are
x = i
(2.3)
each root appearing twice.
with X = + i are
a = (1/I2)(Px ix) = 1/J2
The normalized eigenvectors associated
\i
b = (l/f2)(Py iy) = 1/F2
0
(2.4)
(2.4)
and those associated with X = i are
a* = (I/F2)(Px + ix) = 1/W2
0o b* = (l/F2)(Py + iy) = 12 i
(2.5)
where the column matrices are the matrix representations of the
functions on the left. One may easily verify this representation,
for example
(H, a) = (H, Px ix) = i(Px ix) (2.6)
while multiplying the matrix in Eq. (2.2) into the column matrix repre
senting a also gives i times the same column matrix. The quantum me
chanical analogues of these four quantities are the well known raising
and lowering operators for the harmonic oscillator. With these four
quantities ten linearly independent quadratic monomials can be con
structed which are eigenfunctions of H and have as eigenvalues 0, + 2i
or 2i. Only four of these are of interest here, namely those with
eigenvalue zero. They are the four linearly independent constants of
the motion (since they commute with H under Poisson bracket) listed
below:
aa* = (P2 + x2)/2 (2.7)
x
bb = (P2 + y2)/2 (2.8)
y
ab" = [(PxPy + xy) + i(yPx xPy) /2 (2.9)
a*b = [(PxPy + xy) i(yPx xPy)]/2 (2.10)
For the purpose of physical interpretation it is more convenient to
deal with the real and imaginary parts of these constants separately,
as well as to separate the Hamiltonian from them. Accordingly the
following four quantities are introduced:
H = aa* + bb* = (P2 + P2 +2 + y2)/2 (2.11)
x y
D = aa bb = (P2 + x2 P2 y2)/2
x y
(2.12)
L = i(a*b ab*) = yPx xPy (2.13)
K = a*b + ab* = xy + PxPy (2.14)
These four constants have the following simple interpretations.
a) H, the Hamiltonian, is the total energy of the system.
b) D is the energy difference between the two coordinates.
c) L is the angular momentum of the system and the generator of
rotations in the xy plane.
d) K is known as the correlation and is a peculiar feature of the
harmonic oscillator. As a generator of an infinitesimal con
tact transformation, it generates an infinitesimal change in
the eccentricity of the orbital ellipse while preserving the
orientation of the semlaxes, and preserving the sum of the
squares of their lengths. The energy of a harmonic oscillator
depends only on the sum of the squares of the semiaxes of its
orbital ellipse, which remains constant under such a transfor
mation.
The set of functions (K, L, D) is closed under the Poisson bracket
operation. Explicitly, their Poisson bracket table is
K L D
K 0 2D 2L (2.15)
L 2D 0 2K
D 2L 2K 0
Aside from the factor 2 these are the Poisson bracket relations of the
generators of the three dimensional rotation group, or the three com
ponents of the angular momentum in three dimensions. Since the con
figuration space of the two dimensional isotropic oscillator is appar
ently only two dimensional and has only rotations in the xy plane as
an obvious symmetry, the occurrence of the three dimensional rotation
group is rather anomalous.
By performing a series of geometrical transformations on a, a b
and b* one can explicitly demonstrate the spherical symmetry of the
system. This series of transformations is sometimes called the Hopf
mapping.
The Hopf Mapping [31]
Perhaps the best analytic representation of the Hopf mapping is
obtained by introducing polar coordinates (a, r, o, p) for the complex
variables a and b. Explicitly, these are
a = (Xeip cos T)A/2 (2.16)
b = (\eio sin 1)AT2, (2.17)
where X does not denote the eigenvalue of H as previously. The first
step in the Hopf mapping is to form the ratio w defined as the quotient
S= X2 a/b (2.18)
The formation of this ratio may be regarded as a gnomonic projection of
the four dimensional space, regarded as having two complex dimensions,
onto a two dimensional space (Fig. 2.1), having one complex dimension.
This maps the point (a,b) lying on the two dimensional complex circle
of radius X2 into the point
\2 a/b = \2 ei(OP) cot r (2.19)
on the complex line.
The second step is to regard the complex point w = X2 ei(0a) cot
as not lying on a complex line; but rather as a point in a real two
dimensional plane. One now performs an inverse stereographic projection
onto the surface of a three dimensional sphere whose south pole is tan
gent to the plane at its origin. By choosing the origin of the three
dimensional space to be at the center of the sphere, a point on the
surface of the sphere will be specified by giving the radius r, the
azimuth 4 and the colatitude 0. For convenience the sphere is chosen
to have diameter X2. The azimuth is measured in a plane parallel to the
original plane and hence
0 = o p (2.20)
Thus to determine the location of the projected point P on the sphere,
only the colatitude 0 is left to be determined. From Fig. 2.2 one has
tan (A 2) = cot T (2.21)
2 2
but
tan ( = cot (2.22)
hence
0 = 2T .
(2.23)
26
N
3
0
27
0
u
c/C
"4
a
L.
O.
U
t
O
0
4/
(O
C)
1I
LC
;)
I (
c&
o
Thus the coordinates of the point P projected onto the surface of
the sphere from the point w = X2 ei(op) cot T are
0 = 2Tr
S= O p
r = \2/2 .
(2.24)
(2.25)
(2.26)
Expressing the quantities H,
these new coordinates gives
K, L and D in Eq. (2.4) in terms of
H = X2/2
K = (X2/2) sin 0 cos 0 = X
L = (X2/2) sin 0 sin P = Y
D = (X2/2) cos 0 = Z .
Hence, K, L and D simply determine a point on the surface of this sphere.
It should be noted that K generates infinitesimal rotations about X, L
about Y, and D about Z.
Since these three constants of the motion are then directly the
coordinates of the fixed point corresponding to the orbit of the oscil
lator, it is apparent that they generate rotations of the sphere, which
transform one orbit into another of the same energy.
(2.27)
(2.28)
(2.29)
(2.30)
Canonical Coordinates
Whenever two functions satisfy the Poisson bracket relation
(A, B) = LB (2.31)
where 4 is a scalar, a corresponding equation,
(A, (1i ) In B) = 1 (2.32)
may always be written, which is the rule satisfied by a coordinate and
momentum which form a canonically conjugate pair. The first of these
two equations simply states that B is an eigenfunction of A with respect
to the Poisson bracket operation. Consequently it is to be expected
that such eigenfunctions be a ready source of canonical variables. A
system in which the Hamiltonian is one of the moment can be obtained
in the present case. Consider the following diagram in which the two
headed arrows imply that the Poisson bracket of the two functions is
zero and the oneheaded arrows imply that the Poisson bracket is + 1
in the direction indicated:
H = a a* + bb* ( D = a a* bb
(1/4i) In (a*b*/ab) < (1/4i) In (ab*/a*b).
(2.33)
If these four quantities are rewritten in terms of the four parameters
introduced in the Hopf mapping and p + o is denoted by T, then Eq. (2.33)
becomes
H = X2/2 > D = (X2/2) cos 0
SxI
T/2 < 0/2 (2.34)
Hence the Hamiltonian and the energy difference are the two moment
and the angles F/2 and 0/2 are the corresponding canonical coordinates.
Since '/2 is conjugate to H, T increases linearly with time.
The final result of this transformation is to map an entire orbit
into a single point. Since the radius of the image sphere determines
the energy, the only information lost as a consequence of the mapping is
the phase of the point representing the harmonic oscillator in the great
circle comprising its orbit. The lost phase may be recovered by attach
ing a flag to the point representing the orbit, which will then rotate
with constant velocity. The missing angle may then be taken as the
angle which the flag makes with its local meridian. In this way the
motion of the harmonic oscillator in its phase space is shown to be
equivalent to the motion of a rigid spherical rotator, and in fact the
parameters 0, 0 and Y are just the Euler angles describing the orien
tation of the rotor. These same parameters may also be used for the
description of spinors, and a particularly lucid geometrical interpre
tation of this transformation may be found in a paper of Payne [32]
describing various methods of representation of twocomponent spinors.
The Anisotropic Harmonic
Oscillator in Two Dimensions
The anisotropic oscillator, both with commensurable frequencies
and with noncommensurable frequencies has been discussed by Jauch and
Hill [10] and also by Hill [331. These authors were able to show that,
In the case of commensurable frequencies, the symmetry group of the
system is again the three dimensional rotation group. However, they
were unable to come to any conclusion concerning the case with incom
mensurable frequencies. Recently R.L. HudsonI has been investigating
the incommensurable case also.
It can be shown though, that even with noncommensurable frequen
cies the symmetry group is also a three dimensional rotation group.
The proof is sketched below and follows closely that used in the iso
tropic case, except where noted.
If units are chosen such that the frequency in the x coordinate is
unity, the Hamiltonian is
H = (P2+ P2 + x ++ 2y2)/2 (2.35)
x y
where w is the frequency in the y coordinate and is an irrational number.
The eigenfunctions and eigenvalues of H under the Poisson bracket opera
tion are listed below:
IPrivate communication.
Eigenfunction Eigenvalue
a = (l1/2)(Px ix) i
a* = (1I/2)(Px + ix) i
b = (l/A2)(Py iwy) iw
b*= (lJ2)(Py + iwy) i (2.36)
Although the most obvious four constants of the motion are aa*, bb ,
(a) b*, and (a*)Wb, a more convenient set is the following collection of
independent combinations of them:
H = aa* + bb* (2.37)
D = aa* bb* (2.38)
K = [(a)Wb* + (a*) b)/(aa*)(w1)/2 (2.39)
L = i[(a)b* (a*)wb]/(aa*)(w)/2 (2.40)
Again the set (K, L, D) is closed under the Poisson bracket operation
and forms a Lie group, the three dimensional rotation group. The
Poisson bracket table is listed below:
K L D
K 0 2wD 2wL
L 2wD 0 2wK
D 2wL 2wK 0 .(2.41)
A mapping to polar coordinates as in the isotropic case can once
again be performed, however it is not a Hopf mapping. The form of the
final equations for the constants is the same as in Eqs. (2.27) to
(2.30), however, for the anisotropic case
0 = up o (2.42)
A set of functions which proves to be convenient to use as the new
canonical coordinates is diagrammed below with the same convention as
before:
H D
I I
I I(a) wb%*\ 1 (aI*) (a)w
4it (a)(b 4iw (a)b* (2.43)
If the parameters from the mapping are substituted into Eq. (2.43)
this table becomes identical to Eq. (2.34) with the conditions that 0
is defined by Eq. (2.42) and
T = wp + o (2.44)
It is rather interesting to note that the entire effect of the
possible incommensurability of the frequencies in the two coordinates
is completely absorbed by the power mapping (a)W, which is well defined
as exp (w In a) in the case of an irrational w. In this latter case
one has a mapping with infinitely many branches rather than a finite
number as in the case of a rational w. It is of course on this account
that the orbita Lissajous figure in the xy planewould then be open,
and would not correspond to periodic motion.
Commensurability and Constants of the Motion
In reviewing the analysis of the theory of the plane harmonic
oscillator, it is probably appropriate to draw attention to a number
of statements made, both in the literature of classical mechanics,
and the lore of accidental degeneracy. First of all, the existence of
accidental degeneracy is often considered to be allied with the exist
ence of bounded, closed orbits. When this argument is used, it is
applied to the existence of algebraic integrals of the motion, and in
deed it should be noted here that as long as the frequencies of the two
coordinates are commensurable, bounded closed orbits result and the in
tegrals are algebraic, as required by the theory. However, the above
results also verify something which has been generally known [34], namely
that nonalgebraic integrals can exist, and as would be expected, are
associated with spacefilling motion.
Another point regards the separability of the HamiltonJacobi
equation in several coordinate systems. It is pointed out that if a
HamiltonJacobi equation is separable in more than one coordinate system,
it is necessarily degenerate [351. The argument usually given concerns
the fact that orbits with incommensurable frequencies in one coordinate
system will be spacefilling curves, the boundaries of which are co
ordinate arcs. Thus, in whatever coordinate system the motion is de
scribed, the system of separation is uniquely defined as that one form
ing an envelope of the various orbits. Of course, the option left in
the enunciation of this theory is the fact that changes of scale may
still be made, in which the new coordinates are functions each of only
one of the old coordinates, and not several or all. So it is in the
polar mapping, where one forms a ratio of powers of functions of only
a single coordinate. However, this change of "scale" seems to eliminate
the problem of the incommensurable frequencies at the same time, and
thus circumvents the "proof'of uniqueness. This result has only been
secured at a price, namely the inverse of the transformation is infin
itely multiple valued [36]. Otherwise, one could easily reduce the
anisotropic oscillator to a static point on the Riemann sphere, return
from this to an isotropic oscillator, and separate in whatever system
desired.
These considerations are not unrelated to the problem of finding
a Lie group generated by the constants of the motion. In this respect,
a Lie group has always been regarded as being more important quantum
mechanically than classically, because its irreducible representations
would be associated with accidental degeneracies of the Schroedinger
equation. Since the anisotropic harmonic oscillator has no accidental
degeneracies in the case of incommensurable frequencies, it is some
what surprising to find that the unitary unimodular group is still its
symmetry group, at least classically.
In the case of commensurable frequencies, proportional say to m
and n, the operator amb*n commutes with the Hamiltonian, and is inter
preted as creating m quanta in the first coordinate and annihilating n
quanta in the second coordinate. Thus, m quanta in one coordinate are
equal to n quanta in the second coordinate, the energy is restored to
its original level, and the system is left In another eigenstate.
Thus, it is possible to expect and explain degeneracy in the case of
commensurable frequencies.
In the classical incommensurable case, the symmetry group still
exists. After all, the total energy of the oscillator is the sum of
the energies in each coordinate, which is reflected in the fact that
it is the sum of the squares of the semiaxes of the orbit (A Lissajous
figure, like an ellipsoid, has principal axes), and not the individual
semiaxes, which determine the energy. Thus, the dimensions of the
bounding rectangle of a spacefilling Lissajous figure may be changed
subject to this constraint. This, in addition to an adjustment of the
relative phases in the two coordinates, is the effect of the constants
of the motion acting as the generators of infinitesimal canonical trans
formations.
The problem of carrying this result over to quantum mechanics seems
to be in finding a proper analogue of the power functions (a)' and b,
which as ladder operators may not possess fractional powers, even though
they did so when considered as functions of a complex variable.
The Isotropic Oscillator in n Dimensions
Using units such that both the mass and force constant are unity,
the Hamiltonian for this system is
H = (P? + x?)/2. (2.45)
i=1
The eigenfunctions of H under the Poisson bracket operation are
aj = (Pj ixj)/2 (2.46)
with eigenvalues
X = i (2.47)
as well as
a= (Pj + ix.)/T2 (2.48)
with eigenvalues
\ = i (2.49)
Hence one can immediately construct n2 linearly independent constants
of the motion of the form aia', i = 1, ...n and j = 1, ...n, although
as before, it is somewhat more convenient to deal with their real and
imaginary parts separately; in particular,
n
H = aia (2.50)
i=1
Di = ala ai;a i = 2, ...n (2.51)
K.. = a.a + a.a. j > i (2.52)
Ij i J I J
Lj = i[a'aj aiaj] j > i .(2.53)
The latter three sets of constants are closed under Poisson bracket
and satisfy the following commutation relations:
(Kj, Kk) = i[(aam aiam)6k + (aa akaj)5im
(aak aia)m + (a*a ama*) (254)
k jm m j m j ik (2.54)
(Li L il[(aa a al)jk + (aak a.a )
I ma jk + k a)k im
(aak aiak)bjm (ajam aja)bik] (2.55)
(Di, Dj) = 0 (2.56)
(Kij Lkm) = (aam + a.aj k (ak + aia*)jm
+ (aa + a a*) (a a + a a*)bm (2.57)
1m i m jk j k j k im
(Lij, Dk) = (aiak + aj) k (a ak + aia ) k
(a a + aja)5il (2.58)
(D, Kkm) = i[(aak a a)bJm + (aa aa)J
(aa ala ) lk] (2.59)
Depending upon the values of i, j, k, and m all of the relations will
give linear combinations of the functions Di, Kij, or Lij as a result.
The relations also identify the symmetry group of the problem [37],
namely SUn where the n2 1 constants are considered as the generators
of the group.
The Hopf Mapping
In general only 2n1 independent constants of the motion including
2
the Hamiltonian are expected to be found and not n2. Hence there should
exist n2(2n1) = (n1)2 functional relationships among the n2 constants.
As a consequence, the following set of 2nl constants have been chosen
for the Hopf mapping:
the (n1) constants Klj where j = 2, ...n,
the (n1) constants Llj where j = 2, ...n,
and the following linear combination of H and the Di
n
D = BIH + z ;Di (2.60)
i=2
where
2n
1 2 (2.61)
and
B 2 (i = 2, ...n) (2.62)
n
In terms of the ai and ai
n
D = ala aia (2.63)
i=2
The (n1)2 relationships between the constants are
LijKli LIIKIJ
Lij = H +(2.64)
LliLl + KllKl(2.65)
J = H + D
(H+D)2 K 2 L 2
Di = li (2.66)
2(H + D)
n n
where n n
where H2 = D2 + KI + Li (2.67)
i=2 i=2
as may be verified by substitution from the definitions given in Eqs.
(2.50) through (2.53).
For the Hopf mapping let the 2n variables ai and af be redefined as
a1 =i 2(cos T)e"i (2.68)
a = e.(sin T)eij (j = 2, ...n) (2.69)
,J 2 J
where the Ej are direction cosines such that
n
E = 1. (2.70)
j=2
Forming the ratios
 = Ej(cot )e ) (2.71)
aI
gives a complex projective space of nI complex dimensions. Regarding
the space as having 2n2 real dimensions and performing an inverse
stereographic projection a real space of dimension 2n1 is obtained.
As in the case of the two dimensional oscillator the sphere of
2n1 dimensions is considered tangent to the 2n2 dimensional space
at its south pole. By defining 0 as the angle between a radius vector
from the center of the sphere and the diameter between the north and
south poles 2T can be identified with 0. The angles in the planes
41
perpendicular to the northsouth diameter projected from the 2n2 space
parallel to this line are identical with the angles in the space and
are denoted
Pi = 7i p (2.72)
The substitution of Eqs. (2.68) and (2.69) in Eq. (2.63) and into
the equations defining H, Klj, and Lij gives
H = (1/2)%2 (2.73)
D = (1/2)%2 cos 0 = Z (2.74)
Klj = (1/2)X2 sin OEj cos Vj = Xj (2.75)
Llj = (1/2)%2 sin Ocj sin qj = Yj (2.76)
where the Kli, Lli, and D have been identified with the 2nl coordi
nates in this space. As in the case for n = 2 these 2n1 constants
determine a point on the surface of a sphere of radius H in this space.
At the same time the remaining constants still define infinitesimal
transformations in the respective coordinate planes such as rotations
or changes in the correlation or relative phases. It should be noted
that only in three dimensions is it possible to identify planes and
their perpendicular coordinates in a one to one fashion. In general
for n dimensions there are (1/2)n(n1) coordinate planes.
I_ ~XXI I
Other Canonical Coordinates
In finding canonical sets of coordinates and moment a wide choice
of momenta,even when the n moment are selected from amongst the con
stants of the motion,is available. The final decision will depend upon
the particular application desired. However, there are some simple
linear combinations which are often sufficiently useful, to make their
display worthwhile.
As a first selection for the canonical moment consider the follow
n
ing:
P1 = H =) ala (2.77)
i=2
j1 n
Pj a= a aia (j = 2, ...n), (2.78)
i=l i=j
where P. is the energy difference between the first (jl) coordinates
J
and the last n(jl) coordinates. The coordinates conjugate to these
moment must then have the form
1 ala
S= In n (2.79)
Q1 = ala n
I a la
Q. = In a laj (j = 2, ...n) (2.80)
J 4i a. a
In terms of the parameters used in the Hopf mapping these become
H = X2/2
(2.81)
P2 (X2/2) cos 0
Pj =
(X2/2) [cos2 T
j1
+i.
i=2
j1
= (%2/2) [cos a + Efsin2 0/2]
i=2
2 2 2 ]
eisl n ]
n
E Sin T 
i=j
(j = 3, ...n) (2.83)
Ql = 1/2(p on)
Q2 = 1/2(p 02)
Qj = 1/2(ao. 0.)
J J1 J
(j = 3, ...n).
One might also consider taking the total energy and the energy in the
last n(jl) coordinates for the canonical moment,
n
PI = H = aia'
i=l
(2.87)
n
PJ = a;a
(j = 2, ...n)
with the corresponding coordinates
1 a.a.
S naj a2+
2: aa+
(j = 1, ...n1)
1 a
n
S= In 
2i a,
and
(2.84)
(2.85)
(2.86)
(2.88)
(2.89)
(2.90)
(2.82)
which in terms of the Hopf parameters become
H = X2/2
n
P = (2/2)Z ? sin2 1
i=j
Q = p 02
Qj = oj aj+
Qn= on
(j = 2, ...n)
(j = 2, .. .n1)
Finally, one might consider using the energy and(nl) successive
energy differences,
n
PI = H aiaYi
i=l
P = .a j = (2, ...n)
The corresponding coordinates are
S a1a2 ..an
Q In n
2ni ala2 ...an
I k1 nk+1 n
Qk In ( (1) ) ( i
2ni j=l aj j=k
kI
a
J
(2.96)
(2.97)
(2.98)
(2.99)
k = (2, ...n).
(2.91)
(2.92)
(2.93)
(2.94)
(2.95)
In terms of the Hopf coordinates
P, = H = X2/2 (2.100)
P2 = (2/2)(cos2 T E2 sin2T) (2.101)
P = (\2/2) sin2 E e:) (j=3 ...n) (2.102)
n
Q = (l/n)(p + oi) (2.103)
i=2
and J
Qj= (/n)[(nj+l)(p + aO)(j1) ok] (2.104)
k=2 k=j
In all three cases, the transformation has been such as to give
the actionangle variables. The effect has been to map the orbit into
a point on the surface of a sphere of 2nl dimensions whose radius is
the energy H. Since the mapping is many to one, one sees that as in
the two dimensional case an entire orbit is mapped into a single point.
The orbital phase may still be recovered by attaching a flag to the
point representing the orbit which will then rotate with constant
angular velocity Q1, since Q, is conjugate to H and hence increases
linearly with time.
The second case also has another interesting property, namely D2
is the Hamiltonian for an isotropic oscillator in (n1) dimensions,
and the other Di are constants of its motion. Hence one may perform
another Hopf mapping and continue in this manner until one has the 2
dimensional oscillator again.
In general, any of the 2n1 independent constants from the n21
constants could be chosen as generators, and used as coordinates. It
should be noted that the expressions for the coordinates could just
as well be expressed in terms of an arctangent by using the well known
relation
tanI x = i n ix (2.105)
2 l+ix
which would exhibit the similarity between these coordinates and those
obtained by Goshen and Lipkin [23].
The Anisotropic Oscillator in n Dimensions
In units where the mass is unity, the Hamiltonian for the n
dimensional anisotropic harmonic oscillator is
n
H = (P + xi), (2.106)
i=1
where wi is the frequency in the ith coordinate. The eigenvectors are
all of the form
a = (Pj i w. xj)AF2 (2.107)
a = (Pj + i ~j xj)/2 (2.108)
with eigenvalues
x. = i J (2.109)
J J
There also exists n2 constants similar to those obtained in the iso
tropic case, namely
n
H aial
i=1
Di = aia ai+1 aTi+ (i = 1, ...n1)
(a J a' + ai J a ))
Kij w.l uil
(aiaf)2 (aj ay) 2
i[a. J a. a. J a ']
Lij i 1
(aia ) 2 (aja) 2
(2.110)
(2.111)
(2.112)
(2.113)
Once again this set is closed under Poisson bracket and satisfies
the commutation relations of the generators of SUn*
It should be noted that every normal quadratic Hamiltonian is
equivalent to some anisotropic harmonic oscillator, because they can
be put in the canonical form given in Chapter 1.
For nonnormal Hamiltonians such as the free particle
1 2
H = (2.114)
2m
there is no equivalence with the harmonic oscillator. However, symmetry
groups still exist for these cases; for example the above expression is
invariant under rotations. The symmetry of such Hamiltonians has been
discussed in considerable detail by Moshinsky [38] and F.T. Smith [39,
40].
For the Hamlltonian given In Eq. (2.114) (where m 1)
ing set of nine quadratic constants of the motion exists,
H = (P2 + p2 + p2)/2
x y z
2 2
A = P2 P2
A2 z2
=PP
A3 = PXPy
A4 = PxPz
A5 PyPz
L1 = yPz zPy
L2 = zPx xPz
L3 = xPy yPz
The relations between these constants are summarized in the
Poisson bracket table:
the follow
(2.115)
(2.116)
(2.117)
(2.118)
(2.119)
(2.120)
(2.121)
(2.122)
(2.123)
following
H Al A2 A3 A4 A5 L1 L2 L3
0
0
0
0
0
0
2A5
2A4
4A3
0
0
0
0
0
0
4A5
2A4
2A3
0
0
0
0
0
0
A3
AI+A2
A5
0
0
0
0
0
0
A2
A3
A4
0
2A5
4A5
A4
A3
A2
0
L3
L2
0 0
2A4 4A3
2A4 2A3
A5 Al
AIA2 A5
A3 A4
L3 L2
0 Li
Ll 0
(2.124)
These define a group which is the
dimensional rotation group and the
semidirect product (41] of the three
five dimensional abellan group formed
by the five symmetric second rank tensors. In general, a free particle
in n dimensions will have n21 linearly independent constants of which
n(nl)/2 will be the generators of Rn. The others, n(n+l)/2 1 In
number, commute among themselves and are the generators of the abellan
group. These constants are second rank tensors formed by taking pro
ducts of the moment.
The harmonic oscillator Hamiltonlan will be generalized in the
next chapter and consideration given to the problem of a plane harmonic
oscillator in a uniform magnetic field.
CHAPTER 3
THE HARMONIC OSCILLATOR IN A UNIFORM MAGNETIC FIELD
Introduction
In the previous chapter it has been shown how the accidental de
generacy of the harmonic oscillator may be accounted for by the methods
presented in Chapter 1. Although historically, attention has most often
been focused upon the accidental degeneracies of the harmonic oscillator
[9] and the Kepler problem [4, 51, the problem of the cyclotron motion
of a charged particle in a uniform magnetic field possesses certain
unique features, among which is that it has two linear constants of the
motion. It also has accidental degeneracies of Its own [42]. Its
Hamiltonlan,
H = ( A)2/2m, (3.1)
where, for one choice of gauge
A = ( By, + BoX, 0) (3.2)
Is quadratic.
As shown In the previous chapter, the constants of the motion of
the two dimensional harmonic oscillator K, L, and D arose in a natural
fashion as products of linear elgenfunctions of the operator TH. The
eigenvalues of TH occurred in negative pairs, so that their products,
belonging to the sums of the corresponding eigenvalues, were constants
50
of the motion. Moreover, there was a convenient mapping, the Hopf mapping,
which yielded a set of canonical angular coordinates 4 and T, whose con
jugate moment were respectively, H and D.
In the present instance, one finds by expanding Eq. (3.1) and using
the gauge In Eq. (3.2) that
1 e28~2 eBo
H (P2 2) o (x2 + y2) + (yx xPy). (3.3)
2m 8mc2 2mc
Effectively this is the sum of a harmonic oscillator Hamiltonlan and a
term proportional to the angular momentum.
There is one main difference between the problem of cyclotron motion
and the harmonic oscillator. The representation of the harmonic oscil
lator Hamlltonian, applied to the homogeneous space H(1) has negative
pairs of nonzero elgenvalues. All constants of the motion are generated
by products of the corresponding eigenfunctions, for which the eigenvalue
sum Is zero. Thus the lowest order constants of the motion are quadratic,
belonging to H(2), and all others are linear combinations of these pro
ducts. With the cyclotron Hamiltonlan there are zero eigenvalues, and
therefore there are elements of H(l) which are constants of the motion.
Whereas the Poisson bracket of two quadratic constants Is quadratic, the
Polsson bracket of two linear constants is a scalar. A rather different
kind of symmetry group results in the two cases; for the harmonic oscil
lator one finds the unitary unimodular group, while for cyclotron motion
one obtains a group generated by the harmonic oscillator ladder operators.
Classical Zeeman Effect for the Harmonic Oscillator
It is somewhat more instructive to consider a plane isotropic har
monic oscillator in a uniform magnetic field rather than the cyclotron
motion exclusively. Its Hamiltonian,
(Pe 2 1 22
H = ( A) + m r
= I (p2 + p2) + (w + ()(x2 + y2) + (yx xP)
Tm x 2 '
(3.4)
Is also quadratic. In Eq. (3.4) w is the Larmor frequency eB/2mc and
o is the natural frequency of the oscillator. This Hamiltonian reduces,
in the limit as w +0 (B .0) to that of the plane harmonic oscillator,
while as (. i0 It reduces to that of cyclotron motion.
The matrix representation of this Hamlltonian considered as an
operator under Poisson bracket and denoted by TH Is
0 W m((w2 + 2) 0
( 0 0 m(w2 +
TH 1 0 0
m
0 J 0 (3.5)
m
where the basis is composed of the monomials (x, y, Px' Py). The
elgenvalues and elgenvectors of TH are
Elgenvector Eigenvalue
u i[(w2 + w2)1/2 +
u* il[(2 )12 + / ]
v 1[(w2 + 2)1/2 w]
v* i[(12 + 12 ] (3.6)
By defining
r = x ly (3.7)
SP = Px iPy (3.8)
the eigenvectors u and v can be written as follows:
u [m(w2 + w2)]1/2 r. + P+ (3.9)
v = [m(w2 + w2)]1/2 r + P (3.10)
where u* and v* are simply the complex conjugates of u and v respec
tively. These four elgenfunctions satisfy the following relation
(u*, u) = (v*, v) = 41(2 + w2)1/2 (3.11)
The constants of the motion will be products of eigenfunctlons,
the sum of whose eigenvalues is zero. Hence one establishes the
quantities uu*, vv*, u*Rv, and uRv* as constants of the motion where
R is a number such that
R\I = \2 ,
(3.12)
and
\ = (w2 + w2)1/2 + W
\2 ( + U )1/2 .
In order to
following linear
display the symmetry group in a convenient
combinations are taken as the constants of
H = [l uu* + 2 vv*]/4/2 + 02
Rl
K [uRv* + u*Rv]/R1/2(uu*)
Rl
L = i[uRv* u*Rv]/R1/2 (uu*)f
D = [uu* Rvv*]/R ,
where H has been written in accordance with
quantities satisfy the following relations:
Eq. (1.29). These four
TH(K) TH(L) = TH(D) = 0
TK(L) = aD
TL(D) = aK
TD(K) = CZL,
(3.19)
(3.20)
(3.21)
(3.22)
where
where a 8/2+ j (3.23)
The first of these equations is simply a statement of the fact that K,
L, and D are constants of the motion, while the latter three equations
show that the symmetry group of the system is SU2.
(3.13)
(3.14)
form the
the motion:
(3.15)
(3.16)
(3.17)
(3.18)
If in Eq. (3.12) R is a rational number, then the classical system
has bounded closed orbits and a quantum mechanical analog exists for
the operators K and L. However if R happens to be irrational, then the
orbits are space filling.
Limiting Cases
In taking the limit as w +0, i.e. as the magnetic field is turned
off, Eq. (3.4) becomes the Hamiltonian for the two dimensional Isotropic
harmonic oscillator, the elgenvectors become,
u Nm wo r + mP (3.24)
v m Wo r" + P" (3.25)
.vm
while the four eigenvalues degenerate Into two, namely iwo, whence
R = 1. Expressed in terms of the new u and v, the four constants are
H (uu* + vv*)/4 (3.26)
K uv* + u*v (3.27)
L i(uv* u*v) (3.28)
D = uu* vv* (3.29)
which are the constants previously obtained for the Isotropic oscil
lator. The commutation rules for K, L, and D also still hold.
In considering the limit as w~ ~0 one finds that
u =Nm u r*+ + P+
v = m m r" + P~
J^m
and
(3.30)
(3.31)
The eigenvalues of u and u* approach 21w in this limit while those
belonging to v and v* both approach zero. Hence, in order to satisfy
Eq. (3.12) R must also approach zero. The Hamiltonian In this limit
approaches that for pure cyclotron motion as in Eq. (3.3). From the
values of the eigenvalues one immediately has two linear constants of
the motion, v and v*, and one quadratic constant uu*.
The constants can be explicitly derived by considering the com
mutation rules of K, L, and D in the limit as wo 0.
Rewriting the commutation relations explicitly gives
SuRy* + u*Rv
R '
Ri/2(uu*)
i uRv* u*Rv
Rl '
R1/2 T) 2
uRv* u*Rv
Rl
R1/2 (uu 2
uu* Rvv*
R
uu* Rvv*,(3.32)
,(3.32)
R
= uRv* + u*Rv (3.33)
Rl
Ri/2(uu*) 2
Suu* Rvv*,
R
uRv* + u*Rv
RL
RI/2(uu ) 2
i uRv* u*Rv (3.34)
R2 R1(
R1/2(uu ) 2
Multiplying the first of these by R and the latter two by R3/2 and then
and
taking the limit as wo 0 (R +0), results in the equations
I ((uu*)1/2 (v*+v), (uu*)1/2 (v*v)) = 8 uu* (3.35)
I ((uu*)1/2 (v*v), uu*) = 0 (3.36)
(uu*, (uu*)1/2 (v* + v)) = 0 (3.37)
Since uu* Is simply twice the cyclotron Hamiltonlan, it follows from
the last two equations that both the real and imaginary parts of v are
constants of the motion. Dividing Eq. (3.35) by uu* gives
(2 v = 8w (3.38)
Expanding Eq. (3.1) results in
H = (P2 + P)/2m + mw2 (x2 + y2)/2 + (yPx xP ).
x y y ( 3 3 9)
(3.39)
Thus the cyclotron Hamiltonlan splits into two parts, one being the
harmonic oscillator Hamiltonlan Ho and the second proportional to the
Zcomponent of the angular momentum L. Both of these terms commute
with the total Hamiltonian H = Ho + L. Hence another quadratic con
stant of the motion Is
D = Ho L (3.40)
For convenience in notation define
S m (vv*) my Px (3.41)
Q Sm (vv*) mW + Py (3.42)
With these definitions the following commutation rules hold:
(H, D) = (H, S) (H, Q) J 0 (3.43)
(S, Q)} 2mw (3.44)
(D, S)} 2wQ (3.45)
(D, Q) = 2wS (3.46)
These commutation relations coincide with those obtained by Johnson
and Lippmann [42]. These authors have discussed the two constants S
and Q in considerable detail and show that they are simply related to
the location of the center of the circular orbit and to its diameter,
which can readily be seen In the following manner. The canonical
momentum expressed in terms of the mechanical momentum is
S= my + E (3.47)
c
Substituting for the canonical moment in Eqs. (3.41) and (3.42) gives
S = m(2wy vx) (3.48)
Q m(2w + vy) (3.49)
Evaluating S when vx = 0 and Q when vy = 0 gives the center of the
orbit as
(xc, Yc) = (Q/2mw, S/2mw). (3.50)
These constants also determine the diameter of the orbit. Since the
orbit is a circle only one of the constants must be considered. For
example, consider Q. When y takes on its maximum positive value, x
takes on its minimum value and when vy takes on Its maximum negative
value, x is a maximum and hence the diameter d is
d = xmax Xmin m (3.51)
w
Because of the continuum of points available for the center for a
given energy the degeneracy of this problem is infinite.
Canonical Coordinates
As in the case of the plane harmonic oscillator a set of canonical
coordinates can be found such that the Hamiltonian becomes a canonical
momentum. In fact, two moment for the problem are:
S= uu* + X2 vv*
H = (3.52)
4( 2 2 1/2
and
D uu* Rvv* (3.53)
R
while the coordinates conjugate to these moment are
Ql 2[ln + %i In v]/4i Xl,2 (3.54)
u v
Q2 2[In u* X In v]/16i ( 2 + 2)1//2 (3.55)
Q %2 u I v X2 (w +0) (.5
respectively.
A mapping similar to the Hopf mapping may also be performed where
u = [% /2 cos T e"p] (3.56)
and
v = [l/2 sin T eI ] (3.57)
Under this mapping the moment become
H = k/4 (2 + 2) 1/2 (3.58)
D = %I cos 9 (3.59)
where
O 27 (3.60)
and the coordinates are
QI = Y/2X2 (3.61)
and
and 2 2 1/2
Q2 = */8R(2 + W /2, (3.62)
where
T = Rp + o (3.63)
and
= Rp a (3.64)
Performing the mapping on the other two constants of the motion gives
K = Xi sin 0 cos 0 (3.65)
and
L = Xs sin 0 sin 0. (3.66)
There also exists another set of canonical coordinates which were
originally defined by Goshen and Llpkin [24]. Written In terms of the
cartesian coordinates the moment are taken to be
P = [(p + P2)/2m + m(W2 + 2)(x2 + y2)/2]/(w2 + 2 .67)
q x y 0 0
P = XPy yP .
The corresponding coordinates are
o = tan'
2
1
q = I tan1 2
2 (2
[P P /m + m(i2 + w2) xy]
2 2y 2 2 (3.69)
(Px + P2)/2m + m(u2 + o2)(x + y2)/2
(o2 + 2)1/2 (xP + ypy)
P)/2 m( T + )(x +
In terms of these variables the Hamiltonlan has the particularly simple
form
H = (WZ + U) I/2 Wp0
(3.71)
from which It follows that both Pq and PO are constant In time.
Rotating Coordinates
The problem of the plane harmonic oscillator in a uniform magnetic
field has a certain uniqueness when viewed from a rotating coordinate
system. However, the problem will first be solved In plane polar
coordinates. An extremely lucid description and tabulation of these
orbits has been given by E.R. Harrison [431.
Assuming the direction of the magnetic field to be in the negative
Z direction, the Lagranglan is
and
(3.68)
and
2)/2.(3.70)
y2)/2
I
m 112
L = (r2 + r262) m r2 mwr2b
2 2
where the following gauge has been chosen,
Ar = Az = 0
and
AO = Bor
The moment conjugate to r and 6 are
Pr = mr
and
SP = mr2( w).
Since L Is cyclic in ,. Pn is a constant of the motion.
(3.72)
(3.73)
(3.74)
(3.75)
(3.76)
The Hamiltonian
is
H Prr
2
= (Po
and the equations of
+ P,9 L
+ Pg/r2)/2m +
motion are:
2
P PO
mr3
PO = 0
= Pr/m
6 +P
mr
In general, the effect of imposing
m(w2 + w2)r2/2 + uPo
S2 2
m(W + +w)r
a uniform magnetic field on a system
(3.77)
(3.78)
(3.79)
(3.80)
(3.81)
with a central potential Is to add two terms to the Hamiltonian, namely
a harmonic oscillator potential which Is often neglected for small
fields [44], and a term proportional to Pg, the angular momentum.
The solutions for the orbit equations are:
0 0. a sin1 br22a2 )+ sin1 [ 2W2r2b )
21aI r2[b2_4a2W21/2 2W [b24a2W2]1/2
a br2 22 w 2W2r2 b
 sin' sinl
21a r[b224a2 1/2 2W [b24aW2 1/2
(3.82)
where 2 22 2
b = 2 2r + (3.83)
r2
0
a = Pe/m (3.84)
W = ( +2 )/2 (3.85)
The subscripts on the coordinates and velocities denote Initial
values. These orbits are plotted In Fig.3.1 through 3.4. In all cases
units have been chosen such that m q c = 1. Fig. 3.1 and 3.2
show the high and low field orbits respectively, for the same set of
Initial conditions and wo. Fig. 3.3 and 3.4 show the orbits for a
fixed magnetic field, but for different values of the initial tan
gential velocity. In general the orbit will be a processing ellipse
for PO > 0 and will be a hypotrachold similar to Fig. 3.1 for Pg < 0.
64
Parameters Initial Conditions
= 63/4 r = 1 r 0
o o
w 4 =0 6 =1
o o o
0=0
e 1
Fig. 3.1 Particle orbit in a strong uniform magnetic field
and harmonic oscillator potential.
E
U 3 0
0. 3
O
0
. I
 0
U I
0 0
IL C
O
0
0.
41
u
0
o
u
0
U
41
0
L
0
L.
t
C
0,
(U
u
0
6I
E
0
U
*i
0>
L'
H
0
u >
4
c 
0U
EO
L
0 >
o 
Cr0
10
o .
U
c
*J
0
E 
 C
CO
U
LU7
(gIn
I
0
o<
Cu
U
4
0)
4 E
C
u
CU
E 
04
Sc4
 0
.00
OL.
0&
 i.
U
4 O
w
C)
3 HI
OUJ
IU *
u0
^v
(U <
 m
In a rotating coordinate system defined by
S= 0 wt (3.86)
S= r (3.87)
where w is the Larmor frequency, the Lagrangian is
m 2 2 2 2 m(o 2
L = [r + r ( r (3.88)
2 2
The canonical moment in the rotating coordinates system are:
PF = mr (3.89)
P = m;2 (3.90)
e
L is also cyclic in 8 and hence P is a constant. The Hamiltonian in
the rotating coordinate system is
H PrF+ P 0 L
0
S(Pr2 + P2 )/2m + m(2 + 2)r /2 (3.91)
and hence the equations of motion are
= m(w2 + 2)? (3.92)
r mF3 o
P = 0 (3.93)
r Pp/m (3.94)
P
S=  (3.95)
m:2
which are the equations of motion for a plane Isotropic harmonic oscil
lator with force constant m(w2 + w(). In general the transformation of
a Hamiltonian with a central force from a stationary coordinate system
to a rotating coordinate system will subtract a term proportional to PO.
Generalizing then, It can be stated that a Hamiltonian with a central
potential and with a uniform magnetic field present, when viewed from
a rotating coordinate system, has a harmonic oscillator potential added,
that Is
H T + V(r) + 1 mw2 2 (3.96)
2
where T is the kinetic energy In the rotating coordinate system and w
is the rotation frequency of the coordinate system with respect to the
fixed system. In the present case since V(r) is a harmonic oscillator
potential one simply has a harmonic oscillator with a larger force con
stant.
The solution for the orbit defined by the Eqs. (3.92) through
(3.95) is
a b a2 a b2a2a2
60 sin] (n si) sin) ( o
21al F2[b2 4a2 2 2 1/2 2a s 1 2b24a 2 2 1/2
(3.97)
where W is defined in Eq. (3.85) and
a = P5/m (3.98)
b a2/F2 + w2 F2 + (3.99)
Equation (3.97) Is the equation for an ellipse, regardless of the
magnitude of w. The orbits corresponding to Fig. 3.1 and 3.2 in a
rotating coordinate system are shown in Fig. 3.5 and 3.6.
Discussion of the Constants of the Motion
Of the four constants of the motion defined in Eqs. (3.15) through
(3.18) one is H, the energy of the system. Another constant is the
angular momentum PO which is a linear combination of H and D
PO = H + R (3.100)
XI2 4(w2 2o)1/2
Since only three of the constants are independent one would like to
discover one other constant, independent of the two above, which has
a physical meaning. This can be done by considering the constants of
the motion in the rotating coordinate system and then transforming back
to the original coordinate system.
In the rotating coordinate system two constants of the motion are
aa* and bb*, where
a = PR i(2 + o)1/2 x (3.101)
and
b P 1(w2 + ) 2)1/2 (3.102)
and where the mass is taken to be unity. By a simple rotation the axis
can be oriented so that the semimajor axis of the ellipse lies along
x and the semiminor axis lies along y. In this case then aa* Is pro
portiona to the maximum value which 2 attains, because since
portional to the maximum value which x attains, because since
71
o 44
0o
u
CC
o
0
o
0
U
E
0
L
4
o
3
0,
C
o
r)
L
tjE
72
o
0
E
U)
4)
O
L
0
0
U
0
'.
lu
r
E
O
0
L
/I I 8
,
,10
3
Lo
0
u
r
a
aa* P2 + (w2 + (2) g2
x o
is a constant, it can be evaluated when P. = 0. Similarly bb* Is pro
x
portional to y max. Alternatively, one can say that aa and bb' give
the boundaries of the orbit i.e. aa* gives the maximum radial distance
the particle can attain and bb* is the minimum radial distance.
These two radii are also constants of the motion In the stationary
system since the transformation only involved the angle. However, it
can explicitly be shown that these are constants.
By transforming from the rotating coordinate system to the station
ary one aa* can be written in terms of the eigenvectors and their con
jugates as
aa* =(1/4)(uu* + vv* + uv* e2iWt + u*v e21t) .
In the stationary system the time rate of change is
(3.104)
d (aa*) = (aa*, H) + (aa*)
dt /t
=(1/4)[21wuxjve21wt + 21 u*ve21wt
+21iuv*e2iwt 21iu*ve21tt] = 0,
(3. 105)
and hence aa* is a constant of the motion in the stationary system also.
A similar proof also shows bb* to be a constant. However aa* and bb*
are not Independent, but are related to one another through H and D.
Even though a good physical interpretation cannot be given to the
constants K, L, and D their effect on the orbit under Poisson bracket
can be calculated i.e., the infinitesimal change each produces.
(3.103)
Before calculating the effects of the constants on the orbit, it
is desirable to write the orbit equation in terms of the constants of
the motion and the time. This can be done by defining a complex vector
r = x + ly (3.106)
which can be written In terms of the eigenvectors by inverting Eqs.
(3.9) and (3.10)
r = (u + v*)/2c (3.107)
where
c = [m(w2 + W21/2. (3.108)
Employing the polar forms of both the eigenvectors in Eqs. (3.56) and
(3.57) and of the constants in Eqs. (3.59), (3.65) and (3.66), u and v*
may be expressed in terms of the constants as
u = (R/2)1/2 [(K2 + L2 + D2)1/2 + D]1/2 ei' (3.109)
v* = (R/2)1/2 ((K2 + L2 + D2)1/2 D]1/2 ei' (3.110)
However from Eqs. (3.63) and (3.64)
p = (1/2R)(Y + P) (3.111)
and
o = (1/2)(T P) (3.112)
Hence
u = (R/2)1/2 [(K2 + L2 + D2)1/2 + D]1/2 e1i/2R eI~/2R
(3.113)
and
S (R /2 [(2 2 2)1/2 1/2 1I/2 1T/2
v =(11/2) [(K + L + D ) DI e e
(3.114)
From Eqs. (3.65) and (3.66) It is apparent that
0 = tan' L (3.115)
K
The time dependence is in Y since Y1/2\2 in Eq. (3.61) Is the coordinate
canonical to H and hence from
d/dt (~/2%2) = ({/2%2, H) = 1 (3.116)
it follows that
Y = 2%2(t to) (3.117)
or
S= 2R l(t to) (3.118)
Substituting Eqs. (3.115), (3.117), and (3.118) Into Eqs. (3.113) and
(3.114) and then substituting these Into Eq. (3.107)gives
r = 1 (.E)1/2([(K2 + L2 + D2)1/2+D]1/2e(i/2R)tan1 (L/K) ei'l(tto)
2c 2
+ [(K2 + L2+ D2 1/2D /2(i/2)tan (L/K) ei.2(tto))
= r + r (3.119)
where
r 1 ()/2 [(K2 + L2 + D2) 1/2D1/2 (i/2R)tan1 (L/K) ell(tto)
S2c 2
(3. 120)
and
+ = ( R /2[(K2 + L2 + D2)1/2D]1/2e(i/2)tan' (L/K) eIX2(tt).
+ 2c 2
(3.121)
Equation (3.119) expresses the orbit equation as the sum of two contra
rotating vectors, whose frequencies of rotation are %1 and X2, and
whose amplitudes are functions of the constants of the motion.
The Poisson bracket of the constants K, L, and D with r are:
(K, T") = L I KD
2 [(K2 + L2 + D2)1/2_D] K2 + L2
L L IKD 
2 [(K2 + L2 + D2 1/2+D] R(K2 + L2) (3.122)
S7 K 1LD
L, ) = [ 2 2 K +ILD
2 [(K2 + L2 + D2 D/22 K2 + L
a K ILD
2 2 2 1/2 2 72 .
2 [(K + L + D ) +D] R(K2 + L )
(3.123)
i(, la ia
(D, 1J = r.. r+ (3.124)
where a is a constant defined in Eq. (3.23).
The infinitesimal change induced by the constants is that each
changes both the amplitude and phase of r+ and r_, while preserving
the sum of the squares of their sum and difference. If the orbit were
an ellipse, this last statement would be equivalent to stating that the
sum of the squares of the semi axes of the ellipse is a constant. The
proof of the statement goes as follows:
Before the infinitesimal change the quantity is
Ir+ + r I2 + r+ r12 = 2(r +12 + Ir_12)
(3.125)
and after the change it is
jr+(1 + E) + r_(l + 8) 2 + r( + r ) r.(1 6) 2
2[r 12 + jr 12 + (E + *) Ir 2 + (5 + 8*) r 2] (3.126)
where second order terms have been Ignored. The difference between the
two terms is
2[(E + E*) Ir+12 + (8 + 5*) IrI2]
which is zero In all three cases. For example, if the infinitesimal
change is induced by K, one has for the difference
2[(e + E*)Ir+12 + (5 + 8*)Ir_ 2]
2a 2L jrj2 2L r12
2 [(K2 + L2 + D2)1/2D [(K2 + L2 + D2) /2+D]
2aR L[(K2 + L2 + D2)1/2D L[(K2 + L2 + D2)1/2+Di
8c [(K2 + L2 + D2)1/2D] [(K2 + L2 + D2) /2+D]
0. (3.127)
Gauge Transformations in Uniform Magnetic Fields
in uniform magnetic fields and in the absence of other external
potentials the following fact concerning a charged particle in this
field Is true: The change in momentum In going from one point to
another Is Independent of the path taken between the points. The
proof Is as follows: Let P denote the canonical momentum, A the mechani
cal momentum, B the uniform field, and A the vector potential such that
B Vx A (3.128)
and since B is uniform
A B x r (3.129)
2
Since
d e V x (3.130)
dt c
and
v dr (3.131)
dt
one has 
rI r2
f d f f ( x B) dt
c dt
F2
= .f dr x
Trl
T2
2e f dA (3.132)
c 
rI
Hence
S(3.133)
c
and
aP  (3.134)
c
Equation (3.134) states that the difference in moment between any two
points is proportional to the difference in the vector potential evalu
ated at these points. From this it Is Inferred that the choice of a
certain gauge Is equivalent to picking the zero of momentum.
One can also see why a translation of coordinates must be accompanied
by a gauge transformation. When the translation occurs the zero of momen
tum changes and hence this change must be subtracted from A to give the
the same zero of momentum. However, the change In A is simply a gauge
transformation.
In order to show that a gauge transformation can induce a translation,
first consider the case of a free particle In a uniform field with the
following gauge A B (y, x, 0) (3.135)
2
In units where m and w eBo/2mc are one, the constants of the motion are
and
Using the fact that
P y a
P + x B .
p +x1
P v+ ,
c
(3.136)
(3.137)
(3.138)
the location of the center of the orbit is found to be
(Xc, yc) 2 (2, 2)
(3.139)
However with the gauge
B
A = 2 ( (X + y), + x, 0),
2
(3.140)
even though the two constants are the same, one finds the center of the
orbit to be
(xc yc) ( __2 a3 X )
2 2
(3.141)
Hence the center of the orbit has been translated 1t/2 units in the x
direction and X/2 units in y.
The center of the cyclotron orbit may also be rotated, however
this is not accomplished by a simple gauge transformation. Instead
one must rotate the vector A and also rotate the coordinate system.
For the vector potential defined in Eq. (3.135) the center of the orbit
is given in Eq. (3.139). A rotation of the vector field defined by 1
gives
cos 6 sin 0 0 y
B,
AR = sin 0 cos O 0 x
2
0 0 1 0
B
= ( y cos 6 + x sin 0, y sin 0 + x cos 6, 0). (3.142)
2
The rotation of the coordinates defined by
x U x' cos 6 y' sin 6 (3.143)
and
y x' sin 0 + y' cos 6 (3.144)
transforms AR into A' where
B ,
A' = (y', x', 0) (3.145)
Hence, the center of the orbit is
(xc ) ( , ) (3.146)
and
B B cos 0 + a sin 0
2 2 2
(3.147)
a B' a'
S  sin 0 + cos 0, (3.148)
which shows that the center of the orbit has been rotated with respect
to the original coordinate system, while the diameter of the orbit has
remained unchanged.
The vector field can also be rotated without rotating the coordi
nates and still preserve the curl if at the same time the vector is
dilated by a factor of sec 6. The new vector potential is
Bo
A = sec 6 (y cos e + x sin 6, y sin 6 + x cos 0, 0).
2 (3.149)
The transformation leaves the center of the orbit invariant and in fact,
Is equivalent to a gauge transformation where the vector AD
B
AD = tan O(x, y, 0), (3.150)
whose curl is zero, has been added to the original vector potential A.
CHAPTER 4
THE KEPLER PROBLEM
The Two Dimensional Case in Polar Coordinates
The Hamiltonian for the two dimensional Kepler problem in polar
coordinates Is
H = (1/2)(P2 + p2/r2) 1/r (4.1)
where the mass and force constant have been chosen equal to unity.
The corresponding HamiltonJacobi equation is
(S )2 + (,)2 2(W+ ) = 0. (4.2)
or r2 ao r
where W is the total energy and
S(r, 4) = Sr(r) + S(~Q) (4.3)
Hence
p S (4.4)
Pr 7r
and
p Bs
P ; (4.5)
Since the partial differential equation (4.2) has no P dependence
bS/ba Pp Is a constant and is, in fact, the angular momentum of the
system. Solving Eq. (4.2) for 6S/or = Pr one finds
2_2 1/2
= Pr = [2W + 2 (46)
dr r r2
The action variables Jr and JO are given by the integrals
and
(4.7)
(4.8)
Jr P rdr
J, = 5 PVd$ .
Equation (4.8) can be Integrated immediately since the integral over
one cycle is simply an Integral from 0 to 2A, hence
J = 2P .
(4.9)
The Integral in Eq. (4.7) can
as described in Born [35] and
be evaluated by a contour Integration
the result is
2*
Jr = 1/2 J
[2W]
(4.10)
The angle variables wr and wo
jugate to Jr and Jo are given
which are canonical coordinates and con
by
6S Pr
Wr = j" dr
r r
and
sS Pr 3P$
r = f dr + JZ d
Upon performing the integrations
wr = (u E sin u)/2xr
and
w = [u 6 sin u tanl (y/x) + <]/2n,
(4.11)
(4.12)
(4.13)
(4.14)
where the elliptic anomaly u is defined through
r = a(l E cos u),
a being the semimajor axis and E being the eccentricity of
a = W/2
2 2
By choosing o to be the azimuth of the perihelion,
x r cos ( 0 0)
y = r sin ( () ,
(4.15)
the ellipse,
(4.16)
(4.17)
and Eq. (4.14) reduces to
wo = (u E sin u + 0 )/2x (4.20)
Solving Eq. (4.10) for the energy in terms of Jr and J
22
w = 2(4.21)
(Jr + 2
and hence the frequencies vr and v are degenerate i.e.,
V v & = (4.22)
T* a ab
Thus the quantities aa bb ab, and ab are constants of the motion
where J 1/2 2riw
a = () e r (4.23)
and b = ()1 e2nw (4.24)
2n
These quantities satisfy the following Polsson bracket relations
(a, a*) (b, b*) i (4.25)
where one uses the fact that the J's and w's form a set of canonical
coordinates and moment.
Taking the real and imaginary parts of these constants and the sum
and difference of the real ones gives the following four constants of
the motion,
H = aa* + bb* (4.26)
D = aa* bb* (4.27)
K = ab* + a*b (4.28)
L = i(ab* + a*b). (4.29)
Of these constants H commutes with the other three and is a function of
the Hamiltonian, namely, it is inversely proportional to the square
root of the energy. The three constants K, D, and L satisfy the com
mutation relations of the generators of the group SU2 or R3.
(K, D) 2L (4.30)
(D, L) 2K (4.31)
(L, K) = 20 (4.32)
Although they have no simple Interpretation as constants of the motion,
one can make a mapping of these constants which preserves the commuta
tion rules. The new constants can be interpreted as the angular
momentum and the two components of the Runge vector [Il].
is
K K' = K(2 + H+D)1/2
HD
L +L' = L(2 + H+D)1/2
HD
D D' = H D
These constants satisfy the rules
(K', D')
(D', L')
(L', K'}
= 2L'
= 2K'
= 2D'
Otherwise they are a direct consequence of
K, L, D, so that this mapping is valid for
must be made of the identity
H2 = K2 + L2 + D2
the commutation rules for
any SU2 group where use
(4.39)
K' is essentially the x component of the vector P where
P = R 2E = (e//2E) Cos 0
sin 0
0
(4.40)
and R is the Runge vector, while L' is the y component and D' is the
angular momentum.
The mapping
(4.33)
(4.34)
(4.35)
(4.36)
(4.37)
(4.38)
The Kepler Problem in Three Dimensions
Using the same units as in the previous case, the Hamiltonian for
the Kepler problem in three dimensions is
2 2
1 2 i' P
H = (Pr +  )  (4.41)
2 r2 r2 sin2 0 r
Setting up the HamiltonJacobi equation and separating gives
P = S= a = constant (4.42)
S a 1/2
0 sin2 )
and
S 2 1a /2
Pr (2W + (4.44)
ar r r2
where a9 and aO are separation constants, aO being the magnitude of
the Z component of the angular momentum and a0 the magnitude of the
total angular momentum. Since the potential is central the orbit is
planar and the normal to the plane is necessarily parallel to the
angular momentum. Denoting by B the dihedral angle which the orbital
plane makes with the xy plane one can deduce that
cos B  (4.45)
a0
Solving for the action variables gives
JO = 2r0( = 2 nP
(4.46)
88
Je 2(C(a ca) (4.47)
Jr JO JO (4.48)
and from Eq. (4.48) It follows that
2
(Jr + Je + j)2 (4.49)
In the plane of the orbit rectangular coordinates 4 and v can be
introduced, where the 4 axis is the major axis of the ellipse and the
origin coincides with the center of force and the origin of the polar
system. In this system
4 = r cos (T To) (4.50)
v = r sin (Y To), (4.51)
where Y is an angle measured from the 4 axis to the particle and To
is the angle to the perihelion.
The angle variables are
oP
wr Pr dr = (u e sin u)/2n (4.52)
6Jr
we J dr + j dJ
PJO r r
= (u e sin u tan1 V + sln1 cos O)/2, (4.53)
p. sin B
w= dr + f dO + f dP
= (u E sin u tan1 V + sin1 Cs sin' cot + +)/2
i. sin B tan B
(4.54)
where u is defined as in Eq. (4.15) and E as in Eq. (4.17). Using Eqs.
(4.50) and (4.51) It follows that
tan1 v = T o "(4.55)
It follows from Fig. 4.1 that sin1 cot 6 is the projection of the
tan B
angle T onto the plane in which 0 is measured I.e., the xy plane, which
shall be denoted . It also follows that sin1 cos 0 Is the angle T,
sin B
where the line of nodes has been chosen to coincide with the major axis of
the ellipse. Hence, the line of nodes makes an angle (0 ) with the
xaxis. Rewriting Eqs. (4.53) and (4.54) results in the following set
of equations,
we = (u e sin u + To)/2i (4.56)
we = (u e sin u + \o + ( ^ )/2n (4.57)
By defining
a = ()/2 e2iwr (4.58)
2ir
b (Je)1/2 e2niv (4.59)
2n
J= 1/2 27iw
e (4.60)
2c
Fig. 4.1 Three dimensional coordinate system for
the Kepler problem.
one finds the following nine constants of the motion: aa*, bb*, cc*,
ab a b, bc b c, ac and a c. These constants are the generators
of the group SU3. The constant angles B, To and 4  and also the sum
"+ To can be expressed In terms of the constants,
sin Yo M i(ab* a*b) (4.61)
2(aa*bb*)1/2
ab* + a*b
cos Yo + 1/2 (4.62)
2(aa bb*)
sin( ') = i(bc* b*c) (4.63)
2(bb*cc*)1/2
cos(4  (bc + b c) (4.64)
2(bb*cc*)1/2
sin( + ) = i(ac* ac) (4.65)
2(aa*cc*)1/2
cos(  + ) ac* + a*c(4.66)
2(aa*cc*)1/2 66
and
cos B = cc (4.67)
bb* + cc*
Generally, the rotation group In four dimensions R4, is taken to
be the symmetry group of the Kepler problem. Its generators are the
three components of the angular momentum vector L and the three compo
nents of the normalized Runge vector P which Is defined in terms of
in Eq. (4.40). The rectangular components of P are
cos Yo cos (0
P E cos Yo sin (0
sin B sin To
) cos B sin yo sin (0 )
) cos B sin T' cos ( )
(4.68)
For L one has
sin B sin ( )
sin B cos (P )
cos B
Using Eqs. (4.61) through (4.67) one can express
in terms of the generators of SU3 as follows:
S(aa* + 2bb* + 2cc)1/2 [(ab* + a*b)(bc*
4bb*(cc*)/2(bb* + cc*)
the generators of R4
+ b*c)(bb* + cc*)
+ (ab* a*b)(bc* b*c) cc*]
(4.70)
p l(aa* + 2bb* + 2cc*)1/2
S4bb*(cc*)I/2(bb* + cc*)
[(ab* + a*b)(bc* b*c)(bb* + cc*)
+ (ab* a*b)(bc* + b*c) cc*]
i(aa* + 2bb* + 2cc*)1/2
z 2(bb* + cc*)
Lx = i(bc* b*c)(bb* + 2cc*)
2(cc*)1/2
L (bc* + b*c)(bb* + 2cc)l1
y 2(cc*)1/2
(bb* + 2cc*) /2(ab* a*b)
1/2
JO +
2x
(4.69)
(4.71)
(4.72)
(4.73)
(4.74)
and Lz = cc* (4.75)
These six quantities satisfy the relations
(Li, P)  Eij Pk (4.76)
(Li, Lj} Eljk Lk (4.77)
(Pi' Pj) ijk Lk (4.78)
which define R4, once again directly from the commutation rules for SU3.
The process can be inverted and the generators of SU3 can be ex
pressed in terms of the components of P and L. These are
aa* (L2 + L2)1/2 L (4.79)
S x2 y
bb* = (L2 + L2)1/2 (4.80)
cc* L (4.81)
ab* (L2 + L2)/2 L1/22 + L2/2(cos o I sin )
(4.82)
bc* (Lz(L2 + L2)1/21/2 (cos () 1 sin (0 )) (4.83)
z x y
ac* [ (L2 + L2)l/2Lz]/2Lz /2[cos( ~ + Y ) sin(4 + )]
=2W x y z 1
(4.84)
where L
4 ^= tan" x (4.85)
Ly
and
(4.86)
 t 1 PzLy(L + L2 + L2)1/2
SP(L2 + L2) + PzLxL
2 2 1/2
S tan1 (L + L)
B = tan'VX^^ 
(4.87)
The Two Dimensional Kepler Problem in Parabolic Coordinates
in parabolic coordinates, defined by
x =L v (4
y = (L2 v2)/2 (4
and where m = e = 1, the Kepler Hamiltonian is
P2 + p2 4
H = 0 (4
2(42 + 2)
From the HamiltonJacobi equation
P[ = (2wL2 + 2a1)1/2 (4
PV = (2Wv2 + 2a2)1/2 (4
where
al + a2 = 2. (4
The action variables defined by
Jl = Pidq (4
give 2nal
= (4
GL 2
.88)
.89)
.90)
.91)
.92)
.93)
.94)
.95)

Full Text 
xml version 1.0 encoding UTF8
REPORT xmlns http:www.fcla.edudlsmddaitss xmlns:xsi http:www.w3.org2001XMLSchemainstance xsi:schemaLocation http:www.fcla.edudlsmddaitssdaitssReport.xsd
INGEST IEID EVBD08PTA_2YMF18 INGEST_TIME 20170719T20:26:26Z PACKAGE UF00097926_00001
AGREEMENT_INFO ACCOUNT UF PROJECT UFDC
FILES
PAGE 1
A STUDY OF ACCIDENTAL DEGENERACY IN HAMILTONIAN MECHANICS By VICTOR AUGUST DULOCK, JR. A DISSERTATION PRESENTED TO THE GRADUATE COUNQL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY UNIVERSITY OF FLORIDA June, 1964
PAGE 2
In memory of John F. Kennedy
PAGE 3
ACKNOWLEDGMENTS The author wishes to express his sincere appreciation to his supervisory cocmittee and in particular to Professor H.V. Mcintosh, who directed this research. He would also like to acknowledge the counsel and guidance given him by Dr. B.S. Thomas throughout his graduate program. The author also wishes to thank the Department of Physics and the Graduate School of the University of Florida for their financial support during his graduate career. Finally, the author gratefully acknowledges the help and encouragement given him by his wife, Maryann.
PAGE 4
TABLE OF CONTENTS Page ACKNOWLEDGMENTS ill LIST OF FIGURES vl INTRODUCTION 1 Purpose I Historical Sketch 2 CHAPTER 1. MATHEMATICAL FOUNDATIONS 7 Introduction 7 Lie Algebra Defined by Poisson Brackets 8 2. THE HARMONIC OSCILLATOR 20 Introduction 20 Tlie Plane Isotropic Oscillator 20 Tine Hopf Mapping 2k Canonical Coordinates 29 The Anisotropic Harmonic Oscillator in Two Dimensions 31 Commensurabi 1 i ty and Constants of the Motion .... 3^ The Isotropic Oscillator in n Dimensions 36 The Hopf Mapping 38 Other Canonical Coordinates k2 The Anisotropic Oscillator in n Dimensions ^6 3. THE HARMONIC OSCILLATOR IN A UNIFORM MAGNETIC FIELD . . 50 Introduction 50 Classical Zeeman Effect for the Harmonic Oscillator 52 Limiting Cases 55 Canonical Coordinates 59 Rotating Coordinates 61 Discussion of the Constants of the Motion 70 Gauge Transformations in Uniform Magnetic Fields . . 77
PAGE 5
TABLE OF CONTENTSContinued Page CHAPTER 4. THE KEPLER PROBLEM 82 The Two Dimensional Case in Polar Coordinates .... 82 The Kepler Problem In Three Dimensions 87 The Two Dimensional Kepler Problem in Parabolic Coordinates 94 5. SUMMARY 98 LIST OF REFERENCES IO5 BIOGRAPHY 108
PAGE 6
LIST OF FIGURES Figure Page 2.1 Gnomon I c Projection 26 2.2 Inverse Stereographic Projection 27 3.1 Particle orbit In a strong uniform magnetic field and harmonic oscillator potential 6k 3.2 Particle orbit in a weak uniform magnetic field and harmonic oscillator potential 65 3.3 Particle orbit in a fixed uniform magnetic field and harmonic oscillator potential with a large Initial tangential velocity 66 3.^ Particle orbit in a fixed uniform magnetic field and harmonic oscillator potential with a small Initial tangential velocity 67 3.5 Particle orbit shown In Fig. 31 as viewed from a rotating coordinate system 7' 3.6 Particle orbit shown in Fig. 32 as viewed from a rotating coordinate system 72 4.1 Three dimensional coordinate system for the Kepler problem 90
PAGE 7
INTRODUCTION Purpose One ordinarily defines the accidental degeneracy of a system as that which does not follow from an analysis of the obvious geometrical symmetry of the system in configuration space. Frequently, nevertheless it can be shown that there is a higher order symmetry group in the phase space of the system, formed with the help of certain hidden symmetries, which is adequate to account for this accidental degeneracy. A well known example of a system possessing this property, and one of the first to be investigated historically, both classically and quantum mechanically, is the nonrelat i vist ic Kepler problem, whose quantum mechanical analogue is the problem of a oneelectron atom. In the quantum mechanical case, the invariance of the Hamiltonian with respect to the three dimensional rotation group accounts for the degeneracy in the "m" quantum number, but not for that in the "1" quantum number. The purpose of this work is to review the more important results obtained on the accidental degeneracy of harmonic oscillators and the Kepler problem and to find some general principles with which the symmetry groups of Hamiltonian systems can be found. In Chapter 1 a basic mathematical foundation is laid. This foundation is based on the fact that all quadratic functions over the phase space of a system form a Lie algebra. On this basis a method is developed which enables 1
PAGE 8
one to find, rather easily, the constants of the motion for systems with quadratic Hami 1 tonians or for those which can be made quadratic through a canonical transformation. This method is then applied to various systems in the succeeding two chapters. The harmonic oscillator is treated in Chapter 2, both in the two and in the n dimensional case. Chapter 3 deals with a general Hamiltonian which may be interpreted as that of a charged mass point in a plane harmonic oscillator potential and uniform magnetic field. Finally, in Chapter 4 the Kepler problem is discussed and found to have an SU3 symmetry besides the Rl^ generally associated with it. All calculations are carried out classically. The results obtained seem to be valid in the quantum mechanical case except as noted. The transition to the quantum mechanical version is carried out in the usual manner by replacing Poisson bracket relations with commutator braci
PAGE 9
aphelion and a length equal to the eccentricity of the classical orbit and he also gives the quantum mechanical operator definition. The comutatlon rules between the components of this vector and the angular momentum are also stated explicitly. As a consequence of the the commutation rules he shows that the system is degenerate in "1" also. However, L. Hulthen [3] observed the fact that these commutation rules are the same as those of the generators of the four dimensional rotation group. 0. Klein was given credit for first recognizing this fact. V. Fock [k] , in 1935. using a stereographic projection in momentum space, solved the Schrodinger equation in momentum space and showed that the solutions are spherical functions of the four dimensional sphere and hence that the symmetry group is R^. Later V. Bargmann [5], in a discussion of Fock's results, showed that Fock's group was generated by the components of the two vector constants R and L, the angular momentum. It was about this time that 0. Laporte [6] obtained a differential operator identical to the Runge vector and presented a method for obtaining all the eigenfunct ions of the hydrogen atom by differentiation, similar to the method used to find the eigenfunct ions of L by using the raising and lowering operators L.,. and L_. He made extensive use of the stereographic parameters, developed by him and G.Y. Rainich [7] about the same time, in momentum space. The wave functions were found to be proportional to hyperspher leal surface harmonics. In 1939 and 1940 E.L. Hill and his student, J.M. Jauch [8, 9, 10] published two articles and a Ph.D. thesis dealing with accidental degeneracy. These authors considered both the two and three dimensional
PAGE 10
Kepler problem and the harmonic oscillator in two and three dimensions. The Kepler problem was found to have R, and R^ symmetry groups in two and three dimensions respectively, while the harmonic oscillators had the symmetry of SU2 and SUt, the groups being generated by the constants of the motion in all cases. The transformations induced by the generators were found to have an analogy in classical mechanics, namely that they represented transformations of one orbit of phase space into another of the same energy. They also showed the existence of a correspondence between transformations in classical mechanics and quantum mechanics. In 19^9, A.W. Saenz [11] wrote a Ph.D. thesis under Laporte, in which the symmetry groups explaining the accidental degeneracy of various problems were found by reducing the problems to force free motion on hyperspherical surfaces of various dimensions. In this manner he considered the Kepler problem, harmonic oscillator and rigid rotor. More recently, the n dimensional isotropic oscillator has been discussed by G.A. Baker [12]. He found the synmetry group to be SU^ and demonstrated that this group could account for all the degeneracy found. This had been demonstrated earlier by Y.N. Demkov [I3] from a slightly different point of view. Whereas Baker attacked the problem by finding the most general group which left the Hamiltonian invariant, Demkov used the generators of infinitesimal transformations which commuted with the Hamiltonian and then showed that this defined an SU^^ group. In a later article Demkov [l4] demonstrated the equivalence of his and Baker's works.
PAGE 11
The connection between accidental degeneracy and hidden symmetry was confirmed by A. P. Alliluev [15] in 1958. He considered two examples, the two dimensional oscillator and the Kepler problem in n dimensions, which reduced to Fock's result in the case where n = 3. In the general case, the hidden symmetry of the Kepler problem in n dimensions was found to be that of the n + 1 dimensional hypersphere. The one dimensional Kepler problem has been treated by R. Loudon [16] and was found to be doubly degenerate in all levels except the ground state, an apparent violation of the theorem that the energy levels of one dimensional systems are nondegenerate. Loudon discusses his results in the light of this theorem and shows that it holds only when the potential is free of singularities. In 1959. H.V. Mcintosh [17] wrote a qualitative review of the work up to that time and stated the results. The most recent article to appear is one by Demkov [18] in which the concept of excessive symmetry groups is introduced. He gives a prescription for finding the minimal symmetry group and as an example applies it to the anisotropic two dinnensional oscillator with incommensurable frequencies. In addition to the works cited above dealing with thestudy of accidental degeneracy itself, there has been a considerable amount of work in using the symmetry groups for various problems and applying them to the many body problem, particularly to the shell model of the nucleus. One of the first of these was by J. P. Elliott [19, 20] who undertook to study collective motion in the shell model utilizing a coupling scheme associated with the degeneracy of the three dimensional
PAGE 12
harmonic oscillator. The direct result was a classification of the states according to the group SUo. Following this treatment was a series of papers by V. Bargmann and M. Moshinsky [21, 22] on the group theory of harmonic oscillators, in which they developed a classification scheme for states of n particles in a common harmonic oscillator potential. They were primarily concerned with finding a scheme which would exhibit explicitly the collective nature of the states. Similar work was also done by S. Goshen and H.J. Lipkin [23, m] with the exception that the particles were in a one dimensional f>otentIal well. This system could also be described from a collective viewpoint and these states could be classified using the group of transformations. In recent times, L.C Biedenharn has extensively studied the Kepler problem and its symmetries, both in the nonrelat i vist ic and relativistic cases. In the nonrelat i vist ic [25] case he has obtained relations for the representations of R. exploiting the fact that Rk is locally isomorphic to Ro x R^. The nonrelat ivist ic problem has also been considered by Moshinsky [26], where use has been made of the accidental degeneracy in the coulomb problem to obtain correlated wave functions for a system of n particles in this potential. In the relativistic case, Biedenharn [27] and Biedenharn and Swamy [28] have obtained operators which are analogous to the angular momentum and Runge vector in the nonrelat ivist ic problem.
PAGE 13
CHAPTER 1 MATHEMATICAL FOUNDATIONS Introduction In classical and quantum mechanics "quadratic" Hamiltonians play an exceptional role, a fact which is traceable in the end to the fact that the Poisson bracket operation satisfies axioms whose expression is tantamount to saying that the Poisson bracket is an alternating bilinear functional which is a derivative in each of its arguments. If H^"' designates the vector space of homogeneous polynomials of degree n, whose basis may be chosen to be composed of homogeneous monomials, the Poisson bracket operation maps the cartesian product ^(m) X h(") into H^"^""^) , as will be shown later. The Poisson bracket operation is a bilinear functional, which means that it defines a linear transformation when either argument is held fixed, and so it always has a matrix representation. For the case m = 2, one has a mapping from h(") into itself, representable by a square matrix. In this latter case, the extensive lore of matrix theory is available to help discuss the transformations. In particular, the mapping may be described in terms of its eigenvalues and eigenvectors. Thus the mapping Tq(p) = (q. p}. (1.1) for a fixed q Â€ h' ' , and argument p e h'"^, can be expected to have 7
PAGE 14
specialized properties which become all the more significant when n = 2. In this case the elements of H^ ' form the basis of a Lie algebra under the Poisson bracket operation. In any event, the general mapping T^ of Eq. (I.I) yields representations of this algebra as it acts on polynomials, homogeneous of various degrees as shall be shown later for the case n = 2. Thus the allusion to the exceptional role of quadratic Hamiltonians refers to the expectation that the algebraic structure of the Lie algebra will have its reflection in the dynamical properties of the system which it represents. Lie Algebra Defined by Poisson Brackets The set of 2n linearly independent coordinates and momenta of an arbitrary physical system is said to form its phase space 0. These variables may be paired in such a way that they can be indexed by fj, i = n, ...I, 1, ...n, where the positive indices indicate coordinates and the negative ones the corresponding conjugate momenta. The Poisson bracket of any two functions of these variables is defined by ''Â•"''Zl^srrs^lTr<'^' i = l It follows from Eq. (1.2) that the Poisson bracket operation satisfies the axioms of a Lie algebra [29, 3?]: (f.g) = [g.f) (alternating rule), (13)
PAGE 15
{f, ag + Bh} = a[f , g} + B(f, h} (linear), {] .k) if. ig. h)) + {g. (h, f}) + {h, {f, g}) = (Jacobi identity). (15) The Poisson bracket also acts like a derivative in that it satisfies the relation {f. gh) = g(f, h} + [f, g)h . (1.6) Although all second degree polynomials form a Lie algebra, the linear monomial functions of the coordinates in phase space itself also form a linear vector space 0, where the basis may be chosen to be the 2n fj. Denoting the Poisson bracket operation { , } as the symplectic inner product in i> , it is noted that this product does not satisfy the usual axioms for an inner product since it is antisymmetric under interchange of the arguments. This difficulty may be remedied in the following manner. Using the definition of the Poisson bracket it is found that 1 i = j > (f.. f_.) =0 i ^ j I i = j < . (1.7) From this it is observed that if i > 0, the dual vector f. to f. may be taken to be f_j. Similarly for i < the dual vector is taken to be f .. Thus the dual basis to r= (f^. Â•Â•Â• f,. f.,. ... f.n) (l8a)
PAGE 16
10 is then r*= (f.Â„. .. f... fp ... V(1.8b) with these definitions the inner product may be defined as (f. g) = {f. g}. (1.9) where f and g are linear combinations of the fj and = )^atfi . (1. 10) f i It is evident then that (f,. fj) = (f,. fj} =6,j (1.11) and (f. f) = (f. f) =)^i^il^ Â• <^'2) The space also has the peculiar property of being selfdual in a special way, namely the mapping Tj: is a linear mapping and is in a one to one correspondence to the vectors in the dual space $*". Considering tensor products of $, it is observed that (f, g) also defines a mapping on these products. If one defines h(") = Â® = (^^, (113) and if f Â€ h("') and g e h(") . then (f, g) 6 h("^""^). In particular, when m= 2 i.e., f e Q s h' ' , a mapping is obtained which is representable by a square matrix.
PAGE 17
n Such a mapping may be expected to have eigenfunct ions and eigenvalues. Let h be a fixed element of Q (a quadratic Hamiltonian) and let g; and \be the eigenfunct ions and eigenvalues such that (h. g;) = Th(g;) = ^jg;. (1.14) where the gj are taken to be linear I.e. g; Â€ H^ '. If h is representable as a normal matrix, the gj form a basis for the space h' ', and are simply linear combinations of the f. In particular, since the gj Â€ <& s h' ', they form a basis for all the tensor spaces as well where all possible products of the gj are taken i.e., ((gi)) = H<') ((9,9^)) = H^2) ((gj.... 9j) = H<"). (1.15) On account of the derivative rule Eq. (1.6) the following relation is satisfied {h, gjgj) = g;(h, g j ) + {h, gj) gj = (K; + Xj) gjgj, (1.16) and hence the monomials in gj are eigenfunct ions of h in every H^"' and form a basis even in the case of degeneracy, for all possible eigenfunct ions. In cases where the operator Tf is normal, a complete set of orthogonal eigenvectors exists which will be denoted g.. Hence, if the
PAGE 18
12 eigenfunctions are indexed such that gj = g_j then the eigenvectors of any normal operator satisfy the rule (gj. 9j} = &j,.j. (1.17) The eigenvalues \ corresponding to the eigenvectors gj of the normal operator h occur in negative pairs. This follows from a consideration of the Jacobi identity lh,{g_., g.)} + {g_j, {g., h}) ^ {g..{h. g_.)) = 0. thus (h, 1} + (\.) [g_.. g.} + \_. {g., g_.) = 0. and (^i + ^0 {gj. gi) = 0, (1.18) and hence, \. = \ . since {q., q .} / 0. Since h e H^ ' , h may be written as n = ) C. .g.g. , h = ) C,.g.g.. (1.19) i,j=n Â• = J where the g; are eigenfunctions of h belonging to eigenvalues \j and are indexed so that K^ = \_i . (1.20) Assuming that the eigenvalues of h are distinct and using the identity {h, h) = one has, with the aid of Eq. (119)
PAGE 19
13 =\ C..[h. g.g.} Â«.J i2y \ C..(\. + \.) g.g.. (1.21) '.J Hence each coefficient must vanish and either Cj j = or {\+ >.:) = 0. However X.? = \j only if j = i. Therefore n h=^ Cj^.j g;g.;. (1.22) i=i To further evaluate these coefficients recall that {h, g^) = \^g^. (1.23) thus A H9k =2. ^!.' ^^'^'' ^'<^ iÂ«l n Using the fact that the coefficients C.. are symmetric in the indices the coefficient of g^ on the right side is \=Ck,k^^k' ^k^' <'25)
PAGE 20
]k and hence C. .^^ ('26) The final conclusion is that h = Iji^''''<'"' l=i I I Consider now the case where h may have degenerate eigenvalues. Using the general form for h, Eq. (1.19), and calculating {h, h) as before Cj j = unless X. = \. . However in the present instance it is possible to have X. = \ , k ?^ m and hence Xi = X Â„. Hence "^ km k m g.g . (1.28) , m ^k^m * n
PAGE 21
15 + ) C, [g, {g , g.) + g [g. , g.}] . / ^ k,m ^k ^m ^j ^m^k* ^j k,nT=l ^k=Nn k it m Consider the first term in the second sum. If m=j then k?^ implies k.^} and hence C. _. must be zero. Therefore, as before h = i = I II y r^ rgigi(130) It should also be noted that the transformation to eigenvector coordinates for h is a canonical transformation since the eigenvectors satisfy (gj, gj) = Sj _., which is a necessary and sufficient condition for g. and g . to be canonical variables. Consider now the canonical form of a nonnormal quadratic operator h, of the Jordan canonical form. Again n n "=11 C,.g,g.. (,.3.) i = n jÂ«=n where the gj are the principal vectors of h and indexed such that the same orthogonality conditions hold as before. Using the property that (h, g.) = g;^,. (1.32) it is found that
PAGE 22
16 i=n j=l i i J = ) ) C [g {g., g } + g.{g., g, )] l=n j=n i i J 'Ck+l,k3k+l '<>0 .Ck.k+l Sk+l Â•<" . (1.33) Hence all C: ; = 5. _. , and n h=)^ gig.i^, . (1.3^) i=l The sympletic norm also produces certain symmetries and antisymmetries. For example, the dual mapping R(fj), which maps the coordinates into momenta and the momenta into the negatives of the coordinates , i > (1.35) R(f .) = f. preserves the inner product; i.e. (R(fi). R(fj)} = {f;. fj) . (136) On the other hand the mapping P(f:) which simply exchanges coordinates and momenta
PAGE 23
17 P(fi) = fi. (1.37) or the time reversal mapping T(fj) Kfj) = f, i > (138) T(f.i) = f.i reverses the sign of the symplectic norm {P(fj). P(fj)) = {f,. fj] (T(f.). T(f.)) = (f.. f.} . (1.39) These operators may be expressed in terms of the eigenvectors, g. of the normal operator h, previously defined. Recall also that the gr form a set of canonical coordinates and momenta. Using the convention that positive indices represent coordinates and negative indices the corresponding canonical momenta, these operators have the following canonical forms: n i=i
PAGE 24
18 These operators have the following effects on eigenvectors of h: r^k Â•< > 1^(9^) " (T. g^) = j (1.43) Tp(g,^) = {P.g,^} = g_,^ all k, (1.^) f g.k Â•< > L9k ^ " Â° Â• The composite of any two of the operations on an eigenfunct ion is equivalent to operating with the third operator: Mk ^^ 0] VV9,)).9., = yg,) (I.V) It is also interesting to note that the three operators satisfy the following Poisson bracket relations among themselves: (P. T) = 2R (1.^9) (R. T) = 2P (1.50) {p. R) = 2T (1.51) The composite of the normal operator h with each of these operators on the eigenvectors is also worth noting. It is relatively easy to show
PAGE 25
19 that and hence the matrices representing h and T commute. However one finds that P and R anticommute with h, i.e. VWVVV> ('Â•") and T (T (g )) = T (T (g )). (1.54) R h k h R k The existence of either one of these last equations is sufficient to insure that the eigenvalues appear in negative pairs [30]. In the next two chapters the method described above will be applied to harmonic oscillator Hamiltonians in two and n dimensions and to variations on the oscillator Hamiltonian.
PAGE 26
CHAPTER 2 THE HARMONIC OSCILLATOR Introduct ion The accidental degeneracy and symmetries of the harmonic oscillators has been discussed by several authors. A particularly thorough account has been given in the thesis by J.M. Jauch [9], in which he considers the isotropic oscillator in two, three and n dimensions, and finds the symmetry group to be SU^, i.e., the Hamiltonian is invariant under transformations induced by the generators of the group. A short discussion of the anisotropic oscillator in two dimensions with commensurable frequencies is also given. The main results of this work are also contained in the notes of E.L. Hill [33]More recent discussions of the n dimensional isotropic oscillator have been given by G.A. Baker [12] and Y.N. Demkov [13, 14. 18]. The plane Isotropic Oscillator In units such that the mass and spring constant are unity, the Hamiltonian of the plane isotropic harmonic oscillator is H = (P? + Pj + x^ + y2)/2 . (2.1) X y With respect to the basis (x,y,P ,P ), the matrix representation of H 20
PAGE 27
2] as an operator under the Poisson bracket relation is.
PAGE 28
22 {H. a] = (H. P^ ix} = i(P^ ix) (2.6) while multiplying the matrix in Eq. (2.2) into the column matrix representing a also gives i times the same column matrix. The quantum mechanical analogues of these four quantities are the well known raising and lowering operators for the harmonic oscillator. With these four quantities ten linearly independent quadratic monomials can be constructed which are eigenfunct ions of H and have as eigenvalues 0, + 2i or 2i. Only four of these are of interest here, namely those with eigenvalue zero. They are the four linearly independent constants of the motion (since they commute with H under Poisson bracket) listed be 1 ow : aa* = (p2 + x2)/2 (2.7) bb*= (p2 + y2)/2 (2.8) ab"'= [(P^Py + xy) + i(yP^ xP^) ]/2 (2.9) a*b = [(P^Py + xy) i(yP^ xPy) ]/2 . (2.10) For the purpose of physical interpretation it is more convenient to deal with the real and imaginary parts of these constants separately, as well as to separate the Hamiltonian from them. Accordingly the following four quantities are introduced: H = aa'^ + bb* = (P^ + p2 + x^ + y2)/2 (2.11) X y D = aa''^ bb'^ = (p2 + x^ pj y2)/2 (2.12) X y
PAGE 29
23 L = i(a*b ab*) = yP^ xP, K = a""b + ab'^ = xy + P^P y (2.13) (2.14) These four constants have the following simple interpretations. a) H, the Hamiltonian, is the total energy of the system. b) D is the energy difference between the two coordinates. c) L is the angular momentum of the system and the generator of rotations in the xy plane. d) K is known as the correlation and is a peculiar feature of the harmonic oscillator. As a generator of an infinitesimal contact transformation, it generates an infinitesimal change in the eccentricity of the orbital ellipse while preserving the orientation of the semiaxes, and preserving the sum of the squares of their lengths. The energy of a harmonic oscillator depends only on the sum of the squares of the semiaxes of its orbital ellipse, which remains constant under such a transformat ion. The set of functions {K, L, D} is closed under the Poisson bracket operation. Explicitly, their Poisson bracket table is (2.15)
PAGE 30
2k generators of the three dimensional rotation group, or the three components of the angular momentum in three dimensions. Since the configuration space of the two dimensional isotropic oscillator is apparently only two dimensional and has only rotations in the xy plane as an obvious symmetry, the occurrence of the three dimensional rotation group is rather anomalous. By performing a series of geometrical transformations on a, a*, b and b' one can explicitly demonstrate the spherical symmetry of the system. This series of transformations is sometimes called the Hopf mapping. The Hopf Mapping [31] Perhaps the best analytic representation of the Hopf mapping is obtained by introducing polar coordinates (\, t, a, p) for the complex variables a and b. Explicitly, these are a = (\e'P cos t)//2 (2.16) b = (\e'<' sin i)Ar2, (2.17) where \ does not denote the eigenvalue of H as previously. The first step in the Hopf mapping is to form the ratio to defined as the quotient u) = X^ a/b . (2.18) The formation of this ratio may be regarded as a gnomonic projection of the four dimensional space, regarded as having two complex dimensions,
PAGE 31
25 onto a two dimensional space (Fig. 2.1), having one complex dimension. This maps the point (a,b) lying on the two dimensional complex circle of radius \ into the point X^ a/b= \2 e'^*^^^ cot T (2.19) on the complex line. The second step is to regard the complex point u = \2 gi (ap) ^q^ ^ as not lying on a complex line; but rather as a point in a real two dimensional plane. One now performs an inverse stereographic projection onto the surface of a three dimensional sphere whose south pole is tangent to the plane at its origin. By choosing the origin of the three dimensional space to be at the center of the sphere, a point on the surface of the sphere will be specified by giving the radius r, the azimuth 4> and the colatitude 9. For convenience the sphere is chosen to have diameter \ . The azimuth is measured in a plane parallel to the original plane and hence 4> = o p . (2.20) Thus to determine the location of the projected point P on the sphere, only the colatitude Q is left to be determined. From Fig. 2.2 one has tan (f ) = cot T (2.21) but tan (f ) = cot I (2.22) hence e = 2t . (2.23)
PAGE 32
26
PAGE 33
27 a. n 1. O 4)
PAGE 34
28 Thus the coordinates of the point P projected onto the surface of the sphere from the point w = X^ e'(^'P) cot t are e = 2t (2.24) = a p (2.25) r = \2/2 . (2.26) Expressing the quantities H, K, L and D in Eq. (2.4) in terms of these new coordinates gives H = \^n (2.27) K = (X^/2) sin cos = X (2.28) L = (\^/2) sin e sin * = Y (2.29) D = (\2/2) cos e = Z . (2.30) Hence, K, L and D simply determine a point on the surface of this sphere. It should be noted that K generates infinitesimal rotations about X, L about Y, and D about Z. Since these three constants of the motion are then directly the coordinates of the fixed point corresponding to the orbit of the oscillator, it is apparent that they generate rotations of the sphere, which transform one orbit into another of the same energy.
PAGE 35
29 Canonical Coordinates Whenever two functions satisfy the Poisson bracket relation U. B} HB (2.31) where p. i s a scalar, a corresponding equation, (A, (1/n) In B) = 1 , (2.32) may always be written, which is the rule satisfied by a coordinate and momentum which form a canonical ly conjugate pair. The first of these two equations simply states that B is an eigenfunct ion of A with respect to the Poisson bracket operation. Consequently it is to be expected that such eigenfunct ions be a ready source of canonical variables. A system in which the Hamiltonian is one of the momenta can be obtained in the present case. Consider the following diagram in which the twoheaded arrows imply that the Poisson bracket of the two functions is zero and the oneheaded arrows imply that the Poisson bracket is + 1 in the direction indicated: H = a a + b b'^ < > D = a a"^ b b (lAi) In {aWab)< > (lA') In (ab"Va*b) . (2.33) If these four quantities are rewritten in terms of the four parameters introduced in the Hopf mapping and p + o is denoted by ^, then Eq. (2.33) becomes
PAGE 36
30 ^/2 D = (X /2) cos 9 ^ 4)/2 (2.3^) Hence the Hamiltonian and the energy difference are the two momenta and the angles Y/2 and 4>/2 are the corresponding canonical coordinates. Since ^/2 is conjugate to H, "i increases linearly with time. The final result of this transformation is to map an entire orbit into a single point. Since the radius of the image sphere determines the energy, the only Information lost as a consequence of the mapping is the phase of the point representing the harmonic oscillator in the great circle comprising its orbit. The lost phase may be recovered by attaching a flag to the point representing the orbit, which will then rotate with constant velocity. The missing angle Y may then be taken as the angle which the flag makes with its local meridian. In this way the motion of the harmonic oscillator in its phase space is shown to be equivalent to the motion of a rigid spherical rotator, and in fact the parameters 4), and â€¢ are just the Euler angles describing the orientation of the rotor. These same parameters may also be used for the description of spinors, and a particularly lucid geometrical interpretation of this transformation may be found in a paper of Payne [32] describing various methods of representation of twocomponent spinors.
PAGE 37
31 The Anisotropic Harmonic Oscillator in Two Dimensions The anisotropic oscillator, both with commensurable frequencies and with noncommensurable frequencies has been discussed by Jauch and Hill [10] and also by Hill [33]These authors were able to show that, in the case of commensurable frequencies, the symmetry group of the system is again the three dimensional rotation group. However, they were unable to come to any conclusion concerning the case with incommensurable frequencies. Recently R.L. Hudson has been investigating the incommensurable case also. It can be shown though, that even with noncommensurable frequencies the symmetry group is also a three dimensional rotation group. The proof is sketched below and follows closely that used in the isotropic case, except where noted. If units are chosen such that the frequency in the x coordinate is unity, the Hamiltonian is H = (p2 + p2 + x2 + u)2y2)/2 , (2.35) where oj is the frequency in the y coordinate and is an irrational number. The eigenfunct ions and eigenvalues of H under the Poisson bracket operation are listed below: Private communication.
PAGE 38
32 Eigenfunction
PAGE 39
33 (2.30), however, for the anisotropic case Â» Cijp o . (2.42) A set of functions which proves to be convenient to use as the new canonical coordinates is diagrammed below with the same convention as before: In * * < *Â• In * 'Â— X kiui (a)*^ 4iu) (a)'^* (2.i*3) If the parameters from the mapping are substituted into Eq. (2.43) this table becomes identical to Eq. (2.34) with the conditions that 4> Is defined by Eq. (2.42) and ^ = up + a (2.44) It is rather interesting to note that the entire effect of the possible incommensurability of the frequencies in the two coordinates is completely absorbed by the power mapping (a)*^, which is well defined as exp (u) In a) in the case of an irrational u). In this latter case one has a mapping with infinitely many branches rather than a finite number as in the case of a rational w. It is of course on this account that the orbit Â— a Lissajous figure in the xy planewould then be open, and would not correspond to periodic motion.
PAGE 40
34 Conmensurabil ity and Constants of the Motion In reviewing the analysis of the theory of the plane harmonic oscillator, it is probably appropriate to draw attention to a number of statements made, both in the literature of classical mechanics, and the lore of accidental degeneracy. First of all , the existence of accidental degeneracy is often considered to be allied with the existence of bounded, closed orbits. When this argument is used, it is applied to the existence of algebraic integrals of the motion, and indeed It should be noted here that as long as the frequencies of the two coordinates are commensurable, bounded closed orbits result and the integrals are algebraic, as required by the theory. However, the above results also verify something which has been generally known [3^]. namely that nonalgebraic integrals can exist, and as would be expected, are associated with spacefilling motion. Another point regards the separability of the Ham i 1 tonJacob i equation in several coordinate systems. It is pointed out that if a HamiltonJacobi equation is separable in more than one coordinate system, it is necessarily degenerate [35]The argument usually given concerns the fact that orbits with incommensurable frequencies in one coordinate system will be spacefilling curves, the boundaries of which are coordinate arcs. Thus, in whatever coordinate system the motion is described, the system of separation is uniquely defined as that one forming an envelope of the various orbits. Of course, the option left in the enunciation of this theory is the fact that changes of scale may still be made, in which the new coordinates are functions each of only
PAGE 41
35 one of the old coordinates, and not several or all. So it is in the polar mapping, where one forms a ratio of powers of functions of only a single coordinate. However, this change of "scale" seems to eliminate the problem of the incommensurable frequencies at the same time, and thus circumvents the "proof 'of uniqueness. This result has only been secured at a price, namely the inverse of the transformation is infinitely multiple valued [36]. Otherwise, one could easily reduce the anisotropic oscillator to a static point on the Riemann sphere, return from this to an isotropic oscillator, and separate in whatever system desired. These considerations are not unrelated to the problem of finding a Lie group generated by the constants of the motion. In this respect, a Lie group has always been regarded as being more important quantum mechanically than classically, because its irreducible representations would be associated with accidental degeneracies of the Schroedinger equation. Since the anisotropic harmonic oscillator has no accidental degeneracies in the case of incommensurable frequencies, it is somewhat surprising to find that the unitary unimodular group is still its symmetry group, at least classically. In the case of commensurable frequencies, proportional say to m and n, the operator a'"b'''" commutes with the Hamiltonian, and is interpreted as creating m quanta in the first coordinate and annihilating n quanta in the second coordinate. Thus, m quanta in one coordinate are equal to n quanta in the second coordinate, the energy is restored to its original level, and the system is left in another eigenstate.
PAGE 42
36 Thus, it is possible to expect and explain degeneracy in the case of commensurable frequencies. In the classical incommensurable case, the symmetry group still exists. After all, the total energy of the oscillator is the sum of the energies in each coordinate, which is reflected in the fact that it is the sum of the squares of the semiaxes of the orbit (A Lissajous figure, like an ellipsoid, has principal axes), and not the individual semiaxes, which determine the energy. Thus, the dimensions of the bounding rectangle of a spacefilling Lissajous figure may be changed subject to this constraint. This, in addition to an adjustment of the relative phases in the two coordinates, is the effect of the constants of the motion acting as the generators of infinitesimal canonical transformations. The problem of carrying this result over to quantum mechanics seems to be in finding a proper analogue of the power functions (a) and b, which as ladder operators may not possess fractional powers, even though they did so when considered as functions of a complex variable. The Isotropic Oscillator in n Dimensions Using units such that both the mass and force constant are unity, the Hamiltonian for this system is H = V (P? + x?)/2. (2.45) i = l
PAGE 43
37 The eigenfunct ions of H under the Poisson bracket operation are aj = (Pj iXjO/Tz (2.46) with eigenvalues as wel I as with eigenvalues >^= i (2.47) a*= (Pj + ix.)//2 (2.48) (2.49) 2 Hence one can immediately construct n linearly independent constants of the motion of the form a;aV. i = 1, ..n and j = 1, ...n, although as before, it is somewhat more convenient to deal with their real and imaginary parts separately; in particular, n H =y aja* (2.50) Dj = a,a'f ajaf i = 2, ...n (2.51) K.J = a.av + a*aj. j > i (2.52) Ljj = i[afaj a.ap j > i . (2.53) The latter three sets of constants are closed under Poisson bracket and satisfy the following commutation relations:
PAGE 44
38 im > *VC . /**> (afak ajat)6.^ (aja, aja*)6.,] (2.55) {Oj. Dj) = (2.56) Mata,Ha.a*)6.,(ata,.a.ap6,^ (2.57) (L.j, D^) = (a^a^ + a.a*)6,^ (aja^ + aja^)8j,^ (aa, .a.at)5., (258) tDj. ^kn,^ = 't(^?k ^"k)^Tn * <^K ^j<)^k (ata, a,a;)6,,] . (259) Depending upon the values of 5, j, k, and m all of the relations will give linear combinations of the functions Dj, Kjj . or L;j as a result. The relations also identify the symmetry group of the problem [3?]. namely SUp where the n^ 1 constants are considered as the generators of the group. The Hopf Happing In general only 2nl independent constants of the motion including the Hamiltonian are expected to be found and not n2. Hence there should exist n^(2nl) = (n1)^ functional relationships among the n constants.
PAGE 45
39 As a consequence, the following set of 2nl constants have been chosen for the Hopf mapping; the (nl) constants K,j where j = 2, ...n, the (nl) constants Lj Â• where j = 2, ...n, and the following linear combination of H and the Dj n D= B,H +y B;Dj (2.60) i=2 where 6,=^ (2.61) I n and 2 n In terms of the a; and a'j B. = i (i = 2, ...n) . (2.62) I i=2 The (n1)^ relationships between the constants are D = a^a', ^ j, i=2 (2.63) 'J H + 'Â• H + D 2 2 2 (H+D)^ K,i L,. Dj = TTTT^T^ (2.66) 2(H + D)
PAGE 46
40 where H^' = D^ +) K,f +) L,f (2.67) iÂ«2 i=2 as may be verified by substitution from the definitions given in Eqs. (2.50) through (2.53). For the Hopf mapping let the 2n variables aj and a* be redefined as a, =jY^cos T)e'P (2.68) a. = ^ Â€.(sin T)e'<'j (j = 2, ...n) (2.69) J 72 J where the e: are direction cosines such that n j=2 (2.70) Forming the ratios = Â€j(cot T)e'^^"''J^ (2.71) gives a complex projective space of n1 complex dimensions. Regarding the space as having 2n2 real dimensions and performing an inverse stereographic projection a real space of dimension 2nl is obtained. As in the case of the two dimensional oscillator the sphere of 2nl dimensions is considered tangent to the 2n2 dimensional space at its south pole. By defining d as the angle between a radius vector from the center of the sphere and the diameter between the north and south poles 2t can be identified with 0. The angles in the planes
PAGE 47
perpendicular to the northsouth diameter projected from the 2n2 space parallel to this line are identical with the angles in the space and are denoted Â«!>!= aj p . (2.72) The substitution of Eqs. (2,68) and (2.69) in Eq. (2,63) and into the equations defining H, K) : , and L : gives H = (l/2)\2 (2.73) D = (l/2)\^ cos e = Z (2.7^) K,j = (l/2)\2 sin eej cos Oj = Xj (2.75) L,j = (l/2)\2 sin 0ej sin
PAGE 48
42 Other Canonical Coordinates In finding canonical sets of coordinates and momenta a wide choice of momenta, even when the n momenta are selected from amongst the constants of the motion, is available. The final decision will depend upon the particular application desired. However, there are some simple linear combinations which are often sufficiently useful, to make their display worthwhile. As a first selection for the canonical momenta consider the follown ing: p, = H =) ajaf (2.77) i=2 Pj =y ajat V aiaf (j 2. ...n). (2.78) i=l i=j where P. is the energy difference between the first (jI) coordinates and the last n(jl) coordinates. The coordinates conjugate to these momenta must then have the form (2.79) (j =2, ...n) . (2.80) In terms of the parameters used in the Hopf mapping these become H = \2/2 (2.81) Ql
PAGE 49
43 P2 = (\2/2) cos d (2.82) j1 n /^2,^^ r 2 \7 2 . 2 V 2 . 2 1 p. = (\ /2) [cos T+^ Â€Sin T\ eisin tJ 1=2 lj > 1 = (\^/2) [cos 9 +y esin^ 0/2] (j = 3. ...n) (2.83) iÂ»2 a, = l/2(p oÂ„) (2.84) Q2 = l/2(p 02) (2.85) d. = l/2(a._, a.) (j = 3, ...n). (2.86) One might also consider taking the total energy and the energy in the last n(jl) coordinates for the canonical momenta, n P, = H =y a. at (2.87) 1=1 and n P J . =^ ajaf (j = 2. ...n) (2.88) with the corresponding coordinates Q: = Â— In ^i^^ (j = 1. ...n1) (2.89) J 2, a^aj^, Qn = 1"(2.90) 2 1 an
PAGE 50
kk which in terms of the Hopf parameters become H = >^/2 (291) n p. = {\^/2)Y e? sin^ T (j = 2. ...n) (2.92) Q, = p 02 (2.93) Qj = o. Oj.^i (j = 2, ...n1) (2.9^) Finally, one might consider using the energy anci(nl) successive energy differences. n Pi = H^ aja; 1=1 n (2.96) P. = a.., a*, a, at j = (2. ...n) . (297) The corresponding coordinates are 1 a)a2 ...a^ Ql = 77 ^" 2ni aa2 ^n k1 ^* nk+1 n a"^"^ (2.98) Q,, = In ( ir (1) ) ( t (4) ) o 99) 2n' j = i ^j j = k ^j k = (2. ...n).
PAGE 51
^5 In terms of the Hopf coordinates p, = H = x2/2 (2.100) P2 = (\2/2)(cos2 T Â€2 sin^T) (2.101) p. = (\2/2) sin2T(ej., Â€j) (j=3. ...n) (2.102) n Ql = (I/n)(p *Y o.) (2.103) iÂ»2 Qj = (l/n)[(nj + l)(p^^ a,^)(J0^ ^k^ Â• (2.10A) k=2 kpj In all three cases, the transformation has been such as to give the actionangle variables. The effect has been to map the orbit into a point on the surface of a sphere of 2nl dimensions whose radius is the energy H. Since the mapping is many to one, one sees that as in the two dimensional case an entire orbit is mapped into a single point. The orbital phase may still be recovered by attaching a flag to the point representing the orbit which will then rotate with constant angular velocity Qi , since Q.i is conjugate to H and hence increases 1 inear ly with t ime. The second case also has another interesting property, namely D2 is the Hamiltonian for an isotropic oscillator in (n1) dimensions, and the other 0; are constants of its motion. Hence one may perform another Hopf mapping and continue in this manner until one has the 2dimensional oscillator again.
PAGE 52
1*6 In general, any of the 2nl independent constants from the n^1 constants could be chosen as generators, and used as coordinates. It should be noted that the expressions for the coordinates could just as well be expressed in terms of an arctangent by using the well known relat ion tan"' X = i i In LÂ±i (2.105) 2 1 + ix which would exhibit the similarity between these coordinates and those obtained by Goshen and Lipkin [23]. The Anisotropic Oscillator in n Dimensions In units where the mass is unity, the Hamiltonian for the ndimensional anisotropic harmonic oscillator is H = iÂ£ (P? + Â«^? x2). (2.106) i = l where uj is the frequency in the i coordinate. The eigenvectors are all of the form aj = (p. i w. x.)Af2 (2.107) a^ = (p. + i u,. Xj)//2 (2.108) with eigenvalues \j = Â± i ojj . (2.109) There also exists n^ constants similar to those obtained in the isotropic case, namely
PAGE 53
^7 II 5=1 * (2.110) Di = a..ai aj^, af^, (i = 1. ...nl) (2.111) (a; J a; ' Â•Â»Â• a, J aj ') (2.112) Uj1 Wjl *^2 /o =.*^ 2 (ajaf)^ (a^ap i[a*'^J a."' a"J af ' ] Â•I J Ujl U)j1 (a.af) 2 (a^aj)"^ Once again this set is closed under Poisson bracket and satisfies the commutation relations of the generators of SU^. It should be noted that every normal quadratic Hamiltonian is equivalent to some anisotropic harmonic oscillator, because they can be put in the canonical form given in Chapter 1. For nonnormal Hamiltonians such as the free particle H = l.P^ . (2.114) 2m there is no equivalence with the harmonic oscillator. However, symmetry groups still exist for these cases; for example the above expression is invariant under rotations. The symmetry of such Hamiltonians has been discussed in considerable detail by Moshinsky [38] and F.T. Smith [39, ko].
PAGE 54
48 For the Hamlltonlan given In Eq. (2.114) (where m 1) the following set of nine quadratic constants of the motion exists, H = (P^ + P^ + P^)/2 (2.115) X y z \ 'I Ai Pj Pj (2.116) A2 = Pj Pz (2.117) A3 = P^Py (2.118) A4 = P^P^ (2.119) A5 =Â• PyP^ (2.120) Li = yPz HPy (2.121) L2 = zPx xP^ (2.122) L3 = xpy yP^ . (2.123) The relations between these constants are summarized In the following Poisson bracket table:
PAGE 55
^9
PAGE 56
CHAPTER 3 THE HARMONIC OSCILLATOR IN A UNIFORM MAGNETIC FIELD Introduct ion In the previous chapter it has been shown how the accidental degeneracy of the harmonic oscillator may be accounted for by the methods presented in Chapter 1. Although historically, attention has most often been focused upon the accidental degeneracies of the harmonic oscillator [9] and the Kepler problem [k, 5], the problem of the cyclotron motion of a charged particle in a uniform magnetic field possesses certain unique features, among which is that It has two linear constants of the motion. It also has accidental degeneracies of its own [kl]. Its Ham 1 1 ton! an, _ H = (P I A)2/2m. (3.1) where, for one choice of gauge A= (^B^y, ^iB^x. 0) . (3.2) is quadratic. As shown in the previous chapter, the constants of the motion of the two dimensional harmonic oscillator K, L, and D arose in a natural fashion as products of linear eigenfunct ions of the operator T^. The eigenvalues of Tm occurred In negative pairs, so that their products, belonging to the sums of the corresponding eigenvalues, were constants 50
PAGE 57
51 of the motion. Moreover, there was a convenient mapping, the Hopf mapping, which yielded a set of canonical angular coordinates and ^, whose conJugate momenta were respectively, H and D. In the present instance, one finds by expanding Eq. (3.1) and using the gauge in Eq. (32) that H J. (p2 ^ p2) + liÂ£ (x2 + y2) ^ f?o (yp^ . ^p^j (33) 2m 8mc 2mc Effectively this is the sum of a harmonic oscillator Hamiltonian and a term proportional to the angular momentum. There is one main difference between the problem of cyclotron motion and the harmonic oscillator. The representation of the harmonic oscillator Hamiltonian, applied to the homogeneous space H^ ' has negative pairs of nonzero eigenvalues. All constants of the motion are generated by products of the corresponding eigenfunct ions , for which the eigenvalue sum is zero. Thus the lowest order constants of the motion are quadratic, belonging to h' ', and all others are linear combinations of these products. With the cyclotron Hamiltonian there are zero eigenvalues, and therefore there are elements of h' ' which are constants of the motion. Whereas the Polsson bracket of two quadratic constants is quadratic, the Poisson bracket of two linear constants Is a scalar. A rather different kind of symmetry group results In the two cases; for the harmonic oscillator one finds the unitary unlmodular group, while for cyclotron motion one obtains a group generated by the harmonic oscillator ladder operators.
PAGE 58
52 Classical Zeeman Effect for the Harmonic Oscillator It Is somewhat more instructive to consider a plane isotropic harmonic oscillator in a uniform magnetic field rather than the cyclotron motion exclusively. Its Hamlltonian, H i_ (p Â£ A)^ + imuj ^2 2m 2m (p2 + p2) + El (0^ + co2)(x2 ^ y2) + ^yp^ . ^p ) (3.^) Is also quadratic. In Eq. (3^) w Is the Larmor frequency eB/2mc and Wq is the natural frequency of the oscillator. This Hamlltonian reduces, In the limit as cj Â» (B Â» 0) to that of the plane harmonic oscillator, while as co^ Â» It reduces to that of cyclotron motion. The matrix representation of this Hamlltonian considered as an operator under Polsson bracket and denoted by T^ Is
PAGE 59
53 Eigenvector
PAGE 60
5^ ^"'^ 2 2 1/2 In order to display the symmetry group in a convenient form the following linear combinations are tatcen as the constants of the motion: H = [\, uu* + \2 vv*]/i+ yj+ J^ (3.15) R1 K = [u^v* + u*'^v]/r'/2(uu*) 2 (3.16) R1 L = i[uRv* u*Rv]/r'/2 (iju*)"2~ (3.17) D = [uu'^ Rvv*]/R , (3.18) where H has been written in accordance with Eq. (1.29). These four quantities satisfy the following relations: Tj^(K) T^(L) = T^(D) (3.19) Tk(L) = aD (3.20) TjD) aK (3.21) T[,(K) = ai, (3.22) a Â« 8^u2 + 0)2 . (3.23) where The first of these equations Is simply a statement of the fact that K, L, and D are constants of the motion, while the latter three equations show that the symmetry group of the system is SU2.
PAGE 61
55 If in Eq. (312) R is a rational number, then the classical system has bounded closed orbits and a quantum mechanical analog exists for the operators K and L. However if R happens to be irrational, then the orbits are space filling. Limiting Cases In taking the limit as uj Â»0, i.e. as the magnetic field is turned off, Eq . (3.^) becomes the Hamiltonian for the two dimensional isotropic harmonic oscillator, the eigenvectors become, u Â«"/m Uq r"*" + LP* 0.2k) V m V NTm cOq r+ Lp, (3.25) V m while the four eigenvalues degenerate into two, namely Â± IUq, whence R = 1 . Expressed in terms of the new u and v, the four constants are H (uu* + vv*)A (3.26) K Â» uv* + u*v (3.27) L i(uv* u*v) (3.28) D uu* vv* , (3.29) which are the constants previously obtained for the isotropic oscillator. The commutation rules for K, L, and D also still hold.
PAGE 62
56 In considering the limit as u)q Â»0 one finds that and u = NTm w r* + U P+ (3.30) V m V = fm a)r~ + LP". (3.31) The eigenvalues of u and u* approach Â± 2iw in this limit while those belonging to V and v* both approach zero. Hence, in order to satisfy Eq. (3.12) R must also approach zero. The Hamiltonian in this limit approaches that for pure cyclotron motion as in Eq. (33). From the values of the eigenvalues one immediately has two linear constants of the motion, v and v*, and one quadratic constant uu*. The constants can be explicitly derived by considering the commutation rules of K, L, and D in the limit as uJq Â»0. Rewriting the commutation relations explicitly gives , ruRv* + u*^ u^* u''^ 1 ^ uu^ Rvv^'^ ,. Â„. ' uRy* u*K uu* Rvv* ^ Â„ a u^* + u^'^v (333) Rr ( Rr and r'/2(,,*) 2 , J R>/2(,,*) 2 uu* Rvv* u'^v* + u*'^v I _ . ii^v* 11*'^ R r'/2(iju ) 2 J r1/2(uu ) 2 Multiplying the first of these by R and the latter two by R^^ and then
PAGE 63
57 taking the limit as cJq +0 (R 0), results in the equations I ((uu*)'/2 (v'Vv), (uu*)'/2 (v*v)} o 8 a)uu* (335) I {(uu*)>/2 (v*v). uu*} = (336) {uu*. (uu*)'/2 (v* + v)] . (3.37) Since uu is simply twice the cyclotron Hamiltonian, it follows from the last two equations that both the real and imaginary parts of v are constants of the motion. Dividing Eq. (335) by uu* gives ^^J1, 1^^ . 8U, . (3.38) Expanding Eq. (31) results in H = (P^ + p5)/2m + mJ(x^ + y2)/2 + u, (yp xP ) . ^ (3.39) Thus the cyclotron Hamiltonian splits into two parts, one being the harmonic oscillator Hamiltonian Hq and the second proportional to the Zcomponent of the angular momentum L. Both of these terms commute with the total Hamiltonian H = Hq + LHence another quadratic constant of the motion is D ' Ho L . (3.^0) For convenience in notation define S "^ (vv*) . muy P^ (3.m) r Q ^ (vv*) mux + Py . (3.^2)
PAGE 64
58 With these definitions the following commutation rules hold: (H. D) = {H, S} (H, Q} (3.^3) {S. Q} = 2maj (3./+4) {D. S) 2coQ (3.Zi5) (D. Q) Â«= 2a)S . (3.^6) These commutation relations coincide with those obtained by Johnson and Lippmann [kl] . These authors have discussed the two constants S and Q in considerable detail and show that they are simply related to the location of the center of the circular orbit and to its diameter, which can readily be seen In the following manner. The canonical momentum expressed in terms of the mechanical momentum is P = m V + Â® A . (3.^7) Substituting for the canonical momenta in Eqs. (3^1) and (3.42) gives S = m(2wy v>^) (348) Q m(2wx + vy) . (3.49) Evaluating S when Vj^ Â™ and Q when Vy = gives the center of the orbit as (xc. Yc) " (Q/2m(J, S/2mw). (350) These constants also determine the diameter of the orbit. Since the orbit is a circle only one of the constants must be considered. For
PAGE 65
59 example, consider Q. When v takes on its maximum positive value, x takes on its minimum value and when Vy takes on its maximum negative value, X is a maximum and hence the diameter d is dWmin = ^^(3.51) u Because of the continuum of points available for the center for a given energy the degeneracy of this problem is infinite. Canonical Coordinates As In the case of the plane harmonic oscillator a set of canonical coordinates can be found such that the Hamiltonian becomes a canonical momentum. In fact, two momenta for the problem are: Xi uu' + \o vv' H = ' (352) ., 2 2,1/2 ^^ and D = HHIjlJRvv: ^ (3.53) while the coordinates conjugate to these momenta are Q, = \2['n Â— + >.i In ^]Ai \,\2 (35^) u V Q2 X2[ln y^ \, In ^]/l6i \2 (u)2 + u,2)'/2 (3.55) respect i vely,
PAGE 66
60 A mapping similar to the Hopf mapping may also be performed where u = [\^/2 cos T e"'P] (3.56) and V = [\]/2 sin T e'*'] . (3.57) Under this mapping the momenta become H = X^X^/k {J+ u)^)'^^ (3.58) D = \] cos (359) where e Â•= 2t (3.60) and the coordinates are Q, = ^/2^2 (3.61) ^"'' 2 2 1/2 Q2 = 'J>/8r(u^ + %) ' (3.62) where Y = Rp + o (3.63) and Rp o . (3.64) Performing the mapping on the other two constants of the motion gives K = X.] sin cos (J> (365) and L = \, sin sin
PAGE 67
61 and P = [(p2 + p2)/2m + m(u)2 + J) {k^ + y^)/2]/{J+ u)^) ' /2(3 .67) P0 = xPy yP^ The corresponding coordinates are (3.68) and 0=1 tan"' 2 q = 1 tan"' 2 [P P /m + m(u2 + cj2) xy] X y 2_ o (Px + Py)/2m + m(uj2 + '^o)^^^ + V^^^^ (0)2 ^ 0,2)1/2 (^p^ ^ yp^) ^ I (Px + P5)/2m m(uj2 + ui2j(^2 ^ ^2^/2 (3.69) (3.70) In terms of these variables the Hamiltonian has the particularly simple form H ((^ + 0,2)1/2 p^ . ^^^ (3.7,) from which it follows that both P^ and Pq are constant in time. Rotating Coordinates The problem of the plane harmonic oscillator in a uniform magnetic field has a certain uniqueness when viewed from a rotating coordinate system. However, the problem will first be solved in plane polar coordinates. An extremely lucid description and tabulation of these orbits has been given by E.R. Harrison [^3 ] Â• Assuming the direction of the magnetic field to be in the negative Z direction, the Lagrangian is
PAGE 68
62 2 L = (r2 + r202) . !I!^ r2 mwr^Q (3.72) where the following gauge has been chosen, Ar =Â» A^ = (3.73) and The momenta conjugate to r and are Pr = "If(3.75) and . . Pq mr (6 u). (3.76) Since L Is cyclic in 9, Pq Is a constant of the motion. The Hamiltonlan Is = (Pj + p/r2)/2m + m(u)2 + wj) r^/2 + uPq . (3.77) and the equations of motion are: Â„ 2 Pe 2 2 p^ = _Â£_ . ni(cj + u)o)r (3.78) mr^ (3.79) r = Pr/m (3.80) e = _^+ u . (3.81) In general, the effect of imposing a uniform magnetic field on a system
PAGE 69
63 with a central potential is to add two terms to the Hamlltonlan, namely a harmonic oscillator potential which is often neglected for small fields [MÂ»], and a term proportional to Pg, the angular momentum. The solutions for the orbit equations are: _a sin"' [ br 2a i + Â±_ sin"' ( ^^ ^ ^ ] 2 la I r2[b2.4aV]'/2 2W [b^^aVj^/^ a br^ 2I "* . 1 ^^''o '^ sin"' { ^ Â— s 701/9 ) sin"' { Â— 5 55 \ Jn ^ 2Â« rJ[b2.4aV]'/2 2W [b^UaVl'^^ (3.82) where 2 Â« Â« b = Â£_+ w2r2 + r2 (383) a = P0/m (3.8^) W= (0,2 + 0,2)1/2 . (3.85) The subscripts on the coordinates and velocities denote initial values. These orbits are plotted in Fig. 31 through 3^. I" al ' cases units have been chosen such that m <Â» q = c = 1. Fig. 31 and 3.2 show the high and low field orbits respectively, for the same set of initial conditions and Wq. Fig. 33 and 3^ show the orbits for a fixed magnetic field, but for different values of the initial tangential velocity. In general the orbit will be a processing ellipse for Pq ^ and will be a hypotrachoid similar to Fig. 31 for Pq < 0.
PAGE 70
64 Parameters Initial Conditions u) = 63A 10 Â« 4 o e = o r o Fig. 31 Particle orbit in a strong uniform magnetic field and harmonic oscillator potential.
PAGE 71
65
PAGE 72
66 o >. Â•Â— *J c Â— o o E O m 0) ^ > o Â— _
PAGE 73
67 c
PAGE 74
68 In a rotating coordinate system defined by e = e cot (3.86) r r , (3.87) where w is the Larmor frequency, the Lagrangian is The canonical momenta in the rotating coordinates system are: P= mr (3.89) P= m?^5 . (3.90) o L is also cyclic in and hence Pis a constant. The Hamiltonian in the rotating coordinate system is H = P^r* PL 2 Â„ 2.2. ,, ,2 2.2 + PS /r )/2m + m(u) + o) ) r o ^ o (P; + Pq /? )/2m + m(u + <)r /2 (391) and hence the equations of motion are p. = _L m(cj2 + cj2)p (3.92) r _9 u Pg = (3.93) ? = P?/m (3.9M e ^ . (3.95) mp2
PAGE 75
69 which are the equations of motion for a plane Isotropic harmonic oscillator with force constant m(u)2 + u^) . In general the transformation of a Hamiltonian with a central force from a stationary coordinate system to a rotating coordinate system will subtract a term proportional to Pq. Generalizing then, it can be stated that a Hamiltonian with a central potential and with a uniform magnetic field present, when viewed from a rotating coordinate system, has a harmonic oscillator potential added, that is H = f + V(r) + 1 mu)^ r^ (396) where T Is the kinetic energy in the rotating coordinate system and u Is the rotation frequency of the coordinate system with respect to the fixed system. In the present case since \/(r) is a harmonic oscillator potential one simply has a harmonic oscillator with a larger force constant. The solution for the orbit defined by the Eqs. (3.92) through (3.95) is . a br^2a^ a , br^2a^ eeo'Â— sin' { , , , ^ ^ } sin' [ , Â° ! . ... ]. 2a r^Cb^^aV]'/^ 2a r^[b^^aV] '^^ where W Is defined in Eq. (3.85) and a = Pg/m (3.98) b = a2/r2 + W^ r^ + f2 ^3 g^j Equation (397) Is the equation for an ellipse, regardless of the
PAGE 76
70 magnitude of o). The orbits corresponding to Fig. 3.1 and 3.2 in a rotating coordinate system are shown in Fig. 35 and 36. Discussion of the Constants of the Motion Of the four constants of the motion defined in Eqs. (315) through (3.18) one is H, the energy of the system. Another constant is the angular momentum Pq which is a linear combination of H and D Pe^ JiLH ' ^ ^ D. (3.100) Since only three of the constants are independent one would like to discover one other constant, independent of the two above, which has a physical meaning. This can be done by considering the constants of the motion in the rotating coordinate system and then transforming back to the original coordinate system. In the rotating coordinate system two constants of the motion are aa* and bb'"', where a = Pjj \{J+ u)2)'/2 (3.101) and b P1(0)2 ^ ^)l/2 ^ (3.102) and where the mass is taken to be unity. By a simple rotation the axis can be oriented so that the semimajor axis of the ellipse lies along X and the semlmlnor axis lies along y. In this case then aa' is pro2 portional to the maximum value which x attains, because since
PAGE 77
71 o I CD
PAGE 78
72 D
PAGE 79
73 aa* = P2 4(uj2 + cj2) ^2 (3.103) is a constant, it can be evaluated when P. = 0. Similarly bb* Is proportional to y ^ . Alternatively, one can say that aa and bb'' give the boundaries of the orbit i.e. aa' gives the maximum radial distance the particle can attain and bb* is the minimum radial distance. These two radii are also constants of the motion in the stationary system since the transformation only involved the angle. However, it can explicitly be shown that these are constants. By transforming from the rotating coordinate system to the stationary one aa''' can be written in terms of the eigenvectors and their conjugates as aa"'^ =(l/i+)(uu'"' + vv"'^ + uv'^ e^''^^ + u*v e'^'^^tj ^ (3.104) In the stationary system the time rate of change is ^ (aa*) (aa*, h} + ^ (aa*) (lA)[2iaxjvV'^t + 2iu)u*ve2''*'t +2iu)uvV'^t 2iuju*ve"2''^t] Â» 0, (3.105) and hence aa* is a constant of the motion in the stationary system also. A similar proof also shows bb* to be a constant. However aa* and bb* are not independent, but are related to one another through H and D. Even though a good physical interpretation cannot be given to the constants K, L, and D their effect on the orbit under Poisson bracket can be calculated i.e., the infinitesimal change each produces.
PAGE 80
7^ Before calculating the effects of the constants on the orbit, it is desirable to write the orbit equation in terms of the constants of the motion and the time. This can be done by defining a complex vector 7 = X + iy (3106) which can be written in terms of the eigenvectors by inverting Eqs. (3.9) and (3.10) r Â» (u + v'0/2c , (3.107) where c = [m(u)2 + u^)]'/2. (3.108) Employing the polar forms of both the eigenvectors in Eqs. (3.56) and (3.57) and of the constants In Eqs. (359), (3.65) and (366), u and v* may be expressed in terms of the constants as u = (R/2)'/2 [(k2 + l2 + 0^)1/2 ^ D]'/2 e'P (3.109) v*Â«. (R/2)'/2 [(k2 ^ L^ + d2)1/2 . q]W2 ^\a ^3,10) However from Eqs. (363) and (3.64) p = (1/2R)(Y + ) . (3.112) Hence u (R/2)'/2 [(k2 ^ l2 ^ d2)'/2 + D]'/^ e'*/2R e'^/^R (3.113) and v*.(Fl/2)'^' Uk' * l' * 0')''' . 01''' e'"^ e""^ . (3.114)
PAGE 81
75 From Eqs. (365) and (366) It Is apparent that D2)l/2+D]'/2e(l/2R)tanl (L/K) e'i^Ktto) 2c 2 + [(K^ + L^ + D2)'/2.p^]/2^(i/2)tanl (L/K) ^i^2(tto), = '^+ ^ . (3.119) where 2c 2 (3.120) 7_ . 1_ (a)'/2[(K2 + L^ ^ D2)'/^D]'/2e('/2R)tan' (L/K) ^U,(tto) and 7^ = i_ (R)'/2[(k2 + L^ + D')'/2.o]'/23('/2)tanl (L/K) J^2(tto). 2c 2 (3.121) Equation (3.119) expresses the orbit equation as the sum of two contra rotating vectors, whose frequencies of rotation are \j and Xj, and
PAGE 82
76 whose amplitudes are functions of the constants of the motion. The Poisson bracket of the constants K, L, and D with 7 are: (K. 7) = ^ 2 IKD [(K^ + L^ + d2)'/2.d] k^ + L^ iKD [(K^ + L^ + D^)'''^D] R(K^ + L^) (3.122) (L, r) a + Â— r/ 2 2 2J/2 , ,2 2, [(K + L + D ) +D] R(K ^ L )_ (3.123) (D. r) \a _ Aa _ 2R ''V 'Â» (3.12^) where a is a constant defined in Eq. (323). The infinitesimal change induced by the constants is that each changes both the amplitude and phase of r and r_, while preserving the sum of the squares of their sum and difference. If the orbit were an ellipse, this last statement would be equivalent to stating that the sum of the squares of the semi axes of the ellipse is a constant. The proof of the statement goes as follows: Before the infinitesimal change the quantity is Ir^+ r.2 + r^ rj^ = 2(rj2 + rj2) (3.125)
PAGE 83
77 and after the change it is r^(l + Â€) + r.(l + 5)2 + r^(] + e) r.(l 5)2 2[r I^ + rj^ + (Â€ + Â€*) lr_^2 + (6 + d'") rj^] (3.126) where second order terms have been ignored. The difference between the two terms is ^ o u o 2[(Â€ + e') r^2 + (5 + 5*) irj2] . which is zero in all three cases. For example, if the infinitesimal change is induced by K, one has for the difference 2[(e + 6^)r^2 + (6 + 6*)rj2] 2a 2 2L frj' 2L r.r [(k2.l2.d2)'/2.d] [(k2.l2.d2)'/^d: Â„ 2aR " 8c 0. [(K^ + L^ + d2)'/2.d] [(^2 + L^ + d2)'/2+oJ (3.127) Gauge Transformations in Uniform Magnetic Fields In uniform magnetic fields and in the absence of other external potentials the following fact concerning a charged particle in this field is true: The change in momentum in going from one point to another Is Independent of the path taken between the points. The
PAGE 84
78 proof is as follows: Let P denote the canonical momentum, ir the mechani cal momentum, B the uniform field, and A the vector potential such that and since B is uniform Since B " Vx A (3.128) A i B X r . (3.129) ^^ B.7 xB . (3.130) at c and one has dt / dn = J / (f X B) dt ^2 _ _ = ^ / d? X B c Â•' ^2 c r Hence and c (3.131) ^/ dA . (3.132) 1 Art " ^ AA (3.133) c AP ^ AA . (3.13M Equation (3.134) states that the difference in momenta between any two points is proportional to the difference in the vector potential evaluated at these points. From this it is inferred that the choice of a certain gauge is equivalent to picking the zero of momentum.
PAGE 85
79 One can also see why a translation of coordinates must be accompanied by a gauge transformation. When the translation occurs the zero of momentum changes and hence this change must be subtracted from A to give the the same zero of momentum. However, the change in A is simply a gauge transformat ion. In order to show that a gauge transformation can induce a translation, first consider the case of a free particle in a uniform field with the following gauge A^(y x 0) 2 (3.135) In units where m and co eBQ/2mc are one, the constants of the motion are Px y = Â« (3.136) and Py + ^ = 8 Â• (3.137) Using the fact that P^+fA. (3.138) the location of the center of the orbit is found to be ^\' Vc) " (f. f) Â• (3.139) However with the gauge D A 2 ( (\ + y). ^i + X, 0). (3.1^0) even though the two constants are the same, one finds the center of the orbit to be ("c Vc)(2^ ^) . (3.141)
PAGE 86
80 Hence the center of the orbit has been translated g,/2 units In the x direction and \/2 units in y. The center of the cyclotron orbit may also be rotated, however this Is not accomplished by a simple gauge transformation. Instead one must rotate the vector A and also rotate the coordinate system. For the vector potential defined in Eq. (3135) the center of the orbit is given In Eq. (3.139). A rotation of the vector field defined by A g Ives h'f cos e sin e sin cos d 1 2. { y cos 9 + X sin d, y sin + x cos 6, 0). (3.1'*2) The rotation of the coordinates defined by and X X ' cos y ' s in y"x' sin0+y' cos (3.143) (3.1M*) transforms Ad into A' where AÂ» _Â° (y', X', 0) . (3.1^5) Hence, the center of the orbit Is and (x' , y ') ( Â— , ) ^ c c 2 2 S Â„ 1 cos + 1 sin 2 2 2 (3.146) (3.147)
PAGE 87
81 1= . 1 sin e + 1 cos 0, (3.1^8) which shows that the center of the orbit has been rotated with respect to the original coordinate system, while the diameter of the orbit has remained unchanged. The vector field can also be rotated without rotating the coordinates and still preserve the curl If at the same time the vector is dilated by a factor of sec 0. The new vector potential is D A = Â— sec 9 (y cos 0+xsin0, ysin0+x cos d, 0). 2 (3.149) The transformation leaves the center of the orbit invariant and in fact, Is equivalent to a gauge transformation where the vector A^ A. = ^ tan e(x. y. 0). (3.150) whose curl is zero, has been added to the original vector potential A.
PAGE 88
CHAPTER k THE KEPLER PROBLEM The Two Dimensional Case in Polar Coordinates Tiie Ham i I ton i an for tlie two dimensional Kepler problem in polar coordinates is H Â« (1/2)(P; * p;/r2) . 1/r {k.\) where the mass and force constant have been chosen equal to unity. The corresponding Hami 1 tonJacobi equation is (^,^ . ' (^)^ 2(W . 1) = 0. (^2) br r^ ^ r where W is the total energy and S(r, ) . (4.3) Hence p Â«^ (^Â•^) and '* ' i ("5) Since the partial differential equation (^+.2) has no dependence bS/^ " P(K Is a constant and is, in fact, the angular momentum of the system. Solving Eq. (42) for 3s/3r Â» P^. one finds S Â„ p [2W+ 2 li] (4.6) or r r'^ 82
PAGE 89
83 The action variables J^ and J^ are given by tlie integrals Jr = 55 P^dr (4.7) J^' P P4,d]/2jr, (4.14)
PAGE 90
m where the elliptic anomaly u Is defined through r = a(l Â€ cos u). ('Â•15) a being the semimajor axis and Â€ being the eccentricity of the ellipse, a Â« W/2 (4.16) Â€^ =. 1 + 2wpJ . (4.17) By choosing to be the azimuth of the perihelion, X = r cos (4" ^^) (4.18) y = r sin ( *)Â• ) , (4. 19) and Eq. (4.14) reduces to w. = (u Â€ sin u + <(> )/2jt . (4.20) (p o Solving Eq. (4.10) for the energy in terms of J and J^ , 2 W= ^ (4.21) (Jr + V and hence the frequencies v and v^ are degenerate i.e., r <*> ^= ^ = ^= ^ Â• ^"^""^ Thus the quantities aa' , bb , ab' , and a b are constants of the motion where J 1/2 2jtiw a = (2^ e ' (4.23)
PAGE 91
85 and b^A'/'e^"'^. (4.2^ 2 It These quantities satisfy the following Poisson bracket relations {a, a*) = (b, b*} I (^.25) where one uses the fact that the J's and w's form a set of canonical coordinates and momenta. Taking the real and imaginary parts of these constants and the sum and difference of the real ones gives the following four constants of the mot ion, H Â» aa* + bb"'^ (^.26) D = aa''' bb''^ (^.27) K = ab* + a'*b (4.28) L f(ab* + a*b). (4.29) Of these constants H commutes with the other three and is a function of the Hamiltonian, namely, it is inversely proportional to the square root of the energy. The three constants K, D, and L satisfy the comtation relations of the generators of the group SUo or R,. mu (K, 0} 2L (4.30) (D, L) " 2K (4.31) (L. K] = 20 . (/4.32) Although they have no simple interpretation as constants of the motion, one can make a mapping of these constants which preserves the commutation rules. The new constants can be interpreted as the angular
PAGE 92
86 momentum and the two components of the Runge vector [11]. The mapping (^.33) K .K' = K(2 + HÂ±D)l/2 HD L .^L' L(2 + HÂ±D)l/2 HD (k.3k) D Â»D' = H D . These constants satisfy the rules {K', D'} = 2L' (D', L') = 2K' (L', K'} = 2D' (^.35) (4.36) (^.37) (^.38) Otherwise they are a direct consequence of the commutation rules for K, L, D, so that this mapping is valid for any SU2 group where use must be made of the identity h2 = K^ Â» L^ + d2 . (4.39) K' Is essentially the x component of the vector P where P = RAT^ = (cATTi) cos <)>, s in . (i+.40) and R is the Runge vector, while L' Is the y component and D' is the angular momentum.
PAGE 93
87 The Kepler Problem In Three Dimensions Using the same units as in the previous case, the Hamiltonian for the Kepler problem in three dimensions Is 2 ? ' 2 Pe P* I "(Pr^7 % \ ) Â• (^.^0 2 r^ r^ sin2 e r Setting up the Hami I tonJacobi equation and separating gives as ' ' ae V ^sln2 6/ and as 2 ai P^ Â» .= (2W + ^)'/2^ (1^ i^) where Oq and o^ are separation constants, a^ being the magnitude of the Z component of the angular momentum and Oq the magnitude of the total angular momentum. Since the potential is central the orbit Is planar and the normal to the plane Is necessarily parallel to the angular momentum. Denoting by B the dihedral angle which the orbital plane makes with the xy plane one can deduce that cos B = ^ . (4.^5) Solving for the action variables gives J4, Â«. 2na^ 2nP^ {k.kS)
PAGE 94
88 Jq 2(00 a^) 2 2W J0 J4, (k.kQ) and from Eq. (4.48) It follows that 2^ (j^ + Je + J({,)' (4.49) In the plane of the orbit rectangular coordinates Mand v can be introduced, where the \l axis is the major axis of the ellipse and the origin coincides with the center of force and the origin of the polar system. In this system H = r cos (Y Yq) (4.50) V = r sin (^ ^o), (4.51) where ^ is an angle measured from the Maxis to the particle and â€¢q Is the angle to the perihelion. The angle variables are 3Pr w_ => / Â— dr = (u e sin u)/2n (u e sin u tan1 y. + sin"' S2^)/2n ^ [I sin B (4.52) (4.53)
PAGE 95
89 (u e sin u tan' ^^ + sin"' SÂ£L^ . sin"' ^^ig + 0)/2Â« ^ Msin B tan B ' (^.54) where u Is defined as In Eq. (^.15) and Â€ as In Eq. (4.17). Using Eqs. (4.50) and (4.51) it follows that tan^ V . ^j, . Y (4 55J It follows from Fig. 4.1 that sin"' SÂ°l2. is the projection of the tan B angle H onto the plane in which is measured i.e., the xy plane, which shall be denoted^ . It also follows that sin"' Â£2LÂ£ Is the angle 1, sin B where the line of nodes has been chosen to coincide with the major axis of the ellipse. Hence, the line of nodes makes an angle (4> ^) with the Xaxis. Rewriting Eqs. (4.53) and (4.54) results in the following set of equations, v^ = (u e sin u + ^o)/2n (4.56) W4, = (u e sin u + ^Q + <^ ^)/2jt . (4.57) By defining a = (L) e "^ (4.58) 2it b Â» (:i^)'/2 e2Â«Iv^ (4.59) In J* 1/2 27tIwij, ^, ^ ^ c = (i) e ^ (4.60) In
PAGE 96
90 Fig. ^.1 Three dimensional coordinate system for the Kepler problem.
PAGE 97
91 one finds the following nine constants of the motion: aa*, bb*, cc*, ab , a b, be , b c, ac , and a'c. These constants are the generators of the group SU^ . The constant angles B, ^^ and <> ^ and also the sum ^ + Hq can be expressed in terms of the constants, sin ^o iii.^Â— ^iM. (4.61) Â° 2(aa*bb*)'/2 ^ ab 4a b ,, ^Â„, cos ^o fTf (^.62) 2(aa*bb*) slnC^. '\'^l\f,% (^.63) 2(bb'*cc''')'/2 cos(Â» ^) = ^^\\^y, (^.64) 2(bb*cc^')'/2 in(^ +YJ = '(^"^^ ^y,l Â° 2(aa*cc*)'/2 (^Â•65) and cos(..^.,Â„). ^Â— ^,, C^^ COS B Â£S ^ . (i,.67) bb* + cc* Generally, the rotation group in four dimensions R^, is tal
PAGE 98
92 Â•2W For L one has cos ^o cos ( cos ^Q sin ( sin B sin ^Â„ ") cos S s i n ^Q sin ( ^) ") cos S sin ^Q cos ( ^) (4.68) In sin B sin ( ^) s in S cos ( ^) cos B (4.69) Using Eqs. (4.61) through (4.67) one can express the generators of R/^ In terms of the generators of SU3 as follows: ^'"!^f?Vy^'^! I(ab* . a*b)(bc* * b*c)(bb* + cc*) "" 4bb*(cc*)'/2(bb'' + cc*) + (ab* a*b)(bc* b*c) cc*] (4.70) p . l(aa*H2bb*.2cc'^)^^^ ^3t^.v ^ a*b)(bc* b*c)(bb* + cc''') y 4bb*(cc*)'/2(bb*+ cc*) + (ab* a*b)(bc* + b*c) cc*] (4.71) 1 /o i(aa* > 2bb* + 2cc*) ,,,k . Â„ Â•>> 1 /o, .* v.v ,, ,,v P Â» i Â— Â— T ;Â— ^ (bb +2cc)'/'^(ab a b) (4.72) ^ 2(bb* + cc*") i(bc* b*c)(bb* * 2cc*)^^^ 2(cc*)'/2 (4.73) _ (bc*^b'c)(bb^H2cc'^)'/2 Â•y 2(cc*)'/2 (4.74)
PAGE 99
93 and L^ cc* . (4.75) These six quantities satisfy the relations ^^l' Pj5 ^ijk Pk ('+Â•76) a,. Lj} Â€,.^ L, (ii.77) ^Pi' Pj^ ^ijk k (^.78) which define R^, once again directly from the commutation rules for SU3, The process can be inverted and the generators of SU3 can be expressed in terms of the components of P and L. These are ..*.^.(,2,,V/2.,^ (..79) N/2W bb'^ML^.L^)'/^ (4.80) cc* L^ (4.81) '* " IT^ " ^'' * '?^"'^^^"'<^x l2)>/2(cos 1^ i sin ^ 2W (4.82) be*tl^d^ + l2)'/2]1/2 (^^5 (^ .^) . , sin (4, ._)) (1^83) ac*[7=i=(L2 + l2)'/2.l^]1/2lJ/2[,,3(, .^.,^^).i sin(* * ^^ where 2w (i*.84) "x ^ ^' tan"' ^ (4.85)
PAGE 100
94 ^Q tan P^LydJ + d ^ l2)'/2 Â•1 ^z'y P (L^ + L^) + P L L ^x^ X y' z X z and S = tan ^2 ^ ,2J/2 1 (Lx + Ly) (^.86) (4.87) The Two Dimensional Kepler Problem in Parabolic Coordinates In parabolic coordinates, defined by X a ^i V y = (n^ v2)/2 and where m =Â» e = 1, the Kepler Hamiltonian is Pu + P^ 4 H [i ^ ry 2(^2 + v2) From the Hami 1 tonJacob i equation where P, = (2Vii2 + 2a,) ^/2 P = (2Wv2 + 20^)^/2 tt) + a2 = 2, The action variables defined by J, ^ Pjdq, give iL 2no!i (4.88) (4.89) 2W (4.90) (4.91) (4.92) (4.93) (4.94) (4.95)
PAGE 101
95 and 2itZ2 (^.96) Solving Eqs. (4.91), (492), (4.95) and (4.96) for a^ and ^2 ^'^^ equating the solutions and Hence "^^ JH (Pj 2Hx^)/2 (^97) 2n 2W J = (p2 2Wv2) p2 J jtsT^ (v^ + Â— i:L) "* (2W) (4.98) (4.99) (4.100) The form of these latter two equations suggests that two new variables be defined a = J. + Â•2W M(4.101) b " V + Â— Py , so that ^ jtN/2W aa' Jy Â» nN/2w bb* . (4.102) (4.103) (4.104) Both a and b are very similar to the harmonic oscillator operators introduced in Chapter 2 and by analogy one would expect to find that quantities proportional to ab* and a^b are also constants of the motion, This is in fact the case. Taking real and imaginary parts of these quantities and normalizing, the constants are
PAGE 102
96 K ^ (ab* + a*b) and L = i "LiW (ab* a*b). 2 (^.105) (4.106) The third independent constant can be chosen to be the difference between i and J^ 2W (aa* bb*) (4.107) It Is not difficult to verify that these constants satisfy the commutation relations of the generators of SU2 or R, i.e., {K. L) = 2D (L, D} = 2K (D, K) = 2L , where (a'^. a} (b*. b) = _ 2i 2W (4.108) (4.109) (4.110) (4.111) It should be noted that in this case the angle variables w^ and w were not used. This occurred only because a similarity was noted between the constants of the motion Jm and Jy and the constants of the motion for the harmonic oscillator. Implicit in the definitions of a and b are two variables, denoted B and r, which are conjugate to Jj^ Jy but these are not the angle variables. These variables are S _ tan1 .^_ In si 2\4x and I r " Â— tan' 1 "^ (4.112) (4.113) 2ff ^/2Wv t is easy to verify that these are conjugate to the action variables
PAGE 103
97 for (r, Jy) Â» {S. J^] = 1 (4.1li*) and (r. 4x) {s. Jv) Â•= . e+.ns) It should not be surprising to find a set of coordinates conjugate to the action variables which are not the angle variables because for a given set of momenta the set of canonical coordinates is not unique. This occurs because one can always add a total time derivative to the Lagrangian which is a function of q, q and t. The constants K, L, and D have remarl
PAGE 104
CHAPTER 5 SUMMARY First of all, the motivation of this dissertation was to find general principles under which symmetry groups of Hamiltonlan systems may be found. One is guided in the search by the i
PAGE 105
99 quadratic operators except in the isotropic case, for which the normal operator is degenerate. The second stage is to consider what Hamiitonians are equivalent to quadratic Hamiitonians by canonical transformations. In this regard one can say almost all are. Here, one reinterprets the HamiltonJacobi equation. If it is thought of as yielding forcefree motion in some (generally noneucl idean) phase space, then a further transformation to polar coordinates will yield circular motion in the phase planes, which presumably implies a quadratic relation between the coordinate and conjugate momentum. This relation is delicate on two accounts. First the relationship in one phase plane may depend upon separation constants arising from other phase planes, and thus irreconcilable conflicts may arise from Inverting all the transformations to the Hami I tonJacobi system. Secondly, with an eye toward quantum mechanical applications, the symmetry operators may well be transcendental functions of the coordinates and momenta, making it impossible to define the proper quantum mechanical functions of noncommuting operators. The organization of the paper was as follows: In the first chapter consideration Is given to the Lie algebra generated by quadratic functions over the phase space under the Poisson bracket operation, and its relation to different homogeneous function spaces also taken over the phase space. It is seen that the eigenfunct ions of these quadratic operators in the phase space itself determine the structure of all the other homogeneous spaces; and In particular there is a simple criterion for finding a set
PAGE 106
100 of operators which cormute with the original operator under Poisson brackets. Namely, it is that they be constructed from products of phase space eigenfunct ions belonging to positive and negative pairs of eigenvalues. These correspond to the "constants of the motion" when the operator Is tal
PAGE 107
101 momenta are found such that the momenta are n1 of the generators of the SUn group which generate rotations in the coordinate planes. The third chapter concerns another specific class of quadratic Hami I tonians , which may be interpreted as belonging to a charged particle moving in a uniform magnetic field as well as a plane harmonic oscillator potential. The problem is of some interest in its own right, and among other things, it is learned how to transfer constants of the motion to a rotating coordinate system. In particular a rule is given for writing Hamiltonians with central potentials when a uniform magnetic field Is added and a transformation is made to a rotating coordinate system with a frequency of rotation equal to the Larmor frequency. The generators of the symmetry group yield a set of intrinsic coordinates, which allow a clean separation of the magnetic field (or angular momentum) from the oscillatory motion; a transformation used by Goshen and Lipl
PAGE 108
102 field and potential field's effects may be seen clearly. A careful treatment is given of the known constants of the motion of the cyclotron problem, known from the work of Johnson and Lippmann [kl]. Particular attention is paid to the influence of the vector potential on the obvious geometric symmetry, and it is seen that canonical rotations and geometrical rotations yield two entirely different tyf>es of symmetry and constants of the motion. Although the specific form of the Hamiltonian has a particularly rich variety of interpretations as a physical problem regarding the motion of a charged particle, mathematically it is essentially nothing but a (possibly anisotropic) harmonic oscillator, a fact which reduces the discovery of the symmetry group to a previous case, and renders the analytic solution of the problem trivial. In general the orbit is a Lissajous figure and the constants of the motion have their interpretation as operators which move the boundaries of the figure or change the phase of the orbit in various ways, preserving the sum of the squares of its semiaxes upon which the energy of the orbit solely depends. In the fourth chapter consideration is given to a nonharmonic problem, the Kepler problem, which is well known to be accidentally degenerate, and equivalent in the bound orbits at least to the harmonic oscillators. As a new result it is found that there is an SU3 symmetry group, in addition to the known R/^ symmetry group, the additional constants defining the line of nodes. In any event the generators of the two groups are functions of oneanother and not subsets. Thus, one has a general functional relation between the generators of these two groups.
PAGE 109
103 As a general rule one can say that any Hamiltonian which is separable will possess an SU^ symmetry even though the frequencies in each angle coordinate may not be equal. Since it is of considerable importance for the further application of this technique, the relation of the SUo group arising from different separations of the Ham i 1 ton Jacob i equation in polar and parabolic coordinates is studied. It is found that considerable attention has to be paid to the problem of finding rational constants of the motion. This work offers many prospects for continuation. There are still many potentials, for which hope of finding rational constants seems slim, but which nevertheless ought to be investigated. Among these are the Stark effect (Redmond recently found a new constant [kS]) , the problem of two fixed attracting centers [46], the motion of a charged particle in the field of a magnetic dipole, as well as the relativistic Kepler problem which Biedenharn has investigated recently [27, 28]. In addition, the relations uncovered between the orthogonal groups and unitary groups might lead to some interesting formulas for the higher dimensional groups, but a point of particular interest is the appearance of the angular momentum operators as a subgroup of SU,. Since they depend on two of the three coordinates of phase space, one can exploit the symmetry of SU^ as a permutation of its generators to obtain SU2 groups involving the remaining coordinate, and hence the radial action. In quantum mechanical terms, this yields another viewpoint toward the known fact that coulomb integrals behave as Wlgner coefficients [47].
PAGE 110
10^ In conclusion it should be pointed out that due to the nonlinearIty of the various relations found, their extension to quantum mechani cal operators will require exceptionally careful scrutiny.
PAGE 111
LIST OF REFERENCES 1. Pauli, W., Z. Physik 36, 336 (1926). 2. Lenz, W. , Z. Physik 24, 197 (1924). 3. Hulthen, L. . Z. Physik 86, 21 (1933). k. Fock, v., Z. Physik 28, 145 (1935). 5. Bargmann, V., Z. Physik g^. 576 (1936). 6. Laporte, 0., Phys. Rev. 50, 400A (1936). 7. Laporte, 0. and Rainich, G.Y., Trans. Am. Math. Soc. 39, 154 (1936). 8. Jauch, J.M., Phys. Rev. 5i, 1132A (1939). 9. Jauch, J.M. , "On Contact Transformations and Group Theory in Quantum Mechanical Problems", Ph.D. Thesis, University of Minnesota (1939). 10. Jauch. J.M. and Hill. E.L., Phys. Rev. 57, 641 (1940). 11. Saenz, A.W. , 'On Integrals of the Motion of the Runge Type in Classical and Quantum Mechanics ", Ph.D. Thesis, University of Michigan (1949). 12. Baker, G.A. , Phys. Rev. J03. 1119 (1956). 13. Demkov, Y.N., J. Exptl. Theoret. Phys. (U.S.S.R.) 26, 757 (195^+). 14. Demkov, Y.N. , Soviet Physics JEPT 263 (1959)Translated from J. Exptl. Theoret. Phys. (U.S.S.R) 36, 88 (1959). 15. Alliluev, A. P., Soviet Physics JEPT 6, I56 (1958). Translated from J. Exptl. Theoret. Phys. (U.S.S.R.) 33, 200 (1957). 16. Loudon. R., Am. J. Phys. 27, 649 (1959). 17. Mcintosh, H.V., Am. J. Phys. 27. 620 (1959). 18. Demkov, Y.N., Soviet Physics JEPT 17. 1349 (1963)Translated from J. Exptl. Theoret. Phys. (U.S. S.R. 7^44, 2007 (1963). 19. Elliott, J. P., Proc. Roy. Soc. (London) A245, 128 (1958). 105
PAGE 112
106 20. Elliott, J. P., Proc. Roy. Soc. (London) A245 . 562 (1958). 21. Bargmann, V. and Moshinsky, M. , Nuclear Phys. j^, 697 (I960). 22. Bargmann, V. and Moshinsky, M. , Nuclear Phys. 23, 177 (1961). 23. Goshen, S. and Lipkin, H.J., Annals of Physics 6, 301 (1959). 2^4. Goshen, S. and Lipkin, H.J., Annals of Physics 6, 310 (1959)25. Biedenharn, L.C., J. Math. Phys. 2, 433 (I960). 26. Moshinsky, M. , Phys. Rev. j26, 1880 (1962). 27. Biedenharn, L.C., Phys. Rev. _[26, 845 (1962). 28. Biedenharn, L.C. and Swamy, N.V.V.J., Phys. Rev. ]^, B1353 (1964). 29. Goldstein, H., "Classical Mechanics ', Addison Wesley, Reading, Mass (1959). 30. Mcintosh, H.V., Journal of Molecular Spectroscopy 8, I69 (1962). 31. Hopf, H., Mathemetishe Anna 1 en ]0k, 637 (1930 , Fundamenta Mathematicae 25^, 427 (1935)The mapping is discussed from a topological point of view. The mappings themselves are presumably of a much older date. 32. Payne, W.T. , Am. J. Phys. 20, 253 (1952). 33. Hill, E.L., "Seminar on the Theory of Quantum Mechanics ', (unpublished). University of Minnesota (1954). 34. Whittaker, E.T. , "Analytical Dynamics , 164, Cambridge University Press (I96I). 35. Born, M., 'The Mechanics of the Atom", Frederick Ungar Publishing Co. , New York (I96O). 36. Wintner, A., 'The Analytical Foundations of Celestial Mechanics", 125, 128, Princeton University Press (1947). 37. Jacobson, N. , "Lie Algebras ', Interscience Publ ishers , New York (1962). 38. Moshinsky, M. , The Three Body Problem and the SU4 Group", (unpublished notes). University of Mexico (1963). 39. Smith, F.T., Phys. Rev. 120, IO58 (I960).
PAGE 113
107 40. Smith, F.T., J. Math. Phys. I. 735 (1962). 41. Mcintosh, H.V., Journal of Molecular Spectroscopy^, 269 (I960). 42. Johnson, M.H. and Lippmann, B.A. , Phys. Rev. J^, 828 (1949). 43. Harrison, E.R., Am. J. phys, 22, 315 (1959). 44. Corben H.C. and Stehle, P., "Classical Mechanics". Second Edition. John Wiley and Sons, Inc., New York (I96O). 45. Redmond, P.J., Phys. Rev. JJ^, B1352 (1964). 46. Landau, L.D. and Lifshltz, E.M. , Mechanics ', 48, Addison Weslev Reading, Mass. (I96O). 47. Melvin, M.A. and Swamy, N.V.V.J., "Journal of Mathematics and Physics 36, 157 (1957).
PAGE 114
BIOGRAPHY Victor August Dulock, Jr. was born February 26, 1939, at Waco, Texas. He was graduated from St. Thomas High School, Houston, Texas in June, 1956. In June, I96O he received the degree of Bachelor of Arts from the University of St. Thomas in Houston, Texas. From September, I96O until the present time he has pursued his work toward the degree of Doctor of Philosophy in the Department of Physics at the University of Florida. During his graduate career, he worked as a graduate assistant and interim instructor in the Department of Physics. Mr. Dulock is married to the former Maryann Welhoelter of Houston, Texas and is the father of three children. He is a member of the American Physical Society. 108
PAGE 115
This dissertation was prepared under the direction of the chairman of the candidate's supervisory committee and has been approved by all members of that committee. it was submitted to the Dean of the College of Arts and Sciences and to the Graduate Council, and vyas approved as partial fulfillment of the requirements for the degree of Doctor of Phi losophy. June 18. 196^ ^^/Dean, College of Arts anw Sciences Dean, Graduate School Supervisory Committee; Cochai rman (xL^y.f/j^^ ;' "Ia :^oJf>i /H.(S <.A^ L^L^ vti.^

