a comprehensive description of electromagnetism in classical theory is given in terms of the electromagnetic field strength acting locally and directly on charged matter .
this is not the case in quantum theory where a state of charged matter is influenced by electromagnetic field even in the situation when the region of the nonvanishing field strength does not overlap with the region accessible to charged matter ; the indispensable condition is that the latter region be non - simply - connected .
in particular , a magnetic field of an infinitely long current - carrying solenoid can be shielded from the penetration of charged matter ; such a field configuration may be denoted as an impenetrable magnetic vortex .
perhaps , w.  ehrenberg and r.  e.  siday were the first to discover theoretically the purely quantum effect which is due to the influence of the vector potential of the impenetrable magnetic vortex on the motion of a charged particle , and they proposed an interference experiment to detect this @xcite .
ten years later the effect was rediscovered by y.  aharonov and d.  bohm who made an important addition by predicting also an effect which is due to the influence of the electric scalar potential @xcite . immediately , in a year , the first interference experiment to verify the magnetic effect was performed @xcite , and this caused a burst of interest opening a long way of an extensive research both in theory and in experiment , see , e.g. , @xcite and references therein ; the effect conventionally bears the name of the authors of @xcite .
although a great progress in understanding the aharonov  bohm ( ab ) effect has been achieved since the time of foundational works @xcite , one may agree with the authors of @xcite that `` the investigation and exploitation of the ab effect remain far from finished '' .
it should be noted that the concern of paper @xcite was in the elaboration of the consistent theory of refraction of electron rays in the framework of electron wave optics .
following thoroughly their flow of reasoning , the authors arrive at the conclusion that the medium in a spatial region outside a magnetic vortex is anisotropic and that there exist wave - optical phenomena which are due to the influence of the vector potential of the vortex . finally , in the very end of the paper , they suggest an interference experiment to verify their findings .
an approach of paper @xcite was quite different : the authors heuristically , with the use of few rather general and obvious formulas ( which are even left unnumbered in the text ) , arrive at their suggestion of interference experiments already in the beginning of the paper .
a powerful theoretical background behind this intuitive suggestion is hereafter given in section 4 .
`` exact solution for scattering problems '' ( which exceeds in its volume two preceding sections ) .
the background relies directly on the formalism of quantum mechanics , namely on the schrdinger wave equation , which allows the authors to consider their effect as a scattering event .
perhaps , this was a reason , why paper @xcite , unlike paper @xcite , had received an immediate response from the physics community .    in this respect
it looks rather paradoxical that , whereas interference experiments aiming at the verification of the ab effect were numerous ( see @xcite ) , a scattering experiment aiming at the same was _ never _ performed .
one of the motivations of the present paper is to provide an explanation for this paradox .
let us commence by recalling the famous formula obtained by the authors of @xcite for the scattering differential cross section per unit length of the magnetic vortex ( see ( 21 ) and ( 22 ) in @xcite ) : @xmath2 } { \sin^2(\varphi/2)},\label{eq1}\ ] ] where @xmath3 is the momentum of a scattered particle , @xmath4 is its charge , @xmath5 is the scattering angle ( with @xmath6 corresponding to the forward direction ) and @xmath7 is the flux of the magnetic vortex directed along the @xmath8-axis . a characteristic feature of ( 1 ) is a periodic dependence in the magnetic flux with the period equal to the london flux quantum , @xmath9 .
also , one could suggest that ( 1 ) describes a purely quantum effect , since it vanishes in the formal limit of @xmath10 , or , more physically , in the limit of large enough values of the particle momentum , @xmath11 .
however , result ( 1 ) corresponds to the idealized case of the magnetic vortex of the vanishing transverse size .
since any real magnetic vortex ( either that produced by a current - carrying solenoid , or that inside a magnet ) is of a nonzero transverse size , the quotient of squared sines in ( 1 ) in the real case becomes a certain function of dimensionless variable @xmath12 , where @xmath13 is a length characterizing the transverse size of the vortex .
the transverse size effects were taken into account in @xcite , and results were found to converge with ( 1 ) in the case of @xmath14 . as to the case of @xmath15 ,
it was only scarcely and insufficiently considered , see @xcite , and , therefore , let us dwell on this case at more length .
it should be noted in the first place that the smallest possible values of @xmath13 can be of order 1 @xmath16 ( ferromagnetic whiskers enclosed in screening shells ) @xcite , while , in a scattering event , the largest possible wavelengths , @xmath17 , of the lightest charged particle , electron , are of order 0.1 nm ( slowly moving electrons of energies 10 - 100 ev , e.g. , such as those in famous diffraction experiments by c.  j.  davisson and l.  h.  germer in 1927 , nobel prize in 1937 ) .
hence , quantity @xmath18 takes values of order @xmath19 at least , and formula ( 1 ) which might be supposed to approximate the case of @xmath20 has no relation to scattering of real particles .
albeit this formula looks essentially quantum , the `` quantum '' factor of @xmath21 appears just due to dimensional reasons because of @xmath22 ; it is a somewhat more explicating to write @xmath21 as product @xmath23 , where the second ( dimensionless ) factor is vanishingly small for all charged particles existing in nature .
this was certainly understood by the authors of @xcite , who never proposed cross section ( 1 ) for the experimental verification .
it may seem quite reasonable to apply the quasiclassical ( called also semiclassical ) method of g.  wentzel , h.  a.  kramers and l.  brillouin ( wkb ) @xcite to solving the problem of scattering of real charged particles by an impenetrable magnetic vortex : the size of the interaction region may be estimated as being of order of @xmath13 and larger , then condition @xmath24 allows one to hope that the use of the wkb method will be adequate , see , e.g , @xcite .
however , a skepticism immediately arises : how can a purely and essentially quantum effect as is the ab one be deduced with the help of a quasiclassical ( semiclassical ) method which is commonly believed to give nothing else but the yield of classical theory ? in this respect it should be noted that the findings of the authors of @xcite indicate that scattering off an impenetrable magnetic vortex in the case of @xmath15 is reduced to classical scattering of point particles off an impenetrable tube with no flux .
nevertheless , in spite of the skepticism substantiated by the results of @xcite , we shall proceed in the present paper by applying the wkb method to the solution of the scattering ab problem .
let us recall that interference experiments aiming at the verification of the ab effect were performed with rather energetic electrons .
we are considering electrons of the same ( as in the interference experiments ) and even smaller ( down to 10 - 20 ev ) energies and study direct scattering of such particles on an impenetrable magnetic vortex , using the wkb method . whereas the authors of @xcite employed the dirichlet boundary condition to ensure the vortex impenetrability , we go further by employing the most general ( unorthodox in the language of @xcite ) boundary condition that is required by the basic concepts of probability conservation and self - adjointness of the hamilton operator . in the course of our study we shall find , in addition to the classical effect of @xcite , a purely quantum effect that can be denoted as the scattering ab effect with real particles .    in the next section
we consider a quantum - mechanical charged particle in the background of an impenetrable magnetic vortex and choose a boundary condition for the particle wave function at the vortex edge .
the wkb method is used to obtain the partial waves of the solution to the schrdinger equation in section 3 .
the scattering amplitude and cross section are obtained in section 4 .
the conclusions are drawn and discussed in section 5 .
we relegate some details of calculation of the scattering amplitude to appendix a , while the use of the hamilton - jacobi equation of motion to derive the classical effect is outlined in appendix b.
defining a scalar product as @xmath25 , we get , using integration by parts ,    @xmath26\chi\right\}= \nonumber \\ = \!\int\limits_{\omega}\!{\rm d}^3x\left\{\left [ -\frac{1}{2m}\left(\hbar{{\mbox{\boldmath $ \partial$}}}-{\rm i}\frac{e}{c}{\bf a}\right)^2\!+\!v\right]\tilde{\chi}\right\}^*\chi- \frac{\hbar}{2m}\int\limits_{\partial\omega}\!{\rm d}{{\mbox{\boldmath $ \sigma$}}}\cdot \left\{\tilde{\chi}^*\left[\left(\hbar{{\mbox{\boldmath $ \partial$}}}-{\rm i}\frac{e}{c}{\bf{a}}\right)\chi\right]\right.\!-\nonumber \\
\left.-\left[\left(\hbar{{\mbox{\boldmath $ \partial$}}}+{\rm i}\frac{e}{c}{\bf{a}}\right)\tilde{\chi}^*\right]\chi\right\}\equiv \left(h^\dag\tilde{\chi},\chi\right)-{\rm i}\hbar\int\limits_{\partial\omega}{\rm d}{{\mbox{\boldmath $ \sigma$}}}\cdot{\bf u } , \label{eq2}\end{aligned}\ ] ]    where @xmath27 is a two - dimensional surface bounding the three - dimensional spatial region @xmath28 , @xmath29 is the formal expression for a general quantum - mechanical hamilton operator in an external electromagnetic field ( electric scalar potential is included in the general potential energy part @xmath30 ) , and @xmath31 -\left[\left(\hbar{\mbox{\boldmath $ \partial$}}+{\rm i}\frac{e}{c}{\bf a}\right)\tilde{\chi}^*\right]\chi\right\}.\label{eq4}\ ] ] operator @xmath32 is hermitian ( symmetric ) , @xmath33 if @xmath34 the latter condition can be satisfied in various ways by imposing different boundary conditions for @xmath35 and @xmath36 .
however , among the whole variety , there may exist a possibility that a boundary condition for @xmath36 is the same as that for @xmath35 ; then the domain of definition of @xmath37 ( set of functions @xmath36 ) coincides with that of @xmath32 ( set of functions @xmath35 ) , and operator @xmath32 is called self - adjoint .
whether such a possibility exists is in general determined by the weyl
von neumann theory of self - adjoint operators , see , e.g. , @xcite .
the action of a self - adjoint operator results in functions belonging to its domain of definition only , and , therefore , a multiple action and functions of such an operator ( for instance , the evolution operator ) can be consistently defined .    in the present case ,
the problem of self - adjointness of operator @xmath32 ( 3 ) is resolved by imposing a boundary condition in the form @xmath38\biggr|_{{\bf
x}\in\partial\omega}= \nonumber \\
\biggl.=\left[\sin(\rho\pi)\,{\bf n}\cdot\left({\mbox{\boldmath $ \partial$}}-\frac{{\rm i}e}{\hbar c}{\bf a}\right)\tilde{\chi}+ \cos(\rho\pi)\frac{1}{\bf n\cdot{\bf x}}\tilde{\chi}\right]\biggr|_{{\bf x}\in\partial\omega}=0 , \label{eq7}\end{aligned}\ ] ] where @xmath39 is the internal unit normal to boundary @xmath27 and @xmath40 is the real parameter ( the self - adjoint extension parameter ) . defining current @xmath41- \left[\left(\hbar{\mbox{\boldmath $ \partial$}}+{\rm i}\frac{e}{c}{\bf a}\right)\chi^*\right]\chi\right\},\label{eq8}\ ] ]
we note that self - adjointness condition ( 7 ) results in the vanishing of the normal component of the current at the boundary : @xmath42    if @xmath35 is a solution to the schrdinger equation , @xmath43 then current ( 8) is the probability current density satisfying the continuity equation @xmath44 thus , self - adjointness of the hamilton operator ensures probability conservation , @xmath45 .    in the following we consider a quantum - mechanical charged particle in the background of a static magnetic field , hence @xmath46 .
the magnetic field is confined to an infinite tube which is classically impenetrable to the charged particle ; thus , spatial region @xmath28 is an infinite three - dimensional space with the exclusion of the tube containing the magnetic flux lines , and boundary surface @xmath27 is a surface of the tube . assuming the cylindrical symmetry of the flux tube , we use cylindrical coordinates @xmath47 with the @xmath8-axis coinciding with the axis of the tube and choose the components of the vector potential in @xmath28 in the form @xmath48 where @xmath7 is the total magnetic flux confined in the tube .
then the hamilton operator takes form @xmath49,\label{eq13}\ ] ] and the solution to schrdinger equation ( 10 ) with @xmath32 ( 13 ) is @xmath50 where @xmath51 is the solution to equation @xmath52\psi(r,\varphi)=0,\label{eq15}\ ] ] and obeys the boundary condition , see ( 7 ) , @xmath53|_{r = r_c}=0;\label{eq16}\ ] ] @xmath13 is the tube radius , @xmath54 is the energy of the stationary scattering state , @xmath55 and @xmath3 are the momenta in the longitudinal and radial transverse directions , respectively
. one could recognize that ( 16 ) is known as the robin boundary condition ; the case of @xmath56 corresponds to the dirichlet condition ( perfect reflectivity of the boundary ) and the case of @xmath57 corresponds to the neumann condition ( absolute rigidity of the boundary ) .
it should be emphasized that parameter @xmath40 is in general dependent on @xmath8 and @xmath5 .
thus , the `` number '' of self - adjoint extension parameters is infinite , moreover , it is not countable but is of power of a continuum .
this distinguishes the case of an extended boundary from the case of an excluded point ( contact interaction ) when the number of self - adjoint extension parameters is finite , being equal to @xmath58 for the deficiency index equal to \{@xmath59 } ( see , e.g. , @xcite ) .
the motion of the particle in the longitudinal direction is free , and our aim is to determine the particle motion in the orthogonal plane , i.e. to find a solution to ( 15 ) and ( 16 ) .
let us start with the decomposition of wave function @xmath51 into partial waves @xmath60 where @xmath61 is the set of integer numbers .
the radial component of each partial wave satisfies equation @xmath62r_n(r)=0,\label{eq18}\ ] ] where @xmath63 and @xmath64 @xmath65 .
note that the amplitude of @xmath51 ( 17 ) is periodic in the value of flux @xmath7 with the period equal to the london flux quantum .
we present the radial component of the partial wave as @xmath66 $ ] , where @xmath67 satisfies the ricatti - type equation @xmath68 expanding @xmath67 into series @xmath69 and collecting terms of the same order in @xmath70 , we obtain the system of equations : @xmath71 ^ 2=p_n^2 , \\    2\left[\partial_r\sigma_n^{(0)}\right]\left[\partial_r\sigma_n^{(1)}\right]=-pr_c\delta_r\sigma_n^{(0 ) } , \\
2\left[\partial_r\sigma_n^{(0)}\right]\left[\partial_r\sigma_n^{(2)}\right]=-pr_c\delta_r\sigma_n^{(1)}- \left[\partial_r\sigma_n^{(1)}\right]^2 , \\    \qquad\qquad \vdots
\\     2\left[\partial_r\sigma_n^{(0)}\right]\left[\partial_r\sigma_n^{(2l)}\right]=-pr_c\delta_r\sigma_n^{(2l-1)}- 2\sum\limits_{l'=1}^{l-1}\left[\partial_r\sigma_n^{(l')}\right]\left[\partial_r\sigma_n^{(2l - l')}\right]- \left[\partial_r\sigma_n^{(l)}\right]^2,\\    2\left[\partial_r\sigma_n^{(0)}\right]\left[\partial_r\sigma_n^{(2l+1)}\right]=-pr_c\delta_r\sigma_n^{(2l)}- 2\sum\limits_{l'=1}^{l}\left[\partial_r\sigma_n^{(l')}\right]\left[\partial_r\sigma_n^{(2l+1-l')}\right ] , \\
\qquad\qquad \vdots    \end{array } \label{eq22}\ ] ] consecutively solving equations beginning from the top , one can obtain @xmath67 ( 21 ) up to an arbitrary order in expansion parameter @xmath72 .
the terms corresponding to odd @xmath73 contribute to the amplitude of @xmath74 , while the terms corresponding to even @xmath73 contribute to the phase of @xmath74 .
however , the power of the wkb method reveals fully itself in those cases when it suffices to take account for the first two terms , @xmath75 and @xmath76 , only . it should be noted that equation ( 18 ) , or ( 20 ) , has its analogue in optics , describing the propagation of waves in media @xcite : quotient @xmath77 corresponds to the index of refraction .
terms @xmath78 @xmath79 are negligible , if the variation of the refraction index along the ray of propagation becomes appreciable on the distances which are much larger than the effective wavelength , @xmath80 where @xmath39 points in the direction of the ray which in our case is orthogonal to the boundary . in view of ( 19 ) this yields condition @xmath81 which is satisfied in the quasiclassical region , @xmath82 as well as in the deeply nonclassical region @xmath83 solving the first two equations in ( 22 ) , we get @xmath84.\label{eq27}\ ] ] it should be noted that the second equation in ( 22 ) is nothing more but the continuity equation , @xmath85 , for the partial wave current , @xmath86,\ ] ] in the quasiclassical region .    actually , the wkb approximation ( i.e. the neglect of @xmath78 with @xmath87 ) is efficient in a much more extended regions than those given by ( 25 ) and ( 26 ) .
the wkb approximation is not valid in the vicinity of @xmath88 which is defined by @xmath89 this point is the turning point of a classical trajectory of the particle . in view of ( 27 ) , the two linearly independent solutions to ( 18 ) in the quasiclassical region are @xmath90 while the two linearly independent solutions to ( 18 ) in the deeply nonclassical region are @xmath91 where @xmath92 the crucial question is how oscillating solutions ( 29 ) in the outer region go over to solutions ( 30 ) in the inner region and vice verse . to answer this question , one has to find a solution in the vicinity of the turning point , where @xmath93 and @xmath94 , and to match it with the solutions in the inner and the outer regions .
but a more instructive way , as stated in @xcite , is to consider formally @xmath74 as a function of complex variable @xmath95 and to perform transitions along paths in the complex plane , where condition ( 24 ) is satisfied .
namely , a path goes along the real positive semiaxis up to the left vicinity of the turning point , circumvents it along a semicircle in the upper or lower half - plane , and then goes again along the real positive semiaxis . for the transition from the inner to the outer region
we have @xmath96 where the ( @xmath97 ) sign corresponds to the path in the upper or lower half - plane , and @xmath98 @xmath99 , in contrast to @xmath100 , is divergent in the inner region . for the transition from the outer to the inner region we have @xmath101 and @xmath102 \rightarrow r_{n,{\rm in}}^{(-)}(r),\label{eq35}\ ] ] where the account is taken for the fact that @xmath103 is real .
boundary condition ( 16 ) in terms of partial waves takes form @xmath104|_{r = r_c}=0.\label{eq36}\ ] ] in the case @xmath105 we choose @xmath106e^{{\rm i}\pi/4},\label{eq37}\ ] ] where the factor of @xmath107 is inserted for future convenience .
coefficient @xmath108 is determined from condition ( 36 ) : @xmath109|_{r = r_c } } { \cot(\rho\pi)+[r\partial_r\ln r_{n,{\rm out}}^{(+)}(r)]|_{r = r_c}}= \nonumber \\
= \exp\left\{-\frac{2{\rm i}}{\hbar}\int\limits_{r_t}^{r_c}{\rm d}r\,p_n-2{\rm i}\arctan\left [ \frac{r_cp_n(r_c)\hbar^{-1}}{\cot(\rho\pi)-\frac{1}{2}\,\frac{p^2}{p_n^2(r_c)}}\right]\right\ } , \qquad r_c > r_t ; \label{eq38}\end{aligned}\ ] ] note that the explicit form of the turning point ( i.e. the solution to ( 28 ) ) is evidently @xmath110 . in the case
@xmath111 we choose @xmath112,\label{eq39}\ ] ] where coefficient @xmath108 is determined from condition ( 36 ) : @xmath113 @xmath114\right|_{r = r_c } } { \cot(\rho\pi)+\left.\left\{r\partial_r\ln \left[r_{n,{\rm in}}^{(-)}(r)+cr_{n,{\rm in}}^{(+)}(r)\right]\right\ } \right|_{r = r_c } } , \,\,r_c < r_t,\ ] ] and constant @xmath115 is to be fixed before long . using correspondence rules ( 33 ) and ( 35 ) , we continue @xmath74 ( 39 ) to the case @xmath116 as @xmath117 @xmath118 + c\left[r_{n,{\rm out}}^{(+)}(r)e^{-{\rm i}\pi/4}+ r_{n,{\rm out}}^{(-)}(r)e^{{\rm i}\pi/4}\right]\right\}.\ ] ] the converging wave , @xmath119 , should disappear from the part with coefficient @xmath108 , cf .
this fixes the constant , @xmath120 , yielding in the case @xmath116 @xmath121 with @xmath122^{-1},\,\,r_c < r_t.\label{eq41}\end{aligned}\ ] ] we can rewrite ( 37 ) in the case @xmath105 as @xmath123r_{n,{\rm out}}^{(+)}(r)e^{{\rm i}\pi/4}.\label{eq42}\ ] ] we conclude this section by stating that ( 40 ) and ( 42 ) give the radial components of partial waves of @xmath51 ( 17 ) in the quasiclassical region in the wkb approximation .
going over to asymptotics @xmath124 , where @xmath125,\qquad \int\limits_{r_t}^{r}drp_n = pr-\frac{\hbar}{2 }    @xmath126 where @xmath127 and @xmath128+\right . \nonumber \\
\left.+\sum\limits_{|n-\mu|>pr_c/\hbar}e^{{\rm i}n\varphi}a_ne^{-\frac{\rm i}{2}|n-\mu|\pi } c_n(r_c,\rho)\right\}.\label{eq46}\end{aligned}\ ] ] to determine coefficient @xmath129 , let us recall the asymptotics of the partial wave decomposition of a plane wave , @xmath130= \nonumber \\
= \sqrt{\frac{2\pi\hbar}{pr}}\left[\delta(\varphi)e^{{\rm i}(pr\hbar^{-1}-\pi/4 ) } + \delta(\varphi-\pi)e^{-({\rm i}pr\hbar^{-1}-\pi/4)}\right],\label{eq47}\end{aligned}\ ] ] where @xmath131 is the angular delta - function , @xmath132 .
the asymptotics of the plane wave is naturally interpreted as a superposition of two cylindrical waves : the diverging one , @xmath133 , going in the forward , @xmath6 , direction and the converging one , @xmath134 , coming from the backward , @xmath135 , direction . namely the converging cylindrical wave should be present without distortions in the asymptotics of wave function @xmath51 ( 44 ) , whereas the diverging cylindrical wave is distorted and differs from that in ( 47 ) ;
actually , this is the condition that @xmath51 be the scattering state solution . equating terms before @xmath134 in @xmath51 to the terms before @xmath136 in the first line of ( 47 ) , we get @xmath137.\label{eq48}\ ] ] as a result , we obtain the following expression for the @xmath13-independent part of the wave function @xmath138 where @xmath139}=\nonumber \\ = \sqrt{\frac{2\pi\hbar}{pr}}\left[e^{-{\rm i}\mu{\rm sgn}(\varphi)\pi } \delta(\varphi)e^{{\rm i}(pr\hbar^{-1}-\pi/4)}+\delta(\varphi-\pi ) e^{-{\rm i}(pr\hbar^{-1}-\pi/4)}\right],\label{eq50}\end{aligned}\ ] ] is the incident wave and @xmath140\!]+\frac 12\right)\varphi } } { \sin(\varphi/2)}\label{eq51}\ ] ] is the scattering amplitude @xcite ; the sign function is @xmath141 at @xmath142 , @xmath143\!]$ ] denotes the integer part of quantity @xmath144 ( i.e. the integer which is less than or equal to @xmath144 ) and it is implied that @xmath145 .
the squared absolute value of ( 51 ) yields differential cross section ( 1 ) .
incident wave ( 50 ) differs from the plane wave : the distortions are due to the long - range nature of the vector potential outside the flux tube .
the apparent divergence of amplitude ( 51 ) in the forward direction is a phantom , because ( 51 ) is valid under condition @xmath146 , whereas , otherwise , the divergence is absent and the discontinuity in incident wave ( 50 ) at @xmath147 is cancelled , yielding @xmath148 which is continuous and differentiable at @xmath6 , see , e.g. , @xcite and references therein .
amplitude @xmath149 ( 51 ) is of order @xmath150 $ ] and is negligible as compared to the contribution from the @xmath13-dependent part of the wave function , which can be presented in the following form : @xmath151 \frac{e^{{\rm i}(pr\hbar^{-1}+\pi/4)}}{\sqrt{r}},\label{eq52}\ ] ] where @xmath152 @xmath153 @xmath154 we show in appendix a ( see ( a.22 ) ) that amplitude @xmath155 ( 55 ) , although exceeding amplitude @xmath149 ( 51 ) , is still negligible , being of order @xmath156 $ ] .
amplitude @xmath157 ( 54 ) is of nonnegligible value which is calculated in appendix a : @xmath158-{\rm i}\pi/4\right\}\times\nonumber \\
\times\exp\left\{-2{\rm i}\arctan\left[\frac{2pr_c\hbar^{-1}|\sin^3(\varphi/2)| } { 2\cot(\rho\pi)\sin^2(\varphi/2)-1}\right]\right\}.\label{eq56}\end{aligned}\ ] ] amplitude @xmath159 ( 53 ) is a finite sum of geometric progression , which is straightforwardly calculated : @xmath160\!]+\frac 12)\varphi } \frac{\sin(s_c\varphi/2)}{\sin(\varphi/2)}\cos(\mu\pi+s_c\varphi/2)\label{eq57}\ ] ] in the case @xmath161\!]-[\![\mu]\!]=[\![pr_c\hbar^{-1}-\mu]\!]+[\![\mu]\!]+1=s_c,\label{eq58}\ ] ] or @xmath162\!]+\frac 12\mp\frac 12)\varphi } \left\{\frac{\sin[(s_c+1/2)\varphi/2]}{\sin(\varphi/2)}\cos\left[\mu\pi+(s_c+1/2)\varphi/2\right]-\right .
\nonumber \\\left.-\tan(\varphi/4)\pm{\rm i}\right\}\label{eq59}\end{aligned}\ ] ] in the case @xmath161\!]-[\![\mu]\!]-\frac
12\pm\frac 12=[\![pr_c\hbar^{-1}-\mu]\!]+[\![\mu]\!]+ \frac 12\mp\frac 12=s_c.\label{eq60}\ ] ] since @xmath163 , amplitude @xmath159 is strongly peaked in the forward direction ( actually as a smoothed delta - function ) , whereas amplitude @xmath157 is vanishing in the forward direction .
therefore , the interference between the amplitudes , @xmath164 , is vanishing in the same manner as @xmath165 .
the differential cross section of the scattering process , @xmath166 , is hence a sum of two terms : @xmath167 and @xmath168,\label{eq62}\ ] ] where we have recalled that @xmath65 and introduced function @xmath169 which at @xmath170 can be regarded as a regularized ( smoothed ) delta - function , @xmath171    term ( 61 ) which is independent of the self - adjoint extension parameter , the magnetic flux and even the particle momentum is obtainable in the framework of classical theory , for instance , with the use of the hamilton - jacobi equation of motion , see appendix b. amplitude @xmath157 ( 54 ) can be rewritten as @xmath172\right\},\label{eq64}\ ] ] where @xmath173 is the wkb phase shift , @xmath174 and @xmath175 are given by ( a.2 ) and ( a.3 ) in appendix a. the derivative of the wkb phase shift , @xmath176 , can be related to the classical reflection angle ( equal to the classical incidence angle ) , @xmath177 , see ( b.10 ) in appendix b : @xmath178 thus , the condition of stationarity of the phase in ( 64 ) , @xmath179 ( @xmath180 is negligible , see appendix a ) , coincides in its form with the classical relation between scattering angle @xmath5 and reflection angle @xmath177 , see ( b.12 ) in appendix b , and the classical impact parameter can be introduced as @xmath181 note also that the first equation in ( 22 ) is equivalent to the equation determining the squared @xmath95-derivative of the hamilton - jacobi action , see ( b.5 ) in appendix b.    however , the merits of the wkb method are not restricted to determining the wkb phase shift , i.e. to the only description of classical elastic reflection .
the use of the method has also allowed us to find term ( 62 ) which is generically independent of the self - adjoint extension parameter but depends on both the magnetic flux and the particle momentum .
this term describes the purely quantum effect of diffraction .
it should be noted that the quantum effect is well separated in the scattering angle from the classical one : namely , diffraction is in the forward direction where classical reflection is absent .
although the range of angles where diffraction is to be observed is quite narrow , the whole ( i.e. integrated over @xmath145 ) contribution of the diffraction effect to the total cross section is the same as that of the reflection effect : @xmath182 and the total cross section is twice the classical one ; the latter is needed for the optical theorem to be maintained ( see @xcite ) .
it should be noted in the first place that the ab scattering amplitude ( see ( 51 ) ) and , consequently , cross section ( 1 ) were obtained in @xcite by imposing the condition of regular behavior for the particle wave function at the location of a singular ( i.e. of zero transverse size ) magnetic vortex .
however , the regularity condition is not the most general one that is required by probability conservation .
the most general condition in the case when a spatial dimension along the vortex is ignored involves four self - adjoint extension parameters ( the deficiency index is \{2,2 } ) @xcite .
the dependence on the self - adjoint extension parameters enters into terms which are added to amplitude @xmath149 ( 51 ) .
therefore , the cross section depends on four arbitrary parameters or , if a longitudinal dimension is taken into account , on four arbitrary functions of @xmath8 .
if an invariance under rotations around the vortex is assumed to be a physical requirement , then the cross section still remains to be dependent on two arbitrary functions .
the problem of seeking some , if any , sense for this arbitrariness is of purely academic interest , since the case of a magnetic vortex of zero transverse size has no relation to physics reality .
turning now to the physical case of a magnetic vortex of nonzero transverse size , which is addressed in the present paper , we would like to emphasize that the classical initial equation ( see ( b.2 ) in appendix b ) , as well as the quantum - mechanical one ( see ( 10 ) with ( 13 ) ) , involves the vector potential of the vortex ( see ( 12 ) ) .
however , in classical theory , this potential has no physical consequences , as is demonstrated , for example , in appendix b. one may agree with the authors of @xcite that `` the ab effect was already implicit in the 1926 schrdinger equation '' , but it takes more than eight decades before a physical consequence of the vector potential from the schrdinger equation is at last pointed out , see ( 62 ) .
this explains why a scattering experiment aiming at the verification of the ab effect is still lacking .
the main motivation of the present work is to make such an experiment plausible from a theorist s point of view .
that is why we have taken care to prevent from any objection of theoretical kind : we have employed the most general boundary condition that is required by probability conservation and have shown how to grasp the ab effect even with the use of the 1926 wkb method .    the ab effect , as well as quantization of the flux trapped in superconductors , the van hove singularities in the excitation spectra of crystals etc , explicates an importance of topological concepts for quantum physics @xcite .
the ab phase , @xmath183 , which is acquired by the charged - particle wave function after the particle has encircled a magnetic vortex is certainly the same for the particles of different energies .
nonetheless , a question always is how to observe this phase in real experiments , and what are physical restrictions on the particle energy in these experiments .
for instance , interference experiments to verify the ab effect involve electrons of energies restricted to the range of 10 - 100 kev @xcite .
interference is immanently a wave - optical phenomenon which is indirectly described in quantum theory through reference to wave optics ( a quantum - mechanical particle possesses de broglie s wavelength and thus is subject to the laws of wave optics )
. meanwhile , quantum theory provides a direct description of scattering phenomena and , in particular , of the ab effect as a scattering event . a scattering experiment with particles of the wavelength exceeding
the vortex radius is hardly realistic , and that is why we are considering a scattering event with particles of the wavelength which is less than the vortex radius .    by using the wkb method in the present paper ,
we obtain the scattering differential cross section , @xmath166 , consisting of two parts : one ( 61 ) describes the classical effect of elastic reflection from the vortex and another one ( 62 ) describes the quantum effect of diffraction on the vortex .
the effects are equal to each other , when integrated over the whole range of the scattering angle , see ( 69 ) .
it should be noted that , although the amplitude yielding the classical effect depends on the self - adjoint extension parameter ( @xmath40 ) , the flux of the vortex ( @xmath7 ) and the particle momentum ( @xmath3 ) , see ( 56 ) , all the dependence disappears in the differential cross section , see ( 61 ) ; certainly , this cross section can be obtained with the use of purely classical methods ( the hamilton - jacobi equation of motion ) .
the amplitude yielding the quantum effect is generically independent of the self - adjoint extension parameter , see ( 57)-(60 ) .
this signifies that diffraction ( to be more precise , the fraunhofer diffraction , i.e. the diffraction in almost parallel rays ) is independent of the choice of different boundary conditions which are compatible with the conservation of probability . owing to diffraction the ab effect
persists in scattering of real particles , e.g. , electrons of energies in the range from 10 ev to 100 kev .
the classical effect opens a window in the forward direction , where the ab effect is to be observed as a fringe shift in the diffraction pattern exhibited by the scattering differential cross section , see ( 62 ) .
we present the window for the observation of the ab effect on fig .
1 , where the scattering differential cross section and the scattering angle are normalized in such a way that the plot remains actually the same for @xmath184 ; recall that the present - day realistic restriction is @xmath185 .
the classical effect given by ( 61 ) is unobservable at @xmath186 , whereas the diffraction oscillations given by ( 62 ) are unobservable on the classical background at @xmath187^{1/3}<|\varphi|<\pi$ ] ; at @xmath188^{1/3}$ ] , both effects are unobservable being indistinguishable from the @xmath189-background which is of order of @xmath190^{1/3}$ ] , see ( a.22 ) in appendix a. a gate for particles propagating in the strictly forward ( @xmath6 ) direction is opened when the vortex flux equals an integer times the london flux : more than 90% of the diffraction cross section is accumulated in a peak centred at @xmath6 and having width @xmath191 , see the dashed line on fig .
the gate for the strictly forward propagation of particles is closed when the vortex flux equals a half - integer times the london flux : a dip at @xmath6 is surrounded by two symmetric peaks accumulating more than 85% of the diffraction cross section ( each peak having width @xmath191 ) , see the solid line on fig .
thus , the classically forbidden propagation of particles in the strictly forward direction is controlled by altering the amount of the vortex flux .
is along the ordinate axis and @xmath192 is along the abscissa axis ( @xmath193 ) . dashed and solid lines correspond to cases @xmath194 and @xmath195 , respectively ( @xmath196 ) .
the area under one central peak by the dashed line ( @xmath197 ) is 0.4514119@xmath970.0000002 , and the area under two central peaks by the solid line ( @xmath198 ) is 0.4278549@xmath970.0000013.,title="fig : " ] +    this scattering ab effect was never observed experimentally , since a detailed proposal for its observation was just recently elaborated @xcite .
such an experiment is in our opinion of fundamental significance and seems to be quite feasible with the present - day facilities providing for sources of bright and yet coherent electron beams .
the work of yu .
a.  s. was supported by the national academy of sciences of ukraine ( project no.0112u000054 ) , by the the program of fundamental research of the department of physics and astronomy of the national academy of sciences of ukraine ( project no.0112u000056 ) and by the ictp  seenet - mtp grant prj-09 `` strings and cosmology '' . the work of n.  d.  v. was supported by the science and technology center in ukraine and the national academy of sciences of ukraine within the framework of the targeted research & development initiatives ( trdi ) program under grant no.5716 .
we acknowledge the support from the state agency for science , innovations and informatization of ukraine under the sffr - brrfr grant f54.1/019 .
we present amplitude @xmath157 ( 54 ) in the form @xmath199\right\},\end{aligned}\ ] ] where , see ( 38 ) , @xmath200 @xmath201,\ ] ] and we have introduced dimensionless variable @xmath202 . to calculate the sum in ( a.1 ) in asymptotics @xmath203
, we use the poisson summation formula @xmath204\right\}+\frac
12e^{{\rm i}\eta(s_+,s)}+ \frac 12e^{{\rm i}\eta(-s_-,s)},\ ] ] where @xmath205\!]$ ] .
if @xmath206 and @xmath207 is convex upwards , @xmath208 , on interval @xmath209 , then only a finite number of terms in the series on the right - hand side of ( a.4 ) contributes to the leading asymptotics at @xmath203 , and one can use the method of stationary phase for its evaluation .
namely , if equation @xmath210\right|_{n = n_j}-2\pi l_j=0\ ] ] determines a stationary point inside the interval , @xmath211 , for some values of @xmath73 denoted by @xmath212 , then @xmath213\right\}\left\{\frac{2\pi e^{-{\rm i}\pi/2 } } { \left.-\left[\frac{\partial^2}{\partial n^2}\eta(n , s)\right]\right|_{n = n_j}}\right\}^{1/2}+o(1).\ ] ]    in the present case we have @xmath214.\ ] ] differentiating ( a.2 ) and ( a.3 ) over @xmath215 , we obtain @xmath216 and @xmath217^{1/2 } \left\{\left[1-\left(\frac{n-\mu}{s}\right)^2\right]\cot(\rho\pi)-\frac 32\right\ } } { \left\{\left[1-\left(\frac{n-\mu}{s}\right)^2\right]\cot(\rho\pi)-\frac 12\right\}^2+s^2\left [ 1-\left(\frac{n-\mu}{s}\right)^2\right]^3}.\ ] ] ( a.9 ) is negligible as compared to ( a.8 ) , being of order @xmath218 at @xmath219 , @xmath203 , where the wkb approximation is valid , see ( 25 ) . therefore equation ( a.5 )
takes form @xmath220 choosing the range for @xmath5 as @xmath145 , we find that a solution to ( a.10 ) exists at @xmath221 only .
denoting the solution by @xmath222 , we obtain @xmath223 and , hence , @xmath224\right|_{n = n_0}=-\frac { 2}{s|\sin(\varphi/2)|}.\ ] ] we see that the term corresponding to @xmath225 in the right - hand side of ( a.4 ) is of order @xmath226 , whereas all other terms are of order @xmath227 . as a result , we obtain expression ( 56 ) .
amplitude @xmath155 ( 55 ) is presented in the form @xmath228}\times \nonumber \\
\times\exp{\left[2\sqrt{(n-\mu)^2-s^2}-2|n-\mu|{\rm arccosh}(|n-\mu|/s)\right]}\times\nonumber \\
\times\left\{\frac{\cot(\rho\pi)+\frac 12\frac{s^2}{(n-\mu)^2-s^2}-\sqrt{(n-\mu)^2-s^2 } } { \cot(\rho\pi)+\frac 12\frac{s^2}{(n-\mu)^2-s^2}+\sqrt{(n-\mu)^2-s^2}}+\right.\nonumber \\\left.+ \frac{\rm i}2\exp\left[2\sqrt{(n-\mu)^2-s^2}-2(n-\mu){\rm arccosh}\left(\frac{|n-\mu|}{s}\right)\right]\right\}^{-1},\end{aligned}\ ] ] where ( 41 ) is taken into account .
the sum in ( a.13 ) at @xmath6 is evaluated asymptotically at @xmath229 by converting it into integral : @xmath230 + o(1/\sqrt{s})\right\},\ ] ] where @xmath231 \frac{\cot(\rho\pi)+\frac 12\frac{1}{v^2 - 1}-s\sqrt{v^2 - 1}}{\cot(\rho\pi ) + \frac 12\frac{1}{v^2 - 1}+s\sqrt{v^2 - 1}}\times\nonumber \\ \times\left\{1 + 4\exp\left[4s\left(v{\rm arccosh}v-\sqrt{v^2 - 1}\right)\right ] \left[\frac{\cot(\rho\pi)+\frac 12\frac{1}{v^2 - 1}-s\sqrt{v^2 - 1}}{\cot(\rho\pi)+\frac 12 \frac{1}{v^2 - 1}+s\sqrt{v^2 - 1}}\right]^2\right\}^{-1}\end{aligned}\ ] ] and @xmath232 \left[\frac{\cot(\rho\pi)+\frac
12\frac{1}{v^2 - 1}-s\sqrt{v^2 - 1}}{\cot(\rho\pi)+\frac 12 \frac{1}{v^2 - 1}+s\sqrt{v^2 - 1}}\right]^2\right\}^{-1 } .\end{aligned}\ ] ] in the cases of @xmath233 and @xmath234 , the integrand in ( a.15 ) is maximal at @xmath235 , and one can use the laplace method for the estimation of @xmath236 in this case , see , e.g. , @xcite .
namely , let @xmath237 be continuous function with @xmath238 and @xmath239 be continuous differentiable function obeying conditions @xmath240 at @xmath241 and @xmath242 at @xmath243 @xmath244 , then @xmath245 where @xmath246 is the euler gamma - function . in the present case @xmath247 , @xmath248 , @xmath249 , @xmath250 ,
and we obtain @xmath251 a similar estimate is obtained for @xmath252 ( a.16 ) in the case of @xmath233 : @xmath253 in the case of @xmath254 , the integrand in ( a.16 ) is maximal at @xmath255 where @xmath256 .
the integrand in the vicinity of the maximum is approximated as @xmath257(v - v_0)^2\right\},\ ] ] and the integral is estimated as @xmath258 the integrand in ( a.15 ) in the case of @xmath254 changes sign at @xmath255 .
similarly to ( a.20 ) , we obtain the estimate for integral @xmath236 in this case : @xmath259 in the case of @xmath234 , the numerical analysis of @xmath252 yields an estimate of order @xmath260 .
hence , we obtain the estimate for amplitude @xmath261 ( a.13 ) : @xmath262 where constant @xmath263 is independent of @xmath264 , @xmath265 and @xmath5 .
in classical theory , we consider scattering of a point charged particle by an impenetrable tube containing magnetic flux @xmath7 .
the hamilton function corresponding to the motion outside of the tube is @xmath266,\ ] ] where @xmath267 , and the hamilton - jacobi equation ( see , e.g. , @xcite ) is @xmath268=0.\ ] ] a general solution to ( b.2 ) is @xmath269 primary conservation laws are @xmath270=e\ ] ] with @xmath271 , @xmath272 and @xmath273 being arbitrary constants ( @xmath274 ) .
taking the conservation laws into account , we obtain action @xmath275 ( b.3 ) in the form @xmath276 secondary conservation laws are @xmath277 with @xmath278 , @xmath279 and @xmath280 being arbitrary constants . thus we obtain relations @xmath281 @xmath282 @xmath283 which determine a trajectory of the particle , i.e. the functional dependence of @xmath284 , @xmath8 and @xmath285 on @xmath95 .
one can conclude that the classical trajectory is independent of the enclosed magnetic flux , since @xmath286 can be absorbed into constant @xmath271 : @xmath287 .
the trajectory is symmetric with respect to the point of reflection , @xmath288 .    as a consequence of ( b.7 ) ,
we obtain the following expression for the incidence angle which is equal to the reflection angle : @xmath289 where @xmath290 @xmath291 . defining the impact parameter as @xmath292 and using elementary tools of plane geometry
, we determine the angle of deflection ( called also the scattering angle ) in range @xmath145 as @xmath293 hence , the scattering differential cross section is @xmath294


there is large observational evidence of the presence of supermassive black holes ( smbh ) in the centres of most nearby galaxies  @xcite . the dynamical evidence related to the orbital motion of stars in the cluster surrounding sgr
@xmath7 indicates the presence of a smbh with mass @xmath8  @xcite .
in addition , the observed correlation between the central black hole masses and the stellar velocity dispersion of the bulge of the host galaxies suggests a direct connection between the formation and evolution of galaxies and smbh  @xcite .
the observation of luminous quasars detected at redshifts higher than 6 in the sloan digital sky survey ( sdss ) implies that smbh with masses @xmath9 , which are believed to be the engines of such powerful quasars , were formed within the first billion years after the big bang ( e.g. fan 2006 for a recent review ) .
however , it is still an open question how smbh seeds form and grow to reach such high masses in such a short amount of time  @xcite .
a number of different routes based on stellar dynamical processes , hydrodynamical processes or a combination of both have been suggested ( e.g. volonteri 2010 for a recent review ) .
one of the theoretical scenarios for smbh seed formation is the gravitational collapse of the first generation of stars ( population iii stars ) with masses @xmath10 that are expected to form in halos with virial temperature @xmath11 at @xmath12 where cooling by molecular hydrogen is effective . as a result of the gravitational collapse of such pop iii stars ,
very massive bhs would form and then grow via merger and accretion  @xcite .
another possible scenario proposes that if sufficient primordial gas in massive halos , with mass @xmath13 , is unable to cool below @xmath14 , it may lead to the formation of a supermassive object @xcite , which would eventually collapse to form a smbh .
this route assumes that fragmentation , which depends on efficient cooling , is suppressed , possibly by the presence of sufficiently strong uv radiation , that prevents the formation of molecular hydrogen in an environment with metallicity smaller than a given critical value  @xcite .
furthermore , fragmentation may depend on the turbulence present within the inflow of gas , and on the mechanism redistributing its angular momentum  @xcite .
the `` bars - within - bars '' mechanism  @xcite is a self - regulating route to redistribute angular momentum and sustain turbulence such that the inflow of gas can proceed without fragmenting as it collapses even in a metal - enriched environment .
depending on the rate and efficiency of the inflowing mass , there may be different outcomes .
a low rate of mass accumulation would favor the formation of isentropic supermassive stars ( sms ) , with mass @xmath15 , which then would evolve as equilibrium configurations dominated by radiation pressure  @xcite .
a different outcome could result if the accumulation of gas is fast enough so that the outer layers of sms are not thermally relaxed during much of their lifetime , thus having an entropy stratification  @xcite .
a more exotic mechanism that could eventually lead to a sms collapsing into a smbh is the formation and evolution of supermassive dark matter stars ( sdms ) ( spolyar et al 2008 ) .
such stars would be composed primarily of hydrogen and helium with only about @xmath16 of their mass in the form of dark matter , however they would shine due to dark matter annihilation .
it has recently been pointed out that sdmss could reach masses @xmath17  @xcite .
once sdmss run out of their dark matter supply , they experience a contraction phase that increases their baryon density and temperature , leading to an environment where nuclear burning may become important for the subsequent stellar evolution .
if isentropic sms form , their quasi - stationary evolution of cooling and contraction will drive the stars to the onset of a general relativistic gravitational instability leading to their gravitational collapse  @xcite , and possibly also to the formation of a smbh .
the first numerical simulations , within the post - newtonian approximation , of @xcite concluded that for spherical stars with masses greater than @xmath18 thermonuclear reactions have no major effect on the collapse , while less massive stars exploded due to hydrogen burning .
later @xcite performed the first relativistic simulations of the collapse of a sms in spherical symmetry .
they were able to follow the evolution until the formation of a bh , although their investigations did not include any microphysics . @xcite revisited the work of @xcite and performed simulations of non - rotating sms in the range of @xmath19-@xmath18 with post - newtonian corrections and detailed microphysics that took into account an equation of state ( eos ) including electron - positron pairs , and a reaction network describing hydrogen burning by the cno cycle and its break - out via the rapid proton capture ( rp)-process .
they found that sms with zero initial metallicity do not explode , while sms with masses larger than @xmath20 and with metallicity @xmath21 do explode .
more recently @xcite carried out general relativistic hydrodynamic simulations of the spherically symmetric gravitational collapse of sms adopting a spacetime foliation with outgoing null hypersurfaces to solve the system of einstein and fluid equations .
they performed simulations of spherical sms with masses in the range of @xmath22 using an eos that accounts for contributions from baryonic gases , and in a tabulated form , radiation and electron - positron pairs .
they were able to follow the collapse from the onset of the instability until the point of bh formation , and showed that an apparent horizon ( ah ) enclosing about 25% of the stellar material was formed in all cases when simulations stopped .
@xcite carried out general relativistic numerical simulations in axisymmetry of the collapse of uniformly rotating sms to bhs .
they did not take into account thermonuclear burning , and adopted a @xmath23-law eos , @xmath24 with adiabatic index @xmath25 , where @xmath26 is the pressure , @xmath27 the rest - mass density , and @xmath28 the specific internal energy .
although their simulations stopped before the final equilibrium was reached , the bh growth was followed until about 60@xmath29 of the mass had been swallowed by the smbh .
they estimated that about 90@xmath29 of the total mass would end up in the final smbh with a spin parameter of @xmath30 .
the gravitational collapse of differentially rotating sms in three dimensions was investigated by @xcite , who focused on the post - bh evolution , and also on the gravitational wave ( gw ) signal resulting from the newly formed smbh and the surrounding disk .
the gw signal is expected to be emitted in the low frequency lisa band ( @xmath31 hz ) .    despite the progress made ,
the final fate rotating isentropic sms is still unclear . in particular
, it is still an open question for which initial metallicities hydrogen burning by the @xmath32-limited hot cno cycle and its break - out via the @xmath33o@xmath34ne reaction ( rp - process ) can halt the gravitational collapse of rotating sms and generate enough thermal energy to lead to an explosion . to address this issue
, we perform a series of general relativistic hydrodynamic simulations with a microphysical eos accounting for contributions from radiation , electron - positron pairs , and baryonic matter , and taking into account the net thermonuclear energy released by the nuclear reactions involved in hydrogen burning through the pp - chain , cold and hot cno cycles and their break - out by the rp - process , and helium burning through the 3-@xmath35 reaction .
the numerical simulations were carried out with the nada code  @xcite , which solves the einstein equations coupled to the general relativistic hydrodynamics equations .
greek indices run from 0 to 3 , latin indices from 1 to 3 , and we adopt the standard convention for the summation over repeated indices . unless otherwise stated we use units in which @xmath36 .
next we briefly describe how the system of einstein and hydrodynamic equations are implemented in the nada code .
we refer to  @xcite for a more detailed description of the main equations and thorough testing of the code ( namely single bh evolutions , shock tubes , evolutions of both spherical and rotating relativistic stars , gravitational collapse to a bh of a marginally stable spherical star , and simulations of a system formed by a bh surrounded by a self - gravitating torus in equilibrium ) .
we follow the 3 + 1 formulation in which the spacetime is foliated into a set of non - intersecting spacelike hypersurfaces . in this approach
, the line element is written in the following form @xmath37 where @xmath35 , @xmath38 and @xmath39 are the lapse function , the shift three - vector , and the three - metric , respectively .
the latter is defined by @xmath40 where @xmath41 is a timelike unit - normal vector orthogonal to a spacelike hypersurface .
we make use of the bssn formulation  @xcite to solve the einstein equations .
initially , a conformal factor @xmath42 is introduced , and the conformally related metric is written as @xmath43 such that the determinant of the conformal metric , @xmath44 , is unity and @xmath45 , where @xmath46 .
we also define the conformally related traceless part of the extrinsic curvature @xmath47 , @xmath48 where @xmath49 is the trace of the extrinsic curvature .
we evolve the conformal factor defined as @xmath50 @xcite , and the auxiliary variables @xmath51 , known as the _ conformal connection functions _ , defined as @xmath52 where @xmath53 are the connection coefficients associated with @xmath44 .    during the evolution
we also enforce the constraints @xmath54 and @xmath55 at every time step .
we use the cartoon method @xcite to impose axisymmetry while using cartesian coordinates .
in addition to the bssn spacetime variables ( @xmath56 ) , there are two more quantities left undetermined , the lapse , @xmath35 and the shift vector @xmath57 .
we used the so - called `` non - advective 1+log '' slicing @xcite , by dropping the advective term in the `` 1+log '' slicing condition . in this case
, the slicing condition takes the form @xmath58 for the shift vector , we choose the `` gamma - freezing condition '' @xcite written as @xmath59 @xmath60 where @xmath61 is a constant that acts as a damping term , originally introduced both to prevent long term drift of the metric functions and to prevent oscillations of the shift vector .      the general relativistic hydrodynamics equations , expressed through the conservation equations for the stress - energy tensor @xmath62 and the continuity equation are    @xmath63    where @xmath27 is the rest - mass density
, @xmath64 is the fluid four - velocity and @xmath65 is the covariant derivative with respect to the spacetime metric .
following @xcite , the general relativistic hydrodynamics equations are written in a conservative form in cylindrical coordinates . since the einstein equations are solved only in the @xmath66 plane with cartesian coordinates ( 2d ) , the hydrodynamic equations are rewritten in cartesian coordinates for @xmath66 .
the following definitions for the hydrodynamical variables are used @xmath67 @xmath68 @xmath69 @xmath70 @xmath71 where @xmath72 and @xmath73 are the lorentz factor and the specific fluid enthalpy respectively , and @xmath26 is the pressure .
the conserved variables are @xmath74 , @xmath75 , @xmath76 .
we refer to @xcite for further details .
isentropic sms are self - gravitating equilibrium configurations of masses in the range of @xmath77 , which are mainly supported by radiation pressure , while the pressure of electron - positron pairs and of the baryon gas are only minor contributions to the eos .
such configurations are well described by newtonian polytropes with polytropic index @xmath78 ( adiabatic index @xmath25 ) .
the ratio of gas pressure to the total pressure ( @xmath32 ) for spherical sms can be written as  @xcite    @xmath79    where @xmath80 is the mean molecular weight .
thus @xmath81 for @xmath82 .
since nuclear burning timescales are too long for @xmath83 , evolution of sms proceeds on the kelvin - helmholtz timescale and is driven by the loss of energy and entropy by radiation as well as loss of angular momentum via mass shedding in the case of rotating configurations .
although corrections due to the nonrelativistic gas of baryons and electrons and general relativistic effects are small , they can not be neglected for the evolution .
firstly , gas corrections raise the adiabatic index slightly above @xmath84    @xmath85    secondly , general relativistic corrections lead to the existence of a maximum for the equilibrium mass as a function of the central density . for spherical sms
this means that for a given mass the star evolves to a critical density beyond which it is dynamically unstable against radial perturbations @xcite :    @xmath86    the onset of the instability also corresponds to a critical value of the adiabatic index @xmath87 , i.e. configurations become unstable when the adiabatic index drops below the critical value    @xmath88    this happens when the stabilizing gas contribution to the eos does not raise the adiabatic index above @xmath84 to compensate for the destabilizing effect of general relativity expressed by the second term on the righ - hand - side of eq .
( [ gcri ] ) .
rotation can stabilize configurations against the radial instability .
the stability of rotating sms with uniform rotation was analyzed by @xcite .
they found that stars at the onset of the instability have an equatorial radius @xmath89 , a spin parameter @xmath90 , and a ratio of rotational kinetic energy to the gravitational binding energy of @xmath91 .      to close the system of hydrodynamic equations ( eq .
[ hydro eqs ] ) we need to define the eos .
we follow a treatment which includes separately the baryon contribution on the one hand , and photons and electron - positron pairs contributions , in a tabulated form , on the other hand .
the baryon contribution is given by the analytic expressions for the pressure and specific internal energy    @xmath92    @xmath93    where @xmath94 the universal gas constant , @xmath95 the temperature , @xmath96 the baryon specific internal energy , and @xmath97 is the mean molecular weight due to ions , which can be expressed as a function of the mass fractions of hydrogen ( @xmath98 ) , helium ( @xmath99 ) and heavier elements ( metals ) ( @xmath100 ) as @xmath101 where @xmath102 is the average atomic mass of the heavy elements .
we assume that the composition of sms ( approximately that of primordial gas ) has a mass fraction of hydrogen @xmath103 and helium @xmath104 , where the metallity @xmath105 is an initial parameter , typically of the order of @xmath106 ( see table 1 details ) .
thus , for the initial compositions that we consider the mean molecular weight of baryons is @xmath107 ( i.e. corresponding to a molecular weight for both ions and electrons of @xmath108 ) .
effects associated with photons and the creation of electron - positron pairs are taken into account employing a tabulated eos . at temperatures above @xmath109 ,
not all the energy is used to increase the temperature and pressure , but part of the photon energy is used to create the rest - mass of the electron - positron pairs . as a result of pair creation , the adiabatic index of the star decreases , which means that the stability of the star is reduced .
given the specific internal energy , @xmath28 and rest - mass density , @xmath27 , as evolved by the hydrodynamic equations , it is possible to compute the temperature @xmath95 by a newton - raphson algorithm that solves the equation @xmath110 for @xmath95 ,    @xmath111    where @xmath112 is the iteration counter .      in order to avoid the small time steps , and cpu - time demands connected with the solution of a nuclear reaction network coupled to the hydrodynamic evolution , we apply an approximate method to take into account the basic effects of nuclear burning on the dynamics of the collapsing sms .
we compute the nuclear energy release rates by hydrogen burning ( through the pp - chain , cold and hot cno cycles , and their break - out by the rp - process ) and helium burning ( through the 3-@xmath35 reaction ) as a function of rest - mass density , temperature and mass fractions of hydrogen @xmath98 , helium @xmath99 and cno metallicity @xmath100 .
these nuclear energy generation rates are added as a source term on the right - hand - side of the evolution equation for the conserved quantity @xmath113 .
the change rates of the energy density due to nuclear reactions , in the fluid frame , expressed in units of [ @xmath114 @xmath115 @xmath116 are given by :    pp - chain @xcite :    @xmath117    where @xmath118 , and @xmath119 is given by @xmath120    3-@xmath35 @xcite :    @xmath121    where @xmath122 .
cold - cno cycle @xcite :    @xmath123    hot - cno cycle @xcite :    @xmath124    rp - process @xcite :    @xmath125    since we follow a single fluid approach , in which we solve only the hydrodynamics equations  eq .
( [ hydro eqs ] ) ( i.e. we do not solve additional advection equations for the abundances of hydrogen , helium and metals ) , the elemental abundances during the time evolution are fixed
. nevertheless , this assumption most possibly does not affect significantly the estimate of the threshold metallicity needed to produce a thermal bounce in collapsing sms .
the average energy release through the 3-@xmath35 reaction is about 7.275 mev for each @xmath126c nucleus formed . since the total energy due to helium burning for exploding models is @xmath127 ergs ( e.g. @xmath128 ergs for model s1.c ) ; and even considering that this energy is released mostly in a central region of the sms containing @xmath129 of its rest - mass  @xcite , it is easy to show that the change in the metallicity is of the order of @xmath130 .
therefore , the increase of the metallicity in models experiencing a thermal bounce is much smaller than the critical metallicities needed to trigger the explosions .
similarly , the average change in the mass fraction of hydrogen due to the cold and hot cno cycles is expected to be @xmath131 for exploding models .    [ cols="<,^,^,^,^,^,^,^,^,^,^,^ " , ]      after each time iteration the conserved variables ( i.e.  @xmath132 ) are updated and the _ primitive _ hydrodynamical variables ( i.e.  @xmath133 ) have to be recovered .
the recovery is done in such a way that it allows for the use of a general eos of the form @xmath134 .
we calculate a function @xmath135 , where @xmath136 and @xmath137 depend only on the conserved quantities and the pressure guess @xmath138 .
the new pressure is computed then iteratively by a newton - raphson method until the desired convergence is achieved .
the eos allows us to compute the neutrino losses due to the following processes , which become most relevant just before bh formation :    pair annihilation ( @xmath139    ): most important process above @xmath109 . due to the large mean free path of neutrinos in the stellar medium at the densities of sms
the energy loss by neutrinos can be significant . for a @xmath140 sms most of the energy release in the form of neutrinos originates from this process .
the rates are computed using the fitting formula given by @xcite .    photo - neutrino emission ( @xmath141    ): dominates at low temperatures @xmath142 and densities @xmath143 @xcite .
plasmon decay ( @xmath144    ): this is the least relevant process for the conditions encountered by the models we have considered because its importance increases at higher densities than those present in sms .
the rates are computed using the fitting formula given by @xcite .
the evolution equations are integrated by the method of lines , for which we use an optimal strongly stability - preserving ( ssp ) runge - kutta algorithm of fourth - order with 5 stages  @xcite .
we use a second - order slope limiter reconstruction scheme ( mc limiter ) to obtain the left and right states of the primitive variables at each cell interface , and a hlle approximate riemann solver  @xcite to compute the numerical fluxes in the @xmath145 and @xmath146 directions .
derivative terms in the spacetime evolution equations are represented by a fourth - order centered finite - difference approximation on a uniform cartesian grid except for the advection terms ( terms formally like @xmath147 ) , for which an upwind scheme is used .
the computational domain is defined as @xmath148 and @xmath149 , where @xmath150 refers to the location of the outer boundaries .
we used a cell - centered cartesian grid to avoid that the location of the bh singularity coincides with a grid point .
since it is not possible to follow the gravitational collapse of a sms from the early stages to the phase of black hole formation with a uniform cartesian grid ( the necessary fine zoning would be computationally too demanding ) , we adopt a regridding procedure @xcite . during the initial phase of the collapse we rezone the computational domain by moving the outer boundary inward , decreasing the grid spacing while keeping the initial number of grid points fixed .
initially we use @xmath151 grid points , and place the outer boundary at @xmath152 where @xmath153 is the equatorial radius of the star .
rezoning onto the new grid is done using a polynomial interpolation .
we repeat this procedure 3 - 4 times until the collapse timescale in the central region is much shorter than in the outer parts . at this point , we both decrease the grid spacing and also increase the number of grid points @xmath154 in dependence of the lapse function typically as follows : @xmath155 if @xmath156 , @xmath157 if @xmath158 , and @xmath159 if @xmath160 .
this procedure ensures the error in the conservation of the total rest - mass to be less than @xmath161 on the finest computational domain .      to deal with the spacetime singularity from the newly formed bh we use the method of excising the matter content in a region within the horizon as proposed by @xcite once
an ah is found .
this excision is done only for the hydrodynamical variables , and the coordinate radius of the excised region is allowed to increase in time . on the other hand ,
we do neither use excision nor artificial dissipation terms for the spacetime evolution , and solely rely on the gauge conditions .
here we define some of the quantities listed in table 1 .
we compute the total rest - mass @xmath162 and the adm mass @xmath163 as @xmath164 @xmath165,\end{aligned}\ ] ] where @xmath166 ( @xmath41 being the unit normal to the hypersurface ) and @xmath167 is the scalar curvature associated to the conformal metric @xmath44 .
the rotational kinetic energy @xmath168 and the gravitational potential energy @xmath72 are given by @xmath169 where @xmath170 is the angular velocity .
@xmath171 where the internal energy is computed as    @xmath172    in axisymmetry the ah equation becomes a nonlinear ordinary differential equation for the ah shape function , @xmath173 @xcite .
we employ an ah finder that solves this ode by a shooting method using @xmath174 and @xmath175 as boundary conditions .
we define the mass of the ah as @xmath176 where @xmath177 is the area of the ah .     with zero metallicity ) for three different eos ( @xmath23-law and the microphysical eos with and without the electron - positron pair creation).,width=283 ]     and zero metallicity ) with a @xmath23-law and the microphysical eos with electron - positron pair creation .
the middle and lower panels display the ah mass and the disk mass as a function of time for the collapse simulation with a @xmath23-law eos.,width=283 ]
the initial sms are set up as isentropic objects .
all models , except model d2 , are chosen such that they are gravitationally unstable , and therefore their central rest - mass density is slightly larger than the critical central density required for the onset of the collapse of a configuration with given mass and entropy . a list of the different sms we have considered is provided in table 1 .
models s1 and s2 represent a spherically symmetric , nonrotating sms with gravitational mass of @xmath178 and @xmath179 , respectively , while models r1 and r2 are uniformly rotating initial models again with masses of @xmath178 and @xmath179 , respectively . the rigidly and maximally rotating initial models r1 and r2 , and the differentially rotating models d1 and d2 are cconstructed with a polytropic eos with the lorene code ( url http://www.lorene.obspm.fr ) .
we obtain temperatures for our microphyscial models by inverting the corresponding energy density with our eos of section 3.2 .
we also introduce a perturbation to trigger the gravitational collapse by reducing the pressure overall by @xmath180 .    in order to determine the threshold metallicity required to halt the collapse and produce an explosion we carry out several numerical simulations for each initial model with different values of the initial metallicity .
the initial metallicities along with the fate of the star are given in table 1 .
axisymmetric calculations without rotation ( i.e. models s1 and s2 ) retain the spherical symmetry of the initial conditions .
there are no physical phenomena like convective or overturn instabilities to produce asphericity , i.e. we can directly compare our 2d non - rotating models with those computed in spherical symmetry ( 1d calculations ) by  @xcite and @xcite .
the main differences with respect to the results obtained by  @xcite are most likely due to two reasons .
first , we apply a fully general relativistic treatment while they used a post - newtonian treatment of gravity , and second , there are differences in the treatment of nuclear burning ( fuller et al .
1986 solved the relevant nuclear network without the approximations adopted in our work , see section 3.3 for further details ) . despite these differences ,
the results agree fairly well . as discussed in detail in section 7.1
the initial metallicities required to produce an explosion are similar , and a thermal bounce can be produced only if sufficient energy is liberated during the phase when the hcno cycle is active .
the main difference with respect to the work of  @xcite resides in the formulation of einstein s field equations ; in particular , in the foliation of the spacetime ( foliation into a set spacelike hypersurfaces versus a foliation with outgoing null hypersurfaces ) . in order to compare with the results of @xcite , we computed the redshifted total energy output for a model having the same rest - mass , doing the time integration until approximately the same evolutionary stage as in @xcite .
we find that the total energies released in neutrinos differ by less than @xmath181 ( for mode details see section 7.3 ) .
cm      previous simulations of sms collapse to bh in general relativity have been performed with a @xmath23-law eos with @xmath182 ( with the only exception being the work of linke et al .
2001 ) . in order to elucidate the influence of the eos on the dynamics of collapsing sms
, we performed three simulations of the same initial model ( model r1.0 , a marginally unstable uniformly rotating sms with zero initial metallicity ) without nuclear burning effects , and with three eos : a @xmath23-law eos with @xmath182 ( i.e. a similar set - up as in shibata & shapiro 2002 ) and the microphysical eos with and without including electrons and the e@xmath183 pairs , i.e. for the last eos case we consider eqs .
( 19 ) and ( 20 ) for the baryons plus @xmath184 and @xmath185 for the photons ( where @xmath186 is the radiation density constant ) .
we denote the @xmath23-eos as eos-0 , the full microphysical eos , our canonical one for the studies of this work as eos-1 , and the reduced microphysical case as eos-2 .    in figure  [ fig1 ]
we show with a dotted line the time evolution of the central density of model r1.0 with eos-0 , and with a dashed ( solid ) line the time evolution of the central density with eos-2 ( eos-1 ) .
the first thing to note is that the collapse timescale obtained with the @xmath23-law eos is shorter than that obtained with the microphysical eos ( both eos-1 and eos-2 ) , because the ion pressure contribution to the eos raises the adiabatic index above @xmath84 ( see eq .
[ ion_gamma ] ) .
this increase in the adiabatic index helps to stabilize the star against the gravitational instability , and therefore delays the collapse .    on the other hand
, the effect of pair creation reduces the adiabatic index below @xmath84 at @xmath187k .
this explains the differences between the solid and dashed lines in figure  [ fig1 ] at central densities @xmath188 g/@xmath189 , which correspond to central temperatures @xmath190k .
once the collapse enters this regime , pair creation becomes relevant enough to reduce the adiabatic index below @xmath84 , which destabilizes the collapsing star .
compared to previous works ( shibata & shapiro 2002 ) the use of a microphysical eos instead of a @xmath23-law eos delays the collapse ( mostly due to the baryons while e@xmath183 destabilize ) of an initially gravitationally unstable configuration .
we also performed 2d axisymmetric simulations of differentially rotating sms .
first , we investigated the influence of the eos on _ gravitationally unstable stars _ using model d1 as a reference .
this model corresponds , within the accuracy to which the initial conditions can be reproduced , to a differentially rotating unstable sms discussed by  @xcite ( i.e. their model i ) .
results are displayed in figure  [ fig2 ] . in the upper panel of this figure
, we show the time evolution of the central density for model d1 with eos-0 ( dashed line ) , and with the microphysical eos-1 ( solid line ) .
opposite to the behavior in the case of the uniformly rotating sms r1.0 , the collapse timescale is longer with the @xmath23-law than with the microphysical eos .
the reason for this difference is that the initial central temperature ( @xmath191k ) of model d1 is an order of magnitude higher than the initial central temperature in r1.0 .
therefore , electron - positron pair creation reduces the stability of the star ( by reducing @xmath23 ) already during the initial stages of the collapse .
this behavior is expected to be present also in 3d , and since the collapse timescale is reduced when using the microphysical eos , nonaxisymmetric instablilities would have even less time to grow before the formation of a bh .
it reinforces the conclusions of  @xcite , who showed that the three dimensional collapse of rotating stars proceeds in an approximately axisymmetric manner .
the lower two panels of figure  [ fig2 ] display the growth of the ah mass and the disk mass ( defined as the rest - mass outside the ah of the newly formed bh ) as a function of time ( in units of the gravitational mass for comparison with figures 9 and 10 of saijo & hawke 2009 ) , respectively .
the values of both quantities at the end of the simulation agree , within a @xmath192 difference , with those obtained in 3d by  @xcite .
we note that there is also good agreement ( less than @xmath192 difference ) regarding the time at which an ah is first detected .
these observed small differences are likely due to differences in the initial models and numerical techniques rather than to the influence of nonaxisymmtric effects .
this suggests that our collapse simulation with the same treatment of physics yields good agreement with the 3d simulations of  @xcite    we also investigated the influence of electron - positron pair creation on the evolution of _ gravitationally stable differentially rotating _ sms using model d2 which is similar to the stable differentially rotating model iii of @xcite .
we performed three simulations of model d2 varying the eos . in figure  [ fig3 ]
we show the central rest - mass density as a function of time for a @xmath23-law ( dashed line ) , and for the microphysical eos-1 ( solid line ) and eos-2 ( dotted line ) . in agreement with the results obtained by @xcite we find that model d2 represents a stable differentially rotating sms when a @xmath23-law is used .
a persistent series of oscillations is triggered by the initial perturbation in the pressure .
this is also the case with the microphysical eos-2 without the inclusion of electrons and the e@xmath193 pairs .
however , the time evolution of d2 is completely different when e@xmath193 pairs are taken into account .
the influence of pairs is large enough to destabilize model d2 against gravitational collapse .
we note that unlike all other sms considered in this paper , which are @xmath25 models initially unstable to gravitational collapse , model d2 is an initially stable @xmath25 model which becomes gravitationally unstable only by the creation of electron - positron pairs at high temperatures .
hence , using a microphysical eos with electron - positron pairs is crucial to determine the stability of differentially rotating sms .
we note that the central temperature of the initial models d1 and d2 is of the order of @xmath194k . at this temperature , the main source of thermonuclear energy is hydrogen burning via the rp - process .
it is however expected that such sms would previously experience a phase of hydrogen burning via the cold and hot cno cycles which would significantly affect the evolution of the models such that configurations with high @xmath195 as in models d1 and d2 might never be reached .
therefore , models d1 and d2 are not particularly well suited to investigate the existence of a thermal bounce during collapse ( see section 7 ) . exploring in detail the parameter space for the stability of differentially rotating sms with the microphysical eos , and
the existence of a thermal bounce during the collapse phase depending on the initial stellar metallicity , is a major task on its own , which is beyond the scope of this paper .
for these reasons we do not consider differentially rotating sms in this work .
first we consider a gravitationally unstable spherically symmetric sms with a gravitational mass of @xmath178 ( s1.a , s1.b and s1.c ) , which corresponds to a model extensively discussed in  @xcite , and therefore allows for a comparison with the results presented here .
@xcite found that unstable spherical sms with @xmath178 and an initial metallicity @xmath196 collapse to a bh while models with an initial metallicity @xmath197 explode due to the nuclear energy released by the hot cno burning .
they also found that the central density and temperature at thermal bounce ( where the collapse is reversed to an explosion ) are @xmath198 g/@xmath189 and @xmath199 k , respectively .
the left panels in figure  [ fig4 ] show the time evolution of the central rest - mass density ( upper panel ) and central temperature ( lower panel ) for models s1.a , s1.c , r1.a and r1.d , i.e. , non - rotating and rotating models with a mass of @xmath200 .
in particular , the solid lines represent the time evolution of the central density and temperature for model s1.c ( @xmath201 ) and r1.d ( @xmath202 ) . as the collapse proceeds ,
the central density and temperature rise rapidly , which increases the nuclear energy generation rate by hydrogen burning .
since the metallicity is sufficiently high , enough energy can be liberated to increase the pressure and to produce a thermal bounce .
this is the case for model s1.c . in figure  [ fig4 ]
we show that a thermal bounce occurs ( at approximately @xmath203 s ) entirely due to the hot cno cycle , which is the main source of thermonuclear energy at temperatures in the range @xmath204 .
the rest - mass density at bounce is @xmath205 g/@xmath189 and the temperature @xmath206 .
these values , as well as the threshold metallicity needed to trigger a thermonuclear explosion ( @xmath201 ) , are higher than those found by  @xcite ( who found that a spherical nonrotating model with the same rest - mass would explode , if the initial metallicity was @xmath197 ) .
on the other hand , dashed lines show the time evolution of the central density and temperature for model s1.a ( @xmath197 ) . in this case , as well as for model s1.b , the collapse is not halted by the energy release and continues until an ah is found , indicating the formation of a bh .
we note that the radial velocity profiles change continuously near the time where the collapse is reversed to an explosion due to the nuclear energy released by the hot cno burning , and an expanding shock forms only near the surface of the star at a radius @xmath207 ( i.e. @xmath208 ) where the rest - mass density is @xmath209 .
we show in figure  [ fig5 ] the profiles of the @xmath145-component of the three - velocity @xmath210 along the @xmath145-axis ( in the equatorial plane ) for the nonrotating spherical stars s1.a ( dashed lines ) and s1.c ( solid lines ) at three different time slices near the time at which model s1.c experiences a thermal bounce .
velocity profiles of model s1.c are displayed up to the radius where a shock forms at @xmath211 and begins to expand into the low density outer layers of the sms .
the evolutionary tracks for the central density and temperature of the rotating models r1.a and r1.d are also shown in figure  [ fig4 ] .
a dashed line corresponds to model r1.a , with an initial metallicity @xmath202 , which collapses to a bh .
a solid line denotes model r1.d with @xmath196 , which explodes due to the energy released by the hot cno cycle .
we find that model r1.c with a lower metallicity of @xmath212 also explodes when the central temperature is the range dominated by the hot cno cycle .    as a result of the kinetic energy stored in the rotation of models r1.c and r1.d , the critical metallicity needed to trigger an explosion decreases significantly relative to the non - rotating case .
we observe that rotating models with initial metallicities up to @xmath213 do not explode even via the rp - process , which is dominant at temperatures above @xmath214 k and increases the hydrogen burning rate by @xmath215 times relative to the hot cno cycle .
we also note that the evolution time scales of the collapse and bounce phases are reduced because rotating models are more compact and have a higher initial central density and temperature than the spherical ones at the onset of the gravitational instability .
the right panels in figure  [ fig4 ] show the time evolution of the central rest - mass density ( upper panel ) and the central temperature ( lower panel ) for models s2.a , s2.b , r2.a and r2.d , i.e. , of models with a mass of @xmath216 . we find that the critical metallicity for an explosion in the spherical case is @xmath217 ( model s2.b ) , while model s2.a with @xmath218 collapses to a bh .
we note that the critical metallicity leading to a thermonuclear explosion is higher than the critical value found by  @xcite ( i.e. , @xmath219 ) for a spherical sms with the same mass .
the initial metallicity leading to an explosion in the rotating case ( model r2.d ) is more than an order of magnitude smaller than in the spherical case . as for the models with a smaller gravitational mass
, the thermal bounce takes place when the physical conditions in the central region of the star allow for the release of energy by hydrogen burning through the hot cno cycle .
overall , the dynamics of the more massive models indicates that the critical initial metallicity required to produce an explosion increases with the rest - mass of the star .
figure  [ fig6 ] shows the total nuclear energy generation rate in erg@xmath220s for the exploding models as a function of time during the late stages of the collapse just before and after bounce .
the main contribution to the nuclear energy generation is due to hydrogen burning by the hot cno cycle .
the peak values of the energy generation rate at bounce lie between several @xmath221 for the rotating models ( r1.d and r2.d ) , and @xmath222 for the spherical models ( s1.c and s2.b ) .
as expected the maximum nuclear energy generation rate needed to produce an explosion is lower in the rotating models .
moreover , as the explosions are due to the energy release by hydrogen burning via the hot cno cycle , the ejecta would mostly be composed of @xmath223he .
-component of the three - velocity @xmath210 along the @xmath145-axis ( in the equatorial plane ) for the nonrotating spherical stars s1.a ( dashed lines ) and s1.c ( solid lines ) at three different time slices near the time at which model s1.c experiences a thermal bounce .
velocity profiles of model s1.c are displayed up to the radius where a shock , that expands into the low density outer layers of the sms , forms.,width=302 ]    s for the exploding models ( s1.c , r1.d , r1.c , s2.b and r2.d ) as a function of time near the bounce .
the contribution to the nuclear energy generation is mainly due to hydrogen burning by the hot cno cycle .
the peak values of the energy generation rate at bounce lie between @xmath224 $ ] for the rotating models ( r1.d and r2.d ) , and @xmath225 $ ] for the spherical models ( s1.c and s2.b).,width=302 ]    as a result of the thermal bounce , the kinetic energy rises until most of the energy of the explosion is in the form of kinetic energy .
we list in the second but last column of table 1 the radial kinetic energy after thermal bounce , which ranges between @xmath226 for the rotating star r1.c , and @xmath227 for the spherical star s2.b .
due to the lack of resolution at the surface of the star , it becomes difficult to compute accurately the photosphere and its effective temperature from the criterion that the optical depth is @xmath228 .
therefore , in order to estimate the photon luminosity produced in association with the thermonuclear explosion , we make use of the fact that within the diffusion approximation the radiation flux is given by      where @xmath230 is the energy density of the radiation , and @xmath231 is the opacity due to electron thompson scattering , which is the main source of opacity in sms .
the photon luminosity in terms of the temperature gradient and for the spherically symmetric case can be written as @xmath232 where @xmath186 is the radiation constant , and @xmath233 the speed of light .
as can be seen in the last panel of figure  [ fig7 ] the distribution of matter becomes spherically symmetric during the phase of expansion after the thermonuclear explosion . in this figure ( fig .
[ fig7 ] ) we show the isodensity contours for the rotating model r1.d .
the frames have been taken at the initial time ( left figure ) , at @xmath234s ( central figure ) just after the thermal bounce ( at @xmath235s ) , and at @xmath236 s when the radius of the expanding matter is roughly 4 times the radius of the star at the onset of the collapse .
the photon luminosity computed using eq.([photon_lum ] ) for model r1.d is displayed in figure  [ fig8 ] , where we also indicate with a dashed vertical line the time at which the thermal bounce takes place .
the photon luminosity before the thermal bounce is computed at radii inside the star unaffected by the local dynamics of the low density outer layers which is caused by the initial pressure perturbation and by the interaction between the surface of the sms and the artificial atmosphere .
once the expanding shock forms near the surface , the photon luminosity is computed near the surface of the star .
the lightcurve shows that , during the initial phase , the luminosity is roughly equal to the eddington luminosity @xmath237 until the thermal bounce .
then , the photon luminosity becomes super - eddington when the expanding shock reaches the outer layers of the star and reaches a value of @xmath238 .
this value of the photon luminosity after the bounce is within a few percent difference with respect to the photon luminosity @xcite found for a nonrotating sms of same rest - mass .
the photon luminosity remains super - eddington during the phase of rapid expansion that follows the thermal bounce .
we compute the photon luminosity until the surface of the star reaches the outer boundary of the computational domain @xmath239 after the bounce . beyond that point
, the luminosity is expected to decrease , and then rise to a plateau of @xmath240 due to the recombination of hydrogen ( see fuller et al .
1986 for a nonrotating star ) .
the outcome of the evolution of models that do not generate enough nuclear energy during the contraction phase to halt the collapse is the formation of a bh .
the evolutionary tracks for the central density and temperature of some of these models are also shown in figure  [ fig4 ] .
the central density typically increases up to @xmath241 and the central temperature up to @xmath242 just before the formation of an ah .
three isodensity contours for the rotating model r1.a collapsing to a bh are shown in figure  [ fig9 ] , which display the flattening of the star as the collapse proceeds .
the frames have been taken at the initial time ( left panel ) , at @xmath243 ( central panel ) approximately when model r1.d with higher metallicity experiences a thermal bounce and , at @xmath244 , where a bh has already formed and its ah has a mass of @xmath245 of the total initial mass .    at the temperatures reached during the late stages of the gravitational collapse ( in fact at @xmath246 )
the most efficient process for hydrogen burning is the breakout from the hot cno cycle via the @xmath33o@xmath34ne reaction .
nevertheless , we find that models which do not release enough nuclear energy by the hot cno cycle to halt their collapse to a bh , are not able to produce a thermal explosion due to the energy liberated by the @xmath33o@xmath34ne reaction .
we note that above @xmath109 , not all the liberated energy is used to increase the temperature and pressure , but is partially used to create the rest - mass of the electron - positron pairs . as a result of pair creation
, the adiabatic index of the star decreases , which means the stability of the star is reduced .
moreover , due to the presence of e@xmath183 pairs , neutrino energy losses grow dramatically .
figure  [ fig10 ] shows ( solid lines ) the time evolution of the redshifted neutrino luminosities of four models collapsing to a bh ( s1.a , r1.a , s2.a , and r2.a ) , and ( dashed lines ) of four models experiencing a thermal bounce ( s1.c , r1.d , s2.b , and r2.d ) .
the change of the slope of the neutrino luminosities at @xmath247 denotes the transition from photo - neutrino emission to the pair annihilation dominated region .
the peak luminosities in all form of neutrino for models collapsing to a bh are @xmath248 .
neutrino luminosities can be that important because the densities in the core prior to bh formation are @xmath249 , and therefore neutrinos can escape .
the peak neutrino luminosities lie between the luminosities found by @xcite for the collapse of spherical sms , and those found by @xcite ( who only took into account the luminosity in the form of electron antineutrino ) .
the maximum luminosity decreases slightly as the rest - mass of the initial model increases , which was already observed by  @xcite .
in addition , we find that the peak of the redshifted neutrino luminosity does not seem to be very sensitive to the initial rotation rate of the star .
we also note that the luminosity of model r1.a reflects the effects of hydrogen burning at @xmath250 .
the total energy output in the form of neutrinos is listed in the last column of table 1 for several models .
the total radiated energies vary between @xmath5 ergs for models collapsing to a bh , and @xmath6 ergs for exploding models .
these results are in reasonable agreement with previous calculations .
for instance , @xcite obtained that the total energy output in the form of electron antineutrinos for a spherical sms with a mass @xmath251 and zero initial metallicity was @xmath252 , although their simulations neglected general relativistic effects which are important to compute accurately the relativistic redshifts . on the other hand , @xcite , by means of relativistic one - dimensional simulations , found a total radiated energy in form of neutrinos of about @xmath253 for the same initial model , and about @xmath254 when redshifts were taken into account . in order to compare with the results of @xcite
, we computed the redshifted total energy output for model s1.a , having the same rest - mass , until approximately the same evolution stage as @xcite did ( i.e. when the differential neutrino luminosity @xmath255 ) .
we find that the total energy released in neutrinos is @xmath256 .
the neutrino luminosities for models experiencing a thermonuclear explosion ( dashed lines in fig .  [ fig10 ] ) peak at much lower values @xmath258 , and decrease due to the expansion and disruption of the star after the bounce .      the axisymmetric gravitational collapse of rotating sms with uniform rotation is expected to emit a burst of gravitational waves  @xcite with a frequency within the lisa low frequency band ( @xmath259 ) .
although through the simulations presented here we could not investigate the development of nonaxisymmetric features in our axisymmetric models that could also lead to the emission of gws ,  @xcite have shown that the three dimensional collapse of rotating stars proceeds in an approximately axisymmetric manner .        where @xmath264 refers to the second time derivative of the quadrupole moment .
the gravitational wave quadrupole amplitude is @xmath265 .
following @xcite we compute the second time derivative of the quadrupole moment by finite differencing the numerical results for the first time derivative of @xmath266 obtained by          where d is the distance of the source , and @xmath269 the spectral energy density of the gravitational radiation given by @xmath270 with @xmath271 we have calculated the quadrupole gravitational wave emission for the rotating model r1.a collapsing to a bh .
we plot in figure  [ fig11 ] the characteristic gravitational wave strain ( eq.[h_char ] ) for this model assuming that the source is located at a distance of 50 gpc ( i.e. , @xmath272 ) , together with the design noise spectrum @xmath273 of the lisa detector @xcite .
we find that , in agreement with @xcite , @xcite and @xcite , the burst of gravitational waves due to the collapse of a rotating sms could be detected at a distance of 50 gpc and at a frequency which approximately takes the form @xcite          furthermore , @xcite have recently investigated , by means of three - dimensional general relativistic numerical simulations of equilibrium tori orbiting bhs , the development of the nonaxisymmetric papaloizou - pringle instability ( ppi ) in such systems  @xcite , and have found that a nonaxisymmetric instability associated with the @xmath277 mode grows for a wide range of self - gravitating tori orbiting bhs , leading to the emission of quasiperiodic gws .
in particular ,  @xcite have pointed out that the emission of quasiperiodic gws from the torus resulting after the formation of a smbh via the collapse of a sms could be well above the noise sensitivity curve of lisa for sources located at a distance of 10gpc .
such instability appears for tori whose angular velocity in the equatorial plane expressed as @xmath278 has @xmath279 where @xmath280 corresponds to the keplerian limit , i.e. @xmath281 in newtonian gravity .
we find that the torus ( defined as the rest - mass outside the ah ) that forms after the collapse to a bh of the uniformly rotating model r1.a ( when the mass of the ah exceeds @xmath245 of the adm mass ) does not fulfill the above condition for the development of the ppi .
however , we find that the torus that forms when the differentially rotating model d1 collapses to a bh has a distribution of angular momentum such that @xmath278 with @xmath282 .
this suggests that the torus may be prone to the development of the nonaxisymmetric ppi , which would lead to the emission of quasiperiodic gws with peak amplitude @xmath283 and frequency @xmath284 during an accretion timescale @xmath285 .
we have presented results of general relativistic simulations of collapsing supermassive stars using the two - dimensional general relativistic numerical code nada , which solves the einstein equations written in the bssn formalism and the general relativistic hydrodynamic equations with high resolution shock capturing schemes .
these numerical simulations have used an eos that includes the effects of gas pressure , and tabulated those associated with radiation pressure and electron - positron pairs .
we have also taken into account the effects of thermonuclear energy release by hydrogen and helium burning .
in particular , we have investigated the effects of hydrogen burning by the @xmath32-limited hot cno cycle and its breakout via the @xmath33o@xmath34ne reaction ( rp - process ) on the gravitational collapse of nonrotating and rotating sms with non - zero metallicity .
we have presented a comparison with previous studies , and investigated the influence of the eos on the collapse .
we emphasize that axisymmetric calculations without rotation ( i.e. models s1 and s2 ) retain the spherical symmetry of the initial configurations as there are no physical phenomena to produce asphericity and numerical artifacts associated with the use of cartesian coordinates are negligibly small .
overall , our collapse simulations yield good agreement with previous works when using the same treatment of physics .
we have also found that the collapse timescale depends on the ion contributions to the eos , and electron - positron pair creation affects the stability of sms .
interestingly , differentially rotating stars that are gravitationally stable with a @xmath25 eos can become unstable against gravitational collapse when the calculation is performed with the microphysical eos including pair creation .
we have found that objects with a mass of @xmath286 and an initial metallicity greater than @xmath1 explode if non - rotating , while the threshold metallicity for an explosion is reduced to @xmath2 for objects which are uniformly rotating .
the critical initial metallicity for a thermal explosion increases for stars with a mass of @xmath3 .
the most important contribution to the nuclear energy generation is due to the hot cno cycle .
the peak values of the nuclear energy generation rate at bounce range from @xmath287 for rotating models ( r1.d and r2.d ) , to @xmath288 for spherical models ( s1.c and s2.b ) . after the thermal bounce , the radial kinetic energy of the explosion rises until most of the energy is kinetic , with values ranging from @xmath289 for rotating stars , to up @xmath290 for the spherical star s2.b .
the neutrino luminosities for models experiencing a thermal bounce peak at @xmath291 .
the photon luminosity roughly equal to the eddington luminosity during the initial phase of contraction . then , after the thermal bounce , the photon luminosity becomes super - eddington with a value of about @xmath292 during the phase of rapid expansion that follows the thermal bounce .
for those stars that do not explode we have followed the evolution beyond the phase of black hole formation and computed the neutrino energy loss .
the peak neutrino luminosities are @xmath4 .
sms with masses less than @xmath3 could have formed in massive halos with @xmath293 . although the amount of metals that was present in such environments at the time when sms might have formed is unclear , it seems possible that the metallicities could have been smaller than the critical metallicities required to reverse the gravitational collapse of a sms into an explosion .
if so , the final fate of the gravitational collapse of rotating sms would be the formation of a smbh and a torus . in a follow - up paper
, we aim to investigate in detail the dynamics of such systems ( collapsing of sms to a bh - torus system ) in 3d , focusing on the post - bh evolution and the development of nonaxisymmetric features that could emit detectable gravitational radiation .
we thank b. mller and p. cerd - durn for useful discussions . work supported by the deutsche forschungsgesellschaft ( dfg ) through its transregional centers
sfb / tr 7 `` gravitational wave astronomy '' , and sfb / tr 27 `` neutrinos and beyond '' , and the cluster of excellence exc153 `` origin and structure of the universe '' .