arxiv:hep-th/ v1 23 Aug 1993

Transcription

1 hep-th/930809, IMA Preprint Series no.53, July 993 QUADRICS ON COMPLEX RIEMANNIAN SPACES OF CONSTANT CURVATURE, SEPARATION OF VARIABLES AND THE GAUDIN MAGNET arxiv:hep-th/930809v 23 Aug 993 E.G. KALNINS, V.B. KUZNETSOV*and WILLARD MILLER, Jr. Department of Mathematics and Statistics University of Waikato Hamilton, New Zealand Department of Mathematics and Computer Science University of Amsterdam Plantage Muidergracht TV Amsterdam The Netherlands. School of Mathematics and Institute for Mathematics and its Applications, University of Minnesota, Minneapolis, Minnesota 55455, USA. Abstract. We consider integrable systems that are connected with orthogonal separation of variables in complex Riemannian spaces of constant curvature. An isomorphism with the hyperbolic Gaudin magnet, previously pointed out by one of us, extends to coordinates of this type. The complete classification of these separable coordinate systems is provided by means of the corresponding L-matrices for the Gaudin magnet. The limiting procedures (or ǫ calculus) which relate various degenerate orthogonal coordinate systems play a crucial result in the classification of all such systems. PACS: b, m, p, Fd. Classical Integrable Systems on Complex Constant Curvature Spaces and the Complex Gaudin Magnet. Separation of variables in the Hamilton Jacobi equation n (.) H(p,..., p n ; x,..., x n ) = g αβ p α p β = E, α,β= p α = W x α, α =,..., n *This author was supported by the National Dutch Science Organization (NWO) under the Project # Work supported in part by the National Science Foundation under grant DMS

2 2 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. amounts to looking for a solution of the form n (.2) W = W α (x α ; h,..., h n ), h n = E The solution is said to be a complete integral if det( 2 W/ x i h j ) n n 0. The solution then describes that of free motion on the corresponding Riemannian space with contravariant metric g αβ. Indeed, if we require b j = W h j h n tδ nj, j =,..., n for parameters b j, we find that the functions x(b,h), p(b,h) satisfy Hamilton s equations (.3) ẋ α = H p α, ṗ α = H x α In this article we allow the Riemannian space to be complex and we consider variable separation of (.) on the following two classes of spaces. () The n dimensional complex sphere S nc. This is commonly realised by the set of complex vectors x = (x,..., x ) which satisfy x2 α = and have infinitesmal distance dx dx = dx2 α (2) The n dimensional complex Euclidean space E nc. This is the set of complex vectors x = (x,..., x n ) with infinitesmal distance dx dx = n dx2 α. A fundamental problem from the point of view of separation of variables on these manifolds is to find all inequivalent coordinate systems. As yet, this is an unsolved problem, principally because many such coordinate systems are intrinsically nonorthogonal. For orthogonal coordinate systems the problem is completely solved and in this case the constants h i occuring in the complete integral can be chosen to be the values of an involutive set of constants of motion (.4) A j = where n α,β= {F(x α, p β ), G(x α, p β )} = a (j) αβ p αp β, j =,..., n, A n = H, {A j, A k } = 0, n ( F G F G ), p γ x γ x γ p γ γ= α, β =,..., n is the Poisson bracket. These constants of the motion are such that () each of the tensors a (j) αβ is a Killing tensor and satisfies Killing s equations (α a (j) βγ) = 0, [], and (2) the A j can be represented as a sum of quadratic elements of the enveloping algebra of the Lie algebra of symmetries of each of these two considered spaces. The Lie algebras of these spaces have respectively bases of the form () SO(n + ) : M αβ = x α p β x β p α, α, β =,..., n +. Here, {M αβ, M γδ } = δ αγ M δβ + δ αδ M βγ + δ βγ M αδ + δ βδ M γα (2) E(n) : M αβ, P γ = p γ, α, β, γ =,, n, α β. Here, {M αβ, P γ } = δ βγ P α δ αγ P β, {P α, P β } = 0.

3 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 3 In this section the separable coordinate systems classified in [,5,0] are given an algebraic interpretation. This is done using the complex analogue of the isomorphism between all integrable systems connected with all possible separable systems and the m site SO(2, ) Gaudin magnet, [2,4]. The m-site complex Gaudin magnet can be realised as follows. Consider the direct sum of Lie algebras, each of rank, (.5) A = m so α (3,C). The generators s α C 3, α =,, m of A satisfy the Poisson bracket relations (.6) {s i α, sj β } = δ αβǫ ijk s k α. The following metric will be used subsequently when norms and scalar products are calculated: s 2 α = (s α,s α ) = (s α) 2 + (s 2 α) 2 + (s 3 α) 2, (s α,s β ) = s αs β + s 2 αs 2 β + s 3 αs 3 β. If for each α, s 2 α = c2 α then the variables s α lie on the direct product of n complex spheres in C 3. The complex Gaudin magnet is the integrable Hamiltonian system described by the n integrals of motion H α which are in involution with respect to the Poisson bracket, (.7) H α = 2 m (sα,s β ), {H α, H γ } = 0. e α e β β= (We will give a simple proof of this involution property later.) Here the e α are taken to be pairwise distinct. This integrable Hamiltonian system is called an m- site SO(3,C) XXX Gaudin magnet. The H α are all quadratic functions in the generators of the A algebra and the following identities are satisfied: (.8) m H α = 0, m m e α H α = J 2 c 2 α, where we have introduced the variables m (.9) J = s α, J 2 = (J,J), the total sum of the momenta s α. Indeed (.0) {J i, J j } = ǫ ijk J k, {J k, H α } = 0, i, j, k =, 2, 3, α =,, m. The complete set of involutive integrals of motion is provided by H α,j 2 and, for example, J 3. The integrals are generated by the 2 2 L matrix [3,4,6,7,8], (.) L(u) = m ( s 3 α 2 3 s α ) is2 α = ( A B ),

4 4 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. where (.2) det L(u) = (A 2 + BC) = m H α m Furthermore, L(u) satisfies the linear r matrix algebra [2,2] c 2 α ( ) 2. (.3) { L(u), 2 L(v)} = i u v [P, L(u) + 2 L(v)], i = where (.4) P =, L(u) = L(u) I, 2 L(v) = I L(v). The algebra specified by (.3) contains all the necessary Hamiltonian structure of the problem in question. Note that (.3) is equivalent to the easily proved relations (.5) {A(u), A(v)} = {B(u), B(v)} = {C(u), C(v)} = 0, {A(u), B(v)} = {B(u), C(v)} = i i (B(v) B(u)), {A(u), C(v)} = (C(u) C(v)), u v u v 2i (A(v) A(u)). u v From (.5) and the Leibnitz property of the Poisson bracket it is straightforward to deduce that {det L(u), detl(v)} = 0. In his article [4] Kuznetsov has explicitly given the nature of the isomorphism between the XXX Gaudin magnet models and the separation of variables on the n dimensional real sphere S n. The purpose of this article is to extend these ideas to complex orthogonal coordinate systems on the complex n sphere S nc and, of course, as a consequence complex Euclidean space E nc. Following Kuznetsov [4,7] in the case of the sphere, we set c α = 0, α =,, n +. The coordinates on the resulting cones are parametrised by (.6) s α = 4 (p2 α + x2 α ), s2 α = i 4 (p2 α x2 α ), s3 α = i 2 p αx α. It follows from (.6) that {x α, x β } = {p α, p β } = 0, {p α, x β } = δ αβ. Introducing the new variables M αβ = x α p β x β p α which are the generators of rotations we have (s α,s β ) = 8 M2 αβ. These generators satisfy the commutations relations given previously. This equality establishes the simple quadratic connection between the generators s α of A and the

5 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 5 subsequent discussion) transform into the following integrals for the free motion on the n sphere H = h α H α = h α h β Mαβ 2 4 e α e. β For h γ = e γ we obtain the Casimir element of the so(n + ) algebra α<β M2 αβ. The total momentum J takes the form α<β J = 4 (p p + x x), J2 = i 4 (p p x x), J3 = i 2 p x where the scalar product for the vectors x and p in C is Euclidean.The quantities M and J form the direct sum so(n + ) so(3) as a result of the commutation relations {M αβ, J i } = 0, {J i, J j } = ǫ ijk J k. Therefore, in addition to the involutive set of integrals H α we can choose J 2 = α<β M2 αβ and 2(J +ij 2 ) = x x = c, which gives the equation of the n sphere Generic Ellipsoidal Coordinates on S nc and E nc. Critical to the separation of variables on the n sphere S nc is the system of ellipsoidal coordinates graphically pictured by the irreducible block (2.) (S nc e e 2 e. where in general e α e β for α β. The separation variables are defined as zeros of the off diagonal element B(u) of the L matrix, i.e., B(u j ) = 0, j =,..., n. It follows that (2.2) x 2 α = 0 for u = u j and x 2 α u e = c Πn j= (u j e α ) α Π β α (e β e α ). Each vector of momentum s α is associated with a cell e α of the block. Note that x 2 α = 2(s α + is 2 β ). For each u j the conjugate variable v j is defined according to (2.3) v j = ia(u j ) = 2 x α p α u j e α. From (.5) one can show that the u i, v i satisfy the canonical relations [4]: (2.4) {u j, u i } = {v j, v i } = 0, {v j, u i } = δ ij. The change to the new variables v j, u i, c and J 3 is effectively the procedure of variable separation of (.) in ellipsoidal coordinates on the n-sphere. Writing the L matrix in terms of the new variables we obtain L(u) = ( A(u) B(u) ), B(u) = c Π n j= (u u j),

6 6 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. (2.5) A(u) = 2i n c B(u) ij 3 v j Π γ= + (u j e γ ), u u j Π i j (u j u i ) j= where A(u j ) = iv j, j =,, n, and A(u) (/u)j 3 + as u. To obtain C(u) we first notice that equating residues at e α on the right and left hand side of A(u) gives p α = 2x α c ij 3 + n v j Π e j= α u j γ= (u j e γ ). Π i j (u j u i ) This together with the expression (2.2) for x 2 α in terms of u j gives C(u) in the new variables. Three other useful formulae are (2.6) M αβ = 2x αx β c (e α e β ) n j= ˆΠ γ= v (u j e γ ) j Π i j (u j u i ), where the hat in (2.6) means that the product terms with γ = α and γ = β are omitted, (2.7) J 2 = n j= v 2 j Π γ= (u j e γ ) n Π i j (u j u i ) = j= cv 2 j ( x 2 α ), (u j e α ) 2 H α = x2 α c n vj 2 Π u j e α j= γ= (u j e γ ) Π i j (u j u i ). These relations together with (2.2) establish the explicit connection between the two sets of 2n+2 variables p α, x α and u i, v j, c, J 3. The equation for the eigenvalue curve Γ : det(l(u) iλi) = 0 has the form λ 2 A(u) 2 B(u)C(u) = 0. If we put u = u j into this equation then λ = ±v j. Thus variables u j and v j lie on the curve Γ: (2.8) v 2 j + n H α u j e α v 2 j det L(u j ) = 0. Equations (2.8) are the separation equations for each of the n degrees of freedom connected with the values of the integrals H α. For the sphere c α = 0 these have the form H α = M 2 αβ, Hα = 0, eα H α = M 2 αβ = J 2.

7 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 7 The Hamilton Jacobi equation (.) when parametrised by these variables has the form (2.9) 4 Mαβ 2 = J 2 = α<β n ( W 2Π γ= ) (u j e γ ) u j Π i j (u j u i ) = E, j= v j = W u j, which can be solved by the separation of variables ansatz (2.0) W = n W j (u j, H,..., H n ) = j= n j= v j du j. It will be convenient to employ an alternative form of the L matrix. If we use the vector (2.) L(u) = where s = we see that L satisfies s α, L(u) = (s,l(u)), det L(u) = L(u) L(u) ( ) 0, s 2 = 0 ( ) ( ) 0 i, s 3 0 =, i 0 0 (2.2) {L i (u), L j (v)} = ǫ ijk ( L k (u) L k (v) ). u v At this point we must consider a crucial difference between the real sphere and its complex counterpart. In the case of the complex sphere the generic ellipsoidal coordinates can admit multiply degenerate forms: the restriction e α e β, for α β can be lifted. The resulting coordinates can be denoted by the block form (S nc e λ eλ 2 2 eλ q q, λ + + λ q = n +, where the λ α denote the multiplicities of the e α. To understand how the previous analysis applies to these types of coordinates we first illustrate with an example corresponding to the coordinates with diagram In this case we write (S nc e 2 e 3 e. (2.3) L(u) = a s + a 2s 2 s α +. u e u e 2 Putting α=3 (2.4) x 2 x + ǫx 2, x x, p 2 p + ǫp 2, p p,

8 2 2 8 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. then in the limit as ǫ 0 we find where L(u) = z (u e ) + z u e α=3 z =< 4 (p 2 + x 2 i ), 4 (p 2 x 2 i ), 2 p x > z α z 2 =< 2 (p p 2 + x x 2 ), i 2 (p p 2 x x 2 ), i 2 (p x 2 + p 2 x ) > (2.5) z α = s α, α = 3,, n +. The components of z,z 2 satisfy the E(3,C) algebra relations {z i, zj } = 0, {zi, zj 2 } = ǫ ijkz k, {zi 2, zj 2 } = ǫ ijkz k 2 where the Poisson bracket (expressed in the primed coordinates) is {F, G} = j,k= B jk ( F p j G x k F x j G ) p k and B = (B jk ) = It can easily be verified from these relations that relations (2.2) are again satisfied. Thus, {L L(u),L L(v)} = 0 so the coefficients of the various powers (u e j ) k in the expansion of L L(u) form an involutive set of integrals of motion. We obtain L L = (u e ) 4z z + + α 3 2 (u e ) 3z z 2 + 2z 2 z α (u e )( ) + α,β 3, α β α=3 (u e ) 2z 2 z 2 + α 3 ( )(u e β ) (z α,z β ) (2.6) = 4(u e ) 2(p x 2 p 2x ) 2 + [ 4(u e ) 2 ( ) (p x α p αx ) 2 + 2(u e )( ) (p x α p α x )(p 2 x α p α x 2 )] n + 8( )(u e β ) (p β x α p α x β )2. α,β=3,α β In particular, 2z z α (u e ) 2 ( )

9 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 9 z α z 2 = (p x α p α x )(p 2 x α p α x 2 )/4. To relate this to the projective coordinates on the complex n sphere we recall that under the transformation (2.4) the fundamental quadratic forms X = a x 2 +a 2x α=3 x2 α, P = a p 2 + a 2p p2 α transform to X = 2x x 2 + α=3 x 2 α, P = 2p p 2 +. Therefore if we take the coordinates α=3 p 2 α x = 2 (X + ix 2 ), x 2 = 2 (X ix 2 ), x α = X α, α = 3,, n + p = 2 (P + ip 2 ), p 2 = 2 (P ip 2 ), p α = P α, α = 3,, n + we then can write p x α p αx = 2 (M α + im 2α ), p 2x α p αx 2 = 2 (M α im 2α ), p x 2 p 2 x = i 2M 2, X = Xj 2, P = Pj 2, where M jk = X j P k X k P j. The integrals of motion H α, Z in this case have, using partial fractions, the form L L = + j= j= (u e ) [ 2 4 (p x 2 p 2 x )2 + 4(e e α ) (p x α p α x )2 ] α=3 α=3 [ 4(e e α ) 2(p x α p α x )2 2(e e α ) (p x α p α x )(p 2 x α p α x 2 ) + (p α x β p β x α )2 ] + [ (p 4(e α e β ) (u e ) 4(e β α=3 e α ) x α p αx ) (e e α ) (p x α p α x )(p 2 x α p α x 2 )] = α=3 H α + Z (u e ) + Y 2 u e where Y = α=3 H α. The analysis presented so far could have been deduced from Kuznetsov [4] where the double root is essentially contained in the s systems of type C on the real hyperboloid. Furthermore, the threefold root is contained in Kuznetsov s type D systems. The question we now answer is how to use these techniques on the case of ellipsoidal coordinates corresponding to multiply degenerate roots. For this we use the limiting procedures developed by Kalnins, Miller and Reid [5]. We recall that the process of using these limiting procedures amounts to altering the elementary divisors of the two quadratic forms x 2 α = 0, x 2 α = c2.

10 0 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. Theorem. Let u i be the generic ellipsoidal coordinates coordinates on the n sphere viz. (2.8) a α x 2 α = c2 Π n j= (u j e α ) Π β α (e β e α ), α =,, n + with corresponding infinitesmal distance ds 2 = 4 n (du i ) 2 Π j i (u i u j ) Π j= (u i e j ) i= and coordinate curves (2.9) a α x 2 α u e j = 0, a α x 2 α = c2. Then the degenerate ellipsoidal coordinates having the infinitesimal distance (2.20) ds 2 = 4 n (du i ) 2 [Π j i (u i u j )] Π p j= (u i e j ) N j i= can be obtained from generic ellipsoidal coordinates via the transformations (2.2) e J j ej +J ǫ j, j =,, N J, J =,, p N J p J j pj + i=2 N J x J j xj + i=2 J ǫ i j+ i pj i J ǫ i j+ i xj i where J ǫ i j+ i = Πi l=2 (J ǫ j J ǫ l 2 ), aj j = /[Π k j( J ǫ j J ǫ k )], k =,, N J and N + + N p = n +. (We require J ǫ i 0 = 0 and take the limit as the J ǫ h 0 for h =,, N J.) In particular, (2.22) L(u) = where p N J J= j= z J j (u e J ) N J j+ z J j =< (p J i 4 pj j+ i +xj i xj j+ i ), i (p J i 4 pj j+ i xj i xj j+ i ), i p J i 2 xj j+ i >. i i i

11 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET (2.24) {(z J j ) l, (z J k ) m } = δ JJ (z J j+k N J ) s ǫ lms and we also have (2.25) L L(u) = 4 (x J j p L k xl k pj j )(xj m pl n xl n pj m ) (u e J ) p (u e L ) q where p N J, q N L, j, m N J, k, n N L, subject to j +m = N J +2 p, k +n = N L +2 q and the summation extends over the indices J, L, p, q, j, m, k, n subject to these restrictions. From expressions (2.24) we can verify that relations (2.2) are satisfied, so that {L L(u),L L(v)} = 0. The constants of the motion can be read off from the partial fraction decomposition of (2.25). This clearly illustrates the compactness of the r matrix formulation for the operators describing the integrable systems examined so far. In dealing with the case of Euclidean space E nc the most transparent way to proceed is as follows. The generic ellipsoidal coordinates in n dimensional complex Euclidean space are given by [,9] (2.26) x 2 α = Πn j= (u j e α ) Π β α (e β e α ), j, α, β =,, n with coordinate curves n x 2 α (2.27) ( ) =, u = u j, j =,, n. Proceeding in analogy with (2.2)-(2.0) we can obtain the r matrix algebra. The corresponding L(u) operator is n s α (2.28) L(u) = +. u e α 4 0 (Indeed the equation L (u) + il 2 (u) = 0 is just (2.27). Moreover it is obvious, due to the fact that expressions (2.) satisfy (2.2), that expressions (2.28) also satisfy (2.2).) The conjugate variables v j are defined by v j = il 3 (u j ) and they must satisfy the canonical relations (2.4). The integrals of motion H α are determined from n (2.29) L 2 H α (u) = where H α = 2 n (sα,s β ) e α e β 2 (s α + s2 α ) = n 4 ( β= β= M 2 αβ e α e β p 2 α ) with α H α = p 2 /4. The separation equations are of the form v 2 j + L2 (u j ) = 0, j =,, n,

12 2 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. 3. Cyclidic Coordinates. Associated with the separation of variables problem for the Hamilton Jacobi equation (.) with E 0 is the corresponding E = 0 problem. In this case the equation is (3.) n α,β= g αβ p α p β = 0, p = W x α, α =,, n and we consider only complex Euclidean space. While it is true that all the coordinate systems discussed for E nc with E 0 will provide a separation of variables of this equation, there are coordinates that provide an additive separation of variables only when E = 0. This is related to the fact that the E = 0 equation admits a conformal symmetry algebra [,,3]. The most convenient way to proceed is to introduce hyperspherical coordinates {x,, x n+2 }, n n x = t 2 ( zj 2 ), x 2 = it 2 ( zj 2 + ), x k+2 = 2z k t 2, k =,, n, j= j= related to the usual Cartesian coordinates {z,, z n } according to z k = x k+2 /( x ix 2 ), k =,, n n+2 x 2 j = 0. j= We consider the system of in n + 2 dimensions, where J = 0 and the c α = 0. The general (separable) cyclidic coordinates are specified by Ω = Furthermore, n+2 j= x 2 n+2 j = 0, Φ = λ e i j= x 2 j = 0, λ = u,, u n, e k e j. σx 2 α = Πn j= (u n+2 j e α ) Π β α (e β e α ), σ = [ e j x 2 j]. The Hamilton Jacobi equation is given by n ( W (3.2) J 2 = σ 2 u j j= j= ) 2 Π n+2 γ= (u j e γ ) Π i j (u j u i ) = 0. The quadratic forms Ω and Φ have elementary divisors [ ], see [5,]. It is known that the geometry of these fourth order coordinate curves is unchanged under birational transformations of the form e k αe k + β γe k + δ, u j αu j + β γu j + δ, λ αλ + β γλ + δ,

13 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 3 Now we can mimic the exposition given for the Gaudin magnet integrable systems using the hyperspherical coordinates and the Poisson bracket n+2 (3.3) {F(x i, p j ), G(x i, p j )} h = ( F G F G ), p k x k x k p k k= the x i, p j now being regarded as independent. The analysis then proceeds much as in the construction (2.2)-(2.2), but with n replaced by ; indeed the coordinates are defined as zeros of the off diagonal element B(u) = n+2 x 2 α subject to the restriction Φ = 0. The crucial difference is that J = 0. The expressions for p α and M αβ are altered by a factor σ: p α = 2σx α [ n j= v j Π n γ=(u j e γ ) e α u j Π i j (u j u i ) ] which together with x 2 α in terms of u j gives C(u) in the new variables. Another useful formula is n ˆΠn+2 γ= M αβ = 2σx α x β (e α e β ) v (u j e γ ) j Π i j (u j u i ). In fact the Poisson bracket {, } h can be identified with the Poisson bracket {, } for functions defined in the n dimensional space spanned by z,, z n. This can readily be seen by noting that F(z,, z n, p z,..., p zn ) = F( x 3 /(x + ix 2 ),, x n+2 /(x +ix 2 ), (x +ix 2 )p x3 +x 3 (p x +ip x2 ),, (x +ix 2 )p xn+2 + x n+2 (p x +ip x2 )) from which the equality of the Poisson brackets follows identically. The infinitesmal distance for general cyclidic coordinates is ( (3.4) ds 2 = n ) 4 (x + ix 2 ) 2 (du i ) 2 Π j i (u i u j ) Π n+2 j= (u. i e j ) i= j= We denote the coordinates defined by this graph as [E nc, E = 0 e... e n+2. Other coordinates of this type are those corresponding to the graphs (S pc e... e p+ (S qc f... f q+ where p + q = n + 2. These coordinates are given by (3.5) σx 2 k = Πp j= (u j e α ) Π β α (e β e α ), α, β =,, p +, k =,, p σx 2 k = Πq j= (v j f α ), α, β =,, q +, k = p +,, p + q + 2 Π β α (f β f α ) p+ q+ σ = [ ei x 2 i + fj x 2 j ].

14 4 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. We note here that the e i the f j are pairwise distinct, for if they were not then for any of these quantities which occured with multiplicity more than a birational transformation could transform it to and hence to a graph corresponding to E nc. It can happen just as in the case of generic ellipsoidal integrable systems on the sphere that some of the e i in (3.5) are equal to some of the f j also. In this case the rules for obtaining the corresponding L matrix are summarised in the following theorem. Theorem 2. Denote the generic ellipsoidal coordinates by the graph [E nc e... e n ], e i e j and generic cyclidic coordinates by the graph {CE nc e... e n+2 }, e i e j. Separable coordinates for the Hamilton Jacobi equation (3.) corresponding to generic graphs with multiplicities [E nc e n en p p ], n + + n p = n, {CE nc e n en q q }, n + + n q = n + 2, [5] can be obtained via the transformations of Theorem applied to the quadratic forms m a i x 2 m i Ω = = 0, Φ = a i x 2 i λ e = 0 i i= where m = n for generic ellipsoidal coordinates in E nc and m = n + 2 for the corresponding cyclidic coordinates. For coordinates corresponding to the direct sum of two spheres viz. (3.5) we can merely apply the result of the previous theorem for spherical coordinates to each of the pairs of quadratic forms i= p+ Ω = a α x 2 α, Φ u e = p+ a α x 2 α α n+2 Ω 2 = α=p+2 a α x 2 n+2 α, Φ 2 u e = α α=p+2 a α x 2 α. Indeed the freedom to subject the coordinates u j and the e i to birational transformations in the expressions for generic cyclidic coordinates allows us to let e. (In this particular case the resulting coordinates can be identified with generic elliptical coordinates on the n sphere.) The process described in Theorem enables one to pass from the elementary divisors [...] to [N, N 2,..., N p ], N i, i =,..., p, see [5,]. It is then possible (via Theorem 2) to take e in which case [ N, N 2,..., N p ] corresponds to the various generic coordinate systems in Euclidean

15 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 5 consider the quadratic forms ˆΩ, ˆΦ corresponding to elementary divisors [2 ] viz. ˆΩ = x 2 (λ e ) 2 + 2x x 2 λ e + + x2 n+2 λ e n+2 = 0, ˆΦ = 2x x 2 + x x2 n+2 = 0. Putting λ /λ, e i /e i we find (with the use of ˆΦ = 0) or, in the limit as e, ˆΩ = λ 2 e 2 [ x2 n+2 (e λ) 2 + x 2 n+2 k=3 k=3 (e e k ) (e λ) x 2 k λ e k = 0. x 2 k (e k λ) ], Now if we perform the limiting procedure (2.3) - (2.5) on L(u) (in n + 2 dimensions), let u /u, e i /e i, let e and then evaluate at p = 0, x = we find (up to a common factor u) L(u) = n+2 α=3 s α Here, we have made use of the fact that J = 0. This agrees with (2.28) and with the L(u) operator given by Kuznetsov [4], and corresponds to elementary divisors [ 2,,, ]. This process can be generalised for quadratic forms corresponding to elementary divisors [ N,, ], N >. We have that. ˆΩ = u N 2 x 2 +un 3 2x x 2 + i x i x N i n+2 k=n+ x 2 k u e k, ˆΦ = N x i x N+ i + i= n+2 k=n+ x 2 k. The corresponding L(u) operator is where L(u) = n+2 α=n+ s α u e k + s s = 4 [ˆΩ(0, p 2,..., p n+2 ) ˆΩ(, x 2,..., x n+2 )], s 2 = 4 [ˆΩ(0, p 2,..., p n+2 ) + ˆΩ(, x 2,..., x n+2 )], s 3 = 2 [S (0, p 2,..., p n+2 ;, x 2,..., x n+2 )], S = N j=2 u N j l+m=j x l p m.

16 x 2 = (u e )(u 2 e ) (ū 3 f 4 )(ū 4 f 4 ). 6 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. 4. Branching Rules for the Construction of Orthogonal Non Generic Complex Integrable Systems on S nc and E nc. To deal with the non generic separable coordinate systems in E nc and S nc we must combine coordinate systems for both manifolds [5]. (See [3,4] for tabulations of all cases for small values of n.) The branching laws for graphs on these manifolds are summarised below () (S nc e i S pi C (2) (3) (S nc e λ i i, λ i > E pi C (E nc e i S pi C (4) (E nc e λ i i, λ i >. E pi C As an example consider the coordinate system given by the graph on S 4C, (S 2C e 2 e 5 [E 2C f 3 f 4 The coordinates for this graph can be obtained from those of the generic graph via the limiting transformations (S 4C e 2 e 3 e 4 e 5 e j = e + ǫ + ǫ 2 f j, j = 3, 4, f 3 f 4 u j = e + ǫ + ǫ 2 ū j, j = 3, 4 where ǫ 0. The corresponding coordinates are then given implicitly by x 2 = (u e )(u 2 e ), x 2 5 e 5 e = (u e 5 )(u 2 e 5 ) (e e 5 ) 2 2x x 2 = (u e ) (u 2 e ) + (u e )(u 2 e ) e 5 e e 5 e (e 5 e ) 2 (u e )(u 2 e ) e 5 e (f 3 + f 4 ū 3 ū 4 ) x 2 3 = (u e )(u 2 e ) (e 5 e ). (ū 3 f 3 )(ū 4 f 3 ), (f 3 f 4 )

17 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 7 Using the nomenclature given previously we have that L(u) has the two different forms z L (u) = (u e ) + z 2 + z 3 + z 4 + z 5, u = u 2, u 2 u e u e u e 5 L 2 (u) = z + z 3 u f 3 + z 4 u f 4, u = ū 3, ū 4. As usual, the separation variables are the zeros of the equation L λ (u)+il2 λ (u) = 0 with u = u, u 2 when λ = and u = ū 3, ū 4 when λ = 2. Each of the L λ separately satisfy (2.2) for λ =, 2. In addition we have {L L (u),l 2 L 2 (v)} = 0. This example illustrates how to derive the substitutions that enable the various branching laws to be obtained from a generic form. For the sphere S nc and generic coordinates (S nc... e λ+q s 0 e e s q > s, λ >, if we make the substitutions e k = e 0 + ǫ q s+ + ǫ λ+q s f k, k = s, u j = e 0 + ǫ q s+ + ǫ λ+q s ū j, j = q where ǫ 0, then the graph illustrated transforms into (S nc e λ 0. [E qc f f s The remaining branching rule is obtained from a graph of the type By means of the substitutions e j = e 0 + ǫ(f 0 f j ), j =,, p + (S nc... e 0 e... e s. u k = e 0 + ǫ(ū k f 0 ), k =,, p we obtain the coordinate system coming from the graph (S nc e 0. (S pc f 0 f p The corresponding substitutions for the analogous Euclidean space branching rules are essentially identical. To completely specify the coordinate systems on these manifolds we need a few more substitution rules. Firstly consider the graph [E pc h... h p + [E qc f... f q. This graph can be obtained from the generic graph for E (p+q)c via the substitutions u i = K ǫ + ū i, e i = K ǫ + h i, i =,, p u i = K 2 ǫ + ū i, e i = K 2 ǫ + f i p, i = p +,, p + q, K K 2. The graph [E qc h... h q, q < n can be obtained from [E nc h... h n via the substitutions u i = e + ǫ q n + ǫ q n ū i, i =,, n e j = e + ǫ q n + ǫ q n h j, j =,, q e k = h k, k = q +,, n. Given these substitution rules all the corresponding graphs for E nc and S nc can

18 8 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. 5. Quantum Integrable Systems on Complex Constant Curvature Spaces and the Quantum Gaudin Magnet. To deal with the quantum version of this description of separation of variables we consider the Schrödinger or Helmholtz equation [3]. (5.) HΨ = [ g g αβ g y α ]Ψ = EΨ, y β where E 0 for the moment. Separation of variables means (roughly) the solution of this equation of the form (5.2) Ψ = Π n ψ α(y α ; h,, h n ) where the quantum numbers h j are the eigenvalues of mutually commuting operators n A i = a (j) 2 n (5.3) αβ + b β, j =,, n y α y β y β α,β= β= [A i, A j ] =0, i, j =,, n, A n = H, A j Ψ = h j Ψ. Furthermore these operators can be represented as symmetric quadratic elements in the enveloping algebra of the symmetry algebra of the equation (5.) in the case of S nc and E nc. The standardised representations of these symmetry algebras are (5.4) (5.5) SO(n + ) : M αβ = ˆx αˆp β ˆx β ˆp α, ˆp α =, α, β =,, n +. x α [M αβ, M γδ ] = δ αγ M δβ + δ αδ M βγ + δ βγ M αδ + δ βδ M γα E(n) : M αβ, P γ = x γ [M αβ, P γ ] = δ βγ P α δ αγ P β, [P α, P β ] = 0 where [, ] is the commutator bracket. Much of the analysis goes through as it did in the classical case with, of course, some critical differences. For the quantum Gaudin magnet one considers the sum of rank Lie algebras A = m so α (3) where the generators of the algebra satisfy the commutation relations (.6) and the inner product is defined as in. The Casimir elements of A have the form (5.6) (s α,s α ) = k α (k α + ) where k α is a constant when the generators of so α (3) determine an irreducible representation,[2]. The quantum Gaudin magnet is the quantum integrable Hamiltonian system on A given by m commuting integrals of motion : H α = 2 m (sα,s β ), [H α, H β ] = 0. e α e β β= This is the m site SO(3, C)-XXX quantum Gaudin magnet, [3]. These integrals (operators) satisfy (5.7) m Hα = 0, m eα H α = J 2 m kα (k α + )

19 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 9 where J = m s α is the total momentum operator. The complete set of commuting operators consists of H α,j 2 and J 3. The integrals are generated by the 2 2 operator L(u) given by (.) understood in the operator sense. The quantum determinant is (5.8) q-detl(u) = A(u) 2 2 {B(u), C(u)} = m with L-operator satisfying the r matrix algebra H α (5.9) [ L(u), 2 L(v)] = i u v [P, L(u) + 2 L(v)], i =, (compare with (.3)). The total momentum J has components m k α (k α + ) ( ) 2 (5.0) J = 4 (ˆp ˆp + ˆx ˆx), J 2 = i 4 (ˆp ˆp ˆx ˆx), J 3 = i (ˆp ˆx + ˆx ˆp) 4 The Lie symmetries M αβ and J form the direct sum so(m) so(3) with commutation relations (5.) [M αβ, J i ] = 0, [J i, J j ] = iǫ jkl J l, subject to the constraints (5.2) J 2 = 4 α<β M 2 αβ + 6 m(m 4), 2(J + ij 2 ) = ˆx ˆx. Considering again separation of variables using the coordinates of the irreducible block (2.), the separation variables are defined, as before, as the zeros of the off diagonal elements B(u) of the L matrix. The q-determinant is the generating function of the commuting integrals of motion (5.3) [q-detl(u), q-detl(v)] = 0. On the n sphere the algebra A is realised by taking the canonical operators (5.4) s α = 4 (ˆp2 α + ˆx2 α ), s2 α = i 4 (ˆp2 α ˆx2 α ), s3 α = i 4 {ˆp α, ˆx α }, α =,, n + where (s α,s α ) = 3 6. Furthermore (5.5) (s α,s β ) = 8 (M2 αβ + 2 ) and this establishes the relationship between the s α of A and the M αβ of so(). The integrals transform into the family of integrals (5.6) H = hα H α = h α h β (M 2 αβ + )

20 20 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. and the coordinates are given by B(u j ) = 0 : x 2 α = 0. For each u j variable there is defined the conjugate variable (operator) (5.7) v j = ia(u j ) = 4 u j e α {ˆp α, ˆx α }. Separation of variables is then the process of changing to the new variables v j, u i, c and J 3. In fact we have that (5.8) 4 {ˆp α, ˆx α } = ˆx2 α c [ ij 3 n j= e α u j D j v j ], D j = Π (u j e α ) Π i j (u j u i ). The separation equations can be obtained by substituting u = u j into the q-detl(u), making the choice (5.9) v j = i( + ), u j 4 u j e α and looking for the solutions of the spectral problem H α Ψ = h α Ψ. The separation equations are then (5.20) d 2 du 2 j ψ j + 2 ( n ) d h α ψ j = ψ j u j e α du j u j e α where Ψ = Π j= ψ j(u j ), h α = h α +. 8 e α e β β α These are the separation equations for the Helmholtz equation in generic ellipsoidal coordinates. Note that for this choice of v j we have taken (s α,s α ) = 3/6. The L operator is given in direct analogy with the classical case, (5.2) L(u) = s α. The treatment of generic ellipsoidal coordinates on E nc follows similar lines, with the L operator given in this case by (5.22) L(u) = n s α The limiting procedures described in the classical case work also in the quantum.

21 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 2 In the case of cyclidic coordinates we can adopt the same strategy as in 3: we impose the conditions J = 0 and proceed as before. The natural setting in this case is again to use hyperspherical coordinates. The total momentum has components as in (5.0). If ˆp α, ˆx α are now vectors in hyperspherical coordinates then we can derive the standard quantum r matrix algebra as above with the same formulas valid. The constraints are now 2(J + ij 2 ) = ˆx ˆx = 0, J 2 = 8 n(n 4). Coordinates u j and their conjugate operators v j can be chosen as before. In particular if we make the choice (5.23) v j = + n+2 u j 4 u j e α and look for solutions of HΨ = 0 of the form Ψ = σ (n 2)/2 Π n i= ψ i(u i ), then we obtain the separation equations which coincide with the the equations given by Bocher,[], viz. d 2 du 2 ψ j + n+2 j 2 u j e α d du j ψ j = 6 (n2 4)u n j 6 (2n n2 )( n+2 e α)u n j + n 2 β=0 λ βu β j Π n+2 γ= (u ψ j j e γ ) for suitable γ β. Note that the solutions Ψ are not strictly separable in this case but are what is termed R separable, i.e., separable to within a non separable factor σ (n 2)/2, []. We conclude this work by pointing out that one can obtain a complete set of constants of the motion associated with an orthogonal separable coordinate system {u, u n } for the Schrödinger equation (5.) directly from (5.) itself, [5, 6]. It is well known that all orthogonal separable coordinate systems for the Schrödinger equation on a Riemannian space are obtainable via the Stäckel construction, e.g., [7, 8]. Thus if {u} is a separable orthogonal coordinate system for (5.) there exists an n n nonsingular matrix S = (S γ,β (u α )) such that uγ S αβ = 0 if γ α, and such that the nonzero components of the contravariant metric tensor in the coordinates {u} are g αα (u) = T nα (u), α =,, n, where T is the inverse matrix to S, (5.24) n T αβ (u)s βγ (u β ) = δγ α. β= The constants of the motion are n (5.25) A β = T βα (u) ( u 2 ) α + f α (u α ) uα, An = H. Here, f α (u α ) = uα ln( gg αα ). (The fact that αα

22 22 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. follows for any space of constant curvature by noting that it is equivalent to the statement that the off diagonal elements of the Ricci tensor must vanish for an orthogonal coordinate system: R αβ = 0, α β, [7].) One can show that [A α, A β ] = 0, [7, 8]. Furthermore, if Ψ = Π n ψ α (y α ; h,, h n ) satisfies the separation equations (5.26) 2 u α Ψ + f α (u α ) uα Ψ = then n S αβ (u α )h β Ψ, α =,, n β= (5.27) A β Ψ = h β Ψ, β =,, n. Now fix γ, γ n, and denote by y the coordinate choice y = (u,, u γ, τ, u γ+,, u n ) where τ is a parameter. We see from (5.26) that if Ψ satisfies the separation equations then (5.28) ψ γ ( 2 uγ ψ γ + f α (u γ ) uγ ψ γ ) uγ =τψ(u) = On the other hand, from (5.) we have where ψ γ ( 2 uγ ψ γ + f α (u γ ) uγ ψ γ ) uγ =τψ(u) = H γ (y) = g(y) This suggests the operator identity (5.29) α γ g γγ (y) (H(u) Hγ (y)) = n S γβ (τ)a β Ψ(u). β= g γγ (y) (H(u) Hγ (y))ψ(u) uα ( g(y)g αα (y) uα ). n S γβ (τ)a β, i.e., that the left-hand side of (5.29) is a one-parameter family of constants of the motion and that as τ runs over a range of values and γ =,, n the full space of constants of the motion associated with this separable system is spanned, [5, 6]. Now (5.29) is equivalent to the conditions T nα (u) T nγ (y) T nα (y) T nγ (y) ( δα γ ) = β= n S γβ (τ)t βα (u), β= and these conditions are easily seen to follow from (5.24) and from n β=,β γ T nβ (y)s βξ (u β ) + T nγ (y)s γξ (τ) = δ n ξ. Similarly, for the classical case an expression for the Hamiltonian analogous to

23 SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 23 References. E.G. Kalnins, Separation of Variables for Riemannian Spaces of Constant Curvature, Pitman Monographs and Surveys in Pure and Applied Mathematics 28, Longman Scientific and Technical, Essex, England, M. Gaudin, La fonction d onde de Bethe, Masson, Paris, V.B. Kuznetsov, Equivalence of two graphical calculi, J. Phys. A Math. Gen. 25 (992), V.B. Kuznetsov, Quadrics on real Riemannian spaces of constant curvature. Separation of variables and connection with Gaudin magnet, J. Math. Phys. 33 (992), E.G. Kalnins, W. Miller Jr. and G.J. Reid, Separation of variables for complex Riemannian spaces of constant curvature. Orthogonal separable coordinates for S nc and E nc, Proc. Roy. Soc. Lond. A394 (984), I.V. Komarov and V.B. Kuznetsov, Quantum Euler-Manakov top on the 3-sphere S 3, J. Phys. A: Math. Gen. 24 (99), L E.K. Sklyanin, Separation of variables in the Gaudin model, J. Soviet Math. 47 (989), I.V. Komarov and V.B. Kuznetsov, Kowalewski s top on the Lie algebras o(4), e(3), and o(3,), J. Phys. A: Math. Gen. 23 (990), W.D. Niven, VI. On ellipsoidal harmonics, Phil.Trans. CLXXXII (89), E.G. Kalnins and W.Miller Jr., Separation of variables on n-dimensional Riemannian manifolds.. The n-sphere S n and Euclidean n-space R n, J. Math. Phys. 27 (986), 72.. M. Bôcher, Die Reihentwickelungen der Potentialtheorie, B.G. Teubner, Leipzig, N.J.Vilenkin, Special Functions and the Theory of Group Representations, Translations of Mathematical Monographs Vol22 American Math.Society, Providence, Rhode Island, E.G. Kalnins and W. Miller Jr., The wave equation and separation of variables on the complex sphere S 4, J. Math. Anal. Appl. 83 (98), E.G. Kalnins and W. Miller Jr., Lie theory and the wave equation in space time. 2. The group O(4, C), SIAM J. Math. Anal. 9 (978), N.M.J. Woodhouse, Killing tensors and the separation of the Hamilton-Jacobi equations, Commun. math. Phys. 44 (975), W. Dietz, Separable coordinate systems for the Hamilton-Jacobi, Klein-Gordon and wave equations in curved spaces, J. Phys. A(GB) 9 (976), L.P. Eisenhart, Separable systems of Stäckel, Ann. Math. 35 (934), W. Miller Jr., Mechanisms for variable separation in partial differential equations and their relationship to group theory, Proceedings of School on Symmetry and Nonlinear Phenomena, Paipa, Colombia, February 22-26, 988 in Symmetries and Nonlinear Phenomena, D. Levi and P. Winternitz eds., World Scientific, London (989). 9. E.G. Kalnins and W. Miller, Jr., The theory of orthogonal R-separation for Helmholtz equations, Advances in Mathematics 5 (984), E.G. Kalnins and W. Miller, Jr., Intrinsic characterization of orthogonal separation of one coordinate in the Hamilton-Jacobi equation, J. Phys. A: Math. Gen. 5 (982), L.D. Faddeev and L.A. Takhtajan, Hamiltonian methods in the theory of solitons, Springer Verlag, Berlin, 986.