arxiv:math-ph/ v2 23 Oct 2003

Transcription

1 arxiv:math-ph/03008v 3 Oct 003 Hamiltonians separable in cartesian coordinates and third-order integrals of motion Simon Gravel Département de phsique et Centre de recherche mathématiques Université de Montréal, C.P.618, Succursale Centre-Ville Montréal, Québec H3C 3J7 Canada Simon.Gravel@UMontreal.CA June 14, 018 Abstract We present in this article all Hamiltonian sstems in E() that are separable in cartesian coordinates and that admit a third-order integral, both in quantum and in classical mechanics. Man of these superintegrable sstems are new, and it is seen that there eists a relation between quantum superintegrable potentials, invariant solutions of the Korteweg-De Vries equation and the Painlevé transcendents. 1

2 I Introduction In classical mechanics, an n-dimensional Hamiltonian sstem is called Liouville integrable if it allows n functionall independent integrals of motion in involution (including the Hamiltonian), that is {H,X i } = 0, {X i,x j } = 0, i,j. The Hamiltonian H = H( 1,... n,p 1,...,p n ) and the integrals of motion X i = X i ( 1,... n,p 1,...,p n ) must be well defined functions on phase space ([1, 10]). The sstem is superintegrable if it allows more than n functionall independent integrals, n of them in involution. It is called maimall superintegrable if it allows n 1 integrals of motion. The best known superintegrable sstems in n dimensions are the harmonic oscillator V = ωr and the Coulomb potential V = α r, and the are indeed maimall superintegrable. This ma be closel related to Bertrand s theorem ([1, ]) which states that these are the onl rotationall invariant sstems for which all finite trajectories are closed. In quantum mechanics, a Hamiltonian sstem is said to be integrable if there eists a set {X i } of n well defined, algebraicall independent operators (including the Hamiltonian) that commute pairwise. It is superintegrable if it possesses further independent operators, {Y j } that commute with the Hamiltonian. The Y j do not necessaril commute with each other, nor with the X i. The independence of operators in quantum mechanics remains to be defined rigorousl [1, 13, 16, 8]. Since we are dealing here onl with polnomial differential operators, we can proceed b analog with the classical case, keeping in mind that a rigorousdefinition will be needed as soon as we will wantto make some more general statements. This choice of a definition will be used onl for discussion purposes, since we will find all potentials that admit third-order integrals and all their integrals. The results obtained will therefore hold for an definition of the independence of operators. Integrable and superintegrable sstems, both in quantum and in classical mechanics, attracted considerable interest in the last ears. Etensive literature eists about sstems with second-order integrals of motion, either in euclidian space [6, 7, 8, 9, 1, 9], or in spaces with nonzero constant [18, 4] or nonconstant curvature [19]. As long as there was no magnetic field in the Hamiltonian, the quantum and classical integrals of motion obeed the same determining equations, and therefore quantum and classical integrabilit were ver similar. Both properties were related to separation of the Hamilton-Jacobi or Schrödinger equations, and also to eact solvabilit [6] and generalized smmetries [5]. Sstems with higher-order integrals have been studied and classified as earl as 1935 in a well-known paper b Drach. This paper considered classical Hamiltonians in comple Euclidian space. Efforts were made recentl to understand and classif more completel sstems with higher-order integrals in classical

3 [1, 7] and quantum mechanics [1]. In spite of these efforts still relativel few such sstems are known. This paper is the logical sequel of a sstematic search for superintegrable sstems with higher order integrals started in [1]. Here we consider two-dimensional real Euclidian space with a one-particle Hamiltonian; H = 1 ( p +p ) +V(,). We request the eistence of two additional integrals of motion, one of second order in the momenta and the other of third order. The condition of eistence for second-order integral implies, both in classical and quantum mechanics, that the Hamiltonian be separable in cartesian, polar, parabolic or elliptic coordinates. In this paper we consider potentials that are separable in cartesian coordinates; H = 1 ( p +p ) +V1 ()+V (). We found all such sstems which admit third-order integrals. Quantum and classical mechanics will be treated simultaneousl, for the conditions of eistence of integrals of motion, even though not equivalent, are quite similar in both cases. II Eistence of a third-order integral In quantum and classical mechanics, the general third-order commuting operator X = 3 i+j=0 can be reduced to a much simpler form, P ij (,)p i pj X = i,j,k i+j+k=3 L 3 =p p. A ijk {L i 3,p j p k }+{g 1 (,),p }+{g (,),p } (1) where the A ijk are real constants, and the g i real functions. The bracket is the anticommutator. It is not needed in classical mechanics, but in quantum mechanics it allowsus to get rid of terms with even powersof the p i and to make sure the operator is self-adjoint. Furthermore, its use allows us to see clearl the relations between the quantum and the classical case. Indeed, it was found in [1] that the requirement that the operator commutes (or Poisson-commutes) with the Hamiltonian H = 1 ( p +p ) +V(,). 3

4 implies equations that behave well in the classical limit. Namel, commutativit implies ( 0 = g 1 V +g V h f 1 V +f V +f 3 V +f 4 V 4 ) +8A 300 (V V )+(A 10 V +A 01 V ), (g 1 ) = 3f 1 ()V +f (,)V, (3) (g ) = f 3 (,)V +3f 4 ()V, (4) (g 1 ) +(g ) = (f (,)V +f 3 (,)V ), (5) in quantum mechanics, where () f 1 () = A A 10 A 10 +A 030, f (,) = 3A 300 A 10 +A 01 +A 10 A 111 +A 01, f 3 (,) = 3A 300 +A 10 A 01 +A 111 A 10 +A 01, f 4 () = A A 01 +A 10 +A 003. The equations in classical mechanics are obtained b setting h = 0. We ma notice from the quantum equations, or directl from the condition [H,X] = 0, that we can epress all quantum integrable potentials as h Ṽ(,) where Ṽ does not depend on h. It is often more natural and interesting though to choose arbitrarparameters contained in Ṽ to be depending on h, so that the potential V does not vanish in the classical limit. One ma alwas verif though that through an appropriate transformation of the arbitrar parameters one can write the potential as being proportional to h The three last equations, identical in the quantum and classical cases, ield a linear compatibilit condition for V, which reads 0 = f 3 V +(f 3f 4 )V +( 3f 1 +f 3 )V f V +(f f 3 )V +( 3f 1 +f +f 3 3f 4 )V +( f +f 3 )V +( 3f 1 +f f 3 )V +( f +f 3 3f 4 )V. (6) Further nonlinear compatibilit conditions can be obtained from () to (5) for the potential, and these are listed in [1]. The are quite complicated though and were not used for the results stated in this paper. 4

5 III Potentials separable in cartesian coordinates If we set V = V 1 ()+V () in equations () to (5), we find ( 0 = g 1 V 1 +g V h f 1 V 1 +f 4 V 4 +8A 300 (V V 1 )+(A 10 V 1 +A 01 V ) ), (7) (g 1 ) = 3f 1 ()V 1 +f (,)V, (8) (g ) = f 3 (,)V 1 +3f 4 ()V, (9) (g 1 ) +(g ) = (f (,)V 1 +f 3 (,)V ), (10) with h = 0 in the classical case. Equations (8) and (9) are readil integrated, so in the cartesian case two equations remain to be solved. The compatibilit condition (6) allows us to find ODEs for V 1 and V. If we set alternativel = 0 and = 0, we find (A 10 +A 111 +A 01 )V (3) 1 ()+4(A 10 +A 111 )V 1 ()+1A 10V 1 () = a+b (11) (A 01 A 111 +A 01 )V (3) ()+4(A 01 A 111 )V ()+1A 01 V () = c +d (1) The solutions to the homogeneous part of these equations are easil found and brought to a simple form b translations in and. If we take the first one for definiteness, we have four different tpes of solution. When A 10 0, we have two possible solutions, V 1hom = c 1 (+α) + c ( α) V 1hom = c 1 + c 3. If A 10 = 0 and A we get V 1hom = c 1 +c. Finall if onl A 01 0, the solution ma be brought to the form V 1hom = c +c 1. Special solutions are also simple. If A 10 0, we have V 1part = α +β. Otherwise, when A 111 0, V 1part = α 3 + β, finall, if onl A 01 0, V 1part = α 4 +β 3. Provided that (11) or (1) do not vanish triviall, we can choose V 1 or V, respectivel, amongst the following functions: (A.1) f 1 = c1 (+α) + c ( α) +c 3 +c 4 (A.) f = c1 + c 3 +c 3 +c 4 5

6 Superintegrable potentials Leading-order terms of the integrals V a =a( + ) L 3 ;{L,p p }; {L,p }; {L,p } V b =a( + )+ b + c {L,p p } V c =a( + )+ h + h L 3 ; {L,p p } V d =a( + )+ h L 3 ; {L,p p }; {L,p } V e = h 8α ( + )+ h 4 ( α) + h (+α) L 3 3α {L,p }; {L,p } V f = h 8α ( + )+ h 4 + h (+α) + h ( α) L 3 3α {L,p } V g = h 8α ( + ) 4 L 3 3α ({L,p }+{L,p }) + h ( α) + h ( α) + h (+α) + h (+α) V h =a(4 + )+ b +c p p V i =a(9 + ) {L,p } V j =a(9 + )+ h {L,p } V k = h 8α (9 + )+ h 4 (+α) + h ( α) {L,p } V l = h + a {L,p }; {L,p p }; p 3 V m = h + h L 3 ; {L,p }; {L,p } {L,p p }; p 3 ; p3 V n =a+ h {L,p }; p 3 ; p p V o = h +V() p 3 Table 1: Superintegrable potentials that satisf linear compatibilit conditions for nonzero parameters. (A.3) f 3 = c1 +c 3 +c 3 +c 4 (A.4) f 4 = c 1 4 +c +c 3 (A.5) f 5 = c 1 3 +c (A.6) f 6 = c 1 (A.7) f 7 = c 1, and then solve (7) to (10). These long but rather straightforward calculations ield the 15 superintegrable potentials included in Table 1. Some of them are obviousl particular cases of others, but we listed them separatel to account for their additional integrals. Onl the third-order integrals are listed, some of them being trivial consequences of lower-order ones. With the eception of the harmonic oscillator, potentials that have first-order integrals are not listed here for the were alread presented in [1] with all their third-order integrals. The complete integrals of motion can be found in Appendi I. Man of these potentials were not known. The onl classical potentials among these are indicated with a. These are well-known superintegrable 6

7 potentials (see e.g. [9]), and all of them, ecept V i, are in fact quadraticall superintegrable. All the potentials are superintegrable in the quantum case. We therefore notice that classical nontrivial potentials can have man different quantum equivalents. The classical harmonic oscillator V a can be seen as a limiting case of the quantum potentials V a, V c and V d, and also, if we set α = h/ω, of V e, V f and V g, not to mention the similar potentials that can be obtained b permutations of and. The anisotropic harmonic oscillator with ratio 1 : 3 also admits man quantum deformations but, interestingl, the anisotropic oscillator with ratio 1 : does not admit such deformations. Notice also that if we want to deal with real potentials onl, α must be either real or purel imaginar in potentials V e, V f, V g and V k. Therefore these have as a classical limit harmonic oscillators with a > 0. All quantum superintegrable potentials reduce to classical ones when the classical limit is considered, sometimes in more than one wa. For eample, potentials V e, V f and V g give the free motion potential instead of the harmonic oscillator if α remains constant as h 0. These potentials all satisf the linear equations (11) and (1), and can be epressed as sums of simple superintegrable potentials. Let us now set A 10 = A 111 = A 01 = 0 so that (11) vanishes triviall. We ma also assume that V 1 does not take one of the forms (A.1) to (A.7), for we have alread worked these cases out. This is quite useful, for if we set = 1 in (6), we obtain for V 1 an equation of the same form as (11) with different coefficients. These coefficients must therefore vanish, so A 300 = A 01 = A 10 = 0. This is a significant simplification that allows us to restrict our attention, when considering equation (1), to the following three cases: i) V = a ; ii) V = a; iii) A 10 = A 01 = 0 Before we consider each case separatel, it is worth noticing that potentials of the form V = V 1 () that admit third-order integrals independent of and p should appear as solutions here, for the integral remains if we add a function V () to these potentials. These potentials were found in [1] and [16] to satisf equation h V 1 = 4V 1 3 g V 1 g 3, (13) and can therefore be written as V 1 = h P(), (14) where P() is the Weierstrass elliptic function. Since the variable plas no role here and these potentials admit integrals with leading-order terms proportional to p 3, these solutions will appear in cases (i) to (iii). 7

8 A. Case i: V = V 1 ()+a When V = V 1 ()+a and a 0, we find that A 300 = 0, and the following two equations must be satisfied: 0 = A 030 ( h V (3) 1 6(V 1 ) ) +γ 1 V 1 (15) 0 = A 10 ( h V (4) 1 4a(V 1 ) +6(V 1 ) 4a V 1 +8a ) +8aA 01 (a (V 1 ) V 1 )+4η(a V 1 ). (16) where γ 1 and η are arbitrar constants. When A 030 0, equation (15) is equivalent to (13) (up to a translation of V 1 to get rid of γ 1 ), hence its solutions are of the form (14). These potentials cannot satisf simultaneousl equation (16) for nontrivial parameters. This can be observed b epanding (13) in series around = 0 and substituting the result in (16). Therefore solutions given b A are of no special interest here. Let us now set A 030 = γ 1 = 0. Equation (16) can be greatl simplified. We assume that A 10 0, for otherwise equation (16) can be solved to give potential (III). Then b an appropriate translation of and V, we can get rid of the terms involving A 01 and η and finall divide b A 10 ; 0 = h V (4) 1 4a(V 1 ) +6(V 1 ) 4a V 1 +8a. (17) This equation admits a first integral, namel k = h (V 1 V 1 )+4( a 3V 1 ) V 1 +6V 1 +1a V 1 a 4. (18) Both (17) and (18) can be simplified b setting V 1 = W()+a /3. Then h W (4) = 1WW +1(W ) +bw +bw 1 6 b (19) with b = 8a 0 for (17), and k = 3 h (W W ) 18(W ) +(a +3W). (0) for the first integral. Equation (19) is well known. It is equivalent to equations (3.16) in [4] and (.17) in [0], which were obtained b nonclassical reduction of the Boussinesq Equation. It was also shown in [3] to be a nonclassical reduction of the Kadomtsev-Petviashvili equation. It has the Painlevé propert, and, when b 0, its solution, given in [5] (equation.88) ma be written in terms of the fourth transcendant function of Painlevé, namel W = h b 1P 4 (, b h ) 1 bp 4 (, b h ) 1 bp b 4(, h ) 1 6 (b + h K 1 hb 1 ) (1) 8

9 where b 1 ± b = ± 8a and P 4 (, b ) = P h 4 (, b,k h 1,K ) is the fourth transcendant function of Painlevé, and therefore satisfies equation P 4 (,α) = (P 4(,α)) P 4 (,α) 3α P 4(,α) 3 αp 4 (,α) ( α +K 1 )P 4 (,α)+ K P 4 (,α). K 1 and K are integration constants. The potential therefore reads () V(,) =a( + )+ h b 1P 4(, 8a h )+4aP 4(, 8a h ) +4aP 4 (, 8a h )+ 1 (3) 6 ( h K 1 + hb 1 ). This potential admits as special cases two anisotropic harmonic oscillators, V = a( + ) when K = 0 (and P 4 = 0), and V = a( /9+ ) when K 1 = 0 and K = 1/18 (and P 4 = /3), as well as all their quantum deformations that have the form V = a(p + ) +f(), that is potentials V d, V e, V j and V k (up to a permutation of and ). The constant term in the potential, ( h K 1 + hb 1 ), can be set to zero, but we will keep it in order to be able to write the quantum and classical integrals in a unified wa. In classical mechanics, the equation (18) with h = 0 admits a first integral, which reads c = (9V 1 a )(V 1 a ) 3 + k 4 k(v 1 a )(3V 1 +a ). We ma therefore write the solution for V 1 implicitl as c d +d(v 1 a )(3V 1 +a ) = (9V 1 a )(V 1 a ) 3, (4) where c and d are arbitrar constants. If c = d = 0, we find either the familiar anisotropic harmonic oscillators, or a potential obtained b joining at = 0 two halves of anisotropic harmonic oscillators with different ratios. Even though this potential does not have a continuous second derivative, it can be obtained as a limiting case of the famil of smooth superintegrable potentials (6). In the general case the potential ma be epressed as the root of a fourthorder polnomial with three arbitrar parameters (a, c and d), although we can set a = 1 and one of the other two coefficients to ±1 b scaling and V 1. Since equation(4) ma describe new two-dimensional potentials with bounded motion, it is worth studing for interesting special cases. Indeed, if we assume a > 0 and d 0, we can consider the case c = 18a 4 d 3 /7 and d = 4a d /7. Four solutions eist in that case, two of which are translated harmonic oscillators. The remaining two potentials are V = a + a ( ) d+5 ±4 ( d+ 9 ). (5) 9

10 Up to an additive constant, potentials (5) are equal to V = a + a 9 ( ) ± d+. (6) Those potentials are smooth interpolations between anisotropic harmonic oscillators with ratio 1 : 1 and 1 : 3. When d = 0 (and therefore c = d = 0) the are the junctions of the two halves of harmonic oscillators mentioned above. The integral of motion is similar in quantum and classical mechanics, and reads X = {L,p }+{a 3V 1 (),p } 1 a { h 4 V 1+(a 3V 1 )V 1,p }, (7) where V 1 is a solution to equation (4). Therefore the integral for (6) must be slightl modified to take into account the constant term we removed. B. Case ii: V = a +V 1 () Here we find again two equations: 0 = h A 10 V (3) 1 6A 10 (V1 ) +1a A 003 +γ 1 V 1, 0 = A ( ) 030 6(V1 ) + h V (3) a 1 4 A ( 10 6(V1 ) + V 1 aa 01 ((V 1) +V 1) γ V 1 +aγ 1. ) (8) This time wecannottreatthe equationsseparatel,but wecan, when A 10 0, translate V to annihilate γ 1. Then we can substitute the first equation in the second to get rid of A 030, and finall translate to get rid of A 01. We can then solve the remaining sstem, first b solving the second, linear equation, and then b substituting the result in the first one. The onl solution remaining is then potential V n = a + h. Hence we can set A 10 = 0, and therefore also A 003 = γ 1 = 0 (as V 1 b and a 0). In the second equation we can set γ = 0. If then A 030 = 0, we find V 1 = a/ which was alread classified. If A and A 01 = 0, though, the potential V 1 () is solution to h V 1 = 6V 1 +λ+k (9) If λ 0 we can set k = 0, and we find that V 1 can be epressed in terms of the first Painlevé transcendent: Comme λ est une constante arbitraire, il est plus pratique d crire ω 5 = λ h 4, with ω 5 = λ h 4. The integral of motion is V = a + h ω P 1 (ω), (30) 10

11 X = p 3 +3{V (),p }+{ h4 ω 5 4a,p }. Notice that in order to consider the limiting case a = 0, we can multipl the integral b a, and thus we find a potential that depends onl on with a trivial p integral. We can also set ω 5 = aω 5 first, in which case we find the potentials (14). If we look for the classical limit, we find with the integral V = a +b. (31) X = p 3 +3b{,p } { 3b a,p }. Returning to equation (9) and assuming λ = 0, we find again that V 1 has the form (14), and the integral depends onl on and p. Finall, if we set A 030 A 01 0, we have to solve the equation 0 = ( 3 ) (V 1 ) + h 4 V (3) 1 +b((v 1 ) +V 1 ), b = aa (3) A 030 It can be integrated once to give C 1 = 3 (V 1 ) + h 4 V (3) 1 +b((v 1 )+V 1). We can set C 1 = 0 b translations in and V. The equation admits a first integral which reads b h (V 1 () b)v 1 ()+b h (b V 1())V 1() 8bV 1 ()(V 1 () b) = k 1. (33) In the classical case, h = 0, this is enough to solve. The solution, which can be written implicitl in a more compact form, is given b d = V 1 (V 1 b), (34) where d is an arbitrar constant. When d = 0 we find the familiar case V = b + a which admits a first-order integral. We ma notice that the implicit form of the solution is somewhat similar to (4). In the quantum case, we notice that the transformation W() = V 1 () b, that preserves the Painlevé propert, simplifies equation (33) that becomes h ( WW W ) 8(W +b)w = k. (35) or W = W W + 4W h + 4bW h + k h W (36) 11

12 We can also substitute Y() = W() to find h Y = ( Y +b ) Y + k 4Y 3. If k = 0, we can set h = 1 and b = 1/ b the change of variables Y = ( hb) 1 3 Z, = ( h The solution for Z(ξ) is then a special case of the second Painlevé transcendant, defined b the equation The solution for V is b ) 1 3 P (,α) = P (,α) 3 +P (,α)+α. V(,) = b+a +( hb) 3 P ( (b h ξ. )1 3,0 ), (37) with the particular case V = b+a. Now if k 0, we use (36), that can be normalized b a change of variables, We then find ( h ) 1 3 = (4b) k ξ, W = Y. (4 hb) 1 3 Y = Y Y +4βY ξy 1 Y, (38) where β = k /(4 hb) 0 is an arbitrar constant. Equation (38) corresponds to case XXXIV, p. 340, in [17]. The solution for Y(ξ) reads βy = P (ξ, β 1 )+P (ξ, β 1 ) + ξ (39) wherep is onceagainthe second Painlevétranscendant. Since β isan arbitrar constant, we can set κ = β 1. Back to the original variables, we get V(,) = a + ( h b ) 1 ) 3 (P ( (4b/ h ) 1 3,κ)+P ( (4b/ h ) 1 3,κ) (40) This potential admits V = a and V = a + h / as particular cases. In all cases with V = a+v 1 () and A 01 A the integral of motion is ap 3 bp p +a{3v 1 () b,p } b{v 1 (),p }. (41) Limiting values for a, b and h give either known potentials or trivial integrals. 1

13 C. Case iii: A 10 = A 01 = 0 Since we set here A 10 = A 01 = 0, we assume A 030 A (otherwise we would find potentials of the form (14)). The conditions for the eistence of a third-order operator then read h V 1 () = 6V 1 ()+A 003 σ h V () = 6V () A 030 σ. (4) If σ vanishes, we find the potential V = h (P()+P()), with two integrals that each depend on onl one variable, as could have been predicted from the results of equation (13). If σ does not vanish, let us set b 1 = A 003 σ and b = A 030 σ. If h = 0, we find that V = ± β 1 ± β is superintegrable, where the β i = b i /6 are arbitrar constants. This defines similar real potentials on each quadrant of the plane. We can again patchwork the pieces to find real continuous potentials defined everwhere, e.g. V = c 1 +c, although in that case neither the Hamiltonian nor the integral are differentiable at the equilibrium point. Finall, if h 0, we find that V 1 and V can be both written using the first Painlevé transcendent. V = h ω 1P 1 (ω 1 )+ h ω P 1 (ω ) (43) where the ω i = (b i / h 4 ) 1/5 are arbitrar constants. The integral of motion is X = b p 3 b 1 p 3 +3b {V 1 (,b 1 ),p } 3b 1 {V (,b ),p }, (44) both in quantum and classical mechanics. IV Conclusion We have found all sstems in two-dimensional euclidian space that admit separation of variables in cartesian coordinates and at least one third-order integral, both in quantum and classical mechanics. Man new superintegrable potentials were found, and, interestingl, all the quantum superintegrable potential are found as solutions of equations having the Painlevé propert, and this is probabl not accidental. Man of the quantum potentials can in fact be written in terms of different transcendent functions of Painlevé, and man are related to group invariant solutions of the Korteweg-de Vries equation, and to reductions 13

14 of the Boussinesq and KP equations. All classical integrable potentials were found to be limiting cases of quantum ones. The do not obe to equations having the Painlevé propert, though, for man classical integrable potentials havemovable branch points of the form b. Thus, in that respect, quantum integrable potentials behave more regularl than classical ones. A natural question is what does the Painlevé propert, in quantum mechanics, tell us about classical integrable potentials. Since this paper deals mostl with the classification of superintegrable sstems, the consideration of such questions regarding their properties and their solutions will be postponed to a future article. Our investigation provided interesting new eamples of the differences between quantum and classical integrabilit. A sstematic search for sstems with higher-order integrals is therefore a useful task since we still know little about such sstems, and the are likel to share interesting properties. Note added in proof: The future article mentioned earlier([11]) indeed shows that superintegrabilit of separable sstems tells a lot about the phsical properties of the separated sstems. Even though the sstems considered here can all be separated in two independent one-dimensional sstems, two-dimensional superintegrabilit provides information on them that is far from obvious from a one-dimensional perspective. Equations () to (5) are the necessar and sufficient conditions for the eistence of third-order integrals. The general solution to these equation is highl nontrivial, and it is not likel that a direct approach will lead to such a solution for integrals of order higher than four, without the use of new methods. In order to develop such methods, it would be useful to have a rigorous definition of quantum integrabilit. Acknowledgements The author thanks gratefull P. Winternitz for his ver helpful comments and suggestions, and also J. Bérubé for useful discussions. The author benefits from a NSERC ES-A fellowship. V Anne I Here is the complete list of two-dimensional Hamiltonians separable in cartesian coordinates that admit at least one third-order integral. We also gave all their third-order integrals. Potentials depending on onl one cartesian coordinates are not listed here, for the were classified in [1]. Man of the potentials listed below were alread known to be superintegrable, but we listed them here for completeness and in some cases to take into account additional third-order integrals. Potentials in boes are those who do not admit enough first- or second-order integrals to make them superintegrable. Most of them were not known before. Some potentials (namel V p and V q ) are not superintegrable according to 14

15 the usual definition, for their integrals have simple though nontrivial relations with the Hamiltonian. For definitions of parameters the reader is referred to the bod of the article. In order to write the integrals of motion in a compact form we often use the notation V(,) = V 1 ()+V () = V 1 +V. A. Quantum potentials (Q.1) V = a( + ) X 1 = L 3 X = {L,p p }+a{,p } a{,p } X 3 = {L,p }+a( {,p } { 3,p } ) X 4 = {L,p } a ( {,p } { 3,p } ) (Q.) V = a( + )+ b + c X 1 = {L,p p }+{ ( b 3 +a ),p } { ( ) c +a,p 3 } (Q.3) V = a( + )+ h + h X 1 = L 3 h { ,p }+ h { ,p } ( ) ( ) X = {L,p p }+{ h +a,p 3 } { h +a,p 3 } (Q.4) V = a( + )+ h X 1 = L 3 h { + 3,p }+ h { 33 +,p } X = {L,p p } {a h,p }+a{,p } X 3 = {L,p } {a 3 + h 1,p }+{ h 3 +a,p } ( ) 1 (Q.5) V = h 8α 4( + 1 )+ ( α) + 1 (+α) X 1 = L 3 3α {L,p h 3 }+ 4 {( 8+ α 4 ( +α ) ( α) (+α) ),p } + h 4 {(8 3 ( 4 10α 4α 4 ) α ( α) (+α) ),p } X = {L,p }+ h {( 4α 4α 6( +α )) 4 ( α ) ),p } + h { ( 4α ) 4α 4 α + 4( +α ) ( α) (+α),p } ( 1 (Q.6) V = h 8α 4( + )+ 1 ) + 1 (+α) + 1 ( α) X 1 = L 3 3α {L,p }+ h { 33 4α + 63 ( +α ) ( α) (+α) 3( α ),p }+ 3 h {( 3α 3 8α 1α α + 4 ( +α ) ( α) (+α) ),p } ( ) 1 (Q.7) V = h 8α 4( + 1 )+ ( α) + 1 ( α) + 1 ( +α) + 1 (+α) X 1 = L 3 3α ({L,p }+{L,p h })+ 4 {(14+ 3( + ) α + 4( 5 ) α 15

16 144 α + 4(3 )( +α ) ( α) (+α) + 48( )( +α ) ( α) (+α) 4(5 ) α ),p } h 4 {(14+ 3( + ) α 144 α 4( 3 )( +α ) ( α) (+α) + 48( )( +α ) ( α) (+α) ),p } (Q.8) V = a(4 + )+ b +c X 1 = p p +{ a + b,p }+{8a +c,p } (Q.9) V = a(9 + ) X 1 = {L,p }+ 3 a{3,p } 6a{,p } (Q.10) V = a(9 + )+ h X 1 = {L,p }+{a3 3 h,p }+{3 ( 1 (Q.11) V = h 8α 4(9 + )+ X 1 = {L,p }+ h { α 4 ),p } (Q.1) V = h + a ( ) a + h,p } ) 1 ( +α) + 1 ( α) ( ) 1α 8α 4 ( α ) α,p }+ 3 h +α ) 4 {(8( ( α ) X 1 = {L,p } h {,p }+{3 h +a + h,p } X = {L,p p } h {,p }+a{,p } X 3 = p 3 +{3 h,p } (Q.13) V = h + h X 1 = L 3 h { ,p }+ h { ,p } X = {L,p } h {,p }+{3 h + h + h,p } X 3 = {L,p } h {,p }+{3 h + h + h,p } X 4 = {L,p p } h {,p }+ h {,p } X 5 = p 3 +{ 3 h,p } X 6 = p 3 +{3 h,p } (Q.14) V = a+ h X 1 = {L,p } { h,p }+{3 h a,p } X = p 3 +{3 h,p } X 3 = p p + h { 1,p }+a{,p } 16

17 (Q.15) V = h P()+V() X 1 = p 3 +{3 h P(),p } (Q.16) V = h (P()+P()) X 1 = p 3 +{3 h P(),p } X = p 3 +{3 h P(),p } (Q.17) V = h ω 1P 1 (ω 1 )+ h ω P 1 (ω ) X 1 = ω 5 p 3 ω 5 1p 3 +3ω 5 {V 1 (,ω 1 ),p } 3ω 5 1{V (,ω ),p } (Q.18) V = a( + )+ h b 1P 4 8a 8a (, h )+4aP 4 (, h )+4aP 4(, 8a h )+ 1 6 ( h K 1 + hb 1 ) X 1 = {L,p }+{a 3V 1,p } 1 a { h 4 V 1 +(a 3V 1 )V 1,p } (Q.19) V = a + h ω P 1 (ω) X 1 = p 3 +3 h ω {P 1 (ω),p }+{ ω5 h 4 4a,p } (Q.0) V = b+a +( hb) 3 P ( (b h )1 3,0 ) X 1 = ap 3 bp p +a{3v 1 () b,p } b{v 1 (),p } (Q.1) V = a + ( h b ) 1 ) 3 (P ( (4b/ h ) 1 3,κ)+P ( (4b/ h ) 1 3,κ) X 1 = ap 3 bp p +a{3v 1 () b,p } b{v 1 (),p } B. Classical potentials (C.1) V = a( + ) X 1 = L 3 X = {L,p p }+a{,p } a{,p } X 3 = {L,p }+a( {,p } { 3,p } ) X 4 = {L,p } a( {,p } { 3,p } ) (C.) V = a( + )+ b + c X 1 = {L,p p }+{ ( b 3 +a ),p } { (C.3) V = a(4 + )+ b +c X 1 = p p +{ a + b,p }+{8a+c,p } (C.4) V = a(9 + ) X 1 = {L,p }+ 3 a{3,p } 6a{,p } ( ) c +a,p 3 } 17

18 (C.5) V = ± β 1 ± β X 1 = β p 3 β 1p 3 +3β {V 1 (,6β 1 ),p } 3β 1 {V (,6 b ),p } (C.6) V = a +V 1 where V 1 satisfies equation (4) X 1 = {L,p }+{a 3V 1,p } 1 a {(a 3V 1 )V 1,p } (C.7) V = a +b X 1 = p 3 +3b{,p } { 3b a,p } (C.8) V = a +V 1 (), where (V 1 b) V 1 = d X 1 = ap 3 bp p +a{3v 1 () b,p } b{v 1 (),p } References [1] V.I. Arnold, Mathematical Methods in Classical Mechanics, (Springer,New York,1980). [] J. Bertrand, Théorème relatif au mouvement d un point attiré vers un centre fie, C.R. Ac. Sci. 77, (1873). [3] P. A. Clarkson and P. Winternitz, Nonclassical smmetr reductions for the Kadomstev-Petviashvili equation, Phsica D, 49, 57-7 (1991). [4] P. A. Clarkson and M.D. Kruskal, New similarit solutions of the Boussinesq equations, J. Math. Phs. 30, (1989). [5] C.M. Cosgrove, Higher-order Painlevé equations in the polnomial class I. Bureau smbol P, Studies in Applied Mathematics 104, 1-65 (000). [6] N.W. Evans, Superintegrabilit in classical mechanics, Phs. Rev. A41, (1990). [7] N.W. Evans, Group theor of the Smorodinsk-Winternitz sstem, J. Math. Phs. 3, (1991). [8] N.W. Evans, Superintegrabilit of the Winternitz sstem, Phs. Lett. A147, (1990). [9] I. Friš, V. Mandrosov, J. Smorodinsk, M. Uhlíř and P. Winternitz, On higher-order smmetries in quantum mechanics, Phs. Lett. 16, (1965). [10] H. Goldstein, Classical Mechanics, (Addison-Wesle, Reading, MA,1990). [11] S. Gravel, Superintegrabilit, isochronicit, and quantum harmonic behavior, Ariv: math-ph/ (003) 18

19 [1] S. Gravel et P. Winternitz, Superintegrable sstems with third-order integrals in classical and quantum mechanics, J. Math. Phs, 43, (00). [13] J.Hietarinta, Pure quantum integrabilit, Phs. Let. A46, (1998). [14] J.Hietarinta, Classical vs quantum integrabilit, J. Math. Phs. 5, (1984). [15] J.Hietarinta and B. Grammaticos, On the h -correction in quantum mechanics, J. Phs. A: Math. Gen., (1989). [16] J.Hietarinta, Solvabilit in quantum mechanics and classicall superfluous invariant, J. Phs. A: Math. Gen., L143-L147 (1989). [17] E.L. Ince, Ordinar Differential Equations, (Dover, New York, 1944) [18] E.G. Kalnins, W. Miller Jr and G.S. Pogosan, Completeness of multiseparabilit in two-dimensional constant curvature space, J. Phs. A34, (001). [19] E.G. Kalnins, J.M. Kress and P. Winternitz, Superintegrabilit in a twodimensional space of non-constant curvature, J. Math. Phs. 43, (00). [0] D. Levi and P. Winternitz, Non-classical smmetr reduction: eample of the Boussinesq equation, J. Phs. A, (1989) [1] A.Makarov, J. Smorodinsk, Kh. Valiev and P. Winternitz, A sstematic search for non-relativistic sstems with dnamical smmetries, Nuovo Cim. A5, (1967). [] M.F. Rañada, Superintegrable n = sstems, quadratic constants of motion and the potentials of Drach, J. Math. Phs. 38, (1997). [3] M.F. Rañada et M. Santander, Comple Euclidian super-integrable potentials, potentials of Drach and potentials of Holt, J. Phs. Lett. A, 78, (001). [4] M.F. Rañada and M. Santander, Superintegrable sstems on the twodimensional sphere S and the hperbolic plane H, J. Math. Phs. 40, (1999). [5] M. B. Sheftel,P. Tempesta and P. Winternitz, Superintegrable sstems in quantum mechanics and classical Lie theor, J. Math. Phs. 4, (001). [6] P. Tempesta, A.V. Turbiner and P. Winternitz, Eact solvabilit of superintegrable sstems, J.Math.Phs. 4, (001). 19

20 [7] A. V. Tsiganov, The Drach superintegrable sstems, J. Phs. A: Math. Gen. 33, (000). [8] S. Weigert, The problem of quantum integrabilit, Phsica D 56, (199). [9] P. Winternitz, J. Smorodinsk, M. Uhlíř and I. Friš, Smmetr groups in classical and quantum mechanics, Yad. Fiz. 4, (1966) [Sov. J. Nucl. Phs. 4, (1967)]. 0