Extensions of Hamiltonian systems dependent on a rational parameter Claudia M. Chanu, Luca Degiovanni, Giovanni Rastelli arxiv:30.5690v [math-ph] 2 Oct 203 Dipartimento di Matematica G. Peano, Università di Torino. Torino, via Carlo Alberto 0, Italia. e-mail: claudiamaria.chanu@unito.it luca.degiovanni@gmail.com giorast.giorast@alice.it July 7, 208 Abstract The technique of extension allows to build(n+)-dimensional Hamiltonian systems with a non-trivial polynomial in the momenta first integral of any given degree starting from a n-dimensional Hamiltonian satisfying some additional properties. Until now, the application of the method was restricted to integer values of a certain fundamental parameter determining the degree of the additional first integral. In this article we show how this technique can be generalized to any rational value of the same parameter. Several examples are given, among them the anisotropic oscillator and a special case of the Tremblay-Turbiner-Winternitz system. Introduction The extensions of natural Hamiltonians were introduced in [4] as a tool for building Hamiltonian systems with polynomial first integrals of any given degree. The technique of extension was developed as the generalization of an iterative procedure to generate a first integral of degree λ N {0} in the momenta [3] for the system 2 p2 r + ( r 2 2 p2 ψ + ) a sin 2, () (λψ) This superintegrable system was already considered in [2], where it was remarked the existence of additional polynomial first integrals for several values of λ Q and a general form for odd positive integers was conjectured. Then, a generalization of the system was introduced by Tremblay, Turbiner and Winternitz (TTW) [] and, subsequently, studied by many authors (see for instance [8, 7, 9] and references therein), proving the superintegrability of the TTW system, and of several related new systems, for any rational value of λ.
The technique of extensions was generalized in [5] and employed to build new superintegrable systems from existing ones in [6]. The extension theory can be easily generalized to complex manifolds, see [7] for details. In this paper we develop a more general iterative procedure that in particular allows us to construct, for any positive rational λ = m/n, a polynomial first integral of degree m+n of the system (), as well as of the two uncoupled harmonic oscillators with rational ratio λ of the frequencies 2 p2 y + 2 p2 x +ω2 x 2 + ω2 λ 2y2. (2) For both systems, the case of integer λ is recalled in Sec. 2, where the extension technique is briefly summarized. Its generalization is exposed in Sec. 3, while Sec. 4 contains several detailed examples. In the Appendix, the trigonometric tagged functions, employed extensively in the article, are defined and some of their properties are shown. 2 Extension of Hamiltonian systems The extension procedure has been characterized in[4, 5]. Given a Hamiltonian L on a 2n-dimensional Poisson manifold Q, we construct the (2n + 2)-dimensional Poisson manifold M = T Q, where T is the cotangent bundle of a onedimensional manifold with the canonical symplectic form dp u du. The main result (Proposition in [5]) states that a Hamiltonian L on the Poisson manifold Q, with Hamiltonian vector field X L, admits an extension of the form 2 p2 u +α(u)l+β(u), with a polynomial first integral of the form U m (G), where U = p u +γ(u)x L, if and only if there exists a function G on Q that satisfies the relation X 2 L (G) = 2m(cL+L 0)G, (3) for some m N {0} and c,l 0 R. If a solution is found for (mc,l 0 ) with c 0, then, without loss of generality we can set L 0 =0. When the relation (3) holds, the function γ is any solution of γ +2cγγ = 0 (4) and the functions α and β are calculated directly from γ through the relations α = mγ, β = ml 0 γ 2. (5) The form of γ, and therefore of the operator U, depends on the value of c. If c 0, then for any κ R a solution of (4) and the corresponding operator are γ = T κ (cu), U = p u + T κ (cu) X L. 2
The choice L 0 = 0 leads to β = 0, hence, up to inessential constants, the extended Hamiltonian is 2 p2 u + mc Sκ(cu) 2 L where the trigonometric tagged functions S κ, T κ are defined in the Appendix. Conversely, if c = 0, then for any A R γ = Au, U = p u AuX L and the extended Hamiltonian, written up to inessential constants, becomes 2 p2 u +mal+ml 0A 2 u 2. The parameters in equation (3) and (4) are chosen in order to obtain a function γ, and therefore an operator U, not depending on m. This choice, going back to [4], gives a simpler form for γ but has a drawback: G seems to depend on the three independent parameters m, c and L 0 instead that on the two parameters mc and ml 0, as in fact it is. By setting instead c = mc and L 0 = ml 0, equations (3) and (4) become respectively X 2 L(G) = 2( cl+ L 0 )G, (6) mγ +2 cγγ = 0. (7) Within this notation, G explicitly depends on two parameters only and it is clear that G remains the same for different values of m, but γ and U become dependent on m. A simple calculation shows that if γ(u;c) satisfies (4) then mγ(u; c) satisfies (7); moreover, relations (5) hold unchanged. In the following, the choice of parameters made in equations (6) and (7) will be used and therefore, for c 0, κ R, γ = m T κ ( cu), and the extended Hamiltonian becomes Conversely, for c = 0, A R, we have and the extended Hamiltonian becomes U = p u + m T κ ( cu) X L, (8) 2 p2 u + m2 c (9) ( cu)l. S 2 κ γ = mau, U = p u maux L (0) 2 p2 u +m2 AL+m 2 L0 A 2 u 2. () When L is a natural Hamiltonian, and therefore Q is the cotangent bundle of a Riemannian manifold, the configuration manifolds of the extended Hamiltonians are warped manifolds. Several explicit examples of extensions are exposed into details in [3, 4, 6]. 3
Example. (See also[3]) We apply the extension procedure to the Hamiltonian L = 2 p2 φ + a sin 2 (φ). (2) One of the possible extensions (with κ = 0) is H cal = ( 2 p2 u + λ2 u 2 2 p2 φ + a ) sin 2, (3) φ which coincides with the Calogero-type system () after the rescaling φ = λψ. In this example c 0 and λ 2 = m 2 / c. The relation (6) imposes a very strong constraint on the parameters: if the function G does not depend on the momenta, then (6) has a solution only for c = and L 0 = 0. Hence, λ 2 = m 2 and λ is necessarily an integer. Example 2. (See also [6]) We consider the system of two uncoupled harmonic oscillators (2). By rescaling u = λy and dividing the Hamiltonian by the constant factor λ 2, we get the equivalent Hamiltonian H osc = 2 p2 u +λ2 ( 2 p2 x +ω2 x 2 )+λ 2 ω 2 u 2. (4) Here we have an extension of the one-dimensional harmonic oscillator with c = 0, λ 2 = m 2 A and ω 2 = A L 0. Even in this case, the equation (6) admits a solution not depending on the momenta only for L 0 = ω 2. This constraint implies A = and λ 2 = m 2, forcing λ to be an integer. Both systems of the above examples admit the maximal number of functionally independent first integrals, namely, H, L and U λ G, for λ integer, so that they are superintegrable. However, it is well known that Hamiltonians (3) and (4) are superintegrable for any rational value of the parameter λ ([7, 9, 8]). In the following we show an effective procedure to construct a third independent first integral for non integers values of λ by means of a more general extension procedure. 3 Extensions with rational parameters It is useful to slightly modify the expression of the extensions introduced so far and rename them as follows Definition. The Hamiltonian H m = 2 p2 u +m 2 α(u)l+m 2 β(u), (5) defined on the Poisson manifold T Q, is a m-extension of L generated by G, if K m = (p u +m γ(u)x L ) m (G) (6) is a first integral of H m. 4
The comparison between (8)-() and (5)-(6) shows that α, β, γ are independent of m and that in particular the following relations hold α = α m 2, β = β m 2, γ = γ m, (see Table for their expressionsaccording to the value of c). Moreover, G must satisfy (6). c = 0 c 0 α = γ = A c Sκ 2( cu) β = L 0 γ 2 = L 0 A 2 u 2 0 γ = Au T κ ( cu) Table : Functions involved in the m-extension of L In order to write as extensions the Hamiltonians (3), (4) with λ = m n, one, naively, should divide α, β, γ by n 2. This corresponds, roughly speaking, to a m-extension of L/n 2. More precisely, Definition 2. The Hamiltonian H m,n = 2 p2 u + m2 m2 α(u)l+ n2 n β(u), 2 (7) defined on the Poisson manifold T Q, is called a (m,n)-extension of L generated by a function G n if ( K m,n = p u + m m(gn L) n 2 γ(u)x ) (8) is a first integral of H m,n, where α, β, γ are defined as in Table. Lemma. A Hamiltonian L admits a (m, n)-extension generated by a function G n, if and only if the function G n satisfies X 2 L(G n ) = 2n 2 ( cl+ L 0 )G n. (9) Proof. The Hamiltonian L (n) = n2l admits a m-extension if and only if there exists a function G n satisfying X 2 L (n) (G n ) = 2 The m-extension of L (n) is then given by ( cl (n) + ) (n) L 0 G n. with first integral 2 p2 u +m2 αl (n) +m 2 β(n), (20) (p u +m γx L (n)) m (G n ), 5
where β (n) = implies L (n) 0 γ 2 and α, β, γ are given by Table. The definition of L (n) Hence, after setting β = L 0 γ 2, one has ( ) c n 4X2 L (G (n) n) = 2 n2l+ L 0 G n. β (n) = L (n) 0 γ2 = n 2 L 0 γ 2 = n 2 β and the extension (20) becomes 2 p2 u + m2 m2 αl+ n2 n β 2 which coincides with the (m, n)-extension of L. Remark. The Proof of Lemma shows that it is not restrictive to assume in Definition 2 that α, β, γ are those given in Table. Consequently, the search for extensions with rational parameter reduces to the search for solutions of (9). Theorem 2 shows how to construct iteratively solutions G n of (9) starting from a known solution G of (6). Theorem 2. Let G be a function on Q satisfying with X L (Λ) = 0, then the recursion satisfies, for any n N {0}, X 2 L(G) = ΛG G = G, G n+ = X L (G)G n + n GX L(G n ), (2) X 2 L (G n) = n 2 ΛG n. Proof. By induction on n. For n =, one has G 2 = 2GX L (G) and the relation X 2 L (G 2) = 4ΛG 2 is straightforward. Then, let us assume that G n satisfies X 2 L (G n) = n 2 ΛG n. Recalling that for two functions A and B the formula holds, we get X 2 L (AB) = X2 L (A)B +2X L(A)X L (B)+AX 2 L (B) X 2 L (G n+) = X 2 L [ X L (G)G n + ] n GX L(G n ) = XL 3 (G)G n +2XL 2 (G)X L(G n )+X L (G)XL 2 (G n)+ [ ] X 2 n L(G)X L (G n )+2X L (G)XL(G 2 n )+GXL(G 3 n ) = ΛX L (G)G n +2ΛGX L (G n )+n 2 ΛX L (G)G n + [ ] ΛGX L (G n )+2n 2 ΛX L (G)G n +n 2 ΛGX L (G n ) n = (+2n+n 2 )ΛX L (G)G n + +2n+n2 ΛGX L (G n ) n = (n+) 2 ΛG n+. 6
Corollary 3. A Hamiltonian L admits a m-extension H m if and only if it admits a (m,n)-extension H m,n for any positive rational m/n. Proof. If L admits m-extensions, then there exists G satisfying(6). By Theorem 2 with Λ = 2( cl+ L 0 ), (22) we can construct G n verifying the condition (9) for any n. Hence, we can construct the (m, n)-extension (7). The converse is straightforward. Proposition 4. The closed form for G n satisfying the recursion (2) is G n = [ n 2 ] where [ ] denotes the integer part. Proof. Since ( ) n Λ k G 2k+ (X L G) n 2k, (23) 2k + X L (G 2k+ (X L G) n 2k ) = (2k +)G 2k (X L G) n 2k and using the identity n+2k+ n ( n 2k+ +(n 2k )ΛG 2k+2 (X L G) n 2k 2, ) + n 2k+ n ( n ) ( 2k = n+ 2k+), by applying the recursion (2) we have that for n = 2i G 2i+ = X L (G)G 2i + 2i GX L(G 2i ) = = i i + ( )( 2i + 2k+ ) 2k + Λ k G 2k+ (X L G) 2i 2k 2i ( ) 2i 2k 2i Λ k+ G 2k+3 (X L G) 2i 2k 2 2i 2k + i ( ) 2i 2i+2h+ Λ h G 2h+ (X L G) 2i 2h 2h+ 2i h=0 i ( ) 2i 2h+ 2i + Λ h G 2h+ (X L G) 2i 2h n 2h h= = (2i+)G(X L G) 2i +Λ i G 2i+ i [ ) + + 2i 2h+ = h= i ( 2i+2h+ 2i 2i 2h+ ( ) 2i+ Λ k G 2k+ (X L G) 2i 2k, 2k+ 2i ( 2i ) ] 2h Λ h G 2h+ (X L G) 2i 2h 7
that is (23) with n = 2i+. The case n = 2i+ is analogous: G 2i+2 = X L (G)G 2i+ + 2i+ GX L(G 2i+ ) i ( )( ) 2i+ 2k + = + Λ k G 2k+ (X L G) 2i+ 2k 2k+ 2i+ = + i i h=0 + 2i 2k 2i+ ( ) 2i+ Λ k+ G 2k+3 (X L G) 2i 2k 2k + ( ) 2i+ 2i+2h+2 Λ h G 2h+ (X L G) 2i+ 2h 2h+ 2i+ i ( ) 2i 2h+ 2i Λ h G 2h+ (X L G) 2i 2h n 2h h= = (2i+2)G(X L G) 2i+ i [ ) + + 2i 2h+2 = h= i 2i+2h+2 2i+ that is (23) with n = 2i+2. ( 2i+ 2h+ 2i+ ( ) 2i+2 Λ k G 2k+ (X L G) 2i+ 2k, 2k+ ( 2i+ ) ] 2h Λ h G 2h+ (X L G) 2i+ 2h It follows that, if L is quadratic in the momenta, then G n is polynomial in the momenta of degree at most n(+deg(g)). From [5], we can derive the explicit expression of the first integrals K m,n of any (m, n)-extension Proposition 5. Let G n be a recursion with G = G satisfying (6), that is with Λ given by (22). We have with D m,n = n K m,n = P m,n G n +D m,n X L (G n ), (24) P m,n = [(m )/2] [m/2] ( ) m (m ) 2k 2k n γ p m 2k u Λ k, ( ) m (m ) 2k+ 2k + n γ p m 2k u Λ k, m >, where [ ] denotes the integer part and D,n = m n 2 γ. We remark that in Theorem 3 of [5] the upper limit of the sum in D m is misprinted as [m/2] instead of [(m )/2]. 8
Proof. Since the first integralk m,n of a (m,n) extension oflis the first integral of a m extension of L/n 2, by Lemma it follows that Theorem 3 of [5] can be easily adapted, once we set γ γm/n 2, 2m(cL+L 0 ) n 2 Λ and G G n. The result is straightforward. It follows that the degree of the first integral K m,n, for L quadratic in the momenta, is m+deg(g n ) m+n(+deg(g)). When m and n are not reciprocally prime, one can expect that the polynomial K m,n factorizes into a number of factors, some of them again in the form K r,s, where m = a 0 r, n = a 0 s, a 0,r,s N {0}. Indeed, the computation of several examples suggests that K r,s is then a divisor of K m,n. 4 Examples of (m, n)-extensions Once a m-extension of a Hamiltonian L is known, it is straightforward to build any (m, n)-extension of L, as the following examples show. 4. Two uncoupled oscillators Let us consider the two uncoupled harmonic oscillators described by the Hamiltonian (4), extension of L = 2 p2 x +ω 2 x 2, for integer values of the parameter λ. In this case the vector field X L is given by X L = p x x 2ω2 x. p x WecaneasilyfindaG(x)notdependingonthemomentasatisfyingthecondition (6), with c = 0, L 0 = ω 2 ; indeed (6) reduces to p 2 x (G ) 2ω 2 (xg G) = 0. Thus, G = x up to an inessential multiplicative constant and X L (G) = p x. Then, by setting A =, we can construct the (m,n)-extension of L H m,n = ( ( m ) ) 2p p 2 2 u + 2 n x +ω 2( m ) 2(x 2 +u 2), n for any m,n positive integers, which is equivalent (up a rescaling and a constant factor) to (2) with λ = m/n. The first terms in the recursion (2) are G = x, G 2 = 2x p x, G 3 = 3x p 2 x 2ω 2 x 3, G 4 = 4x p 3 x 8ω2 x 3 p x, G 5 = 5x p 4 x 20ω 2 x 3 p 2 x +4ω 4 x 5. 9
Alternatively, formula (23) becomes in this case G n = [ n 2 ] ( ) n ( 2ω 2 ) k x 2k+ px n 2k. 2k+ For some values of (m,n) we get the following first integrals K, = xp u up x, K,2 = 2x p x p u u ( ) 2 p2 x ω 2 x 2, K 2,2 = 2(xp u up x )(p x p u +2ω 2 xu), K 3,2 = 2x p x p 3 u 9 2 u p2 x p2 u +9ω2 x 2 u p 2 u 27ω2 xu 2 p x p u, + 27 4 ω2 u 3 p 2 x 27 2 ω4 x 2 u 3. We remark that K 2,2 is a multiple of K,. This factorization is related to the fact that (,) and (2,2) represent the same rational number. 4.2 Calogero-type systems We consider now the (m, n)-extensions of the Hamiltonian (2). In this case the vector field X L is given by X L = p φ φ + 2acosφ sin 3. φ p φ WecaneasilyfindaG(φ)notdependingonthemomentasatisfyingthecondition (6), with c =, L 0 = 0; indeed (6) splits into (G +G) = 0, G + sinφ cosφ G = 0. Thus, up to an inessential multiplicative constant, a solution is G = cosφ and X L (G) = p φ sinφ. Then, by (2) we can construct the sequence of G n, whose first terms are G = cosφ, G 2 = sin2φ p φ, G 3 = cos3φ p 2 φ 2acos3 φ sin 2 φ, G 4 = sin4φ p 3 φ 8acos3 φ sinφ p φ, G 5 = cos5φ p 4 φ + 4a(6cos2 φ 5)cos 3 φ sin 2 φ We can construct the (m,n)-extension of L H m,n = 2 p2 u + m2 n 2 S 2 κ (u) 0 ( 2 p2 φ + p 2 φ + 4a2 cos 5 φ sin 4 φ. ) a sin 2, (φ)
for any m,n positive integers. According to the sign of the parameter κ, the extended Hamiltonian is defined on the cotangent bundle of a sphere (κ > 0), of a plane (κ = 0) orof a pseudo-sphere(κ < 0). Moreover,for κ = 0 the extension coincides with (3) that, after the rescaling φ = λψ, becomes the Hamiltonian () H = 2 p2 u + ( ) u 2 2 p2 ψ + a sin 2, (λψ) for any rational λ = m/n. This Hamiltonian generalizes the Calogero three particle chain without harmonic term (obtained for λ = 3, see [2], also known as Jacobi system) and it was the starting point of our work in [3]. This is a particular case of the TTW system []. As an example, we compute a first integral of the extension (3) for (m,n) = (2,3) i.e., λ = 2/3, K 2,3 = cos3φ p 2 up 2 φ + 4 3u sinφ(4cos2 φ ) p u p 3 φ + 4cos3φ 9u 2 p 4 φ 2acos3 φ sin 2 φ p2 u + 8acos2 φ usinφ p up φ + 8a(5cos2 φ 3)cosφ 9u 2 sin 2 p 2 φ + 6a2 cos 3 φ φ 9u 2 sin 4 φ. 4.3 Three-sphere Any example of m-extension taken from [3, 4, 5, 6] can be easily transformed into a (m, n)-extension. For example, let us consider from [4] the m-extension of the geodesic Hamiltonian ( ) L = 2 p 2 η + p2 ξ sin 2 η + p2 ξ 2 cos 2 η on the three-sphere S 3, with coordinates (η,ξ,ξ 2 ) where 0 < η < π/2, 0 ξ i < 2π and the parametrization in R 4 given by x = cosξ sinη, y = sinξ sinη, z = cosξ 2 cosη, t = sinξ 2 cosη. Equation (6) admits a complete solution (i.e., a solution depending on the maximal number of parameters a i ) if and only if mc = c equals the curvature K = of the three-sphere [4].The complete solution is G = (a 3 sin(ξ )+a 4 cos(ξ ))sin(η)+(a sin(ξ 2 )+a 2 cos(ξ 2 ))cos(η) = a 4 x+a 3 y +a 2 z +a t, with a i constants. We look for compatible potentials V, i.e., functions that can be added to L in such a way that G satisfies (6) for the natural Hamiltonian L+V. It is easier to find compatible V when some of the parameters a i in G are assigned. For example, if a 2 = a 3 = a 4 = 0, a compatible potential is V = sinξ cosξ 2 cosηsinη.
With this choice of (a i ), the first three G n obtained from G = sinξ 2 cosη are G = sinξ 2 cosη, G 2 = sin2ηsin 2 ξ 2 p η +2sinξ 2 cosξ 2 p ξ2, G 3 = cos3ηsin 3 ξ 2 p 2 η 6sinηsin2 ξ 2 cosξ 2 p η p ξ2 cos3 η sin 2 η sin3 ξ 2 p 2 ξ + sinξ 2(3cos 2 ξ 2 cos 2 ηsin 2 ξ 2 ) cosη the (m,n)-extension of L+V is H m,n = ( m ) 2 L+V 2 p2 u + n Sκ 2(u), and, for example, a first integral of H,2 is K,2 = T κ (u) cos2 ηsin 2 ξ 2 2sin 2 η p 2 ξ 2 2 sin3 ξ 2 sinξ cos 2 η, cosξ 2 sinη ( 2 cos2ηsin2 ξ 2 p 2 η sinηcosξ 2sinξ 2 p η p ξ2 cosη p 2 ξ + cos2 ξ 2 cos 2 ηsin 2 ξ 2 2cos 2 η 2sin2ηsin 2 ξ 2 p u p η +2sinξ 2 cosξ 2 p u p ξ2. 5 Conclusions ) p 2 ξ 2 sin2 ξ 2 sinξ cosη sinηcosξ 2 In this article we show that the technique of extensions can be modified in order to be applied successfully also to a class of Hamiltonian systems (including some known superintegrable ones) depending on rational values of a parameter. Many properties of the first integrals obtained, for instance their factorization, need a deeper analysis and the possibility of the use of the technique to build new extended systems with rational parameters is still unexplored. This short exposition is certainly not a complete theory of extensions with rational parameters, but it represents a solid ground for building such a theory. 6 Appendix The trigonometric tagged functions sin κx κ κ > 0 S κ (x) = x κ = 0 sinh κ x κ < 0 κ cos κx κ > 0 C κ (x) = κ = 0 cosh κ x κ < 0, T κ (x) = S κ(x) C κ (x), are employed, explicitly or not, by several scholars (see [0], []) and appear in several branches of mathematics. The trigonometric tagged functions satisfy a number of properties analogous to those of ordinary trigonometric functions. Their main advantage is to unify 2
trigonometric and hyperbolic functions in a homogeneous way. From the definition it follows that S κ (x) is an odd function while C κ (x) is an even function but, if κ 0 and x R they are no more periodic. The basic properties are C 2 κ(x)+κs 2 κ(x) =, S κ (x±y) = S κ (x)c κ (y)±c κ (x)s κ (y), C κ (x±y) = C κ (x)c κ (y) κs κ (x)s κ (y), from these the duplication formulas can be obtained S κ (2x) = 2S κ (x)c κ (x), C κ (2x) = C 2 κ (x) κs2 κ (x) = { 2C 2 κ (x), 2κS 2 κ (x). The bisection formula Cκ 2 (x) = +C κ(2x) 2 is always true, while one has to set κ 0 in order to obtain Sκ 2 (x) = C κ(2x). 2κ The functions S κ (x) and C κ (x) are related through differentiation: d dx S κ(x) = C κ (x) d dx C κ(x) = κs κ (x). Hence, the linear combinations of S κ (x) and C κ (x) provide the general solution of the differential equation F +κf = 0. Acknowledgement This work has been partially (CC) supported by PRIN 200/20 Research project Teorie geometriche e analitiche dei sistemi Hamiltoniani in dimensioni finite e infinite grant no 200JJ4KPA 007. References [] D. Bao, S. Chern and Z. Shen, An introduction to Riemann-Finsler Geometry (Springer-Verlag, New York, 2000). [2] C. Chanu, L. Degiovanni and G. Rastelli, Superintegrable three-body systems on the line, J. Math. Phys. 49, 290 (2008). [3] C. Chanu, L. Degiovanni and G. Rastelli, Polynomial constants of motion for Calogero-type systems in three dimensions, J. Math. Phys. 52, 032903 (20). [4] C. Chanu, L. Degiovanni and G. Rastelli, First integrals of extended Hamiltonians in (n + )-dimensions generated by powers of an operator, SIGMA 7, 038 (20). 3
[5] C. Chanu, L. Degiovanni and G. Rastelli, Generalizations of a method for constructing first integrals of a class of natural Hamiltonians and some remarks about quantization J. Phys.: Conf. Ser. 343, 020 (202). [6] C. Chanu, L. Degiovanni and G. Rastelli, Superintegrable extensions of superintegrable systems, SIGMA 8, 070 (202). [7] A.J. Maciejewski, M. Przybylska and H. Yoshida, Necessary conditions for super-integrability of a certain family of potentials in constant curvature spaces, J. Phys. A: Math. Theor. 43, 38200 (200). [8] E.G. Kalnins, J.M. Kress and W. Miller, Jr., Tools for Verifying Classical and Quantum Superintegrability, SIGMA 6, 066 (200). [9] W. Miller, Jr., S. Post and P. Winternitz, Classical and quantum superintegrability with applications, arxiv:309.2694, (203). [0] M.F. Rañada, M. Santander, Superintegrable systems on the twodimensional sphere S2 and the hyperbolic plane H2, J. Math. Phys. 40, 5026 5057 (999). [] F. Tremblay, A.V. Turbiner and P. Winternitz An infinite family of solvable and integrable quantum systems on a plane, J. Phys. A: Math. Theor. 42, 24200 (2009). 4