Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the semiclassical quantisation of non - integrable systems using properties of their periodic orbits was triggered by m. gutzwiller @xcite and extended by several groups @xcite . it allows one to express the oscillating part @xmath7 of the quantum - mechanical density of states , given exactly in terms of a quantum spectrum @xmath8 and separated into two terms by g(e ) = _ n ( e - e_n ) = g(e)+g(e ) , [ gqm ] through a semiclassical trace formula of the form g_sc(e ) _ po a_po(e ) [ s_po(e)/-_po/2 ] . [ dgsc ] the sum is over all periodic orbits @xmath9 of the classical system , @xmath10 are their action integrals , the amplitudes @xmath11 depend on their stabilities and degeneracies , and @xmath12 are some phases called maslov indices . the average part @xmath13 of the density of states , which by definition varies smoothly with energy , is obtained by the extended thomas - fermi ( etf ) model ( see , e.g. , @xcite , chapter 4 ) . then , the sum @xmath14 usually turns out to be a good approximation to the exact quantity , although the sum over the @xmath15 in is only an asymptotic one , correct to leading order in @xmath16 , and in chaotic systems is hampered by convergence problems @xcite . for integrable systems , berry and tabor @xcite and later creagh   @xcite showed how a trace formula of the form can quite generally be obtained from the einstein - brillouin - keller ( ebk ) quantisation @xcite . this method will be used in the present paper for a spherical system . the semiclassical theory , often referred to as periodic orbit theory ( pot ) , has been very successful not only in ( partial ) quantisation of a given hamiltonian , but also in interpreting many experimentally observable quantum oscillations in finite fermion systems  so - called `` shell effects ''  in terms of classical mechanics . examples are atomic nuclei @xcite , metallic clusters @xcite , semiconductor quantum dots @xcite , and metallic nanowires @xcite . an overview over many aspects of the pot and illustrative applications are given in @xcite . one problem that remains with all the trace formulae developed in the work quoted above is that the amplitudes @xmath17 diverge in situations where a continuous symmetry is broken or restored under the variation of a system parameter ( which may also be the energy @xmath18 ) , or where bifurcations of periodic orbits occur . in such situations one has to go beyond the stationary phase approximation which is underlying the semiclassical approach . this has , besides @xcite , been developed most systematically by ozorio de almeida and hannay @xcite for both situations , leading to local uniform approximations with finite amplitudes @xmath17 . so - called global uniform approximations , which yield finite amplitudes at symmetry - breaking and bifurcation points and far from them go over into the standard ( extended ) gutzwiller trace formula , were developed for the breaking of u(1 ) symmetry in @xcite , for some cases of u(2 ) and so(3 ) symmetry breaking in @xcite , and for various types of bifurcations in @xcite . ( details and further references may be found in @xcite , chapter 6.3 . )    in this paper we investigate an anharmonically perturbed three - dimensional spherical harmonic oscillator ( ho ) hamiltonian ( for a particle with mass @xmath19 ) : h ( , ) = 12p^2 + 12 ^ 2 r^2 + 14r^4 = h_0(r , p ) + h(r ) , [ hpert ] where @xmath20 and @xmath21 are the absolute values of the three - dimensional momentum and radius vectors . @xmath22 is the strength of the perturbation which first is assumed to be small , but later may assume arbitrary positive values . the unperturbed ho system has u(3 ) symmetry @xcite which is broken by the anharmonic term , leading to the so(3 ) symmetry of the perturbed spherical system . the hamiltonian suggested itself in a recent study @xcite of harmonically trapped fermionic atoms with a short - range repulsive two - body interaction , treated self - consistently in the hartree - fock approximation . as a result , very pronounced shell effects in the total energy @xmath23 of the interacting system as a function of the number @xmath24 of atoms were found , which remind about the so - called super - shells predicted @xcite and observed @xcite in metallic clusters . in a first attempt to interpret these shell effects semiclassically @xcite , we parameterised the self - consistent mean field of the interacting system by the hamiltonian and applied the semiclassical perturbation theory of creagh @xcite to explain qualitatively the shell structure of the hf results .    in the present paper we describe some of the mathematical details of the semiclassical perturbation theory for the symmetry breaking u(3 ) @xmath25 so(3 ) and develop a uniform trace formula valid for arbitrary strengths @xmath1 in the hamiltonian . in section [ secpert ] we present the perturbative trace formula for the density of states which is valid in the limit of small @xmath1 and has been presented shortly in @xcite . it already puts the dominance of the shortest periodic orbits , namely the circles and diameters , into evidence . although their individual contributions to the perturbative trace formula diverge in the limit @xmath26 , their sum restores to the unperturbed ho trace formula with u(3 ) symmetry in this limit , which at the same time is the limit of zero energy . in section [ secuni ] we develop the uniform trace formula that includes the contributions of the diameter and circle orbits for arbitrary values of @xmath1 . for this purpose , we start from the ebk quantisation of an arbitrary system with spherical symmetry and apply the poisson summation formula . we obtain a one - dimensional trace integral for the semiclassical density of states , whose end - point contributions correspond to the diameter and circle orbits , valid for an arbitrary spherical one - body potential . for the hamiltonian , the gross - shell structure of the density of states is at low energies always dominated by the families of diameter and circle orbits , although these orbit families only have a two - fold degeneracy . at sufficiently high energies and repetition numbers , `` rational tori '' with frequency ratios @xmath27 with @xmath28 bifurcate from the circle orbits with repetition numbers @xmath29 , as discussed in section [ sectori ] . these rational tori have the highest ( three - fold ) degeneracy possible in a three - dimensional spherical system , but due to their length they only affect the finer quantum structures of the density of states . this situation is completely different from that of a spherical billiard @xcite where the contribution of the diameter orbits to the density of states is practically negligible .    in the limit @xmath30 , becomes a purely quartic oscillator which in many respects is easier to handle . in particular , in this limit the system acquires the `` scaling property '' that its classical dynamics becomes independent of the energy which can be absorbed by a rescaling of coordinates and time . this limit will be discussed separately in section [ secquart ] . one interesting classical aspect is that no bifurcations of the circle orbits occur . the same rational tori , which bifurcate from the circle orbit in the perturbed ho system , exist here at all energies . but , again , they affect only the finer quantum structures of the density of states , while its gross - shell structure is largely dominated by the circles and diameters .    in the appendices [ app1 ] - [ app3 ] , we collect some mathematical details about the integration over the manifolds @xmath2p@xmath3 and s@xmath31 and some explicit analytical expressions for action integrals and periods in terms of elliptic integrals . finally , in the appendix [ appothers ] , we re - derive from our general trace integral the known trace formulae for two popular spherical systems : the spherical billiard and the coulomb potential . we first write down the exact density of states of the unperturbed three - dimensional harmonic oscillator ( ho ) , given by the hamiltonian h_0(r , p ) = 12 p^2 + 12 ^ 2 r^2 . [ hho ] its quantum - mechanical density of states can be written in the form g_0(e ) = _ n=0^d_n ( e - e_n ) = g_0^etf(e ) + g_0(e ) , [ gsum ] using the spectrum @xmath32 with the degeneracy factor @xmath33 . here @xmath34 is the smooth part given by the etf model g_0^etf(e ) = , [ getfho ] and @xmath35 is the oscillating part given by the following exact trace formula @xcite g_0(e ) = 2g_0^etf(e)_k=1^(-1)^k ( 2k ) = 2g_0^etf(e)_k=1^(-1)^k . [ dg3ho ] it can be understood as a sum over all periodic orbits of the system . hereby @xmath36 represents the repetition number of the primitive classical periodic orbit which has the action @xmath37 . keeping only the leading tf term of @xmath34 in , the density of states can also be written in the form g_0(e ) = e _ -^(-1)^ke^i2ke/ + o(^-1 ) , [ g3hoc ] since the imaginary parts cancel upon summation . note that the smooth tf part in comes from the contribution with @xmath38 in . next we follow creagh @xcite by writing the perturbed trace formula in the form g_pert(e ) = e _ -^(-1)^k m(k/)e^i2ke/ , [ gpert ] where @xmath39 is a modulation factor that takes into account the lowest - order perturbation of the action of the ho orbits . here @xmath40 is a small action that , quite generally , depends on @xmath1 and the energy @xmath18 . the factor in front of the sum in takes into account only the leading term of the smooth part @xmath34 ; this is consistent with the fact that the perturbation theory only deals with the terms of leading order in @xmath41 of the semiclassical trace formula . as we will see in section [ secqm ] , the contributions of terms of relative order @xmath42 in @xmath34 are , indeed , negligible . as described in detail in @xcite , one obtains @xmath39 from an integration of a perturbative phase function over the manifold that describes the classical degeneracy of the ho orbits . the unperturbed ho has u(3 ) symmetry , so that @xmath39 is formally obtained by the average ( k/ ) = e^iks(o)/_ou(3 ) , [ mav ] where @xmath43 is an element of the group u(3 ) characterising a member of the unperturbed ho orbit family , and @xmath44 is the lowest - order primitive action shift brought about by the perturbation @xmath45 in . in the present case @xmath44 is nonzero already in first - order perturbation theory and therefore @xmath46 is proportional to @xmath1 . we need , however , not integrate over the full nine - dimensional space of group elements of u(3 ) . it is sufficient to consider only the subset of symmetry operations which transform the orbits within a given degenerate family into each other ( without changing their actions ) . the dimension of this subset is the degree of degeneracy @xmath47 of that family . for the three - dimensional spherical ho we have @xmath48 . this can most easily be seen by the following argument which also will allow us to find a suitable parametrisation for the four - fold integration . the full phase space is six - dimensional ; due to energy conservation it is reduced to a five - dimensional energy shell which has the topology of a five - sphere s@xmath31 . of the five parameters that specify a point on s@xmath31 , one can be chosen as a trivial time shift along the periodic orbits corresponding to a simple phase factor @xmath49 . this parameter forms a subgroup u(1 ) , so that its elimination restricts us to the four - dimensional manifold ( see appendix [ app1 ] for its mathematical definition ) s^5/u(1 ) = ^2 , [ coset ] which is neither a four - torus nor a four - sphere @xcite . ( note that for the two - dimensional ho one is led to the manifold @xmath2p@xmath50 which happens to be homomorphous to the two - sphere s@xmath3 , see @xcite . ) a suitable parametrisation of @xmath2p@xmath3 is given in appendix [ app1 ] in terms of four angles @xmath51 whose meaning will become clear in a moment . the modulation factor therefore becomes ( k/ ) = d_^2 e^iks/ = _ 0^ d _ 0^ ^3d _ 0 ^ 2 d_2 _ 0 ^ 2 d_3 e^iks(,,_2,_3)/. [ mfac ] in the zero - perturbation limit @xmath26 , where the action shift @xmath52 and hence also @xmath46 is zero , the modulation factor becomes unity , as it should for to approach . next we have to determine the action shift @xmath52 of a primitive periodic orbit of the ho , caused by the perturbation @xmath53 . the harmonic solutions of the classical equations of motion of the unperturbed system are well known : x_i = ( t+_i ) , i=1,2,3 ( x_1,x_2,x_3)=(x , y , z ) . [ xi ] hereby @xmath54 are the three conserved energies in the three dimensions . depending on the values of @xmath54 and the phases @xmath55 , describes circles , ellipses or librations through the origin ; we will in the following call the latter the `` diameter orbits '' . we now re - parameterise the @xmath56 in the following way : x(t ) = rn_1(t),y(t ) = rn_2(t+_2),z(t ) = rn_3(t+_3 ) , [ xi3d ] where @xmath57 , and @xmath58 is given by r=/,e = e_1+e_2+e_3 . since the @xmath59 must lie on the sphere @xmath60 , due to the conservation of energy , we can also write them as : n_1 = , n_2 = , n_3 = . [ s2par ] actually , since the @xmath61 must be restricted to positive definite values in order to cover once all classically allowed values of the @xmath62 , we only have to integrate over the first octant of s@xmath3 , so that @xmath63 $ ] . the four angles @xmath64 , with their ranges of definition , are precisely the parameters describing the manifold @xmath2p@xmath3 as explained in appendix [ app1 ] , and the correct integration measure is that used in . we might have kept the phase angle @xmath65 of @xmath66 in , too ; integration over it is tantamount to a trivial integration over time along the orbits . the resulting five - dimensional integral becomes equivalent to that over the s@xmath31 sphere discussed in appendix [ app2 ] .    according to creagh @xcite , the first - order action shift brought about by a perturbation @xmath67 is given by _ 1 s = - _ po h(t ) dt , where @xmath15 stands for a particular member of the unperturbed periodic orbit family , and the perturbating hamiltonian @xmath68 has to be evaluated along the periodic orbit . inserting - into the perturbation @xmath69 leads to elementary integrals over powers of the trigonometric functions . the result is _ 1 s(,,_2,_3 ) & = & -(3r^4 - 4l^2/^2 ) = -(3e^2-^2l^2 ) + & = & -[3-^2(2 ) ( ^2 ^ 2_2+^2 ^ 2_3 ) -^4 ^ 2(2)^2(_2-_3 ) ] , [ ds3d ] where @xmath70 is the conserved squared angular momentum , and the first - order action unit @xmath46 is given by = e^2/4 ^ 5 . [ sigma ]    the integral , with the function given in , has been integrated numerically and found to be identical with the corresponding five - dimensional integral over the 5-sphere s@xmath31 , expressing @xmath71 in terms of the five hyperspherical angles of the six - dimensional polar coordinates ( cf . appendix [ app2 ] ) . however , even the four - dimensional integration in is more than needed . in fact , the same result is obtained if we replace the phase function under the @xmath2p@xmath3 integration as follows : _ 1 s(,,_2,_3 ) _ 1 s ( ) = -(3-^2 ) . [ ds3da ] the integrals over @xmath72 , @xmath73 and @xmath74 then become trivial and the result is found ( with @xmath75 ) to be ( k/ ) & = & 2 _ 0^ d e^-ik(3-^2)/ = _ 0 ^ 1 dze^-ik(3-z)/ = .   [ modfex ] the replacement is not obviously justifiable in a direct way , but a suitable reduction of the five - dimensional s@xmath31 integral yields exactly the integral , as demonstrated in appendix [ app2 ] . identically the same result is also obtained independently from the ebk quantisation of the perturbed hamiltonian using poisson summation and keeping only first - order terms in @xmath1 ; it is given in eq .   at the end of sect . [ secpois ] . the form of the modulation factor suggests a simple physical interpretation : the resulting two terms correspond to two separate families of orbits that survive the breaking of the u(3 ) symmetry . as we will see below , these are the circle and the diameter orbits with maximal and minimal ( i.e. , zero ) angular momentum , respectively , at fixed energy . for @xmath76 these are actually the shortest periodic orbits found in the perturbed system . in contrast to the four - fold degenerate unperturbed ho orbits , they only have a two - fold degeneracy , since only rotations about two of the euler angles change their orientation . in general , the most degenerate periodic orbits in a three - dimensional system with spherical symmetry so(3 ) can be rotated about three euler angles . however , for each of the circle and diameter orbits one of the three rotations is redundant in the sense that it maps these orbits onto themselves : for the circles , it is the in - plane rotation around their centre , and for the diameters , it is the rotation about their own axis of motion . the loss of two degrees of degeneracy of these two orbit types with respect to the unperturbed ho orbits appears in the form of the factor @xmath77 in their amplitudes in , in agreement with the general counting of the powers of @xmath77 in semiclassical amplitudes @xcite . we will see in the next section that for values @xmath78 and large enough values of @xmath1 , there exist fully three - fold degenerate orbits , namely orbits with rational ratios @xmath79 of radial and angular frequencies with @xmath80 . they are created at bifurcations of the @xmath36-the repetitions of the circle orbits with @xmath81 . this happens , however , only at finite values of @xmath1 , so that these orbits do not contribute in the limit @xmath26 . inserting into , we find the following perturbed trace formula for the oscillating part of the level density : g_pert(e ) & = & _ k=1^ + & = & _ k=1^ , [ dgp ] where @xmath82 and @xmath83 are the perturbed primitive actions of the diameters and the circle , respectively : s_d = s_0 - 3 , s_c = s_0 - 2,s_0=2e/. [ sorbits ] these values will be confirmed from the general expressions derived in the next section . the perturbative trace formula uniformly restores the oscillating part of the exact ho trace formula in the limit @xmath26 . in section [ secqm ] we will test its range of validity comparing against quantum - mechanical results . note that  as is usual in perturbation theory  the semiclassical amplitudes and actions of the orbits contributing to are correct only in the small-@xmath1 limit . in general , one has to generalise the perturbative trace formula to a uniform version that , in the limit of large perturbation , goes over into the corresponding ( extended ) gutzwiller trace formula . this uniformisation has been done for the breaking of u(1 ) symmetry in @xcite , and for some cases of u(2 ) symmetry breaking in @xcite ; one of the latter result applies , with suitable changes , also to some cases of so(3 ) @xmath25 u(1 ) symmetry breaking . the symmetry breaking u(3 ) @xmath25 so(3 ) under study here has not been treated in the literature so far . however , the simplicity of our above results makes the uniformisation particularly easy in the present case : since the modulation factor already has its own asymptotic form  or , inversely speaking : since the asymptotic expansion of the integral in in the limit of large @xmath46 happens to be exact also for @xmath84  it will be sufficient to replace in the perturbed actions @xmath82 , @xmath83 and their semiclassical amplitudes by those valid for all values of @xmath1 , which will be derived in the next section . in this section we will calculate the full actions and semiclassical amplitudes of the diameter and circle orbits of the hamiltonian , valid for arbitrary values of @xmath1 . the general trace formula for an arbitrary spherical potential in three dimensions has been given by creagh and littlejohn @xcite . we choose a different approach here , which is in spirit that of berry and tabor @xcite , but goes beyond their leading - order approximation in taking into account the end - point corrections to a trace integral which yield precisely the circle and diameter orbit contributions . the classical hamiltonian of the system ( with mass @xmath19 ) h ( , ) = e = 12p^2 + v(r ) , v(r ) = 12 ^ 2 r^2 + 14 r^4 , [ hrad ] is integrable due to the spherical symmetry of the potential , so that we can apply the standard ebk ( or radial wkb ) approximation to it . writing @xmath85 in terms of polar coordinates @xmath86 and the associated canonical momenta @xmath87 , we have the usual form involving an effective potential @xmath88 that includes a centrifugal term : e = h(r , p_r , l ) = 12 p_r^2 + v_eff(r ) , v_eff(r ) = v(r ) + , [ heff ] where @xmath89 is the radial momentum . the three independent ( and poisson - commuting ) constants of the motion are the energy @xmath18 , the total angular momentum @xmath90 , and its @xmath91 component @xmath92 . the momentum @xmath93 can be expressed in terms of the latter two as p_= . [ pthet ] before we specialise to our particular potential @xmath94 in , we briefly recall the radial ebk quantisation and derive from it a general trace formula for an arbitrary spherical potential , starting from the hamiltonian .      in the standard radial ebk method @xcite , one quantises the three following action integrals s_r & = & p_rdr = 2(n_r+1/2 ) , n_r = 0,1,2,  [irad ] + i _ & = & p_d= ( n_+1/2 ) , n_= 0,1,2,  [ithet ] + i _ & = & l_z = m , m=0,1,2,  [iphi ] where the maslov index 1/2 in is correct only for smooth potentials . since @xmath93 in does not depend on the potential , the integral for @xmath95 can be done once for all and yields i_= l-|l_z| 0 , [ ithet1 ] so that the quantisation condition for @xmath96 is given by l = ( n_+|m|+1/2 ) = ( + 1/2 ) , = n_+|m|=0,1,2 ,  [ lquant ] in terms of a single angular momentum quantum number @xmath97 . the relation includes the so - called langer correction ; the quantised squared angular momentum @xmath98 agrees with the exact quantum - mechanical value @xmath99 in the limit of large @xmath97 . solving for @xmath89 , we can write the radial action as s_r(e , l ) = dr = 2(n_r+1/2 ) , [ irad1 ] showing that the quantised energies will only depend on the radial and angular momentum quantum numbers @xmath100 and @xmath97 ; they have , of course , the usual @xmath101-degeneracy @xmath102 , since @xmath103 . inverting the relation , we may rewrite the hamiltonian in the form e = ( s_r , l ) . [ ham ] inserting the right - hand side of and into , we obtain the ebk - quantised eigenenergies : e_n_r^ebk = ( 2(n_r+1/2),(+1/2 ) ) , , n_r = 0,1,2,  [eebk ] they can in general only be obtained by numerical iteration after doing the radial action integral over @xmath104 within the classical turning points .      before continuing , we simplify the situation by a scaling of the energy . in principle we have to vary the three parameters @xmath18 , @xmath105 , and @xmath1 to study the dynamics of our present system . however , we can introduce a scaling of the energy in such a way that the classical dynamics only depends on one single parameter , a dimensionless scaled energy @xmath106 . if we multiply by the factor @xmath107 , we can write the r.h.s .  in terms of scaled coordinates @xmath108 and momenta @xmath109 and a scaled time @xmath110 : q_i = x_i,_i = p_i , = t , [ scal ] so that the scaled energy @xmath106 becomes e = e = 12**^2+v(q ) = 12q^2+v(q)+ , v(q)=12q^2 + 14q^4 . [ escal ] where the dimensionless scaled angular momentum @xmath111 is given by l = l / s , s = . [ sscal ] in @xmath112 is the scaled radial variable and @xmath113 the scaled radial momentum , and the dot means the derivative with respect to the scaled time @xmath110 . @xmath114 is the action unit . this brings about a considerable simplification of the classical dynamics : we only need to vary one parameter @xmath106 ; at the end of our calculations we just have to remember that energies are measured in units of @xmath115 , angular momentum and actions in units of @xmath116 , times in units of @xmath117 , etc . in the following we shall give all quantities as functions of the dimensionless scaled variables @xmath106 and @xmath111 . ( other scaled dimensionless quantities such as actions , periods etc . will be denoted by lower - case letters . ) the scaling in is specific for our present hamiltonian . it can , however , easily be modified for ho potentials perturbed by arbitrary central potentials which are pure power laws in @xmath112 . the results of this subsection can therefore easily be generalized to the corresponding suitably scaled potentials @xmath118 . we now take the spectrum as a starting point to write down the density of states in the ebk approximation : g_ebk(e ) = _ n_r=0^_=0^ ( 2 + 1)(e - e_n_r^ebk ) . [ gebk ] next we apply poisson summation @xcite to convert the sums over @xmath100 and @xmath97 into integrals : g_ebk(e ) = _ n=-^_m=-^ _ 0^d ( 2 + 1 ) _ 0^dn_r ( e - e_n_r^ebk ) e^i2(nn_r+m ) +  [ gebk1 ] due to the finite lower limits of the summations in , there are boundary corrections to , which we have indicated by the dots . we shall comment on them in a moment . using and , we substitute the variables @xmath97 and @xmath100 by the classical actions @xmath96 and @xmath119 , using @xmath120 and @xmath121 . then @xmath122 becomes , expressed in terms of the dimensionless scaled variables , g_ebk(e ) = _ n=-^ _ m=-^_0^l_m(e ) dl2l _ 0^s_r(e , l ) ds_r(e-(s_r , l ) ) e^is[ns_r+2ml]/-i(n+m ) +  [ gebk2 ] hereby @xmath123 is the scaled radial action integral and @xmath124 the scaled hamiltonian . note that the integration limits have been imposed by the energy conservation ; @xmath125 is the maximum scaled angular momentum at fixed energy . the lower limits are , strictly speaking , @xmath126 for the integral over @xmath111 and @xmath127 for the integral over @xmath128 . replacing them by zero corresponds to neglecting corrections of higher order in @xmath77 . in the following , we keep terms up to order @xmath77 with respect to the leading factor @xmath129 in . there are exactly two corrections of relative order @xmath77 in what is indicated by the dots above : one is a boundary correction from @xmath130 to , and the second is coming from the lower limit @xmath131 of the @xmath128 integral in . a short calculation shows that these two corrections cancel identically . all other terms neglected in correspond to corrections of relative order @xmath42 or higher . we now use the relations ( e-(s_r , l ) ) = |s_r(e , l)|(s_r - s_r(e , l ) ) and s_r(e , l ) = = t_r(e , l ) . here @xmath132 is the scaled period of the classical motion at fixed values of @xmath106 and @xmath111 , which is the energy derivative of the corresponding scaled action integral @xmath133 : s_r(e , l ) & = & _ q(e , l)dq = 2_q_1^q_2 dq , [ srad ] + t_r(e , l ) & = & = 2_q_1^q_2 dq . [ trad ] here @xmath134 and @xmath135 are the scaled lower and upper turning points of the radial motion , respectively , which both depend on @xmath106 and @xmath111 . the above integrals can in many cases be expressed in terms of complete elliptic integrals . ( exceptions are the ho and coulomb potentials where they become simple algebraic functions of @xmath106 and @xmath111 . ) for our present potential @xmath118 in , the analytical expressions for , and other quantities of interest are given in the appendix [ app3 ] . one important result derived there is the taylor expansion of @xmath133 around @xmath136 , whose first terms are s_r(e , l ) = s_r(e,0 ) - l + a(e)l^2 + o(l^3 ) , [ srltay ] where @xmath137 is given in . the structure of  but not the explicit form of @xmath137  appears to be a general result valid for all regular central potentials @xmath118 with a minimum at @xmath138 ( and hence not for the coulomb potential , see appendix [ appcoul ] ) , but we were not able to prove it in the general case . the relevance of the linear term @xmath139 in will become clear in the next subsection .    using the above relations we can now do the integral over @xmath128 in exactly , due to the delta function , and obtain the following `` ebk trace integral '' : g_ebk(e ) = _ n=-^ _ m=-^(-1)^n+m _ 0^l_m(e ) dllt_r(e , l ) e^i[ns_r(e , l)+2msl]/ + o(^-1 ) . [ gebkint ]    the phase of the integrand of can be interpreted as the full action @xmath140 , divided by @xmath77 , of a given classical orbit labelled by @xmath141 and @xmath24 , consisting of the radial part @xmath142 and the angular part @xmath143 . when doing the full double summation over all @xmath141 and @xmath24 and the integration over @xmath111 exactly , yields the ebk spectrum to leading order in @xmath77 . note that the limits of the @xmath111 integral are the cases of zero angular momentum , which corresponds to the diameter orbits , and its maximum value @xmath125 at a given energy , which corresponds to the circle orbits . the latter have the radius @xmath144 at which the effective scaled potential @xmath145 has its minimum ; for this motion the radial action is zero : @xmath146 , since all energy is in the angular motion . all contributions with @xmath147 to the integral correspond to motion which has both radial and angular components . the formula , valid for arbitrary ( but correctly scaled ) spherical potentials @xmath118 with smooth walls , is in principle a trace formula , but it does not yet have its characteristic form . that form is obtained by evaluating the integral over @xmath111 in the semiclassical limit @xmath148 by evaluating its leading contributions from the critical points of the phase function @xmath149 in the exponent of the integrand . to leading order in @xmath77 , these are stationary points  as far as they exist . using the standard stationary - phase evaluation of the integral around the stationary points , one obtains contributions to @xmath122 with amplitudes proportional to @xmath150 , as expected for the fully three - fold degenerate orbits in a spherical system ( cf .  the discussion at the end of the previous section ) . next to leading order in @xmath77 , one obtains contributions from the end points of the integral ( see , e.g. , @xcite ) , sometimes referred to as `` edge corrections '' . as already announced above and shown explicitly below , these correspond here to the diameter and circular orbits , giving contributions with amplitudes proportional to @xmath151 . a special contribution to @xmath122 comes from @xmath152 . as generally proved by berry and mount @xcite , this must be the thomas - fermi ( tf ) value of the density of states , which is the leading contribution to its smooth part . from we obtain with @xmath152 g_ebk^(0)(e ) = _ 0^l_m(e)dl^2 t_r(e , l ) . [ tfint ] the integral of over @xmath153 is straightforward . note that the contribution from the upper limit @xmath125 in is zero , since the turning points coincide : @xmath154 . the result is g_ebk^(0)(e ) = _ 0^q_2 ^ 0q^2dq = g_tf(e ) . [ gtf ] here @xmath155 is the upper turning point for @xmath136 . ( the lower one is zero : @xmath156 . ) that really is equal to the tf density of states follows from its general definition g_tf(e ) = d^3p d^3r ( e- ( , ) ) = d^3 d^3q ( e- ( , * ) ) . using polar coordinates for @xmath157 and @xmath158 and doing the @xmath158 integration leads to . the analytical expression valid for the potential is given in . in the introduction we have stated that the average part of the density of states is generally given by the etf approximation which contains @xmath77 corrections to the tf limit . for the spherical three - dimensional ho , these corrections are of relative order @xmath42 , as seen in . for the present perturbed ho potential we will see that the tf approximation is sufficient to reproduce the average part of the quantum - mechanical density of states , at least up to the energies for which the quantum spectrum is numerically available . * the stationary condition for the phase function in at a point @xmath159 reads [ ns_r(e , l)+2ml]|_l_nm = n|_l_nm+2m= 0 . [ stat ] due to energy conservation we may write de = d(s_r , l ) = |_l ds_r + |_s_r dl= 0 , so that = -2 , [ dsrdl ] where the frequencies of angular and radial motion are defined as usual by _ ( e , l ) = , _ r(e , l ) = 2 . with this , the stationary condition becomes = = , [ tori ] which is the periodicity condition for the `` rational tori '' as the most degenerate classical orbits @xcite . clearly , @xmath24 and @xmath141 must have the same sign , since frequencies and periods are positive quantities . because of , the lower integration limit @xmath136 in is a stationary point corresponding to the diameter orbit with @xmath160 . however , the stationary - phase integration of for @xmath161 yields a zero contribution due to the factor @xmath111 in the integrand , so that the only semiclassical contribution due to the diameter orbits is the end - point correction for @xmath136 discussed in the next section . as we shall see further below , solutions of with @xmath162 do not exist for all energies @xmath106 and for all pairs @xmath163 of integers . therefore we postpone the contributions of the rational tori to and concentrate first on the diameter and circle orbits .    before developing the corresponding trace formula for our present system , valid to all orders in @xmath1 , we want to establish here the connection to the first - order perturbative approach used in sect . [ secpert ] . from the taylor expansions given in appendix [ app3 ] , we obtain for the radial action integral the following first terms : 2s_r(e , l ) = 2(-l ) -(3e^2-^2l^2 ) + o(^2 ) . [ srad1 ] in the second term we recognise precisely the first - order action shift given in . inserting into it @xmath164 and the leading - order expression for @xmath165 given in , we obtain the first - order action shifts of the diameter and circle orbits , respectively , given in . we now insert into the expression for the density of states , neglecting the terms of higher order in @xmath1 and keeping only the zero - order terms in the amplitude of @xmath166 , given in , and of the upper integration limit @xmath165 . noting furthermore that for the unperturbed ho orbits the ratio @xmath167 becomes equal to 2 , so that the resonance condition implies @xmath168 , we assume for the moment that all other combinations of @xmath24 and @xmath141 in this limit may be neglected . writing @xmath169 , we obtain the following result for the density of states up to first order in @xmath1 ( note that the linear terms in @xmath96 cancel in the exponent of the integrand ) : g_ebk^(1)(e ) = _ k=-^(-1)^k e^i2ke/ _ 0^e/ dl^2 e^-ik(3-^2l^2/e^2)]/ , [ gebkp ] where @xmath46 is the quantity defined in . using the substitution @xmath170 , the above integral becomes identically equal to the integral for the perturbative modulation factor , and the result is precisely the perturbed density of states defined in with the oscillating part given in . rather than proving at this stage that all contributions to with @xmath171 either cancel or are negligible to the leading order in @xmath77 , we now go on to derive the full trace formula for the diameter and circle orbits , doing the summations over all @xmath24 and @xmath141 exactly in the semiclassical limit . equipped with the above results , we are now in a position to evaluate the full trace formula for the contribution of diameter and circle orbits for our present hamiltonian . the smooth tf part of the level density , @xmath172 , is given explicitly in the appendix [ app3 ] . the contributions of the rational tori , which bifurcate from the circle orbits only for high enough energy and repetition numbers , will be discussed in section [ sectori ] . we therefore now evaluate the end - point contributions to the integral in in the semiclassical limit @xmath148 . the lower end point @xmath136 in yields the asymptotic contribution g_d(e ) = t_r(e,0 ) _ n=-^ ( -1)^n e^i[ns_r(e,0)/+/2]u_n(e ) , [ gd1 ] where we have left out the contribution from @xmath173 which contributes to the tf part , and defined the quantity u_n(e ) = _ l0 . using the identity @xcite = _ m=-^ , ( zn , n ) [ polesum ] we find u_n(e ) = _ l0 . [ une ] we now exploit the result from which we find = -+ 2a(e ) l + o(l^2 ) . for odd values of @xmath24 , the sin function in gives always a nonzero denominator in the limit @xmath136 , so that @xmath174 becomes zero . for even @xmath175 we get in the limit @xmath176 ( ) = ( -k+2kla(e ) ) = ( -1)^k(2kla(e ) ) ( -1)^k 2kla(e ) , so that we obtain u_2k(e ) = . inserting this into using , leaving out the contribution @xmath38 which is contained in @xmath177 , we obtain the final contribution of the diameter orbits to the density of states : g_d(e ) = -_k=1^ ( ks_d(e)/ ) , [ dgd2 ] where the actions @xmath178 are given explicitly in . in the low - energy limit we obtain with exactly the diameter contribution to the perturbative trace formula . note that the presence of the linear term @xmath179 in the quantity is instrumental in annihilating the contributions with odd @xmath24 to the sum in and hence in establishing the fact that only an even number of radial oscillations yields a physical periodic orbit with angular momentum @xmath164 . the upper end point @xmath180 in yields the asymptotic contribution g_c(e ) = l_m(e)t_r(e , l_m ) _ m=-^(-1)^m e^i[2l_m(e)/-/2 ] _ n=-^. [ gc1 ] using , and noting that @xmath181 , we get |_l_m = - . employing again to perform the @xmath24 summation , writing @xmath182 and omitting the @xmath38 term , we obtain after some manipulations the contribution of the circle orbits to the density of states g_c(e ) = _ k=1^ ( ks_c(e)/ ) , [ dgc2 ] where @xmath183 is given explicitly in , and the quantity @xmath184 is defined by w(e ) = = . [ wofe ] in the limit @xmath185 this quantity is found with the r.h.s .  of the analytical results and to go like w(e ) = 2(1+e/4 +  ) = 2+e +  expanding the denominator in , using the above @xmath184 , up to first order in @xmath1 brings the amplitude exactly into that of the circle orbit contribution to the perturbative result . note that the denominator under the sum in looks exactly like that in the gutzwiller trace formula @xcite for an isolated stable orbit , whereby @xmath186 corresponds to one - half of the stability angle . indeed , the semiclassical amplitude in diverges when @xmath187 becomes an integer @xmath188 . at the corresponding nonzero energies , the circle orbit bifurcates ; the condition for this to happen is exactly the resonance condition for the rational tori with @xmath189 , taken at @xmath190 . ( we come back to this point in section [ sectori ] . ) using and , we find that @xmath184 is restricted to the following range : 2 < w(e ) = < 0 . [ wlim ] the smallest @xmath36 for which we can have @xmath191 with @xmath192 is @xmath193 with @xmath194 , so that @xmath195 , which happens at the scaled energy @xmath196 ( see below ) . this means that the shortest new orbit is the 7:3 torus , bifurcating from the 3rd repetition of the circle orbit . the next bifurcation ( from @xmath197 ) is that of the 9:4 torus at @xmath198 . still longer orbits bifurcate at lower energies , but from higher repetitions of the circle orbit . therefore , we can ignore the bifurcations at low energies and small repetition numbers @xmath36 and concentrate on the contributions of the diameter and circle orbits to the density of states . let us finally write down the trace formula which we have obtained for the oscillating part of the density of states : g(e ) _ k=1^ , [ dguni ] where the semiclassical amplitudes are given by _ k^c(e ) = , _ k^d(e ) = . [ ampcd ] in the following we shall test this formula against the quantum - mechanical density of states . as we have shown above , goes over into the perturbative trace formula in the limit @xmath26 . this confirms the assumption , made at the end of sect . [ secpois ] , that the contributions with @xmath171 to do not contribute in this limit .      in this section we compare our semiclassical results with those obtained from the exact quantum - mechanical spectrum ( obtained numerically by solving the radial schrdinger equation on a discrete mesh ) . in order to focus on the _ coarse - grained _ shell structure , we convolute both results over the energy @xmath18 with a normalised gaussian of width @xmath199 . the quantum - mechanical density of states then becomes g_qm(e ) = _ n e^-(e - e_n)^2/^2 . [ gqmgam ] in order to obtain its oscillating part , we subtract from it the tf expression @xmath172 which we can calculate analytically ( see appendix [ app3 ] ) . the semiclassical trace formula becomes , after doing the convolution in stationary - phase approximation ( cf . @xcite , sec . 5.5 )    g_sc(e ) _ po a_po(e ) e^-[t_po(e)/2]^2 [ s_po(e)/-_po/2 ] , [ dgscgam ] so that orbits with longer periods @xmath200 will be exponentially suppressed .    in we show a comparison of the results obtained with the uniform trace formula for the case @xmath201 with those obtained from the exact numerical quantum spectrum . the width of the gaussian smoothing was chosen to be @xmath202 ; in the semiclassical result the harmonics @xmath203 then do not contribute noticeably . the agreement is seen to be perfect ; tiny differences can only be noted near the beat minima on the amplified scale below .    in we show a similar comparison for @xmath204 , exhibiting much fine structures . we now start to resolve the spectrum with a higher resolution ; note that we still only plot the oscillating part @xmath7 of the density of states . in the top panel , the semiclassical result includes harmonics ( i.e. , repetitions of the two orbits ) with @xmath205 . this is evidently sufficient to reproduce the quantum - mechanical oscillations up to @xmath206 with a high accuracy . in the centre panel , we have added two more harmonics , going up to @xmath207 . here we recognise the appearance of divergences at the scaled energies @xmath208 and @xmath209 which correspond to the bifurcations of the 19:9 and 17:8 tori from the repetitions of the circle orbit with @xmath210 and @xmath211 , respectively ( cf .   below ) . the small wiggles close to the divergences , with decreasing amplitudes when going away from them , are due to the missing contributions from those torus orbits in the trace formula ( see the following section ) . = 13.5 cm    in order to anticipate the results of a suitable uniform trace formula including the contributions of the tori , we show in the bottom panel of the result of the ebk trace integral obtained by a numerical integration over the variable @xmath111 . hereby the summations over @xmath24 and @xmath141 have been done in the limits @xmath212 , @xmath213 using @xmath214 . [ the period to be used in the gaussian damping factor of is in this case given by @xmath215 . ] now the divergences and all other discrepancies have disappeared ; the difference to the quantum - mechanical result can hardly be recognised . the integral therefore is also a uniform expression for the semiclassical density of states , containing the contributions of all periodic orbits . the numerical evaluation becomes , however , very time consuming for even finer energy resolutions ( i.e. , still smaller values of @xmath199 ) . our main result here is that up to a rather fine resolution , the shell structure of the present perturbed ho system is dominated by the two - fold degenerate families of diameter and circle orbits , although these are not the most degenerate orbits . this is a rather unusual situation , due to the fact that the shortest tori with the highest degeneracy are considerably longer ( by a factor @xmath21610 ) than the primitive diameter and circle orbits . let us finally test the range of validity of the perturbative trace formula . in we compare its results to that of the full trace formula for the two values @xmath217 ( upper panel ) and @xmath201 ( lower panel ) . as expected , the two curves agree in the limit @xmath185 . however , already for @xmath217 the amplitude modulation of the perturbative result deviates so much that the position of the first beat node is predicted by about 10% too low in energy . for @xmath201 , the perturbative result becomes so bad that it predicts the second beat node approximately at the position of the correct first one . these results underline the fact that in @xcite , the quantum results for weak interactions could be qualitatively reproduced by the perturbative trace formula , but not quantitatively with correct positions of the beat nodes . a more detailed analysis of the selfconsistent hf results of @xcite using the uniform trace formula is in progress @xcite .      as we have seen in the previous section , periodic orbit with both angular and radial motion appear as the rational tori fulfilling the periodicity condition . for the hamiltonian , this condition can only be fulfilled in the limits given in . this means that stationary points @xmath218 of the @xmath96 integral in satisfying and hence only exist if @xmath219 . hereby @xmath220 corresponds to the repetition number @xmath36 of the circle orbit that undergoes a bifurcation exactly when @xmath218 is equal to the maximum angular momentum @xmath221 . in we list some of the lowest energies @xmath222 at which this happens for the smallest integers @xmath189 . each of the @xmath223 tori only exists above its bifurcation energy , i.e. , at energies @xmath224 . [ cols="^,^,<,<,<",options="header " , ]     the contributions of the bifurcated tori to the semiclassical trace formula is obtained by a stationary - phase evaluation of the @xmath111 integral in , leading to g_tori(e ) = _ n>2m>0 a^t_nm(e)(-1)^n+m ( + ) . [ dgt ] the summation goes only over those pairs @xmath163 of positive integers for which the stationary condition can be fulfilled . ( note that the contribution of the corresponding pairs of negative integers is included in the amplitudes by a factor 2 . ) the actions of the tori are given by s^t_nm(e ) = ns_r(e , l_nm)+m2l_nm(e ) , [ stori ] and their amplitudes by ^t_nm(e ) = t_r(e , l_nm)l_nm(e ) ( n| |_l = l_nm(e))^-1/2(e - e_n : m ) , [ amptori ] where @xmath114 is the action unit given in . the heaviside step function in is defined as ( x ) = 0 x < 0,(x ) = 1/2 x = 0,(x ) = 1 x > 0,where the value 1/2 for @xmath225 accounts for the fact that one obtains only half a fresnel integral when the stationary point is the upper end point of the integral in . the power @xmath150 in the amplitude @xmath226 is characteristic of the three - fold degenerate orbits in a three - dimensional spherical system @xcite . .   corresponds to the standard berry - tabor trace formula @xcite for the most degenerate orbits in an integrable system . here , however , each contribution to is only valid if the energy @xmath106 is sufficiently larger than the corresponding bifurcation energy @xmath222 . in the neighbourhood of the bifurcation energies , we have to replace the sum of the diverging circle orbit contribution and the corresponding torus contribution with @xmath227 by a common uniform contribution . as it is well - known from the semiclassical theory of bifurcations @xcite , one has to include hereby also the contributions of so - called `` ghost orbits '' , which are the imaginary ( or complex ) continuations of the bifurcated orbits ( here : the tori ) to the other side of the bifurcation ( here : @xmath228 ) . it is , however , not the purpose of our present paper to discuss in more detail the uniform approximation suitable for this situation . the present scenario is actually almost identical to that described in @xcite for the bifurcations of tori in the integrable hnon - heiles system . although there the orbit from which the tori bifurcate is an isolated one , while it here is a family of circle orbits , the uniform trace formula given in @xcite can be applied to the present system in a straightforward manner . we refer the interested reader to this paper for all the detailed formulae as well as for a numerical demonstration showing how diverging discrepancies of the type seen in the centre part of disappear when a proper uniform trace formula is used . its result would here be practically indistinguishable from the result shown in the bottom part of , in which we have evaluated the integral in numerically instead of using the stationary phase approximation with end - point corrections . in this section we discuss briefly a purely quartic oscillator potential , i.e. , we start from the hamiltonian h ( , ) = 12p^2 + r^4 = e. [ hr4 ] here both the energy and @xmath1 can be scaled away . after dividing the above equation by @xmath18 and introducing the scaled variables q_i = ( ) ^1/4r_i , = ( e)^1/4t , l = ^1/4e^-3/4l , we obtain the `` unit energy '' equation 1 = 12q^2 + 14q^4 + , where the dot again means derivative with respect to the scaled time @xmath110 . hence the classical equations of motion become independent of energy and @xmath1 . after obtaining all the classical results , we can reintroduce @xmath18 and @xmath1 just by remembering that lengths scale with @xmath229 , momenta with @xmath230 , times ( periods ) with @xmath231 , actions and angular momenta with @xmath232 . we obtain the semiclassical trace formula from the results of the previous sections and of the appendix [ app3 ] simply by taking the limit @xmath30 , which means practically by extracting the leading terms for @xmath233 . in this way , we obtain the following primitive periods and actions for the circle orbit : t_c(e ) = 2 ^ -1/4(4e/3)^-1/4,s_c(e ) = 2l_m(e ) = 2 ^ -1/4(4e/3)^3/4 , [ circr4 ] and for the diameter orbit : t_d(e ) = 2 ^ -1/4k_0e^-1/4,s_d(e ) = 2 ^ -1/4(4/3)k_0e^3/4 . [ diagr4 ] here @xmath234 is the constant elliptic integral k_0=k(_0)=1.8540746773  ,_0 = 1/. furthermore we get t_r(e , l_m ) = ( 3e)^-1/4 . the resonance condition becomes independent of the energy = w(l_nm ) = , [ wofl ] since both periods scale with the same power of the energy . for the frequency ratio we find the limits 2 w(l_nm ) = n : m < 0 l_nm < l_m = ( 4/3)^3/4 = 1.240806  [tori4 ] note that the upper limit @xmath235 never becomes a stationary point , so that no bifurcations of the circle orbits occur in this system . the values of the stationary points @xmath236 ( in scaled dimensionless units ) for the lowest values of @xmath24 and @xmath141 are given in the rightmost column of ; as we just have seen , they do not depend on the energy @xmath18 . the result @xmath237 is a special case of the general result , given in @xcite , that for a potential @xmath94 = @xmath238 one has @xmath239 .    for the amplitudes of the circle and diameter orbits in the trace formula , which keeps the same form , we obtain here _ k^c(e ) = , _ k^d(e ) = . [ ampr4 ] the actions @xmath83 and @xmath82 to be used in are those given in and . finally , the tf level density is found to be g_tf(e ) = . the tori , which in the perturbed ho system discussed in the last section occur only at sufficiently high energies after their bifurcations from the repeated circle orbits , exist in the present system at all energies . they affect again only the finer quantum structures of the density of states since they only exist for @xmath240 due to the selection by their resonance condition ( cf .  also ) . this is demonstrated in the following two figures .    in , we compare the results of the semiclassical trace formula for the diameters and circles using the above quantities to the exact quantum - mechanical density of states , both gaussian - averaged over a width @xmath241 ( in energy units ) . the agreement is again perfect even on the enlarged scale in the lower part of the figure . the reason is that for the chosen energy resolution with @xmath241 , only repetition numbers @xmath242 contribute , for which there exist no tori ( see ) . = 9.2 cm    = 4.5 cm    shows the same comparison for @xmath243 . here some small discrepancies appear which become more visible with increasing energy . they are due to the missing contributions from the tori with @xmath244 . since there are no bifurcations involved in the present system , the contributions of the tori are given by the berry - tabor trace formula , employing the appropriate expressions for the actions and periods used therein . the analytical expression of the action integral @xmath245 in terms of elliptic integrals is the same as for the perturbed ho potential , as given in eqs .   - , except that here the term @xmath246 in the potential and the corresponding contribution @xmath247 must be omitted .    even at the rather fine resolution obtained in , our main result is the observation that the circle and diameter orbits dominate the shell structure also for the purely quartic oscillator . motivated by recent hartree - fock ( hf ) calculations of harmonically trapped fermionic gases with a repulsive interaction @xcite , in which pronounced super - shell effects were obtained , we have developed a semiclassical trace formula for the density of states of the anharmonically perturbed three - dimensional isotropic harmonic oscillator ( ho ) defined in . we find that its gross - shell structure is dominated by the periodic orbits of diameter and circle shape , whose interference leads to super - shell beats which explain the findings from the numerical hf results .    in a first step , we have used the perturbative approach of creagh @xcite to describe the symmetry breaking u(3 ) @xmath25 so(3 ) for weak anharmonicities @xmath1 , resulting in the perturbative trace formula . it uniformly restores the u(3 ) limit , yielding the exact trace formula of the unperturbed ho system in the limit @xmath26 . for the derivation of the perturbative result , one must in principle integrate over the 4-fold degenerate families of unperturbed ho orbits which cover the manifold @xmath2p@xmath3 . as shown in app . [ app2 ] , however , it turns out to be easier to integrate over the energy shell in phase space , which is a five - sphere s@xmath31 , whereby the integration can be reduced , under exploitation of the so(3 ) symmetry , to the one - dimensional integral given analytically in . an interesting result hereby is the fact that its two contributions , corresponding to the diameter and circle orbits , already have the characteristic form occurring in the general trace formula . usually , this form is obtained from a perturbative trace formula only asymptotically in the semiclassical limit @xmath148 ; this happens , e.g. , also when one investigates the corresponding hamiltonian in two dimensions or other perturbed 2-dimensional systems @xcite . next we have used the standard ebk quantisation and the poisson summation formula to derive a general semiclassical trace formula for arbitrary central potentials in terms of a one - dimensional integral over the angular momentum @xmath96 of the classical orbits , given in . exact integration and summation over all @xmath248 yields , to leading order in @xmath77 , the ebk spectrum . its asymptotic evaluation in the limit @xmath148 yields the standard - type contributions to the trace formula . the end - point contributions yield the diameter orbits with @xmath164 and the circle orbits with maximum angular momentum ; both contributions are of next - to - leading order in @xmath77 , whereas the leading - order contributions corresponding to the typical rational tori of integrable systems according to the berry - tabor theory @xcite come from the stationary points inside the integration interval . one interesting mathematical aspect is the mechanism by which the ratio @xmath249 of radial to angular frequency of the diameter orbit is naturally selected by a property of the general radial action integral @xmath250 as a function of angular momentum @xmath96 which appears to be universal for regular spherical potentials . it holds also for a particle in a spherical box with specular reflection , for which we have re - derived in appendix [ appbb ] the trace formula given by balian and bloch @xcite from our general formula , taking into account the changes of the maslov indices due to hard - wall reflections ) . the coulomb potential , for which @xmath251 , is discussed in appendix [ appcoul ] where we also re - derive the exact trace formula for the rydberg spectrum given in @xcite .    for the hamiltonian , we have derived a uniform trace formula for the diameter and circle orbit contributions that is valid for arbitrary @xmath1 . we find that for low energies and low repetition numbers no tori exist , so that the gross - shell structure of this system is dominated by the circle and diameter orbits and their lowest repetitions . the tori bifurcate , at sufficiently high energy , out of the repetitions with @xmath78 of the circle orbits , as was also observed recently for homogeneous power - law potentials in @xcite . their contributions sufficiently high above the respective bifurcation energies are given by the usual berry - tabor type trace formula . we have not discussed here the common uniform treatment of the circle orbits and the tori bifurcating from them , since this scenario is identical to the one discussed in detail in @xcite . the same qualitative results are found also in the limit @xmath252 in which the potential becomes a purely quartic oscillator . although the tori with @xmath80 here exist at all energies , they are so much longer ( by a factor @xmath253 ) than the shortest diameter and circle orbits that they only affect the finer details of the quantum spectrum .    we thus have found the interesting and quite atypical situation of a three - dimensional system in which the periodic orbits with highest , i.e. , three - fold degeneracy are only responsible for finer details of the quantum spectrum , whereas its gross - shell structure is to an astonishing degree dominated by the orbits of next - to - leading order in @xmath77 , i.e. , the diameter and circle orbits occurring in two - fold degenerate families . this is totally different from the situation , observed first by balian and bloch @xcite , of a spherical cavity potential ( with ideally reflecting walls ) which has turned out to be a realistic model for large spherical metal clusters . ( we re - derive the semiclassical trace formula of balian and bloch for the spherical cavity in appendix [ appbb ] . ) there the famous super - shells @xcite come from the interference of the shortest three - fold degenerate tori ( the triangle and square orbits ) , whereas the two - fold degenerate diameter orbit family has virtually no effect on the observable shell structure of these systems . our results explain qualitatively the numerical quantum - mechanical hf results @xcite for harmonically trapped fermionic atom gases . a more quantitative comparison , in which the strength @xmath1 of the perturbation in is determined directly from the numerically obtained self - consistent hf potentials , is in progress @xcite . lazzari _ et al . _ @xcite have studied the periodic orbits in a spherical woods - saxon potential in connection with the physics of metal clusters and thereby focused on the contributions of the diameter orbits . we find agreement of their results for the tori and the diameters with our results , and , respectively , if we interpret their function @xmath254 as being _ the radial period @xmath255 at the corresponding stationary value @xmath218 of the angular momentum (= 0 for the diameters ) . they have , however , neglected the contributions from the circle orbits ; furthermore they did not compare their semiclassical results with quantum - mechanical results . we suspect that the neglect of the circle orbits might be justifiable due to the chosen steepness of their woods - saxon potential , for which the triangle - like tori are perhaps more important than the circles . ozorio de almeida _ et al . _ @xcite have discussed the summation of all periodic orbits in integrable systems in connection with spectral determinants . they conclude that edge corrections corresponding to lower - degenerate orbits may be neglected because they are implicitly approximated by the nearest - lying rational tori ( characterised by very large numbers @xmath248 in our notation ) . while their argument applies to the _ full quantisation _ of the exact ( or ebk - approximated ) quantum spectrum obtained by summing over _ all orbits _ , it can not be used when the spectrum is coarse - grained with a finite energy window @xmath199 , such as we have used the periodic - orbit sum in our present investigation ( cf .  also the discussion in the appendix of ref . @xcite ) which is more practically oriented . we acknowledge enlightening discussions with s. bechtluft - sachs , s. creagh , and n. sndergaard , and helpful comments by an anonymous referee . we are particularly grateful to k. jnich for providing us with the elegant reduction of the s@xmath31 integration given in appendix [ app2 ] . we thank a. magner for valuable comments on our manuscript and for bringing the refs . @xcite to our attention . this work was financially supported by the swedish foundation for strategic research and the swedish research council . m. .  acknowledges financial support from the deutsche forschungsgesellschaft ( graduiertenkolleg 638 : _ nichtlinearitt und nichtgleichgewicht in kondensierter materie _ ) , and m.b . acknowledges the warm hospitality at the lth during several research visits . we reproduce here the parametrisation of @xmath2p@xmath3 given in @xcite . we start with three complex numbers @xmath256 : ( z_1,z_2,z_3 ) = ( n_1,n_2e^i_2,n_3e^i_3 ) , with @xmath257 , @xmath258 , and @xmath259 , @xmath260 , whereby the @xmath261 are restricted to n_1 ^ 2+n_2 ^ 2+n_3 ^ 2=1 . [ ensh ] the @xmath262 define the complex projective space @xmath2p@xmath3 . the @xmath261 are in one - to - one correspondence with the points on the first octant of the 2-sphere @xmath60 , and the @xmath263 form a 2-torus . at the edges of the octant , the torus contracts to a circle or to a point , and one or both of the @xmath263 are not defined . physically , we interpret the real and negative imaginary parts of the @xmath262 as coordinates and momenta , respectively , of a starting point @xmath264 r_(0 ) = e z_,p_(0 ) = - m z _ , in the 6-dimensional phase space @xmath265 of the 3-dimensional isotropic harmonic oscillator . the hamiltonian , which has u(3 ) symmetry @xcite , is given by h ( , ) = 12 _ = 1 ^ 3 = e. the @xmath266-periodic solutions of the equations of motion are then given by r_(t ) = e ( z_e^it ) , p_(t ) = -m ( z_e^it ) . the equation defines the energy shell , so that the six quantities @xmath267 are points on a 5-sphere , s@xmath31 ( see also appendix [ app2 ] below ) . by taking the total energy to be @xmath268 and setting @xmath269 , we fix the radius of the 5-sphere to be unity . taking out the phase factor @xmath270 takes us from s@xmath31 to the projective space @xmath2p@xmath3 .    using the fubini - study metric , the authors of @xcite find the squared distance @xmath271 between two points on @xmath2p@xmath3 to be ds^2 = dn_1 ^ 2+dn_2 ^ 2+dn_3 ^ 2 + n_2 ^ 2(1-n_2 ^ 2)d_2 ^ 2+n_3 ^ 2(1-n_3 ^ 2)d_3 ^ 2 - 2n_2 ^ 2n_3 ^ 2d_2d_3 . [ ds2 ] the first three terms correspond to the standard metric on @xmath60 , and the last three to the metric on a flat 2-torus whose shape depends on where we are on the octant . we now wish to re - parametrise the @xmath261 in terms of the two polar angles @xmath272 on @xmath60 by writing n_1 = , n_2 = , n_3 = , [ coord ] with @xmath273 $ ] . collecting our four angles in a vector @xmath274 by defining = ( _ 1,_2,_3,_4 ) = ( , , _ 2,_3 ) , we can write the squared distance @xmath271 as ds^2 = _ i , j=1 ^ 4 g_ijd_id_j , [ ds2a ] and find the elements @xmath275 of the metric tensor to be @xmath276 , @xmath277 , @xmath278 , @xmath279 , @xmath280 , and @xmath281 . from this , we get the integration measure on @xmath2p@xmath3 to be d_^2 = |det(g_ij)|^1/2 d_1d_2d_3d_4 = ^3d dd_2d_3 . the integrated volume of @xmath2p@xmath3 , consistent with a general formula for @xmath2p@xmath282 , is _ ^2 = d_^2 = _ ^n = . here we give a strict proof of eq .   of the modulation factor @xmath283 for the anharmonically perturbed ho system . the most straightforward parametrisation of the six - dimensional phase space of the unperturbed ho is given in terms of the coordinate vector @xmath284 and the momentum vector @xmath285 @xmath286 . since a periodic solution is uniquely determined by their initial values at @xmath287 : @xmath288 and @xmath289 , and the total energy @xmath290 is a constant of motion , the six components of @xmath291 and @xmath292 must cover a 5-sphere with radius @xmath293 . we introduce dimensionless coordinate and momentum vectors @xmath294 and @xmath295 , respectively , by = _ 0 , = _ 0 , [ unvec ] so that in these variables the energy shell becomes the unit sphere s@xmath31 : ^2+^2 = 1 . [ ens5 ] from eq . , the first - order action perturbation @xmath296 is given in terms of energy @xmath18 and angular momentum @xmath96 by s_1 = -(3-^2l^2/e^2 ) , [ ds1 ] with @xmath46 defined in . the modulation factor for the first - order perturbed ho system becomes therefore ( k/ ) = d_s^5 e^iks_1/ , [ modfac5 ] where @xmath297 is the five - dimensional solid angle and @xmath298 is the integration measure on s@xmath31 with @xmath299 . one may now express @xmath296 in directly in terms of the five polar angles of six - dimensional hyperspherical coordinates , which becomes a very complicated function , and integrate over all five angles . we have done this numerically and verified that it yields exactly the same result as the numerical 4-dimensional @xmath2p@xmath3 integration of with the phase function . it is , however , not necessary to perform the full five - dimensional integral for . we follow a much more economical route @xcite , exploiting the rotational symmetry of the system and the so(3 ) invariance of the expression for @xmath296 . due to the restriction , we can write the s@xmath31 sphere as s^5 = \{(*e*_,*e * _ ) | ; _ , * e*_s^2 } , [ s5new ] where @xmath300 and @xmath301 are the unit vectors in the directions of @xmath294 and @xmath295 , respectively . the square of the conserved total angular momentum then is l^2 = * l*^2 = ( _ 0_0)^2 = ^2 ^ 2 ( * e*__)^2 , so that the action perturbation becomes , after the substitution @xmath302 , s_1 = - , . [ ds1a ] the integration measure for is d_s^5 = ^2d_s^2(*e * _ ) ^2d_s^2(*e*_)d = 18 ^ 2d d_s^2(*e * _ ) d_s^2(*e * _ ) , and the modulation factor becomes ( k/ ) = _ 0^d^2 d_s^2(*e * _ ) d_s^2(*e * _ ) e^-ik/. [ mf5a ] we now introduce the angle @xmath303 between the unit vectors @xmath300 and @xmath301 , so that ( * e*__)^2 = ^2 , and choose @xmath300 as the direction of the north pole ( @xmath304 ) of polar angles @xmath305 for the vector @xmath301 . the integrand then becomes independent of @xmath300 , so that the corresponding s@xmath3 integral just gives @xmath306 . for the other s@xmath3 integral we have @xmath307 since the integrand does not depend on @xmath308 . we thus obtain , after the standard substitution @xmath309 and an obvious reduction of the integration limits , ( k/ ) = _ 0^/2 d^2 _ 0 ^ 1 du e^-ik/. next , we make the substitution @xmath310 and then a further substitution @xmath311 to obtain ( k/ ) = _ 0^/2d _ 0^ dt e^-ik(3-^2+t^2)/ = _ 0 ^ 1 ds _ 0^ dt e^-ik(2+s^2+t^2)/. since the last integral goes exactly over the first quadrant of a unit disk in the @xmath312 plane , we can use polar coordinates @xmath313 and furthermore @xmath314 to obtain our final result ( k/ ) = _ 0^/2 d _ 0 ^ 1 drre^-ik(2+r^2)/ = _ 0 ^ 1 dze^-ik(3-z)/ , which is identical to that given in or in . in this appendix we derive analytical expressions for the particular scaled potential @xmath118 in . rewriting the radial action integral using the substitution @xmath315 , we obtain s_r(e , l ) = 2_q_1^q_2 dq = _ z_1^z_2 , [ srint ] where g(z ) = 4ez-2z^2-z^3 - 2l^2 = ( z_1-z)(z_2-z)(z_3-z ) . [ gz ] the classical turning points are defined by the roots of the cubic equation @xmath316 , which are always real : the discriminant can easily be shown to be negative except for @xmath180 where it becomes zero due to the double root @xmath317 . selecting the roots such that z_3 < z_1 z_2 , one finds that @xmath318 is always negative , whereas @xmath319 and @xmath320 are positive definite and represent the squares of the real classical turning points @xmath134 and @xmath135 . the integrals in can not be found in most tables . however , after an integration by parts we may write the r.h.s . as s_r(e , l ) = s_r^(1)(e , l)+ s_r^(2)(e , l)-l^2 s_r^(-1)(e , l ) , [ srltot ] where we have defined s_r^(n)(e , l ) = _ z_1^z_2 dz . in byrd and friedman @xcite , these integrals are found to be : s_r^(1)(e , l ) & = & [ ( ^2-^2)k()+^2e ( ) ] , + s_r^(2)(e , l ) & = & [ ( 3 ^ 4 - 6 ^ 2 ^ 2 + 2 ^ 4+^2 ^ 4 ) k()+2 ^ 2(3 ^ 2-^2-^2 ^ 2 ) e ( ) ] , + s_r^(-1)(e , l ) & = & ( ^2 , ) . [ sradn ] here @xmath321 , @xmath322 and @xmath323 are the complete elliptic integrals of first , second and third kind with modulus @xmath324 , and = , ^2 = 1- . [ kap ] note that @xmath325 , which corresponds to the circular case of the elliptic integral @xmath323 @xcite . this function becomes singular like @xmath326 , i.e. , for @xmath327 which happens when @xmath176 . since @xmath319 goes to zero like @xmath328 in this limit , @xmath323 has a first - order pole at @xmath136 . using the laurent expansion of @xmath323 given in eq . 906.04 of @xcite , we find s_r^(-1)(e , l ) = -2a(e ) + o(l ) = -2a(e ) + o(l ) , [ srlm1 ] with a(e ) = - = ( 1-e +  ) . [ aofe ] here @xmath329 and @xmath330 are the @xmath136 values of the quantities given in eqs .   and below . the last equality in follows because @xmath331 , as seen from . the functions @xmath332 with @xmath333 are both regular in @xmath136 . altogether we obtain the following taylor expansion for the total radial action integral : s_r(e , l ) = s_r(e,0 ) - l + a(e)l^2 + o(l^3 ) , [ srlexp ]    the radial period is easily found to be t_r(e , l ) = s_r^(0)(e , l ) = k ( ) . [ trad1 ]    in the limit @xmath136 , which is relevant for the diameter orbits and the tf density of states , the above results simplify considerably . we obtain s_r(e,0 ) & = & 23(1 + 4e)^1/4 , [ sr0 ] + t_r(e,0 ) & = & 2(1 + 4e)^-1/4k(_0 ) [ tr0 ] . note that in this limit we get z_1 ^ 0=0,z_2 ^ 0=(q_2 ^ 0)^2=-1,z_3 ^ 0=--1 , [ zi0 ] and the modulus of the elliptic integrals becomes _ 0 = = = ( 1-e +  ) . [ kap0 ] the tf expression for the density of states becomes g_tf(e ) & = & ( 1 + 4e)^1/4 [ 2(1 + 3e)e(_0)-(1 + 3e+)k(_0 ) ] . [ gtf0 ] it is of interest to see the first terms in the taylor expansion of the above results around @xmath334 , whose leading terms are those of the pure harmonic oscillator : s_r(e , l ) & = & ( e - l)-18(3e^2-l^2)+  , [ srtay ] + t_r(e , l ) & = & -34e+(7e^2-l^2)+  , [ trtay ] + g_tf(e ) & = & e^2(1 - 5e/4 +  ) .[gtftay ] the unscaled leading ho terms are s_r^ho(e , l ) = ( -l ) = , t_r^ho = , g_tf^ho(e ) = . the diameter orbit with @xmath164 makes two full radial oscillations during its primitive period . therefore its primitive action and period are given by s_d(e ) = 2s_r(e,0),t_d(e)=2t_r(e,0 ) . [ stdia ] for the circle orbit , the primitive action integral @xmath183 is just @xmath266 times the maximum value of the angular momentum , @xmath335 , which for the present potential is found to be l_m(e ) = s , s_c(e ) = 2l_m(e ) = ( 1-e/4 +  ) . [ scirc ] the primitive period of the circle orbit is given by its energy derivative t_c(e ) = = ( 1-e/2 +  ) . [ tcirc ] we also give here the value of the radial period @xmath336 at @xmath180 . the corresponding frequency @xmath337 is easily obtained from the second derivative of the effective potential at its minimum @xmath144 , which for the present potential is given by q_0 ^ 2=z_0=2(-1)/3 , using @xmath338 , and the result becomes t_r(e , l_m ) = = ( 1 + 3e)^-1/4 = ( 1 - 3e/4 +  ) . [ trlm ] this agrees with the result obtained from , noting that for the maximum value of @xmath180 , @xmath339 and @xmath340 , so that @xmath324 becomes zero and @xmath341 . in order to illustrate the validity and the usage of our general ebk trace integral , we present in this appendix briefly its application to two popular spherical potentials . we consider a particle with mass @xmath19 in a spherical box with radius @xmath58 and ideal specular reflection from the boundary , i.e. , a three - dimensional spherical billiard . quantum - mechanically , one obtains the spectrum by solving the schrdinger equation with dirichlet boundary conditions . the semiclassical trace formula for this system has been derived by balian and bloch @xcite using a multi - reflection expansion of the green function . to derive it from our formula , we need the radial action integral for arbitrary energy @xmath342 and angular momentum @xmath343 . working throughout with unscaled variables , we find ( with @xmath344 ) s_r(e , l ) = 2_r_1^r dr & = & 2 - 2l(l / r ) + & = & 2r - l + l^2/r + o(l^3 ) , [ srlbox ] thus fulfilling our general relation with a(e ) = 1/r . with this , the stationary condition gives 2n(l / r ) = 2 m with the solutions l_nm = r(m / n ) . as in the potentials discussed in sect . [ secuni ] , the solution with @xmath164 corresponding to n=2 m contributes only via the end - point correction of to the diameter orbit with the primitive action s_d(e ) = 2s_r(e,0 ) = 4r and , according to the lower equation in , with the amplitude of its @xmath36-th repetition _ k^d(e ) = . ( k=|m| ) the upper end - point correction with @xmath345 can , for @xmath346 , only be reached formally in the limit @xmath347 . it corresponds to the `` whispering gallery mode '' with amplitude @xmath348 , cf . below , and does therefore not contribute to the trace formula . the tori with @xmath349 have the actions s^t_nm(e ) = ns_r(e , l_nm)+2 m l_nm = 2nr_nm , _ nm = m / n . [ stbb ] this corresponds exactly to the actions of the polygonal orbits given in @xcite with winding number @xmath350 and @xmath351 corners ( i.e. , reflections from the boundary ) . the case @xmath352 corresponds to the @xmath36-th repetition of the diameter orbit . for the polygons , @xmath353 is half the polar angle covered by one of their segments . the amplitudes of the tori with @xmath354 become , according to , ^t_vk(e ) = ( 2_vk ) . [ ampbb ]    before we can write down the trace formula with the correct phases , we have to correct the radial quantisation condition , because the maslov index changes from 1/4 to 1/2 per turning point if a reflection happens there . hence the quantisation condition is @xmath355 for the diameter orbit and @xmath356 for the tori . this changes the phases in and all the results derived from it . taking this into account , we obtain the semiclassical trace formula for the spherical billiard g(e ) = _ v>2k>0^t_vk(e ) - _ k>0 a_k^d(e ) , [ dgbb ] which is identically the result of balian and bloch @xcite .      here we derive the trace formula for the coulomb potential as a special central potential which is not regular in @xmath357 : v(r ) = - . [ vcoul ] we need not use scaled variables here , but it is useful to introduce the positive energy @xmath358 since all bound states in the potential have @xmath359 . the maximum angular momentum and the action of the circle orbit are l_m(e ) = , s_c(e)= , from which the period of the circle orbit is found to be t_c(e ) = 2 = . the radial action integral becomes elementary ( cf . @xcite ) : s_r(e , l ) = _ r_1^r_2 dr = 2(-l ) = 2 . [ srcoul ] from it we find the radial period t_r(e ) = . the difference from the other potentials treated in this paper is that the term linear in @xmath96 in here has the coefficient @xmath360 . furthermore , we find from the general relation = -2 _ ( e ) = _ r(e ) , t_r(e ) = t_(e ) . [ trat ] thus , the radial and angular frequencies are identically the same for all values of the angular momentum @xmath96 . this is a special property of kepler s ellipses . [ note that here , in contrast to the spherical ho , the angular momentum @xmath96 is defined with respect to one of the focal points and not to the symmetry centre of the ellipse ! ]    inserting the above results into the general ebk trace integral for the density of states , we obtain g_coul(e ) = _ m , n=-^ ( -1)^m+n_0^l_m(e ) ldl e^i2/. [ gcoul ] due to , the general resonance condition implies @xmath361 for all periodic orbits . we therefore only use the terms with @xmath362 in the double sum of ( cf .  the discussion at the end of sect .   [ secpois ] ) . we then find g_coul(e ) = . writing this in terms of atomic units ( note that we have put @xmath19 ) e_at= = ry , we finally obtain the quantum - mechanically exact trace formula g_coul(e ) = _ n=1^n^2 ( e+e_at / n^2 ) = , [ gcoul1 ] which has been given in @xcite . the eigenenergies e_n = -,(n = n_r++1 ) with degeneracies @xmath363 are , as is well known , also found directly from the ebk quantisation conditions - using and . m. c. gutzwiller : _ chaos in classical and quantum mechanics _ ( springer verlag , new york , 1990 ) ; + h .- j . stckmann : _ quantum chaos : an introduction _ ( cambridge university press , cambridge , uk , 1999 ) ; + f. haake : _ quantum signatures of chaos _ ( springer , 2nd edition , 2001 ) .
Abstract:

Answer: we develop a uniform semiclassical trace formula for the density of states of a three - dimensional isotropic harmonic oscillator ( ho ) , perturbed by a term @xmath0 . this term breaks the u(3 ) symmetry of the ho , resulting in a spherical system with so(3 ) symmetry . we first treat the anharmonic term for small @xmath1 in semiclassical perturbation theory by integration of the action of the perturbed periodic ho orbit families over the manifold @xmath2p@xmath3 which is covered by the parameters describing their four - fold degeneracy . then we obtain an analytical uniform trace formula for arbitrary @xmath1 which in the limit of strong perturbations ( or high energy ) asymptotically goes over into the correct trace formula of the full anharmonic system with so(3 ) symmetry , and in the limit @xmath1 ( or energy ) @xmath4 restores the ho trace formula with u(3 ) symmetry . we demonstrate that the gross - shell structure of this anharmonically perturbed system is dominated by the two - fold degenerate diameter and circular orbits , and _ not _ by the orbits with the largest classical degeneracy , which are the three - fold degenerate tori with rational ratios @xmath5 of radial and angular frequencies . the same holds also for the limit of a purely quartic spherical potential @xmath6 .