A New Superintegrable Hamiltonian P.E. Verrier and N.W. Evans Institute of Astronom, Madingle Rd, Universit of Cambridge, CB3 HA, UK Dated: Ma 29, 18) arxiv:712.3677v1 [nlin.si] 21 Dec 7 Abstract We identif a new superintegrable Hamiltonian in 3 degrees of freedom, obtained as a reduction of pure Keplerian motion in 6 dimensions. The new Hamiltonian is a generaliation of the Keplerian one, and has the familiar 1/r potential with three barrier terms preventing the particle crossing the principal planes. In 3 degrees of freedom, there are functionall independent integrals of motion, and all bound, classical trajectories are closed and strictl periodic. The generalisation of the Laplace-Runge-Len vector is identified and shown to provide functionall independent isolating integrals. The are quartic in the momenta and do not arise from separabilit of the Hamilton-Jacobi equation. A formulation of the sstem in action-angle variables is presented. PACS numbers: Valid PACS appear here Kewords: Suggested kewords Electronic address: pverrier@ast.cam.ac.uk Electronic address: nwe@ast.cam.ac.uk 1
I. INTRODUCTION The Kepler problem is well known to be superintegrable that is, it has five functionall independent integrals of motion. The are the energ and the components of angular momentum and the Laplace-Runge-Len vectors, obtainable b separating the Hamilton-Jacobi equation in spherical polar and rotational parabolic coordinates [1]. For Hamiltonians with three degrees of freedom, the eistence of five integrals of motion implies that ever bound trajector is closed. Sommerfeld and Born, in the das of the old quantum theor, appear to have been the first to realie that if a potential is separable in more than one coordinate sstem, it possesses additional isolating functionall independent integrals [2, 3]. The first sstematic inquir into this problem was begun b Winternit, Smorodinsk and co-workers, who found ever potential in two degrees of freedom for which the Hamilton-Jacobi equation is separable in more than one wa [4]. Subsequentl, the etended this work to three degrees of freedom b finding ever potential separable in spherical polars and at least one additional coordinate sstem []. Evans [6] then completed this work b investigating all the remaining possibilities. A useful introduction to the subject of superintegrabilit, as well as summar of recent work, is given in the conference proceedings of Tempesta et al. [7]. If the Hamilton-Jacobi equation separates, then the corresponding integral of motion is necessaril linear or quadratic in the canonical momenta. Consequentl, all the superintegrable sstems listed in [4,, 6] have integrals that are quadratic in the momenta. As an eample, let us consider the Keplerian Hamiltonian H = 1 2 ṙ 2 k r, 1) where k is a real positive constants. Letting r denote the position vector, then the integrals are the energ E, the components of the angular monetum vector L = r ṙ and the Laplace-Runge-Len vector all of which are at most quadratic in the velocities. A = ṙ L k r, 2) r Superintegrable sstems with higher-order integrals are known [8], although the are etremel scarce. Eamples include the anisotropic harmonic with rational frequenc ratio 2
l : m : n where l+m+n, i.e., H = 1 2 ṙ 2 +l 2 2 +m 2 2 +n 2 2. 3) The potential separates in rectangular cartesians and possesses two commuting quadratic integrals. There are two additional integrals which ma be taken as polnomials of degree l+m 1 and l+n 1 [9]. The Calogero potential in a harmonic well H = 1 2 ṙ 2 +k 2 + 2 + 2 )+ k 1 ) + k 1 2 ) + k 1 2 ) 2, 4) is known to be super-integrable and possesses an integral of the motion that is cubic in the velocities [, 11]. In this paper, we introduce a new superintegrable Hamiltonian, namel H = 1 2 ṙ 2 k r + k 1 + k 2 2 + k 3 2 2. ) This is recognied as a generaliation of the familiar Keplerian Hamiltonian. The constants k 1,k 2 and k 3 are taken as positive, so the Hamiltonian is perfectl phsical and motion is confined to, sa, the octant >, > and >. The terms involving the k i correspond to repulsive barriers preventing the orbit crossing an of the principal planes. In Section II, we present numerical integrations showing that the orbits give closed curves in all cases. This motivates a search for the isolating integrals, one of which is found to be quartic in the momenta in Section III. Finall, an action angle formalism is given in Section IV and the relation to the Kepler problem discussed. II. EVIDENCE OF SUPERINTEGRABILITY A. Analtic Proof Let us recall that Keplerian motion in N degrees of freedom alwas possesses 2N 1 functionall independent integrals of motion [12]. Specialising to 6 degrees of freedom, we have the Hamiltonian H = 1 2 p 2 k s 6) 3
where s has Cartesian coordinates s 1,s 2,s 3,s 4,s,s 6 ). Now, let us introduce coordinates,,,θ,θ,θ ) according to s 1 = cosθ, s 3 = cosθ, s = cosθ, s 2 = sinθ s 4 = sinθ s 6 = sinθ The Hamiltonian becomes H = 1 2 p2 1 +p2 2 +p2 3 + p2 θ 1 2 + p2 θ 2 2 + p2 θ 3 2 ) k 2 + 2 + 2 ) 1 2 The coordinates θ 1,θ 2,θ 3 ) are ignorable, so we obtain a new Hamiltonian H = 1 2 p2 1 +p2 2 +p2 3 ) k r + k 1 2 + k 2 2 + k 3 2 7) where k 1,k 2 and k 3 are the constant values of the momenta conjugate to the ignorable coordinates. In the original 6 degrees of freedom Hamiltonian 6), ever bound trajector is closed. Consequentl, in the reduced 3 degrees of freedom Hamiltonian 7), ever bound trajector is also closed. Evidentl, the proof can be readil generalised to N degrees of freedom. B. A Sampler of Orbits It is interesting to investigate characteristic orbits corresponding to the Hamiltonian 7). Using a standard Burlisch-Stoer code [13], the motion of a unit mass particle was followed in the octant with >, > and >. Stepsie and tolerances were set to maintain accurac to a level of around 12 relative energ change, and integration lengths were tpicall tens of periods. Starting with initial conditions that would give a circular orbit in a true Keplerian potential, the effect of relativel weak barriers can be seen in Figure 1. The orbit appears similar in shape to the Keplerian ellipse but reflects off the three aes planes. In this case, as a consequence of the k i being equal and the smmetr in the initial phase space position, the orbit is confined to a plane. The effect of larger barriers can be seen through increasing the k i b a factor of ten, as shown in Figure 2. The orbit is now further distorted, and takes the form of a figure-of-eight. If the k i are not equal, the orbit is still closed, as shown in 4
2 E=-.49 2 2 1 1 1 1 2 1 2 1 2 =.8,.8,.8) p =.,.3,.) k i =.1,.1,.1) k = 1. FIG. 1: A perturbed circular Keplerian orbit in the potential. The three panels show the projections in the, and planes respectivel from left to right. The initial conditions are given to 2 s.f.) below the plot and the energ in the top corner. Note that although in this case the orbit lies in a plane this is not generall true. 2 E=-.41 2 2 1 1 1 1 2 1 2 1 2 =.8,.8,.8) p =.,.3,.) k i =.,.,.) k = 1. FIG. 2: As for Figure 1 but now the centrifugal barriers are an order of magnitude larger. Figures 3 and 4. The latter of these two cases has a different set of initial conditions, which would place it on an initiall elliptical orbit in the true Keplerian problem. Man more initial conditions and combinations of parameter values were investigated. In all bound cases, ever orbit is closed and strictl periodic, which is reassuring confirmation of the eistence of a fifth isolating integral.
2 E=-.24 2 2 1 1 1 1 2 1 2 1 2 =.8,.8,.8) p =.,.3,.) k i =.,.,.) k = 1. FIG. 3: As for Figure 2 but now the centrifugal barriers are different and the motion is no longer confined to a plane. 2 E=-.28 2 2 1 1 1 1 2 1 2 1 2 =.,.,.) p =.1,., -.2) k i =.,.,.3) k = 1. FIG. 4: This orbit in the potential has initial conditions that would place it on an elliptic orbit in the original Kepler problem. III. THE INTEGRALS OF MOTION Makarov et al. [] and Evans [6] alread showed that all Hamiltonians of the form H = 1 2 ṙ 2 +Fr)+ k 1 + k 2 2 + k 3 2 2. 8) possess four isolating integrals of motion, arising from separabilit of the Hamilton-Jacobi equation in the spherical polar and conical coordinate sstems. Here, Fr) is an arbitrar function of the spherical polar radius. The four isolating integrals are the energ E and 6
three generaliations of the angular momentum components, namel E = 1 2 p 2 +Fr)+ k 1 2 + k 2 2 + k 3 2 9) I 1 = 1 2 L2 1 + k 2 2 2 + k 3 2 I 2 = 1 2 L2 2 + k 1 2 2 + k 3 2 I 3 = 1 2 L2 3 + k 1 2 2 + k 2 2 2 ) 2 11) 2 12) where p and L are the linear momentum and angular momentum vectors. As Eq ) is of this form, four of the integrals of motion are alread known and arise from separabilit. The pule is that there is a fifth integral whose form is unknown and which does not arise from separabilit. In the case where one of the barriers, sa that in the = = plane, vanishes, the fifth integral is known to be [6] I 4 = L 1 p 2 p 1 L 2 2 k 2r + k 1 + k ) 2 2 2 and follows from separabilit in the rotational parabolic coordinate sstem. If the fifth integral for the general problem ) is quartic, it must reduce to the above integral in the limit k 3. This suggests taking the ansat I 4 = L 1 p 2 p 1 L 2 2 k 2r + k 1 + k )) 2 2 +k 2 2 3 g,p) 14) where g,p) is a function of both position and momentum et to be determined. Requiring the Poisson bracket of I 4 with the Hamiltonian to vanish leads to a solution for g. So, we arrive at an isolating integral of the form I 4 = L p) 3 2 k 2r + k 1 + k 2 2 + k 3 2 2 13) )) 2 + 2k 3 2 r p)2 1) It is not et proven that this is a functionall independent integral, as it is possible to construct an infinite number of quartic integrals from combinations of the four eisting quadratic integrals. To test for functional independence, the 6 Jacobian E,I 1,I 2,I 3,I 4 ) i,p i ) can be constructed and shown to be of rank. Thus, the integral given in Eq 1) is the fifth functionall independent isolating integral of motion for the Hamiltonian. 7 16)
In the case when the k i are all ero, this integral reduces to the -component of the Laplace-Runge-Len vector. In fact, if we work through the same derivation but make the cclic permutations, two more integrals are obtained, The are the equivalents of Eq 1) with the coordinates permuted, and reduce to the and components of the Laplace-Runge-Len vector. As is epected the are not functionall independent, and this is easil demonstrated b including them in the Jacobian and noting that it remains of rank. The three components are related through I 4 +I 4 +I 4 = 4EI 1 +I 2 +I 3 +k 1 +k 2 +k 3 )+k 2 17) where I 4 is the integral given b Eq 1) and I 4 and I 4 the cclicl permuted versions. IV. ACTION-ANGLE VARIABLES It is also possible to solve this general problem in action-angle variables. This is worthwhile as it provides insight into the relationship our of our new superintegrable Hamiltonian with the Kepler problem. Following [14], the actions can be shown to be J φ = p φ dφ = 2 2π I2 k 1 ) k 2 J θ = p θ dθ = 2π 2I1 2I 2 ) 2k 3 J r = p r dr = 2π 2I 1 k ) 2E and hence E = 18) 2k 2 π 2 J r +J θ +J φ +2 2π k 1 + k 2 + k 3 ) 2 19) Asepected, the Hamiltoniandepends ontheactions onlthroughthecombination J r +J θ + J φ, impling that the three frequencies of the classical motion are the same. If w r,w θ,w φ ) are the angles conjugate to 19), then we can make a canonical transformation to new action-angle coordinates J 1,J 2,J 3,w 1,w 2,w 3 ), using the generating function F = w φ w θ )J 1 +w θ w r )J 2 +w r J 3 ) The new actions are related to the old via J 1 = J φ, J 2 = J θ +J φ, J 3 = J r +J θ +J φ, 21) 8
and so the Hamiltonian becomes H = E = 2k 2 π 2 J 3 +2 2π k 1 + k 2 + k 3 ) 2 22) and depends on onl one of the new actions. Using Hamilton s equations, we see that the anglew 3 increases linearl withtime, whilst the anglesw 1 andw 2 arethe additionalintegrals of motion. The can be found b eplicit construction of Hamilton s characteristic function S S = p φ dφ+ p θ dθ+ p r dr 23) followed b use of the equations w i = S/ J i see [14] for the equivalent calculation for the Keplerian potential). We find that where w 1 = 1 4π arccos 1 4π arccos w 2 = 1 4π arcsin 1 2π arcsin ) A 2 cos2φ k 1 +k 2 A2 k 1 +k 2 ) 2 4A 2 k 2 ) 2A 2 cot 2 θ B 2 +A 2 +k 2 B2 A 2 k 2 ) 2 4A 2 k 2 ) B 2 cos2θ+a 2 k 3 A2 +B 2 k 3 ) 2 4A 2 B 2 kr 2B 2 r k 2 +4B 2 E ) 24) 2) A = J 1 2 2π + k 1 + k 2 = I 2 26) B = J 2 2 2π + k 1 + k 2 + k 3 = I 1 27) Note that in the case that k 1 = k 2 = k 3 = the Kepler problem), it is usual to introduce the inclination i = arccosa/b) of the orbital plane. Tnen, the angle w 1 reduces to w 1 = 1 φ arcsincotθcoti)) 28) 2π which is the longitude of the ascending node. The second angle w 2 is easiest evaluated in the orbital plane with polar coordinates r, ϕ) and becomes w 2 = 1 2π ϕ ϕ lan) ϕ ϕ peri )) = ω. 29) The first integral is therefore the angular difference between the orbital position and the longitude of the ascending node, the second the angular difference between the orbital position and the periapse. Thus, w 2 reduces to ω, the longitude of the periapse in the Kepler problem. 9
V. SUMMARY AND CONCLUSIONS We have found a new superintegrable Hamiltonian, which is a generaliation of the wellknown Kepler problem. There are five isolating integrals of the motion, namel the energ and generaliations of the componenets of the angular momentum and Laplace-Runge-Len vectors. Intriguingl, some of the integrals of motion are quartic in the momenta and do not arise from separabilit of the Hamilton-Jacobi equation. There are three interesting questions which merit further research. First, it is clear that the N degrees of freedom Hamiltonian H = 1 2 N i=1 p 2 i k r + k i 2 i is also superintegrable. It would be interesting to find the complete set of 2N 1 functionall independent integrals of motion in this case. Second, it is well-known [1, 16] that the additional integrals of motion in the Kepler problem arise from the eistence of the dnamical smmetr group SO4). It would be interesting to understand the group theoretic interpretation of the integrals of motion discussed in this paper. Third, although the reduction technique we used to generate the superintegrable potential in Section II is simple, it can be made to do some more work. For eample, it is also clear that the Hamiltonian [c.f., eq 3)] H = 1 2 p2 1 +p2 2 +p2 3 )+l2 2 +m 2 2 +n 2 2 + k 1 2 + k 2 2 + k 3 2. alwas has independent integrals of motion as well, whose form remains to be established. Perhaps all superintegrable potentials in three degrees of freedom can be viewed as projections of higher dimensional Keplerian or harmonic oscillator motion? 3) [1] Landau, L., & Lifshit, E.M., 1968, Mechanics, Pergammon Press, Oford, p. 11ff [2] Sommerfeld, A., 1923, Atomic Structure and Spectral Lines, Methuen, London, p. 118ff [3] Born, M., 1927, Mechanics of the Atom, G. Bell, London, p. 26ff [4] Fris, J., Mandrosov, V., Smorodinsk, Y. A., Uhlí, M., & Winternit, P. 196, Phsics Letters, 16, 3 [] Makarov, A. A., Smorodinsk, Y. A., Valiev K., & Winternit, P., 1967 Nuovo Cimento 2, 61.
[6] Evans, N. W. 199, Phs. Rev. A, 41, 666 [7] Tempesta, P., Winternit, P., Harnad, J., Miller Jr, W., Pogosan, G., Rodrigue, M.,, Superintegrabilit in Classical and Quantum Sstems, American Mathematical Societ [8] Gravel, S., Winternit, P. 2, J Math Phs, 43, 92 [9] Jauch, J.M., & Hill, E.L., 194, Phs. Rev., 7, 641 [] Adler, M., 1977, Comm. Math. Phs.,, 19 [11] Wojciechowski, S. 1983, Phsics Letters A, 9, 279 [12] Moser, J., 197, Comm. Pure and Applied Maths, 23, 69 [13] Press, W. H., Teukolsk, S. A., Vetterling, W. T., & Flanner, B. P. 2, Numerical recipes in C++ : the art of scientific computing, Cambridge Universit Press [14] Goldstein, H., 198, Classical Mechanics, Addison-Wesle, Reading, p. 613-61 [1] Fock, V., 193, Z. Phs., 98, 14 [16] Bargmann, V., 193, Z. Phs., 99, 76 11