the naked reissner - nordstrm ( r - n ) singularity is a classical general relativistic solution in electrovacuum . the solution is expected to have a very limited meaning , due to the fact that such singularities can not be created neither by a gravitational collapse , nor by dropping a charge into the black hole . ( according to the weak cosmic censorship conjecture general naked singularities should be prohibited in general theory of relativity , although there are indications that by including quantum effects the violations of the conjecture could be considered @xcite . ) moreover , a naked singularity created from some exotic initial data conditions should become quickly neutralized ( classically , or via quantum pair production ) .
some results also indicate that if one considers electro - gravitational perturbations the r - n naked singularity becomes linearly unstable @xcite .
however it was discovered that the scalar field scattering problem on such a singular background can be still well defined @xcite , since the waves remain regular at the origin .
( however , the back reaction of the given scalar field configuration might still excite some of the unstable electro - gravitational modes and this would eventually lead to a breakdown of the perturbation approach . ) despite the nice regularity property of the scattering problem , the spacetime is non - globally hyperbolic and the time evolution of the fields is not unique @xcite .
this means one has to specify an additional boundary condition at the singularity to obtain a fully unique time evolution .
another way of seeing the problem is through the language of operators : one can understand the spatial part of the wave operator as a positive symmetric operator acting on a @xmath4 hilbert space , and then obtain the scalar field dynamics through a suitable positive self - adjoint extension of such a symmetric operator @xcite .
( one `` preferred '' way in which such a self - adjoint extension can be always realized is through the so called friedrich s extension @xcite , which will also be the case of this paper . ) anyway , after uniquely specifying the dynamics , one should be able to characterize the scattering by a set of characteristic oscillations , the quasi - normal modes .
low damped quasi - normal modes are in general used as a possible source of information about potential astrophysical objects ( such as neutron stars , black holes ) , and the highly damped modes are potentially interesting from the point of view of quantum gravity @xcite . since a lot of work was devoted to the problem of quasinormal modes of the reissner - nordstrm black hole , it might be interesting to observe what happens if one transits from the r - n black hole case to the r - n naked singularity case ( with a reflective boundary condition ) .
information about `` what happens '' shows how many features of the quasinormal modes of the black hole spacetimes are specific to the black holes themselves and what features survive much more general conditions .
let us also give one concrete example why asking what happens with the quasinormal modes of the naked r - n singularity might be interesting : in the black hole case the behavior of the asymptotically highly damped modes is widely suspected to be linked to the properties of the black hole spacetime horizon(s ) , ( or more specifically it is considered to carry information about _ quantum _ black holes ) .
what does then happen with the asymptotically highly damped modes in case there are _ no _ horizons present ?
( if the behavior of the highly damped modes in case of r - n naked singularity would , for instance , resemble the behavior of the highly damped modes of the black hole , it would be a disturbing fact from the point of view of the popular conjectures linking the modes to the horizon s quantum area spacing .
moreover , investigating the r - n naked singularity from this point of view is especially attractive , as the r - n naked singularity is obtained by a continuous transition in the @xmath5 parameters from the r - n black hole spacetime . )
thus , briefly , we hope that despite the fact that most likely the r - n naked singularity model does not correspond to a realistic physical situation , there are still many interesting things one can learn from such a model .
the structure of this paper is as follows : in the second section we analyse the problem of the uniqueness of the time evolution of scalar fields on a r - n naked singularity background . in the sections
three and four we analyse the properties of the effective potential for the scalar fields scattering and the geometrical optics ( eikonal ) limit of such a scattering problem . in the fifth section
we define analytically solvable potentials that can give good approximations to the problem of the low damped qnm frequencies ( such that characterize the given scattering problem ) . in the sixth section we use those analytical approximations to derive semi - analytical results for the qnm frequencies in the eikonal limit . in the seventh section we use the numerical characteristic integration to obtain the low damped frequencies for the low values of @xmath3 . in the eighth section we analyse what happens with
the asymptotically highly damped modes and we suggest that in case of naked r - n singularity such modes do _ not _ exist ( as one might expect considering some presently popular conjectures @xcite ) .
we give the final conclusions in the section nine .
we provide also appendices with more detailed results and some further suggestions for the analytical approximations of the problem .
in this section we will follow the standard analysis of the scalar field evolution in a curved background .
( as an example of such an analysis see the treatment of schwarzschild black hole perturbations in @xcite . for a review that presents also such techniques
see for example @xcite . )
take the klein - gordon equation for the complex ( charged ) scalar field : @xmath6 with the metric line element given as @xmath7 for the reissner - nordstrm ( r - n ) singularity the function @xmath8 is in planck units given as : @xmath9 take the decomposition of the field into the spherical harmonics @xmath10 after we separate the variables we obtain the following reduced equation @xmath11   \nonumber \\    & -&\frac{l(l+1)f(r)}{r^{2}}\psi_{l}(t , r ) .
\end{aligned}\ ] ] we are interested only in compactly supported data initial value problem : first assume that @xmath12 is everywhere bounded and hence the laplace transform of @xmath12 exists : @xmath13 the laplace transformed equation gives the following equation : @xmath14 \nonumber \\
& -&\frac{l(l+1)f(r)}{r^{2}}\tilde\psi_{l}(s , r)+i_{l}(s , r),\end{aligned}\ ] ] where @xmath15_{|t=0}.\ ] ] the solution that corresponds to the initial data term @xmath16 is obtained by the inverse laplace transform of @xmath17 and the function @xmath17 is given as @xmath18 here @xmath19 is a green s function satisfying @xmath20   \\ & + & \left[s^{2}+\frac{l(l+1)f(r)}{r^{2}}\right ] g_{l}(s , r , r')=\delta(r - r ' ) .
\nonumber\end{aligned}\ ] ] since @xmath12 was bounded , its laplace transform @xmath17 must be also bounded in @xmath21 .
this translates to the boundedness of the green s function in @xmath21 .
the unique solution of the given initial data problem is obtained if the condition of boundedness of the green s function leads to a _
unique _ way to construct green s function from the two linearly independent solutions @xmath22 of the homogeneous equation @xmath23 \\ & + & \left[s^{2}+\frac{l(l+1)f(r)}{r^{2}}\right ] u_{l1,2}(s , r)=0.\nonumber\end{aligned}\ ] ] if @xmath8 goes to 1 at spatial infinity ( the metric is asymptotically flat ) , then there is only one solution of that stays bounded as @xmath24 .
if there exists only one solution of the equation , such that it is linearly independent from the solution bounded at infinity and in the same time it is bounded at @xmath25 then these two solutions uniquely define the green s function .
if all the solutions of are singular at 0 , there is no solution ( bounded in @xmath21 ) of the given initial value problem such that it can be laplace transformed ( this can be taken as an indication that there is no solution at all ) .
if both of the linearly independent solutions are regular at 0 , and ( at least ) two different green functions lead to a function in the domain of the inverse laplace transform , then there is no _ uniquely _ defined solution to the initial value problem . in the `` worst '' case
there are infinitely many solutions , given by arbitrary linear combination of @xmath22 , that are linearly independent to the solution bounded at @xmath24 .
in such a case the problem is underdetermined and one needs one more condition at @xmath26 that selects a unique green s function between the different green s functions marking different time evolutions . for each one of the choices of the green s function
, one can reproduce the calculation from @xcite and see that the quasinormal modes defined by    1 .   the choice of the green s function close to 0 , 2 .   the outgoing radiation condition ,    characterize the time evolution of the field at a fixed point within some specific time interval .
unfortunately for the case of r - n naked singularity ( @xmath27 ,   @xmath28 ) both of the linearly independent solutions @xmath22 are regular at 0 and the problem is underdetermined .
it is easy to show that both of the solutions are regular at the origin .
write as : @xmath29\frac{d u_{l1,2}(s , r)}{dr}\nonumber\\ & + & \left[s^{2}+\frac{l(l+1)f(r)}{r^{2}}\right ] u_{l1,2}(s , r)=0.\end{aligned}\ ] ] now taking the @xmath25 limit of the equation one obtains the following @xmath30\frac{d u_{l1,2}(s , r)}{dr}\nonumber \\&+&\frac{l(l+1)}{q^{2 } } u_{l1,2}(s , r)=0.\end{aligned}\ ] ] this means the solutions @xmath31 behave close to 0 as @xmath32 with @xmath33 hence both are regular .
the fact that the problem is underdetermined is not surprising , since the space - time is _ not _ globally hyperbolic and anything can fall out at any time from the singularity .
this means the singularity has `` hair '' ( carries some other information beyond the metric ) and the quasi - normal modes obviously depend on the  hair " . let us add here that , as we already mentioned in the introduction , the `` hair '' of the singularity relates to the existence of many different self - adjoint extensions of the `` hamiltonian '' operator in the equation .
( in @xcite one can find a nice analysis of the uniqueness of the self - adjoint extensions of such operators for many different types of naked spacetime singularities including the r - n naked singularity . )
is there any intuitive physical condition that we can further impose on the fields , that will uniquely select the appropriate green s function ?
at least to get the geometrical optics continuous extension of the black hole case one can impose the condition that nothing falls in or out of the singularity .
this means there is neither absorption nor superradiation in the scattering and the s - matrix of the k - g field is a unitary operator .
what does this condition mean ?
the conserved current 4-vector for the complex klein - gordon field is given by @xmath34 let us integrate the 4-current along a cyllindrical hypersurface given by @xmath35 and @xmath36 $ ] .
then since is a conserved current it holds that @xmath37    here @xmath38 is the integral along a hypersurface given by the interior of the cylinder at the constant time ( @xmath39 or @xmath40 ) , with a future orientated surface normal vector .
@xmath41 is a cylindrical hypersurface given by @xmath35 and @xmath42 is an induced density given as @xmath43 .
furthermore @xmath44 is a normal vector to @xmath41 given as @xmath45 .
since we want @xmath46 to remain constant with respect to time as we take the limit @xmath47 ( nothing flows out or into the singularity ) , the equation reduces to the following    @xmath48 = 0.\ ] ]    let us analyse the equation in decomposition into spherical harmonics , hence @xmath49 then we can rewrite the equation as @xmath50_{|r = r_{0 } } = 0,\end{aligned}\ ] ] for arbitrary @xmath51 .
since @xmath52 as @xmath47 , in order to fulfill the equation we impose for every @xmath3 and every @xmath53 @xmath54_{|r=0}=0.\ ] ] but this means that the functions @xmath12 should be always constrained either by the condition @xmath55 or @xmath56 .
( one might argue that it will be enough to claim that the fields and their first @xmath21 derivatives should be real at zero , but that does not put any general constraint on the normal modes . )
now we obtain the function @xmath12 from the normal modes as @xmath57,\ ] ] and take the condition @xmath55 for arbitrary time @xmath53 .
this translates to @xmath58=0.\ ] ] but since both @xmath59 ( see ) and fourier - like transform given by should not map non - zero functions to zero , one obtains the condition @xmath60 .
this means we compose the relevant wave packet only from the following modes : @xmath61 the same line of reasoning applies to the condition @xmath56 and the wave - packets that fulfill such condition must be composed entirely from modes given as ( @xmath62 ) : @xmath63 for @xmath64 , we have the coefficient @xmath65 and the solution @xmath66 behaves as a constant for @xmath67 ( with the value set to 1 ) , which means that the condition does not make strictly sense ; in such case the modes are given simply by the function @xmath66 .
this shows that the condition of s - matrix being a unitary operator gives additional constraints on the green s function .
the previous conditions mean that we shall consider only a linear space of wave - packets formed purely from modes that vanish at zero , or a linear space of wave - packets formed purely from modes whose first @xmath21 derivatives vanish at zero .
( of course one can not superpose wave - packets formed from modes having the vanishing radial derivative at 0 with wave packets formed from modes that are vanishing at 0 , since the 4-current does not linearly depend on a wave - function . )
so let us pick one of those two types of modes , fix the particular green s function and employ the following reasoning : since @xmath68 is conserved and the coefficients in the equation are everywhere outside 0 a smooth function of @xmath21 , the function @xmath69 , arising from compactly supported initial data could be unbounded in time only in the case the wave - packet becomes concentrated around the singularity and slowly growing asymptotically with time into delta function @xmath70 but this scenario is prevented by the boundary condition that `` nothing flows in or out of the singularity at any time '' .
this means we expect the solution to be bounded with respect to both time and space and all the reasoning based on the assumption of the existence of laplace transform is justified .
further in the text we will employ the field vanishing condition at 0 .
( this boundary condition at the singularity corresponds to what is known as friedrich s extension of a symmetric operator . )
thus the quasinormal modes will relate to the scattering problem following from the time evolution determined by the boundary condition @xmath71 . but everything that we will do in the following text can be repeated for any other meaningful normal modes boundary condition giving another time evolution of the fields . the general dependence of the quasi - normal modes ( characteristic oscillations ) on such a boundary condition ( given by some linear combination @xmath72 ) describes the way the quasi - normal modes depend on the  hair " of the singularity . since the quasi - normal modes are the ones that carry astrophysical information about the astrophysical sources , their dependence on the  hair " of the singularity ( in some simplifying sense given by the complex parameter @xmath73 ) has a potential astrophysical importance ( of course in case naked singularities have any astrophysical importance ) .
using @xmath74 defined as @xmath75 and @xmath76 the tortoise coordinate given by the condition : @xmath77 one can rewrite the equation into the following form @xmath78 with @xmath79 f\big(r(x)\big).\ ] ] and , for the normal modes  @xmath80 , we can write @xmath81\phi_{l}(x)=0.\ ] ]    if @xmath82 , we can see that @xmath8 given by eq .
( [ tortoise ] ) has no zeros for real arguments , but eq .
( [ x ] ) can still be integrated to give @xmath83 where @xmath84 is an integration constant .
we remark here @xmath85 as @xmath86 , but for @xmath25 we have @xmath87 as can be seen in fig .
( [ fig : rstar ] ) .
further in the text we take the tortoise coordinate with the boundary condition @xmath88 .
it means that with respect to the usual ( @xmath89 ) tortoise coordinate such tortoise coordinate is shifted to the origin by the transformation @xmath90 .    ) for a spacetime with @xmath91 and different values of @xmath92 , taking @xmath93 . ]
the potential has for the ratio @xmath94 less than approximately @xmath95 and the relevant @xmath76 ( in the naked singularity case the domain of @xmath76 is constrained ) 3 extrema , one smaller `` outer '' maximum , one dominant `` inner '' maximum and minimum in the potential valley between them .
( for @xmath25 the function @xmath96 . ) for @xmath94 more than approximately @xmath95 the potential has only one maximum ( thus only one peak ) .
these features of the potential can be seen in the figure [ fig : pot ] .
moreover , in the case @xmath94 less than approximately @xmath95 and @xmath97 , the inner maximum becomes completely dominant , making the outer peak negligible as compared to the size of the inner peak .
for the outer peak and for @xmath97 the terms @xmath98 represent only a small correction ( as compared to the term proportional to @xmath99 ) and the second peak will vanish in this approximation at @xmath100 .
( this is because the outer peak lies always at @xmath101 .
it is also quite obvious that the inner peak exists only due to the fact that the terms @xmath98 become dominant for @xmath21 close to 0 . ) the features described in this paragraph can be observed in the figure [ fig : pot_l ] .
further in the text we will use the following notation related to the potential parameters : by @xmath102 we mean the heights of the two peaks .
( @xmath103 is the first larger peak , @xmath104 the second smaller peak . ) by @xmath105 we call the curvature of the first peak and by @xmath106 the curvature of the second peak .
furthermore , by @xmath107 we mean the point of the location of the top of the first peak and by @xmath108 the point of the location of the second peak .
many of these parameters can be exactly calculated in the @xmath97 limit .
the results of these calculations are given in the appendix [ b ] .
given by eq .
( [ potential ] ) with @xmath109 , @xmath91 for @xmath110 and 0.56 .
( the curves from left to right correspond to the increase of charge . )
note that the dashed part of the potential ( for @xmath111 and 0.5 ) is inside the black hole horizon . ]
first let us have a look at the geometrical optics ( eikonal ) limit of the klein - gordon equation , which is valid for @xmath97 .
in particular what does the geometrical optics limit tell us about the fundamental mode ?
the effective potential for the massless particle is @xmath112 the effective potential goes to plus infinity as @xmath25 and its extremes are given as @xmath113 this means the potential has one minimum and one maximum for @xmath114 for @xmath115 there are no extrema .
the maximum represents the unstable orbit of a massless particle , and for the black hole case it is only the maximum that is relevant ( because it is located above the horizon ) . for the naked singularity with @xmath116
there exists also a stable orbit of a massless particle .
however the fundamental mode in the geometrical optics limit is related to the unstable orbit as @xmath117 ( for a very good paper that discusses this topic see @xcite . )
note that here @xmath118 gives the unstable orbit decay rate and furthermore holds the following @xmath119 ( again @xmath76 is the tortoise coordinate . ) particularly for the reissner - nordstrm black hole / naked singularity ( with @xmath120 ) equation can be expressed as :    @xmath121    with @xmath122 , given as before .     given by eq .
( [ potential ] ) with @xmath91 for @xmath123 ( @xmath124 ) and different values of @xmath3 ( same caption as on the right plot ) .
note in this case the existence of the secondary peak , reminiscent of the black hole potential .
below : the same , but this time for @xmath125 ( @xmath126 ) .
there is no secondary peak in this case.,title="fig : " ]   given by eq .
( [ potential ] ) with @xmath91 for @xmath123 ( @xmath124 ) and different values of @xmath3 ( same caption as on the right plot ) .
note in this case the existence of the secondary peak , reminiscent of the black hole potential .
below : the same , but this time for @xmath125 ( @xmath126 ) .
there is no secondary peak in this case.,title="fig : " ]    from the black hole qnm - s boundary conditions one can see that such massless particles decaying from their unstable orbit and eventually escaping either to the b - h horizon , or the infinity represent the fundamental qnm - s in the @xmath97 ( geometrical ) limit
. one can observe ( see the discussion later ) that the same holds in the naked singularity case . in the @xmath97 case the fundamental mode of the massless perturbations
is described by the picture of a particle decaying from its unstable `` photon '' orbit .
( note that in the naked singularity case the stable photon orbit makes the particle decaying inwards oscillate around the stable orbit . )
this should suggest that in the geometrical optics limit    1 .
one shall expect continuity of the fundamental mode ( as a function of @xmath127 ) when turning from the black hole to the naked singularity case , 2 .   for @xmath128
there is a clear indication that there do not exist low damped quasi - normal modes .
the way how to understand the continuity ( for @xmath97 ) of the fundamental mode as one transfers from the black hole case to the naked singularity case could be the following : consider the partial s - matrix @xmath129 . within
the s - matrix one can look for poles on the complex plane : obviously any of the poles must be necessarily non - real . for the relativistic time dependence convention chosen as  @xmath130   the poles with @xmath131 must be purely imaginary ( as a result of the given hamiltonian being a symmetric operator ) and they represent bound states with the non - relativistic energies given as @xmath132 .
the poles with @xmath133 represent quasinormal frequencies .
some of the quasinormal frequencies correspond to the resonances in the phase factor related to the partial s - matrix , and the wave packets formed out of the resonance energies represent quantum particles tunnelling out of the potential valley ; for the non - relativistic quantum particles ( for example ) with the decay rate given as @xmath134 . now in the naked singularity case the fundamental mode is represented by a resonance that has a low enough energy ( at the level of the smaller peak ) and gets trapped in the potential valley in between the peaks for some time and then radiated away .
the larger peak behind the smaller peak is hugely dominant ( large @xmath3 ) and effectively acts to the wave packet as an infinite barrier .
there are probably much higher resonance energies related to higher qnm frequency overtones that are determined by the higher peak , but left unaffected by the details of the small peak .
so there are two different effective regimes for the quasinormal modes , in the first regime they are sensitive only to the details of the smaller peak , in the second regime they are sensitive only to the details of the larger peak .
this splitting disappears in the case of small @xmath3-s since the two peaks are of comparable heights and effectively interfere .
let us finish this section with one more remark : in the black hole case one can relate the fundamental mode to the peak of the potential ( independently of @xmath97 ) following the way of thinking in @xcite .
the ratio of the amplitudes of the reflected and transmitted waves @xmath135 is for quasinormal modes 1 .
the same happens for real @xmath136-s , when the energies @xmath137 are close ( `` almost at '' ) to the peak of the potential .
so for `` almost '' real @xmath137 ( fundamental mode ) the continuity of the ratio of reflected and transmitted amplitudes analytically continued to the complex plane suggests that the real part of the frequency should be near the square root of the peak of the potential .
( this intuition is then to some extend confirmed also semi - analytically , such as by the ferrari and mashhoon approach using the poeschl - teller potential @xcite , or by using the inverted harmonic oscillator potential . ) in the case of a naked singularity life is not so easy ( in the reduced 1d problem we have no transmitted amplitudes ) , so this simple logic fails .
( one can say luckily it fails as if it remained valid one would also for @xmath97 expect a discontinuity in the fundamental mode as one jumps from the black hole to the naked singularity .
this is because the second , dominant peak suddenly appears in the domain of the @xmath76 coordinate , due to the fact that the domain of @xmath76 discontinuously jumps when passing from the black hole to the naked singularity case . )
certainly one can say that any of the quick fits , such as were done in the black hole case do not lead to anything close to the numerical data obtained in this paper .
in this section we want to suggest some ways to model the naked singularity scattering ( with a general time evolution ) by solutions of the equation with @xmath138 replaced by analytically solvable potentials .
the reasons are the following : first , this section serves as a basis for the analysis of the  @xmath97   cases provided in the section [ largelsection ] .
second , it demonstrates that the naked singularity scattering is _ in principle _ treatable via analytical approximations and thus shows the general power of analytical techniques .
( also the approximations obtained here might produce some future results for the quasinormal frequencies .
for some further calculations see the appendix [ a ] .
furthermore , for a very nice overview of the results for the quasi - normal frequencies and related transmission resonances of the analytically tractable potentials see @xcite . )
third , it might bring more insights into the physics obtained through the exact scalar field potential .
the key point is to split the scalar field potential into different domains and approximate it on each of those domains ( or directly the solution of the equation ) by a different , analytically solvable potential .
then one has to impose the standard procedure : the logarithmic derivatives of the solutions on different domains must be glued on the domain s boundary .
we might add here that since the analytically solvable potentials are typically fitted by the parameters of the peaks one expects such approximation to work for the modes that are _ not _ too damped ( the _ low _ damped qnm modes ) .
the approximation _ can not _ be taken too seriously for the _ highly _ damped modes .
let us suggest modeling all the solutions of the relevant cases as follows :      for this case it might be interesting to consider the morse potential .
( for the definition and origins of all the potentials used in this section see also @xcite . ) on the other hand the infinite valley behind the one and only peak may not be well modelled by only a finite valley that is behind the peak of the morse potential .
thus we consider the solution in the region ` behind' morse potential region to be given as for @xmath139 and the logarithmic derivative gluing condition originates from :    1 .
@xmath140  : use function given as ( for the details see equation ): @xmath141 2 .
@xmath142  : use the appropriate solution of the equation with the morse potential : @xmath143    here @xmath144 is a point most conveniently chosen where the morse potential is 0 ( for any arbitrary @xmath3 there is such point for @xmath145 ) .
this means @xmath144 can be analytically given as @xmath146 let us make here one remark : instead of approximating the original scalar field equation for @xmath67 one can consider taking the equation in the @xmath25 approximation .
this leads to the following : @xmath147 the equation has solutions @xmath148 .
the solution such that fulfills @xmath149 must have @xmath150 .
this approximation of the solution we think to be less exact as the approximation given as , as it is in fact effectively only the first non - zero term of the power series expansion of taken at 0 .
( for such reasons we decided to use in our analytical approximations the approximate solution . )
let us consider now the limit of large @xmath3 . in this case
, since the peak already lies close to 0 and the change in the highest power of inverted @xmath76 becomes more and more dominant ( with higher @xmath3-s ) , one might simplify the approximation given in case [ one ] , by using the following potential :    1 .
@xmath151 : use function given as 2 .
@xmath152 : use the appropriate solution of the equation with the potential ( poeschl - teller potential ) : @xmath153.\ ] ]    here we also used the poeschl - teller potential instead of the morse potential , because the logarithmic derivative gluing condition is much easier to solve for the solutions of the poeschl - teller than for the solutions of the morse potential .
the expense of the simplicity is that the approximation by such potential might be slightly less exact than by using the morse potential , but one still expects it to be accurate enough .
if one does not want to neglect the second smaller peak ( so one is interested in resonances related to the valley between the two peaks ) , the analytically treatable potential describing all the resonances ( and also higher damped qnm frequencies ) could be :    1 .
@xmath140 : use function given as , 2 .
@xmath155  : solutions of the equation with the morse potential @xmath156 3 .
@xmath157 : solution of the equation with the poeschl - teller potential : @xmath158.\ ] ]    here @xmath144 is best chosen as in the case [ one ] and @xmath159 is chosen to be such that the resulting potential is continuous .
generally @xmath159 has to be obtained numerically , but for larger @xmath3-s ( certainly @xmath160 is more than enough , as we checked ) , the curvature of the second peak is very small comparing to the curvature of the first peak and also the curvature of the first peak becomes ( for large @xmath3-s ) very large comparing to the scales of the potential , so that one can calculate ( for @xmath160 with a good approximation at least to the 6 decimal places ) @xmath159 just by    1 .   taking the poeschl - teller potential
to be constant and given by the hight of the outer peak 2 .
simplifying the morse potential by the following approximation : @xmath161    the resulting formula then becomes @xmath162 then one can easily show that @xmath163 as @xmath164 .
now consider what happens with the resonances related to the valley between the two peaks in the large @xmath3 limit .
in such case the dominant peak grows to infinity as compared to the smaller peak , the infinite valley behind the dominant peak shrinks to 0 .
also @xmath165 and the valley between the two peaks becomes flat as compared to the difference between the height of the first peak and the value of the potential at the bottom of the valley .
if the resonances related to the smaller peak locate close to the top of the peak ( they should since the first peak grows with @xmath3 an makes the tunnelling harder ) , then one might effectively approximate the case [ three ] by the poeschl - teller potential : @xmath158,\ ] ] for @xmath166 .
( one can also derive this approximation straight from the logarithmic derivative gluing condition in the appendix [ section ] by taking the @xmath164 limit . )
in this section we want to analytically confirm the results obtained directly from the eikonal limit .
the case [ four ] from the previous section is easily solvable .
the solution in the region @xmath168 must fulfil the outgoing radiation condition . ( to be exact after analytically extending the solution to the complex plane it must give asymptotically the outgoing waves on the line @xmath169 . )
such solution is given as follows    @xmath170    where @xmath171 is the standard hypergeometric function ( of the type 2 - 1 ) and @xmath172 the boundary condition at 0 gives @xmath173 , leading to @xmath174    unless @xmath94 is not too close to the upper limit given for @xmath3 large as 9/8 , one can reliably approximate @xmath175 and the condition turns to be : @xmath176    now can be rewritten through gamma functions and reduces to the simple problem of finding poles of the product of gamma functions @xmath177\gamma\left[-\frac{i\omega+i\sqrt{v_{2}}}{\alpha_{2}}+\frac{1}{2}\right]=0.\ ] ] the poles are located at @xmath178 for @xmath179 being a natural number and this gives @xmath180 this is precisely the formula for the lowest damped black hole qnm frequencies .
this means in the @xmath97 limit we see ( at least for the lowest modes ) a continuous transition from the black hole case to the naked singularity case , as expected .
( at least in the case where it holds that @xmath175 , which becomes a less accurate approximation for @xmath94 close to @xmath95 .
on the other hand for such ratios of @xmath94 the peak is almost vanished and one can assume that the decaying circular orbit eikonal picture has already quite limited sense . )    all this means that in the regime @xmath97 the lowest qnm frequencies should be located close to the first peak of the potential and are ( up to certain @xmath179 ) given by the formula .
they represent resonances such that they are insensitive to the details of the larger peak as in their case the larger peak can be already seen to act effectively as an infinite potential barrier .
hence we can conclude that in the @xmath97 case we are able to match the geometrical ( eikonal ) limit of the original problem considered .
take the case [ two ] from the previous section . on the left hand side of the logarithmic derivative gluing condition
we need a logarithmic derivative of the function : @xmath181 on the right hand side of the logarithmic derivative gluing condition we need the logarithmic derivative of the solution on the interval @xmath182 $ ] .
it is given as in , only the parameters are related to the first and only peak in this case .
then the logarithmic derivative gluing condition leads to the following :    @xmath183.\ ] ]    if one takes the @xmath164 limit then the equation becomes : @xmath184.\ ] ]    here @xmath185 is defined as @xmath186 . for not highly damped frequencies in the @xmath97 limit it must hold that @xmath187 . but by plotting the right side of the equation minus the left side of the equation in mathematica it can be shown that the equation does not have any real solutions . for more than large intervals of values of @xmath185
the absolute value of the left hand side of the equation minus the right hand side of the equation always seems to grow ( almost ) linearly with respect to @xmath185 from the value approximately given as 0.35 at @xmath188 to infinity .
it might not be hard to prove also analytically that does not have real solutions by using many relatively simple properties of gamma functions .
this means there are _ no _ low damped modes in the case of @xmath189 and @xmath97 , as expected from the eikonal limit / particle picture .
( it means the damping of the fundamental mode grows to infinity as @xmath164 and grows no less rapidly than @xmath105 , hence no less rapidly than cubically with @xmath3 . )
in @xcite , the quasi - normal modes of a scalar field in an electrically charged vaidya background were studied using a similar numerical setup to what we are going to use here .
an interesting investigation done in @xcite tried to determine what would happen with the quasi - normal modes as the time - dependent background approached a naked singularity , but the numerical code used was not suitable for following the field evolution after the extremal configuration ( @xmath190 ) was reached . in our present work , we have a static configuration that describes a naked a singularity in order to study the properties of the quasi - normal modes .
our objective in this section is to solve eq .
( [ wave ] ) with potential  ( [ potential ] ) numerically , in the case where @xmath191 as described in the last section . to do this ,
we rewrite eq .
( [ wave ] ) in terms of the light - cone variables @xmath192 and @xmath193 , where @xmath76 corresponds to the tortoise coordinate ( [ x ] ) , as @xmath194 that can be integrated with the boundary conditions @xmath195 where condition  ( [ bc1 ] ) is a necessary condition on the field @xmath196 near the origin ( see the discussion on fig .
[ fig : phi_r ] below ) and condition  ( [ bc2 ] ) defines an `` arbitrary '' relevant initial signal to be propagated .
we use the algorithm @xmath197 see where @xmath198 and @xmath199 are the integration steps in @xmath200 and @xmath201 , respectively , and the definitions of @xmath202 etc can be seen on fig .
[ fig : grid ] . note that here @xmath203 is the potential ( [ potential ] ) evaluated at the same @xmath21 coordinate as @xmath204 ( and @xmath202 ) .    ) with boundary conditions ( [ bc1 ] ) and ( [ bc2 ] ) .
the points marked with `` x '' are out of our domain , values of @xmath196 at the positions marked with filled circles are given by the boundary conditions and the values at the empty circles are obtained with the algorithm ( [ algo ] ) . ]    as we can see in fig .
[ fig : phi_r ] , the boundary conditions ( [ bc1 ] ) and ( [ bc2 ] ) ensure the necessary conditions on the fields @xmath196 and @xmath205 near the center . as we discussed previously in section [ sec2 ] ,
the physically correct boundary condition for @xmath205 is @xmath206 . from this
we must have for @xmath207 that @xmath208 and @xmath209 .
in fact , the quantity @xmath210 dx\ ] ] is invariant along the @xmath53-evolution governed by ( [ waveuv ] ) with the boundary conditions @xmath211 . taking into account ( [ x ] )
, the integrand can be expressed near @xmath26 as @xmath212 with the prime denoting the derivative with respect to @xmath21 . any initial condition with finite @xmath213 , as those ones obeying ( [ bc1 ] ) and ( [ bc2 ] ) , is such that @xmath214 for @xmath25 , with @xmath215 , implying that @xmath216 for initial condition with finite e obeying @xmath217 .
moreover , we have checked that our boundary conditions for @xmath196 also reproduce the @xmath209 condition , by testing the code with the condition @xmath218 , where @xmath219 is a very small value , and we found no differences in the late time evolution of the fields .     with @xmath220 as a function of @xmath21 near the center @xmath26 for a late time @xmath221 , shown here in order to exemplify the effect of conditions ( [ bc1 ] ) and ( [ bc2 ] ) in the numerical integration , for a spacetime with @xmath222 and @xmath123 .
below : the same as in the left plot , but this time for the function @xmath223.,title="fig : " ]   with @xmath220 as a function of @xmath21 near the center @xmath26 for a late time @xmath221 , shown here in order to exemplify the effect of conditions ( [ bc1 ] ) and ( [ bc2 ] ) in the numerical integration , for a spacetime with @xmath222 and @xmath123 .
below : the same as in the left plot , but this time for the function @xmath223.,title="fig : " ]    in the left plot of fig .
[ fig : field_qm ] we present some typical time evolutions of @xmath196 , for a @xmath220 and different @xmath224 ratios . in the right plot
we present the obtained frequencies of the qnms ( fundamental mode ) in the @xmath225 plane .
we can see a discontinuity in the frequencies as @xmath226 , as was expected from the discussion of the potential @xmath227 ( see fig .
[ fig : pot ] ) .
we also point here that we see no significant changes , but rather a smooth behavior as @xmath228 ( @xmath229 in the plot ) .
but we see a point of inflection in @xmath230 at @xmath231 , for which we did not find an analytical explanation .
we remark here some numerical issues that prevented us from extending the results shown in fig .
[ fig : field_qm ] for lower and higher @xmath224 rations . as @xmath226 and @xmath196 becomes less damped , different modes subsist for longer times and and longer evolutions are needed in order to obtain the clear frequency of the fundamental mode . on the other hand , as we increase @xmath224 , @xmath196 is damped so quickly that we can not observe enough oscillation cycles to obtain the frequencies .
these are the issues that have limited our results .
with @xmath220 at @xmath232 for a spacetime with @xmath233 and different values of @xmath234 . below : frequencies of the fundamental mode with @xmath220 in the @xmath225 plane , parametrized by the @xmath224 ration.,title="fig : " ]   with @xmath220 at @xmath232 for a spacetime with @xmath233 and different values of @xmath234 .
below : frequencies of the fundamental mode with @xmath220 in the @xmath225 plane , parametrized by the @xmath224 ration.,title="fig : " ]    finally , in fig .
[ fig : field_l ] we explore how the frequencies of the qnms change with @xmath3 . as usual in black hole scattering problems , we see that the oscillation frequency @xmath230 increases with @xmath3 . but the qualitative behaviour of @xmath235 changes significantly with @xmath224 . in the upper plots ( @xmath236 )
, @xmath237 decreases with @xmath3 , that is , the damping time is longer . in the lower plots ( @xmath2 )
, we have the opposite tendency .
this behaviour is connected to the potential @xmath227 shown in fig .
[ fig : pot_l ]
. it might be also interesting to mention that in case @xmath2 , there is a qualitative similarity in the behaviour of the imaginary part of the frequencies as a function of @xmath3 , between the case when @xmath3 is small and the @xmath3 large limit .     for @xmath91 and @xmath123 ( @xmath124 ) .
below : same as above , but this time for @xmath91 and @xmath125 ( @xmath126 ) . ]
let us explore the following limit , @xmath3 will be fixed and @xmath238 , assuming that there exist an infinite number of qnms .
assume further that @xmath239 diverge as @xmath240 .
in such a limit , to have a good understanding of what the qnm boundary conditions are , it is convenient to take an analytic continuation of the solution into the complex @xmath76 plane and impose the outgoing wave boundary conditions on the ( stokes ) line @xmath241 .
this is for almost imaginary asymptotic frequencies ( @xmath242 as @xmath240 ) the line @xmath243 .
the key point is to realize that , in the asymptotic case , @xmath244 holds everywhere apart of a tiny region around @xmath245 , where the potential @xmath246 can be already approximated by @xmath247 ( there is one subtlety : the potential , and also later the solutions , are multi - valued functions in the complex variable , but after specifying a branch cut together with the branch , this does not cause any problems to our analysis . ) so let us make the following statement : the solutions of the equation @xmath248 are _ everywhere _ good approximations to the solutions of the scalar field equation for the asymptotically highly damped modes .
( the region where the approximations @xmath249 and @xmath250 cease to hold is the region where we can already neglect the potential as a whole with respect to the @xmath137 term . )
the equation has a general solution given as @xmath251~,\ ] ] where @xmath252 is the bessel function .
the solution for which @xmath253 fulfils the vanishing boundary condition at 0 is given by @xmath254 .
our qnm boundary condition says that this solution taken along the stokes line @xmath241 , should give purely outgoing radiation for @xmath255 .
but in the @xmath256 limit we know that we can approximate the solution ( via bessel functions approximations in such a limit ) as @xmath257 this linear combination is independent on @xmath136 and does not give the purely outgoing radiation ( @xmath258 ) , which suggests that our qnm boundary conditions can not be fulfilled for asymptotically highly damped modes .
this can be taken as an evidence that the asymptotically highly damped modes ( @xmath259 ) do _ not _ exist for the r - n naked singularity and this can be taken as a confirmation that the highly damped qnms could link to the black hole spacetime horizon s properties ( non - existence of the horizon leads to non - existence of the modes in the asymptotic limit ) .
one can see also some link to the fact that for the extremal r - n black hole the spacing between the asymptotic frequencies goes to zero as proportional to the surface gravity ( this might make a whole infinite tower of arbitrarily highly damped modes eventually `` collapse '' to a single mode ) .
moreover , we assume that the non - existence of the highly damped modes can be proven along the same lines also for another naked singularity spacetimes ( like negative mass schwarzschild spacetime etc @xcite ) .
in this paper we analysed the problem of the scalar field scattering on a r - n naked singularity background from the point of view of quasi - normal modes .
the evolution on the r - n naked singularity is non - unique unless one specifies an additional boundary condition representing a `` hair '' of the singularity .
the quasi - normal modes then carry information about the `` hair '' .
we applied a particular boundary condition , that nothing comes out , or in from the singularity and analysed analytically , as well as numerically , the characteristic oscillations of the scalar field perturbations ( low damped quasi - normal modes ) .
we analysed the eikonal @xmath97 case via the analytical approach confirming the intuition obtained through the massless particle viewpoint , and showed that an approach based on analytical approximations can be useful also for the small @xmath3 wave mode numbers .
for the small @xmath3-s we calculated the frequencies numerically via the characteristic integration method .
we also suggested arguments showing that the asymptotically highly damped modes ( limit @xmath3 fixed and @xmath260 ) do _ not _ exist in case of r - n naked singularity .
this might confirm the intuition one has about such modes from the black hole physics ( and from the presently popular conjectures @xcite ) .
the basic results can be summarized as follows : for the low modes and the large @xmath3 there is a continuous transition in the low damped qnm modes between the r - n black hole and the r - n naked singularity . however ,
when the ratio @xmath94 becomes larger than approximately @xmath95 then the picture becomes significantly different and the low damped modes do _ not _ exist for large @xmath3-s .
( this is a very different picture from the bh based intuition . ) for the small @xmath3 numbers the modes face a discontinuous transition when transiting from the black hole to the naked singularity .
furthermore , the @xmath3 dependence @xmath237 ( for small @xmath3 ) changes as @xmath94 becomes larger than approximately @xmath95 : @xmath237 decreases for @xmath261 and increases for @xmath262 .
it might be interesting to notice that for @xmath262 the increase of @xmath263 as a function of @xmath3 ( for small @xmath3-s ) matches the behaviour of @xmath263 for large @xmath3-s . in the case of large @xmath3-s and @xmath262
we have shown that @xmath263 of the fundamental mode grows at least cubically with @xmath3 and thus , as we already mentioned , the low damped modes do not exist .
for the asymptotically highly damped modes our results seem to suggest that they do not exist , which means that the imaginary parts of the frequencies are bounded .
the authors are grateful to fapesp , cnpq , and the max planck society for the financial support .
in this appendix section we provide a list of some of the quantities ( also some of their derivations ) that occur in the calculations relevant to the problem analysed in this paper .
write @xmath227 in a convenient form : @xmath264 with @xmath265,\\ c&\equiv & q^{2}l(l+1)-(2m)^{2}-2q^{2},\\ d&\equiv & 6mq^{2},\\ e&\equiv & -2q^{4}.\end{aligned}\ ] ] now let us calculate the second derivative of the potential with respect to the tortoise coordinate : @xmath266 at @xmath267 the derivative @xmath268 is zero , then necessarily @xmath269_{|x_{max}}\ ] ] and since @xmath270 we obtain @xmath271 for the global extremum condition holds : @xmath272    this equation can have maximally 4 roots . from it
can be easily seen that for @xmath24 becomes for any @xmath273 positive ( dominant term is @xmath274 ) and goes to 0 . in case of naked singularity for @xmath25 the potential goes always to @xmath275 ( dominant term is @xmath276 ) , in case of black hole for @xmath277 it always goes to 0 ( dominant term is @xmath8 ) .
this means , together with , that for any @xmath278 ( black hole or not ) there are always either two local maxima and one local minimum , or one local maximum without local minima .
for @xmath97 the following approximations hold : for large @xmath3 the maximum of the inner peak can be approximately found analytically and behaves as : @xmath279 also holds the following : for @xmath280 @xmath281 the parameters @xmath282 turn in the large @xmath3 limit to be : @xmath283^{3/2}}{|q|},\ ] ] @xmath284^{3/2}}{3^{3/2 } |q|}.\ ] ] this implies also the following results : for @xmath164 @xmath285
write the two linearly independent solutions on the domain @xmath286 in the form convenient for the logarithmic derivative gluing condition : @xmath287 . if the domain extends to the infinity and @xmath288
is taken to be the asymptotically incoming wave solution and @xmath289 the asymptotically outgoing wave solution , then @xmath290 is the s - matrix ( ratio of the coefficients of the outgoing and incoming waves ) .
anyway , the logarithmic derivative gluing condition at the boundary of the regions @xmath286 and @xmath291 turns to be :      this is the condition connecting the @xmath293 coefficient with @xmath294 coefficient through the values of the solutions and their derivatives at the boundary of the regions .
on one side consider asymptotically free - wave region , whereas on the other side consider infinite potential barrier - like boundary condition for @xmath295 . in such case
one can through equations of the type connect ( after finite number of steps ) the s - matrix of the external ( asymptotic ) region with the boundary condition at the origin ( where we can conveniently put @xmath296 to be zero ) .
then one gets the s - matrix expressed through a function of the boundary values ( at all the boundaries ) of the solutions and their derivatives
. this will be algebraically more difficult condition than the one coming just from posing the outgoing wave condition in the external region , as we did in the previous sections ( case of large @xmath3 and @xmath94 larger than approximately 9/8 ) .
the reason for posing the more complicated condition of the type is that it might be easier ( and also brings more insight ) to look for a real energy resonances ( by plotting phase change and looking for rapid phase shifts , hence rapid change of the logarithm of the s - matrix ) , than to try to numerically solve the more simple outgoing wave condition ( with the help of appropriate computer software ) . such resonances should then 1 to 1 correspond to the low damped quasi - normal frequencies .
let us apply these ideas and try to calculate the expressions for the partial s - matrix .
this is the case [ one ] .
we are gluing only two regions .
the solution on the left side of the boundary @xmath298 , such that fulfils @xmath299 , is again given as .
on the right side ( the morse potential side ) the solution is given as :      by @xmath301 we mean asymptotically outgoing and incoming wave solutions in the sense that @xmath302 and @xmath303 exactly .
furthermore @xmath304 is being the kummer ( confluent ) hypergeometric function and @xmath305 also @xmath306 is the s - matrix .
the logarithmic derivatives gluing condition at @xmath307 is then @xmath308=   \frac{\alpha_{1}z(a_{1})}{2}+i\omega+\frac{s'_{1}}{s'_{2}}\frac{m\left(s'_{1}+1 , s'_{2}+1,z(a_{1})\right)}{m\left(s'_{1 } , s'_{2},z(a_{1})\right)}.\ ] ] considering that the following holds : @xmath309 one can rewrite into : @xmath310=   i2\sqrt{v_{1}}+i\omega+\frac{s'_{1}}{s'_{2}}\frac{m\left(s'_{1}+1 , s'_{2}+1,\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}{m\left(s'_{1 } , s'_{2},\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}.\ ] ]    the partial s - matrix can be obtained from ( @xmath311 ) as : @xmath312 since for real @xmath136 holds that @xmath313 and @xmath298 is a real valued function , one can immediately observe that ( for @xmath314 ) @xmath315 , hence partial s - matrix is unitary .
@xmath316-\left[i2\sqrt{v_{1}}-i\omega+\frac{s_{1}}{s_{2}}\frac{m\left(s_{1}+1 , s_{2}+1,\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}{m\left(s_{1 } , s_{2},\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}\right]}{f(a_{1})\left[\frac{1}{a_{1}}+\frac{\beta_{1}e^{\beta_{1}a_{1}}-\beta_{2}e^{\beta_{2}a_{1}}}{e^{\beta_{1}a_{1}}-e^{\beta_{2}a_{1}}}\right]-\left[i2\sqrt{v_{1}}+i\omega+\frac{s'_{1}}{s'_{2}}\frac{m\left(s'_{1}+1 , s'_{2}+1,\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}{m\left(s'_{1 } , s'_{2},\frac{i4\sqrt{v_{1}}}{\alpha_{1}}\right)}\right ] } \nonumber\end{aligned}\ ] ]      the two peak case is a more complicated problem .
anyway , one can write the logarithmic derivative gluing condition as @xmath317^{-1}\}}{f_{21}\{g_{1},g_{2},g_{3 } , [ 1+\exp(2\alpha_{2}(a_{2}-x_{2max}))]^{-1}\}}. \nonumber\end{aligned}\ ] ] here @xmath318 is the s - matrix from the case [ one ] and is given by the formula . the s - matrix related to the case
[ three ] ( call it @xmath319 ) is related to @xmath318 through the condition . in the gluing formula @xmath319 stands for @xmath293 and @xmath318 for @xmath294 .
also @xmath320 are the solutions of the poeschl - teller case @xmath321^{-1}\right],\ ] ] and @xmath322^{-1}\right].\ ] ]    @xmath323 are the morse potential solutions given as .
we see that if there is a regime in which the modes do not feel the smaller peak , then they are given by the poles of @xmath318 and must be in the same time poles of @xmath319 .
the formula for @xmath319 in this case ( for frequencies that represent poles of @xmath318 ) simplifies to @xmath324_{|_{a_{i}}}.\ ] ] there is a possibility that for the case of the low modes in the two peak model one might be able to neglect the infinite depth of the valley between the infinite barrier and the morse peak . in such case
one can take a simpler model where the morse potential directly follows after the infinite barrier .
the @xmath318 term appearing in the formula will be in such case given by a much more simple expression : @xmath325