arxiv: v2 [math-ph] 9 May 2018

Fourth-order superintegrable systems separating in Polar Coordinates. II. Standard Potentials Adrian M. Escobar-Ruiz, 1,a J. C. López Vieyra,,b P. Winternitz, 1,c and İ. Yurduşen 1,,d 1 Centre de recherches mathématiques and Département de mathématiques arxiv:1804.05751v [math-ph 9 May 018 et de statistique, Université de Montreal, C.P. 618, succ. Centre-ville, Montréal QC HC J7, Canada Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Apartado Postal 70-54, 04510 México, D.F., Mexico Department of Mathematics, Hacettepe University, 06800 Beytepe, Ankara, Turkey Dated: 11 May 018 Superintegrable Hamiltonian systems in a two-dimensional Euclidean space are considered. We present all real standard potentials that allow separation of variables in polar coordinates and admit an independent fourth-order integral of motion. The general form of the potentials satisfies a linear ODE. In the classical case, the standard potentials coincide with the Tremblay-Turbiner-Winternitz TTW or Post- Winternitz PW models. In the quantum case new superintegrable systems are obtained, in addition to the TTW and PW ones. Their classical limit is free motion. Keywords: Superintegrability, separation of variables a Electronic mail: escobarr@crm.umontreal.ca b Electronic mail: vieyra@nucleares.unam.mx c Electronic mail: wintern@crm.umontreal.ca d Electronic mail: yurdusen@hacettepe.edu.tr 1

I. INTRODUCTION The purpose of this paper is to obtain and classify all classical and quantum superintegrable Hamiltonians that allow separation of variables in polar coordinates and admit a fourth-order integral of motion Y. We focus on standard potentials, namely those that satisfy not trivially a linear ODE. This linear ODE is nothing but the compatibility condition for the existence of the integral Y. Superintegrable systems separating in polar coordinates, allowing fourth-order integrals of motion and involving exotic potentials were studied in an earlier article 1. Roughly speaking, a Hamiltonian with n degrees of freedom is called integrable if it allows n independent well defined integrals of motion in involution. It is minimally superintegrable if it allows n+1 such integrals, and maximally superintegrable if it admits n 1 integrals where only subsets of n integrals among them can be in involution. The best known superintegrable systems are the harmonic oscillator with its sun + 1 algebra of integrals, and the Kepler-Coulomb system with its on + 1 algebra when restricted to fixed bound state energy values. A more recent review article provides precise and detailed definitions, general settings and motivation for studying superintegrable systems 9. A systematic search for superintegrable classical and quantum systems in E and E established an interesting connection between second-order superintegrability and multiseparability in the Schrödinger or Hamilton-Jacobi equation 5,6,15,16,1,,4. An extensive literature exists on second-order superintegrability in spaces of, and n dimensions, Riemannian and pseudo-riemannian, real or complex, 4 8,. The systematic study of higher order integrability is more recent. Pioneering work is due to Drach 10,11. For more recent work see,7 9,14,19,6 8,41,44 47,49. The present article is a contribution to a series 1,1,18,19,5,8,40,4,48 devoted to superintegrable systems in E with one integral of order n and one of order n. In this article we restrict ourselves to the space E. The Hamiltonian has the form H = 1 p x +p y +Vx,y, 1 Vx,y is a scalar potential and in classical mechanics p x and p y are the momenta conjugate to the Cartesian coordinates x and y. In quantum mechanics they are the corresponding

operators p x = i, p x y = i. In polar coordinates x,y rcosθ, rsinθ, the y classical Hamiltonian reads H = 1 p r + p θ r + Vr,θ, wherep r andp θ aretheassociatedcanonicalmomenta. Thecorrespondingquantumoperator takes the form H = r + 1 r r + 1 r θ + Vr,θ. In paper 1 we found all the exotic potentials, i.e. those which do not satisfy any linear differential equation. Here we solve the complementary problem, namely we find all the standard potentials which do satisfy a linear ODE. In all equations we keep the Planck constant explicitly. Classical standard potentials will be obtained in the limit 0. Before going into the question of standard superintegrable potentials with a fourth-order integral of motion, let us recall that two families of standard superintegrable potentials in E are already known. Both allow separation of variables in polar coordinates and hence admit the second-order integral X = p θ + Sθ, p θ = L z = i θ. 4 They also admit a further polynomial integral of arbitrary order N. One of them is the TTW potential V TTW = br + 1 [ α r cos kθ + β sin, 5 kθ where k = m/n and m and n are two integers with no common divisors 46,47. The other is the PW potential 4 where again k = m/n is rational. [ V PW = a r + 1 µ r cos kθ + ν sin kθ, 6 In the TTW case 5 the lowest order of an additional integral Z was shown to be 17,9,0,44 both in classical and quantum mechanics. N = m+n 1, 7 The TTW and PW potentials are related by coupling constant metamorphosis 4,,,4. This transformation takes integrable and superintegrable systems into other systems that

are also integrable or superintegrable, respectively. However, in general the order of the integrals of motion is not preserved, nor is their polynomial character. Special conditions must be satisfied for polynomial integrals to be transformed into polynomials of the same order. It has been conjectured and verified for k = 1 and that these conditions are satisfied for the TTW and PW systems 4. Thus we may expect that in any study of superintegrable systems separating in polar coordinates in E with an additional integral of order N with N < the TTW and PW systems will appear as standard potentials. A further comment is in order here and should be taken into account in any classification of superintegrable systems. If a certain potential Vx,y appears as superintegrable in E with two integrals X 1 and X that are polynomials of order m and n, respectively, then the same potential will reappear at higher orders. Indeed X 1 X +X X 1 will be an integral of order m + n for the same Hamiltonian. The commutator [X 1,X or Poisson bracket {X 1,X } P.B. willbeapolynomialintegraloforderm+n 1. Such trivially superintegrable systems should be weeded out of lists of new superintegrable systems. For fourth order superintegrable potentials, in addition to the Hamiltonian H, we have two more conserved quantities. A second-order integral 4, associated with the separation of variables in polar coordinates for which the potential takes the form and a fourth-order one Vr, θ = Rr + Sθ r, 8 Y = A 1 {L z, p x} + A {L z, p y} + A {L z, p x p x +p y } + A 4{L z, p y p x +p y } 9 + B 1 p 4 x p4 y + B p x p y p x +p y + B {L z, p x p y } + B 4{L z, p xp y } + C 1 {L z, p xp y p y} + C {L z, p x p yp x } + D 1 p 4 x +p 4 y 6p xp y + 4D p x p y p x p y + {g 1 x,y, p x} + {g x,y, p x p y } + {g x,y, p y} + g 4 x,y, where the A i, B i, C i, D i are real constants. The bracket {, } denotes an anticommutator and Rr, Sθ, g 1,,,4 x,y are real functions such that [H, Y = [H, X = 0. 10 Under rotations around the z-axis, each of the six pairs of parameters in 9 A 1,A, A,A 4, B 1,B, B,B 4, C 1,C, D 1,D, 4

forms a doublet. In general, Y may contain the O singlets but by linear combinations of the form Y + u 1 X + u H + u XH, where the u i are constants, we eliminated all such trivial terms. Under rotationsthroughtheangleθthedoublets A i, B i, C i, D i rotatethrough θ, θ, θ, 4θ, respectively. We have [Y,X = C 0 where C is in general a 5 th -order linear operator. We thus obtain a finitely generated polynomial algebra of integrals of motion. We are looking for fourthorder superintegrable systems, so at least one of the constants A i, B i, C i, D i is different from zero. The operator Y is the most general polynomial expression for a fourth-order Hermitian operator of the required form. The commutator [H, Y contains derivatives of order up to three. The operator Y in 9 is given in Cartesian coordinates for brevity. Putting p x = i cosθ r sinθ r θ, p y = i sinθ r + cosθ r θ, we obtain the corresponding expression in polar coordinates. Explicitly, the leading terms of the integral Y are Y = 4 B 1 cosθ+b sinθ+d 1 cos4θ+d sin4θ 4 r + 1 r 4 [D sin4θ +D 1 cos4θ r A 1 r +A 4 sinθ B +B 4 r sinθ+r A r +A 4 cosθ B 1 +B r cosθ rc 1 cosθ C sinθ [D r 1 cos4θ+d sin4θ r B cosθ +B 4 sinθ +ra sinθ A 4 cosθ C 1 cosθ C sinθ r θ [ B 1 sinθ B cosθ+d 1 sin4θ D cos4θ rc 1 sinθ+ C cosθ r +A cosθ+a 4 sinθ r θ r [ B 1 sinθ B cosθ D 1 sin4θ+d cos4θ ra cosθ+ A 4 sinθ C 1 sinθ+c cosθ +r B sinθ B 4 cosθ r A 1 cosθ+a sinθ r θ +lower order terms. For convenience let us introduce the functions G 1 r, θ = g 1 cos θ+g sin θ+g cosθsinθ, 5 4 θ 11

G r, θ = g 1 sin θ+g cos θ g cosθsinθ r, G r, θ = g 1 sinθ g sinθ g cosθ r, G 4 r, θ = g 4. 1 Then, we can write {g 1 x,y, p x} + {g x,y, p x p y } + {g x,y, p y} + g 4 x,y = {G 1 r,θ, r } + {G r,θ, r θ } + {G r,θ, θ } + G 4r,θ. 1 II. DETERMINING EQUATIONS FOR A FOURTH-ORDER INTEGRAL A. Commutator [H, Y Since Y is a fourth-order operator the commutator [H, Y would be a fifth-order one, i.e. we have [H, Y = N+1 k+l=0 Z k,l r, θ k+l r k θ l, N = 4 14 and we require Z k,l = 0 for all k and l. The terms of order k +l = 5 and k +l = 4 vanish automatically. It was already shown that for arbitrary N the terms of order N +1 and N vanish 4. Moreover, only terms with k+l of the opposite parity than N provide independent determining equations Z k,l = 0. Those with the same parity provide equations that are differential consequences of the first ones. Finally, we obtain [H,Y = A rrr r + A rrθ +C r r + C θ r θ + A rθθ + A r θ θθθ θ θ = 0, 15 here the coefficients A rrr,a rrθ,..., are real functions of r and θ. The determining equations are obtained by setting all coefficients in 15 equal to zero. We thus obtain a total of 6 6

PDEs for the functions G 1, G, G, G 4, R and S. Four equations are independent of 1 G 0,1 r + 1 r G G 1,0 1 = F 1 R F 1 S + F, 16 r r S = F R F r S + F 4 r S, 17 1 r G0,1 1 +G 1,0 = F R 6F S + F 5, 18 r r S r G 1 +G 1,0 + 1 r G0,1 = F 5 R F 5 S + F r r S. 19 They are the same for classical or quantum systems, and two equations contain quantum corrections proportional to 1 G r S +G 1 R r S 1 G 4 1,0 [ + 0, G1 G 0,1 + G 0, G 1,0 1 r 4r 4r r G 1,0 + 5G 1,1 4r + G 1 1, + 1 r G 1, + G,0 1 1 r rg,0 + 4 G,1,0 +G 1 [ = 6F6 r 4F 1 F 11 6F10 S + 4F 7 + F 1 S 4 r 5 r r 4 r r F8 + r F 5 S + F r r S +F 11 R +F 6 R +F 1 R, 0 G R r S + G 1 r S G 0,1 4 [ r + G 4r G 0,1 + 5G 0, + G 0, G 1,0 + G 1,1 1 r 4r r 4r r + G 1, + G,0 + G 1,1 + 1 4r r r G,1 + 1 4 G,0 [ = 6F7 4F F 1 6F5 S + 4F 8 + F 1 S r 4 r 5 r r 4 r r F9 + r 6F S + F 4 r r S +F 1 R +F 7 R +F R. 1 Thefunctions F 1...F 1 arecompletely determined bythe constants A i,b i,c i,d i figuring in the leading part of the integral Y. They are presented in the Appendix. The system of equations 16-1 is the main object of study of the present paper. 7

B. Linear compatibility condition The compatibility equation for 16-19 is a linear equation 0 = 1 [ 16 r A cosθ +A 4 sinθ 9C 1 sinθ +C cosθ S [ 0 A 4 cosθ A sinθ +9C sinθ C 1 cosθ S [ +80 C 1 sinθ + C cosθ S [ +5 A sinθ A 4 cosθ+c sinθ C 1 cosθ S [ A cosθ +A 4 sinθ +C 1 sinθ +C cosθ S 4 + 1 [ 96 r B cosθ B 1 sinθ +8D 1 sin4θ D cos4θ S [ +40 B 1 cosθ +B sinθ 0D sin4θ +D 1 cos4θ S [ +0 B cosθ B 1 sinθ+14d cos4θ D 1 sin4θ [ +10 B sinθ +B 1 cosθ+4d 1 cos4θ +D sin4θ S [ B cosθ B 1 sinθ +D cos4θ D 1 sin4θ S 4 + [4r A 1 cosθ+a sinθ 1r B 4 cosθ B sinθ r A 4 sinθ +A cosθ+18c 1 sinθ +C cosθ +B cosθ B 1 sinθ+0d cos4θ D 1 sin4θ R + [6r 4 A 1 cosθ+a sinθ 1r B sinθ B 4 cosθ +r A cosθ+a 4 sinθ+18c 1 sinθ +C cosθ S +r B 1 sinθ B cosθ 0D cos4θ D 1 sin4θ R + [1r 5 A 1 cosθ+a sinθ +14r 4 B 4 cosθ B sinθ +r 5A cosθ+a 4 sinθ C 1 sinθ +C cosθ +1r D cos4θ D 1 sin4θ R + [r 6 A 1 cosθ+a sinθ +r 5 B 4 cosθ B sinθ +r 4 A 4 sinθ+a cosθ C 1 sinθ +C cosθ r B 1 sinθ B cosθ +D cos4θ D 1 sin4θ R 4, for the functions Rr and Sθ in the potential. It is the same in classical and quantum mechanics as are equations 16,..., 19. 8

We will henceforth define exotic potentials as those for which the compatibility condition is satisfied trivially. For N = 4 this means all the coefficients in the integral Y of 9 except A 1,A,B and B 4 vanish. Standard potentials must satisfy the linear PDE nontrivially. Hence at least one of the constants A i,b i,c i or D i in must be nonzero. Differentiating with respect to r we eliminate Sθ for the equation. Following paper 1 we obtain 8 ODEs for Rr by setting the coefficients of coskθ and sinkθ for k = 1,,,4 equal to zero. These together with the determining equations imply that the only possible forms of Rr are Rr = 0, Rr = a r, Rr = br. The result is actually a consequence of Bertrand s theorem and is valid for integrals of motion Y of any order 1. Let us now consider each of the cases listed in separately. III. NON CONFINING POTENTIAL Vr,θ = Sθ r Here we address the case of a non confining potential, namely Rr = 0. It turns out that this case provides the most general form of the function Sθ. The same function Sθ will reoccur for Rr 0. The equations 16, 18 and 19 corresponding to the determining equations A rrr = A rrθ = A rθθ = 0, respectively, take the following form G 1,0 1 r,θr A cosθ+a 4 sinθ +C 1 sinθ+c cosθ r B 1 sinθ B cosθ+ D 1 sin4θ D cos4θ S +4 B 1 cosθ+ B sinθ+ D 1 cos4θ+ D sin4θ S = 0, 4 G 1,0 r,θr 4 +G 0,1 1 r,θr B cosθ+ B 4 sinθ A sinθ A 4 cosθ C 1 cosθ C sinθ D 1 cos4θ+ D sin4θ S +6 A cosθ+a 4 sinθ+c 1 sinθ+c cosθ r r B 1 sinθ B cosθ+d 1 sin4θ D cos4θ r S = 0, 5 9

G 1,0 r,θr 5 +G 0,1 r,θr +G 1 r,θr A 1 cosθ+ A sinθr +B 4 cosθ B sinθr +A cosθ+ A 4 sinθ C 1 sinθ+ C cosθr B 1 sinθ+ B cosθ+ D 1 sin4θ D cos4θ +4 B cosθ+ B 4 sinθr S 6 +A 4 cosθ A sinθ+ C 1 cosθ C sinθr D 1 cos4θ+ D sin4θ S = 0. The above equations can be integrated to find the functions G 1 r,θ, G r,θ and G r,θ in terms of Sθ. They are given by the following formulas: 1 G 1 r,θ = B r 1 sin θ B cos θ+d 1 sin 4θ D cos 4θ 1 A cos θ+a 4 sin θ +C 1 sin θ+c cos θ S r + B r 1 cos θ +B sin θ +D 1 cos 4θ +D sin 4θ S +γ 1 θ, 7 G r,θ = r B cos θ +B 4 sin θ S + 1 [ A r cos θ +A 4 sin θ +C 1 sin θ +C cos θ S + A sin θ A 4 cos θ C 1 cos θ C sin θ S 1 A cos θ +A 4 sin θ+ C 1 sin θ +C cos θ S + 1 [ 10 B r 1 sin θ B cos θ+d 1 sin 4θ D cos 4θ S + B 1 cos θ +B sin θ +4D 1 cos 4θ +D sin 4θ S + 1 B 1 sin θ B cos θ +D 1 sin 4θ D cos 4θ S 6 + 1 r γ 1 θ+γ θ, 10 8

G r,θ = r A 1cosθ+ A sinθs + 1 [B r cos θ +B 4 sin θs +5B sinθ B 4 cosθs B cosθ+b 4 sinθs + 1 [ 7 A r sinθ A 4 cosθ C 1 cosθ C sinθ S 1 A cosθ + A 4 sinθ 47C 1 sinθ+ C cosθ S 6 + A sinθ A 4 cosθ C 1 cosθ C sinθ S 1 A cosθ + A 4 sinθ+c 1 sinθ+ C cosθ S 6 + 1 [ B r 4 1 cosθ+ B sinθ+1d 1 cos4θ+ D sin4θ S 1 B 1 sinθ B cosθ+0d 1 sin4θ D cos4θ + 1 B 1 cosθ+ B sinθ+4d 1 cos4θ+ D sin4θ B 1 sin θ B cosθ+d 1 sin4θ D cos4θ + 1 4 + 1 r γ θ+ 1 r γ 1θ+ 1 r γ 1 θ+γ θ, where γ 1 θ,γ θ and γ θ are functions of θ still to be determined. S S S The determining equation A θθθ = 0 17 has not been used yet. We now substitute the expressions 7-9 into 17 in order to obtain further information on functions of θ, namely Sθ andγ k θ, k = 1,,. The variablerstill figures intheobtained result but only explicitly as r k, k = 0,1,,,4 no unknown functions of r remain. Thus, the coefficient of r k in A θθθ = 0 17 must vanish for each k and we obtain a system of 5 equations for the 4 remaining functions of θ. They are 9 0 = 4 A 1 cos θ+ A sinθ S +14 6 A 1 cosθ+ A sinθ 0 = γ, 0 A 1 sin θ A cosθ S S +γ +γ, 1 0 = 16 B sinθ B 4 cosθ S +8 B cosθ+ B 4 sinθ +14 B sinθ B 4 cosθ S B cosθ+ B 4 sinθ S S +4γ 1 +γ 1, 11

0 = 96 B cos θ B 1 sin θ+8d 1 sin 4θ D cos 4θ S + 40 B 1 cos θ+ B sin θ 0D 1 cos 4θ+ D sin 4θ S + 0 B cos θ B 1 sin θ+14d cos 4θ D 1 sin 4θ S + 10 B 1 cos θ+ B sin θ+4d 1 cos 4θ+ D sin 4θ S + B 1 sin θ B cos θ+d 1 sin 4θ D cos 4θ S 4, 0 = 16 A cos θ+a 4 sin θ 9C 1 sin θ+ C cos θ S + 0 A sin θ A 4 cos θ+9c 1 cos θ C sin θ S + 80 C 1 sin θ+ C cos θ S + 5 A sin θ A 4 cos θ C 1 cos θ C sin θ S A cos θ+a 4 sin θ+c 1 sin θ+c cos θ S 4. 4 Let us now discuss these equations. 1. The equation 0 determines γ = constant.. Equations 1 and must be used to determine γ 1 θ and γ θ in terms of the function Sθ figuring in the potential.. Ofthe1 doubletcoefficients figuringintheleadingpartoftheintegraly,fourfigure onlyin, namelyb 1, B, D 1, D andfouronlyin4, namelya, A 4, C 1, C. The remaining ones A 1, A and B, B 4 appear only in 1 and, respectively. By definition, the linear compatibility condition must be satisfied trivially for exotic potentials at Rr = 0. Hence, for these potentials the only possible constants figuring in Y are A 1, A and B, B 4. This case was analyzed conclusively and in detail in our previous article 1. Let us now turn to the problem of finding all standard superintegrable potentials. To do this we proceed as follows I. Assume that at least one of the constants B 1, B, D 1, D is not vanishing and solve for Sθ. The general solution for Sθ will depend on these constants and on additional integration constants. 1

II. We proceed similarly with equation 4, i.e. solve it assuming that at least one of the A, A 4, C 1, C is non-zero. In both cases it is still necessary to solve the determining equations C r = 0, C θ = 0, 5 with C r and C θ figuring in 15, see 0-1. The determining equation C r = 0 defines G 4 r,θ, C θ = 0 will provide further information and make it possible to determine Sθ completely and hence also the integrals X and Y. III. It is also possible for both of the above sets of constants to contain nonzero constants. The corresponding potentials will then allow more than one fourth-order integral. However, at most integralsh, X and Y can be polynomially independent. In any case, the function S III θ obtained in this case will be a special case of both S I and S II and we shall not pursue this issue further. A. General form of the angular part Sθ Case I. B 1, B, D 1, D 0, 0, 0, 0 We concentrate in equation. We first define T I θ by putting T I θ S Iθ. 6 The general solution of the fifth-order linear ODE for T I obtained from is T I θ = 1 M I s 1 +s sinθ+s cosθ+s 4 sin4θ+s 5 cos4θ, 7 where the s i are integration constants and M I = B 1 sinθ B cosθ+d 1 sin4θ D cos4θ. 1

Using 6 we obtain S I θ as S I θ = 1 M I s 1 [B 1 cosθ+b sinθ+4d 1 cos4θ+d sin4θ + s [B 1 +D 1 cos6θ+d sin6θ+d 1 cosθ+d sinθ + s [B D cos6θ+d 1 sin6θ+d cosθ D 1 sinθ 8 + s 4 [8D 1 B 1 cos6θ B sin6θ+b 1 cosθ B sinθ + s 5 [8D +B cos6θ B 1 sin6θ+b cosθ+b 1 sinθ. There seem tobe toomanyintegration constants in8, since is afourth-order ODE for Sθ. Indeed the five solutions corresponding to s 1,...,s 5, are not independent and one can be eliminated in terms of the others. We shall not go into this here since constraints on the constants s i will be imposed by the last still unsolved determining equation, namely C θ = 0 and these will always leave less than 5 free constants. Case II A, A 4, C 1, C 0, 0, 0, 0 Similarly, the solution to the second linear equation 4 is T II θ = 1 M II s 1 +s sinθ +s cosθ +s 4 sinθ+s 5 cosθ, 9 where M II = A cosθ+a 4 sinθ+c 1 sinθ+c cosθ. Consequently S II θ = 1 M II s 1 [A sinθ A 4 cosθ+c sinθ C 1 cosθ + s [A C cos4θ C 1 sin4θ+c cosθ+c 1 sinθ s [A 4 +C 1 cos4θ C sin4θ+c 1 cosθ C sinθ 40 + s 4 [C +A cos4θ+a 4 sin4θ+a cosθ A 4 sinθ s 5 [C 1 A 4 cos4θ+a sin4θ+a 4 cosθ+a sinθ. Again the five solutions corresponding to the 5 constants s i, are not linearly independent. As mentioned above, the determining equation C r = 0 0 defines the function G 4 r, θ in both cases. The last equation C θ = 0 1 constrains the coefficients s i in both cases. We again treat cases I and II, separately. 14

Case I B 1, B, D 1, D 0, 0, 0, 0 In this case, from the equation C r = 0 0 we obtain the function G 4 r, θ 48r 4 G 4 r,θ = 4 SS B 1 sinθ B cosθ+d 1 sin4θ D cos4θ +48S B 1 cosθ+b sinθ +D 1 cos4θ+d sin4θ +S [ S B cosθ B 1 sinθ D 1 sin4θ D cos4θ 6S B 1 cosθ +B sinθ+4d 1 cos4θ+d sin4θ + [07S D 1 cos4θ+d sin4θ +64S B 1 sinθ B cosθ+6d 1 sin4θ D cos4θ 41 96S B 1 cosθ+b sinθ+0d sin4θ+d 1 cos4θ +S 5B cosθ B 1 sinθ+440d cos4θ D 1 sin4θ +1S 4 B 1 cosθ+b sinθ+4d sin4θ+d 1 cos4θ +S 5 B 1 sinθ B cosθ+d 1 sin4θ D cos4θ. 15

The determining equation C θ = 0 reads 4 56S B 1 sinθ B cosθ+d 1 sin4θ D cos4θ 40SS B sinθ+b 1 cosθ+4d 1 cos4θ+d sin4θ +60S B cosθ B 1 sinθ 8D 1 sin4θ D cos4θ +40SS B cosθ B 1 sinθ D 1 sin4θ D cos4θ +0S S B 1 cosθ+b sinθ+4d 1 cos4θ+d sin4θ +S B 1 sinθ B cosθ+d 1 sin4θ D cos4θ +S S B 1 sinθ B cosθ+d 1 sin4θ D cos4θ + [84S B cosθ B 1 sinθ 88D 1 sin4θ D cos4θ S 1B 1 cosθ+b sinθ 1676D 1 cos4θ+d sin4θ +16S 5B cosθ B 1 sinθ+114d 1 sin4θ D cos4θ 4 80S B 1 cosθ +B sinθ+16d 1 cos4θ+d sin4θ 16S 4 11B 1 sinθ B cosθ+118d 1 sin4θ D cos4θ +4S 5 B 1 cosθ+b sinθ +4D 1 cos4θ+d sin4θ +S 6 B 1 sinθ B cosθ+d 1 sin4θ D cos4θ = 0. 16

First, by an appropriate rotation we can always set B 1 = 0. Then, we substitute the function 8 in the remaining equation 4. As a result we obtain algebraic equations for the s i i = 1,,,4,5 in 8 by setting the coefficients of coskθ and sinkθ for k = 0,,4,...,18 equal to zero. We obtain the following solutions s l 1 = B D 1 D B +40D 1 q l D 1 q l [ B D 8D 1 B B D 1 +D +8D 1 4D 1 +D 6 4D1 6 B 8D1 D + 8B D 5 1 6D1 D B 4 D 1 +D D 1 +64 D 1 +D D 5 1 4 q l B +8D1 q l, s l = q l, s l = q l D 1 D 6 B 8D 1 + B D 1 +4D D 1 s l 4 = q l D 6 B 8D 1 [ D B 8D1 B D 1 +D +8D 1 D 1 +D 6 B D1 5 4 q l + D1 D 7B +8D1 ql + B D1 q l [ 4B D D 1 +D B 8D1 6 + D 1 B D 1 +6D 16D 1 D 1 +D 4 q l + 8B D1D ql + D1q l,, s 5 = 0, where q l, l = 1,,,4, are the four roots of the quartic equation 4 D1 4 q4 + 4B D D1 q +B D 4 1 +6D D 1 16D 4 1 D 1 +D 4 q B D 1 8D 1 B D D 1 +D 6 q + B 1 8D D D 1 +D 8 = 0. 44 The discriminant Γ of the quartic equation 44 is given by Γ = 56 4 D 4 1 D B 8D1 [ B 4 60D 48D1 +768B D 1 +D +B 6 4096 D1 +D. This discriminant iszero if andonlyif at least two rootsof44 areequal. If thediscriminant is negative there are two real roots and two complex conjugate roots. If it is positive the roots are either all real or all non-real. From a physical point of view only the real solutions are admitted. Let us analyze the zeros of the discriminant 45. They will correspond to the TTW model 46,47. For Γ 0, by substituting 4 into 8 we obtain an angular component S I θ proportional to with no classical analog, it cannot be transformed or reduced to that of the TTW model. 17 45

1. Case = 0 For = 0, the discriminant 45 vanishes. The corresponding coefficients take the values s 1 = s 1, s = 0, s = 0, s 4 = s 4, s 5 = s 5, which yields the potential S I θ = 4D 1cos4θ+ D sin4θs 1 +4D 1 s 5 + D s 4 D1 D cos 8θ+D 1 D sin 8θ D 1 +D. 46 For D = 0, the angular part of the potential becomes S I θ = s 1 cos 4θ+ s 5 cos 8θ 1, 47 where s 1 = 4s 1 /D 1, s 5 = 4s 5 /D 1. This potential 47 coincides with that of the TTW Hamiltonian 46,47. The angular component of the TTW potential is S TTW θ = αk cos kθ + βk sin kθ = 4k α βcoskθ 4k α + β cos4kθ 1. 48 Then for k = and α = s 1 s 5 /,β = s 1 + s 5 /, the potential 47 actually corresponds to a rotated TTW model with no radial component Rr = 0 which is a superintegrable system in both the classical and quantum cases.. Case D 1 = 0 For D 1 = 0, the discriminant 45 vanishes as well. The corresponding coefficients vanish, s 1 = s = s = s 4 = s 5 = 0, which gives the trivial solution S I θ = 0.. Case D = 0 For D = 0, the discriminant 45 vanishes again. The corresponding coefficients are given by s 1 = s 1, s = 0, s = B s 4 8D 1s 1 +s 4 B D 1, s 4 = s 4, s 5 = 0, 18

thus Sθ = B s 4 +D 1 s 1 sinθ+d 1 s 4 sinθ B D 1 1+cos4θ This solution corresponds to the angular component of the TTW model 48 with k = 1. The zeros of the discriminant 45 correspond to a TTW model. Case II: A, A 4, C 1, C 0, 0, 0, 0 Now, we proceed to the case II. First, from the determining equation C r = 0 0 we obtain the function G 4 r, θ 48r G 4 r,θ = 16S [S A 4 cosθ A sinθ C sinθ C 1 cosθ. + S 14S A cosθ +A 4 sinθ +C 1 sinθ+c cosθ + [ 16S A 4 cosθ A sinθ+15c sinθ C 1 cosθ 4S 7A cosθ+a 4 sinθ 17C 1 sinθ+c cosθ +4S A sinθ A 4 cosθ +9C 1 cosθ C sinθ +S 9A cosθ+a 4 sinθ+11c 1 sinθ+c cosθ +6S 4 A sinθ A 4 cosθ+c sinθ C 1 cosθ S 5 A cosθ+a 4 sinθ+c 1 sinθ+c cosθ, 49 19

and correspondingly C θ = 0 1 takes the form [ 4 6S A cosθ+a 4 sinθ+c 1 sinθ+c cosθ +15S A cosθ +A 4 sinθ +9C 1 sinθ +C cosθ +0SS A cosθ +A 4 sinθ +C 1 sinθ+c cosθ S A cosθ+a 4 sinθ+c 1 sinθ+c cosθ +S 0S A 4 cosθ A sinθ+c 1 cosθ C sinθ +5S A sinθ A 4 cosθ+c sinθ C 1 cosθ S A cosθ +A 4 sinθ +C 1 sinθ+c cosθ + [ 96S A cosθ +A 4 sinθ+9c 1 sinθ+c cosθ +4S A sinθ A 4 cosθ+0c sinθ C 1 cosθ +8S 1A cosθ+a 4 sinθ 719C 1 sinθ+c cosθ +0S A sinθ A 4 cosθ+79c 1 cosθ C sinθ +8S 4 A cosθ+a 4 sinθ+67c 1 sinθ +C cosθ +1S 5 A sinθ A 4 cosθ+c sinθ C 1 cosθ S 6 A cosθ+a 4 sinθ+c 1 sinθ+c cosθ = 0. 50 0

Now we focus on the potential S II θ 40. By an appropriate rotation we can always set A = 0. Then, we substitute the expression 40 into the remaining determining equation 4. As a result we obtain algebraic equations for the coefficients s i i = 1,,,4,5 by setting the coefficients of coskθ and sinkθ for k = 1,,,...,16 equal to zero. In general, we obtain the following three solutions s 1 = 0, s = s, s l = q l, s l 4 = 4 C 1 +C 1 A 4 +C q l +C 1 C s C 1 +A 4 q l q l A 4 C 4, s l 5 = 4 4C 1 +C 1A 4 +6C A 4 q l +C 1 +A 4 C s C 1 +A 4 q l q l A 4 C 1 +A 4 4, where q l, l = 1,,, is the solution to the cubic equation 51 q +A 4 q 4 q C1 +A 4C 1 +C + C 1 +A 4 C1 +C C 1 +A 4 C 6 = 0. 5 Again, substituting 51 into 40 we obtain the corresponding angular component S II θ of the potential. Notice that there exist potentials S II θ proportional to that have no classical analog. The discriminant Ω of the cubic equation 5 is given by Ω = 4 1 C [ A 4 4 + 8A 4 C 1 + 18A 4 C 1 +C 7C 1 +C. 5 The discriminant 5 is zero if and only if at least two roots are equal. It is positive if the roots are all distinct real numbers, and negative if there exist one real and two complex conjugate roots. For = 0, the discriminant 5 vanishes. The corresponding coefficients take the values s 1 = s 1, s = 0, s = 0, s 4 = s 4, s 5 = s 5, which yields the potential S II θ = 6s 5 +s 1 cosθ C 1 1 cos6θ, 54 1

A 4 = C = 0 it corresponds to that of the TTW potential with k =. Similarly we can show that for C = 0, thus Ω = 0, the TTW angular potentials with k = 1 and k = 1 occur. IV. CONFINING POTENTIALS Deformed Coulomb Potential Here we address the case of the deformed Coulomb potential where a 0 is a real constant. Vr, θ = a r + Sθ r, 55 Similarly to the non confining case Rr = 0, from 16, 18 and 19 corresponding to the determining equations A rrr = A rrθ = A rθθ = 0, respectively, we determine the functions G 1 r,θ, G r,θ, and G r,θ up to arbitrary additive functions γ 1 θ, γ θ and γ θ, respectively. Then, we substitute such functions G 1,, into the determining equation A θθθ = 0 17. By doing so, such equation yields again a set of five linear ODE 1. Two fourth-order ODE for the angular component Sθ. The first ODE contains the doublets B 1,B,D 1,D, while the second ones depends on the doublets A,A 4,C 1,C and the parameter a figuring in the Coulomb term.. Three linear ODE that define the functions γ 1 θ, γ θ and γ θ, respectively. The equation for γ 1 depends on the doublet B,B 4 and the parameter a. For the function γ the doublet A 1,A and the parameter a are involved and the equation for γ contains only the parameter a. For a 0 the most general form for the angular component Sθ is given by the Case II 40 for which at least one of the constants C 1, C, A, A 4 is not vanishing. Therefore, the function Sθ is not new with respect to the non-confining case Rr = 0 and it was already described in detail in the Section III.

The corresponding functions G 1,, are given by G 1 r,θ = G 0 1 r,θ, [ G r,θ = G 0 r,θ + 5a A4 cosθ A sinθ+c 1 cosθ C sinθ, r [ G r,θ = G 0 A cosθ+a 4 sinθ + C cosθ+c 1 sinθ r,θ + a, r where the G 0 i r,θ, i = 1,,, coincide with those of the non-confining potential given by 7, 8 and 9, respectively, putting B 1 = B = D 1 = D = 0 and γ 1 = γ = γ = 0. From the determining equation C r = 0 given in 0 we calculate the corresponding function G 4 r,θ [ G 4 r,θ = G 0 1 1 4 r,θ +a r 4 7A 4 cosθ A sinθ+99c sinθ C 1 cosθ +5S θa 4 sinθ+a cosθ+c 1 sinθ+c cosθ, here G 0 4 r,θ is given in 49. The angular component Sθ is given by 40 with the coefficients s i 51 presented in the Case II. 56 57 Deformed Harmonic Oscillator Potential Here we address the case of the confining potential with b 0 a real constant. Vr, θ = br + Sθ r, 58 Following the same strategy as for the previous case, we obtain the following functions G 1 r,θ = G 0 1 r,θ + br B 1 cosθ+b sinθ+d 1 cos4θ+d sin4θ, G r,θ = G 0 r,θ, G r,θ = G 0 r,θ + brb cosθ B 1 sinθ+d cos4θ D 1 sin4θ, G 4 r,θ = G 0 4 r,θ, where the G 0 i r,θ, i = 1,,, coincide with those of the non-confining potential Rr = 0, given in 7, 8 and 9, respectively, putting A = A 4 = C 1 = C = 0 and γ 1 = γ = γ = 0. The function G 0 4 r,θ is given in 41. For b 0, the angular component Sθ is given by 8 with the coefficients s i 4 presented in the Case I for which at least one of the constants B 1, B, D 1, D is not vanishing. 59

V. CONCLUSIONS We considered superintegrable systems in a two-dimensional Euclidean space. Classical and quantum fourth-order superintegrable standard potentials separating in polar coordinates were derived. We can summarize the main results via the following theorems Theorem 1. In classical mechanics, the standard superintegrable confining systems, ab 0, correspond to the TTW potential V TTW r, θ = br + 1 [ α r cos θ + β sin, 60 θ α, β real constants, and the PW potential V PW r, θ = a r + 1 [ µ r cos θ + ν sin θ, 61 µ, ν real constants. The corresponding leading terms of the integral Y in 9 are proportional to p 4 x +p 4 y 6p xp y, and {L z, p xp y p y}, respectively. These terms are independent of a and b. Theorem. In quantum mechanics, the new confining superintegrable systems correspond to Vr, θ = br + 1 r S Iθ, and Vr, θ = a r + 1 r S IIθ, where S I θ is given by 8 with 4, and S II θ takes the form 40 with 51. In general, the functions S I and S II are proportional to and cannot be reduced to a TTW or PW potential. The corresponding leading terms of the integral Y in 9 are B p x p y p x +p y + D 1p 4 x +p4 y 6p x p y + 4D p x p y p x p y, B,D 1,D 0,0,0 and A 4 {L z, p y p x +p y} + C 1 {L z, p xp y p y} + C {L z, p x p yp x }, 4

A 4,C 1,C 0,0,0 respectively. We emphasize that these are pure quantum potentials. As particular cases, both potentials V TTW and V PW appear in the quantum case as well. Work is currently in progress on a general and unified study of N th -order exotic and standard potentials separating in polar coordinates 1. Within this study 1 the TTW and PW models correspond to standard classical potentials. For some cases we plan to present the polynomial algebra generated by the integrals of motion and to use it to calculate the energy spectrum and the wave functions in the quantum case. VI. ACKNOWLEDGMENTS The research of PW was partially supported by a research grant from NSERC of Canada. JCLV thanks PASPA grant UNAM, Mexico and the Centre de Recherches Mathématiques, Université de Montréal for the kind hospitality while on sabbatical leave during which this work was done. The research of AME was partially supported by a fellowship awarded by the Laboratory of Mathematical Physics of the CRM. The research of İY was partly supported by Hacettepe University Scientific Research Coordination Unit. Project Number: FBI-017-1405. He also thanks the Centre de Recherches Mathématiques, Université de Montréal and especially Professor Pavel Winternitz for kind hospitality during his sabbatical leave. 5

APPENDIX Explicitly, the functions F 1,...,F 1 figuring in 16-1 are given by F 1 = B 1 cosθ+ B sinθ+ D 1 cos4θ+ D sin4θ F = 1 B cosθ B 1 sinθ D 1 sin4θ D cos4θ r + A cos θ+a 4 sin θ +C 1 sinθ+c cosθ, F = 1 B r cosθ B 1 sinθ+d 1 sin4θ D cos4θ + 1 A r 4 sinθ +A cosθ C cosθ+c 1 sinθ B sinθ B 4 cosθ +A sinθ+a 1 cosθ, r F 4 = D r 4 1 cos4θ+ D sin4θ B 1 cosθ B sinθ + 4 A r 4 cosθ A sinθ C 1 cosθ+c sinθ 4 B r cosθ+b 4 sinθ + 4 A cosθ A 1 sinθ, r F 5 = 6 D r 1 cos4θ+ D sin4θ + A 4 cosθ A sinθ r + C 1 cosθ C sinθ +B cosθ+b 4 sinθ, F 6 = A 4cosθ A sinθ+c 1 cosθ C sinθ 9 r D 1cos4θ+D sin4θ, F 7 = 1 r D 1sin4θ D cos4θ+ 1 r A 4sinθ +A cosθ 15C 1 sinθ+ C cosθ B sinθ B 4 cosθ, F 8 = r B 1cosθ+ B sinθ 5D 1 cos4θ+ D sin4θ + r A 4cosθ A sinθ 9C 1 cosθ C sinθ 5 r B cosθ+b 4 sinθ+ A cosθ A 1 sinθ, F 9 = 9 r 4 B 1sinθ B cosθ D 1 sin4θ D cos4θ + 15 r C 1sinθ +C cosθ A cosθ A 4 sinθ + 1 r B sinθ B 4 cosθ 9 r A 1cosθ+A sinθ, F 10 = r B cosθ B 1 sinθ+d cos4θ D 1 sin4θ + A cosθ+a 4 sinθ+c 1 sinθ+c cosθ, F 11 = r B 1cosθ+B sinθ 5D 1 cos4θ+ D sin4θ 4 B cosθ+b 4 sinθ + 1 r A 4cosθ A sinθ 9C 1 cosθ C sinθ, 6,

F 1 = r B 1sinθ B cosθ 14D 1 sin4θ D cos4θ + r 11C 1sinθ+ C cosθ A cosθ A 4 sinθ + 6 r B sinθ B 4 cosθ A 1cosθ+A sinθ, F 1 = 4 r 4 B 1cosθ+ B sinθ 11D 1 cos4θ+ D sin4θ r A 4cosθ A sinθ 17C 1 cosθ C sinθ + 1 r B cosθ+b 4 sinθ r A cosθ A 1 sinθ. REFERENCES 1 Abouamal I and Winternitz P 018 Fifth-order superintergrable quantum system separating in Cartesian coordinates: Doubly exotic potentials J. Math. Phys. 59 0104 Bermudez D, Fernández C D J and Negro J 016 Solutions to the Painlevé V equation through supersymmetric quantum mechanics J. Phys. A-Math. Theor. 49 50 Bertrand J L F 187 Théorème relatif au mouvement d un point attiré vers un centre fixe C. R. Acad. Sci. 77 849 85 4 Boyer C P, Kalnins E G and Miller W Jr 1986 Stäckel-equivalent integrable Hamiltonian systems SIAM J. Math. Anal. 17 778-797 5 Capel J J, Kress J M and Post S 015 Invariant classification and limits of maximally superintegrable systems in D SIGMA 11 17 6 CariñenaJF,HerranzFJandRañadaMF017Superintegrablesystemson-dimensional curved spaces: Eisenhart formalism and separability J. Math. Phys. 58 0701 7 Celeghini E, Kuru Ş, Negro J and del Olmo M A 01 A unified approach to quantum and classical TTW systems based on factorizations Ann. Phys. 7 7 8 Chanu C, Degiovanni L Rastelli G 011 First integrals of extended Hamiltonians in n+1 dimensions generated by powers of an operator SIGMA 7 1 9 Chanu C, Degiovanni L Rastelli G 01 Superintegrable extensions of superintegrable systems SIGMA 8 1 10 Drach J 195 Sur l integration logique des équations de la dynamique à deux variables: Forces constructives. Intégrales cubiques. Mouvements dans le plan Comptes Rendus Acad. Sci. 00 7

11 Drach J 195 Sur l integration logique et sur la tranformation des équations de la dynamique à deux variables: Forces conservatives. Intégrales cubiques. Mouvements dans le plan Comptes Rendus Acad. Sci. 00 599 1 Escobar-Ruiz A M, López Vieyra J C and Winternitz P 018 Fourth order superintegrable systems separating in Polar coordinates. I. Exotic Potentials J. Phys. A-Math. Theor. 50 49506 1 Escobar-Ruiz A M, Winternitz P and Yurduşen İ General Nth -order superintegrable systems separating in Polar coordinates in preparation 14 Fernández C D J and Morales-Salgado V S 016 SUSY partners of the truncated oscillator, Painlevé transcendents and Bäcklund transformations J. Phys. A-Math. Theor. 49 1950 15 Friš J, Mandrosov V, Smorodinsky Y A, Uhlíř M and Winternitz P 1965 On higher symmetries in quantum mechanics Phys. Lett. 16 54 6 16 Friš J, Smorodinsky Y A, Uhlíř M and Winternitz P 1966 Symmetry groups in classical and quantum mechanics Yadra Fiz. 4 65 5 17 Gonera C 01 On the superintegrability of the TTW model Phys. Lett. A 76 41 18 Gravel S 004 Hamiltonians separable in Cartesian coordinates and third-order integrals of motion J. Math. Phys. 45 100 19 19 Gravel S and Winternitz P 00 Superintegrability with third-order integrals in quantum and classical mechanics J. Math. Phys. 4 590 1 0 Güngör F, Kuru Ş, Negro J and Nieto L M 017 Heisenberg-type higher order symmetries of superintegrable systems separable in cartesian coordinates Nonlinearity 0 1788 808 1 Hakobyan T, Nersessian A and Shmavonyan H 017 Lobachevsky geometry in TTW and PW systems Phys. Atomic Nuclei 80 598 604 Hietarinta J 1998 Pure quantum integrability Phys. Lett. A 46 97 104 Hietarinta J, Grammaticos B, Dorizzi B and Ramani A 1984 Coupling-constant metamorphosis and duality between integrable Hamiltonian systems Phys. Rev. Lett. 5 1707 10 4 Kalnins E G, Kress J M and Miller W Jr 005 Second-order superintegrable systems in conformally flat spaces. I. Two-dimensional classical structure theory J. Math. Phys. 46 05509 5 Kalnins E G, Kress J M and Miller W Jr 005 Second-order superintegrable systems in conformally flat spaces. II. The classical two-dimensional Stäckel transform J. Math. Phys. 46 05510 8

6 Kalnins E G, Kress J M and Miller W Jr 005 Second order superintegrable systems in conformally flat spaces. III. Three-dimensional classical structure theory J. Math. Phys. 46 10507 7 Kalnins E G, Kress J M and Miller W Jr 006 Second order superintegrable systems in conformally flat spaces. IV. The classical D Stäckel transform and D classification theory J. Math. Phys. 47 04514 8 Kalnins E G, Kress J M and Miller W Jr 006 Second-order superintegrable systems in conformally flat spaces. V. Two- and three-dimensional quantum systems J. Math. Phys. 47 09501 9 Kalnins EG, Kress JMandMiller W Jr 010Superintegrability andhigher order integrals for quantum systems J. Phys. A-Math. Theor. 4 6505 0 Kalnins E G, Miller W Jr, Pogosyan G S 011 Superintegrability and higher-order constants for classical and quantum systems Phys. Atomic Nuclei 74 914 8 1 Kalnins E G, Miller W Jr and Winternitz P 1976 The Group O4, Separation of Variables and the Hydrogen Atom SIAM J. Appl. Math. 0 60 64 Kalnins E G, Miller W Jr, Post S 010 Coupling constant metamorphosis and N th -order symmetries in classical and quantum mechanics J. Phys. A-Math. Theor. 4 050 Kalnins E G, Williams G C, Miller W Jr and Pogosyan G S 00 On superintegrable symmetry-breaking potentials in N-dimensional Euclidean space J. Phys. A-Math. Gen. 5 4755 7 4 Makarov A A, Smorodinsky J A, Valiev K and Winternitz P 1967 A systematic search for nonrelativistic systems with dynamical symmetries Nuovo Cimento A 5 1061 84 5 Marchesiello A, Post S and Šnobl L 015 Third-order superintegrable systems with potentials satisfying only nonlinear equations J. Math. Phys. 56 10104 6 Marquette I 009 Superintegrability with third order integrals of motion, cubic algebras, and supersymmetric quantum mechanics. II. Painlevé transcendent potentials J. Math. Phys. 50 0950 7 Marquette I 010 Superintegrability and higher order polynomial algebras J. Phys. A- Math. Gen. 4 150 8 Marquette I, Sajedi M and Winternitz P 017 Fourth order Superintegrable systems separating in Cartesian coordinates I. Exotic quantum potentials J. Phys. A-Math. Theor. 50 1501 9

9 Miller W Jr, Post S and Winternitz P 01 Classical and quantum superintegrability with applications J. Phys. A-Math. Gen. 46 4001 40 Popper I, Post S and Winternitz P 01 Third-order superintegrable systems separable in parabolic coordinates J. Math. Phys. 5 06105 41 Post S, Vinet L and Zhedanov A 011 An infinite family of superintegrable Hamiltonians with reflection in the plane J. Phys. A-Math. Theor. 44 50501 4 Post S and Winternitz P 010 An infinite family of superintegrable deformations of the Coulomb potential J. Phys. A-Math. Gen. 4 01 4 Post S and Winternitz P 015 General N th -order integrals of motion in the Euclidean plane J. Phys. A-Math. Theor. 48 40501 44 Quesne C 010 Superintegrability of the Tremblay-Turbiner-Winternitz quantum Hamiltonians on a plane for odd k J. Phys. A-Math. Theor. 4 08001 45 Rañada M F 01 Higher order superintegrability of separable potentials with a new approach to the Post-Winternitz system J. Phys. A-Math. Theor. 46 1506 46 Tremblay F, Turbiner A V and Winternitz P 009 An infinite family of solvable and integrable quantum systems on a plane J. Phys. A-Math. Theor. 4 4001 47 Tremblay F, Turbiner A V and Winternitz P 010 Periodic orbits for an infinite family of classical superintegrable systems J. Phys. A-Math. Theor. 4 0150 48 Tremblay F and Winternitz P 010 Third-order superintegrable systems separating in polar coordinates J. Phys. A-Math. Theor. 4 17506 49 Tsiganov A V 000 The Drach superintegrable systems J. Phys. A-Math. Gen. 7407 0