arxiv:math-ph/ v1 27 Feb 2007

Transcription

1 Goldfish geodesics and Hamiltonian reduction of matrix dynamics arxiv:math-ph/07009v 7 Feb 007 Joakim Arnlind, Martin Bordemann, Jens Hoppe, Choonkyu Lee Department of Mathematics, Royal Institute of Technology, Stockholm, Sweden. Laboratoire de MIA, 4, rue des Frères Lumière, Univ. de Haute-Alsace, F Mulhouse. Dept. of Physics and Center for Theoretical physics, Seoul National Univ., Seoul 5-74, Korea. Abstract We relate free vector dynamics to the eigenvalue motion of a time-dependent real-symmetric N N matrix, and give a geodesic interpretation to Ruijsenaars Schneider models. Despite of more than two decades of extensive work on Ruijsenaars-Schneider models [, ], their interpretation as describing geodesic motions seems to have gone unnoticed. Apart from wishing to fill this gap, the related) second topic of this paper is the Hamiltonian reduction of free matrix dynamics to, as we will show, free vector dynamics). To start with the latter, let Ẍt) = 0, Xt) a real symmetric N N matrix depending on time ), be described by Writing H[X,P] := TrP. ) Xt) = Rt)Dt)R t) ) Ẋt) = R Ḋ +[M,D] R = RLR, ) the symplectic form TrdX dp becomes cp. [3]), with da := R dr, dq i dp i + df ij da ij i<j f ij da jk da ik +f ik da ij da jk +f jk da ik da ij ), i<j<k 3) which inverting 3)) gives the non-trivial Poisson-brackets {q i,p j } = δ ij 4) {f ij,f kl } = δ jkf il + δ ikf jl + δ jlf ik δ ilf jk 5) {f i<j,a k<l } = δ ikδ jl 6) resp. {r ij,f kl } = δ jkr il δ jl r ik ) 7)

2 for the eigenvalues of) X, their time derivatives q i, and f ij = f ji := R Ṙ q i q j ). ) becomes ij H = N p i + fij q i q j ), 8) i= which, as H = H[ q, p;f ij ], is known under the name Euler Calogero-Moser Hamiltonian [4, 5]), with equations of motion i j q i = k i f ij = k i,j f ik q i q k ) 3 f ik f kj q ik ) qkj. 9) While it is well known see e.g. [6, 7]) that both types of) equations in 9) consistently reduce to q i = q i q j 0) q i q j upon setting j i) f ij = q i q j ) q i q j, ) the Hamiltonian reduction of ) resp. 8)) to 0) has remained open for many years. Simply counting the degrees of freedom H[ q, p;f i<j ] has at most N + NN ) ; to obtain a nice phase-space one would have to eliminate the degrees of freedom corresponding to the Casimirs of the SON) generated by the f ij ), together with {G ij,g kl } = δ jk G il +δ ik G jl +δ jl G ik δ il G jk G ij := f ij +q i q j ) p i p j ) ) shows that one can not simply fix the gauge. Viewing 8) as originating from ), however, one has the degrees of freedom corresponding to the orthogonal matrix R resp. the antisymmetric matrix A, which due to 6) naturally provides NN )/ variables that are canonically conjugate to the G ij ). Note that 8), when written in terms of the variables G ij, reads H = N ) p i + G ij 8 q i= i q j ) G ij pi p j, 3) 4 q i q j i j i j clarifying in the Hamiltonian framework) that it is consistent to put the G ij = 0, due to Ġij = {G ij,h} being weakly zero i.e. using G kl = 0 after computing the Poisson-brackets, according to ) and 4)/5); note that while the q s and p s do not commute with the G ij, the total momentum, P := p i, does, as the G ij contain only the differences of the q s). In order to find out what the reduced phase space is consisting of functions that Poisson-commute with all the G ij ), the reduced Hamiltonian cp. 3)) being H = ) pi = P, 4)

3 it is useful to bring the orbital angular momentum, L ij := q i q j ) p i p j 5) satisfying ), with G replaced by L) into canonical form Q i P j Q j P i ), by making the canonical transformation P i = p i, Q i = q i pi. 6) While the first order form, q i = p i, ṗ i =..., of the) goldfish-equations 0) does) not simplify at all, Q i = P 3 i +Q i P i P i = Pi Pj 3 Q i P j Q j P i j j P 3 j Q i P j Q j P i 7) with Q i P i Q i P i = Pi 4) Q and P now transform as ordinary vectors Q S Q, P SP) under the rotations generated by the G ij : {G ij,q k } = δ ik Q j δ jk Q i {G ij,p k } = δ ik P j δ jk P i. 8) Together with cp. 7)) {G ij,r kl } = δ il r kj δ jl r ki, 9) corresponding to R RS, the N independent variables u := R Q, v := R P 0) are therefore invariant, i.e. natural coordinates for the reduced phase-space. The original free matrix-dynamics with the initial condition that Ẋ0) has positive eigenvalue, and N zero) is thereby reduced to free vector dynamics, u = v v ), v = 0, ) governed by H[ u, v] = v ) = RP) ) = P ) = ) P. ) Using 7), and the time-evolution of the r ij cp. 7)/8)) given according to ) R Ṙ = M ij = f ij ij qij pi p j = q ij P i P j = 3) Q i P j Q j P i one can check the time-evolution of the invariant variables u = R Q, v = R P: d R dt ) Q = RP) P, 3 d RP dt ) = 0. 4)

4 Let us know focus on 0), independent of any previous considerations, as coupled) second order ODE s. Their most natural, and extremely simple though apparently unnoticed), interpretation is that of geodesic equations, q i +Γ i jk q j q k = 0 5) δ i Γ i jk = j δk i) q i q k + δi k δi j ) ) q i q j. 6) The simple expression 6), which is of the form Γ i jk = δi j w ik +δ i k w ij, 7) with w ik = wqi q k ) δ ik ), makes the calculation of the curvature-tensor Radb c = dγ c ab +Γe ab Γc de b d) straightforward: R c adb =δc a δc d w ab δc a δc b w ad +w ca δ abδ c d δc b δ da) +δ c dw ab w ca +w ba w cb w ca w cb ) δ c b w adw ca +w da w cd w ca w cd ) 8) is identically zero for wx) = /x i.e. 6)), due to + q c q a q a q d + q a q d q d q c = 0 q d q c q c q a a c d a), 9) and the w term being cancelled by those ww terms, for which an additional equality of indices holds; alternatively, one could e.g. bring the last 3 ww-terms in 8), provided w x) = wx), into the form wx)wy) wx+y)wx)+wy)), 30) x = q c q a,y = q a q d, and then note the linearity of /w. In order to find the change of variables that will make the Christoffel-symbols vanish, it is easiest to observe that q i = j q i q j q i q j directly implies that the functions b n t) := n )! qiq i q in q in n =,,... 3) with indicating the indices to all be different) do not depend on time: n )! + n )! qiq i q in q in qiq i q in j i n q in q j q in q j = 0, 3) 4

5 as the sum over j / i,...,i n,i n ) gives zero, while q n times the sum over j i,...,i n ) can be replaced by n q in q in. The searched transformation q xq,...,q n ) non-singular, as long as the q i are all different) is therefore provided by The natural induced metric giving 6)) is with x n [ q] := qi q in, 33) n! ) xn det = i q j ), 34) q i i<jq ẍ n [ qt)] = ḃn = 0. 35) g ij = i x j x = J T J ) ij, 36) J n j = xn q j = n )! q j q jn, 37) now writing the q coordinates with upper indices) and the goldfish-equations can therefore be described by the geodesic flow) Hamiltonian j H[ q, π] := π ig ij [ q]π j, 38) g ij = J J T ) ) ij +q i q j + +q i q j ) N = k qi q k ) l qj q l ), 39) q i = g ik [ q]π k, π k = π g ij i q k π j 40) giving 5)/6). To the best of our knowledge, this geodesic interpretation and natural quadratic Hamiltonian structure) of the Ruijsenaars-Schneider model has not been observed before. The inverse of the Jacobian is J ) im = )m q i ) N m j i qi q j ), 4) and the conserved quantities B n [ q, π] := n )! q i q in g inj [ q]π j 4) all linear in the momenta!) Poisson-commute, as the transformation x, p) q, π), with p n = ẋ n = B n trivially gives {p m,p n } = 0. Hamiltonian relations between zeroes of polynomials and their coefficients have been considered before cp. [8, 0]), but to our suprise apparently not for the original goldfish equation. We also became aware of chapter 7 in [9], where just as in [0]) the time-independence of the quantities 3) is thanks to E. Langmann! 5

6 stated and proved), but with the standard exponential cp. 43)) Hamiltonian structure, and not the geodesic structure with respect to which the conserved quantities Poisson-commute). Note that, rewriting [9] the canonical Poissonbracket, {q i, p j } = δ ij, in terms of the exponential variables suggested by the standard RS Hamiltonian one has N i= e pi N j= q i q j ) =: π i = P, 43) {q i,q j } = 0, {q i, π j } = δ ij π i, { π i, π j } = π i π j q i q j δ ij ) 44) which interpreting them as Dirac-brackets, in the context of our Hamiltonian reduction to H = p i ), with p i playing the role of π i ) tells one how to obtain the goldfish equations directly from the reduced Hamiltonian 4): using 44), with π i replaced by p i, and 43) by H = P, gives 0). Finally, let us make some remarks about the hyperbolic case: Starting, as was done for the Calogero-Moser case in [], with d ẊX +X Ẋ) = 0, 45) dt X a positive definite matrix which we take to be real, with eigenvalues e aλjt) a real), Xt) = Rt)e aλt) R t), 46) R T = R, one obtains with, as before, M := R Ṙ) from 45), when defining L via L ij = δ ij λi sinhaλ i λ j )) M ij, a 47) L = [L,M] 48) Rt)Lt)R t) = ẊX +X Ẋ). 49) 4a The solution of 45), on the other hand, can be written as Xt) = e aλ0 e tv0 e aλ0 50) with Λ 0 = Λ0) when choosing R0) = ). The crucial point is that 48) consistently reduces to λ i = j i) a λ i λj sinhaλ i λ j )) 5) when setting a M ij = λ i λj 5) sinhaλ i λ j )) 6

7 corresponding to the initial conditions ci c j V 0 ) ij = a coshaa i a j ), 53) where a i = λ i 0), c i = λ i 0) geometrically speaking, the matrix geodesics 46) project down to geodesics in the space of positive diagonal matrices, if the initial conditions are chosen perpendicular to the action of the rotationgroup ). Except for their intrinsically geodesic interpretation made explicit below) the equations 5), and variants thereof, are well-studied see e.g. [0, ]). In particular one can show, that the solutions of q i = j q i q j cothq i q j ) 54) are the N roots of fq) := N i= c i P tanhpt) = 0, 55) tanhq a i ) resp. see e.g. [9]) the eigenvalues of Zt) = e Λ0 e tl0, with P = q i = c i, L 0 ) ij = c j avoiding the somewhat artificial positivity restrictions in 53)). Rewriting 55) as d e tp d dt dt det Zt) e q )) = 0, 56) the crucial step then not taken in [9]) is to note that the symmetric functions of the eigenvalues e qit), s n t) = n! e q i + +qin ) n =,...,N, 57) inanalogywith33), evolveintimealmostuncoupledasp = ṡ N /s N = const.) as s n Pṡ n = 0 n =,...,N. 58) Acknowledgement We would like to thank the Swedish Research Council, the Brainpool program of the Korea Research Foundation and the Korean Federation of Science and Technology Societies, R , and the Marie Curie Training Network ENIGMA, for support. 7

8 References [] F. Calogero. Nuovo Cimento B43 978) 77. [] S. Ruijsenaars, H. Schneider. Ann. Phys ) 370. [3] J. Arnlind, J. Hoppe. Eigenvalue-Dynamics off the Calogero-Moser system.... Proceedings of the June 004 Les Houches Summer School Applications of Random Matrices in Physics. [4] J. Gibbons, T. Hermsen. Physica D 984) 337. [5] S. Wojciechowski. Phys. Lett. A 985) 0. [6] F. Calogero. J. Math. Phys ) 66. [7] J. Arnlind, J. Hoppe. Lett. Math. Phys ). [8] F. Calogero, J.P. Francoise. J. Phys. A ). [9] Y. B. Suris. The problem of integrable discretization: Hamiltonian approach. Birkhäuser Verlag, 003. [0] F. Calogero. Classical many-body problems amenable to exact treatments. Lecture Notes in Physics, Berlin 00). [] M. Olshanetsky, A. Perelemov. Lett. Nuovo Cimento 6 976) 333. [] H. Braden, R. Sasaki. Prog. Theor. Phys )