closed - orbit theory was first introduced by du and delos  @xcite and bogomolny  @xcite some twenty years ago to interpret the modulations observed in the photo - absorption spectra of hydrogenic rydberg atoms in a magnetic field close to the ionization threshold .
since that time , it turned out to be a powerful and flexible tool for the semiclassical interpretation of a variety of spectra .
it has been used to describe atoms in electric  @xcite as well as parallel  @xcite or crossed  @xcite electric and magnetic fields . in the case of non - hydrogenic atoms ,
the influence of the ionic core can be modelled either by means of an effective classical potential  @xcite or in terms of quantum defects  @xcite . recently
, closed - orbit theory has even been shown to be applicable to the spectra of simple molecules in external fields  @xcite .
a complete description of photo - absorption spectra requires the calculation of the energies @xmath2 of the excited atomic states and the strengths of the spectral lines , which is characterized by the dipole matrix elements @xmath3 between the initial state @xmath4 and the rydberg state @xmath5 , where @xmath6 is the component of the dipole operator describing the polarization of the exciting laser field .
these quantities are neatly summarized in the response function @xmath7 where @xmath8 denotes the retarded green s function .
closed - orbit theory provides a semiclassical approximation to the quantum response function  ( [ gdef ] ) , which splits into a smooth part and an oscillatory part of the form @xmath9 where the sum extends over all classical closed orbits starting from the nucleus and returning to it after having been deflected by the external fields , @xmath10 is the classical action of the closed orbit , and the amplitude @xmath11 describes its stability and its starting and returning directions .
its precise form depends on the geometry of the external fields . in section  [ sec : smatrix ] , it will be specified for systems with and without a rotational symmetry .
although the closed - orbit sum  ( [ resoscgen ] ) appears to provide a straightforward means of calculating the response function from the classical closed orbits , this is actually not the case because the sum usually diverges for real energies @xmath12 .
thus , the quantal information can not be extracted directly from the semiclassical expansion .
one particular and widely applicable method to overcome the convergence problems of the closed - orbit sum is semiclassical quantization by harmonic inversion @xcite . for the hydrogen atom in a magnetic field , this method has been shown @xcite to be capable of extracting semiclassical eigenenergies and transition matrix elements from a closed - orbit sum .    in the present paper we will investigate how these results can be generalized to the hydrogen atom in crossed electric and magnetic fields .
this problem is considerably harder than the treatment of the diamagnetic hydrogen atom , which possesses a rotational symmetry around the field axis . due to that symmetry , in classical mechanics the angular momentum around the field axis
is conserved .
so is , in quantum mechanics , the magnetic quantum number @xmath13 . in crossed fields ,
the rotational symmetry is broken . as a consequence
, the selection rules for the @xmath13-quantum number no longer hold , and a multitude of additional lines appears in the quantum spectrum . at the same time , the determination of classical closed orbits gets significantly more difficult because three non - separable degrees of freedom have to be dealt with .
a detailed description of the intricate pattern of closed orbits and their bifurcations was given in an accompanying paper @xcite .
that data forms the basis of the present work , where the semiclassical treatment of the crossed - fields hydrogen atom will be dealt with , and we will freely use the nomenclature introduced in @xcite .    after the essential properties of the crossed - fields hamiltonian have been summarized in section  [ sec : classham ] , we start , in section  [ sec : smatrix ] , with a derivation of the closed - orbit formula  ( [ resoscgen ] ) in the context of the @xmath0-matrix formulation of closed - orbit theory introduced recently by granger and greene @xcite .
we show that the novel framework can be extended to the crossed - fields situation , and we clarify some misleading conclusions arrived at in @xcite .
section  [ sec : xqm ] describes the quantum spectrum under study , and section  [ sec : scllo ] compares it to a semiclassical spectrum in low resolution . in section  [ sec : sclhi ] , the results of a high - resolution semiclassical quantization using the technique of harmonic inversion are presented .
the semiclassical spectrum correctly identifies the strongest spectral lines , but it fails to describe finer details of the quantum spectrum . in section  [ sec : sclrec ] , we compare a quantum recurrence spectrum to the classical data to show that the principal source of this difficulty lies in the abundance of closed - orbit bifurcations .
uniform approximations provide a tool to cope with the divergences introduced into semiclassical spectra by bifurcations of classical orbits .
a general technique for their construction is described in section  [ sec : uniform ] , and uniform approximations for the two types of generic codimension - one bifurcations identified in @xcite are derived .
finally , in section  [ ssec : unifrecur ] we demonstrate how uniform approximations can be incorporated into recurrence spectra , thus paving the way for their inclusion into the high - resolution semiclassical quantization by harmonic inversion .
throughout this work , we will assume the magnetic field to be directed along the @xmath14-axis and the electric field to be directed along the @xmath15-axis . in atomic units ,
the hamiltonian describing the motion of the atomic electron then reads @xmath16 where @xmath17 and @xmath18 denote the magnetic and electric field strengths , respectively , @xmath19 , @xmath20 , and @xmath21 is the @xmath14-component of the angular momentum vector . by virtue of the scaling properties of the hamiltonian  ( [ specham ] ) , if all classical quantities are multiplied by suitable powers of the scaling parameter @xmath22 the dynamics can be shown not to depend on the energy @xmath12 and the field strengths @xmath17 and @xmath18 separately , but only on the scaled energy @xmath23 and the scaled electric field strength @xmath24 . in particular , classical actions scale according to @xmath25 .
thus , the semiclassical limit of large classical actions corresponds to the limit of large @xmath26 .
the way of recording a quantum spectrum which is best suited for semiclassical investigations is scaled - energy spectroscopy .
a spectrum then consists of a list of the scaling parameters @xmath27 characterizing the quantum states for given scaled energy @xmath28 and scaled electric field strength @xmath29 .
scaled - energy spectroscopy offers the advantage that the underlying classical dynamics does not change across the spectrum .
it will be adopted throughout this work .
the basic observation fundamental to all of closed - orbit theory is a partition of space into physically distinct regions . in the core region
close to the nucleus , the rydberg electron interacts in a complicated manner with all electrons of the ionic core .
this interaction is manifestly quantum mechanical in nature , it can not be described in the framework of semiclassical theories . on the other hand ,
the interaction of the rydberg electron with the external fields is much weaker in the core region than its interaction with the core , so that the fields can safely be neglected .
therefore , a description of the core obtained in the field - free case can be used . in particular
, the initial state of the photo - absorption process is assumed to be localized in the core region and not to be influenced by the external fields .    in the long - range region far away from the nucleus , on the other hand , the external fields play a dominant role , whereas there is no interaction with the ionic core except for the coulomb attraction of its residual charge . in this region , the dynamics of the rydberg electron is well - suited for a semiclassical description .
it is independent of the details of the ionic core .    in order to establish a link between the dynamics in the core and long - range regions ,
a matching region is assumed to exist at intermediate distances from the nucleus where both the external fields and the interaction with the core are negligible .
thus , in the matching region the simple physics of an electron subject to the residual coulomb field of the core is observed .
recently , granger and greene  @xcite proposed a novel formulation of the theory based on ideas borrowed from quantum - defect theory .
their formulation achieves a clear separation between properties of the external field configuration and the ionic core , which are encoded in separate @xmath0-matrices .
suitable approximations to the core and the long - range @xmath0-matrices can be derived independently .
therefore , the formalism can be expected to allow a generalization of closed - orbit theory to atoms with ionic cores exhibiting more complicated internal dynamics than have been treated so far .
the derivation given by granger and greene treated the case of an atom in a magnetic field only .
it will now be extended in such a way that it holds for combined electric and magnetic fields with arbitrary field configurations . to this end ,
the ansatz and basic formulae of granger and greene s theory will be summarized in this section .
a more detailed treatment can be found in their paper  @xcite . in subsequent sections
, we will then turn to a discussion of the long - range scattering matrices pertinent to different external field configurations .    to lay the foundation for a definition of the @xmath0-matrices ,
we pick a basis set @xmath30 and @xmath31 of wave functions of the rydberg electron valid in the core and long - range regions , respectively , and expand in terms of spherical harmonics @xmath32 the channel index @xmath33 is to be understood as a double index @xmath34 characterizing the spherical harmonics . when studying a complicated atom with more than one relevant state of the core
, additional information can be included in the channel functions @xmath35 .    in the matching region , the radial function matrices @xmath36 and @xmath37 can both be expressed in terms of radial coulomb functions .
we use the functions @xmath38 and @xmath39 satisfying outgoing and incoming wave boundary conditions , respectively , given by  @xcite and choose the radial functions to be of the form and @xmath40 used by granger and greene  @xcite differ from those used by robicheaux  @xcite in that they are energy - normalized in rydberg rather than in hartree units .
the radial function matrices given here agree in normalization with those adopted by robicheaux , whereas the matrices used by granger and greene are inconsistent with their equation ( 12 ) . ]
@xmath41\;,\ ] ] @xmath42\;.\ ] ] physically , these choices mean that the basis function @xmath30 is a superposition of a single incoming wave in channel @xmath33 and the outgoing waves in different channels produced from it by scattering off the core .
similarly , @xmath31 consists of an outgoing wave in channel @xmath33 and the returning waves generated by scattering off the external fields . the scattering matrices @xmath43 and @xmath44 thus summarize the physical properties of the core and the external fields , respectively .
they are determined by the condition that the radial functions obey suitable boundary conditions , i.e. @xmath36 is regular at the origin , whereas @xmath37 vanishes or satisfies outgoing - wave boundary conditions at infinity for bound and free states , respectively . for hydrogen , @xmath43 is the identity matrix .    following previous work by robicheaux  @xcite , granger and greene
derive the following expression for the response function  ( [ gdef ] ) : @xmath45 where the vector @xmath46 comprises the energy - dependent dipole matrix elements @xmath47 between the initial state and the core - region channel wave functions .
for hydrogen they can be computed explicitly ( see , e.g. , @xcite or @xcite ) .
the terms of the series  ( [ resseries ] ) embody contributions from paths where the rydberg electron takes zero , one , two , etc
.  trips out into the long - range region and back to the core before interfering with the initial outgoing wave . in the semiclassical approximation ,
@xmath44 will be given in terms of closed orbits .
a returning wave is associated with each returning classical orbit . by a general ionic core
, it is scattered into all directions .
the parts of the wave scattered into the outgoing direction of a closed orbit will then follow this orbit until they return to the core again .
thus , core scattering leads to a concatenation of different closed orbits  @xcite . in hydrogen
, the coulomb center scatters the incoming wave back into its direction of incidence , so that there is no coupling of closed orbits .
terms describing repeated scattering off the external fields are therefore absent from the sum , and the hydrogen response function can be decomposed into a smooth part @xmath48 which is the same as in the field - free case and contains `` direct '' contributions where the electron does not scatter off the external fields at all , and an oscillatory part @xmath49 generated by the electron going out into the long - range region and being scattered back to the nucleus .
it is this part which describes the impact of the external fields .
the basis for a semiclassical approximation is provided by the retarded green s function @xmath50 describing the propagation of the electron from @xmath51 to @xmath52 at the energy @xmath12 .
it can be expanded in terms of the channel functions as @xmath53 with @xmath54 the long - range scattering matrix is related to the green s function matrix by @xcite @xmath55^{-1 }       \underline{g}(r_0,r_0 ) [ \underline{f}^-(r_0)]^{-1 } \;,\ ] ] where @xmath56 is the matching radius , @xmath57 is the diagonal matrix @xmath58 comprising the radial wave functions , and @xmath59 denotes the part of @xmath60 satisfying incoming - wave boundary conditions at the final radius @xmath61 .
the latter condition ensures that only electron paths approaching the matching radius from the long - range region contribute to @xmath44 , whereas paths that traverse the core region are omitted .      to obtain a semiclassical approximation to the long - range scattering matrix
, we make use of the semiclassical green s function derived by gutzwiller  @xcite @xmath62 where the sum extends over all classical trajectories leading from @xmath51 to @xmath52 at the energy @xmath12 , @xmath1 is the number of degrees of freedom , @xmath0 is the classical action along the trajectory , @xmath63 the number of caustics along the trajectory , and @xmath64          \frac{\partial^2 s}{\partial e \partial{\boldsymbol}x ' } &          \frac{\partial^2 s}{\partial e^2 }        \end{array }      \right)\ ] ] is the amplitude for the contribution of the trajectory . by ( [ greenchannel ] ) , we obtain a semiclassical approximation to the green s function matrix @xmath65 as usual in semiclassics , the integrals will be evaluated in the stationary - phase approximation . it yields a sum over all classical trajectories leaving the matching sphere at a direction given by @xmath66 and returning to it at @xmath67 .
the condition that @xmath68 obeys incoming - wave boundary conditions at the final radius translates into the condition that only orbits going out from the matching sphere into the long - range region and then returning to @xmath56 are to be included , whereas orbits passing through the core region are omitted .
if all factors in the integrand except for the exponential are assumed to vary slowly , the stationary - phase approximation reads @xmath69 where @xmath70 is the number of negative eigenvalues of the hessian matrix of @xmath0 occurring in the prefactor .
because the initial state is assumed to be well localized , it is clear that the outgoing waves generated by the photo - excitation originate in the immediate neighborhood of the nucleus .
therefore , only trajectories leaving the matching sphere radially need to be included in ( [ sclgreensp ] ) . by the same token
, the trajectories can be assumed to return to the matching radius radially .
thus , they are parts of closed orbits starting precisely at the nucleus and returning there .    by transforming ( [ greenddef ] )
to spherical coordinates and making use of the relations @xmath71 the amplitude factor @xmath6 for radial trajectories can be simplified to @xmath72 the determinants occurring in ( [ sclgreensp ] ) combine to @xmath73    = & \det\frac{\partial(p_\vartheta',p_\varphi',p_\vartheta , p_\varphi ) }               { \partial(\vartheta,\varphi , p_\vartheta , p_\varphi ) } \cdot      \left(\det\frac{\partial(-p_\vartheta',-p_\varphi',p_\vartheta , p_\varphi ) }                     { \partial(\vartheta',\varphi',\vartheta,\varphi ) }      \right)^{-1 } \\[.5ex ]    = & \det\frac{\partial(\vartheta',\varphi ' ) }               { \partial(p_\vartheta , p_\varphi ) } \;.
\end{split}\ ] ] with these results , the green s function matrix assumes the form @xmath74    the determinant in the denominator of ( [ sclgreenchannel2 ] ) measures the dependence of the final angular momenta of the trajectory upon the starting angles .
as it stands , it suffers from the singularities of the spherical coordinate chart : at the poles , neither the angle @xmath75 nor the angular momenta @xmath76 and @xmath77 are well defined , so that close to the poles , the calculation of the determinant becomes numerically unstable .
the determinant can be rewritten in the form @xcite @xmath78 with a @xmath79-determinant @xmath80 devoid of any singularities .
the parameter @xmath80 was already used in @xcite to study the bifurcations of closed orbits .
we showed there that a closed orbit bifurcates if and only if @xmath81 . with the form  ( [ mprime ] ) of the stability determinant , the semiclassical green s function matrix reads @xmath82 which is free of any singularities introduced by the spherical coordinates .    by virtue of ( [ sgreen ] )
, the semiclassical long - range scattering matrix reads @xmath83 this expression can be further simplified if , for excited states close to the ionization threshold , the radial wave functions @xmath84 are approximated by the zero - energy wave functions , and the hankel functions are replaced with their asymptotic forms for large arguments  @xcite @xmath85 this approximation has proven accurate in many cases of interest , but it was called into question by granger and greene @xcite . it will be discussed further in section [ sec : corot ] , where we will show that there is no reason to doubt its reliability .
it leads to @xmath86 because , due to the conservation of energy , @xmath87 if @xmath88 . in equation  [ scrossed0 ] , the channel indices @xmath89 are finally written out explicitly .    for a radial trajectory in a hydrogen atom
going out from the nucleus to @xmath90 at zero energy , the action is @xmath91 , so that @xmath92 is the action of a closed orbit , measured from its start at the nucleus to its return .
the semiclassical long - range @xmath0-matrix finally reads @xmath93 both the action @xmath10 and the stability determinant @xmath80 are evaluated at the nucleus rather than on the matching sphere .
the response function is given by @xmath94 where the maslov index @xmath95 was increased by 1 to absorb an additional phase , and the function @xmath96 with the core - region matrix elements @xmath97 given by  ( [ channelmat ] ) , characterizes the initial state and the exciting photon . through the @xmath97 ,
the function @xmath98 is energy - dependent . in accordance with the choice of zero - energy radial wave functions in the @xmath0-matrix elements ,
@xmath98 will be evaluated at zero energy .
this approximation has proven accurate in all applications of closed - orbit theory considered in the literature so far . however , from the @xmath0-matrix theory derivation it is obvious that the energy - dependence of both the dipole matrix elements @xmath97 and the @xmath0-matrix elements can easily be included should the need arise . the semiclassical response function ( [ gcrossed ] ) has the anticipated form ( [ resoscgen ] ) with @xmath99      an atom in a single ( electric or magnetic ) external field possesses a rotational symmetry around the field axis , which must be taken into account in the derivation of the closed - orbit formulae .
the symmetry gives rise to a conserved magnetic quantum number @xmath13 , so that the angular momentum quantum number @xmath100 remains the only relevant channel index .
the semiclassical scattering matrix reads @xcite @xmath101 where @xmath102 @xmath103 is the number of poles of @xmath104 encountered along the trajectory , and the sum includes all classical trajectories with azimuthal angular momentum @xmath13 joining the circles given by polar angles @xmath105 and @xmath106 on the matching sphere . if the radius of the matching sphere is much larger than the extent of the initial state , the trajectories can again be assumed to leave the sphere and return to it radially . strictly speaking
, this condition can only be met if @xmath107 , which we will assume in what follows .
if @xmath108 , the initial angular velocity @xmath109 must be non - zero , but it will be small if the matching radius is large . in this case ,
the trajectory will not actually close at the nucleus , but swing by at a short distance .    using , as above , the radial wave functions at zero energy ,
we obtain the semiclassical scattering matrix @xmath110 and the response function @xmath111 with @xmath112 and @xmath113 this result has the form ( [ resoscgen ] ) with @xmath114 it differs from the result obtained previously by du and delos @xcite in that in their work the amplitude factor @xmath104 of ( [ agg ] ) is replaced with @xmath115 this discrepancy was noted and numerically investigated by granger and greene @xcite .
they attribute it to the approximation of using zero - energy wave functions , which can easily be avoided in the @xmath0-matrix theory , but is an integral part of the derivation given by du and delos .    ) , solid line ) and after du and delos @xcite ( equation  ( [ add ] ) , dashed line ) for the closed orbit perpendicular to the magnetic field as a function of the scaled energy .
the scaled matching radius is @xmath116 . ]    for the closed orbit perpendicular to the field in the diamagnetic kepler problem and a scaled matching radius of @xmath116 , the amplitudes ( [ agg ] ) and ( [ add ] ) are plotted in figure [ ampefig ] .
this figure is similar to figure 1 in @xcite , although for the latter the matching radius is not given .
the agreement is excellent at scaled energies close to zero , but becomes poor if the energy decreases .
however , contrary to their conclusions , the lack of agreement is not due to the zero - energy approximation , but rather to the dependence of the amplitudes on the matching radius .
this statement can be verified most conveniently if the motion is described in semiparabolic coordinates @xmath117 if the trajectory is recorded as a function of a parameter @xmath118 related to the time @xmath119 by @xmath120 and a prime denotes differentiation with respect to @xmath118 , for trajectories with vanishing azimuthal angular momentum the equations of motion in the coulomb region read @xmath121 these equations are devoid of any singularities , so that they can conveniently be used to discuss the motion close to the nucleus . the transformation inverse to ( [ sparcoord ] )
is given by @xmath122 the momenta transform according to @xmath123 note that the transformation from semiparabolic to cartesian coordinates is not one - to - one , but that @xmath124 and @xmath125 are fixed up to the choice of sign only .    to evaluate ( [ agg ] ) and ( [ add ] ) , the derivatives @xmath126 and @xmath127 , must be calculated and their dependence on the matching radius @xmath61 must be determined .
as the radial trajectory specified by a starting angle @xmath105 is independent of the radius where the angle is measured , the @xmath61-dependence of the derivatives is determined by the returning trajectories only .
it can be evaluated as follows :    we arbitrarily fix the returning time of a closed orbit at @xmath128 , so that @xmath129 . the solution to ( [ speqns ] ) describing a trajectory returning at an angle @xmath106
is given by @xmath130 where the coefficients were chosen to satisfy the conservation of energy and to give the correct returning angle after a transformation to cartesian coordinates .
the second expression in each line follows from @xmath131 , whence for @xmath132 @xmath133    equations of motion for the derivatives @xmath134 and @xmath135 are obtained by linearizing ( [ speqns ] ) . since ( [ speqns ] ) is already linear , the derivatives satisfy the same equations of motion as the coordinates themselves as long as the electron moves in the coulomb region .
there the solutions read @xmath136 and @xmath137 equation ( [ scr ] ) yields @xmath138 so that the coefficients @xmath139 can be identified with the values of the derivatives obtained at @xmath140 .
analogous expressions hold for @xmath135 .    from ( [ spmomenta ] ) ,
the amplitude ( [ agg ] ) @xmath141 can be evaluated .
it is independent of @xmath61 , as could have been anticipated from the fact that @xmath76 is a component of the total angular momentum and thus is conserved along the trajectory once the electron has entered the coulomb region .
the amplitude @xmath142 can also , up to an immaterial choice of sign , be identified with the monodromy matrix element @xmath143 introduced by bogomolny  @xcite to describe the semiclassical amplitudes , so that the amplitudes derived by granger and greene from the @xmath0-matrix theory agree with bogomolny s .
similarly , the amplitude ( [ add ] ) used by du and delos reads , by ( [ spinverse ] ) , @xmath144 \\      & = \frac{1}{a } + { \cal o}\left(\sqrt{r}\right ) \;.    \end{split}\ ] ] thus , the amplitudes @xmath104 and @xmath145 agree in the limit of vanishing matching radius , but the amplitude @xmath145 proposed by du and delos exhibits a strong dependence on @xmath61 , whereas the amplitude @xmath104 given by granger and greene does not .
these findings can also be confirmed numerically .
figure [ amprfig ] shows the two amplitudes for the closed orbit perpendicular to the magnetic field at a scaled energy of @xmath146 as a function of the scaled matching radius @xmath147 .
the dependence of @xmath145 on @xmath147 is considerable .    ) , solid line ) and after du and delos @xcite ( equatiion  ( [ add ] ) , dashed line ) for the closed orbit perpendicular to the magnetic field as a function of the matching radius at @xmath146 . ]
we have thus shown that , contrary to the conclusion reached by granger and greene , the discrepancy between their semiclassical amplitude and that obtained by du and delos is not due to the zero - energy approximation , but rather due to the choice of a finite matching radius .
in addition , the amplitude derived by granger and greene is not specific to the @xmath0-matrix formulation , it agrees with the result derived earlier by bogomolny in the context of a semiclassical wave function formalism .
nevertheless , as it eliminates the need to specify a finite matching radius and allows one to calculate all classical quantities at the nucleus , it seems more appropriate than the amplitude given by du and delos , which introduces a certain arbitrariness in the choice of a matching radius .
if schrdinger s equation for the crossed - fields hydrogen atom is rewritten in terms of the scaled energy and the scaled electric field strength , a quadratic eigenvalue problem for the scaling parameter @xmath26 is obtained .
an exact numerical method of solution for the quadratic eigenvalue problem has become available only recently @xcite .
we resort to the method introduced by main @xcite , which relies on an approximate linearization of the eigenvalue problem to compute eigenvalues in a small spectral interval .
the accuracy of the linearization can be verified by comparing results calculated using different overlapping intervals .
the eigenvalues are obtained to a relative accuracy of at least @xmath148 , which is far beyond the typical accuracy of semiclassical approximations , so that the algorithm is well suited to this study .     and the scaled electric field strength @xmath149 .
the initial state is @xmath150 , the light is polarized along the magnetic field axis .
the plot shows the squared dipole matrix elements , which for graphical reasons are multiplied by @xmath26 .
the strengths of the extraordinarily strong lines of the lowest @xmath1-manifolds at @xmath151 are scaled down by a factor of 0.2 . ]    in the following we will discuss quantum and semiclassical photo - absorption spectra obtained for the scaled energy @xmath152 and the scaled electric field strength @xmath149 with the initial state @xmath150 and light linearly polarized along the magnetic field axis . a quantum spectrum for these parameter values is shown in figure  [ quantumfig ] . as for a semiclassical analysis ( see section  [ sec : sclrec ] ) it is essential to have as many eigenvalues available as possible ,
the calculation was extended up to @xmath153 .
the spectrum shown in figure  [ quantumfig ] contains nearly 30,000 lines , many of which are too weak to be discernible in the plot .
the eigenenergies of the field - free hydrogen atom satisfy @xmath154 so that in the scaled spectrum the unperturbed @xmath1-manifolds appear equidistantly spaced at @xmath155 these spacings can clearly be recognized in figure  [ quantumfig ] . at low values of @xmath26 ,
neighboring @xmath1-manifolds are isolated .
furthermore , in this region the magnetic quantum number @xmath13 is nearly conserved .
this is apparent from the fact that each @xmath1-manifold contains a central group of strong levels corresponding to @xmath107 , which can be excited even at @xmath156 , and adjacent groups of considerably weaker levels with @xmath157 .
levels with higher magnetic quantum numbers are too weak in this region to be seen in the figure . at higher values of @xmath26 , the conservation of @xmath13 is violated , and individual @xmath1-manifolds acquire strong side bands . at even higher @xmath26 , different @xmath1-manifolds strongly overlap . throughout the spectral range shown , groups of strong lines indicating the centers of different @xmath1-manifolds are clearly discernible .
a semiclassical approximation to a scaled photo - absorption spectrum is obtained if the closed - orbit theory formulae of section  [ sec : cocrossed ] are rewritten in terms of scaled quantities , viz .
@xmath158 with @xmath159    when low - resolution photo - absorption spectra are to be calculated from ( [ goscscal ] ) , a method of cut - off must be adopted to deal with the divergence of the semiclassical closed - orbit sum . for this section ,
we choose a gaussian cut - off , i.e. ( [ goscscal ] ) is replaced with @xmath160 so that orbits with scaled actions larger than the cut - off action @xmath63 are smoothly suppressed .
this smoothing corresponds to a convolution of the quantum signal with a gaussian of width @xmath161 .     and ( b ) @xmath162 . ]    to facilitate the comparison of ( [ goscsmooth ] ) with the convoluted quantum spectrum
, we added the smooth part of the spectrum to @xmath163 , which was calculated by convoluting the quantum spectrum with a gaussian of width @xmath164 .
this function is broad enough to wipe out the distinction between neighboring principal quantum numbers .
results obtained for @xmath165 and @xmath162 are shown in figure  [ xlofig ] . in both cases
it is apparent that the large - scale structure of equidistant principal quantum numbers is well reproduced by the semiclassical approximation . in the quantum spectra
, the substructure of the individual @xmath1-shells can be discerned to a certain degree , given by the smoothing width @xmath161 . in the case of @xmath165 ,
much of this fine structure is also present in the semiclassical spectrum , but often the agreement is not good quantitatively . in particular , the peaks corresponding to the lowest @xmath1-manifolds are considerably wider in the semiclassical than in the quantum spectrum .
if the cut - off action is increased to @xmath162 , finer details are resolved in the quantum spectrum . at the same time , the semiclassical closed - orbit sum becomes more oscillatory to reproduce this fine structure .
it appears , however , to be somewhat over - oscillatory , developing structures absent from the quantum spectrum .
this type of behavior is typical of closed - orbit sums in non - integrable systems .
thus , it can be questioned if the low - resolution closed - orbit sum can meaningfully be extended to even longer orbits . a high - resolution quantization based on the present semiclassical approximation
will be presented in the following section .
for the calculation of a scaled semiclassical spectrum , the method of semiclassical quantization by harmonic inversion of @xmath166 function signals @xcite can be applied .
this technique requires the inclusion of closed orbits up to a maximum scaled action , i.e. it replaces the gaussian cut - off used for the low - resolution semiclassical spectra presented in the previous section with a rectangular cut - off . a rough estimate for the required cut - off action
can be obtained by means of perturbation theory @xcite .
@xmath167 for the case @xmath168 and @xmath169 , i.e. @xmath170 , this estimate yields @xmath171 .    according to ( [ smaxpert ] ) , to compute levels at high quantum numbers @xmath1 a long semiclassical signal
is needed , which can be hard or even impossible to obtain .
we calculated closed orbits up to @xmath172 , so that the orbital data is available for nearly 18,000 closed - orbit multiplets .
however , for reasons to be described in section  [ sec : sclrec ] a useful semiclassical signal can be constructed up to @xmath173 only , so that , from the above estimate , the semiclassical calculation can not reach manifolds much higher than @xmath174 . on the other hand
, the semiclassical approximation must be expected to yield more accurate results for higher quantum numbers .
thus , when a high - resolution semiclassical spectrum is to be calculated , a compromise must be made between the contradictory requirements of describing a spectral region at sufficiently high quantum numbers and with a sufficiently low spectral density .
@r|rrr@ @xmath1 & @xmath175 ( scl . ) & @xmath175 ( qm . ) & @xmath176 + & & 9.88321 & 1.3617 + & & 9.91431 & 3.1145 + & & 9.97747 & 1.7474 + & 10.05366 & 10.05912 & 51.0512 + 6 & 10.09551 & 10.09621 & 20.9313 + & 10.15461 & 10.15378 & 7.0060 + & & 10.24076 & 0.9608 + & & 10.26612 & 2.0777 + & & 10.31803 & 1.9385 +    & & 11.56497 & 2.5663 + & 11.60898 & 11.60820 & 2.5875 + & 11.66889 & 11.67341 & 2.3104 + & 11.72048 & 11.73128 & 32.8808 + 7 & & 11.75121 & 16.7278 + & & 11.78850 & 10.0092 + & & 11.84856 & 5.6249 + & & 11.92188 & 1.9229 + & & 11.95821 & 1.7923 + & & 12.01338 & 2.4821 +    & & 13.23441 & 1.3668 + & & 13.25629 & 2.5141 + & & 13.30255 & 1.9971 + & 13.36921 & 13.36913 & 2.8189 + & 13.40177 & 13.40568 & 30.8875 + 8 & 13.44313 & 13.43744 & 16.0829 + & 13.48737 & 13.48146 & 4.8263 + & & 13.54340 & 4.3111 + & & 13.59258 & 1.0747 + & & 13.61133 & 1.9475 + & & 13.65111 & 1.4081 + & & 13.70866 & 2.9676 +    & & 14.91192 & 2.1880 + & & 14.94654 & 2.9922 + & & 14.99711 & 1.4563 + 9 & [ 2ex]15.06960 & 15.06470 & 3.2226 + & & 15.07888 & 25.1866 + & & 15.10074 & 8.4317 +    [ cols=">,>,>,>",options="header " , ]     for the harmonic analysis of the closed - orbit sum we applied the method of @xmath166 function decimated signal diagonalization @xcite , which yields not only semiclassical eigenvalues and amplitudes , but also an error parameter estimating the precision of the eigenvalues .
results obtained for @xmath152 and @xmath149 with a signal length of @xmath177 are compiled in table  [ xscltab ] .
the table contains the quantum eigenvalues of @xmath26 and their dipole matrix elements for levels satisfying @xmath178 .
it is obvious at a glance that out of the multitude of spectral lines with intensities varying over many orders of magnitude ( most of which are not contained in the table ) only the strongest lines were detected in the semiclassical spectrum .
the semiclassical eigenvalues given are characterized by having small imaginary parts , small error parameters and large amplitudes as well as being stable with respect to a variation of numerical parameters .
the calculation operates at the edge of convergence , and in a few cases one can be in doubt whether a level should be included according to these fairly `` soft '' criteria , but in general a clear decision can be made .
semiclassical values for the transition strengths are not given because they are not reasonably well converged and depend strongly on the numerical parameters .
one might expect that in each @xmath1-manifold it is the strongest lines that are detected semiclassically , and in general this expectation is confirmed by the numerical data .
this can clearly be seen , e.g. , in the manifold @xmath179 , which contains the most stably converged lines in the spectrum .
there are , however , a few conspicuous exceptions , e.g. at @xmath180 , where strong lines are missing whereas comparatively weak lines are found . for @xmath181 ,
no lines at all can be computed from the given semiclassical signal .
if the signal length is decreased to @xmath182 , the three strongest lines appear in the spectrum in this manifold .    at higher
@xmath1 , the number of strong lines in the quantum spectrum increases .
so does the number of lines in the semiclassical spectrum until @xmath183 , where only three semiclassical lines are found .
they appear rather arbitrarily scattered across the quantum spectrum , and their convergence is notably worse than in lower manifolds .
it is clear that in this @xmath1-shell the semiclassical quantization with the given signal is about to break down . at @xmath184 ,
no lines can be detected semiclassically .
as , from the above discussion , this failure was to be expected because the required signal length becomes too large , the obvious way to improve convergence seems to be to use a longer signal .
however , if the signal length is increased to @xmath185 , no reasonably converged semiclassical lines can be found in any @xmath1-manifold .
neither are results improved if the technique of harmonic inversion of cross - correlated closed - orbit sums @xcite is applied .
this method has proven powerful in reducing the signal length required in a semiclassical quantization . in the present case , however , because the cross - correlation increases the total number of frequencies obtained from the harmonic inversion , the true eigenvalues are hidden among a multitude of spurious frequencies , and no useful results can be obtained .    for the time being , therefore , the results given in table  [ xscltab ] represent what can be achieved in the semiclassical quantization of the crossed - fields hydrogen atom .
they confirm the applicability of the closed - orbit theory approach in principle , but they also reveal a fundamental problem in its present formulation . from the analysis of the ideal test signal it is clear that the signal length available is sufficient for a stable signal analysis .
thus , if the semiclassical results are not good , the semiclassical signal itself , rather than the signal analysis , must be to blame .
this conclusion is confirmed by the observation that an increased signal length destroys the results rather than improves them .
we therefore start searching for a flaw in the construction of the semiclassical photo - absorption spectrum .
a conspicuous problem lies in the fact that the set of closed orbits available is incomplete . in no series of rotators or vibrators
can arbitrarily long orbits be calculated . in the case of vanishing electric field
there is a critical angle @xmath186 which the starting angles of both rotators and vibrators approach as the orbits get longer .
this convergence indicates that the orbits approach a separatrix between two families of tori in phase space .
if sufficiently long orbits are studied , there are many closed orbits with very similar initial conditions , so that the numerical search for closed orbits must eventually fail .    , @xmath149 .
]    the region of phase space where the unknown orbits are located is lying close to a separatrix , so that it is highly unstable .
the orbits can therefore be expected not to contribute much to the semiclassical signal .
the magnitude of an orbit s contribution to the closed - orbit sum ( [ gcrossed ] ) is determined mainly by its stability determinant @xmath80 .
figure  [ stabmfig ] shows the stability determinants of the vibrator orbits for @xmath152 , @xmath149 as a function of the scaled action .
different series of vibrators can clearly be discerned in the plot .
it is indeed unstable orbits with large @xmath187 that are missing in the data set , but on the other hand the stability determinants of the missing orbits are not large enough to regard the corresponding semiclassical amplitudes as negligibly small .
because a vast majority of orbits has small @xmath187 and was found , one can still hope that useful results can be obtained from the semiclassical signal , at least for quantum states not located in the separatrix region in phase space , but it is clear that the quality of the semiclassical signal is reduced by its incompleteness .    to assess in detail the detrimental effect of the missing orbits and of any other sources of error that may exist , we carry out a semiclassical analysis of the quantum spectrum .
according to equation  ( [ goscscal ] ) , in a scaled photo - absorption spectrum every closed orbit contributes a purely sinusoidal modulation to @xmath188 .
this contribution can be extracted from the spectrum either by a conventional fourier transform or by means of a high - resolution method . the spectral analysis yields information about classical orbits returning to the nucleus .
for this reason , the transformed spectrum is referred to as a recurrence spectrum .
high - resolution methods @xcite extract the scaled actions and scaled semiclassical amplitudes of individual orbits and thus yield more complete information about the semiclassical spectrum than the fourier transform , but they fail if the average density of closed orbits per unit of scaled action is too large . by contrast , due to its linearity the fourier transform can be applied to any part of the recurrence spectrum with equal ease and numerical stability , irrespective of the spectral density . in dense regions , it will not be able to identify individual closed orbits , but it will nevertheless yield a recurrence spectrum that can be compared to the classical data . in this section we will present results obtained by both the fourier transform and a high - resolution method .
the semiclassical recurrence spectra will be compared to classical results in order to identify the reason why the semiclassical signal is only partially suitable to a semiclassical quantization .    using either method , it is essential to note that the semiclassical closed - orbit formula can not be expected to yield accurate results for the lowest levels .
thus , the low @xmath1-manifolds must be excluded from the semiclassical analysis , i.e. the analysis is based on the quantum spectrum given in an interval @xmath189 $ ] instead of @xmath190 $ ] .
furthermore , to minimize the impact of boundary effects due to the finite length of the semiclassical spectrum , a smooth gaussian cut - off with width @xmath70 centered at @xmath191 is introduced .
the smoothing replaces the peaks of the semiclassical recurrence spectrum by gaussians of width @xmath192 .
the recurrence spectra presented here were calculated from the quantum spectrum shown in figure  [ quantumfig ] , for @xmath152 and @xmath149 , with @xmath193 , @xmath194 , and @xmath195 . for the high - resolution recurrence spectra ,
the method of @xmath166 function decimated signal diagonalization was used .     and
@xmath149 . sticks and squares : semiclassical closed - orbit amplitudes , stars : harmonic inversion of the quantum spectrum.,title="fig : " ]   and @xmath149 .
sticks and squares : semiclassical closed - orbit amplitudes , stars : harmonic inversion of the quantum spectrum.,title="fig : " ]    for low scaled actions , where only few closed orbits exist , the high - resolution analysis can be applied .
results are shown in figure  [ recurreimfig ] , which compares both the scaled actions and the real and imaginary parts of the semiclassical amplitudes extracted from the quantum spectrum to the classical results . for most closed orbits ,
the agreement is excellent .
exceptions occur for the shortest orbits , where the actions of rotator and vibrator orbits are too similar to be resolved by the harmonic inversion . at somewhat larger actions ,
the three orbits in each group fall apart into two rotator orbits with similar actions and a vibrator orbit with a slightly larger action .
these observations can be made even more clearly if the absolute values of the amplitudes are considered .
they are shown in figure  [ recurabsfig ] , where the results of the high - resolution analysis are also compared to those of the fourier transform .
notice that for the fourier transform the semiclassical amplitude is given by the area under a peak rather than the peak height , so that an immediate comparison to the high - resolution results is difficult . in figure
[ recurabsfig ] , the fourier transform was arbitrarily scaled so that the peak heights roughly match the values of the high - resolution amplitudes . for isolated orbits identified both in the fourier transform and the high - resolution spectrum
, the agreement between the two methods is excellent .
where several peaks overlap in the semiclassical spectrum , no direct comparison is possible because the peak phases can not be determined from the figure .
figure  [ recurabsfig ] also extends the results shown in figure  [ recurreimfig ] to higher actions . in this region
the density of closed orbits starts to increase because , on the one hand , rotators of the first series exist and , on the other , bifurcations of closed orbits generate additional orbits .
apart from the fact that many orbits can not be identified individually even by the high - resolution method , the most conspicuous feature of figure  [ recurabsfig ] is that for many orbits the semiclassical amplitudes calculated from the classical data are considerably larger than those extracted from the quantum spectrum . in some cases , this is obvious at a glance , but a closer inspection of the figure reveals that this phenomenon is rather common .
some specific cases will be described in detail in section  [ sec : uniform ] .
the occurrence of exceedingly large semiclassical amplitudes is a well - known problem of both closed - orbit and periodic - orbit theory .
it is associated with bifurcations of classical orbits , where , in the case of closed orbits , the stability determinant @xmath80 vanishes and the closed - orbit amplitude ( [ acrossed ] ) diverges .
close to the bifurcation , @xmath80 is small .
the semiclassical amplitude of the bifurcating orbit is therefore large and exceeds the value determined from the quantum spectrum . in a classical context , we have shown previously @xcite that vanishing @xmath80 is both a necessary and sufficient condition for a bifurcation of closed orbits . in the context of semiclassical closed - orbit theory ,
it is necessary to overcome the divergence of the closed - orbit formula occurring close to a bifurcation .
this problem will be addressed in section  [ sec : uniform ] , after the impact of the bifurcations on the semiclassical signal at hand has been investigated further .
whereas , in figure  [ recurabsfig ] , the vibrator orbits are sufficiently isolated to be resolved by both the harmonic inversion and the fourier transform across the entire range of actions , the rotators occur in groups of several orbits having nearly identical actions .
they are not resolved properly by either method . instead
, the fourier transform produces peaks describing the collective contribution of the orbits in a group .
the harmonic inversion fits this contribution with fewer actions and amplitudes than the actual number of orbits .
although the results can be expected to reproduce the quantum spectrum fairly well , the principal virtue of the high - resolution analysis  that it is capable of giving individual rather than collective contributions  is lost .
it is therefore pointless to extend the high - resolution analysis to higher actions unless a significantly longer quantum spectrum can be obtained , and only the fourier transform will be used henceforth .     of the recurrence spectrum with @xmath195 ( see text ) .
upper part : fourier transform of the quantum spectrum , lower part ( inverted ) : smoothed semiclassical recurrence spectrum . ]
figure  [ recurlongfig ] displays the fourier recurrence spectrum with smoothing @xmath195 for scaled actions up to @xmath196 and compares it to the semiclassical spectrum .
these results extend the semiclassical analysis of quantum spectra to significantly longer orbits than investigated in previous studies .
they allow a verification of closed - orbit theory all the way up to the long orbits .
it is immediately apparent from the figure that the quantum recurrence spectrum retains its pronounced peak structure .
this is to be expected from closed - orbit theory , and indeed the peak locations are given by the actions of closed orbits for long as well as for short orbits .
the basic idea of closed - orbit theory that recurrence peaks are related to classical closed orbits is therefore confirmed in principle even for very long orbits .    even for the largest actions considered ,
the quantum and semiclassical recurrence spectra agree quantitatively for some peaks . for most peaks , however , the peak heights in the quantum and semiclassical spectra disagree .
there are quantum peaks that are smaller in the semiclassical spectrum or even completely absent .
they can be attributed to missing orbits . on the other hand , in many cases
the semiclassical peaks are significantly higher than the quantum peaks , sometimes by several orders of magnitude .
exceedingly high peaks can be traced back to bifurcations of closed orbits if the possibility is ignored that a quantum peak can be small because orbits missing in the semiclassical spectrum interfere destructively with the orbits present .
this latter mechanism becomes the more implausible the larger the semiclassical peak is in comparison to the quantum peak .
taken together , the effects of missing orbits and of bifurcating orbits distort the semiclassical recurrence spectrum to the point where it can no longer be expected to provide a suitable basis for a quantization .
a close inspection of the recurrence spectrum suggests that the problem posed by bifurcating orbits is more severe .
exceedingly high peaks do not only occur frequently , but in addition the very fact that they are high increases their detrimental effect on the semiclassical photo - absorption spectrum .
unless a suitable scheme for dealing with bifurcating orbits can be devised , no improvement of the semiclassical signal can be expected .
we therefore turn to a description of the semiclassical treatment of bifurcating orbits by means of uniform approximations .
exceedingly large contributions of single orbits to a semiclassical spectrum arise when the orbits are too close to a bifurcation to be regarded as isolated , as is implicitly assumed by the stationary - phase approximation used in the derivation of the closed - orbit formula .
uniform approximations furnish a collective contribution of all orbits involved in a bifurcation .
this solution was first suggested by ozorio de almeida and hannay @xcite in the context of periodic - orbit theory .
their original approach was extended by different authors @xcite , so that today uniform approximations are a well - established tool of semiclassical physics . in reference
@xcite , we identified two types of generic closed - orbit bifurcations of codimension one .
the pertinent uniform semiclassical approximations will be derived in what follows .    in most cases of interest ,
a bifurcation destroys real orbits and turns them into complex ghost orbits that exist in the complexified classical phase space .
ghost orbits can yield palpable contributions to semiclassical spectra @xcite .
in particular , their knowledge is essential for the construction of uniform approximation . for the generic closed - orbit bifurcations ,
the ghost orbits were described along with the real orbits in reference  @xcite .    of particular importance
is the observation that bifurcations of codimension higher than one are relevant to semiclassics , although on a classical level they are not generically encountered .
they appear as sequences of generic bifurcations , which , if the individual bifurcations are sufficiently close , must be described collectively by a single uniform approximation .
several examples of uniform approximations for these complicated bifurcation scenarios have been described in the literature @xcite .
the principal requirement a uniform approximation must satisfy is to asymptotically reproduce the known isolated - orbits approximation when the distance from the bifurcation grows large , because in this limit the stationary - phase approximation can be expected to be accurate . in the following , we will advocate a somewhat heuristic technique for the construction of a uniform approximation , which is easy to handle and yields a smooth interpolation between the asymptotic isolated - orbits approximations on either side of the bifurcation .
it will first be described in general terms .
subsequently , uniform approximations describing the generic types of codimension - one bifurcations of closed orbits will be derived .
a bifurcation scenario is described by a normal form @xmath197 depending on @xmath198 variables  @xmath119 and @xmath199 parameters  @xmath200 such that for any fixed value of the parameters @xmath200 there are stationary points of @xmath197 corresponding to the closed orbits involved in the bifurcation .
the parameters @xmath200 must then depend on the energy @xmath12 to reproduce the bifurcations of the closed orbits .    for the uniform approximation we make the ansatz @xmath201 with @xmath202 here , the functions @xmath203 and @xmath204 as well as the parameter values @xmath205 have to be determined .
all of them must be smooth functions of @xmath12 .    to find the asymptotic behavior of the uniform approximation ( [ unifansatz ] )
far from the bifurcations , ( [ unifi ] ) is evaluated in the stationary - phase approximation , which yields @xmath206 where the sum extends over all stationary points @xmath207 of @xmath197 that are real at the given @xmath200 , @xmath208 is the hessian determinant of @xmath209 , and @xmath210 is the number of negative eigenvalues of @xmath211 .
this expression is supposed to reproduce the isolated - orbits approximation @xmath212 in this case , the sum extends over all closed orbits involved in the bifurcation that are real at the given energy @xmath12 .
if the normal form @xmath197 has been chosen suitably , there is a one - to - one correspondence between these orbits and the stationary points @xmath207 .
a comparison of  ( [ unifsp ] ) to  ( [ unifisol ] ) yields the conditions @xmath213 and @xmath214 these equations must be valid for real orbits . in most bifurcation scenarios ,
all orbits are real at least at certain energies . in these cases
, it appears natural to postulate ( [ unifamp ] ) also to hold for ghost orbits .
the parameter values one obtains are then smooth functions of the energy even at the bifurcations where the orbits become ghosts . in some instances , bifurcations involving only ghost orbits occur @xcite . in these cases , the condition  ( [ unifamp ] ) still produces smoothly varying parameters and enforces the desired asymptotics .
the numbers  @xmath210 of negative eigenvalues change discontinuously at a bifurcation . for
orbits which are real on either side of the bifurcation , so do the maslov indices contained in the semiclassical amplitudes  @xmath215 .
these changes must compensate each other if the values  @xmath216 are to be continuous across the bifurcation . for these orbits , therefore , the change of maslov index occurring in a bifurcation must be equal to the change in @xmath210 and can be determined from the normal form . for ghost orbits ,
maslov indices are not well defined classically .
they must be chosen such as to make @xmath216 continuous .
the normal form parameters  @xmath200 and the action @xmath203 can be determined from  ( [ unifaction ] ) .
they usually turn out to be unique .
the amplitude function  @xmath204 , on the contrary , is unknown .
once the parameters  @xmath200 have been found , ( [ unifamp ] ) specifies its values @xmath216 at the stationary points of @xmath197 .
these values , of course , do not suffice to identify @xmath204 uniquely , so that there is considerable freedom in the choice of @xmath204 .
usually , if there are @xmath33  orbits participating in the bifurcation scenario , we will approximate @xmath204 by a polynomial of degree  @xmath217 .
this choice is justified by the observation that the uniform approximation is needed only close to a bifurcation , where all orbits are close to @xmath218 .
thus , in the spirit of the stationary - phase approximation , the dominant contributions to the integral  ( [ unifi ] ) stem from the neighborhood of  @xmath218 , whereas the regions of large  @xmath119 do not contribute . a suitable approximation to  @xmath204 must therefore be precise close to the origin .
this is achieved by a taylor series expansion , which leads to the polynomial ansatz .
simple as it might appear , however , this choice can bring about a mathematical difficulty : a polynomial  @xmath204 diverges as @xmath219 , so that there is no guarantee that the integral  ( [ unifi ] ) will converge .
if it does not , its divergence is an artefact of the choice of @xmath204 , because by construction the regions of large @xmath119 should not significantly influence the value of the integral . in this case , a suitable regularization scheme must be applied .
it can be justified by verifying that the regularized integral possesses the correct asymptotics .
a slightly simpler form of the uniform approximation is obtained if the function @xmath204 is assumed to be a constant .
this approximation does not exactly reproduce the desired asymptotics , but as the transition across the bifurcation mainly results in a change of the stationary points of @xmath197 rather than essential changes in @xmath204 , it can be expected to capture the principal features .
it is clear from the above description that there is a certain arbitrariness in the procedure .
this arbitrariness can be reduced to the choice of a suitable amplitude function  @xmath204 , because by the splitting lemma and the classification theorems of catastrophe theory @xcite the uniform approximation can always be brought into the form  ( [ unifansatz ] ) by a suitable coordinate transformation , provided a normal form @xmath197 equivalent to the actual action function is given .    in the following sections , uniform approximations for the two generic codimension - one bifurcations described in @xcite
will be derived along the lines given here .
they turn out to be analogous to those for isochronous and period - doubling bifurcations of periodic orbits given by schomerus and sieber @xcite .
the simplest closed - orbit bifurcation is the creation of two orbits in a tangent bifurcation .
it is described by the fold catastrophe @xmath220 this normal form has stationary points at @xmath221 , which are real if @xmath222 .
its stationary values are  ( [ foldval ] ) @xmath223 by  ( [ unifaction ] ) , the actions @xmath224 and @xmath225 of the bifurcating orbits must satisfy @xmath226 for these equations to hold , one must assume @xmath227 if the orbits are real and @xmath228 , @xmath229 if they are ghosts .
these conditions determine how the orbits are to be associated with the stationary points of @xmath197 .
equation  ( [ foldaction ] ) can be solved for @xmath230 and @xmath231 the observation that the bifurcating orbits are real if @xmath222 and ghosts if @xmath232 fixes the sign of @xmath200 . both @xmath203 and @xmath200 have thus be determined .    for the semiclassical amplitudes , ( [ unifamp ] )
yields @xmath233 with the ansatz @xmath234 for the amplitude function @xmath204 , we can solve for the parameters @xmath235 and @xmath236 to obtain @xmath237 the uniform approximation thus takes the form @xmath238 with @xmath239 the integral @xmath240 can be evaluated in terms of the airy function @xcite as @xmath241 whereas @xmath242 is given by its derivative @xmath243 with these results , the uniform approximation ( [ uniffold ] ) can be computed once the classical quantities @xmath244 and @xmath245 are known .
after some rearrangements , ( [ uniffold ] ) can be found to agree with the uniform approximation derived by schomerus and sieber @xcite for isochronous bifurcations of periodic orbits , although its present form is much simpler .      the normal form for the symmetrized cusp catastrophe is given by @xmath246 it has stationary points at @xmath218 and @xmath221 and describes a pitchfork bifurcation , where two asymmetric orbits bifurcate off an orbit invariant under a reflection .
we denote their actions and amplitudes by @xmath247 and @xmath248 , respectively , where @xmath249 is understood to be the cumulative amplitude of both asymmetric orbits .    as @xmath250
, the reference action @xmath203 must be chosen equal to the action of the symmetric orbit .
the action difference is given by the stationary value of @xmath197 , which is @xmath251 , so that @xmath252 and @xmath253 the parameter @xmath200 has to be chosen positive if the asymmetric orbits are real , and negative otherwise .
here , @xmath254 was assumed to be positive .
if it is not , the normal form @xmath197 must be replaced with @xmath255 , which changes the sign of the stationary values .    due to the reflection symmetry
, the amplitude function must be an even function of @xmath119 .
we make the ansatz @xmath256 and solve  ( [ unifamp ] ) for the coefficients @xmath257    the complete uniform approximation reads @xmath258 with @xmath259 the integral @xmath240 can be evaluated analytically in terms of bessel functions @xcite : @xmath260 \;.    \end{split}\ ] ] although it is not apparent at first sight , @xmath240 is a smooth function of @xmath200 .
this can be verified if the series expansion @xcite @xmath261 with @xmath262 a power series in @xmath263 is used . in terms of @xmath262 , @xmath264 \ ; ,    \end{split}\ ] ] which is indeed smooth .
the second integral @xmath265 can be evaluated from    @xmath266 \\             & \phantom{{i}\pi\sqrt{|a|}{e}^{-{i}\frac{a^2}8}\bigg\{}+                    \frac a8 \,{e}^{{i}\pi/8 }                         \left[j_{-5/4}\left(\frac{a^2}8\right ) -                                j_{3/4}\left(\frac{a^2}8\right ) \right ] \\
& \phantom{{i}\pi\sqrt{|a|}{e}^{-{i}\frac{a^2}8}\bigg\ { }                    + \operatorname{sign } a \,\frac a8 \,{e}^{-{i}\pi/8 }                         \left[j_{-3/4}\left(\frac{a^2}8\right ) -                                j_{5/4}\left(\frac{a^2}8\right ) \right ]              \bigg\ } \;.    \end{split}\ ] ]    this derivation contains an interchange of differentiation and integration which achieves a regularization of the divergent integral @xmath265 .
it can be justified by verifying that the asymptotic behavior of ( [ i2cusp ] ) for @xmath267 agrees with the stationary phase approximation to ( [ ikcusp ] ) .
the formulae derived in the preceding sections give the uniform approximations directly in terms of the semiclassical actions and amplitudes .
this circumstance makes them easy to apply to scaled spectra : we simply put @xmath25 and @xmath268 .
as @xmath26 is varied , the bifurcation is not encountered because the classical mechanics does not change , so that the isolated - orbits approximation does not actually diverge . however ,
if @xmath26 is small , the action differences between the bifurcating orbits are also small , so that the presence of the bifurcation is felt and the isolated - orbits formula produces exceedingly large contributions . for large @xmath26 , the action differences also grow large , so that the isolated - orbits approximation should be recovered in the limit of large @xmath26 .    ) for the scaled spectrum at @xmath152 and @xmath269 . ]
these findings are illustrated in figure  [ unifscalfig ] for a pitchfork bifurcation taking place in the first series of rotators at a repetition number @xmath270 . at @xmath152 ,
the bifurcation takes place at the scaled electric field strength @xmath271 .
the data shown in figure  [ unifscalfig ] was calculated for @xmath152 and @xmath272 , which is sufficiently far away from the bifurcation for the asymptotic regieme to be reached within the range of @xmath26 shown .
as anticipated , in the limit of @xmath273 the complete uniform approximation agrees with the isolated - orbits formula .
the simple approximation also reproduces the beats correctly , but it has a smaller amplitude .    , @xmath274 and @xmath149 . ]
the scaled uniform approximation can be used to improve the semiclassical recurrence spectrum , but this requires some effort : whereas the isolated - orbits approximation yields @xmath166 function peaks in the recurrence spectrum , which are replaced with gaussians due to the smoothing of the recurrence spectrum ( see section  [ sec : sclrec ] ) , the uniform approximation is a complicated function of @xmath26 .
it must be subjected to a numerical fourier transform in the same was as the quantum spectrum if its contribution to the recurrence spectrum is to be evaluated .
because a bifurcation involves orbits with roughly equal actions , the uniform approximation will produce a recurrence peak at the appropriate action .
an example is shown in figure  [ unifrecfig ] .
it was calculated for the bifurcation already described in figure  [ unifscalfig ] .
the gaussian smoothing used in section  [ sec : sclrec ] was replaced with a rectangular window , so that a number of side peaks appear . in this case , the fourier transform of both the uniform approximation and the isolated - orbits approximation was taken over the rectangular window @xmath275 $ ] .
the bifurcating orbits have the scaled action @xmath276 , which is where the fourier peaks are centered in both approximations .
the peak produced by the uniform approximation is considerably smaller .
if this uniformization procedure is carried out for all excessively high bifurcation peaks , it should be possible to bring the semiclassical recurrence spectrum in figure  [ recurlongfig ] into agreement with its quantum counterpart . in practice , however , several obstacles stand in the way .
first of all , in many cases ghost orbits must be included in the uniform approximation .
they must be found and identified as pertinent to a given bifurcation before the uniformization can be performed .
furthermore , even if all relevant orbits are real , those orbits connected with each other in a bifurcation must be recognized in the data set .
this is by no means an easy task .
for example , if in a given series of rotators and for a given winding number a quartet of orbits appears , there are two different doublet orbits out of which they may have bifurcated , and it is not clear in general which of them must be taken for the uniform approximation . in a single case , this can be found out fairly comfortably by hand .
if many orbits are to be classified , however , it is essential to do the grouping automatically .
we have not yet been able to devise a practical algorithm for this task , so that an automatized uniformization of all bifurcation peaks is presently impossible .
apart from these rather technical difficulties , there are also some obstacles of more fundamental importance .
consider , e.g. , the two high semiclassical peaks at @xmath277 in figure  [ recurabsfig ] .
they are notably too high , and they are well - isolated from neighboring recurrence peaks , so that they may appear to be the ideal testing ground for the uniformization procedure .
these peaks are generated by vibrators with repetition numbers @xmath278 and @xmath279 , respectively .
the pertinent bifurcation scenarios were described in figures 17 and 18 of @xcite .
the `` simple '' scenario taking place at @xmath278 consists of two orbits being generated in the rotational symmetry - breaking at @xmath280 , followed by a tangent bifurcation destroying one of them and a third orbit . to smooth this bifurcation peak , a uniform approximation describing the complete scenario must be found , which requires the construction of a suitable normal form .
although a uniform approximation for the symmetry - breaking is available @xcite , the derivation of the pertinent normal form relies on principles different from the catastrophe theory classification used here , and it is not clear how these two can be united into a single normal form .
thus , the construction of a uniform approximation for this bifurcation scenario , and even more so for the more complicated scenario at @xmath279 , remains an open problem to be solved in the future .
it can be solved within the framework of uniformization presented in section  [ ssec : unifconst ] , but will require a novel way of constructing normal forms .
the approach to high - resolution semiclassical quantization relies on the harmonic inversion of a fourier transformed semiclassical spectrum , i.e. of a recurrence spectrum .
the above method of uniformizing the bifurcation - induced excessively high recurrence peaks in a semiclassical spectrum would therefore , if it could be implemented systematically , also pave the way for the inclusion of uniform approximations into a high - resolution semiclassical quantization , which has not been possible so far .
we were able to demonstrate the feasibility of our method by way of example for the hydrogen atom in an electric field @xcite , which is less demanding classically .
its application to the crossed - fields hydrogen atom , however , remains open for future work .
for the first time , a high - resolution semiclassical quantization of the hydrogen atom in crossed electric and magnetic fields has been presented .
it achieved the identification of the strong spectral lines in different @xmath1-manifolds . by means of a detailed semiclassical analysis of the pertinent quantum spectrum
, it was shown that bifurcations of closed orbits play a crucial role in the semiclassical spectrum and preclude the resolution of finer details in the semiclassical spectrum .
they pose a particular challenge to the semiclassical quantization because they require a special treatment by uniform approximations .
a simple heuristic scheme for the construction of uniform approximations has been proposed .
its simplicity and efficacy was demonstrated by a derivation of the uniform approximations for the codimension - one generic bifurcations of closed orbits .
we have devised a general method for the inclusion of uniform approximations in a high - resolution semiclassical quantization by harmonic inversion . in a recent publication @xcite it was successfully applied to the hydrogen atom in an electric field . in the case of the crossed - fields hydrogen atom ,
the diversity and complexity of the bifurcation scenarios encountered so far hinders the systematic implementation of the uniformization procedure .
the treatment of all relevant bifurcations and the calculation of a detailed semiclassical spectrum thus remain challenging tasks for future studies .