arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

Transcription

1 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/ v [nin.cd] 30 Jan Institut für Theoretische Physik, Universität Stuttgart, Stuttgart, Germany Abstract We extend a theory of first order corrections to Gutzwier s trace formua for systems with a smooth potentia to systems with discrete symmetries and, as an exampe, appy the method to the two-dimensiona hydrogen atom in a uniform magnetic fied. We expoit the C 4v -symmetry of the system in the cacuation of the correction terms. The numerica resuts for the semicassica vaues wi be compared with vaues extracted from exact quantum mechanica cacuations. The comparison shows an exceent agreement and demonstrates the power of the expansion method. Key words: semicassica theories, Gutzwier s trace formua, expansions, discrete symmetries, diamagnetic hydrogen atom PACS: Sq, Mt 1 Introduction Semicassica theories have become very important for a deeper understanding of quantum systems, and Gutzwier s trace formua [1] has become a powerfu too for cassicay chaotic systems. It provides a semicassica approximation of the quantum eve density in terms of cassica periodic orbits. In a systematic expansion of the eve density in powers of it can be considered as the eading order. Higher orders of this asymptotic expansion have been deveoped in severa studies [,3,4], but for a ong time were ony tested for biiard systems, i.e., systems with hard was instead of smooth potentias. By extending an expansion which was derived by Gaspard et a. [3,4], Grémaud [5] deveoped corrections to Gutzwier s trace formua for quantum systems with a smooth potentia. He presented a first-order correction term to Gutzwier s trace Preprint submitted to Esevier Science 8 February 009

2 formua and obtained numerica resuts for the diamagnetic hydrogen atom. Comparisons with vaues extracted from exact quantum cacuations showed that the semicassica resuts can be computed with very high accuracy. However, in the theory presented in [5] some important topics have not yet been considered. For exampe, the corrections cannot be cacuated for cassica orbits which have a turning point (i.e., the veocity vanishes competey at this point). Furthermore, discrete and continuous symmetries have not been incuded. If there are discrete symmetries the eigenstates of the quantum system spit up into severa subspaces. In these subspaces, cassica orbits which are not periodic without a symmetry transformation contribute to the eve density. Preiminary resuts for the diamagnetic hydrogen atom, which is an exampe system with discrete symmetries, were pubished in [6] without expaining which modifications in the numerica cacuations the introduction of discrete symmetries entais. In particuar, an anaysis of the differentia equations which have to be soved to obtain the correction terms is not presented. In this paper we want to give a transparent derivation of the modifications necessary for the cacuation of the correction terms in systems with discrete symmetries. We wi show in detai how such symmetries of the Hamitonian have to be taken into account. It wi be expained in which equations it is necessary to introduce symmetry operations and which transformations ead to the correct boundary conditions of the cassica Greens function, which is an essentia part of the correction terms. Furthermore, we want to give a deeper understanding of the significance of the symmetry of the system for the quantum spectrum and the connection to periodic orbits. We wi appy the method to the diamagnetic hydrogen atom. The threedimensiona hydrogen atom in a uniform magnetic fied has a continuous symmetry, namey the rotationa invariance around the magnetic fied axis. Continuous symmetries have a substantia infuence on the correction terms. Actuay one cannot obtain the correct resuts for the first-order correction for the three-dimensiona hydrogen atom with the formuas presented because the rotationa invariance eads to additiona contributions to the first-order corrections. Therefore, we wi consider the hydrogen atom as a pure twodimensiona system as it was done before in [5]. If one uses semiparaboic coordinates, the potentia of the two-dimensiona diamagnetic hydrogen atom exhibits a discrete C 4v -symmetry. The corrections wi be cacuated for seected periodic orbits taking into account that discrete symmetry. The semicassica resuts wi be compared with the anaysis of exact quantum cacuations. The agreement between the resuts of both methods turns out to be very good. The outine of the paper is as foows. In section we wi first give a summary

3 of the derivation of the corrections without taking into account discrete symmetries. In section 3, we wi introduce the hydrogen atom in a uniform magnetic fied with a aspects reevant for the cacuation of the correction terms. We wi aso discuss the infuence of the rotationa invariance of the three-dimensiona hydrogen atom. Then we wi extend the corrections to discrete symmetries, cacuate resuts for the two-dimensiona hydrogen atom and compare them with exact quantum mechanica cacuations in section 4. Semicassica approximation of the quantum eve density As mentioned in the introduction the corrections to Gutzwier s trace formua for systems with a smooth potentia were derived by Grémaud [5] based on earier work by Gaspard et a. [3,4]. In this section, for the reader s convenience we wi briefy summarize a major steps of this derivation by foowing the way described in [5] in order to be abe to incude discrete symmetries in section 4..1 The trace of the propagator and the first part of the correction The starting point for the derivation of the semicassica eve density which is used in this theory is a discrete version of Feynman s path integra representation of the propagator K( q, q 0, T) = 1 d q (πi t) Nf 1 d q d q N 1 [ ( ) ] i N 1 qn+1 q n exp L, q n t + O ( t), (1) t n=0 for a system with f degrees of freedom and a time independent Hamitonian: H = p + V ( q). () K( q, q 0, T) represents the propagation of a partice going from q 0 to q N = q in time T with t = T/N. L( q, q) denotes the cassica Lagrangian of the system. For the cacuation of the eve density one has to know the trace of the propagator, which can be obtained by setting q 0 = q in (1) and by integrating over 3

4 this variabe: K(T) = 1 K( q 0, q 0, T) d q 0 = d q (πi t) Nf 0 d q 1 d q d q N 1 [ ( ) i N 1 qn+1 q n exp L, q n t + O ( t) t n=0 ]. (3) After expanding the exponentia function and incuding a terms which produce contributions to the eading order in and to the first-order correction, the integra is evauated as described in [4]. The resut in the imit t 0 is given by: K(T) = { K (0) (T) 1 + i Tp T p 0 dt 0 C 1 (T, t 0 ) + O ( )}. (4) In this formua the sum runs over a periodic orbits of the cassica system described by the Hamitonian () with period T. T p represents the time period of the primitive periodic orbit, i.e., if the orbit consists of mutipe repetitions of a periodic orbit, one considers ony the basic traversa in this case. K (0) (T) stands for the eading order of the semicassica approximation. It is identica with Gutzwier s [1] resut, and is given by: K (0) (T) = 1 π T p E T det (m (T) 1) exp [ i W (c) (T) iπ µ + iπ 4 sign ( ET) ]. (5) In this expression one can find the (f ) (f )-dimensiona stabiity matrix m, which represents the stabiity properties of the periodic orbit. The determinant det (m (T) 1) is reated to the f f-dimensiona monodromy matrix. This reation is shown in section.4.. A further term is the cassica action T W (c) (T) = L( q, q) dt, (6) 0 µ is the Masov index, and E T is the derivative of the time period with respect to the energy E. In equation (4) C 1 is the abbreviation of the firstorder correction in. As one can see in (4) one has to use the average of C 1 (T, t 0 ), where t 0 parametrizes the periodic orbit. The integra over t 0 runs from 0 to T p, i.e., one has to integrate over the whoe primitive periodic orbit. The expicit expression of C 1 (T, t 0 ) needs engthy cacuations, which are 4

5 discussed in [3,4,5] in detai. The resut reads: C 1 (T, t 0 ) = 1 8 where T T 0 T dt V,i1 i i 3 i 4 ( q (c) (t) ) G i1 i (t, t; t 0 )G i3 i 4 (t, t; t 0 ) + 1 ( dt dt (c) V,i1 i 4 i 3 q (t) ) ( (c) V,j1 j j 3 q (t ) ) ( 3G i1 i (t, t; t 0 ) 0 0 G i3 j 1 (t, t ; t 0 )G j j 3 (t, t ; t 0 ) + G i1 j 1 (t, t ; t 0 )G i j (t, t ; t 0 )G i3 j 3 (t, t ; t 0 ) ) + V,j(t 0 ) q (c) (t 0 ) T 0 dt V,i1 i i 3 ( q (c) (t) ) G ji1 (0, t; t 0 )G i i 3 (t, t; t 0 ), (7) ( ( (c) V,i1...i n q (t) ) = n (c) V q (t) ) (8) q i1... q in are derivatives of the potentia evauated at the point q (c) (t) on the cassica orbit. G ij are the components of the cassica Green s function, which is a soution of the inearized equation of motion: ( 1 d dt V ( q (c) (t) )) G(t, t ) = 1δ(t t ). (9) q q Due to the factor 1/ q (c) in equation (7), the correction term is singuar for vanishing veocities, and because of the integra over t 0 in equation (4), the correction term diverges for orbits with a turning point.. The cassica Green s function for the trace of the propagator The cassica Green s function is an essentia part of the first-order correction to the trace of the propagator. In order to be abe to ook at the symmetry transformations of C 1 in equation (7), one has to know its structure. In this section we discuss a important parts of this structure. From the derivation of the correction term C 1, it foows that the cassica Green s function has to fufi the boundary conditions [5] G(0, t ) = G(T, t ), P t0 G(0, t ) = P t0 G(T, t ) = 0, Q t0 G(0, t ) = Q t0 G(T, t ), (10) with the projection operator P t0 aong the direction of the cassica orbit at time t 0, which has the form q(t 0 ) q(t 0 ) (P t0 ) ij = q(t 0 ) = q i(t 0 ) q j (t 0 ) q(t 0 ) ; (11) ij 5

6 furthermore, Q t0 = 1 P t0, with the f-dimensiona unity matrix 1. If we use the notation G (t, t ) = G(t, t ) for 0 t t, G + (t, t ) = G(t, t ) for t (1) t T, it is possibe to write the cassica Green s function as a product [5] G ±(t, t ) = M(t) A ±(t ) (13) G ± (t, t ) B ± (t ) with the f f monodromy matrix M(t) and the four f f matrices A ± and B ±. The monodromy matrix is a sympectic matrix, which can be obtained by soving the inearized Hamitonian equations of motion Here, Σ is the matrix Ṁ(t, T) = Σ H X X M(t, T). (14) Σ = and 1 is the f f unity matrix. For a Hamitonian of the form () it has the structure: M(t) = J (t) J 1 (t), (15) J (t) J 1 (t) where J 1 and J are f f matrices. The boundary conditions at time t = t are G (t, t ) = G +(t, t ) + 0. (16) G (t, t ) G + (t, t ) 1 Expoiting the formuation (13) of the cassica Green s function and the sympecticity of the monodromy matrix M(t), one can formuate the boundary condition (16) as foows: A +(t ) = A (t ) M(t ) 1 0 = A (t ) J 1(t ) T. (17) B + (t ) B (t ) 1 B (t ) J (t ) T If this expression is combined with the condition 1 0 A (t ) = 1 0 M(T) A +(t ), (18) 0 Q t0 B (t ) 0 Q t0 B + (t ) 6

7 which is identica with the first and third condition of (10), it eads to the matrix equation: Q t0 [M(T) 1] A (t ) = B (t ) Q t0 M(T) J 1(t ) T J (t ) T. (19) This equation has to be soved in order to obtain the matrices A and B. When A and B are known one has the soution for the cassica Green s function (13). Unfortunatey, the matrix 1 0 [M(T) 1] 0 Q t0 is singuar but, as was shown in [5], there is a soution for cassica trajectories which have no turning point (i.e. the veocity q(t) never vanishes competey aong the orbit), viz. A (t ; t 0 ) 1 q(t 0 ) q(t 0 ) 0 = X = X 0 + B (t ; t 0 ) q(t 0 ) X 0, (0) ṗ(t 0 ) q(t 0 ) 0 where X 0 represents a particuar soution of equation (19), which can be determined by a singuar vaue decomposition (see e.g. [7]). The second condition of equation (10) has been used to derive this resut. In summary, it is possibe to obtain the cassica Green s function by soving the inearized equation of motion (14) of the monodromy matrix and by cacuating A and B from equation (0). Then the correction term C 1 (T, t 0 ) foows from equation (7). In equation (0), the factor 1/ q (c) appears in the construction of C 1 (T, t 0 ) a second time. The reason for this ies in the appication of the projection operator P t0 in the second condition of (10). The integras and doube integras in equation (7) can be transformed into a set of ordinary differentia equations, which can be computed effectivey. The transformation of the integras is discussed in [5] in detai..3 The trace of the quantum Green s function and the second part of the corrections The next step on the way to the eve density is the trace of the quantum Green s function, which can be obtained from the trace of the propagator via 7

8 the semi-sided Fourier transform G(E) = G (E) = 1 [ ] i dt exp i 0 ET K (T). (1) One has to incude a contributions to the eading order as we as the firstorder correction in the semicassica approximation of the integra. The methods which are used to evauate this integra are discussed in [4]. We do not want to repeat the cacuation but present the resuts obtained in [4] and [5], which can be summarized as foows: T 0 p G (E) = 1 i det (m (T 0 ) 1) exp [ i S(c) (T 0 ) iπ µ { 1 + i ( C 1 (T 0 ) + C T E 1 (T 0 ) ) + O ( )}, () where the eading order, beonging to the 1 in the cury brackets, is known from Gutzwier s trace formua [1]. The reduced action S (c) (T 0 ) = W (c) (T 0 ) + ET 0 (3) is used in this expression, and µ represents the Masov index. The first-order correction consists of two terms. The first, C 1 (T 0 ) = 1 T0p C 1 (T 0, t 0 ) dt 0, (4) T 0p 0 is the same as in equation (7) but now with fixed period T 0, and the second is the contribution from integra (1) to the correction, and is given by: C T E 1 (T 0 ) = 1 W () (T 0 ) 1 ( C (1) 0 (T 0) + C () 0 (T 0) ) W (3) (T 0 )C (1) 0 (T 0) 1 W () (T 0 ) 8 W (4) (T 0 ) W () ] (T 0 ) W (3) (T 0 ) (5) W () (T 0 ) 3. W (n) and C (n) 0 are the derivatives W (n) (T) = n W (c) (T), C (n) T n 0 (T) = nc (c) 0 (T) T n of the action W (c) (T) and of the ogarithm of the ampitude of the trace of the propagator respectivey. T p C 0 (T) = n, (6) E T det (m (T) 1) 8

9 .4 Contributions to the correction terms of the trace of the semicassica Green s function In this section we want to summarize the contributions to the second correction term Ct T E which are important for the understanding of the symmetry properties of this term. As we wi see, a parts of equation (5) can be obtained as numerica soutions of differentia equations of some new coordinates..4.1 Derivatives of W (c) (T) Because of the reation W (1) (T 0 ) = (T) = E(T 0 ) = H( X(t, T T 0 )) (7) T=T0 W (c) a higher derivatives of the action W (c) (T 0 ) with respect to the time period of the orbit can be expressed as derivatives of the cassica Hamitonian, where X(t, T 0 ) = ( q(t, T 0 ), p(t, T 0 )) represents the phase space vector of the orbit. Introducing the derivatives X (n) (t, T 0 ) = n X(t, T) T n which appear as coefficients of the Tayor expansion X(t, T) = X(t, T 0 + δt) = n=0, (8) T=T0 1 n! X (n) (t, T 0 )(δt) n (9) of the orbit X(t, T) around a reference orbit X(t, T 0 ) with period T 0, the higher derivatives of W (c) (T 0 ) are: W () (T 0 ) = H,i ( X(t, T 0 ))X (1) i (t, T 0 ), W (3) (T 0 ) = H,ij ( X(t, T 0 ))X (1) i (t, T 0 )X (1) j (t, T 0 ) H,i ( X(t, T 0 ))X () i (t, T 0 ), W (4) (T 0 ) = H,ijk ( X(t, T 0 ))X (1) i (t, T 0 )X (1) j (t, T 0 )X (1) k (t, T 0) 3H,ij ( X(t, T 0 ))X (1) i (t, T 0 )X () j (t, T 0 ) H,i ( X(t, T 0 ))X (3) i (t, T 0 ). (30) From Hamiton s equations of motion for the phase space coordinates Ẋ(t, T 0 ) = Σ H( X(t, T 0 )) X, 9

10 one can infer the equations which govern the motion of the derivatives X (n) (t, T 0 ): Ẋ (1) i (t, T 0 ) = Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0), Ẋ () i (t, T 0 ) = Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0)X (1) (t, T 0 ) + Σ ij H,jk ( X(t, T 0 ))X () k (t, T 0), Ẋ (3) i (t, T 0 ) = Σ ij H,jkm ( X(t, T 0 ))X (1) k (t, T 0)X (1) (t, T 0 )X m (1) (t, T 0) + 3Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0)X () (t, T 0 ) + Σ ij H,jk ( X(t, T 0 ))X (3) k (t, T 0). (31) These equations are inhomogeneous differentia equations. Their soutions are of the type X (1) (t, T 0 ) = M(t, T 0 ) X (1) (0, T 0 ), X () (t, T 0 ) = M(t, T 0 ) X () (0, T 0 ) + F () (t, T 0 ), X (3) (t, T 0 ) = M(t, T 0 ) X (3) (0, T 0 ) + F (3) (t, T 0 ), (3) where M(t, T 0 ) is the monodromy matrix, which was introduced in section.. Together with the condition X(t = 0, T 0 ) = X(t = T 0, T 0 ) (33) for periodic orbits, the soutions (3) ead to the foowing equations for the initia vaues: [1 M(T 0, T 0 )] X (1) (0, T 0 ) = Ẋ(T 0, T 0 ), [1 M(T 0, T 0 )] X () (0, T 0 ) = Ẋ (1) (T 0, T 0 ) + Ẍ(T 0, T 0 ) + F () (T 0, T 0 ), [1 M(T 0, T 0 )] X (3) (0, T 0 ) = 3 Ẋ () (T 0, T 0 ) + 3 Ẍ (1) (T 0, T 0 ) +... X(T 0, T 0 ) + F (3) (T 0, T 0 ). (34) The matrix [1 M(T 0, T 0 )] is singuar. A particuar soution can be found via a singuar vaue decomposition. For a genera unstabe orbit the kerne is one-dimensiona, i.e., the soution space is in genera one-dimensiona (see [5]) but, as aso shown in [5], W () (T 0 ) can directy be computed if a particuar soution X (1) 0 (0, T 0 ) from (34) is known: W () (T 0 ) = H( X(T 0, T 0 )) X (1) 0 (0, T 0 ). (35) Simiar reations can be found for W (3) (T 0 ) and W (4) (T 0 ), in which aso ony a particuar soution of equation (34) is required. 10

11 .4. Derivatives of C 0 (T) If one ooks at the definition of C 0 (T) in equation (6), one can easiy convince onesef that C (1) 0 (T 0) = W (3) (T 0 ) T 0 W () (T 0 ) 1 d n det (m(t) 1), dt T=T0 C () 0 (T 0) = W (4) (T 0 ) T0 W () (T 0 ) 1 W (3) (T 0 ) W () (T 0 ) 1 d n det (m(t) 1). dt T=T0 (36) The derivatives of the action are aready known from section.4.1. In the next step one has to find an expression for the derivatives of n det (m(t) 1). Foowing again the method presented in [5] we can introduce a new matrix N(T) whose determinant is identica with det (m(t) 1). It is given by where and N(T) = M(T) ( 1 P (T) P (T) ), (37) P (T) = ê ê = Ẋ(0, T) Ẋ(0, T) (38) Ẋ(0, T) P (T) = ê ê = Σ Ẋ(0, T) ΣẊ(0, T) = ΣP (T)Σ (39) Ẋ(0, T) are the projection operators parae to the fow of the cassica orbit and perpendicuar to the energy she, respectivey. These reations ead to the expressions d dt ( 1 dn(t) n det (m(t) 1) = Tr N(T 0 ) T=T0 dt T=T0 ), (40) with d ( dt n det (m (T) 1) = Tr N(T 0 ) 1 d N(T) T=T0 dt T=T0 1 dn(t) 1 dn(t) N(T 0 ) N(T dt 0 ) T=T0 dt T=T0 ), (41) 11

12 dn(t) dm(t, T) = dt dt d N(T) = d M(T, T) + d P (T) Σ d P (T) Σ. (43) dt dt dt dt + dp (T) dt Σ dp (T) Σ, (4) dt By using the derivatives X (n) (t, T 0 ) and the definition of the projection operator P (T), one obtains: dp (T) dt = ( 1 Ẋ (1) (0, T) Ẋ(0, T) + Ẋ(0, T) Ẋ (1) (0, T)) Ẋ(0, T) Ẋ(0, T) Ẋ (1) (0, T) P (T), Ẋ(0, T) d ( P (T) 1 = Ẋ dt () (0, T) Ẋ(0, T) + Ẋ(0, T) Ẋ () (0, T) Ẋ(0, T) ( + Ẋ (1) (0, T) Ẋ ) Ẋ(0, T) Ẋ (1) (0, T)) (1) (0, T) Ẋ(0, T) Ẋ(0, T) Ẋ () (0, T) Ẋ (1) (0, T) Ẋ (1) (0, T) Ẋ(0, T) P (T) Ẋ(0, T) Ẋ(0, T) Ẋ (1) (0, T) ( 4 Ẋ (1) (0, T) Ẋ(0, T) Ẋ(0, T) + Ẋ(0, T) Ẋ (1) (0, T) ). The ast pieces which are needed for cacuating C 0 (T) are the first and the second derivative of the monodromy matrix M(T, T) with respect to the period T: dm(t, T) dt d M(T, T) dt M(t = T, T) dt M(t = T, T) = + t dt T = Ṁ(T, T) + M (1) (T, T), = M(T, T) + Ṁ (1) (T, T) + M () (T, T). (44) Simiar to the procedure which was used in section.4.1 to arrive at the equations of motion for the X (n) (t, T 0 ), it is possibe to obtain the differentia equations governing the evoution of the derivatives M (n) (T, T) from differen- 1

13 tia equation (14) of the monodromy matrix. The equations read Ṁ (1) ij (t, T) = Σ ik ( H,km ( X(t, T))M j (t, T)X (1) m (t, T) + H,k( X(t, T))M (1) j (t, T)), Ṁ () ij (t, T) = Σ ik ( with initia vaues H,kmn ( X(t, T))M j (t, T)X (1) + H,km ( X(t, T))M (1) j (t, T)X(1) m (t, T) m (t, T)X(1) n (t, T) + H,km ( X(t, T))M j (t, T)X () m (t, T) + H,k( X(t, T))M () j (t, T) ), (45) M (n) (t = 0, T) = 0. (46) Now we have obtained a resuts necessary for the cacuation of the firstorder corrections to Gutzwier s trace formua without taking into account symmetries. In the next sections we wi appy this theory to the diamagnetic hydrogen atom. 3 The hydrogen atom in a uniform magnetic fied The diamagnetic hydrogen atom was often used as an exampe for a quantum system whose cassica dynamics is chaotic (see e.g. [8] or [9] for an overview). As a rea physica system it was the topic of studies in experimenta physics [10,11]. It has even been used for the numerica test of the correction terms C 1 and C1 T E [5,6]. Because of its simpe scaing property, which is aso fufied for the corrections, it is possibe to compare the semicassica resuts for individua orbits with exact quantum mechanica cacuations. The hydrogen atom in a uniform magnetic fied has a continuous symmetry, namey the rotationa invariance around the magnetic fied axis. This symmetry can be used to formuate the dynamics in a two-dimensiona coordinate system. This method works very we for the eading order but it eads to new difficuties for the first-order correction. This aspect is discussed in section 3.1 before we wi ook at the two-dimensiona hydrogen atom with its discrete symmetries. 3.1 The rotationa invariance and correction terms Written in the cyindrica coordinates = x + y, ϕ and z and atomic units, the cassica Hamitonian of the hydrogen atom in a uniform magnetic 13

14 fied in z-direction has the form: H = p + p z + z 1 r γ, (47) where r = + z, (48) γ = B B 0, B T. (49) In this expression z is the z component of the anguar momentum, which is a constant of motion. The paramagnetic term is not considered in equation (47) because it is constant. If one expoits the cassica scaing property the scaed Hamitonian r = γ /3 r, p = γ 1/3 p, t = γt, (50) γ /3 H = H = p + p z + z 1 r (51) ony depends on the scaed energy ǫ = γ /3 E, but not on the energy E or the magnetic fied strength γ separatey. The reguarization of the Couomb singuarity can be achieved by introducing semiparaboic coordinates [8] and the scaed time τ µ = r + z, ν = r z (5) d t = r dτ = ( µ + ν ) dτ. (53) This transformation eads to the new Hamitonian H = 1 p µ + 1 ( p ν + z 1 µ + 1 ) ǫ ( µ + ν ) + 1 ν 8 µ ν ( µ + ν ) =, (54) where the momenta are defined by p µ = dµ dτ, p ν = dν dτ. (55) The scaed time τ is the parameter for the integration of the differentia equations of the cassica vaues. In the case z = 0, the µ and ν coordinates, which are ony defined as positive coordinates in equation (5), can be extended to negative vaues. Then, p µ and p ν have the same structure as the momenta of Cartesian coordinates. This structure of the momenta is very important because it is often used in the derivation of the corrections (see e.g. equation (15) for the monodromy matrix). However, it cannot be found in the quantum 14

15 mechanica anaogue of equation (54). Starting again with the Hamitonian in cyindrica coordinates H = 1 ( 1 ) z + m 1 r + 1 γ (56) 8 and using the transformations (50) and (5), one arrives at the Schrödinger equation for the (µ, ν)-part of the wave function: { γ/3 ( 1 µ µ µ µ + 1 ν ν ν ) ν m µ m ν µ ν (µ + ν ) ǫ(µ + ν ) } ψ(µ, ν) = ψ(µ, ν). (57) In this quantum mechanica expression, the momentum operators for µ and ν differ from those in Cartesian coordinates even for m = 0, which is the quantum anaogue of z = 0. As was mentioned above, this difference does not have any effect on the eading order of the semicassica eve density. The resuts of the cassica cacuations agree very we with vaues which were extracted from quantum cacuations. But there is no agreement at a between the cassicay evauated correction terms C 1 and C1 T E and quantum mechanica cacuations. This behaviour can be understood, when one introduces the new wave function ψ(µ, ν) = which eads to the foowing Schrödinger equation: { γ/3 ( ) µ + ν + γ /3m 1 4 u(µ, ν) µν, (58) ( 1 µ + 1 ) ν µ ν (µ + ν )+ ǫ(µ + ν ) } u(µ, ν) = u(µ, ν). (59) Now, the momentum operators agree with the structure in the cassica case but an additiona term 1 8 γ/3 (1/µ + 1/ν ) appears, which eads to a nonvanishing centrifuga part of the potentia even for m = 0. The additiona term has aso to be inserted in the cassica Hamitonian: H = 1 p µ + 1 p ν + z γ/3 4 ( 1 µ + 1 ) ǫ ( µ + ν ) + 1 ν 8 µ ν ( µ + ν ) =. (60) Since the new centrifuga part is proportiona to, which was set to one when atomic units were introduced, it contributes to the first-order correction in the expansion of the eve density as one can easiy see when one repeats the cacuations from section. That is the reason why the infuence of the continuous symmetry does not appear in the eading order but in the firstorder corrections. 15

16 The centrifuga term consists of an integra over the cassica periodic orbit with the kerne 1/µ + 1/ν, which diverges at the coordinate axes. This probem does not arise if one regards the hydrogen atom as a purey twodimensiona system because in this case there is no centrifuga term. The momenta have the structure of the momenta of Cartesian coordinates in the cassica dynamics as we as in the quantum case. Therefore, in the foowing sections we ony ook at this two-dimensiona hydrogen atom for the cacuation of the first-order corrections, as it was done in [5] before. 3. Cassica dynamics and Schrödinger equation of the two-dimensiona hydrogen atom The two-dimensiona diamagnetic hydrogen atom can be formuated simiar to the three-dimensiona one. For a magnetic fied B = B ê x, the cassica Hamitonian is given by: H = p x + p y γ y. (61) If one uses the cassica scaing property (50) and semiparaboic coordinates for the two-dimensiona case [1] µ = the transformed Hamitonian reads: + x, ν = x, (6) H = 1 p µ + 1 p ν ǫ ( µ + ν ) µ ν ( µ + ν ) =. (63) The Schrödinger equation associated with the cassica Hamitonian (63) can be obtained by the same procedure as in the three-dimensiona case. The resut is: { + ǫ(µ + ν ) 1 } 8 µ ν (µ + ν ) ψ(µ, ν) { = γ /3 1 ( )} µ + ψ(µ, ν). (64) ν In this case, the momentum operators have exacty the same structure as for Cartesian coordinates. Note that γ 1/3 takes the pace of, which is equa to one in atomic units and which is often caed effective. 16

17 ν (a) (b) ν µ µ Fig. 1. Equipotentia contours of the potentia (65) for different scaed energies ǫ are potted in figure (a). The shadowed area in the coordinate system of figure (b) marks the fundamenta domain. 3.3 The symmetry of the potentia The potentia of the two-dimensiona diamagnetic hydrogen atom in semiparaboic coordinates (see equation (63)) V (µ, ν) = ǫ ( µ + ν ) µ ν ( µ + ν ) (65) has a C 4v -symmetry. This symmetry can be seen in figure 1 (a), in which a few equipotentia contours of the potentia (65) are potted. For the eading order (Gutzwier s trace formua), it is known, how the symmetry of the system can be expoited in the cacuation of the cassica vaues. In section 4 we ook at the infuence of the discrete C 4v -symmetry on cassica vaues which contribute to the first-order correction terms C 1 and C1 T E. One has to integrate differentia equations aong the periodic orbits of the cassica system in order to cacuate the required cassica vaues. It is known that, because of the symmetry, the cassica cacuations can be reduced to a fundamenta domain [13], which is shown in figure 1 (b). It consists of one eighth of the fu coordinate pane. If during the integration aong an orbit one arrives at one of the borders of the fundamenta domain, one refects the orbit at the border. This is possibe because the borders of the fundamenta domain are refection panes of the potentia. In addition, the symmetry of the potentia aows of a compact abeing of the unstabe cassica periodic orbits. In the symmetry reduced fundamenta domain, the orbits are described by a ternary code, which was introduced by Eckhardt and Wintgen [14]. It uses the symbos 0, + and and has its origin in the description of the orbits of the four disk scattering system, which has the same symmetry as the diamagnetic hydrogen atom in semiparaboic coordinates. 17

18 The restriction to the fundamenta domain and the introduction of symmetry operations ead to new periodic orbits, which are ony periodic if one expoits the symmetry properties of the system. The periodic continuation can be achieved by refections at the borders of the fundamenta domain. An exampe is the orbit +, which is shown in figure (a). Starting for exampe on the µ-axis the orbit has to be refected at the ange bisector. After returning to the µ-axis a second refection (this time at the µ-axis) is necessary such that the momenta at the initia and fina points agree. For practica purposes it is often easier not to restrict the cacuation to the fundamenta domain but to find a periodic continuation of the orbit by mapping the fina point of the orbit on its initia point via a symmetry transformation from C 4v, namey rotations by mutipes of 90 degrees (c 4, c 4 = c, c 3 4 ), refections at the coordinate axes (σ v ), and refections at the ange bisectors (σ d ). Then one can obtain the same new periodic orbits as described above. For exampe figure (b) shows the orbit from figure (a) in the case where it is not restricted to the fundamenta domain. The periodic continuation is done by a cockwise rotation by 90 degrees. This method requires ony one symmetry operation to render an orbit periodic, and is in genera easier to impement as the restriction to the fundamenta domain. Figure (c) shows the orbit which corresponds to the doube repetition of the orbit + in the fundamenta domain and figure (d) shows the orbit which consists of four times the same repetition in the fundamenta domain and which is periodic in the pane of semiparaboic coordinates. 3.4 Symmetry properties of the wave functions and cacuation of the quantum spectra Because of the C 4v -symmetry of the Hamitonian (64), the eigenfunctions of the system spit up into subspaces beonging to a representation of the symmetry group. The symmetry group C 4v has four one-dimensiona representations, namey A 1, A, B 1 and B, and a two-dimensiona representation, which is caed E. Appying a symmetry eement from C 4v to a wave function with symmetry E eads in genera to a inear combination of two (energeticay degenerate) wave functions. At fixed scaed energy ǫ, equation (64) can be considered as a generaized eigenvaue probem in the variabe γ /3. The eigenvaues can be cacuated by diagonaizing a matrix representation of the Hamitonian in a compete basis set with the Lanczos agorithm [15,5]. Using ony the wave functions of one of the subspaces, which means that the bock diagona form of the Hamitonian is expoited, reduces the dimension of the eigenvaue probem and eads to separate spectra for each of the subspaces, which is necessary for the anaysis 18

19 (a) (b) 1 1 ν 0 ν µ µ (c) (d) 1 1 ν 0 ν µ µ Fig.. The periodic orbit + in semiparaboic coordinates for a scaed energy of ǫ = 0.1. The dot marks the nuceus at the origin. The orbits shown correspond to the singe (see figures (a) and (b)), the doube (c) and the quadrupe (d) repetition in the fundamenta domain. of the infuence of the symmetry. Up to a maximum γ 1/3 (the eigenvaues are needed in this form) of about 156 one obtains circa 1,000 eigenvaues in each of the one-dimensiona subspaces A 1, A, B 1, and B. In the subspace which beongs to the two-dimensiona representation E one finds roughy 19,000 eigenvaues up to a maximum γ 1/3 of about Anaysis of the quantum spectrum by harmonic inversion If one ooks at the semicassica eve density g osc (E) = 1 π T 0p det (m (T 0 ) 1) { cos ( 1 S(c) (T 0 ) π µ ( 1 C sin S(c) (T 0 ) π ) µ + O ( )}, (66) ) 19

20 which one obtains from the semicassica approximation of the quantum Green s function () via the reation 1 π Im (Tr G(E)) = g(e) = n δ(e E n ), (67) one recognizes that it shoud be possibe to extract the ampitude of the eading order (Gutzwier s trace formua), A = T 0p det (m (T 0 ) 1), (68) and the reduced action S (c) from a Fourier transformation of the quantum eve density. A higher precision in both the ampitude and the action can be achieved with the harmonic inversion method (see e.g. [9]), where the soution of the noninear set of equations g osc (w n ) = Ā e iwn S (69) yieds the parameters Ā and S. The w n are chosen on an equidistant grid. For a spectrum which consists of δ-functions such as that in equation (67), one has to modify it in order to obtain non-vanishing contributions on the grid, e.g., by a convoution of the δ-functions with a Gaussian or by appying a fiter to the signa. We used the atter method, which is described in [16]. Using the harmonic inversion method makes it possibe to gain the cassica vaues from an anaysis of the quantum eigenvaues and to compare these resuts with those from cassica cacuations. The harmonic inversion method can aso be appied to extract higher-order corrections C = C 1 + C T E 1 (70) for individua orbits in the foowing way [5]. The quantum eve density formuated in semiparaboic coordinates depends on the variabe w γ 1/3 : δ(h ) = n w n ψn (µ, ν) 1 w n 1 ( ) + ψn (µ, ν) δ(w w n ), (71) µ ν where H is the quantum Hamitonian which beongs to the cassica Hamitonian (63). Using the cassica scaing property (63) of the system, the semicassica approximation of the quantum Green s function and the reation (67) 0

21 ead to the eve density g osc (E) = w π T ( 0p det (m (T 0 ) 1) cos ws (c) (T 0 ) π ) µ 1 ( w C sin ws (c) (T 0 ) π ) µ + O ( ( ) 1 ) w. (7) Combining the formuas for the quantum eve density (71) and for its semicassica approximation (7) shows that the harmonic inversion of g(w) = n π ψn (µ, ν) 1 w n 1 ( ) + ψn (µ, ν) δ(w w n ) (73) µ ν eads to an expansion of the form (69), where Ā = 1 [ A exp iπ ] µ, S = S (c) (T 0 ) (74) with A from equation (68). This method can provide access to the first-order correction, if one subtracts the eading order and mutipies the resut by w. The fina expression is n πw n ψn (µ, ν) +w 1 w n 1 ( ) + ψn (µ, ν) δ(w w n ) µ ν ( (T 0 ) π µ T 0p det (m (T 0 ) 1) cos ws (c) ), (75) and the harmonic inversion eads to the ampitude Ā = 1 [ i A exp iπ ] µ C, (76) from which one can easiy extract the first-order correction C because a other vaues are known. 4 Infuence of discrete symmetries on the correction terms In this section we investigate the infuence of the discrete C 4v -symmetry on the correction terms C 1 and C1 T E. The corrections are computed for some orbits which are not periodic without a symmetry transformation from C 4v. 1

22 4.1 Symmetry transformations in the cacuations of the correction terms The cassica quantities have to be cacuated for periodic orbits and, as we have seen in section 3, some orbits are ony periodic after the appication of a symmetry operation during or at the end of the integration aong the cassica orbit. It is we known which symmetry operations (refections or rotations of vectors) have to be impemented for the cacuation of the phase space coordinates but, as was mentioned in section, one has to sove a arge number of additiona differentia equations for a new set of coordinates if one wants to cacuate the first-order correction. For exampe the inearized equation of motion (14) has to be soved in order to obtain the monodromy matrix. In addition, the derivatives of the coordinates of the orbit and of the monodromy matrix with respect to the period are required. Thus we face the question of how the symmetry has to be impemented into the correction terms. What are the correct transformations for the additiona coordinates? What is the transformation of the boundary conditions of the cassica Green s function? Correction term C 1 The soution of the inearized equation of motion (14) for the monodromy matrix is essentia for the cacuation of the cassica Green s function by using the formuation (13) and by soving equation (19). The monodromy matrix is one of the cassica vaues which contribute to Gutzwier s trace formua, and its symmetry behaviour has often been used. However, we ook at the transformations of its eements because it is an exampe for a other variabes which are necessary for the first-order corrections. A further eements foow the same scheme. The monodromy matrix can be obtained by soving equation (14) for four inear independent coumn vectors µ M ν X M = M p µm p νm with initia vaues which match the condition M(0) = 1. In our system, the monodromy matrix has the dimensions 4 4 and the inearized equations of

23 motion in semiparaboic coordinates are: µ M = p µm, ν M = p νm, ṗ µm = (ǫ 3 µ ν 1 ) 4 ν4 µ M ( µ 3 ν + µν 3) ν M, ṗ νm = (ǫ 3 µ ν 1 ) 4 µ4 ν M ( µ 3 ν + µν 3) µ M. (77) As an exampe, we ook at the anticockwise rotation by an ange of 90 degrees. In this case the symmetry transformation for the coordinates and momenta of the orbit eads to µ M = p µm, ν M = p νm, ṗ µm = (ǫ ν µ 3 1 ) 4 µ4 µ M + ( ν 3 µ + ν µ 3) ν M, ṗ νm = (ǫ ν µ 3 1 ) 4 ν4 ν M + ( ν 3 µ + ν µ 3) µ M, (78) where µ and ν are the new variabes in the rotated system. Because of the symmetry invariance of the Hamitonian, one has to find a transformation in such a way that µ M = p µm, ν M = p νm, p µm = (ǫ 3 µ ν 1 ) 4 ν4 µ M ( µ 3 ν + µ ν 3) ν M, p νm = (ǫ 3 µ ν 1 ) 4 µ4 ν M ( µ 3 ν + µ ν 3) µ M. (79) In our exampe this condition eads to the resuts: µ M = ν M, ν M = µ M, p µm = p νm, p νm = p µm. (80) The condition is aways fufied if one uses for the eements of X M the same transformations as for the corresponding phase space coordinates of the orbit. Since the matrices A and B can directy be determined from equation (0), there is no need for further symmetry operations in order to obtain the right vaues for these two matrices if the monodromy matrix was cacuated with a symmetry transformations. The same is true for the matrices A + and B +, which foow from equation (17), and for the cassica Green s function, which foows from (13). This is evident if one recas the cacuations from section.. Furthermore, if one appies the correct symmetry transformation to the eements of the monodromy matrix at every time when one transforms the coordinates of the orbit during the integration, a boundary conditions of the 3

24 cassica Green s function are fufied because the monodromy matrix is the ony part that foows from the soution of a differentia equation. A other vaues foow from formuas which incude aready the boundary conditions and are not affected by the symmetry transformations during the integration Correction term C T E 1 In addition to the coordinates considered unti now, one has to obtain the vaues of the derivatives X (n) and M (n) for the second correction term C T E 1. The equations of motion of the components of X (1) Ẋ (1) i (t, T 0 ) = Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0) (81) have exacty the same structure as the differentia equations for the monodromy matrix, i.e., we aready know the symmetry operations which have to be impemented here. Looking at the differentia equations of X () and X (3), namey Ẋ () i (t, T 0 ) = Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0)X (1) (t, T 0 ) + Σ ij H,jk ( X(t, T 0 ))X () k (t, T 0), Ẋ (3) i (t, T 0 ) = Σ ij H,jkm ( X(t, T 0 ))X (1) k (t, T 0)X (1) (t, T 0 )X m (1) (t, T 0) + 3Σ ij H,jk ( X(t, T 0 ))X (1) k (t, T 0)X () (t, T 0 ) + Σ ij H,jk ( X(t, T 0 ))X (3) k (t, T 0), (8) one recognizes that the terms which contain the X () k and the X (3) k, respectivey, have the same form as the equation of motion of the monodromy matrix. Therefore, we have to use the same transformation. The differentia equations of the eements of M (1) and M (), Ṁ (1) ij (t, T) = Σ ik ( H,km ( X(t, T))M j (t, T)X (1) m (t, T) + H,k( X(t, T))M (1) j (t, T) ), Ṁ () ij (t, T) = Σ ik ( H,kmn ( X(t, T))M j (t, T)X (1) m (t, T)X(1) n (t, T) + H,km ( X(t, T))M (1) j (t, T)X(1) m (t, T) + H,km ( X(t, T)) M j (t, T)X m () (t, T) + H,k( X(t, ) T))M () j (t, T) (83) are derived in the same way as those for, e.g., the X () k by differentiating Hamiton s equations of motion, and have the same structure. One can easiy see that the eements of the derivatives M (1) and M () obey the same symmetry transformation as the eements of the monodromy matrix. 4

25 We can summarize the resuts of this section as foows: Under symmetry operations a variabes which contribute to the first-order correction transform in the same way as the corresponding phase space coordinates of the orbit. Whenever a symmetry transformation of the phase space coordinates of the orbit is necessary during the integration of a set of differentia equations in order to obtain a periodic orbit, one has to transform the variabes which were discussed in this section simutaneousy. 4. An exampe for a symmetry reduced orbit As an exampe for an orbit which is periodic ony when a symmetry transformation is appied, we ook at the orbit +, which has been introduced in section 3.3 and is shown in figures (a) and (b) in the fundamenta domain and in the fu (µ,ν)-pane, respectivey. As a first test for the correct impementation of the symmetry in the correction term C 1, the cassica Green s function G(τ, τ ) is considered. In order to iustrate the infuence of the symmetry transformation, the eements of the cassica Green s function of the orbit in the fundamenta domain (figure (a)) is shown in figure 3 (a). The Green s function, which in our case is a matrix, was rotated in such a way that the G 1j (τ, τ ) components represent the direction aong the cassica orbit at the initia point, which is marked by a cross in figure 3 (b). One recognizes the discontinuities at the positions of the refection which appear due to the change in the meaning of the components according to the change of the variabes by using the symmetry transformation. A boundary conditions, which were mentioned in the sections.1 and., are fufied. At the fina point of the orbit, a components have the same vaues as at the initia point. The condition G 1j (0, τ ) = G 1j (T 0, τ ), which is expected from equation (10), is aso fufied. The discontinuities in the eements G j, which appear due to equation (16), are ceary visibe. Figure 3 (c) shows the cassica Green s function for the same orbit, which is not restricted to the fundamenta domain but cosed by a rotation of its fina point (see figure 3 (d)). Therefore, the ony discontinuity appears at the end of the orbit. To show this, the first and the ast point (after the symmetry transformation) of each component are marked with the same symbo. Again, the boundary conditions are fufied. 5

26 (a) G ij (τ,τ, ) G 11 G 1 G 1 G (b) -0. ν τ/t 0 µ (c) G 11 G 1 G 1 G G ij (τ,τ, ) (d) -0.1 ν τ/t 0 µ Fig. 3. The cassica Green s function of the orbit + is shown in two cases for τ = 0.55 T 0. The upper two pictures (a) and (b) represent the orbit in the fundamenta domain, whereas, the two ower ones (c) and (d) stand for the same orbit in the fu pane of semiparaboic coordinates, which becomes periodic by the appication of a rotation of the fina point. In both cases the eements of the Green s function are chosen in such a way that the components G 1j correspond to the direction aong the orbit at the initia point, which is marked by a cross in figures (b) and (d). In the fundamenta domain (figure (a)), where two refections are required, one can ceary see the discontinuities in the eements of the Green s function at the two positions at which the orbit is refected. The discontinuities correspond to the change of the meaning of the components. In figure (b) ony the ast point is affected by the symmetry operation (cockwise rotation by 90 degrees). This effect has been visuaized by marking the ast point of each component with the same symbo as the first point of the same component. The ast point was potted after the rotation has been carried out. In both cases the boundary conditions are fufied. 6

27 Tabe 1 Character tabe of the symmetry group C 4v representation e c 4, c 3 4 c 4 = c σ v σ d A A B B E Infuence of the symmetry behaviour of the quantum wave functions and comparison of semicassica and quantum resuts for the correction terms In the different subspaces, which beong to a representation of the symmetry group, the individua ampitudes A of the orbits can be found with a prefactor, which is given by the character χ (see tabe 1) of the eement from C 4v which provides the periodic continuation of the orbit in the representation of the corresponding subspace. The cause of the prefactor is the projection [13] of the eve density on the subspace beonging to the representation α of the group G (in our case G = C 4v ) with the projection operator P α = d α G X G χ α (X)D α (X). (84) Here, G stands for the order of the group ( C 4v = 8) and d α for the dimension of the representation D α. The eement D α (X) operates on the fina point of the orbit with the matrix representation M(X 1 ). Ony if the symmetry eement X 1 eads to a periodic continuation of the orbit, it wi contribute to the eve density. For the weight of the orbit, one has to take into account that the fundamenta domain appears eight times in the compete coordinate pane and that aways eight identica orbits appear in the fu (µ, ν)-pane. As a consequence of the symmetry (eq. (84)), the orbit +, which requires the rotation c 3 4 in order to become periodic, shoud appear in the subspaces A 1 and A with the prefactor p α = 1, in B 1 and B with p α = 1, and in the subspace beonging to the representation E it shoud not appear at a. Tabe shows the resuts of the harmonic inversion of the eve density which consists ony of the eigenvaues of one of the one-dimensiona subspaces. The quantum vaues are compared with cassica cacuations. The action is independent of the subspace, and its cassica vaue is S (c) = The quantum ampitudes A (qm) incude information from the Masov index µ (see equation (74)). Therefore, they have compex vaues, and the moduus and phase are 7

28 Tabe The actions S (qm) and the ampitudes A (qm) of the orbit + which were extracted from the quantum spectrum are compared with their cassica counterparts. The cassica vaue of the action is S (c) = , independenty of the subspace. Because of the projection to the individua subspaces the moduus of the quantum mechanica ampitudes corresponds to A (qm) = p α A A (c) and the phase is given by arg(a (qm) ) = π µ + arg(p α ) arg(a (c) ). Since the quantum mechanica phases are ony determined moduo π, a arguments are shifted into the standard interva ( π,π]. The orbit + has a vanishing ampitude in subspace E. ( ) ( ) subspace A (c) arg arg S (qm) A (qm) A (qm) A π π/ A π π/ B π π/ B π π/ A (c) Tabe 3 Ampitude of the orbit ++ in different subspaces. See tabe for an expanation. ( ) ( ) subspace A (c) arg arg given by S (qm) A (qm) A (qm) A π π A π π B π π B π π E π 0 A (qm) = p α A A (c), arg ( A (qm) ) π = µ + arg (p α ) arg ( A (c) ), A (c) where p α is the prefactor which appears due to the projection operator P α. In tabe a ampitudes and phases have the expected vaues. The agreement is very good. A arguments are shifted into the standard interva ( π, π]. Tabe 3 shows the same comparison for two repetitions (counted in the fundamenta domain) of the orbit +, which are abeed with the symbo ++. Whie the orbit + vanishes in subspace E, the orbit ++ can be found in this subspace with the expected ampitude prefactor p α = 4. It is important to use the ampitude with the correct prefactor in equation (75) for the extraction of the sum C of the two first-order correction terms. 8

29 Tabe 4 Correction term C of the orbit + in different subspaces. The moduus C (qm) and phase arg(c (qm) ) of the anaysis of the quantum spectra are compared with the cassicay cacuated vaue C (c). subspace C (c) C (qm) ( arg C (qm) ) C (qm) C (c) 1 A π A π B π B π If this condition is taken into account, one obtains the correction terms in every subspace in which the orbit appears. Tabe 4 shows the resuts for the orbit + in a one-dimensiona subspaces of C 4v, which are compared with the cassicay cacuated vaue. In tabes 5 and 6, the semicassica resuts for the correction terms are compared with vaues which were extracted from exact quantum cacuations. Ony the eigenvaues of the subspace beonging to the representation A 1 were used for the harmonic inversion of the quantum spectrum. In this subspace a orbits contribute with a prefactor of 1 independenty of the symmetry eement which is required to find the periodic continuation. Quantum mechanica eigenvaues up to a maximum w of about 156 were used for the anaysis. In tabe 5, cassicay cacuated correction terms are given for some orbits which are ony periodic with a symmetry transformation. The orbits 0 + and 0 + have the same shape but are traversed in opposite directions. In a other cases mentioned in the tabes, the symbo is independent of the direction because the orbit is identica with its time reversed counterpart. The cassica vaues are compared with the resuts from the anaysis of the quantum spectrum in tabe 6. The agreement of the ampitudes from cassica and quantum cacuations is again very good. In most cases the differences are ony of the order In a three cases in which the difference is arger, the actions of these orbits ie cose to those of other orbits. This makes the anaysis of the quantum spectrum difficut. The phases of the quantum mechanicay cacuated correction terms reproduce the correct signs of the cassica vaues but the differences are of the order 10. A simiar behaviour was found in [5] for orbits which are periodic without any symmetry transformation. 9