Citation
A Study of accidental degeneracy in Hamiltonian mechanics

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:
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 )
non-fiction ( 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 non-relat i vist ic Kepler problem, whose quantum mechanical analogue is the problem of a one-electron 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:
Thesis--University of Florida, 1964.
Bibliography:
Bibliography: leaves 105-107.
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 non-profit 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 CONTENTS--Continued

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 non-relativistic Kepler problem, whose

quantum mechanical analogue is the problem of a one-electron 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 non-relativistic and

relativistic cases. In the non-relativistic [251 case he has obtained

relations for the representations of R4 exploiting the fact that R4 is

locally isomorphic to R3 x R3. The non-relativistic 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 non-relativistic 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+n-2), 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 f-1. Similarly for i < 0 the dual vector is taken to

be -f .. Thus the dual basis to
-I


(1.8a)


T = (fn, ... fli f-i .... f-n)










is then
r = (-f- ... -f-1 .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 self-dual 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+n-2). 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 = g-i 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
(X-i + Xi) fgi, g-i) = 0, (1.18)


and hence, i. = --i since (g., g-i) 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 = Ci-i gig-i' 9gk
i=l


n
= ) C [g (9g 9, + 9-i 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 g-k' gk)


(1.24)



indices


(1.25)


h =X Ci,-i gig-i.












and hence
C (1.26)
cI-I = g-i gi)


The final conclusion is that



h = g i9-i- (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 gkg-m, (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 g-i' gj} + g-i gi' gj)]
i=l










n

+ Ck,-m [gk9-m' 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 = -'- gig-i (1.30)
(9-i, 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 non-normal 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


C-k,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= gg-i+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) = f-i, (1.37)

or the time reversal mapping T(fi)


T(fi) fi

i > 0 (1.38)
T(f-i) = -f-i


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 +g-i9-i (1.40)
2(9_ gi


n

i=1
P gig; g ;g ;
=1 2(gi, g-i)


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) = g-k



TR(g = (R, gk = -k
g-k


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 g-k k
TR(TT(gk)) = g-k = 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 = (Xe-ip cos T)A/-2 (2.16)


b = (\e-io 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(O-P) 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(o-p) 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 one-headed 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 two-component spinors.












The Anisotropic Harmonic
Oscillator in Two Dimensions


The anisotropic oscillator, both with commensurable frequencies

and with non-commensurable 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 non-commensurable 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*)(w-1)/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 orbit--a Lissajous figure in the x-y plane--would 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 non-algebraic integrals can exist, and as would be expected, are

associated with space-filling motion.

Another point regards the separability of the Hamilton-Jacobi

equation in several coordinate systems. It is pointed out that if a

Hamilton-Jacobi 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 space-filling 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 space-filling 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 2n-1 independent constants of the motion including
2
the Hamiltonian are expected to be found and not n2. Hence there should

exist n2-(2n-1) = (n-1)2 functional relationships among the n2 constants.












As a consequence, the following set of 2n-l constants have been chosen

for the Hopf mapping:

the (n-1) constants Klj where j = 2, ...n,

the (n-1) 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
2-n
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 (n-1)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)e-ij (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 n-I complex dimensions. Regarding

the space as having 2n-2 real dimensions and performing an inverse

stereographic projection a real space of dimension 2n-1 is obtained.

As in the case of the two dimensional oscillator the sphere of

2n-1 dimensions is considered tangent to the 2n-2 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 north-south diameter projected from the 2n-2 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 2n-l coordi-

nates in this space. As in the case for n = 2 these 2n-1 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(n-1) 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


j-1 n

Pj a= a aia (j = 2, ...n), (2.78)
i=l i=j

where P. is the energy difference between the first (j-l) coordinates
J
and the last n-(j-l) 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


j-1


+i.
i=2


j-1
= (%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 J-1 J


(j = 3, ...n).


One might also consider taking the total energy and the energy in the

last n-(j-l) 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, ...n-1)


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, .. .n-1)


Finally, one might consider using the energy and(n-l) 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 k-1 n-k+1 n
Qk In ( (1) ) ( i
2ni j=l aj j=k


k-I

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)[(n-j+l)(p + aO)-(j-1) ok] (2.104)
k=2 k=j

In all three cases, the transformation has been such as to give

the action-angle variables. The effect has been to map the orbit into

a point on the surface of a sphere of 2n-l 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 (n-1) 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 2n-1 independent constants from the n2-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

relation

tan-I 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, ...n-1)


(a J a' + ai J a ))
Kij w.-l ui-l

(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 non-normal 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

-AI-A2 -A5

-A3 A4

L3 -L2

0 Li

-Ll 0

(2.124)


These define a group which is the

dimensional rotation group and the


semi-direct 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 n2-1 linearly independent constants of which

n(n-l)/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 non-zero 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,

(P-e 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
R-l
K [uRv* + u*Rv]/R1/2(uu*)

R-l
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 =N-m 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
R-l '
R1/2 T) 2


uRv* u*Rv
R-l
R1/2 (uu 2


uu* Rvv*

R


uu* Rvv*,(3.32)
,(3.32)

R


= uRv* + u*Rv (3.33)
R-l
Ri/2(uu*) 2


Suu* Rvv*,

R


uRv* + u*Rv
R-L
RI/2(uu ) 2


i uRv* u*Rv (3.34)
R2 R-1(
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

Z-component 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 -(v-v*) my Px (3.41)


Q Sm (v-v*) 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 e-I ] (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 tan-1 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 sin-1 br2-2a2 )+ sin-1 [ 2W2r2-b )
21aI r2[b2_4a2W21/2 2W [b24a2W2]1/2



a br2 -22 w 2W2r2 -b
--- sin-' sinl
21a r[b2-24a2 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




























4-1














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 *







u-0

^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 b2--a2a2
-60 s-i-n-] (n si)- sin-) ( o
21al F2[b2 -4a2 2 2 1/2 2a s 1 2b2-4a 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 semi-major axis of the ellipse lies along

x and the semi-minor 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* P-2 + (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 e-21t) .


In the stationary system the time rate of change is


(3.104)


d (aa*) = (aa*, H) + (aa*)
dt /t


=(1/4)[-21wuxjve21wt + 21 u*ve-21wt


+21iuv*e2iwt 21iu*ve-21tt] = 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 e-i' (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 e-1i/2R e-I~/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)tan-1 (L/K) e-i'l(t-to)
2c 2


+ [(K2 + L2+ D2 1/2-D /2(-i/2)tan- (L/K) ei.2(t-to))

= r + r (3.119)

where
r 1 ()/2 [(K2 + L2 + D2) 1/2D1/2 (-i/2R)tan1 (L/K) e-ll(t-to)
S2c 2
(3. 120)
and

+ = ( R /2[(K2 + L2 + D2)1/2-D]1/2e(-i/2)tan-' (L/K) eIX2(t-t).
+ 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 |r-12
2 [(K2 + L2 + D2)1/2-D [(K2 + L2 + D2) /2+D]


2aR L[(K2 + L2 + D2)1/2-D L[(K2 + L2 + D2)1/2+Di
8c [(K2 + L2 + D2)1/2-D] [(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 +x-1


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 Hamilton-Jacobi 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 semi-major 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 a-b
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
H-D


L -+L' = L(2 + H+D)1/2
H-D


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 Hamilton-Jacobi 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 x-y 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 tan-1 V + sln-1 cos O)/2, (4.53)
p. sin B










w= -dr + f dO + f- dP


= (u E sin u tan-1 V + sin-1 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


tan-1 v = T o "(4.55)



It follows from Fig. 4.1 that sin-1 cot 6 is the projection of the
tan B
angle T onto the plane in which 0 is measured I.e., the x-y plane, which

shall be denoted --. It also follows that sin-1 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

x-axis. 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 x-2 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/2-Lz]/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 Hamilton-Jacobi 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 UTF-8
REPORT xmlns http:www.fcla.edudlsmddaitss xmlns:xsi http:www.w3.org2001XMLSchema-instance xsi:schemaLocation http:www.fcla.edudlsmddaitssdaitssReport.xsd
INGEST IEID EVBD08PTA_2YMF18 INGEST_TIME 2017-07-19T20: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 CONTENTS--Continued 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. 3-1 as viewed from a rotating coordinate system 7' 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

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 non-relat i vist ic Kepler problem, whose quantum mechanical analogue is the problem of a one-electron 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 non-relat i vist ic and relativistic cases. In the non-relat 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 non-relat 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 non-relat 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^srr-s^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), (1-3)

PAGE 15

{f, ag + Bh} = a[f , g} + B(f, h} (linear), {] .k) if. ig. h)) + {g. (h, f}) + {h, {f, g}) = (Jacobi identity). (1-5) 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) (l-8a)

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 self-dual 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(") = ® = (^^, (1-13) 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. g-i) = 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. (1-19)

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-^ rgig-i(1-30) 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 non-normal 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 .C-k.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) . (1-36) On the other hand the mapping P(f:) which simply exchanges coordinates and momenta

PAGE 23

17 P(fi) = f-i. (1.37) or the time reversal mapping T(fj) Kfj) = f, i > (1-38) 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 •< > L-9-k ^ " ° • The composite of any two of the operations on an eigenfunct ion is equivalent to operating with the third operator: M-k ^^ 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. VW--VVV> ('•") 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 (a-p) ^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 one-headed 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 two-component spinors.

PAGE 37

31 The Anisotropic Harmonic Oscillator in Two Dimensions The anisotropic oscillator, both with commensurable frequencies and with non-commensurable 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 non-commensurable 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 x-y plane--would 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 non-algebraic integrals can exist, and as would be expected, are associated with space-filling 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 Hamilton-Jacobi 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 space-filling 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 space-filling 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[a-faj 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 . /-**> (a-fak 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,^ -(a|a, .a.at)5., (2-58) tDj. ^kn,^ = 't(^?k ^-"k)^Tn * <^K ^j<)^-k (ata, a,a;)6,,] . (2-59) 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 2n-l independent constants of the motion including the Hamiltonian are expected to be found and not n2. Hence there should exist n^-(2n-l) = (n-1)^ functional relationships among the n constants.

PAGE 45

39 As a consequence, the following set of 2n-l constants have been chosen for the Hopf mapping; the (n-l) constants K,j where j = 2, ...n, the (n-l) 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 (n-1)^ 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 n-1 complex dimensions. Regarding the space as having 2n-2 real dimensions and performing an inverse stereographic projection a real space of dimension 2n-l is obtained. As in the case of the two dimensional oscillator the sphere of 2n-l dimensions is considered tangent to the 2n-2 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 north-south diameter projected from the 2n-2 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 (j-I) coordinates and the last n-(j-l) 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) j-1 n /^2,^^ r 2 \7 2 . 2 V 2 . 2 1 p. = (\ /2) [cos T+^ €|Sin T-\ eisin tJ 1=2 l-j > 1 = (\^/2) [cos 9 +y e-sin^ 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-(j-l) 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. ...n-1) (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 (2-91) n p. = {\^/2)Y e? sin^ T (j = 2. ...n) (2.92) Q, = p 02 (2.93) Qj = o. Oj.^i (j = 2, ...n-1) (2.9^) Finally, one might consider using the energy anci(n-l) successive energy differences. n Pi = H^ aja; 1=1 n (2.96) P. = a.., a*, -a, at j = (2. ...n) . (2-97) The corresponding coordinates are 1 a)a2 ...a^ Ql = 7-7 ^" 2ni a|a2 -^n k-1 ^* n-k+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)[(n-j + l)(p^-^ a,^)-(J-0^ ^k^ • (2.10A) k=2 kpj In all three cases, the transformation has been such as to give the action-angle variables. The effect has been to map the orbit into a point on the surface of a sphere of 2n-l 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 (n-1) 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 2n-l 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. ...n-l) (2.111) (a; J a; ' •»• a, J aj ') (2.112) Uj-1 Wj-l *^2 /o =.*^ 2 (ajaf)^ (a^ap i[a*'^J a."' -a"J af ' ] •-I J Uj-l U)j-1 (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 non-normal 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. (3-2) 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 non-zero 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) R-1 K = [u^v* + u*'^v]/r'/2(uu*) 2 (3.16) R-1 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. (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 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. (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 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 (3-33) R-r ( R-r 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* (3-35) I {(uu*)>/2 (v*-v). uu*} = (3-36) {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. (3-35) by uu* gives ^^J1, 1^^ . 8U, . (3.38) Expanding Eq. (3-1) 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 Z-component 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 -"^ (v-v*) . muy P^ (3.m) r Q -^ (v-v*) 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>^) (3-48) 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). (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

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 dW--min = ^^(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 = -' (3-52) ., 2 2,1/2 ^^ and D = HHIjlJRvv: ^ (3.53) while the coordinates conjugate to these momenta are Q, = \2['n — + >.i In ^]Ai \,\2 (3-5^) 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 (3-59) 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> (3-65) 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 5-5 \ Jn ^ 2|«| rJ[b2.4aV]'/2 2W [b^-UaVl'^^ (3.82) where 2 « « b = £_+ w2r2 + r2 (3-83) 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. 3-1 through 3-^. I" al ' 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 Wq. Fig. 3-3 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. 3-1 for Pq < 0.

PAGE 70

64 Parameters Initial Conditions u) = 63A 10 « 4 o e = o r o Fig. 3-1 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 (3-91) 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^ (3-96) 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^ e-eo'-— sin-' { , , , ^ ^ } sin-' [ , ° ! . ... ]. 2|a| r^Cb^-^aV]'/^ 2|a| 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 (3-97) 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. 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 Pq which is a linear combination of H and D Pe^ JiL-H ' -^ ^ 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 semi-major axis of the ellipse lies along X and the seml-mlnor axis lies along y. In this case then aa' is pro-2 portional to the maximum value which x attains, because since

PAGE 77

71 o I CD

PAGE 78

72 D

PAGE 79

73 aa* = P-2 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*ve-2''*'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 (3-106) 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. (3-59), (3.65) and (3-66), 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. (3-63) 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. (3-65) and (3-66) It Is apparent that D2)l/2+D]'/2e(-l/2R)tan-l (L/K) e'i^Kt-to) 2c 2 + [(K^ + L^ + D2)'/2.p^]/2^(-i/2)tan-l (L/K) ^i^2(t-to), = '^+ ^ . (3.119) where 2c 2 (3.120) 7_ . 1_ (a)'/2[(K2 + L^ ^ D2)'/^D]'/2e(-'/2R)tan-' (L/K) ^-U,(t-to) and 7^ = i_ (R)'/2[(k2 + L^ + D')'/2.o]'/23(-'/2)tan-l (L/K) J^2(t-to). 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. (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.|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. (3-135) 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 ton-Jacobi 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. (4-2) 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 semi-major 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 H-D L .^L' L(2 + H±D)l/2 H-D (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 ton-Jacobi 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 x-y 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 tan-1 y. + sin"' S2^)/2n ^ [I sin B (4.52) (4.53)

PAGE 95

89 (u e sin u tan-' ^^ + sin"' S£L^ . sin"' ^^i-g + 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°l-2. is the projection of the tan B angle H onto the plane in which is measured i.e., the x-y 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 X-axis. 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.^— ^-i-M. (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), (4-92), (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 jt-sT^ (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 _ tan-1 .^_ 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 force-free motion in some (generally non-eucl 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 ton-Jacobi 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 non-commuting 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 n-1 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 non-harmonic 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; Co-chai rman (xL^y.f/j^^ ;' "Ia :^oJf>-i /H.(S -<.A^ L^-L^ vti.^