THES-TP 26/94 DEM-NT 94-36 November 1994 arxiv:hep-th/9411218v1 29 Nov 1994 The symmetry algebra of the N-dimensional anisotropic quantum harmonic oscillator with rational ratios of frequencies and the Nilsson model Dennis Bonatsos 1 C. Daskaloyannis 2 P. Kolokotronis and D. Lenis European Centre for Theoretical Research in Nuclear Physics and Related Areas (ECT ) Villa Tambosi Strada delle Tabarelle 286 I-38050 Villazzano (Trento) Italy Institute of Nuclear Physics N.C.S.R. Demokritos GR-15310 Aghia Paraskevi Attiki Greece Department of Physics Aristotle University of Thessaloniki GR-54006 Thessaloniki Greece Abstract The symmetry algebra of the N-dimensional anisotropic quantum harmonic oscillator with rational ratios of frequencies is constructed by a method of general applicability to quantum superintegrable systems. The special case of the 3-dim oscillator is studied in more detail because of its relevance in the description of superdeformed nuclei and nuclear and atomic clusters. In this case the symmetry algebra turns out to be a nonlinear extension of the u(3) algebra. A generalized angular momentum operator useful for labeling the degenerate states is constructed clarifying the connection of the present formalism to the Nilsson model in nuclear physics. 1 e:mail bonat nuclear.ect.unitn.it bonat cyclades.nrcps.ariadne-t.gr 2 e:mail daskaloyanni@olymp.ccf.auth.gr 1
1. Introduction Quantum algebras [1 2] (also called quantum groups) are nonlinear generalizations of the usual Lie algebras to which they reduce in the limiting case in which the deformation parameters are set equal to unity. From the mathematical point of view they have the structure of Hopf algebras [3]. The interest for applications of quantum algebras in physics was triggered in 1989 by the introduction of the q-deformed harmonic oscillator [4 5 6] although such mathematical structures have already been in existence [7 8]. The q- deformed harmonic oscillator has been introduced as a tool for providing a boson realization for the quantum algebra su q (2). Since then several generalized deformed oscillators (see [9 10 11 12 13] and references therein) as well as generalized deformed versions of the su(2) algebra [14 15 16 17 18] have been introduced and used in several applications to physical problems. On the other hand the 3-dimensional anisotropic quantum harmonic oscillator is of interest in the description of many physical systems. Its single particle level spectrum is known [19 20] to characterize the basic structure of superdeformed and hyperdeformed nuclei [21 22]. Furthermore it has been recently connected [23 24] to the underlying geometrical structure in the Bloch Brink α-cluster model [25]. It is also of interest for the interpretation of the observed shell structure in atomic clusters [26] especially after the realization that large deformations can occur in such systems [27]. The 2- dimensional anisotropic quantum harmonic oscillator is also of interest since for example its single particle level spectrum characterizes the underlying symmetry of pancake nuclei [20]. The anisotropic harmonic oscillator in two [28 29 30 31 32 33] and three [34 35 36 37 38 39 40] dimensions has been the subject of several investigations both at the classical and at the quantum mechanical level. The special cases with frequency ratios 1:2 [41 42] and 1:3 [43] have also been considered. The anisotropic harmonic oscillator is a well known example of a superintegrable system both at the classical and at the quantum mechanical level [44 45]. Although at the classical level it is clear that for the N- dimensional anisotropic harmonic oscillator the su(n) or sp(2nr) algebras can be used for its description the situation at the quantum level is not equally clear [28]. In section 2 of this letter the symmetry algebra of the N-dimensional anisotropic quantum harmonic oscillator with rational ratios of frequencies is studied by a method of general applicability to superintegrable systems 2
[46] while in section 3 the case of the 3-dimensional oscillator is considered in more detail because of its importance for the description of nuclei and possibly of atomic clusters. In the case of the 3-dim oscillator the symmetry algebra turns out to be a nonlinear generalization of u(3). In addition a generalized angular momentum operator is constructed suitable for labelling the degenerate states in a way similar to that used in the framework of the Nilsson model [47]. 2. General case: the N-dimensional oscillator Let us consider the system described by the Hamiltonian: H 1 ( ) N p 2 k + x2 k (1) 2 k1 m 2 k where m i are natural numbers mutually prime ones i.e. their great common divisor is gcd(m i m j ) 1 for i j and i j 1...N. We define the creation and annihilation operators [28]: a k ( ) 1 xk 2 m k ip k ( ) a k 1 xk 2 m k + ip k ( ) { U k 1 p 2 2 k + x2 k 1 m 2 ak a } N 2 k H U k. k k1 (2) These operators satisfy the relations (the indices k have been omitted): a U ( U 1 m) a or [U ( a ) m ] ( a ) m a U ( U + m) 1 a or [U (a) m ] (a) m [ a a ] U 1 2m aa U + 1 2m a a 1 ( ) m m a (a) m F(m U) (3) (a) m ( a ) m F(m U + 1) where the function F(m x) is defined by: F(m x) m ( p1 x 2p 1 2m ). (4) 3
Using the above relations we can define the enveloping algebra C defined by the polynomial combinations of the generators { 1 H A k A k U k} and k 1...N 1 where: A k ( a ) mk k (a N ) m N A k (a k ) m ( k a mn N). (5) These operators correspond to a multidimensional generalization of eq. (2): [H A k ] 0 [H A k ] 0 [H U k] 0 k 1... N 1. (6) The following relations are satisfied for k l and k l 1...N 1: while [U k A l ] [ U k A l] [Ak A l ] [ A k A l] 0 (7) U k A k A k (U k + 1) U k A k A k (U k 1) A k A k F ( ) ( m k U k F mn H N 1 U l + 1 ) A k A k F ( m k U k + 1 ) F ( m N H N 1 ) U l. One additional relation for k l can be derived: F ( m N H N 1 ) U l A k A l F ( N 1 m N H (8) U l + 1 ) A l A k. (9) The algebra defined by the above relations accepts a Fock space representation. The elements of the basis E p 1...p N 1 are characterized by the eigenvalues of the N commuting elements of the algebra H and U k with k 1...N 1. The elements A k and A k are the corresponding ladder operators of the algebra. The following relations hold: H E p 1...p N 1 E E p 1...p N 1 U k E p 1...p N 1 p k E p 1...p N 1 A k E p 1...p k...p N 1 F ( ) ( m k p k F mn E N 1 p l + 1 ) E p 1...p k 1... p N 1 A k E p 1...p k...p N 1 F ( m k p k + 1 ) F ( m N E N 1 ) p l E p1...p k + 1...p N 1. (10) 4
Let p min k be the minimum value of p k such that E A k p1...p min k...p N 1 0. (11) From eq. (10) we find that we must have: F ( m k p min ) k 0. (12) Then p k is one of the roots of the function F defined by eq. (4). The general form of the roots is: p min k 2q k 1 q k 1... m k. (13) 2m k Each root is characterized by a number q k. The numbers q k also characterize the representations of the algebra as we shall see. The elements of the Fock space can be generated by successive applications of the ladder operators A k on the minimum weight element E p min 1...p min E [0] N 1 E 0...0 q 1...q k...q N The elements of the basis of the Fock space are given by: E [p min k + n k ] E [n] 1 ( N 1 ) ( ) A nk E [0] C [n] k k1. (14) (15) where [n] (n 1 n 2 n N 1 ) and (q 1 q 2... q N ) while C [n] are normalization coefficients. The generators of the algebra acting on the base of the Fock space give: H E [n] E [n] U k A k E E [n] E [n] ( n k + p min k F ( m k n k + p min k E n 1...n k + 1... n N 1 q 1...q k...q N ) E [n] k 1...N 1 + 1 ) F ( m N E N 1 5 ( nl + p min ) l ) (16) (17)
E [n] A k F ( m k n k + p min ) ( k F mn E N 1 E n 1... n k 1... n N 1. q 1... q k...q N ( nl + p min ) l + 1 ) (18) The existence of a finite dimensional representation implies that after Σ successive applications of the ladder operators A on the minimum weight element one gets zero so that the following condition is satisfied: Therefore: where p min N F(m N E Σ E Σ N 1 N 1 p min l p min l ) 0. p min N is the root of equation F(m N p min N ) 0. Then E Σ + N k1 2q k 1 2m k. (19) In the case of finite dimensional representations only the energies given by eq. (19) are permitted and the elements of the Fock space can be described by using Σ instead of E. The action of the generators on the Fock space is described by the following relations: H Σ [n] U k Σ [n] A k Σ [n] ( Σ + N ) 2q k 1 2m k1 k ) 2m k ( n k + 2q k 1 F ( m k n k + 2q k 1 2m k Σ [n] Σ [n] Σ n 1...n k + 1...n N 1 q 1...q k...q N + 1 ) F ( m N Σ N 1 6 k 1...N 1 n l + 2q N 1 2m N ) (20) (21)
Σ [n] A k F ( ) ( m k n k + 2q k 1 2m k F mn Σ N 1 Σ n 1...n k 1...n N 1. q 1...q k...q N The dimension of the representation is given by ( ) Σ + N 1 d Σ n l + 2q N 1 2m N + 1 ) (Σ + 1)(Σ + 2) (Σ + N 1). (N 1)! (22) It is clear that to each value of Σ correspond m 1 m 2...m N energy eigenvalues each eigenvalue having degeneracy d. Using eqs (3) we can prove that the algebra generated by the generators a l a l N l m l U l 1/2 is an oscillator algebra with structure function [9 12] i.e. with Φ l (x) x/m l Φ l (N l ) U l 1 2m l. This oscillator algebra is characterized by the commutation relations: [ Nl a ] l a l [N l a l ] a l a l a l Φ l (N l ) a l a l Φ l (N l + 1). (23) There are in total N different oscillators of this type uncoupled to each other. The Fock space corresponding to these oscillators defines an infinite dimensional representation of the algebra defined eqs (6-9). In order to see the connection of the present basis to the usual Cartesian basis one can use for the latter the symbol [r] (r 1 r 2...r N ). One then has a l Φ(r [r] l + 1) r 1...r l + 1... r N a l [r] Φ(r l ) r 1... r l 1... r N N l [r] r l [r]. The connection between the above basis and the basis defined by eqs (20-22) 7
is given by: Σ [n] [r] n k [r k /m k ] k 1...N 1 r l n l m l + mod (r l m l ) q l 1... N l mod (r l m l ) + 1 Σ N [r l /m l ] (24) where [x] means the integer part of the number x. Using the correspondence between the present basis and the usual Cartesian basis given in eq. (24) the action of the operators a k on the present basis can be calculated for k 1...N 1: a k Σ [n] n k + q k /m k Σ [n ] [q ] where n l n l q l q l for l k n k n k + [q k /m k ] Σ Σ + [q k /m k ] q k mod (q k m k ) + 1 while for the operator a N one has a N Σ [n] Σ N 1 k1 n k + q N /m N Σ [n ] [q ] where n k n k q k q k for k 1...N 1 Σ Σ + [q N /m N ] q N mod (q N m N ) + 1. Similarly for the operators a k one can find for k 1... N 1: a k Σ [n] n k + (q k 1)/m k Σ [n ] [q ] 8
where while for the operator a N one has Σ [n] a N Σ where n l n l q l q l for l k n k n k + [(q k 2)/m k ] Σ Σ + [(q k 2)/m k ] q k mod (q k 2 m k ) + 1 N 1 k1 n k + (q N 1)/m N Σ [n ] [q ] n k n k q k q k for k 1...N 1 Σ Σ + [(q N 2)/m N ] q N mod (q N 2 m N ) + 1. An important difference between the operators a k a k used here and the operators A (s) i A (s) i used in refs [34 36 37 40] has to be pointed out. The operators used here satisfy the relations [a k a l ] 1 m k δ kl i.e. they represent oscillators completely decoupled from each other while the operators of refs [34 36 37 40] satisfy the relations [A (s) i A (s) j ] δij [A (s) i A (s ) j ] δij i.e. they represent oscillators not completely decoupled. Notice that (s) of refs [34 36 37] ({λ} of ref. [40]) is analogous to the used in the present work. 3. The 3-dimensional oscillator and relation to the Nilsson model In this section the 3-dim case will be studied in more detail because of its relevance for the description of superdeformed nuclei and of nuclear and atomic clusters. The 3-dim anisotropic oscillator is the basic ingredient of the Nilsson model [47] which in addition contains a term proportional to the square of the angular momentum operator as well as a spin-orbit coupling term the relevant Hamiltonian being H Nilsson H osc 2k L S kν L 2 9
where k ν are constants. The spin-orbit term is not needed in the case of atomic clusters while in the case of nuclei it can be effectively removed through a unitary transformation both in the case of the spherical Nilsson model [48 49 50] and of the axially symmetric one [51 52]. An alternative way to effectively remove the spin-orbit term in the spherical Nilsson model is the q-deformation of the relevant algebra [53]. It should also be noticed that the spherical Nilsson Hamiltonian is known to possess an osp(1 2) supersymmetry [54]. In the case of the 3-dim oscillator the relevant operators of eq. (2) form a nonlinear generalization of the algebra u(3) the q-deformed version of which can be found in [55 56 57]. As we have already seen to each Σ value correspond m 1 m 2 m 3 energy eigenvalues each eigenvalue having degeneracy (Σ + 1)(Σ + 2)/2. In order to distinguish the degenerate eigenvalues we are going to introduce some generalized angular momentum operators L i (i 1 2 3) defined by: One can prove that L k iǫ ijk ( (ai ) m i ( a mj j) ( a ) mi i (a j ) m ) j. (25) L 1 i ( A 2 A 2) L2 i ( A 1 A 1 ). The following commutation relations can be verified: [L i L j ] iǫ ijk (F(m k U k + 1) F(m k U k ))L k. (26) It is worth noticing that in the case of m 1 m 2 m 3 1 one has F(1 x) x 1/2 so that the above equation gives the usual angular momentum commutation relations. The operators defined in eq. (25) commute with the oscillator hamiltonian H and therefore conserve the number Σ which characterizes the dimension of the representation. Also these operators do not change the numbers q 1 q 2 q 3 as we can see from eqs (21-22). The eigenvalues of these operators can be calculated using Hermite function techniques. Let us consider in particular the generalized angular momentum projection: L 3 i ( (a 1 ) m ( 1 a m2 2) ( a ) m1 1 (a 2 ) m ) 2. (27) 10
This acts on the basis as follows Σ [n] L 3 i ( F ( ) ( m 1 n 1 + 2q 1 1 2m 1 F m2 n 2 + 2q 2 1 2m 2 + 1 ) Σ n 1 1 n 2 + 1 F ( m 1 n 1 + 2q 1 1 2m 1 + 1 ) F ( ) m 2 n 2 + 2q 2 1 Σ n 1 + 1 n 2 1 ). 2m 2 (28) This operator conserves the quantum number j n 1 + n 2 j 0 1 2 2 1 3 2... In addition one can introduce the quantum number m n 1 n 2. 2 One can then replace the quantum numbers n 1 n 2 by the quantum numbers j µ where µ is the eigenvalue of the L 3 operator. The new representation basis one can label as L 3 Σ j µ µ Σ j µ. (29) This basis is connected to the basis of the previous section as follows Σ j µ j m j c[j m µ] [j + m] 1![j m] 2! where ( [0] k! 1 [n] k! [n] k [n 1] k! [n] k F Σ j + m j m m k n + 2q k 1 2m k (30) and the coefficients c[j m µ] in eq. (30) satisfy the recurrence relation: µc[j m µ] i ([j m] 2 c[j m + 1 µ] [j + m] 1 c[j m 1 µ]). (31) 11 )
These relations can be satisfied only for special values of the parameter µ corresponding to the eigenvalues of L 3. It is worth noticing that in the case of m 1 m 2 which corresponds to axially symmetric oscillators the possible values turn out to be µ 2j 2(j 1)...2(j 1) 2j. In nuclear physics the quantum numbers n n 1 + n 2 and Λ ±n ±(n 2)... ±1 or 0 are used [58]. From the above definitions it is clear that j n /2 and µ Λ. Therefore in the case of m 1 m 2 which includes axially symmetric prolate nuclei with m 1 : m 2 : m 3 1 : 1 : m as well as axially symmetric oblate nuclei with m 1 : m 2 : m 3 m : m : 1 the correspondence between the present scheme and the Nilsson model is clear. In the case of axially symmetric prolate nuclei (m 1 m 2 1) one can easily check that the expression L 2 L 2 1 + L2 2 + (F(m 3 U 3 + 1) F(m 3 U 3 ))L 2 3 satisfies the commutation relations [L 2 L i ] 0 i 1 2 3. 4. Discussion The symmetry algebra of the N-dim anisotropic quantum harmonic oscillator with rational ratios of frequencies has been constructed by a method [46] of general applicability in constructing finite-dimensional representations of quantum superinegrable systems. The case of the 3-dim oscillator has been considered in detail because of its relevance to the single particle level spectrum of superdeformed and hyperdeformed nuclei [19 20] to the underlying geometrical structure of the Bloch-Brink α-cluster model [23 24] and possibly to the shell structure of atomic clusters at large deformations [26 27]. The symmetry algebra in this case is a nonlinear generalization of the u(3) algebra. For labeling the degenerate states generalized angular momentum operators are introduced clarifying the connection of the present approach to the Nilsson model. In the case of the 2-dim oscillator with ratio of frequencies 2 : 1 (m 1 1 m 2 2) it has been shown [59] that the relevant nonlinear generalized u(2) algebra can be identified as the finite W algebra W (2) 3 [60 61]. In the case of the 3-dim axially symmetric oblate oscillator with frequency ratio 1 : 2 (which corresponds to the case m 1 m 2 2 m 3 1) the relevant symmetry is related to O(4) [40 62]. The search for further symmetries related to specific 12
frequency ratios hidden in the general nonlinear algebraic framework given in this work is an interesting problem. One of the authors (DB) has been supported by the EU under contract ERBCHBGCT930467. This project has also been partially supported by the Greek Secretariat of Research and Technology under contract PENED 340/91. 13
References [1] Drinfeld V G 1986 Quantum groups Proceedings of the International Congress of Mathematicians ed A M Gleason (Providence RI: American Mathematical Society) p 798 [2] Jimbo M 1986 Lett. Math. Phys. 11 247 [3] Abe E 1977 Hopf Algebras (Cambridge: Cambridge U. Press) [4] Biedenharn L C 1989 J. Phys. A: Math. Gen. 22 L873 [5] Macfarlane A J 1989 J. Phys. A: Math. Gen. 22 4581 [6] Sun C P and Fu H C 1989 J. Phys. A: Math. Gen. 22 L983 [7] Arik M and Coon D D 1976 J. Math. Phys. 17 524 [8] Kuryshkin V V 1980 Annales de la Fondation Louis de Broglie 5 111 [9] Daskaloyannis C 1991 J. Phys. A: Math. Gen. 24 L789 [10] Arik M Demircan E Turgut T Ekinci L and Mungan M 1992 Z. Phys. C 55 89 [11] Brzeziński T Egusquiza I L and Macfarlane A J 1993 Phys. Lett. 311B 202 [12] Bonatsos D and Daskaloyannis C 1993 Phys. Lett. 307B 100 [13] Meljanac S Milekovic M and Pallua S 1994 Phys. Lett. 328B 55 [14] Bonatsos D Daskaloyannis C and Kolokotronis P 1993 J. Phys. A: Math. Gen. 26 L871 [15] Delbecq C and Quesne C 1993 J. Phys. A: Math. Gen. 26 L127; 1993 Phys. Lett. 300B 227; 1993 Mod. Phys. Lett. A 8 961 [16] Ludu A and Gupta R K 1993 J. Math. Phys. 34 5367 [17] Sun C P and Li W 1993 Commun. Theor. Phys. 19 191 [18] Pan F 1994 J. Math. Phys. 35 5065 14
[19] Mottelson B 1991 Nucl. Phys. A 522 1c [20] Rae W D M 1988 Int. J. Mod. Phys. A 3 1343 [21] Nolan P J and Twin P J 1988 Ann. Rev. Nucl. Part. Sci. 38 533 [22] Janssens R V F and Khoo T L 1991 Ann. Rev. Nucl. Part. Sci. 41 321 [23] Rae W D M and Zhang J 1994 Mod. Phys. Lett. A 9 599 [24] Zhang J Rae W D M and Merchant A C 1994 Nucl. Phys. A 575 61 [25] Brink D M 1966 Proc. Int. School of Physics Enrico Fermi Course XXXVI Varenna 1966 ed C Bloch (New York: Academic Press) p 247 [26] Martin T P Bergmann T Göhlich H and Lange T 1991 Z. Phys. D 19 25 [27] Bulgac A and Lewenkopf C 1993 Phys. Rev. Lett. 71 4130 [28] Jauch J M and Hill E L 1940 Phys. Rev. 57 641 [29] Demkov Yu N 1963 Soviet Phys. JETP 17 1349 [30] Contopoulos G 1960 Z. Astrophys. 49 273; 1963 Astrophys. J. 138 1297 [31] Cisneros A and McIntosh H V 1970 J. Math. Phys. 11 870 [32] Ghosh A Kundu A and Mitra P 1992 Saha Institute preprint SINP/TNP/92-5 [33] Castaños O and López-Peña R 1992 J. Phys. A: Math. Gen. 25 6685 [34] Duimio F and Zambotti G 1966 Nuovo Cimento 43 1203 [35] Maiella G 1967 Nuovo Cimento 52 1004 [36] Vendramin I 1968 Nuovo Cimento 54 190 [37] Maiella G and Vilasi G 1969 Lettere Nuovo Cimento 1 57 [38] Rosensteel G and Draayer J P 1989 J. Phys. A: Math. Gen. 22 1323 15
[39] Bhaumik D Chatterjee A and Dutta-Roy B 1994 J. Phys. A: Math. Gen. 27 1401 [40] Nazarewicz W and Dobaczewski J 1992 Phys. Rev. Lett. 68 154 [41] Holt C R 1982 J. Math. Phys. 23 1037 [42] Boyer C P and Wolf K B 1975 J. Math. Phys. 16 2215 [43] Fokas A S and Lagerstrom P A 1980 J. Math. Anal. Appl. 74 325 [44] Hietarinta J 1987 Phys. Rep. 147 87 [45] Evans N W 1990 Phys. Rev. A 41 5666 [46] Bonatsos D Daskaloyannis C and Kokkotas K 1993 it Phys. Rev. A 48 R3407 [47] Nilsson S G 1955 K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 29 No 16 [48] Castaños O Moshinsky M and Quesne C 1992 Group Theory and Special Symmetries in Nuclear Physics (Ann Arbor 1991) ed J P Draayer and J W Jänecke (Singapore: World Scientific) p 80 [49] Castaños O Moshinsky M and Quesne C 1992 Phys. Lett. 277B 238 [50] DraayerJ P Bahri C and Moszkowski S 1992 Rev. Mex. Fís. 38 Supl 2 26 [51] Moshinsky M and Del Sol Mesa A 1992 Rev. Mex. Fís. 38 Supl 2 146 [52] Castaños O V. Valázquez A. Hess P O and Hirsch J G 1994 Phys. Lett. 321B 303 [53] Del Sol Mesa A Loyola G Moshinsky M and Velázquez V 1993 J. Phys. A: Math. Gen. 26 1147 [54] Balantekin A B Castaños O and Moshinsky M 1992 Phys. Lett. 284B 1 [55] Smirnov Yu F Tolstoy V N and Kharitonov Yu I 1991 Yad. Fiz. 54 721 [1991 Sov. J. Nucl. Phys. 54 437] 16
[56] Smirnov Yu F and Kharitonov Yu I 1993 Yad. Fiz. 56 263 [1993 Phys. At. Nucl. 56 1143] [57] Bonatsos D Faessler A Raychev P P Roussev R P and Smirnov Yu F 1992 J. Phys. A: Math. Gen. 25 L267 [58] Bohr A and Mottelson B R 1975 Nuclear Structure Vol II (Reading: Benjamin) [59] Bonatsos D Daskaloyannis C Kolokotronis P and Lenis D 1994 NCSR Demokritos preprint DEM-NT-94-09 U. Thessaloniki preprint THES- TP 25/94 [60] Tjin T 1992 Phys. Lett. 292B 60 [61] de Boer J and Tjin T 1993 Commun. Math. Phys. 158 485 [62] Ravenhall D G Sharp R T and Pardee W J 1967 Phys. Rev. 164 1950 17