arxiv: v4 [math-ph] 3 Nov 2015

Symmetry Integrability and Geometry: Methods and Applications Examples of Complete Solvability of D Classical Superintegrable Systems Yuxuan CHEN Ernie G KALNINS Qiushi LI and Willard MILLER Jr SIGMA 5) 88 5 pages arxiv:5557v4 [math-ph] 3 Nov 5 School of Mathematics University of Minnesota Minneapolis Minnesota 55455 USA E-mail: yc397@camacuk lixx939@umnedu miller@imaumnedu Department of Mathematics University of Waikato Hamilton New Zealand E-mail: math36@waikatoacnz Received May 5 5 in final form October 7 5; Published online November 3 5 http://dxdoiorg/384/sigma588 Abstract Classical maximal) superintegrable systems in n dimensions are Hamiltonian systems with n independent constants of the motion globally defined the maximum number possible They are very special because they can be solved algebraically In this paper we show explicitly mostly through examples of nd order superintegrable systems in dimensions how the trajectories can be determined in detail using rather elementary algebraic geometric and analytic methods applied to the closed quadratic algebra of symmetries of the system without resorting to separation of variables techniques or trying to integrate Hamilton s equations We treat a family of nd order degenerate systems: oscillator analogies on Darboux nonzero constant curvature and flat spaces related to one another via contractions and obeying Kepler s laws Then we treat two nd order nondegenerate systems an analogy of a caged Coulomb problem on the -sphere and its contraction to a Euclidean space caged Coulomb problem In all cases the symmetry algebra structure provides detailed information about the trajectories some of which are rather complicated An interesting example is the occurrence of metronome orbits trajectories confined to an arc rather than a loop which are indicated clearly from the structure equations but might be overlooked using more traditional methods We also treat the Post Winternitz system an example of a classical 4th order superintegrable system that cannot be solved using separation of variables Finally we treat a superintegrable system related to the addition theorem for elliptic functions whose constants of the motion are only rational in the momenta It is a system of special interest because its constants of the motion generate a closed polynomial algebra This paper contains many new results but we have tried to present most of the materials in a fashion that is easily accessible to nonexperts in order to provide entrée to superintegrablity theory Key words: superintegrable systems; classical trajectories Mathematics Subject Classification: C99; C35; E7 Contents Introduction Examples of D nd degree degenerate systems 4 The D4b) oscillator 5 Contraction to the Poincaré upper half plane 9 3 Contraction to the D3 oscillator 3 Embedding D3 in 3D Minkowski space 3 This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour of Sergio Benenti The full collection is available at http://wwwemisde/journals/sigma/benentihtml

Y Chen EG Kalnins Q Li and W Miller Jr 3 Contraction of the D3 oscillator to the isotropic oscillator 4 4 The Higgs oscillator S3 7 4 Contraction of D4b) to the Higgs oscillator on the sphere 4 Contraction of the Higgs to the isotropic oscillator 3 Examples of D nd degree nondegenerate systems 3 3 The system S7 on the sphere 3 3 S7 in polar coordinates 4 3 Trajectories 4 3 System E6 in Euclidean space as a contraction of S7 8 4 The Post Winternitz superintegrable systems 36 4 The classical system 36 5 Classical elliptic superintegrable systems 37 5 Elliptic superintegrability 37 5 Example 4 53 Another version of the example 44 A Review of some basic concepts in Hamiltonian mechanics 46 A Classical integrability 47 A Classical superintegrability 48 Introduction Classical and quantum mechanical Hamiltonian systems that can be solved explicitly both algebraically and analytically and with adjustable parameters are relatively rare and highly prized Famous classical examples are the anharmonic oscillator Lissajous patterns) and Kepler systems planetary orbits) and the Hohmann transfer for orbital navigation [7] Famous quantum examples are the Coulomb system energy levels of the hydrogen atom leading to the periodic table of the elements) and the quantum isotropic oscillator The solvability of these systems is related to their symmetry not necessarily group symmetry This higher order symmetry is captured by the concept of superintegrability A natural classical Hamiltonian system with the Hamiltonian as kinetic energy plus potential energy) on an n-dimensional Riemannian space is said to be maximally) superintegrable if it admits the maximally possible n functionally independent constants of the motion globally defined usually required to be polynomial or at least rational in the momenta) Similarly a quantum Hamiltonian system H = n +V where n is a Laplace Beltrami operator on an n-dimensional Riemannian manifold and V is a potential function) is maximally) superintegrable if it admits n algebraically independent partial differential operators commuting with H There is a rapidly growing literature concerning these systems eg [ 3 4 6 8 9 3 8 3 3 3 33 38 39 4 4] In this paper we consider only classical systems and n = so all of our systems admit 3 independent constants of the motion By taking Poisson brackets of the classical constants of the motion we generate a symmetry algebra not necessarily a Lie algebra which is never abelian This contrasts with integrable Hamiltonian systems which admit n constants of the motion in involution so that the symmetry algebra is always abelian) It is this nonabelian symmetry algebra and its structure that is responsible for the solvability of superintegrable systems To be more explicit along a specific classical trajectory each constant of the motion L j takes a fixed value l j so the trajectory can be characterized as the intersection of n hypersurfaces L j = l j in the n-dimensional phase space Thus in principle the path of the trajectory can be

Examples of Complete Solvability of D Classical Superintegrable Systems 3 determined algebraically though not how it is traced out in time Since the Poisson bracket of two constants of the motion is again a constant of the motion the nonabelian symmetry algebra gives us relationships between the symmetries In the quantum case each bound state eigenspace of the Hamiltonian is invariant under the action of the symmetry algebra so that a knowledge of the irreducible representations of the symmetry algebra gives useful information about the dimensions of the eigenspaces and of the eigenvalues themselves In this paper we illustrate the value of superintegrabilty by studying and solving several families of classical superintegrable systems Second degree superintegrable systems are those whose generating symmetries are all polynomials in the momenta of degree All of these systems have been classified for dimensions as have the systems with nondegenerate potentials in 3 dimensions) [5 9 3] They occur on constant curvature spaces with 3 Killing vectors admitting the most symmetries) on the 4 Darboux spaces admitting Killing vector) and 6 families of Koenigs spaces admitting no Killing vectors) [5 33] Thus after the constant curvature spaces the Darboux spaces admit the most symmetries The Darboux metrics an be written as D: ds = 4x dx + dy ) D: ds = x + x dx + dy ) D3: ds = ex + e x dx + dy ) D4b): ds = An example of a Koenigs space is ds c = x + y + c x + c ) 3 dx y + c 4 + dy ) cos x + b sin x dx + dy ) We first treat some analogs of the harmonic oscillator A treatment of the Kepler system and its analog on the -sphere from this point of view is contained in Chapter 3 of the article [33]) We start by studying a system on the Darboux space D4b) with a -parameter potential Then by taking limits contractions [5 4]) of this system we obtain systems on the Darboux space D3 on the Poincaré upper half plane the -sphere the Higgs oscillator [4]) and finally the isotropic oscillator on Euclidean space The second class of examples are nondegenerate 3-parameter potential) systems The first S7 in our listing [3] can be regarded as a caged version of an analog to the Coulomb potential on the -sphere It contracts to the caged Coulomb system E6 in Euclidean space The preceding examples can also be studied analytically via separation of variables However the Post Winternitz system [36] cannot see also [9] It is 4th degree superintegrable but nonseparable However we show that the classical orbits can be found exactly via superintegrabilty The last examples are different They are separable in elliptic coordinates and can be derived via an action-angle construction The usual action-angle construction of a superintegrable and separable system requires the addition theorem for trigonometric or hyperbolic functions and leads to polynomial superintegrability eg [] The construction here uses the addition formula for elliptic functions [35] It leads to classical systems where some of the constants are nontrivially rational in the canonical momenta We present two examples where however the systems have polynomial symmetry algebras so we consider them as superintegrable and worthy of study The new feature here is that the symmetry algebra closes polynomially even though the system is only rationally superintegrable These are the first examples of such behavior known to us Again we can determine the trajectories exactly With the exception of the elliptic systems all of the superintegrable systems in this paper have been derived and classified before; we take the structure equations as given and show how superintegrability alone leads to formulas for the trajectories These formulas and their analysis are new The elliptic systems have not been found elsewhere to our knowledge so we demonstrate the procedure to derive them

4 Y Chen EG Kalnins Q Li and W Miller Jr Darboux 4 D4b)) oscillator D3 oscillator Higgs oscillator Poincaré upper half plane isotropic oscillator E5 linear potential) Figure A diagram showing the contractions of oscillators This paper is partly pedagogical the Introduction and the Appendix) but mostly new research Sections 5) In the Appendix we give a brief review of fundamental definitions and results from classical Hamiltonian mechanics adapted from [33] needed as background for our computations The point here is to illustrate how using the structure equations of superintegrable systems alone we can derive and classify the trajectories via algebra Hamilton s equations are used only to determine the periods of orbits for systems with degenerate potential Separation of variable techniques are not used except in the last section; our approach is easy to understand geometrically Analogously the spectra of the corresponding quantum systems can be obtained algebraically from the structure relations though we do not treat this here Some recent papers eg [6 7] adopt a related but different approach by using ladder operators constructed from the structure algebra to compute trajectories and spectra for systems with degenerate potential We know of no prior treatments of the nondegenerate caged systems S7 and E6 or of the use of the structure algebra to call attention to special orbits such as those for which R = We point out the contractions that relate our various superintegrable systems Examples of D nd degree degenerate systems Our first examples are degenerate systems These have -parameter potentials and always admit a symmetry that is a st degree polynomial in the momenta hence a group symmetry that can be interpreted as invariance with respect to rotation or translation corresponds to a constant of the motion which leads to an analog of Kepler s nd Law The only possibilities in dimensions are constant curvature and Darboux spaces [] In [33] there is an example of an analog of the Kepler problem on the -sphere that satisfies Kepler s three laws and then by a limiting process a contraction) goes to the Kepler system in Euclidean space Here we will treat an analog of the harmonic oscillator on the Darboux space D4b) and show that it obeys analogs of Kepler s laws Then by taking contractions to superintegrable systems on the Darboux space D3 on the Poincaré upper half plane equivalent to a system on a hyperboloid in Euclidean space) and finally to the isotropic oscillator system in Euclidean space we will see that using ideas from superintegrability theory alone we can understand the basic properties of these systems: conservations laws explicit trajectories etc and how they are related Fig describes the contraction relationships that we will exploit An important feature of degenerate systems is that they always admit 4 linearly independent symmetries that are nd degree in the momenta whereas only 3 can be functionally independent Thus there must be a relation between these symmetries This relation emerges from the structure algebra obeyed by the symmetries

Examples of Complete Solvability of D Classical Superintegrable Systems 5 The D4b) oscillator We consider the superintegrable system H = sinh x) p coshx) + b x + p α y) + coshx) + b ) with a basis of constants of the motion H J = p y and sinh ) x) Y = cosy) p coshx) + b x + p ) y coshx)p y α cosy) siny) sinhx)p x p y 4 cosh x) + b ) sinh ) x) Y = siny) p coshx) + b x + p ) y coshx)p y α siny) + cosy) sinhx)p x p y 4 cosh x) + b 3) Here the spaces D4 are indexed by a parameter b > In the limit as b the space becomes a -sphere and this system becomes the Higgs oscillator [4]) The variable y can be interpreted as an angle and the space and potential are periodic in y with period π The constants of the motion generate an algebra under the Poisson bracket obeying the structure equations {J H} = {Y H} = {Y H} = 4) {J Y } = Y {J Y } = Y {Y Y } = 4J 3 + bj H αj 5) with Casimir Y + Y = J 4 + H + bj H αj 6) Note that there are 4 linearly independent symmetries J H Y Y as polynomials in the momenta but there can be only 3 functionally independent generators The dependence relation is given by 6) From the first two equations 5) we see that Y Y ) transforms like a -vector under rotations about the 3-axis The sum Y + Y is expressed in terms of constants of motion so the sum is constant We set Y + Y κ where κ is the length of the -vector Thus we can choose a preferred coordinate system such that Y = κ and Y = To determine a trajectory we need to express x and y in terms of the constants of the motion along the trajectory We do this by eliminating the momenta from equations ) 6) We see that Y and Y can be expressed in an alternate form: Y = cosy)h + cosy) coshx)p y siny) sinhx)p x p y 7) Y = siny)h + siny) coshx)p y + cosy) sinhx)p x p y 8) To eliminate p x we multiply equation 7) by cosy) equation 8) by siny) respectively and add them together which gives us Y cosy) + Y siny) = H + coshx)p y Using the facts that Y = κ Y = and J = p y we get the orbit equation coshx)j cosy)κ = H 9) Since coshx) is an even function in x there will be two orbits for positive x and negative x that are symmetric with respect to x-axis Therefore without loss of generality we restrict ourselves to positive x

6 Y Chen EG Kalnins Q Li and W Miller Jr Analog of Kepler s second law of planetary motion We exploit the fact that there is a st degree constant of the motion We consider an orbital motion as taking place in the s s plane with polar coordinates s = r cos θ s = r sin θ where r = coshx) + b)/ sinh x) and θ = y In these coordinates the metric is ds = + br + 6r 4 + 4br + + br ) dr + r dθ If α > the potential is attractive to the origin in the plane; if α < the potential is repulsive and the trajectories are unbounded Theorem Trajectories of the D4b) oscillator sweep out equal areas in equal times with respect to the origin in the s s plane Proof In the interval from some initial time to time t the area swept out by the segment of the straight line connecting the origin and the object is At) = θt) θ) r θ)dθ Thus the rate at which the area is swept out is da dt = da dθ dθ dt = r θt)) dθ dt = r dy dt ) since dy/dt = dy/dθ) dθ/dt) = /) dθ/dt) To get dy/dt) we take the Poisson bracket of H and y which is {H y} = dy dt = p y sinh x) coshx) + b = J sinh x) coshx) + b then plug in this expression for dy/dt) into equation ) to get da/dt) = J / which is a constant Analog of Kepler s third law of planetary motion Theorem For α > the period T of an orbit is ) T = π + b b + ) H bh + α) H bh + α) Proof The total area swept out as the trajectory goes through one period is AT ) = π r θ)dθ ) = 6 J π + b) κ + H + J H + J 4 + b) κ + H J H + J 4 κ + H J H + J 4 ) κ + H + J H + J 4 ) However from the second law we see that AT ) = T/) J and using the fact that κ = J 4 + H + bj H αj the period can be expressed in the form ) Now we begin an analysis of the trajectories We first assume that α > so that the potential is attractive Restriction on H Recall there is a restriction when we try to express r in terms of θ and other constants of the motion That is coshx) = H + cos θκ)/j ) should always be larger than for any θ Then since κ > we have the restriction H κ > J implying that for a closed trajectory H should always be positive Here κ is the length of the -vector Y Y ) To be explicit κ = J 4 + H + bj H αj Then by squaring both sides of an alternate

Examples of Complete Solvability of D Classical Superintegrable Systems 7 S vs S 5 5 - -5 - -5 5-5 - -5 Figure Orbit plot for J = H = 3 κ = b = equal axes) form of the restriction equation H J > κ and plugging in the expression for κ we get α b + )H)J > Noticing the fact that b + H and J are all larger than we reach the conclusion H < α/b + ) for a bounded orbit For larger values of H the trajectory is unbounded Case of bounded trajectory To plot the trajectories on the s s plane it is convenient to write s and s both in terms of θ which is s = J H + cos θ κ + bj ) H + cos θ κ) J 4 cos θ s = J H + cos θ κ + bj ) H + cos θ κ) J 4 sin θ 3) Note: it is useful to divide the r part of s and s by J 4 and introduce new constants H = H/J and κ = κ/j so we can eliminate one constant) Since in s and s J always appears in the form of J so for every positive J there will always be a duplicate case for the corresponding J and without loss of generality we can restrict our discussion to J > A typical plot of the trajectory on the s -s plane is Fig and as one or more of H J and b becomes larger and larger they will dominate the r term and make r less susceptible to the change of cos θ Thus the plots will look more and more circular Furthermore if κ becomes very large the plots will tend to move towards the negative s direction and appear in an elongated form as in Fig 3 The plot will always be symmetric with respect to s -axis because s θ) is even and s θ) is odd in θ Also since the plot will always lean towards the negative s side larger r always occurs at smaller s value To trace out an orbit in time we would need a starting point to integrate Hamilton s equations Conventionally we choose a point that is closest to the origin which we call perihelion From the plots of the orbits and an easy analysis of equations 3) it is obvious that this point is the intersection of the trajectory with the positive s axis The point on the x-y plane that corresponds to this perihelion is x ) where x = κ+h arccosh ) Plugging y = into J equation 8) gives p x p y = Realizing p y = J we know p x must be equal to We see that the perihelion points on the phase space trajectory are uniquely determined by the constants of the motion as follows: coshx) = H + κ J y = p x = p y = J 4)

8 Y Chen EG Kalnins Q Li and W Miller Jr S vs S 4 3-6 -5-4 -3 - - - - -3-4 Figure 3 Orbit plot for J = H = κ = b = axes are not set equal) S vs S s - -5-5 5 s Figure 4 s vs s for J = H = κ = b = and we can see that if we know the perihelion point on the phase space as in equation 4) a set of constants of the motion can be uniquely determined Case of escape velocity This is the case when H = J + κ a bifurcation point on the momentum map In polar coordinates as θ π we have r whereas for s s in equations 3) s and s A plot for s vs s is given in Fig 4 In Fig 4 the two tails will extend to infinity in the direction of negative s -axis and they will get closer and closer to the s -axis for smaller and smaller s value but they will never touch the s -axis Case of unbounded trajectory Here there is more than one value of y ie θ/ that has no corresponding real value of x in equation 9) Example 3 H J ) In this case since H J r will blow up before s goes becomes negative so the entire trajectory will be bounded on the positive-s side of the plane A boundary case H = J ) is plotted in Fig 5 For Fig 5 the constants are J = H = κ = b = so for θ < π/ and 3π/ < θ π corresponding real values of x exist and there is a trajectory For θ = π/ and 3π/ r will goes to infinity and at the same time s will go to infinity and s to which is represented by the two tails in Fig 5 Example 4 J H < J + κ) In this case the trajectory will extend to the negative-s side of the plane A typical plot for an orbit is given in Fig 6

Examples of Complete Solvability of D Classical Superintegrable Systems 9 5 S vs S 4 3 4 6 8 - - -3-4 -5 Figure 5 s vs s for J = H = κ = b = S vs S 5 4 3-3 -5 - -5 - -5 5 - - -3-4 -5 Figure 6 s vs s for J = H = κ = b = Contraction to the Poincaré upper half plane Using the Hamiltonian H ) modified by a constant as H = H [α/ + b)] we let x = ɛy y = ɛx α = [ + b)β]/[ɛ ] J = K/ɛ Y J )/ɛ ) + H X Y X )/ɛ) and go to the limit as ɛ Then H = 4Y p b + X + p ) Y + βy K = p X X = X Y ) p X XY p X p Y βx 4 b + ) X = Xp X + Y p X p Y + The structure equations become βb + ) X 4 {K X } = X {K X } = K + β 4 b + ) {X X } = KX + b + KH and the Casimir is βb + ) X + + K )X b + 4 4 K H = 5)

Y Chen EG Kalnins Q Li and W Miller Jr Circle 8 6 Y 4 Origin 4 6 8 4 6 8 X Figure 7 Circular trajectory on Poincaré upper half plane with β = b = H = K = X = If the potential is turned off this is essentially the Poincaré upper half plane model of hyperbolic geometry; for b = 4 it is exactly that The space here consists of all real points X Y ) with Y > In this limit there is no longer any periodicity Using these identities to eliminate the momenta we arrive at the equation for the trajectories: 4K ) + βb + ) X X + K 4K + βb + ) ) Y H K b + ) = 6) Notice that the orbit equation is even in Y so the trajectory will always be symmetric with respect to the X-axis When restricting to the Poincaré upper half plane with Y > the trajectory will be the upper half of the trajectory traced out by the orbit equation If β > so that the potential is attractive to the boundary Y = this describes the portion of the ellipse x /A ) + y /B ) = in the upper half plane where A = 4H K b + ) 4K + βb + )) B = x = X 4X 4K + βb + ) y = Y H b + ) 4K + βb + ) If β the potential is repulsive If however 4K + b + )β > the trajectory 6) is again a portion of an ellipse There is a special case that when β = we will have A = B such that the trajectory will be upper half of a circle A plot of a trajectory when β = is given in Fig 7 And an example of elliptic trajectory is given in Fig 8 If 4K + b + )β < the trajectory is the upper sheet y > of the hyperbola x /A ) y /B ) = where A = 4H K b + ) 4K + βb + )) B H b + ) = 4K + b + )β 4X x = X 4K + βb + ) y = Y If 4K + b + )β = equation 6) becomes degenerate From 5) we now find X = Y K XY p Y K and have H = [4Y /b+)] p Y Eliminating p Y from these two expressions we find the equation K Y +X ) H K b+)x = for the unbounded trajectories These are upper halves of parabolas with orbit equation Y = X /K) ± H b + )X and we notice that a certain subset of H X K and b would correspond to two half-parabolas which are symmetric with respect to y-axis Interestingly when X /K the two half-parabolas will

Examples of Complete Solvability of D Classical Superintegrable Systems 8 Ellipse 7 6 5 Y 4 3 Origin 3 4 5 6 7 8 9 X Figure 8 Elliptic trajectory on Poincaré upper half plane with β = b = H = K = X = 6 Hyperbola 5 4 Y 3 - -8-6 -4 - Origin 4 6 X Figure 9 Hyperbolic trajectory on Poincaré upper half plane with β = 3 b = H = K = X = 6 A Parabolic Trajectory Y 4-5 - -5 Origin 5 5 X Figure Parabolic trajectory on Poincaré upper half plane with β = b = H = K = X = Y A Intersecting Parabolic Trajectory 7 6 5 4 3-5 - -5 Origin 5 5 X Figure Parabolic trajectory on Poincaré upper half plane with β = b = H = K = X = cross as in Fig A special case occurs when H = so that the trajectory becomes a straight line Since H = 4Y p Y )/b + ) we see this implies that p Y = so the trajectory must be parallel to the x-axis Furthermore if H = and X /K = the trajectory will be exactly the y-axis Moreover for H = and X /K < there will be no trajectory

Y Chen EG Kalnins Q Li and W Miller Jr We mention that this Poincaré upper half plane model can be mapped to the unit disk the Poincaré disk as well as to the upper sheet of the -sheet hyperboloid is 3-space Further this model can be contracted to the flat space system H = p x + p y + αx E5 in our listing [3] Indeed if we set X = y Y = x + b + ɛ β = α b + Hɛ) = ɛ H ɛ ) αb + )3/ 4ɛ then lim ɛ Hɛ) = p x +p y +αx We don t go into detail here but we give another detailed example in the next section 3 Contraction to the D3 oscillator The oscillator in the space Darboux 3 D 3 ) has Hamiltonian H = e x p e x + x + p β y) + e x + 7) with a single Killing vector J = p y and a symmetry algebra basis {H J X X } Here X = ex sin yp x p y + 4 e x + cos yp x e x e x + ) 4 e x + X = ex cos yp x p y + e x 4 e x + sin yp x e x e x + ) 4 e x + The structure relations are e x cos yp y + β cos y e x + sin yp y + {J X } = X {J X } = X {X X } = J H β J β sin y e x + and there is the functional relation 4X + X 4 H J H + β J = Theorem 5 The D3 oscillator system is a contraction of the D4 oscillator system Proof We can get 7) as a limiting case of ) by taking x = x y lnɛ) y = b = ɛ H = 8ɛH α = β 8ɛ 8) e x Then we have H = e x + p x + p y ) J = J = p y and Y 4ɛY Y 4ɛY First it is obvious that J = J = p y Then as ɛ since coshx) = ex lnɛ) + e x +lnɛ) = ex + ɛ e x ɛ sinhx) = ex ɛ e x ɛ plugging these two identities into the Hamiltonian equation ) and also using equation 8) we have H = 8ɛ ex ) ɛ +ɛ 4 e x 4ɛ 4 p x + β 4 p y ) + 8ɛ 8ɛ = e x +ɛ e x ɛ + ɛ e x +ɛ e x ɛ e x ɛ + ɛ 4 e x p e x + ɛ e x + x + p ) β y + e x + ɛ e x + = + ɛ e x e x + p x + p y ) + β e x +

Examples of Complete Solvability of D Classical Superintegrable Systems 3 because the terms with ɛ in front of them will all go to as ɛ This H is exactly the Hamiltonian of the Darboux-3 oscillator Similarly by substituting the identities in equation 8) into equations 7) and 8) we have Y = 4ɛ cos y ) ex ɛ + ɛ 4 e x 4ɛe x + ɛ e x + ) 4 p x + ) 4 p y cos y β 8ɛe x + ɛ e x + ) ) + cos y e x + ɛ e x ɛ e x e x + ) 4 e x + e x = 4 e x + cosy )p x 4 as ɛ In the same sense e x Y = 4 e x + siny )p x 4 4 p y sin y ) ex e x + e x e x + ) e x + e x e x + ) e x + ) e x ɛ e x ɛ 4 p x p y = e x 4 e x + cosy )p x cosy )p y + ex siny )p x p y + β cos y e x + cosy ) p y + ex siny )p x p y + β cos y e x + siny ) p y ex cosy )p x p y + β sin y e x + We notice they are exactly the same vectors as X and X So we have proved that in this limiting case Darboux-4 space contracts to Darboux-3 Using this transformation we can get the orbit equation for the D3 oscillator from the orbit equation for the D4b) oscillator First since κ = Y + Y = 6ɛ Y + Y ) = 6ɛ κ and assuming ɛ > we have κ = 4ɛκ Then putting κ = 4ɛ κ and the identities 8) into 9) we derive e x +ɛ e x ɛ ) J ) cosy ) κ 4ɛ ) = H 8ɛ Thus ex J cosy )κ = H as ɛ 3 Embedding D3 in 3D Minkowski space We can embed D3 as a surface in 3D Minkowski space with coordinates X Y Z in such a way as to preserve rotational symmetry For example let X = e x + e x cos y Y = e x + e x sin y 6 3 Z = ln 6 6 + 5ex ) ) 3 + e x + 5e x + 3 6 6 + 5ex ) e x 3 + e x + 5e x 3 + e x + 5e x Then dx + dy dz = ex +) dx + dy ) Such embeddings are not unique e x Discussion of trajectories To see how the trajectory behaves on this surface we impose the orbit equation e x J cosy )κ = H to X Y and Z From the orbit equation we have e x = H + cos y κ J = H J + κ cos y = a + b cos y 9) J where we ignore all the primes on the letters for convenience Also we treat a = H and b = κ J J as two new constants so the number of constants is reduced by Then we can express X Y and Z all in terms of y We notice that since e x > the trajectory will be closed when a > b and unbounded when a < b Example 6 Plots of an elliptical shape trajectory and its overhead -D view are given in Figs and 3 The size of the cone that appears on the graphics is determined by the minimum value we choose for e x in MATLAB The smaller the minimum value the larger the cone For these two plots the minimum value for e x we choose is and we adjust to different sizes of cones for other examples to make the plots clearer)

4 Y Chen EG Kalnins Q Li and W Miller Jr Figure 3-D View of trajectory with a = b = Figure 3 Overhead view of trajectory with a = b = Example 7 As a becomes smaller the ellipse is prolonged see Figs 4 and 5 Seen from above the shape of the trajectory two merging ellipses The minimum value for e x we choose is for these two plots) Example 8 escape velocity) For a = b we encounter the boundary case between bounded trajectory and unbounded trajectory A plot of the case is given in Fig 6 The shape of the trajectory when looking from above is like two very thin parabolas meeting at the top of the cone The minimum value for e x we choose is for these two plots) Example 9 unbounded trajectory) When a < b there are no real x values corresponding to a range of values of y clearly seen from equation 9) so the trajectory is unbounded Examples are given in Figs 7 and 8 which are a standard 3-D view and a closeup overhead view 3 Contraction of the D3 oscillator to the isotropic oscillator The isotropic oscillator in Euclidean space has the Hamiltonian H = p x + p y + ω x + y ) )

Examples of Complete Solvability of D Classical Superintegrable Systems 5 Figure 4 3-D View of trajectory with a = b = Figure 5 Overhead view of trajectory with a = b = Figure 6 3-D View of escape velocity trajectory with a = b = and basis symmetries K = xp y yp x L = p x p ) y + ω x y ) L = p x p y + ω xy K The structure equations are {L K} = L {L K} = L {L L } = ω K

6 Y Chen EG Kalnins Q Li and W Miller Jr Figure 7 3-D View of trajectory with a = 7 b = Figure 8 Closeup overhead view of trajectory with a = 7 b = The functional relation is L + L = 4 H ω K We can obtain this system as a limit of the D3 oscillator as follows: In 7) we set x = x + ln ɛ ) y = y H = ɛh β = ω )/ɛ Then as ɛ we have K = p y and H = e x p x + p 4ω y ) + e x X = ɛx = sin ex y p x p y + cos y p x ) ) py + ω cos y e x X = ɛx = ex cos y p x p y + sin y p x py ) ) + ω sin y e x In terms of flat space Cartesian coordinates X = r cos θ Y = r sin θ we have e x = 4 r y = θ H = p X + p Y + ω X + Y ) X = 4 p X p Y ) + 4 ω X Y ) X = p Xp Y + ω XY

Examples of Complete Solvability of D Classical Superintegrable Systems 7 Figure 9 Elliptical trajectories with K = Xp Y Y p X ) From the expressions of H and X we get p x = H + 4X ω X p y = H 4X ω Y ) Then from the expression of X we have p xp y = X ω XY Squaring this equation and equating it to the product of the two equations ) we obtain the orbit equation ω X + Y ) H 4X ) 6ω X XY = H 6X 6X The shape of the orbit depends only on 8X /H 4X ) so we are really investigating the equation X + Y axy = where a = 8X /H 4X ) We take the right side to be since we now only care about the shape not the scale When a = ie X = the trajectory is just a circle When < a < ie < 8X /H 4X ) < the trajectory is a tilted ellipse with y = x as the major axis In Fig 9 we show three trajectories for different values of a and also the major axis As a increases the ellipse becomes more elongated in the major axis direction When a = ie H 4X = 4X there is a bifurcation point on the momentum map The trajectory splits into two parallel straight lines For a > ie 8X /H 4X ) > the trajectories are hyperbolas with y = x as the symmetry axis In Fig we present three hyperbolic trajectories shown in darker colors as a increases It is also possible that 8X /H 4X ) < in which case the trajectories will be symmetric about the line y = x 4 The Higgs oscillator S3 The classical system S3 on the -sphere is determined by Hamiltonian [4] H = J + J + J3 + α s + s + ) s 3 ) s 3 where J = s p 3 s 3 p and J J 3 are cyclic permutations of this expression For computational convenience we have embedded the -sphere in Euclidean 3-space Thus we can write H = p + p + p 3 + α s 3 = H + s p + s p + s 3 p 3 ) s + s + s 3

8 Y Chen EG Kalnins Q Li and W Miller Jr Figure Hyperbolic trajectories and use the Euclidean space Poisson bracket {F G} = 3 i= si F pi G + pi F si G) for our computations but at the end we restrict to the unit sphere: s + s + s 3 = and s p + s p + s 3 p 3 = The Hamilton equations for the trajectories s j t) p j t) in phase space are ds j dt = {H s j} dp j dt = {H p j} j = 3 The classical basis for the constants of the motion is L = J + α s s L = J J α s s 3 s X = J 3 3) 3 The structure relations are {X L } = L {X L } = L H + X + α {L L } = L + α)x 4) and the Casimir relation is L + L L H + L X + αx + αl = To analyze the classical trajectories it is convenient to replace the basis elements L L in the algebra with the new basis set S = J J α s )) s s S = J J αs ) s 3 s 5) 3 Thus S = L + X H + α) S = L Now the first two equations 4) become {J 3 S } = S {J 3 S } = S so S S ) transforms as a -vector with respect to rotations about the 3-axis Indeed a rotation through the angle β about the 3-axis rotates the vector by β The remaining relations now become {S S } = X H X + α ) S + S = 4 X H + α ) αx κ

Examples of Complete Solvability of D Classical Superintegrable Systems 9 This verifies that the length κ of the -vector is unchanged under a rotation and shows that S S ) is similar to the Laplace Runge Lenz vector for the Kepler problem in Euclidean space However a better superintegrable -sphere analog of the Kepler problem is the potential αs 3 / s + s called S6 in our listing [3]) In analogy with the choice of periaptic coordinates to simplify the Kepler problem we choose the preferred coordinate system such that the vector S S ) points along the -axis ie S = S = κ Then equations ) and 3) become J J = αs s s 3 J J = κ + α s ) s s 3 J3 = X J + J + J3 = H α s 6) 3 The last 3 equations can be solved to give J = κ + H X α ) αs J = κ + s 3 H X α ) αs s J3 = X 3 Substituting these results into the square of the first equation 6) and simplifying we get the result assuming α ) [ κ + H X ] [ α s + κ + H X ] α s X s 3 = 7) This is the equation of a cone As + Bs + Cs 3 = The trajectories lie on the intersection of this cone and the unit sphere s + s + s 3 = Thus we get conic sections again as in the Euclidean Kepler problem but this time the sections are intersections with the unit sphere rather than planes The possible types of trajectory will depend on the signs of A B C The projection on the s s plane is the curve H + X ) α H + X + κ s ) α + κ s = X 8) The identities [ κ + ] [ H X α) κ + ] H X α) [ κ + ] [ H + X α) κ + ] H + X α) [ κ + ] [ H X + α) κ + ] H X + α) = αx 9) = HX 3) = αh 3) will be important in the analysis of trajectories to follow Analog of Kepler s second law of planetary motion We see from equations 8) and 3) that for the nonzero angular momentum X the projection of the motion on the unit circle in the s s plane is a segment of an ellipse hyperbola or straight line and that none of these curves pass through the center s s ) = ) of the circle Thus as the particle moves along its trajectory s t) s t) s 3 t)) the line segment connecting the projection s t) s t)) to the center of the circle sweeps out an area Introducing polar coordinates s t) = rφt)) cos φt) s t) = rφt)) sin φt) we see that in the interval from some initial time to time t the area

Y Chen EG Kalnins Q Li and W Miller Jr φt) swept out is At) = φ) r φ)dφ Thus the rate at which the area is swept out is da da dφ dφ dt = r φt)) dφ dt Now note from Hamilton s equations that along the trajectory ds X = J 3 = s t)p t) s t)p t) = s dt s ds dt Thus da dt = X a constant = r dφ dt dt = Theorem If the angular momentum is nonzero the projection of the trajectory on the unit circle in the s s plane sweeps out equal areas in equal times Now we begin an analysis of the trajectories There are two cases depending on whether the potential is repulsive α > ) or attractive α < ) Case : α > This is the case of a repulsive potential The equator of the sphere repels the particle Thus motion is confined to a hemisphere We start an analysis of the types of trajectories If X = and using the fact that all trajectories will be periodic for a repulsive force we see that there will necessarily be at least one point on the phase space trajectory for which p = Substituting into the phase space conditions 6) it is easy to see that this is possible only for s = Solving all the equations completely we find that these points on the phase space trajectory are uniquely determined by the constants of the motion as follows: H + X α ) s = s = κ + p = p 3 = H p = κ + H + X α ) 3) This prescription gives us a point on each trajectory comparable to the aphelion for the Kepler system where we can start to trace out the orbit It remains to analyze the possible orbits To see what is the available parameter space for the case X = note that first we must require κ Noting the identity 9) we see that the quantities in brackets must have the same sign However if that sign is negative then the cone 7) will degenerate to a point and not intersect the sphere Thus to obtain trajectories it is necessary that the constants of the motion satisfy κ + H X α) > κ On the other hand if these conditions are satisfied we see from equations 3) that there exist trajectories for which the corresponding constants of the motion are assumed Thus the conditions are necessary and sufficient An Analog of Kepler s 3rd law of planetary motion For nonzero angular momentum X the projection of the motion on the unit circle in the s s plane is an ellipse 8) enclosing the center of the circle The area of this ellipse is easily seen to be πx / H Let T be the period of the trajectory Thus T is the length of time for the projection of the trajectory to trace out the complete ellipse and AT ) = πx / H Since the area of the ellipse is swept out at the constant rate da dt = X da we have AT ) = dt T = X T Equating these two expressions for At) and solving for T we find T = π/ H Theorem For X = the period of an orbit is T = π/ H If X = then from 7) and 8) we see that the projection of the motion in the s s plane is the straight line segment s = Thus this motion takes place in a plane and the trajectory is a portion of a great circle passing through the poles of the sphere Here κ = H α) H α > The motion is periodic and the points of closest approach to the equator are those such that s = H α)/h The period is again T = π/ H There is a special case of equilibrium at a pole when H = α Case : α < This is the case of an attractive potential where the equator of the sphere attracts the particle Again motion is confined to a hemisphere We first consider the case

Examples of Complete Solvability of D Classical Superintegrable Systems where X Without loss of generality we can assume X > Then the projection of the motion is counter-clockwise Now κ > and from the identity 3) we see that there will be 3 classes of trajectories depending on the value of H If H > then the coefficients of s and s in the projection formula 8) must have the same sign necessarily positive for a real trajectory Thus H > : H + X α > κ > H X α κ + H X α > and the projection will be a segment of an ellipse It is straight-forward to check that the ellipse always intersects the circle at 4 points All of these trajectories lead to annihilation at the equator There will necessarily be a perihelion point on each trajectory: s = p = p 3 = s = X / [ κ + H + X α) ] p = κ + H + X α) If H < then the coefficient of s is positive and the coefficient of s is negative in the projection formula 8) Thus H < : H + X α κ < H + X α + κ > and the projections will be segments of hyperbolas Each branch of the hyperbola intersects the circle at exactly points Again these trajectories lead to annihilation at the equator Again there is a perihelion point on each trajectory: s = p = p 3 = s = X / [ κ+ H+X α) ] p = κ + H + X α) If H = then κ = X α) so the coefficient of s is positive and the coefficient of s is in the projection formula 8) Thus H = : X α = κ and the projections will be segments of straight lines s = X /X α) Here p = and p = X α a constant in accordance with the analog of Kepler s Second Law These trajectories lead to annihilation at the equator There is a perihelion point on each trajectory: s = p = p 3 = s = X /X α) p = X α Finally we suppose X = so that the motion takes place in a plane through the poles If H α > then the motion takes place in the plane s = All trajectories annihilate at the equator The projected trajectories each pass through the center of the circle and we have p = H α Hs If H α < then the motion takes place in the plane s = All trajectories annihilate at the equator The projected trajectories do not pass through the center of the circle and we have p = H α Hs If H = α then there are possible trajectories on any plane passing through the poles If we go to new rotated s s coordinates such that the motion takes place in the plane s = then p = and p = αs All trajectories annihilate at the equator except for unstable equilibria at the poles 4 Contraction of D4b) to the Higgs oscillator on the sphere By letting ɛ in b = + ɛ and defining new constant α = 4β the Hamiltonian equation ) becomes H = sinh x) p coshx) x + p 4β y) + coshx) = cosh x p x + p ) βcosh x sinh x) y + sinh x We can ignore the constant term β and write H = cosh x p x + p ) β cosh x y + sinh x sinh x cosh x) = p sinh x) x + p ) 4β y + sinh x) = cosh x p x + p y such that H = H + β This is the Higgs oscillator on the -sphere ) + β cosh x sinh x β

Y Chen EG Kalnins Q Li and W Miller Jr In terms of the Euclidean embedding of the sphere we have s =cos y/cosh x s =sin y/cosh x s 3 = sinh x/cosh x so s +s +s 3 = If we set s = r cos θ s = r sin θ as in polar coordinates then r = /cosh x and θ = y However we notice that in the analog of Kepler s nd law part of the Section r = coshx) + b/ sinh x) and θ = y As ɛ r = / cosh x) Thus r = r and θ = θ/ so this contraction space is a double covering of the original D4b) space We should be careful of this fact in the further computations Expressed in the phase space s s s 3 p p p 3 ) we have H = J + J + J 3 + β where J s = s p 3 s 3 p and J J 3 3 are cyclic permutations of this expression Considering the potential part of the Hamiltonian we have a new basis set Y = J J βs s )/s 3 Y = J J βs s /s 3 ) J = J 3 We will show that in this limit the orbit equation 9) will yield the orbit equation 8) First we transform coshx) and cosy) in equation 9) to express them in terms of s and s s 3 can be expressed in terms of s and s ) From the fact that cosh x = /s + s ) we have coshx) = cosh x = s + s cosy) = cos y sin y = s s s + s Then equation 9) becomes [/s + s ) ]J [s s )/s + s )]κ = H where κ is the length of the two vector Y Y ) and we set Y = such that κ = Y Multiplying the equation by s + s ) plugging in the fact that H = H β and doing some rearrangements we obtain ) ) H + J β + κ s H + J β κ + s = J 33) where we realize that κ = Y which is times the basis S in equation 5) so κ = κ where κ is the one in equation 8) Therefore we can see that 33) is the same orbit equation as 8) Analog of Kepler s second law of planetary motion Following the same procedure as in Section and noticing that r = r and θ = θ/ we have da dt = 4r θt)) ) dθ dt = r dθ dt = da da dt Therefore dt = J which matches the expression for the Higgs oscillator Period of the orbit With the same procedure as in Section and r = r and θ = θ/ we find area swept out as a trajectory goes through one period T is A = π r θ )dθ = π r θ)dθ = 4A where A is the area in equation ) Also from the second law we see that A = J T Hence T = A J = 4A J = T = π + b H bh + α) + ) b H bh + α) where T is the same as equation ) In the limit as ɛ we take b = into the above equation of T and get T 4 = π = 4H + α π 4H β) + 4β = π H which differs from the period in Theorem by a factor of due to the fact that the contraction is a double covering of the original D4 space 4 Contraction of the Higgs to the isotropic oscillator This contraction has a simple geometric interpretation the contraction of a -sphere to a plane We can consider the Higgs oscillator as living in a -dimensional bounded universe of radius

Examples of Complete Solvability of D Classical Superintegrable Systems 3 in some set of units Suppose an observer is situated in this universe near the attractive north pole Our observer uses a system of units with unit length ɛ where < ɛ and we suppose that using these units the universe appears flat to the observer Thus in the observer s units we have s = ɛx s = ɛy s 3 = ɛ X + Y ) = X + Y )ɛ / + Oɛ 4 ) Here ɛ is so small that to the observer it appears that s 3 = Thus to the observer it appears that the universe is the plane s 3 = with local Cartesian coordinates X Y ) We compare the actual system on the -sphere with the system as it appears to the observer and we assume that ɛ is so small that it can be neglected unless we are dividing by it Now we have s = ɛx s = ɛy s 3 p = p X ɛ p = p Y ɛ p 3 Xp X + Y p Y ) We define new constants ω h by α = ω /ɛ 4 E X ω /ɛ 4 = h/ɛ where J 3 = X Then substituting these results into the Higgs Hamiltonian equation H J + J + J 3 + α/s 3 = E we find p X + p Y )/ɛ + ω / ɛ X + Y ) = h/ɛ Multiplying both sides of this equation by ɛ we obtain the Hamiltonian equation for the Euclidean space isotropic oscillator H p X + p Y + ω X + Y ) = h in agreement with ) Using the same procedure we find that the constants of the motion become K = X = xp y yp x and L = lim ɛ ɛ L = p Y + ω Y L = lim ɛ ɛ L = p X p Y ω XY an alternate basis for the symmetries of the isotropic oscillator Similarly the structure equations for the Higgs oscillator go in the limit to the structure equations of the isotropic oscillator 3 Examples of D nd degree nondegenerate systems 3 The system S7 on the sphere The classical system S7 on the -sphere [3] is determined by the Hamiltonian H = J + J + J3 a s + s s + s + a s + a 3s 3 s + s where J = s p 3 s 3 p and J J 3 are cyclic permutations of this expression We have embedded the -sphere in Euclidean 3-space so p + p + p 3 = J + J + J 3 + s p + s p + s 3 p 3 ) s + s + s 3 and we can use the Poisson bracket {F G} = 3 si F pi G + pi F si G) for our computations i= but at the end we restrict to the unit sphere: s + s + s 3 = and s p + s p + s 3 p 3 = The Hamilton equations for the trajectories s j t) p j t) in phase space are ds j dt = {H s j} dp j dt = {H p j} j = 3 The classical basis for the constants of the motion is L = J3 a s + s s + a s + s ) s s

4 Y Chen EG Kalnins Q Li and W Miller Jr L = J J 3 + a s 3 s + s ) + s a s s 3 s + s s + The structure relations are R = {L L } and {R L } = a a 3 + 4L L a 3 s s + s {R L } = 6L L + 4HL a H + 4a L + a 3 a R + 4L 3 4L H + 4L L 4a L + 4a HL a H a 3 a ) L + a a 3 L + a a 3 = 3 S7 in polar coordinates In terms of polar coordinates r θ where s = r cos θ s = r sin θ s 3 = ± r and r θ < π we have H = p θ r + r ) p r + a cos θ r sin θ + a r sin θ ± a 3 r r L = p θ + a cos θ + a 3) cosθ) L = ± r cos θ pθ r ) + sin θp r p θ ± r a + a cos θ + a cos θ) r cos + θ) a 3 cos θ Here ± is interpreted as + in the northern hemisphere and in the southern hemisphere Hamilton s equations give ṙ = r ) p r ṗ r = p θ r 3 θ = p θ r ṗ θ = a + a cos θ + a cos θ r sin 3 θ r r + rp r + a cos θ r 3 sin θ + a r 3 sin θ ± a 3 r ± a 3 Assumptions We require a 3 < and a a > and restrict our attention to trajectories caged in one hemisphere eastern or western): < s 3 < < s < < s < or r > π > θ > Note that the north pole is attractive and the south pole is repulsive 3 Trajectories Eliminating p θ p r in the expressions for H L L we find the implicit equation for the trajectories: = cos θa 3r 4L r cos θ + 4Hr L cos θ 4 cos θa 3 r L 8 cos θl L r r 4L a cos θr a cos θa 3 r r + 4Hr a cos θ r a + 4L 4L a 4a L r r + a + 4L r 4a a 3 r r + 4Hr a 4Hr L + 4a 3 r rl This is a quadratic equation for cos θ as a function of r with solutions cos θ = Nr) ± Sr) 3) Dr)

Examples of Complete Solvability of D Classical Superintegrable Systems 5 in the upper hemisphere where N = a a 3 r + a 3 rl Hra + 4L L r + rl a S = Hr ) r ra 3 L R D = r 4HL 4L + a 3) We note that N 4S = a 3 + 4HL 4L ) 4L r 4a 3 r ra + 4L a 3 r r + 4Hr a 4Hr L 4L a r a 4a L r r + 4L + a ) In the lower hemisphere the trajectories are given by cos θ = Ñr) ± Sr) Dr) Ñ = a a 3 r + a 3 rl Hra 4L L r + rl a 33) S = Hr + r ra 3 L ) R D = r 4HL 4L + a 3) Since R we see that if H L then S Thus for this case the trajectory is never in the lower hemisphere If H L an elementary analysis shows that Sr) > exactly in an interval < α r and is nondecreasing on that interval Case : a > a > All trajectories are closed and periodic We must have L > for trajectories and we initially assume R > The equations for perigee and apogee distance from the projection of the trajectory in the equatorial plane to the origin) are r = a 3 + HL ± a 3 a 3 + 4HL 4L a 3 + H ) Thus in order to have physical trajectories we must have a 3 + 4HL 4L and a 3 + HL Suppose the trajectory touches the equatorial plane at r = θ = θ Then taking a limit in 3) as the trajectory goes to the boundary we find that if the trajectory touches the equatorial plane it does so at an angle cos θ which is a solution of the quadratic equation ) 4 a 3 + H L )L cos θ a 3 L a H L )) cos θ + L H L )L a ) = Indeed cos θ = [ a 3 L + a L H)) ± H L )R ] /[a 3 + 4L H L )] Thus a necessary condition for the trajectory to reach the equatorial plane is that H L If the trajectory just touches the equatorial plane but doesn t go into the lower hemisphere then we must have H = L and a 3 > L > a 3 An example of touching the equatorial plane is Fig The different colors in the accompanying graphs correspond to the different curves 3) 33) that make up the trajectories If the trajectory passes through the equatorial plane then we must have H > L and 4L 4a L + a The angle of crossing lies in the interval a L L 4L 4a L + a cos θ a L + L 4L 4a L + a and for each θ in that interval the possible values of L are L = a 3 cos θ ± L H)cos θ L L + a cos θ + a )

6 Y Chen EG Kalnins Q Li and W Miller Jr 8 6 4 s3 4 6 8 5 8 6 4 5 4 6 8 s s Figure Touching the equatorial plane: a = 3 a = 5 a3 = 6 L = 6 H = 6 L = 8 6 4 s3 4 6 8 5 5 8 6 4 6 4 8 s s Figure Case : a = 3 a = 4 a3 = 5 L = 4 L = H = 6 8 6 4 s3 4 6 8 5 5 5 5 s s Figure 3 Case : a = a = 3 a3 = 4 L = 3 L = 5 H = 5 Note that for each choice of the constants of the motion there are always two crossing angles Thus if a trajectory crosses into the lower hemisphere it must return to the upper hemisphere: No trajectory remains confined to the lower hemisphere Examples are Figs and 3 If the angle is fixed we can use the formula to get L : Examples are Figs 4 and 5 If R = then Sr) The previous analysis is correct except the trajectory is a single arc rather than a loop We call this a metronome orbit The particle moves back and forth along the arc with a fixed period We give no more details here because we will study the analogous systems in our treatment of E6

Examples of Complete Solvability of D Classical Superintegrable Systems 7 8 6 4 s3 4 6 8 8 6 4 4 6 8 s s Figure 4 a = 8 a = a3 = L = 4 H = 5 L = 549 8 6 4 s3 4 6 8 8 6 4 s 4 6 8 s Figure 5 a = 8 a = a3 = L = 4 H = 5 L = 999 8 6 4 s3 4 6 8 5 5 5 5 s s Figure 6 a = a = a3 = 4 L = 5 H = 8 L = Case : a > a Then the portion < s < a /a of the s -axis becomes attractive the trajectories are not closed; each end of the trajectory impacts this interval assuming R > If H > L then one impact occurs in the northern hemisphere and one impact in the southern hemisphere If H L both impacts occur in the northern hemisphere An example for H > L is Fig 6 and for case H < L is Fig 7 If R = then the trajectory is a single arc that impacts the interval < s < a /a of the s -axis in the northern hemisphere

8 Y Chen EG Kalnins Q Li and W Miller Jr Figure 7 a = 4 a = a 3 = L = 5 H = L = 3 System E6 in Euclidean space as a contraction of S7 E6 can be described as a Kepler Coulomb system with barrier in two dimensions [3 3] The system has Hamiltonian H = p x + p y + a x + y + a y x x + y + a 3 x In terms of the rotation generator J = xp y yp x the defining constants of the motion are H and L = J + a y x + y y x + a 3 x L = p x J a y x + y a y x x + y + x + y ) a 3 y x The time derivative of a function F of the phase space parameters along a trajectory is given by F F = {H F } so in Cartesian coordinates we have the equations of motion t ẋ = p x ṗ x = a x 3 x + y ) 3/ x 4 + a y 3x + y ) + a 3 x + y ) 3/) ẏ = p y ṗ y = x + y ) 3/ a y a ) In polar coordinates x = r cos θ y = r sin θ we have p θ = xp y yp x rp r = xp x + yp y p x = sin θ r p θ + cos θp r p y = cos θ r p θ + sin θp r H = p r + p θ r + a r + a sin θ r cos θ + a 3 r cos θ L = p θ + a sin θ cos θ + a 3 sin θ cos θ L = cos θp r sin θ ) r p θ p θ a sin θ a sin θ r cos θ + ) a 3 sin θ r cos θ Setting R = {L L } we can verify that R = 4L + a 3 ) L H L ) + a L + a a L + a H 34) Note also that H = p r + a r + L +a 3 r θ = p θ r ṙ = p r ṗ r = L + a 3 r 3 + a r In polar coordinates the equations of motion are

Examples of Complete Solvability of D Classical Superintegrable Systems 9 ṗ θ = a cos r 3 cos 3 θ ) a 3 sin θ ) θ From these equations we see that points on a trajectory such that θ = are characterized by p θ = If such a point has coordinates θ r p r ) we find that sin θ = a ± a + 4L L + a 3 ) p θ = L + a 3 ) r = a + sin θ ) + a 3 sin θ sin θ )a sin θ L ) p r = H a L + a 3 r Similarly if ṙ = then p r = and if such a point has coordinates θ r p θ ) we find in particular that r = a ± a + 4HL + a 3 ))/H) Assumption : a < a > a 3 > Possible values taken by the constants of the motion depend on the relative sizes of a and a 3 Case : a 3 > a Then L + a 3 > and we must have a + 4HL + a 3 ) for trajectories Unbounded trajectories correspond to H and for H > perigee occurs at p r = r = a + a + 4HL + a 3 ) H p θ = L a sin θ L + a 3 ) sin θ sin θ sin θ = r L + a L + a 3 + a r For H = perigee occurs at p r = r = L + a 3 a p θ = L a sin θ L + a 3 ) sin θ sin θ sin θ = a L + a 3 + L a Bounded orbits correspond to H < in which case perigee occurs at p r = r = a + a + 4HL + a 3 ) H p θ = L a sin θ L + a 3 ) sin θ r L + a sin sin θ = θ L + a 3 + a r and apogee at p r = r = a a + 4HL + a 3 ) H p θ = L a sin θ L + a 3 ) sin θ r L + a sin sin θ = θ L + a 3 + a r We must have a + 4L L + a 3 ) for trajectories There is an equilibrium point at p r = p θ = r = L + a 3 sin θ = a 3 a 3 a a a Period of bounded orbit Since dr dt = H a r L +a 3 we have r dr r = Hr a r L +a 3 ) dt Thus the time needed to go from perigee r to apogee r is lim lim R r + R r R R r dr Hr a r L + a 3 ) = T r

3 Y Chen EG Kalnins Q Li and W Miller Jr Clearly this is the same time as needed to go from apogee to perigee Hence the period T is a [ T = lim H) 3/ arctan α lim arctan β] = a π α α H) 3/ Thus the period depends only on the energy The trajectories Eliminating p θ p r in the expressions for H L L we find the orbit equations r = Nsin θ) ± Ssin θ) N = a a 4a 3 L 4L L ) sin θ a L a L Dsin θ) S = L + a 3 ) sin θ + L a sin θ ) 4L + a 3 )HL L ) + a L + a H + a a L ) D = 4HL + a 3 ) + a ) sin θ + 4a L + 4a H) sin θ 4HL + 4L 35) From structure equation 34) and the fact that R > we see that the nd factor in S is a nonnegative constant Further from 34) we see that the denominator does not vanish for H < that it has a single root sin θ = L /a for H = and a double root for H > For the single root case with H = assuming the plus sign in 35) we set z = sin θ +L /a in 35) and expand r in a Laurent series series to find the asymptotic expression r = 4 R a 3 + a a 8L L 8a 3 L ) z a z 4R L + a 3 ) + 6L L + a 3 ) 8a a L L + a 3 ) + a ) a + 4 a R + Oz) Note that here R = 4L + a 3 )L + a L + a a L If the minus sign is assumed in 35) then we have the finite limit lim r = a [4L z 4R L + a 3 ) + a ] Parametrization of case trajectories Recall that a a a 3 are fixed constants with a < a 3 > a > To parametrize the trajectories we use perigee coordinates l t p r p where l = L + a 3 > < t p < r p > Then H = l rp + a r p L = t p r p l a t p a r p H > l r p > a H = l r p = a H < l r p < a Note that N + S ) N S ) = a + 4l a 3 )l )D 36) For H < apogee occurs for radial distance r a = a rp r p l +a r p l +a r p) For wedge boundaries S S radial lines containing a point on the trajectory where S vanishes) to exist at all we must have a + 4l a 3 )l Then the righthand boundary will always exist but the left hand boundary will exist only if l a However from the R > inequality below the left hand boundary must always exist) The location of the wedge boundaries is defined by t = sin θ where S : t = a a + 4l l a 3 ) R l + a r p = r p S : t = a + a + 4l l a 3 ) l l ) l a 3 a t p l t p)

Examples of Complete Solvability of D Classical Superintegrable Systems 3 Figure 8 Case : positive energy In order for the system to correspond to a trajectory we must have R We will first require that R > the generic case We see that this is violated if either ) l + a r p = or ) t p lies outside the open interval t t ) Thus we have the conditions a + 4l l a 3 ) < t p + a a < + 4l l a 3 ) l + a r p l l l For H > the zeros of the denominator occur at the values t ± of sin θ such that t ± = a r p t p a ± l + a r p )l a 3 a t p l t p) l + a r p For H > if the numerator of the trajectory equation 35) vanishes at sint) = t n then t n = t + or t and t < t n < t Since every zero of the denominator is a zero of one of the numerators we must have t < t < t + < t Further if l < a then D has at most one zero Thus if this case occurs then there is a single asymptote and a single wedge boundary Summary of necessary Case conditions for trajectories with R > ) l = L + a 3 > < t p < r p > ) a + 4l a 3 )l 3) l + a r p 4) a + a + 4l l a 3 ) l < t p < a + a + 4l l a 3 ) l For examples see Figs 8 9 3 3 Case a: R = In this case equation 36) for the trajectories simplifies since S We have N = a + 4l a 3 )l )D so rθ) = a + 4l a 3 )l N Because N is linear in t = sin θ it is easily seen that the trajectories are segments of ordinary conic sections: lines ellipses circles hyperbolas and parabolas Now R = l +a r p r p ) l a 3 a t p l t p) = so there are two possibilities: I) l + a r p = and II) l a 3 a t p l t p =

3 Y Chen EG Kalnins Q Li and W Miller Jr Figure 9 Case : negative energy Figure 3 Case : zero energy Figure 3 Case : zero energy Case ai): Circular arc l + a r p = We can express the constants of the motion in terms of r p t p alone: l = a r p H = a r p < L = a r p

Examples of Complete Solvability of D Classical Superintegrable Systems 33 The equation for the trajectories simplifies to rθ) = a +4l a 3 )l N = r p so the trajectory is a segment of a circle centered at the origin with radius r p The wedge boundaries are t = a a +a r p +a a 3 r p a r p t = a + a +a r p +a a 3 r p a r p Thus for this case to occur r p must satisfy a + a r p + a a 3 r p We conclude that this case always occurs for r p sufficiently small or sufficiently large but that there is a forbidden intermediate band for which no circular trajectory exists Case aii): Conic section arc l a 3 a t p l t p = Again we can express the constants of the motion in terms of r p t p alone: l = a 3 + a t p t H = a 3 + a t p p rp t p) + a L = t pa 3 + a t p ) r p r p t a t p p) a r p The equation for the trajectories becomes where r p a + t p a 3 + a t p) rθ) = ± a t p a t pr p + a r p + a 3 ) sin θ + t p)a a t p r p ) = ± α sin θ + β α = a t p a t pr p + a r p + a 3 ) r p a + t p a 3 + a t β = t p)a a t p r p ) p) r p a + t p a 3 + a t p) and the sign is chosen so that r > Setting y = r sin θ x = r cos θ we see that the trajectory is a segment of the curve αy + βr = ± or βr = αy ± so β x + y ) = α y ± αy + α + β > We have the following possibilities for the curve segments: β > α > : ellipse < β < α : hyperbola β = α : parabola α = : circle β = : horizontal line The wedge boundaries are t p a +a 3 t p a 3 +a t p For the circular arc we have r p = a t p+a 3 a t p ) This coincides with Case ai) For the horizontal line we have t p = a r p a < so trajectories are possible only for r p < a /a Only strictly positive energy solutions are possible for the horizontal line For the parabolic arc we have r p = a 3 a a t p + a t p/a + t p ) t p ) The energy is always negative Only strictly negative energies are possible for the elliptical arc In general orbits are possible for H < if and only if the constants of the motion satisfy the inequality a 3 + a t p / t p) a ) < r p Trajectories are possible in the positive energy case if and only if r p < a 3 + a t p / t p) a ) Case : a 3 < a Now there are no restrictions on the sign of l = L + a 3 The force in the angular direction is always negative In the first quadrant the force in the y direction is always negative It follows from this that there are no periodic orbits contained entirely in the first quadrant The positive y-axis is repelling but the negative y-axis is attracting We assume initially that R > Equations 35) for the trajectories still hold The condition a + 4l l a 3 ) > is now satisfied automatically There is at most a single wedge boundary t = sin θ = a + a + 4l l a 3 ) l < The following general types of behavior occur: If l < the force in the radial direction is always negative The trajectories all impact the lower y-axis and are perpendicular to the axis at impact Since ṗ r is always strictly negative

34 Y Chen EG Kalnins Q Li and W Miller Jr along a trajectory perigee must occur at y p = r p on the negative y-axis Note that R is linear in H always with positive linear coefficient The critical value of H such that R = is E = [4a 3 L + 4L L a a L a L ]/[4L + a + 4a 3L ] The qualitative behavior of the trajectories divides into two basic classes depending on the sign of Nsin θ) for θ = π Note that N ) = L a )l a ) + a a 3 a ) S ) = a a 3 )R Case a: N ) > Then if we increase H from E to E + ɛ for arbitrarily small ɛ > there will be two points of intersection of the trajectory with the negative y-axis one on the curve r = N + S)/D and one on the curve r = N S)/D The trajectory will remain bounded and there will be a wedge boundary If we further increase H leaving L L unchanged this behavior will persist until the critical value E = 4 [L a ) ]/[a a 3 )] E where D ) = Here the intersection point with the negative y-axis and the curve r = N S)/D has moved to y = whereas the intersection point with the curve r = N + S)/D remains finite As H is further increased the behavior of the trajectory is that it is unbounded and asymptotic to a radial line in the nd quadrant that there is a wedge boundary and a single intersection point with the negative y-axis Case b: N ) < Again we increase H from E to E + ɛ for arbitrarily small ɛ > Now there are no points of intersection with the negative y-axis and in fact the trajectory is non-physical until H is increased to ɛ beyond the critical value E The curve is entirely described by the function r = N + S)/D The trajectory is unbounded asymptotic to a radial line in the nd quadrant and intersects the negative y-axis a single time There is no wedge boundary Note that for this case the region {N ) < H E } is nonphysical We do not have a complete proof of this but strong numerical evidence) If l > the behavior is different The trajectories do not necessarily intersect the negative y-axis Some trajectories are disjoint The qualitative behavior of the trajectories divides into two basic classes depending on the sign of Nsin θ) for θ = π Note again that N ) = L a )l a ) + a a 3 a ) S ) = a a 3 )R Case c: N ) > If we increase H from E to E + ɛ for arbitrarily small ɛ > there will be two points of intersection of the trajectory with the negative y-axis one on the curve r = N + S)/D and one on the curve r = N S)/D The trajectory will remain bounded if H < with a wedge boundary) and be unbounded if H > The unbounded trajectory will divide into two disconnected parts If we further increase H leaving L L unchanged this behavior will persist until the critical value E = 4 [L a ) ][a a 3 )] E where D ) = Here the intersection point with the negative y-axis and the curve r = N S)/D has moved to y = whereas the intersection point with the curve r = N + S)/D remains finite As H is further increased the behavior of the trajectory is that it is unbounded and asymptotic to a radial line in the nd quadrant and a single intersection point with the negative y-axis There is a wedge boundary Case d: N ) < Again we increase H from E to E + ɛ for arbitrarily small ɛ > There are no points of intersection with the negative y-axis and the trajectory is non-physical for H < but physical for the H > interval until H is increased to ɛ beyond the critical value E Then the trajectory is unbounded and asymptotic to a radial line in the nd quadrant Physical trajectories have a wedge boundary See Fig 3 3 R = In this case equation 36) for the trajectories simplifies since S We have N = a + 4l a 3 )l )D so rθ) = ± a + 4l a 3 )l N a = ± + 4l a 3 )l a a 4L l )t + a a 3 a L a l ) = ± α sin θ + β

Examples of Complete Solvability of D Classical Superintegrable Systems 35 Figure 3 Case d: N ) < where α = a a 4L l a +4l β = a a 3 a L a l ) a 3 )l a +4l and the sign is chosen so that r > Setting a 3 )l y = r sin θ x = r cos θ we see that the trajectory is a segment of the curve αy + βr = ± or βr = αy ± so β x + y ) = α y ± αy + α + β > We have the following possibilities for the curve segments: β > α > : ellipse < β < α : hyperbola β = α : parabola α = : circle β = : horizontal line Since R = we can express the constants of the motion in terms of α β alone: H = a α 3 a α a 3 α β + a αβ + a β β L = a α a α a 3 αβ + a β β l = a α a β Case e: circular arc a a 4L l = Then H = a 4l The radius is r = ± l a There is a wedge boundary at t = sin θ = a + a +4l l a 3 ) l The trajectory impacts the negative y-axis See Fig 33 Case f: horizontal line segment a L + a l = a a 3 y = a a The energy is H = a l a 3 ) If the energy is negative there is a wedge boundary at t = sin θ = a + a +4l l a 3 ) l < and the trajectory impacts the negative y-axis If the energy is positive the trajectory is unbounded and there is no wedge boundary but the trajectory still impacts the negative y-axis Case g: parabolic arc If α = β we have L = a a a 3 + l ) H = a a 3 )a l a ) l a ) and the energy is positive If α = β we have L = a a + a 3 l ) H = a + a 3 )a l + a ) l + a ) and the energy is negative There is a wedge boundary at t = sin θ = a + a +4l l a 3 ) l and the trajectory impacts the negative y-axis a

36 Y Chen EG Kalnins Q Li and W Miller Jr Figure 33 Case e: circular arc Case h: elliptic arc H = 4l L a a L a l a 3 ) a + 4l l a 3 ) the general expression There is a wedge boundary at t = sin θ = a + a +4l l a 3 ) l and the trajectory impacts the negative y-axis Case i: hyperbolic arc Here again H = 4l L a a L a l a 3 ) a + 4l l a 3 ) the general expression There is a wedge boundary at t = sin θ = a + a +4l l a 3 ) l and the trajectory impacts the negative y-axis There are no radial line trajectories in this case 4 The Post Winternitz superintegrable systems These systems are especially interesting as a pair of classical and quantum superintegrable systems that do not admit separation of variables so they cannot be solved by traditional methods We treat only the classical case here and show that its superintegrability allows us to determine the trajectories exactly The quantum case can also be solved 4 The classical system The classical Post Winternitz system is defined by the Hamiltonian [36] H = p x + p ) cy y + x /3 The generating constants of the motion are of 3rd and 4th degree: L = 3p xp y + p 3 y + 9cx /3 p x + 6 cyp y x /3 L = p 4 x + 4cp y x x /3 cx/3 p x p y c 9x y ) x 4/3

y y Examples of Complete Solvability of D Classical Superintegrable Systems 37 5 5 5 5 3 6 5 4 3 x Figure 34 c = 5 E = L = L = 4 6 4 8 6 4 5 5 5 3 35 x Figure 35 c = E = 4 L = 3 L = 6 Here {L L } = 88c 3 so the constants of the motion generate a Heisenberg algebra Since there is no nd degree constant other than H separation of the Hamilton Jacobi equations is not possible in any coordinate system so separation of variables methods cannot be used to compute the trajectories However we can find a parametric description of the trajectories Choosing constants H = E L = l and L = l we can verify the two identities p 3 y + 9cx /3 p x l + 6Ep y = 3 p4y + 4py E + 4E l 43 l p y ) p x 9 p 3 y l + 6Ep y ) = From these results we can solve for p x and x y as functions of the parameter p y : p 3 ) y l + 6Ep y p x = 6p 4 y + 7Ep y 4l p y + 7E 8l x = 583c 3 6p 4 y + 7Ep y 4l p y + 7E 8l ) 3/ y = 34c 3 l 36E p y + p 6 y + 8l p 3 y 8Ep 4 y + 9l p y + 7E 3 8l E ) 4) and equations 4) enable us to plot the trajectories See Figs 34 35 and 36 There are no closed orbits The trajectories always impact the y-axis from the left or right depending on the signs of the parameters 5 Classical elliptic superintegrable systems 5 Elliptic superintegrability Here we extend the ideas in [4 4] to a system that separates in Jacobi elliptic coordinates [35] We start with a simple example the Higgs oscillator S3 ) [3] expressed in elliptic

y 38 Y Chen EG Kalnins Q Li and W Miller Jr 9 8 7 6 5 4 3 5 5 5 x Figure 36 c = E = L = 5 L = 5 coordinates H = sn α k) sn β k) p α p β + A sn α k) ) A sn β k) The Hamilton Jacobi equation can be taken as H = E for p α = W α p β = W β The separation equations in the coordinates α β are p α E sn α k) + A sn α k) + Λ = p β E A sn β k) + sn β k) + Λ = 5) Here H and L = p α + sn A α k)h sn αk) are constants of the motion and L = Λ on a coordinate hypersurface To find another constant of the motion via the action-angle variable construction we introduce M and N which are computed according to M = dα N = dβ p α p β where p α p β are expressed in terms of coordinates α β by 5) Applying action-angle theory see eg [ ] we see that M N is a constant of the motion The integral M can be written as M = snα k)dα E sn 4 α k) Λ sn α k) A with a similar integral for N This can be calculated as an elliptic integral by choosing the new variable ρ = sn α k) We then have the relations dρ dα = snα k)cnα k)dnα k) cnα k) = ρ dnα k) = k ρ The change of variables is then tantamount to the computation of M = cnα k)dnα k)snα k) dα cnα k)dnα k) E sn 4 α k) Λ sn α k) A = 4k dρ E ρ a)ρ b)ρ c)ρ d) where a = b = and c d are roots of ρ Λ k E ρ A E A similar calculation applies to the substitution µ = sn β k) with exactly the same a b c and d A convenient choice of integral is dt t a)t b)t c)t d)

Examples of Complete Solvability of D Classical Superintegrable Systems 39 ) = sn b d)t a) a c)b d) a d)t b) b c)a d) a c)b d) This integral can be used for t = ρ or µ If we form snk Ea c)b d)m N) κ) for κ = b c)a d)/a c)/b d) and use the addition formula [35] snu + v κ) = snu κ)cnv κ)dnv κ) + snv κ)cnu κ)dnu κ) κ sn u κ)sn v κ) we obtain a constant of the motion In terms of the coordinates it is helpful to write p α = ρ c)ρ d) E p β = µ c)µ d) E ρ µ Then for Z = b d a c a d we have ρ a b a snu κ)cnv κ)dnv κ) = Z µ c)µ d) ρ b µ b with a similar expression for snv κ)cnu κ)dnu κ) We also note that κ snρ κ) snµ κ) = ρ a)µ a)b c)b d) ρ b)µ b)a c)a d) as well as µ c)µ d) = µ µ) k µ) p µ E ρ c)ρ d) = ρ ρ) k ρ) p ρ E ρ )d k snu κ) = ) cnu κ) = dnu κ) = d )ρ k ) ρ c)k ) c ) ρ k ) κ = ρ d)k ) d ) ρ k ) d )c k ) c )d k ) with similar expressions for snv κ) cnv κ) and dnv κ) with ρ replaced by µ From this we deduce that E Z snk Ea c)b d)m N) κ) = X where X is a constant of the motion linear in the momenta viz X = classical Hamiltonian k ρ) k µ) ρ) µ) ρ µ) H = ρ ρ) k ρ ) p ρ µ µ) k µ ) p A ρ µ µ) ρµ ρp ρ µp µ ) in involution with the This is just the first degree symmetry that we expect for this Hamiltonian In this case we choose the variables ρ and µ to obtain the constant In particular cd = A/E c + d = Λ /E This seems to be just a complicated way of deriving a straightforward result However we can use the TTW trick [4 4] by looking at the Hamiltonian H = sn pα k) sn p α p β qβ k) + A sn pα k) ) A sn qβ k)

4 Y Chen EG Kalnins Q Li and W Miller Jr where p and q are arbitrary nonzero numbers The separation equations are p α E sn A pα k) + sn pα k) + Λ = p β E sn qβ k) + A sn qβ k) + Λ = 5) A Here H and L = p α sn pα k)h + sn are constants of the motion To find another pαk) constant of the motion via the action-angle variable construction we introduce M and N which are computed according to M = dα p α N = dβ p β where p α p β are expressed in terms of coordinates α β by 5) Again M N is a constant of the motion The integral M can be written as M = snpα k) dα E sn 4 pα k) Λ sn pα k) A with a similar integral for N This can be calculated as an elliptic integral by choosing the new variable ρ = sn pα k) We then have the relations dρ = p snpα k)cnpα k)dnpα k) dα cnpα k) = ρ dnpα k) = k ρ The change of variables is then tantamount to the computation of M = cnˆα k)dnˆα k)snˆα k) dˆα p cnˆα k)dnˆα k) E sn 4 ˆα k) Λ sn ˆα k) A = 4k dρ Ep ρ a)ρ b)ρ c)ρ d) where ˆα = pα a = b = and c d are roots of ρ Λ k E ρ A E A similar calculation applies to the substitution µ = sn qβ k) with exactly the same a b c and d A convenient choice of integral is dt t a)t b)t c)t d) ) = sn b d)t a) a c)b d) a d)t b) b c)a d) a c)b d) This integral can be used for t = ρ or µ If we now form snk Ea c)b d)pqm N)) and use the addition formula snqu + pv κ) = snqu κ)cnpv κ)dnpv κ) + snpv κ)cnqu κ)dnqu κ) κ sn qu κ)sn pv κ) we obtain a constant of the motion Thus computation yields a rational constant of the motion whenever p/q is a rational number Let us restrict ourselves to p = q = Then u = α/ and v = β where p α = E ρ c)ρ d) p β = E ρ µ c)µ d) b d Y = µ a d

Examples of Complete Solvability of D Classical Superintegrable Systems 4 We form T = snu κ)cnv κ)dnv κ) ρ a)ρ b)ρ c)ρ d)µ c)µ d) = Y µ b)ρ a + b c d) + ρcd ab) cda + b) + abc + d)) T = snu κ)cnv κ)dnv κ) µ a = Y c + b a d)ρ + ad bc)ρ + cba + d) adb + c) ) µ b d + b a c)ρ + ac bd)ρ + cbd a) + adb c) ) We also note that a + b c d)ρ + cd ab)ρ + abc + d) cda + b) ) T 3 = κ snu κ) snµ κ) µ a)ρ a)ρ b)ρ c)ρ d)b d)b c) = 4 µ b) a + b c d)ρ + cd ab)ρ + abc + d) cda + b)) µ c)µ d) = 4µ µ) k µ) p µ E ρ c)ρ d) = ρ ρ) k ρ) p ρ E and ρ = sn α k) and µ = sn β k) From this we deduce that Y sn4k Ea c)b d)m N) κ) = P where P is a constant of the motion rational in the momenta We can calculate the components of this constant as follows: T = n /d where n = 6 ρ) k ρ ) ρµk µp ρ p µ d = k µ ) k Λ E + k) ) ρ + Ak + E ) ρ Λ Ak + ) ) T = n /d where n = µ [ E k 4 + k 4 Λ + 4AE )) ρ 4 + k 4 8AE + Λ E Λ ) + k 8AE 4Λ E Λ ) + Λ E ) ρ 3 + 6 k AE + AE + k Λ + AE) ) ρ + k 4 Λ A + k Λ + Λ A 4AE ) Λ Λ A 8AE ) ρ + k 4 A + 4AE + k ) A Λ )] d = k µ [ k Λ E + k )) ρ + Ak + E ) ρ Λ A k + )] T 3 = n 3 /d 3 where n 3 = 4ρ ) ρk ) µ ) Eρ + Λ ρ + A ) E + Ak 4 + Λ k ) d 3 = k µ ) [ k Λ E + k )) ρ + Ak + E ) ρ Λ A k + )] We could of course use the expressions above given in terms of c + d and cd If we set x = sn α k) = ρ and y = sn β k) = µ then H and L become H = x x) k x ) p x 4y y) k y ) p A x y y) xy 53) L = xy x) k x ) p x 4 y) yk ) p ) y A x y x + ) y E

4 Y Chen EG Kalnins Q Li and W Miller Jr For A = the extra constant of the motion is G = N/D where N = y yxk x k x + x y ) p x 4yy ) + k y ) p y + 4xy ) + k y ) ) p x p y D = y x k + x yk x xyk + y ) p x + 4yy ) + k y ) p y 4y + k y ) ) x )p x p y Setting R = {L G} we obtain the polynomial structure relations R 4 k G ) L + G H L G ) = {G R} = G ) G k ) {L R} = 4G k ) L + k HG k L G H ) Turning to the case when A we have G = N/D where N = y [ x y x yk + x k x xyk + y ) p x + 4x y y ) + k y ) p y 4x y + k y ) p x p y + Ax y) ] D = + k y [ x y xy + y k x + x yk + x ) p x + 4x y y ) + k y ) p y 4x y y ) + k x ) p x p y + Ax y) ] The structure relations are {G R} = G ) k G ) {L R} = 4G [ L k G + Ak 4 G HG L L k Ak + H ] 54) R + 4 [ Ak 4 G L k G + L G 4 Ak G + A L G + HG + L H H 4] = All of the examples that we have constructed using the addition theorem for elliptic functions obey polynomial structure equations but we have no proof that this is true in general 5 Example H = x ) y x p x 4 y ) p ) A y x ) y ) L = x ) y ) y x ) p x 4p y) + Ax + y ) x ) y ) xyx y ) y ) x ) y ) x ) L = yx4 x y + y 3 )p x + 4x y + 4x y 3 )p y + 4x 3 y + 4x 3 )p x p y A x 4 y x + y )p x + 4x 4x y )p y + 4xy + 4xy 3 )p x p y + A x x y ) This is the special case of 53) 54) where we have replaced x y by x y respectively set k = and rescaled Note that H = 4p y + L y ) + A and x ) y ) y x ) p x 4p ) y = L y ) Ax +y ) We have {H L x ) y ) } = {H L } = and with R = {L L } the polynomial structure relations are {L R} 6L L + 3HL 3L 3 H = {L R} 8 L ) = R + 6L 6L L + 3HL 6L 4 H 6H + 6A =

Examples of Complete Solvability of D Classical Superintegrable Systems 43 Using the three constants of the motion to eliminate the momenta terms we find p x = Hx ) + L x ) A x ) p y = Hy ) + L y ) A 4 y ) and the trajectories must satisfy the implicit equation = Hx A + L L x Hx 4 + Hx H)y + 4 H + L + Hx 4 + A ) ) L 3Hx 4 L H x 8 A + AH 3H x 4 + L x A + H x 6 + 3Hx 4 A + AL L Hx 6 + L H + H x L + L x ) H 3L x H Hx A yl + H + L + Hx 4 + A ) y 4 x A + L H) A + L L x Hx 4 + Hx H ) The left-hand side of this equation is quadratic in y so we can solve for y as a function of x We obtain the solutions y = N± S D where N = L + x) 4 x + ) 4 H L Ax 3L x L x 6 + 3L x 4 + 3Ax 4 + A + L ) H + L A + L ) A + L L x ) S = x A + L L x Hx 4 + Hx H ) H + L + Hx 4 + A ) R 6 D = + x) x + ) x + L x + ) x + L x ) H + 4L Ax + L x 4 L + 4L L x A 4L L x 4 + Ax 4) H + A + L ) Let x y ) be a point on the trajectory such that p x = A straightforward calculation gives x = H L ± L + 4AH y = L H p y = A H HL4 + HL + L L L 4L ) L + ) Assumption: A > Note that x L ± ) L = + + 4AH 55) H It follows that if we are in the region x > and H > then since L + 4AH > L only the plus sign is possible in 55) We conclude that the trajectory must be unbounded to the right Examples are Figs 37 and 38 The different colors in the figures correspond to different functions describing the trajectories) We investigate what happens if the trajectory intersects the lines y = ±x If the intersection occurs for y = x then we find that x = L or p x = p y or both If the intersection occurs for y = x then we find that x = L or p x = p y or both In both cases we have p x = 4p y = HL ) +L L ) A L ) Thus in both cases if L the trajectory intersects the lines y = ±x at finite momentum However the velocity is infinite Indeed ẋ = {H x} = p x x )/y x ) which blows up unless p x = This suggests that there are no bounded orbits in the square < x y < However we find some bounded orbits in the region x > for negative energies Examples of these bounded orbits are Figs 39 and 4

y y y y 44 Y Chen EG Kalnins Q Li and W Miller Jr 5 4 y+ y 3 3 4 3 4 5 6 7 8 9 x Figure 37 Case: A = 5 L = L = 4 H = 5 4 y+ y 3 3 4 3 4 5 6 7 8 9 x Figure 38 Case: A = L = L = H = y+ y 5 5 3 35 5 5 3 35 x Figure 39 Case: A = 5 L = L = 3 H = 5 y+ y 5 5 5 x Figure 4 Case: A = L = L = H = 53 Another version of the example We set y = iz where z is real multiply H by and L by i to get the system H = x )p x + 4 + z )p z A x + z + x ) + z )

Examples of Complete Solvability of D Classical Superintegrable Systems 45 L = + z ) x )p x + 4p z) x + z + Ax z ) + z ) x ) Ax zz +x ) +z ) x ) L = zx4 p x 4x 3 p z p x z 4x 3 p z p x x zp x + 4x p zz + 4x p zz 3 p xz 3 x 4 p x 4p zx 4p x zp z x 4p zx z + p xz x p xz 4xz 3 p z p x + Ax x +z ) +z ) x ) ) We will consider this system for A > and resticted to the vertical strip of points x z) where < x < Then H must be > The motion is unbounded there are no bounded orbits and the origin ) and the boundary lines x = ± are repulsive The polynomial structure relations are R = {L L } and {L R} + 8 + L ) = {L R} 6L L 3HL 3L 3 H = R + 6A + 6L + 6L L + 3HL + 6L 4 H + 6H = Expressing the momenta in terms of the coordinates we have p x = A + H + x4 H + L x L Hx x ) p z = Hz + Hz 4 + L z + A + H + L ) 4 + z ) The equations for the trajectories turn out to be y = N± S D where N = L x ) 4 H + L L + 3Ax 4 + A L x 6 + 3L x 4 Ax 3L x ) H + L L + A) L + A L x ) S = x A + H + Hx 4 + L L x Hx ) A + H + L Hx 4) R 6 D = x ) x 4 + 4L x + x + ) H + Ax 4 + L 4L L x 4 + A L x 4 + 4L L x + 4L Ax ) H + L + A) The points x z p x p z ) on the trajectory where p x = satisfy z = L p z = A + H + L ) + L + L L ) 4 + L ) x = + L ± L 4AH H Requiring first that R > we see that L < is a necessary condition for trajectories otherwise S < in the strip < x < If L /H > then a further necessary condition for S > in a strip is L + H + A < If L /H then a necessary condition is L 4AH > Examples here are Figs 4 and 4 Note that here there are 4 different trajectories corresponding to a single choice of constants of the motion Case when R = When L = 9 L = A = H = 8 we have R = such that S = The trajectory in this case is a straight line on the x-axis The two boundary points where p x = for x > are and 9354 The velocity of the particle is ẋ = {H x} = px x ) Substituting x +z L z = and the expression for p x in terms of x into the equation givesẋ = +H H Thus x ẋ = so the particle will go through the origin instead of bouncing back lim x

46 Y Chen EG Kalnins Q Li and W Miller Jr Figure 4 Case: A = L = 5 L = 3 H = Figure 4 Case: A = L = L = 3 H = A Review of some basic concepts in Hamiltonian mechanics Hamiltonian mechanics is a reformulation of classical Newtonian mechanics alternate to Lagrangian mechanics In the Hamiltonian formalism a physical system describing the motion of a particle at time t involves n generalized coordinates qj t) and n generalized momenta pj t) The phase space of the system is described by points pj qj ) Rn The generalized coordinates qj s and generalized momenta pj s are derived from the Lagrangian formulation The dynamics of the system are given by Hamilton s equations [ ] dqj H =+ dt pj dpj H = dt qj A) where H = Hq qn p pn t) is the Hamiltonian which in this paper corresponds to the total time-independent) energy of the system: H = T + V where T and V are kinetic and potential energy respectively Solutions of these equations give the trajectories of the system