arxiv: v1 [math-ph] 31 Jan 2015

Symmetry, Integrability and Geometry: Methods and Applications SIGMA? (00?), 00?,?? pages Structure relations and Darboux contractions for D nd order superintegrable systems R. Heinonen, E. G. Kalnins, W. Miller, Jr. and E. Subag arxiv:150.0018v1 [math-ph] 31 Jan 015 School of Mathematics, University of Minnesota, Minneapolis, Minnesota 55455, USA E-mail: heino07@umn.edu Department of Mathematics, University of Waikato, Hamilton, New Zealand E-mail: math036@waikato.ac.nz School of Mathematics, University of Minnesota, Minneapolis, Minnesota 55455, USA E-mail: miller@ima.umn.edu URL: http://www.ima.umn.edu/ miller/ School of Mathematical Science, Tel Aviv University, Tel Aviv 69978, Israel E-mail: eyalsubag@gmail.com Received???, in final form????; Published online???? Original article is available at http://www.emis.de/journals/sigma/00?/00?/ Abstract. Two-dimensional quadratic algebras are generalizations of Lie algebras that include the symmetry algebras of nd order superintegrable systems in dimensions as special cases. The superintegrable systems are exactly solvable physical systems in classical and quantum mechanics. Distinct superintegrable systems and their quadratic algebras can be related by geometric contractions, induced by Wigner-Inönu type Lie algebra contractions. These geometric contractions have important physical and geometric meanings, such as obtaining classical phenomena as limits of quantum phenomena as 0 and nonrelativistic phenomena from special relativistic as c, and the derivation of the Askey scheme for obtaining all hypergeometric orthogonal polynomials as limits of Racah/Wilson polynomials. In this paper we show how to simplify the structure relations for abstract nondegenerate and degenerate quadratic algebras and their contractions. In earlier papers we have classified all contractions of nd order superintegrable systems on constant curvature spaces and have shown that all results are derivable from free quadratic algebras contained in the enveloping algebras of the Lie algebras e(, C) in flat space and o(3, C) on nonzero constant curvature spaces. The quadratic algebra contractions are induced by generalized Inönü-Wigner contractions of these Lie algebras. As a special case we obtained the Askey scheme for hypergeometric orthogonal polynomials. After constant curvature spaces, the 4 Darboux spaces are the D manifolds admitting the most nd order Killing tensors. Here we complete this theoretical development for D superintegrable systems by showing that the Darboux superintegrable systems are also characterized by free quadratic algebras contained in the symmetry algebras of these spaces and that their contractions are also induced by generalized Inönü-Wigner contractions. We present tables of the contraction results. Key words: contractions; quadratic algebras; superintegrable systems; Darboux spaces; Askey scheme 000 Mathematics Subject Classification: E70; 16G99; 37J35; 37K10; 33C45; 17B60 1 Introduction We define an abstract nondegenerate (quantum) quadratic algebra as a noncommutative associative algebra generated by linearly independent operators H, L 1, L such that H is in the center,

R. Heinonen et al. R = [L 1, L ] 0 and the following relations hold: [L j, R] = M e (j) 1,e,e 3 {L e 1 1, Le }He 3, e k 0, L 0 k = I, (1) 0 e 1 +e +e 3 for some M e (j) 1,e,e 3 C, and where [A, B] = AB BA is the commutator and {L 1, L } = L 1 L +L L 1 is the symmetrizer. Also the operator R is contained in the algebra of symmetrized products: R F R 0 e 1 +e +e 3 3 N e1,e,e 3 {L e 1 1, Le }He 3 = 0 () for some N e1,e,e 3 C. An abstract degenerate (quantum) quadratic algebra is a noncommutative associative algebra generated by linearly independent operators X, H, L 1, L such that H is in the center and the following relations hold: [X, L j ] = 0 e 1 +e +e 3 +e 4 1 P (j) e 1,e,e 3,e 4 L e 1 1 Le He 3 X e 4, j = 1, (3) for some P e (j) 1,e,e 3,e 4 C. The commutator [L 1, L ] is expressed as [L 1, L ] = Q e1,e,e 3,e 4 {L e 1 1 Le, X}He 3 X e 4 (4) 0 e 1 +e +e 3 +e 4 1 for some Q e1,e,e 3,e 4 C. Finally, there is the relation: G S e1,e,e 3,e 4 {L e 1 1, Le, Xe 4 }H e 3 = 0, X 0 = H 0 = I, (5) 0 e 1 +e +e 3 +e 4 for some S e1,e,e 3,e 4 C and where {L e 1 1, Le, Xe 4 } is the 6-term symmetrizer of three operators. For both quantum quadratic algebras there is a natural grading such that the operators H, L j are nd order, X is 1st order and ord ([A, B]) ord (A) + ord (B) 1, ord (AB) = ord (A) + ord (B), (6) ord (I) = 0, ord (A + B) = max {ord (A), ord (B)}, ord (ca) = ord (A), for operators A, B, identity operator I and scalar c, with A B, A, B 0, and c 0. Thus R is usually 3rd order, expression G is 4th order and F is 6th order. The field of scalars can be either R or C. There is an analogous quadratic algebra structure for Poisson algebras. An abstract nondegenerate (classical) quadratic algebra is a Poisson algebra with functionally independent generators H, L 1, L such that all generators are in involution with H and the following relations hold: {L j, R} = R F R 0 e 1 +e +e 3 0 e 1 +e +e 3 3 M e (j) 1,e,e 3 L e 1 1 Le He 3, e k 0, L 0 k = 1. (7) N e1,e,e 3 L e 1 1 Le He 3 = 0 (8) for some M e (j) 1,e,e 3, N e1,e,e 3 C. An abstract degenerate (classical) quadratic algebra is a Poisson algebra with linearly independent generators X, H, L 1, L such that all generators are in involution with H and obey structure equations {X, L j } = 0 e 1 +e +e 3 +e 4 1 P (j) e 1,e,e 3,e 4 L e 1 1 Le He 3 X e 4, j = 1,. (9)

Structure relations and Darboux contractions 3 {L 1, L } = G 0 e 1 +e +e 3 +e 4 1 0 e 1 +e +e 3 +e 4 Q e1,e,e 3,e 4 L e 1 1 Le X He 3 X e 4. (10) S e1,e,e 3,e 4 L e 1 1 Le X e 4 H e 3 = 0, X 0 = H 0 = 1 (11) for some P (j) e 1,e,e 3,e 4, Q e1,e,e 3,e 4, S e1,e,e 3,e 4 C. There is a grading for these quadratic algebras with properties analogous to (6), but with the Poisson bracket instead of the commutator. These quadratic algebra structures arise naturally in the study of classical and quantum superintegrable systems in two dimensions and are key to the exact solvability of these systems, e.g., [1 9]. A quantum D superintegrable system is an integrable Hamiltonian system on an -dimensional Riemannian/pseudo-Riemannian manifold with potential that admits 3 algebraically independent partial differential operators commuting with H, the maximum possible. H = + V, [H, L j ] = 0, L 3 = H, j = 1,, 3. (In dimensions we can always find Cartesian-like coordinates x 1, x such that H = 1 λ(x 1,x ) ( x 1 + x ) + V (x 1, x ) and we adopt these coordinates in the following.) A system is of order k if the maximum order of the symmetry operators L j (other than H) is k; all such systems are known for k = 1,, [10,11]. Superintegrability captures the properties of quantum Hamiltonian systems that allow the Schrödinger eigenvalue problem HΨ = EΨ to be solved exactly, analytically and algebraically. A classical D superintegrable system is an integrable Hamiltonian system on an - dimensional Riemannian/pseudo-Riemannian manifold with potential that admit 3 functionally independent phase space functions H, L 1, L in involution with H, the maximum possible: H = p 1 + p λ(x) + V (x), {H, L j } = 0, L 3 = H, j = 1,, 3, expressed in local Cartesian-like coordinates x 1, x, p 1, p. A system is of order k if the maximum order of the constants of the motion L j, j 3, as polynomials in p 1, p is k. Again all such systems are known for k = 1,, and there is a 1-1 relationship between classical and quantum nd order D superintegrable systems, [11] The possible superintegrable systems divide into four classes: 1. First order systems. These are the (zero-potential) Laplace-Beltrami eigenvalue equations on constant curvature spaces. The symmetry algebras close under commutation to form the Lie algebras e(, R), e(1, 1), o(3, R) or o(, 1). Such systems have been studied in detail, using group theory methods, e.g. [1, 13].. Free triplets. These are superintegrable systems with zero potential and all generators of nd order. The possible spaces for which these systems can occur were classified by Koenigs [10].They are: constant curvature spaces (6 linearly independent nd order symmetries, 3 1st order), the four Darboux spaces (4 linearly independent nd order symmetries, 1 1st order), and eleven 4-parameter Koenigs spaces (3 linearly independent nd order symmetries, 0 1st order). In most cases the symmetry operators will not generate a quadratic algebra, i.e., the algebra will not close. If the system generates a nondegenerate quadratic algebra we call it a free quadratic triplet. 3. Nondegenerate systems, [11]. Here all symmetries are of nd order and the space of potentials is 4-dimensional: V (x) = a 1 V (1) (x) + a V () (x) + a 3 V (3) (x) + a 4 The symmetry operators generate a nondegenerate quadratic algebra with parameters a j.

4 R. Heinonen et al. 4. Degenerate systems, [14] There are 4 generators: one 1st order X and 3 second order H, L 1, L. Here, X is not contained in the span of H, L 1, L. The space of potentials is -dimensional: V (x) = a 1 V (1) (x) + a. The symmetry operators generate a degenerate quadratic algebra with parameters a j. Relation (11) is an expression of the fact that 4 symmetry operators cannot be algebraically independent. Every degenerate superintegrable system occurs as a restriction of the 3-parameter potentials (i.e., 4-dimensional potential space) to 1-parameter ones, such that one of the symmetries becomes a perfect square: L = X. Here X is a first order symmetry and a new nd order symmetry appears so that this restriction admits more symmetries than the original system. Strictly speaking, a nondegenerate D superintegrable system, both classical and quantum, is not a single system but in fact a family of superintegrable systems parameterized by three parameters, a 1, a, a 3. Similarly a degenerate D superintegrable system, both classical and quantum, is a family of superintegrable systems parameterized by one parameter, a 1. For a quadratic algebra that comes from a nondegenerate D superintegrable system (classical and quantum) the constants M (j) e 1,e,e 3 and N e1,e,e 3 are polynomials in the parameters a 1, a, a 3 of degree e 1 e e 3 and 3 e 1 e e 3, respectively. If all parameters a j = 0 the algebra is free. For a quadratic algebra that comes from a degenerate D superintegrable system (classical and quantum) the constants P (j) e 1,e,e 3,e 4, Q e1,e,e 3,e 4 and S e1,e,e 3,e 4 are polynomials in a 1 of degrees 1 e 1 e e 3 e 4, 1 e 1 e e 3 e 4 and e 1 e e 3 e 4, respectively. If all parameters a j = 0 these algebras are free. Basic results that relate these superintegrable systems are the Closure Theorems: Theorem 1. A free triplet, classical or quantum, extends to a superintegrable system with potential if and only if it generates a free quadratic algebra Q. Theorem. A superintegrable system, degenerate or nondegenerate, classical or quantum, with quadratic algebra Q, is uniquely determined by its free quadratic algebra Q. These theorems were proved in [15], except for systems on the Darboux spaces which will be proved in this paper. The proofs are constructive: Given a free quadratic algebra Q one can compute the potential V and the symmetries of the quadratic algebra Q. Thus as far as superintegrable systems on specific spaces are concerned, all information about the systems is contained in the free quadratic algebras. We will refer to quadratic algebras associated with superintegrable systems as geometric, a subset of abstract quadratic algebras. For quadratic algebras associated with quantum superintegrable systems the order of A is its order as a partial differential operator. For quadratic algebras associated with classical superintegrable systems the order of A is its order as a polynomial in the momenta. Although there is a 1-1 relationship between classical and quantum geometric systems the corresponding classical and quantum geometric quadratic algebras are not the same; only the highest order terms in the structure equations agree. Contractions The notion of contractions for quadratic algebras is based on that for Lie algebras, e.g. [16, 17]: Definition 1. Let g be a complex Lie algebra with an underlying vector space V and Lie brackets [, ]. In the following we simply write it as g = (V, [, ]). Suppose that for any ɛ (0, 1], t ɛ : V V is a a linear invertible operator and that lim ɛ 0 + t 1 ɛ [t ɛ X, t ɛ Y ] converge for any X, Y V. We use the notation lim ɛ 0 + t 1 ɛ [t ɛ X, t ɛ Y ] = [X, Y ] 0

Structure relations and Darboux contractions 5 Then [, ] 0 are in fact Lie brackets on V and we denote this Lie algebra by g 0 = (V, [, ] 0 ). We say that g 0 is a contraction of g (that is realized by the family of linear maps {t ɛ } ɛ (0,1] ) and we denote it by g g 0. Thus, as ɛ 0 the 1-parameter family of basis transformations can become singular but the structure constants go to a finite limit. Of particular interest to us are contractions that are induced by ɛ-dependent local analytic coordinate transformations x j (ɛ, x 1,, x n), j = 1,, n on a manifold M such that the Jacobian det( x x ) 0 for ɛ (0, 1], but the Jacobian is undefined or nonsingular in the limit as ɛ 0. Definition. Suppose that g = (V, [, ] 0 ) is a contraction of g = (V, [, ]) that is realized by the family of linear maps {t ɛ } ɛ (0,1]. Let M be a smooth manifold with a local coordinate system (x 1, x,..., x n ) and let ψ : g C (T M) be an embedding of Lie algebras where C (T M) is the space of smooth functions on the cotangent bundle (the phase space) of M which is equipped with its canonical Poisson brackets as Lie brackets. Suppose that x j (ɛ, x 1,, x n), j = 1,, n are ɛ-dependent local analytic coordinate transformations such that the Jacobian det( x x ) 0 for ɛ (0, 1]. If the limit lim ψ(t ɛ(v))(x 1 (ɛ, x 1,, x n),..., x n (ɛ, x 1,, x n)) (1) ɛ 0 + converge for any v V and defines a Lie algebra homomorphism from g 0 into a local expression of a certain space of C (T M ) for some smooth manifold M. Then we say that the contraction g g 0 is implemented by x j (ɛ, x 1,, x n), j = 1,, n and call this procedure a geometric Lie algebra contraction. We give some pertinent examples. Example 1. Consider the complex three dimensional Lie algebra G3 defined by basis elements {P 1, P, D} that satisfy [P 1, P ] = 0, [P 1, D] = P 1, [P, D] = P. This algebra admits an abstract Inönü-Wigner contraction defined by t ɛ (P 1 ) = ɛp 1, t ɛ (P ) = ɛp, t ɛ (D) = D. In this case the contracted Lie algebra, G3 0, coincide with G3 and [P 1, P ] 0 = 0, [P 1, D] 0 = P 1, [P, D] 0 = P. Now considering the complex analytic manifold C with coordinates (x 1, x ) we can realize G3 by ψ(p 1 )(x 1, x, p x1, p x ) = p x1, ψ(p )(x 1, x, p x1, p x ) = p x, ψ(d)(x 1, x, p x1, p x ) = x 1 p x1 +x p x. A geometric implementation of the contraction is obtained by the substitution x 1 = ɛx 1, x = ɛx. Then, taking the limit, we find lim ψ(t ɛ(p 1 ))(ɛx 1, ɛx, p ɛx ɛ 0 + 1, p ɛx ) = lim ψ(t ɛ(p ))(ɛx 1, ɛx, p ɛx ɛ 0 + 1, p ɛx ) = lim ψ(t ɛ(d))(ɛx 1, ɛx, p ɛx ɛ 0 + 1, p ɛx ) = lim ɛp ɛ 0 + ɛx = 1 lim ɛ 0 + ɛ1 ɛ p x 1 = p x 1 (13) lim ɛ 0 + ɛx = lim ɛ 0 + ɛ p x = p x (14) 1 lim ɛ 0 ɛx + 1 ɛ p x + 1 ɛx 1 ɛ p x = x 1p x 1 + x p x (15) Though this Lie algebra contraction acts like the identity map here, we shall see that its action on Darboux quadratic algebras is nontrivial. Example. We again consider the Lie algebra G3 with the same geometric identification. This algebra admits another abstract Inönü-Wigner contraction defined by t ɛ (P 1 ) = P 1, t ɛ (P ) = P, t ɛ (D) = ɛd. The contracted Lie algebra, G3 0, is given by [P 1, P ] 0 = 0, [P 1, D] 0 =

6 R. Heinonen et al. 0, [P, D] 0 = 0. A geometric implementation of the contraction is obtained by the substitution x 1 = x 1 + 1 ɛ, x = x. Then, taking the limit, we find lim ψ(t ɛ(p 1 ))(x ɛ 0 + 1 + 1 ɛ, x, p x 1 + 1, p x ɛ ) = lim p ɛ 0 + x 1 + 1 ɛ lim ψ(t ɛ(p ))(x ɛ 0 + 1 + 1 ɛ, x, p x 1 + 1, p x ɛ ) = lim ɛ 0 + p x = p x lim ψ(t ɛ(d))(p ))(x ɛ 0 + 1 + 1 ɛ, x, p x 1 + 1, p x ɛ ) = lim ɛ 0 + ɛ = lim ɛ 0 + p x 1 = p x 1 ( (x 1 + 1 ) ɛ )p x 1 + 1 + x p x ɛ (16) (17) (18) = p x 1 Hence the resulting map from G3 0 that is given by P 1 p x 1, P p x, D 1 p x 1 is a Lie algebra homomorphism with kernel spanned by P 1 D. Though this geometric implementation is not an isomorphism of the contracted Lie algebra we shall see that its action on Darboux quadratic algebras is isomorphic. There are exactly analogous implementations of geometric contractions in the quantum case. Definition 3. Algebraic Contraction of quadratic algebras. Let A be a classical nondegenerate quadratic algebra with a generating set which consist of a Hamiltonian, H (second order element which lies in the center of A) and two second order constants of motion L 1,L. Let {L 1, L } = R = F(H, L 1, L ) be the Casimir of A. For any ɛ (0, 1] consider a matrix A ɛ GL(3, C) of the form A ɛ = A 1,1 (ɛ) A 1, (ɛ) A 1,3 (ɛ) A,1 (ɛ) A, (ɛ) A,3 (ɛ) 0 0 A 3,3 (ɛ) GL(3). (19) Assume that the map from (0, 1] to GL(3, C) that is given by ɛ A ɛ is continuous. For any ɛ (0, 1] we have another set of generators, {L ɛ 1, Lɛ, Hɛ }, for A that is defined by L 1 L H = A ɛ L ɛ 1 L ɛ H ɛ and satisfies the Casimir relation (0) {L ɛ 1, L ɛ } = (R ɛ ) = R (A 11 (ɛ)a (ɛ) A 1 (ɛ)a 1 (ɛ)) = (1) L ɛ (A 11 (ɛ)a (ɛ) A 1 (ɛ)a 1 (ɛ)) 1 F A ɛ L ɛ = H ɛ α ijk (ɛ)(l ɛ 1) i (L ɛ ) j (H ɛ ) k i+j+k=3,0 i,j,k If lim ɛ 0 + α ijk (ɛ) exists for any i, j, k we denote it by α ijk (0). Then there exists a quadratic algebra, A 0 with a set of generators {L 0 1, L0, H0 } that satisfy {L 0 1, L 0 } = (R 0 ) = i+j+k=3,0 i,j,k we call A 0 the contraction of A with respect to {A ɛ } ɛ (0,1]. α ijk (0)(L 0 1) i (L 0 ) j (H 0 ) k ()

Structure relations and Darboux contractions 7 Note that we can expand {R ɛ, L ɛ 1 },{Rɛ, L ɛ } as quadratic expressions in Lɛ 1, Lɛ, Hɛ. For a contraction it might seem that we must also require these expansion coefficients to have finite limits as ɛ 0. However from the results of 3.1 the convergence of these other structure equations follows from the convergence of the Casimir. There is a completely analogous definition of contraction for quantum quadratic algebras. Just as for abstract classical and quantum Lie algebra contractions there are abstract classical and quantum quadratic algebra contractions that are induced by ɛ-dependent local analytic coordinate transformations x j (ɛ, x 1,, x n), j = 1,, n on a manifold M such that the Jacobian det( x x ) 0 for ɛ (0, 1], but the Jacobian is undefined or nonsingular in the limit as ɛ 0. If an algebraic contraction A B can be implemented by some coordinate transformation, we say that it is a geometric quadratic algebra contraction. In fact, all of the quadratic algebra contractions for Darboux systems discussed in this paper are geometric implementations. We will give many examples in the following sections. In [15] Lie algebra and quadratic algebra contractions for superintegrable systems on constant curvature spaces were related. Theorem 3. Every Lie algebra contraction of G = e(, C) or G = o(3, C) induces a geometric contraction of a free geometric quadratic algebra Q based on G, which in turn induces uniquely a contraction of the quadratic algebra Q with potential. This is true for both classical and quantum algebras. (Here we will demonstrate the analogous result for Darboux spaces, using the conformal symmetry algebra G3 with basis { x, y, x x + y y }.) 3 Structure relations Although the full sets of classical structure equations can be rather complicated, the function F contains all of the structure information for nondegenerate systems and G (only unique up to a nonzero scalar multiple) most of the information for degenerate systems. In particular, it is easy to show that [15], {L 1, R} = 1 F L, {L, R} = 1 F L 1, so the Casimir contains within itself all of the structure equations. For degenerate systems we have [15] {X, L 1 } = K G L, {X, L } = K G L 1, {L 1, L } = K G X, (3) where {H, K} = 0. Here, K is a scalar, unless {X, L 1 } and {X, L } are linearly dependent. In the latter case there would exist 3 algebraically independent elements of the algebra in involution, including the Hamiltonian. This is impossible for a Hamiltonian system. Thus, except for some abstract quadratic algebras unrelated to geometric superintegrable systems, K will always be a scalar that can be normalized to 1. 3.1 The quantum operators F and G The quantum case is similar to the classical case, but more complicated. relation From the Casimir R F R (b 1 L 3 1 + b L 3 + b 3 H 3 + b 4 {L 1, L } + b 5 {L 1, L } + b 6 L 1 L L 1 + b 7 L L 1 L (4) +b 8 H{L 1, L } + b 9 HL 1 + b 10 HL + b 11 H L 1 + b 1 H L + b 13 L 1 + b 14 L + b 15 {L 1, L } +b 16 HL 1 + b 17 HL + b 18 H + b 19 L 1 + b 0 L + b 1 H + b ) = 0,

8 R. Heinonen et al. we want to determine the structure relations. identities Noting that R = [L 1, L ] and using operator L RL = 1 [[L, R], L ] + 1 {R, L }, L 1 RL 1 = 1 [[L 1, R], L 1 ] + 1 {R, L 1}, L 1 RL + L RL 1 = 1 [L, [L 1, R]] 1 [L 1, [L, R]] + 1 {R, {L 1, L }}, (5) [L 1, R ] = {R, [L 1, R]}, [L, R ] = {R, [L, R]}, [L 1, L L 1 L ] = {[L 1, L ], 1 {L 1, L }}, XL 1 L + L L 1 X = 1 {X, {L 1, L }} + 1 [X, [L 1, L ]], (true formally for all operators L 1, L, R, not just for R the commutator), and setting [L 1, R] = A 1 L 1 + A L + A 3 H + A 4 {L 1, L } + A 5 HL 1 + A 6 HL + A 7 L 1 + A 8 L + A 9 H + A 10, [L, R] = B 1 L 1 + B L + B 3 H + B 4 {L 1, L } + B 5 HL 1 + B 6 HL + B 7 L 1 + B 8 L + B 9 H + B 10, we can write [L j, R F ] = 0 in the form {R, Q j } = 0 for some explicit differential operator Q j. This can only hold if Q j = 0. As a result we find: A 1 = b 4 + 1 b 6, A = 3 b, A 3 = 1 b 1, A 4 = b 5 + 1 b 7, A 5 = b 8, A 6 = b 10, A 7 = b 15 3 4 b 1b + b 4 b 5 1 4 b 6b 7, A 8 = b 14 1 b b 4 + b 5 b b 6 + 1 b 5b 7, A 9 = 1 b 17 1 4 b b 9 + 1 b 5b 8 1 4 b 6b 10, A 10 = 1 b 0 1 4 b b 13 1 4 b 1b b 5 1 4 b b 4 + 1 b 5b 15 + 1 b 4b 5 + 1 4 b 1b b 7 + 1 4 b b 6 1 4 b 5b 6 1 4 b 5b 6 b 7 1 4 b 6b 14, B 1 = 3 b 1, B = b 5 1 b 7, B 3 = 1 b 11, B 4 = b 4 1 b 6, B 5 = b 9, B 6 = b 8, B 7 = b 13 + 1 b 1b 5 b 4 + b 1 b 7 1 b 4b 6, B 8 = b 15 + 3 4 b 1b b 4 b 5 + 1 4 b 6b 7, B 9 = 1 b 16 + 1 4 b 1b 10 1 b 4b 8 + 1 4 b 7b 9, B 10 = 1 b 19 + 1 4 b 1b 14 + 1 4 b 1b b 4 + 1 4 b 1b 5 1 b 4b 15 1 b 4b 5 1 4 b 1b b 6 1 4 b 1b 7 + 1 4 b 4b 7 + 1 4 b 4b 6 b 7 + 1 4 b 7b 13. For quantum degenerate systems, in the Casimir relation G(L 1, L, H, X, α) = 0 (6) for nd order superintegrable systems with degenerate potentials we assume that G is given, up to a multiplicative factor, and set G = c 1 L 1 + c L + c 3 H + c 4 {L 1, L } + c 5 HL 1 + c 6 HL + c 7 X 4 + c 8 {X, L 1 } + c 9 {X, L } +c 10 HX + c 11 XL 1 X + c 1 XL X + c 13 L 1 + c 14 L + c 15 H + c 16 X + c 17, [X, L 1 ] = C 1 L 1 + C L + C 3 H + C 4 X + C 5, [X, L ] = D 1 L 1 + D L + D 3 H + D 4 X + D 5, [L 1, L ] = E 1 {L 1, X} + E {L, X} + E 3 HX + E 4 X 3 + E 5 X.

Structure relations and Darboux contractions 9 Using identities (5) we find [X, G] = {[X, L 1 ], A 1 c 11 4 (C 1 + C D 1 ) c 1 4 (C 1D 1 + D 1 D ) + c 13 } + {[X, L ], A + c 11 4 (C 1C + C D ) c 1 4 (C D 1 + D) + c 14 } = 0, (7) [L 1, G] = {[L 1, L ], A + c 8 4 (C 1C + C D ) c 9 4 (C 1 + C 1 D ) + c 1 4 (C 1D C D 1 ) + c 14 } + {[L 1, X], A 3 + k 1 X} + {[L, X], k X} = 0, (8) [L, G] = {[L 1, L ], A 1 + c 8 4 (C 1D + D) c 9 4 (C 1D 1 + D 1 D ) + c 11 4 (C D 1 C 1 D ) c 13 } + {[L 1, X], k 3 X} + {[L, X], A 3 + k 4 X} = 0, (9) where A 1 = c 1 L 1 + c 4 L + c 5 H + (c 8 + c 11 )X, A = c L + c 4 L 1 + c 6 H + (c 9 + c 1 )X, A 3 = c 7 X 3 + (c 8 + c 11 ){L 1, X} + (c 9 + c 1 ){L, X} + c 10 HX, k 1 = c 7 (C 1 + C D 1 ) + c 8 (C 1C 4 + C D 4 C E 1 ) + c 9 (C 4D 1 + C 1 E 1 C 1 D 4 ) + c 1 (C 1E 1 + C 1 D 4 + D 1 E C 4 D 1 ) + c 16, k = c 7 (C 1 C + C D ) + c 8 (C C 4 C E ) + c 9 (C 4D + C 1 C 4 + C 1 E ) + c 1 (C D 4 + C E 1 + D E C 4 D ), k 3 = c 7 (C 1 D 1 + D 1 D ) + c 8 (C 1D 4 + D D 4 D 1 E 1 ) + c 9 (D 1E 1 + D 1 D 4 ) + c 11 (C 4D 1 C 1 D 4 C 1 E 1 D 1 E ), k 4 = c 7 (C D 1 + D ) + c 8 (C D 4 C 4 D D E ) + c 9 (D D 4 + C 4 D 1 + D 1 E ) + c 11 (C 4D C D 4 C E 1 D E ) + c 16. Equating the coefficients of the 4th order terms in (7) and the coefficients of the 5th order terms in (8) and (9) we find [X, L 1 ] = KA + C 5, [X, L ] = KA 1 + D 5, [L 1, L ] = KA 3 + E 5 X, (30) C 1 = Kc 4, C = Kc, C 3 = K c 6, C 4 = K(c 9 + c 1 ), (31) D 1 = Kc 1, D = Kc 4, D 3 = K c 5, D 4 = K(c 8 + c 11 ), E 1 = K(c 8 + c 11 ), E = K(c 9 + c 1 ), E 3 = Kc 10, E 4 = Kc 7, for some constant K. Now we substitute these values back into (7), (8), and (9). We immediately see that k = k 3 = 0 and k 1 = k 4 = c 16 K c 7 (c 4 c 1 c ) + K (c 4 c 9 c c 8 )(c 8 + c 11 ) + K (c 4 c 8 c 1 c 9 )(c 9 + c 1 )

10 R. Heinonen et al. and we obtain ( c13 C 5 A 1 + D 5 A + c ) K 11 4 (c 4 c 1 c ) (KA + C 5 ) (3) ( c14 + c ) K 1 4 (c 4 c 1 c ) ( KA 1 + D 5 ) = 0, ( c14 c ) K 1 4 (c 4 c 1 c ) (KA 3 + E 5 X) C 5 (A 3 + k 1 X) (33) +(E 5 Kk 1 ){X, A } = 0, ( c13 c ) K 11 4 (c 4 c 1 c ) (KA 3 + E 5 X) D 5 (A 3 + k 1 X) (34) (E 5 Kk 1 ){X, A 1 } = 0. These equations split into terms of order 3,,1 and 0. From equation (3) we find C 5 = c 14 K c 1 4 K3 (c 4 c 1 c ), D 5 = c 13 K + c 11 4 K3 (c 4 c 1 c ) (35) except, possibly, for some degenerate cases. The condition that (35) is the unique solution of (3) is exactly that the set ([X, S 1 ], [X, S ]) is linearly independent. Otherwise the solution, though it always exists, is not unique. Substituting (35) into (33) and (34), we have C 5 (E 5 Kk 1 ) X + (E 5 Kk 1 ){X, A } = 0, K whence we find D 5 (E 5 Kk 1 ) X (E 5 Kk 1 ){X, A 1 } = 0, K E 5 = Kc 16 K 3 c 7 (c 4 c 1 c ) + K 3 (c 4 c 9 c c 8 )(c 8 + c 11 ) + K3 (c 4 c 8 c 1 c 9 )(c 9 + c 1 ). (36) We conclude in both the classical and quantum cases that the Casimirs of superintegrable systems determine the structure equations. 4 Free D nd order superintegrable systems As was shown in [11] the free nd order superintegrable system obtained by setting all the parameters in a nondegenerate potential equal to zero retains all of the information needed to reconstruct the potential. Thus we can, in principle, restrict our attention to free systems. First we review from [11, 15] how the structure equations for D nd order nondegenerate classical superintegrable systems are determined. Such a system admits a symmetry L = a ij p i p j + W if and only if {H, L} = 0, i.e., the Killing equations are satisfied and W i = λ j=1 aij V j. Here W i = x1 W with a similar convention for subscripts on V. The equations for W can be solved provided the Bertrand-Darboux equation x1 W = x W 1 holds. For a superintegrable system with independent symmetries L 1 = a ij p i p j + W (1), L = b ij p i p j + W (),we can solve the two independent Bertrand-Darboux equations for the potential to obtain the canonical system Here, V V 11 = A V 1 + B V, V 1 = A 1 V 1 + B 1 V. (37) A 1 = G + D () D, A = G 1 + D (3) D, B1 = G 1 D (0) D, B = G D (1) D, (38) ( a D = det 11 a, a 1 ) ( 3a 1 b 11 b, b 1, D (0) = det, ) a1 3b 1,, b1

Structure relations and Darboux contractions 11 ( 3a 1 D (1) = det, a11 a ) ( 3a 1 3b 1, b11 b, D () = det 1, ) ( a1 3a 1 3b 1 1,, D (3) = det 1, a11 a ) b1 3b 1 1, b11 b, where λ = exp G. If the integrability equations for (37) are satisfied identically then the solution space is 4-dimensional and we can always express the solution in the form V (x) = 3 j=1 a jv (j) (x) + a 4 where a 4 is a trivial additive constant. In this case the potential is nondegenerate and 3-parameter. Another possibility is that the solution space is -dimensional with general solution V (x) = a 1 V (1) (x) + a. For nondegenerate superintegrability, the integrability conditions for the canonical equations must be satisfied identically, so that V, V 1, V, V 11 can be prescribed arbitrarily at a fixed regular point. To obtain the integrability conditions for equations (37) we introduce the dependent variables W (1) = V 1, W () = V, W (3) = V 11, and matrices W (1) 0 0 1 A 1 B 1 0 w = W (), A (1) = A 1 B 1 0, A () = A B 1, (39) W (3) A 13 B 13 B 1 A A 3 B 3 A 1 A 13 = A 1 A 1 + B 1 A + A 1 A 1 B A 1 B 13 = B 1 B1 + A 1 B 1, A 3 = A 1 1 + B 1 A 1, B 3 = B1 1 + B 1 B 1. Then the integrability conditions for system xj w = A (j) w, j = 1,, must hold: A () 1 A (1) = A (1) A () A () A (1) [A (1), A () ]. (40) If and only if (40) holds, the system has a 4D vector space of solutions V. There is a similar analysis for a free nd order superintegrable system obtained by setting the parameter in a degenerate potential equal to zero, [14]: The free system retains all of the information needed to reconstruct the potential. All such degenerate superintegrable systems with potential are restrictions of nondegenerate systems obtained by restricting the parameters so that one nd order symmetry becomes a perfect square, e.g. L 1 = X. Then X is a 1st order constant, necessarily of the form X = ξ 1 p 1 + ξ p, without a function term. Since the degenerate systems are obtained by restriction, the potential function must satisfy the equations (37) inherited from the nondegenerate system, with the same functions A ij, B ij. In addition the relation {X, H} = 0 imposes the condition ξ 1 V 1 + ξ V = 0. By relabeling the coordinates, we can always assume ξ 0 and write the system of equations for the potential in the form V = C V 1, V = V 11 + C V 1, V 1 = C 1 V 1, where C (x 1, x ) = ξ 1 ξ, C (x 1, x ) = A ξ 1 ξ B, C 1 (x 1, x ) = A 1 ξ 1 ξ B 1. To find integrability conditions for these equations we introduce matrices ( ) ( ) ( V v =, B (1) 0 1 = 0 C + C C 1 C, B () 0 C = 0 C 1 V 1 ). (41) Then integrability conditions for system xj v = B (j) v, j = 1,, must hold: B () 1 B (1) = B (1) B () B () B (1) [B (1), B () ]. (4) If and only if (4) holds, the system has a D space of solutions V. Since V = constant is always a solution, (4) is necessary and sufficient for the existence of a nonzero 1-parameter potential system. In this case we can prescribe the values V, V at any regular point x 0 ; there will exist a unique V (x) taking these values.

1 R. Heinonen et al. 4.1 Free triplets Here we review information about free triplets that was presented and proved in [15]. A nd order classical free triplet is a D system without potential, H 0 = p 1 +p λ(x,y) and with a basis of 3 functionally independent second-order constants of the motion L (s) = i,j=1 aij (s) p ip j, a ij (s) = a ji (s), s = 1,, 3, L (3) = H 0. Since the duals of these constants of the motion are nd order Killing tensors, the spaces associated with free triplets can be characterized as D manifolds that admit 3 functionally independent nd order Killing tensors. As mentioned ( above, ) they have been classified in [10]. Since the vectors {h (s) }, h tr (s) (x, y, z) = a 11 (s), a1 (s), a (s) form a linearly independent set, there exist unique 3 3 matrices C (j) such that xj h (s) = C (j) h (s), j, s = 1,. By linearity, any element L = i,j=1 aij p i p j of the space of nd order symmetries spanned by the basis triplet is characterized by matrix equations xj h = C (j) h j = 1,, h tr (x, y, z) = ( a 11, a 1, a ). (43) In particular, at any regular point x 0 we can arbitrarily choose the value of the 3-vector h 0 and solve (43) to find the unique symmetry L of H 0 such that h(x 0 ) = h 0. A normalization condition for the C (j) : (43) is valid for a 11 = a = 1/λ, a 1 = 0, i.e., for H 0. From this and the requirement that the L are free constants of the motion we find C (1) = with the 4 functions C () the requirement G 1, G, 0 1 C() 11, 1 G 1 1 C() 1, 1 C() 11 G 1 C () 1, G C (), C() 1 11, C() 1, C() 1, C(), C () = C () 11, C() 1, G C () 11 C () 1, C(), C() 1 0, G 1, G, free. If we define the functions A1, B 1, A, B by C () 11 = 3 G 3 A1, C () 1 = 1 3 G 1 3 A, C () 1 = 1 3 G 1 1 3 B1, C () = 3 G 1 3 B, then equations (43) agree with the equations that are obtained from a superintegrable system with nondegenerate potential satisfying (37). Thus, for a free system there always exist unique functions A ij, B ij. Then necessary and sufficient conditions for extension to a system with nondegenerate potential V satisfying equations (37) are that conditions (40) hold identically. This analysis also extends, via restriction, to superintegrable systems with degenerate potential. A free triplet that corresponds to a degenerate superintegrable system is one that corresponds to a nondegenerate system but such that one of the free generators can be chosen as a perfect square. For these systems conditions (4) for the potential are satisfied identically. Similarly, we define a nd order quantum free triplet as a D quantum system without potential, H 0 = 1 λ(x) ( 11+ ), and with a basis of 3 algebraically independent second-order symmetry operators L k = 1 λ i,j=1 i (λa ij (k) j)(x), k = 1,, 3, a ij (k) = aji (k), L 3 = H 0 There is a 1-1 relationship between classical and quantum free triplets. 5 Free Darboux systems The Darboux spaces admit a 1-dimensional space of Killing vectors and a 4-dimensional space of nd order Killing tensors. Thus each space can admit at most one superintegrable system with

Structure relations and Darboux contractions 13 degenerate potential, and each space does so. We merely need to check that equations (40) are satisfied. Then we can compute the degenerate potential. Turning off the 1-parameter potential produces a single free degenerate quadratic algebra which we list below. There are no more possibilities. There are a number of possibilities for free triplets to define a nondegenerate quadratic algebra for a Darboux space, however. We classify the possibilities up to conjugacy under the action of the 1-parameter symmetry group of the manifold. Note that the 4-dimensional space of free constants of the motion is not obtained from the enveloping algebra of an underlying symmetry group. We shall see that there is a 1-1 relationship between free quadratic algebras and restrictions of quadratic algebras of nondegenerate superintegrable systems. We adopt the labeling of superintegrable systems and constants of the motion on Darboux spaces introduced in [18], using a tilde to differentiate between a free triplet and its associated superintegrable system. In the following sections, with the aid of MAPLE, we classify all possible free quadratic algebras generated by the nd order Killing tensors, up to conjugacy. Then, using MAPLE again, we verify for each quadratic algebra that the integrability conditions (4) are satisfied and we compute the nondegenerate potentials. Most of the results are presented in lists but in 5. we give more details on the construction of the superintegrable system with potential in one case. Each of the Darboux spaces can be embedded as a surface in 3 dimensions if we regard the ignorable variable as an angle i.e. X = f(x)cos(y), Y = f(x)sin(y), Z = h(x), and this is not unique, [19]. We give an illustrative example for each case. 5.0.1 Free Darboux 1 systems The space Darboux 1 (D 1 ) has free degenerate Hamiltonian D1D H = 1 4x (p x + p y) (44) with a single Killing vector K = p y and a basis, {H, K, X 1, X } for the 4-dimensional space of nd order Killing tensors. Here, X 1 = p x p y y x (p x + p y), The commutation relations are X = p y (yp x xp y ) y 4x (p x + p y). {K, X 1 } = H, {K, X } = X 1, {X 1, X } = K 3, and there is the functional relation 4HX + X1 + K4 = 0. The degenerate potential is V (x, y) = b 1x + b. As shown in [0], a possible embedding of this system in 3-dimensional Euclidean space with Cartesian coordinates X, Y, Z is X = x cos y, Y = x sin y, Z = ( 1 F (φ, ) + ) 4x 3 3 x, where x 1, y is π-periodic and sin φ = x + 1. Here F (φ, k) is an incomplete elliptic integral of the 1st kind. Then ds = 4x(dx + dy ) = dx + dy + dz. A general nd order symmetry, mod H, can be written as L 1 = a 1 X + a X 1 + a 3 K, and the translation group generated by K: x x, y y + α, leaves K and H invariant, but X 1 X 1 αh, X X + αx 1 α H. We classify the distinct free nondegenerate systems under this conjugacy action. We choose one generator L 1 and determine the possibilities for L such that L 1, L, H form a quadratic

14 R. Heinonen et al. algebra. Then we eliminate redundancies. Under conjugacy we can assume that L 1 takes one of the 3 forms X 1 + ak, X + ak, K. 1st case: We choose L 1 = X + ak and try to determine the possibilities for L, up to conjugacy under e αk such that L 1, L, H form a quadratic algebra. (As we go through the cases step by step, we ignore systems that have already been exhibited in earlier steps.) In general L = c 1 X 1 + c 3 K and c 1, c 3 are to be determined. We must require that R = b 1 L 3 1 + b L 3 + b 3 H 3 + b 4 L 1L + b 5 L 1 L + b 6 L 1 L H (45) + b 7 L 1H + b 8 L H + b 9 H L 1 + b 10 H L, for some constants b 1,, b 10. There are possible classes; 1. D1A: L 1 = X + bk, L = X 1 + ik, R = il 3 + 8iL 1L H + 4bL H + 16bH L 1. This class is superintegrable with A 1 = 0, A = x, B1 = 1 5x b + iy x(x b + iy), B = 3i x b + iy. G(x, y) = ln(4x), D = 1 (x b + iy). The potential of the superintegrable system is V (x, y) = b 1(x b + iy) x x b + iy + b x x b + iy + b 3 x + b 4. (This is missing in the tabulation in [0], but pointed out in [18] and [1].). D1B: L 1 = X, L = K, R = 4L 3 16L 1L H. A 1 = 0, A = x, B1 = 1 x, B = 3 y, G(x, y) = ln(4x), D = y. (Listed as a superintegrable system in the tabulation in [0].) The potential of the superintegrable system is V (x, y) = b 1(4x + y ) x + b x + b 3 xy + b 4. nd case: We choose L 1 = X 1. Then the only possibility is 3. D1C L 1 = X 1, L = K, R = 16L H. A 1 = 0, A = x, B1 = 1 x, B = 0, G(x, y) = ln(4x), D = 1. (Listed in the tabulation in [0].) The potential of the superintegrable system is V (x, y) = b 1(x + y ) x + b x + b 3y x + b 4.

Structure relations and Darboux contractions 15 5.0. Free Darboux systems The space Darboux (D ) has free degenerate Hamiltonian DD : H = x x + 1 (p x + p y) (46) with a single Killing vector K = p y and a basis, {H, K, X 1, X } for the 4-dimensional space of nd order Killing tensors. Here, X 1 = xp x p y + y x + 1 (p y x p x), The commutation relations are X = xyp x p y + (y x 4 )p y + x (1 y )p x x. + 1 {K, X 1 } = (K H), {K, X } = X 1, {X 1, X } = 4KX, and there is the functional relation 4HX + X 1 4K X 4H = 0. The degenerate potential is V (x, y) = b 1 x + 1 + b. As shown in [18] The line element ds can be realized as a two-dimensional surface embedded in three dimensions by X = y x + 1 x, Y Z = + 1, Y + Z = (x4 + 5x + 8y ) x + 1 3 x x 8x 8 arcsinhx, in which case, ds = dx + dy dz = x + 1 x (dx + dy ). A general nd order symmetry, mod H, can be written as L 1 = a 1 X + a X 1 + a 3 K and the translation group generated by K: x x, y y + α, leaves K and H invariant, but X 1 X 1 αh + αk, X X + αx 1 + α (K H). We classify the distinct free nondegenerate superintegrable systems under this conjugacy action. We choose one generator L 1 and determine the possibilities for L such that L 1, L, H form a quadratic algebra. Under conjugacy there are 3 possible choices: L 1 = X + ak, X 1, K. 1st case: We choose L 1 = X + ak and try to determine the possibilities for L, up to conjugacy under e αk such that L 1, L, H form a quadratic algebra. (As we go through the cases step by step, we ignore systems that have already been exhibited in earlier steps.) In general L = c 1 X 1 + c 3 K and c 1, c 3 are to be determined. We must require that R = b 1 L 3 1 + b L 3 + b 3 H 3 + b 4 L 1L + b 5 L 1 L + b 6 L 1 L H (47) + b 7 L 1H + b 8 L H + b 9 H L 1 + b 10 H L, for some constants b 1,, b 10. There are possible classes; 1. DC: L 1 = X, L = X 1, R = 4L 1 L + 16L 1 H 16L 1H. A 1 = 0, A = 3x 1 x(x + 1), B1 = (3x4 + x y ) x(x + 1)(x + y ), B = 6y x + y, G(x, y) = ln( x + 1 x ), D = x(x + y ).

16 R. Heinonen et al. The potential of the superintegrable system is ( ) x b V (x, y) = b 1 + x + y (x + 1) y + x + y + b 3 y x + y (Listed as a superintegrable system in [18].). DB: L 1 = X, L = K, R = 16L 1 L 16L 1L H + 16L H. A 1 = 0, A = 3x 1 x(x + 1), B1 = x(x + 1), B = 3 y, G(x, y) = ln( x + 1 x ), D = xy. The potential of the superintegrable system is ( V (x, y) = x x b 1 (x + y ) + b + 1 x + b ) 3 y + b 4. (Listed as a superintegrable system in [18].) nd case: We choose L 1 = X 1. Then the only possibility is 3. DA: L 1 = X 1, L = K, R = 16L 3 3L H + 16L H. A 1 = 0, A = 3x 1 x(x + 1), B1 = The potential of the superintegrable system is ( V (x, y) = x x b 1 (x + 4y ) + b ) + 1 x + b 3y + b 4. (Listed as a superintegrable system in [18].) 5.0.3 Free Darboux 3 systems The space Darboux 3 (D 3 ) has free degenerate Hamiltonian e x + b 4. x(x + 1), B = 0, G(x, y) = ln( x + 1 x ), D = x. D3E : H = 1 e x + 1 (p x + p y) (48) with a single Killing vector K = p y and a basis, {H, K, X 1, X } for the 4-dimensional space of nd order Killing tensors. Here, e x X 1 = 1 ex sin y p x p y + 1 4 e x + 1 cos y p x 1 4 X = 1 ex cos y p x p y + 1 4 e x + 1 sin y p x 1 4 The commutation relations are e x e x (e x + ) e x + 1 e x (e x + ) e x + 1 {K, X 1 } = X, {K, X } = X 1, {X 1, X } = 1 KH, cos y p y, sin y p y.

Structure relations and Darboux contractions 17 and there is the functional relation X 1 + X 1 4 H 1 K H = 0. The degenerate potential is V (x, y) = b 1 e x + 1 + b. As shown in [18], we can embed D3 as a surface in 3D Minkowski space with coordinates X, Y, Z in such a way as to preserve rotational symmetry. Let Then X = e x 1 + e x cos y, Y = e x 1 + e x sin y, Z = 6 1 ln ( 3 6 (6 + 5ex ) ) 3 + e x + 5e x + 1 3 6 (6 + 5ex ) 3 + e x + 5e x 1 dx + dy dz = (ex + 1) e x (dx + dy ). e x 3 + e x + 5e x. An alternate basis is {H, K, Y 1, Y } where ( Y 1 = e x p x p y + i e x e x + 1 p x i e x (e x ) ( + ) e x + 1 p y e iy, Y = e x p x p y i e x e x + 1 p x + i and X 1 = i 4 (Y 1 Y ), X = 1 4 (Y 1 + Y ). The new commutation relations are {K, Y 1 } = iy 1, {K, Y } = iy, {Y 1, Y } = 4iKH, e x (e x ) + ) e x + 1 p y e iy, and the functional relation is Y 1 Y H K H = 0. Returning to the first basis, we note that a general nd order symmetry, mod H, can be written as L 1 = a 1 X +a X 1 +a 3 K. and the translation group generated by K: x x, y y+α, leaves K and H invariant, but X 1 cos αx 1 sin αx, X sin αx 1 + cos αx. We classify the distinct free superintegrable systems under this conjugacy action. We choose one generator L 1 and determine the possibilities for L such that L 1, L, H form a quadratic algebra. Under conjugacy the choices are L 1 = X 1 + ak, X 1 + ix, X 1 + ix K, X 1, K. 1st case: We choose L 1 = X 1 + ak and try to determine the possibilities for a, L, up to conjugacy under e αk such that L 1, L, H form a quadratic algebra. (As we go through the cases step by step, we ignore systems that have already been exhibited in earlier steps.) There are 3 possible classes; 1st class: 1. D3A: L 1 = X 1 1 K, L = X + i K. R = L 3 1 + il3 1 8 H3 + il 1 L L 1 L + 1 (L 1 + L )H + 1 L 1H i 4 L H. Here, A 1 = 0, A = ex ( e x + e x iy e iy + e x ) (1 + e x )(e x + 1)(e x + e iy, B = 3ie iy ) e x + e iy, B 1 = 1 e x ( 1 + e x iy + e x + 4e x iy ) (1 + e x )(e x + e iy, ) G(x, y) = ln((e x + e x )), D = ex 8 (ex e iy ). The potential of the superintegrable system is V (x, y) = b 1 1 + e x + b e x 1 + e x+iy (1 + e x ) + b 3 e x+iy 1 + e x+iy (1 + e x ) + b 4. (This is a superintegrable system listed in [18].)

18 R. Heinonen et al.. D3B: L 1 = X 1, L = X, R = 1 8 H3 + 1 (L 1 + L )H. A 1 = 0, A = 1 e x 1 + e x, B = 0, B 1 = 1 (1 e x ) (1 + e x ), G(x, y) = ln((e x + e x )), D = ex 8. The potential of the superintegrable system is V (x, y) = ex (b e x 1 + e x (b cos y + 1 + b 3 sin y ) ) + b 4. (Listed as a superintegrable system in [18].) 3. D3C: L 1 = X 1, L = K, R = 4L 1 L + L H + L H. A 1 = 0, A = 1 e x 1 + e x, B1 = 1 + e x 1 + e x, B = 3 cot y, G(x, y) = ln((e x + e x )), D = ex 4 sin y. The potential of the superintegrable system is V (x, y) = ( ex e x + 1 b 1 + e x ( b cos y (Listed as a superintegrable system in [18].) + b 3 sin y ) ) + b 4. nd class: We choose L 1 = X 1 + ix. Then the only new possibility is 4. D3D: L 1 = X 1 + ix, L = K, R = 4L 1 L. A 1 = 0, A = 1 e x 1 + e x, B1 = 1 + e x 1 + e x, B = 3i G(x, y) = ln((e x + e x )), D = i 4 ex+iy. The potential of the superintegrable system is V (x, y) = ex 1 + e x (b 1e iy + b e iy ) + b 3 1 + e x + b 4. (Listed as a superintegrable system in [18].) 3rd class: We choose L 1 = X 1 + ix K. There are no new possibilities.

Structure relations and Darboux contractions 19 5.0.4 Free Darboux 4 systems The spaces Darboux 4 (D 4 (b)) have free degenerate Hamiltonian D4(b)D : H = sin x cos x + b (p x + p y), (49) b ±, with a single Killing vector K = p y and a basis, {H, K, X 1, X } for the 4-dimensional space of nd order Killing tensors. Here, X 1 = e y ( H + cos x p y + sin x p x p y ), X = e y ( H + cos x p y sin x p x p y ). The commutation relations are {K, X 1 } = X 1, {K, X } = X, {X 1, X } = 8K 3 4bKH, and there is the functional relation X 1 X K 4 bk H H = 0. The degenerate potential is V (x, y) = b 1 cos x + b + b. A general nd order symmetry, mod H, can be written as L 1 = a 1 X + a X 1 + a 3 K, and the translation group generated by K: x x, y y + α, leaves K and H invariant, but X 1 e α X 1, X e α X. Also the reflection y y, x x leaves H and K fixed but X 1 X. We classify the distinct free superintegrable systems under this conjugacy action. We choose one generator L 1 and determine the possibilities for L such that L 1, L, H form a quadratic algebra. Under conjugacy the choices are L 1 = K + ax, X 1 + X + ak (a = 0,, or a 0, ), X, K. 1st case: We choose L 1 = K + ax and try to determine the possibilities for L, up to conjugacy under e αk such that L 1, L, H form a quadratic algebra. (As we go through the cases step by step, we ignore systems that have already been exhibited in earlier steps.) We first try L = X + cx 1. The class ac 0 doesn t yield a quadratic algebra. There are 3 other possible classes; 1. D4(b)A: L 1 = K, L = X + X 1, R = 64L 3 1 + 16L 1L 64bL 1 H 64L 1H. A 1 = 0, A = ( sin x b cos x ), B 1 = 4( sin x + b cos x + ) sin x ( cos x + b) sin x ( cos x + b) B = 6( 1 + e4y cos x ), G(x, y) = ln( b 1 e4y sin ), D = x The potential of the superintegrable system is ( ) V (x, y) = sin x b1 cos x + b sinh y + b sinh + y (Listed as a superintegrable system in [18].). D4(b)B: L 1 = K, L = X, R = 16L 1 L. A 1 = 0, A = ( sin x b cos x ) sin x ( cos x + b) b cos x sin x G(x, y) = ln( sin ), D = e y. x The potential of the superintegrable system is ( ) V (x, y) = sin x b1 cos x + b sin x + b e 4y + b 3 e y + b 4. (Listed as a superintegrable system in [18].) sin x (e y + e y ). b 3 cos x + b + b 4., B 1 = 4( sin x + b cos x + ), B = 6, sin x ( cos x + b)

0 R. Heinonen et al. nd case: We choose L 1 = K + ax, L = X 1. Then the only possibility is a = 0, which is redundant. 3rd case: We choose L 1 = X 1 + X + ak, L = X + ck. We generate a quadratic algebra for the system 3. D4(b)C: L 1 = X 1 + X + K, L = X + K. R = 16bH 3 + 16L 1 L 16L 1 L + 16bL 1 L H 16bL H 16L 1H. A 1 = 0, A = 4 sin x + b cos x + 4 sin x ( cos x + b), B 6(e 4y 1) = e y cos x + e 4y + 1, B 1 = bey cos x + 6be y + 4b cos x + 4be 4y cos x 8 cos x + 16 cos 4 x sin x (e 4y + e y, cos x + 1)( cos x + b) + 8e4y cos x + 16e 4y cos 4 x 4e y cos 3 x + 16e y cos x sin x (e 4y + e y, cos x + 1)( cos x + b) b cos x sin x G(x, y) = ln( sin ), D = ( cos x + e y + e y ). x The potential of the superintegrable system is e y ( b 1 V (x, y) = Z + (1 e y ) Z + b Z + (1 + e y ) Z + b 3 e y cos x b+ sin x + b cos x where Z = (1 e y ) + 4e y cos x. (Listed as a superintegrable system in [18].) ) + b 4, 4th case: We choose L 1 = X 1, L = X. We do not generate a quadratic algebra. 5.1 Alternate Free Darboux 4 systems There is an alternate form of D4(b) that we shall employ. Set x = ix, y = iy, J = ik = p Y, Y 1 = X 1 + X Then for b ± we have H =, Y = X 1 X. i sinh X cosh X + b (p X +p Y ), Y 1 = cos(y )(H cosh X p Y ) sin(y ) sinh X p X p Y, (50) Y = sin(y )(H cosh X p Y ) + cos(y ) sinh X p X p Y. The structure equations are {J, Y 1 } = Y, {J, Y } = Y 1, {Y 1, Y } = 4J 3 + bj H, Y 1 + Y J 4 bj H H = 0. We can embed D4(b) as a surface in 3D Minkowski space with Cartesian coordinates X, Y, Z. Let X = y b + cosh x, Y Z = b + cosh x, Y + Z = y b + cosh x sinh x + cosh 4x + + 4b cosh x + b b + cosh x (cosh 4x + 3 + b cosh x) dx.

Structure relations and Darboux contractions 1 Then ds = cosh x + b sinh x (dx + dy ) = dx + dy dz. The change of variable u = e x converts the integral into an elliptic integral in u that can be evaluated as a rather complicated sum of elementary functions and the elliptic integrals of types one, two and three. In terms of this alternate form the superintegrable systems can be expressed as: 1. D4(b)A : L 1 = J, L = Y 1, R = 16L 3 1 16L 1L + 16bL 1 H + 16L 1H. A 1 = 0, A = ( sinh X + b cosh X + ) sinh X ( cosh X + b), B cos Y = 6( sin Y ), B 1 = 4(sinh X + b cosh X + ) b + cosh X, G(x, y) = ln( sinh X ( cosh X + b) sinh X ).. D4(b)B : L 1 = J, L = Y 1 iy, R = 16L 1 L. A 1 = 0, A = ( sinh X + b cosh X + ), B = 6i, sinh X ( cosh X + b) B 1 = 4(+ sinh X + b cosh X + ) b + cosh X, G(x, y) = ln( sinh X ( cosh X + b) sinh X ). 3. D4(b)C : L 1 = Y 1 J, L = Y. R = 4bH 3 8L 3 1 8L 1L + 4b(L 1 + L )H + 8L 1H. A 1 = 0, A = 4 sinh X + b cosh X + 4, B = sinh X ( cosh X + b) 6 sin Y cosh X + cos Y, B 1 = b cosh X 6b + (8b + 16) cosh X cos Y sinh X ( cos Y cosh X)( cosh X + b) 3 cosh4 X cos Y + 4 cosh 3 X 16 cosh X + cosh X, G(x, y) = ln(b sinh X ( cos Y cosh X)( cosh X + b) sinh X ). 4. D4(b)D : (free degenerate) J, L 1 = Y 1, L = Y, Y 1 + Y J 4 bj H H = 0. 5. Summary and an example From the results of [15] and the calculations of the preceding sections we see that Theorems 1 and are valid for Darboux spaces. We use DC as an example to give more details about how a nondegenerate superintegrable system with potential is induced from a free system : From Section 5.0. we have the free Darboux system L 1 = X, L = X 1, R = 4L 1 L + 16L 1H 16L 1 H. The potential equations are determined by A 1 = 0, A = 3x 1 x(x + 1), B1 = (3x4 + x y ) x(x + 1)(x + y ), B = 6y x + y, (51)

R. Heinonen et al. and the general solution is x V (x, y) = x + y (x + 1) ( b 1 + The induced classical system has a basis of symmetries H = b y + x + y + b 3 y + x + y ) + b 4. (5) x x + 1 (p x + p y) + V (x, y), L 1 = xyp y p x + (y x 4 )p y + x (1 y )p x x + W (1), + 1 where L = xp x p y + y(p y x p x x + 1 + W (), (53) W (1) = 1 b 1 y(1 x ) + b (( y + x + y ) + 1) b 3 ((y + x + y ) + 1) (x + 1), x + y W () = b 1 x + y b (( y + x + y ) 1) b 3 ((y + x + y ) 1) 4x. + 4 The Casimir is R 4L 1 L 16L 1H + 16L 1 H 4(b + b 3 )H + 16b 4 L 1 + (b + b 3 )L 3b 4 HL 1 + (8b b 4 cb1 + 4cbb 3 + 8b 3 b 4 )H + (16b 4 + b 1)L 1 + b 1 ( b 3 + b )L + ( 1 4 b 1b + 1 4 b 1b 3 + b 1b 4 4b 3 b 4 4b b 3 b 4 4b b 4). (54) We will not work out the details of the quantum case but merely note that the potential terms of the symmetries remain unchanged as do the nd order kinetic energy terms, but there are now 1st order terms: the Hamiltonian kinetic energy is replaced by the Laplace-Beltrami operator on D and the other generating symmetry operators are formally self-adjoint with respect to the D volume measure x dx dy/(x + 1). 6 Contractions of Darboux systems Recall that the scalar curvature of a space with metric ds = e G(x,y) (dx + dy ) where λ(x, y) = e G(x,y) is C = e G ( xx G + yy G). Constant curvature spaces are just those for which C is constant; flat spaces are those for which C = 0. Theorem 4. A Darboux or Koenigs superintegrable system cannot be obtained as a geometric contraction of a nonzero constant curvature or flat space superintegrable system. Proof: Suppose we have a contraction of a nonzero constant curvature system with Hamiltonian H = p x+p y. Then there is a 1-parameter family of Hamiltonians H = p x +p y e G e G (ɛ) where G (1) = G, G depends smoothly on ɛ in the interval 0 ɛ 1, and G (0) defines the metric of the target manifold for the contraction. Further, for ɛ 0 the metric defined by G will be a scalar multiple of a metric on the original constant curvature system. Thus we have C(ɛ) = e G ( x x G + y y G ) for ɛ 0 where C(ɛ) is nonzero and independent of x, y. In the limit we obtain the constant C(0), so the target manifold must be flat or of nonzero constant curvature, hence not a Darboux or Koenigs manifold. Similarly, if the original manifold is flat the target manifold must also be flat..

Structure relations and Darboux contractions 3 Now we investigate contractions involving Darboux superintegrable systems. From Theorem 4 such systems cannot be obtained as contractions of constant curvature systems. (However, they are all Stäckel equivalent to constant curvature space systems.) Thus we limit ourselves to the search for contractions such that the originating manifold is a Darboux space. In distinction to the case of constant curvature spaces as originating manifolds (where all quadratic algebra contractions were induced by Lie algebra contractions of e(, C) and o(3, C)) here the Darboux spaces have only a one-dimensional Lie symmetry algebra, so we have to search for new contraction mechanisms. However, we can restrict our search to free systems and then induce the superintegrable system contractions automatically. To induce the contractions through Lie algebra methods we will employ the conformal symmetry algebra for each free Darboux system generated by functions Q(x, y, p x, p y ) = A(x, y)p x + B(x, y)p y that satisfy the relations {H, Q} = R Q (x, y)h, {K, Q} = K Q (x, y)k, K = p y. A straightforward calculation gives the same algebra in each case: G3 : P x = p x, P y = p y, M = x p x +y p y, {P x, P y } = 0, {P x, M} = P x, {P y, M} = P y. (55) The generalized Inönü-Wigner contractions and their geometric implementations are: 1. P x, P y, ɛm; x = x + 1 ɛ, y = y, or x = x, y = y + 1 ɛ.. ɛp x, ɛp y, M; x = ɛx, y = ɛy. Up to an isomorphism of the Lie algebra that is obtained by contraction of G3 there are no continuous one parametric contractions besides these two, []. Each of these contractions does induce a corresponding contraction of each free Darboux system and we have found no contractions of Darboux systems other than these. This approach is compatible with the use of generalized Inönü-Wigner contractions for the symmetry algebras o(3, C) and e(, C) of constant curvature spaces. Indeed if we compute the conformal symmetry algebra for each free constant curvature system generated by functions Q(x, y, p x, p y ) = A(x, y)p x + B(x, y)p y that satisfy the relations {H, Q} = R Q (x, y)h, {L j, Q} = 3 l=1 K (l) Q,j (x, y)l l, where the L j form a basis for either o(3, C) or e(, C) we find o(3, C) in the case of the sphere and the affine algebra, the semidirect product of e(, C) and the dilation M, for flat space. One can show that the geometrical implementations of generalized Inönü-Wigner contractions of the affine algebra either agree with those of e(, C) or the contractions cannot be implemented geometrically. See for example [3] Table XXIX for the ordinary Inönü-Wigner contractions of the affine algebra. 6.1 D1 contractions We list approaches to finding the contractions (not the contractions themselves). 1. Let x = x + 1 ɛ, y = y, H = 4 ɛ H. Then as ɛ 0 we have H = p x + p y, X 1 p x p y, X p y (y p x x p y ) p y ɛ. Contractions constructed from such limits would have flat space as the target manifold.

4 R. Heinonen et al.. Let x = x, y = y + 1 ɛ, H = H. Then as ɛ 0 we have H = p x + p y 4x, X 1 p x p y y x (p x + p y ) ɛ H, X p y (y p x x p y ) + p x p y ɛ y 4x (p x + p y y ) ( ɛx + 1 4ɛ x )(p x + p y ). Contractions constructed from such limits would have D1 itself as the target manifold. 3. Let x = ɛx, y = ɛy, H = 4ɛH. Then as ɛ 0 we have H = p x +p y x and ɛ X 1 = p x p y y x (p x + p y ), ɛx = p y (y p x x p y ) y 4x (p x + p y ), ɛk = p y. Contractions constructed from such limits would have D1 again as the target manifold. 6. D contractions We list approaches to finding the contractions (not the contractions themselves). 1. Let x = x + 1 ɛ, y = y, H = H. Then as ɛ 0 we have H = p x + p y, ɛx 1 p x p y, ɛ X p y. Contractions constructed from such limits would have flat space as the target manifold.. Let x = ɛx, y = ɛy, H = H. Then as ɛ 0 we have H = x (p x + p y ) and ɛx 1 x p x p y + y p y, ɛx = x y p x p y + y p y + x p x, ɛk = p y. Contractions constructed from such limits would have the complex -sphere as the target manifold. 3. Let x = x, y = y + 1 ɛ, H = H. Then as ɛ 0 we have H = x x +1 (p x + p y ), K = K, X 1(x, y, p x, p y ) = X 1 + ɛ (K H), X (x, y, p x, p y ) = X + 1 ɛ (K H) + 1 ɛ X 1. Contractions constructed from such limits would have D as the target manifold. 6.3 D3 contractions We list approaches to finding the contractions (not the contractions themselves). 1. Let x = x + ln( 1 ɛ ), y = y, H = ɛh. Then as ɛ 0 we have K = p y H = e x (p x + p y ), ɛx 1 ex ɛx ex ( cos y p x p y + 1 sin y (p x py ) ) (sin y p x p y + 1 cos y (p x py ) ), Contractions constructed from such limits would have flat space as the target manifold. In terms of flat space Cartesian coordinates X = r cos θ, Y = r sin θ we have e x = 4 r, y = θ, H = p X+p Y, ɛx 1 1 4 (p X p Y ), ɛx 1 p Xp Y, K = 1 (Xp Y Y p X ) = 1 J.

Structure relations and Darboux contractions 5. Let x = ɛx, y = ɛy, H = 4ɛ H. Then as ɛ 0 we have H = p x + p y and ɛ X 1 1 8 (p x 3p y ), ɛ X 1 p x p y, ɛk = p y. Contractions constructed from such limits would have flat space as the target manifold. 3. Let x = x, y = y + i ln ɛ, H = H. Then as ɛ 0 we have Y 1(x, y, p x, p y ) ɛy 1, Y (x, y, p x, p y ) 1 ɛ Y, K = p y Contractions constructed from such limits would have D3 as the target manifold. 6.4 D4 contractions We list approaches to finding the contractions (not the contractions themselves). 1. Let X = X + ln ɛ, Y = Y, H = 4ɛ H, with b fixed. Then as ɛ 0 we have J = p Y and H = e X (p X + p Y ), J = p Y, 4ɛ Y 1 cos(y )(H e X p Y ) + sin(y ) e X p X p Y p y p x, 4ɛ Y sin(y )(H e X p Y ) cos(y )e X p X p Y p x p y, where x, y are standard Cartesian coordinates. Contractions constructed from such limits would have flat space as the target manifold.. Let X = ɛx, Y = ɛy, H = H, with b fixed. Then as ɛ 0 we have H = 4X +b (p X +p Y ) and Y 1 p Y + H (X Y )p Y 4X Y p X p Y, ɛy Y p Y + X p X p Y, ɛj = p Y. Contractions constructed from such limits would have the complex -sphere as the target manifold. 3. For this case it is most convenient to take the Hamiltonian in the form (49). Let x = x, y = y + 1 ln ɛ, H = H, with b fixed. Then we have X 1(x, y, p x, p y ) = 1 ɛ X 1, X (x, y, p x, p y ) = ɛx. Contractions constructed from such limits would have D4(b) as the target manifold. The manifold D4(b) depends on a parameter, so we can extend the contractions that we consider by allowing the parameter to vary. We find the following new contractions: 4. Let X = X 1 ln(ɛ), Y = Y, b = 1 ɛ, H = ɛ H. Then as ɛ 0 we have H = 1 e X e X +1 (p X + p Y ), J = 1 J = p Y, and ɛ Y 1 cos(y )(H e X p Y ) sin(y )e X p X p Y, ɛ Y sin(y )(H e X p Y )+cos(y )e X p X p Y. Contractions constructed from such limits would have D3 as the target manifold.

6 R. Heinonen et al. 5. Let X = ɛx, Y = ɛy, b = + ɛ, H = ɛ 4 H. Then as ɛ 0 we have H X X +1 (p X + p Y ), J = ɛ J = p Y, and ( ) Y 1 + H J ɛ H X Y p X p Y + (Y X 4 )p Y + X (1 Y )p X, X + 1 Y (X p X p Y + Y ) ɛ X + 1 (p Y X p X ). Contractions constructed from such limits would have D as the target manifold. 6. Let X = ɛx + ln(ɛ), Y = ɛy, b = 1 ɛ, H = ɛ 4 H. Then as ɛ 0 we have H 1 4X (p X + p Y ), J = ɛj = p Y, and ( ) Y 1 + H 1 ɛ J 1 ɛ p Y (Y p X X p Y ) Y 4X (p X + p Y ), Y 1 (p ɛ 3 X p Y Y ) X (p X + p Y ). Contractions constructed from such limits would have D1 as the target manifold. 7. Using Hamiltonian (50) we let b = + ɛ, X = X, Y = Y, H H. Then H = cosh X (p X + p Y ), J = J = p Y, Y 1 = cos(y )(H cosh X p Y ) sin(y ) sinh X p X p Y, Y = sin(y )(H cosh X p Y ) + cos(y ) sinh X p X p Y. Contractions constructed from such limits would have the complex -sphere as the target manifold. Expressed in terms of the symmetries of the -sphere we have Y 1 = J 1 J, Y = J 1 J, J = J 3 6.5 Summary and examples of Darboux contractions From the results of [15] and the preceding sections we have an analog of Theorem 3 for Darboux spaces: Theorem 5. Every Lie algebra contraction of G3 induces a contraction of a free geometric quadratic algebra Q on a Darboux space, which in turn induces a contraction of the quadratic algebra Q with potential. This is true for both classical and quantum Darboux algebras. Some examples follow. 1. We describe how a Lie algebra contraction induces the contraction of DC to S4, including the potential terms. Recall for DC we have the symmetries (53), potential (5), canonical equations (51) and Casimir (54). The coordinate implementation is defined by x = ɛy, y = ɛx. In terms of the coordinates x, y the canonical equations become A 1 = 0, A = 1 y, B 1 = x y y (x + y ), b = 6x (x + y ). The contraction is defined by L 1 = L 1, L = ɛ L 1, H = H, R = ɛ R.

Structure relations and Darboux contractions 7 In the limit we find H = y (p x + p y ) + V = S S 1 S 3 + V, L 1 = S + W 1, L = S 1 S + W, where V = y x + y ( c 1 + c x + y + x + c 3 x + y x ) + c 4 and b 1 = c 1 ɛ, b = c, b 3 = c 3, b 4 = c 4. Here S 1 = p x, S = x p x + y p y, S 3 = 1 (x y )p x + x y p y, are the 1st order generators of the symmetry algebra for the free Hamiltonian H 0 = y (p x + p y ) on the Poincaré upper half plane: x real, y > 0. The Casimir becomes in the limit : R 4L1 L + (c + c 3 )L c 1 4 H + c 1 4 L 1 + c 1 (c c 3 )L + c 1 16 (c + c 3 + 4c 4 ) = 0.. D4(b)A D3E: We again give more details in our nd example, which involves changing the parameter b. For the degenerate system H = we can get sinh (x) cosh(x) + b (p x + p α y) + cosh(x) + b, (56) H = 1 e x e x + 1 (p x + p y ) + β e x + 1 as a contraction case by taking x = x 1 y ln(ɛ), y =, b = 1 ɛ, H = 8ɛH, α = β 8ɛ. (58) Then we have J = J = p y, and Y 1 4ɛY 1, Y 4ɛY. 3. D1A Ẽ9(a = 1): For system D1A we have L 1 = X + bk, L = X 1 + bk, R = il 3 + 8iL 1 L H + 4bL H + 16bH L 1. We choose contraction 1) for D1A: x = x + 1 ɛ, y = y, H = ɛ 4 H and let b = 1 ɛ. Then as ɛ 0 we find L 1 L 1 = p y (y p x x p y ), L L = p y (p x + ip y ), H = p x + p y, and with R = {L 1, L } the Casimir becomes R = il 3 + L H. 4. D1B Ẽ: For system D1B we have L 1 = X, L = K, R = 4L 3 16L 1 L H. We choose contraction 1) for D1B: x = x + 1 ɛ, y = y, H = ɛ 4 H. Then as ɛ 0 we find L 1 L 1 1 ɛ L, L 1 = p y (x p y y p x ), L = L = p y, H = p x + p y, and with R = {L 1, L } the Casimir becomes R = 4L 3 + 4L H. (57)

8 R. Heinonen et al. 5. D1C Ẽ3 : For system D1C we have L 1 = X 1, L = K, R = 16L H. We choose contraction 1) for D1C: x = x + 1 ɛ, y = y, H = ɛ 4 H. Then as ɛ 0 we find L 1 L 1 = p x p y, L = L = p y, H = p x + p y, and with R = {L 1, L } the Casimir becomes R = 0. 6. D4(b)C S7 as b : For system D4(b)C we have L1 = Y 1 J, L 1 = Y, R = 4bH 3 8L 1 3 8L 1 L + 4b(L 1 + L )H + 8L 1 H. In the limit as b we have L 1 = J1 H, L = J 1 J where J 1, J, J 3 are the symmetries of the sphere and H = J1 + J + J 3. The limit system is S7. 7 Tables/figures describing the contractions of the nondegenerate and degenerate Darboux systems We now tabulate the contractions of Darboux systems. The system on the left is the starting system, and an asterisk indicates that a system is degenerate 7.1 D1 contraction table 1: D1 to flat space : D1 to D1 3: D1 to D1 D1A Ẽ3, Ẽ9 D1C, D1A D1C, D1A D1B Ẽ D1C D1B D1C Ẽ3 D1C D1C D1D Ẽ5 D1D D1D 7. D contraction table 1: D to flat space : D to -sphere 3: D to D DA Ẽ3 S1 DA DB Ẽ S DA DC Ẽ3 S4 DA DD Ẽ5 S5 DD 7.3 D3 contraction table 1: D3 to flat space : D3 to flat space 3: D3 to D3 D3A Ẽ3 Ẽ3 D3D D3B Ẽ3 Ẽ3 D3B D3C Ẽ1 Ẽ D3D D3D Ẽ8 Ẽ3 D3D D3E Ẽ3 Ẽ5 D3E

Structure relations and Darboux contractions 9 7.4 D4 contraction table 1: D4 to flat space : D4 to -sphere 3: D4 to D4 D4(b)A Ẽ1 S D4(b)B D4(b)B Ẽ8 S1 D4(b)B D4(b)C Ẽ3 S4 D4(b)B D4(b)D Ẽ3 S5 D4(b)D 4: D4 to D3 5: D4 to D 6: D4 to D1 7: D4 to -sphere D4(b)A D3C DB D1B S9 D4(b)B D3D DA D1C S D4(b)C D3B DC D1C S7 D4(b)D D3E DD D1D S3 7.5 Free nondegenerate Darboux systems to degenerate systems The following are not contractions in the usual sense, but we show which nondegenerate systems restrict to degenerate ones. nondegenerate D1B D1C DA DB D3C D3D D4(b)A D4(b)B degenerate D1D D1D DD DD D3E D3E D4(b)D D4(b)D 7.6 Figures What follows are diagrams illustrating the contractions of Darboux systems. Boxes represent systems (blue boxes being nondegenerate and red boxes being degenerate), and arrows represent contractions. In these diagrams, certain contractions that do not affect the overall hierarchy have been omitted for aesthetic reasons, and a dotted arrow indicates that the limiting process inducing the contraction changes the free Hamiltonian. A diagram illustrating the nondegenerate Darboux systems which restrict to degenerate ones is also presented.

30 R. Heinonen et al. Figure 1. Diagram indicating the contractions of nondegenerate Darboux systems. Figure. Diagram indicating the contractions of degenerate Darboux systems. Figure 3. Diagram indicating restrictions of nondegenerate Darboux systems to degenerate ones. For completeness we give the definitions of the constant curvature superintegrable systems that are targets of Darboux contractions. The following lists contain the defining relations for the free systems and the potentials of the superintegrable systems. Degenerate Euclidean targets: Here the coordinates are x, y and p 1 = p x, p = p y, and J = xp yp 1. 1. Ẽ3: H = p 1 + p, X = J, L 1 = p 1, L = p 1 p, Casimir : L L 1(L 1 H) = 0, potential : V = α(x + y ).