during the past several years considerable interest has been shown in exploring localized nonlinear waves in the variable coefficient nonlinear schrdinger ( nls ) equation and its generalizations @xcite .
the motivation comes from the fact that nls equation and its variants appear in several branches / topics of physics , including nonlinear optics @xcite and bose - einstein condensates ( becs ) @xcite , etc . focusing our attention on becs alone
, it is well known that the gross - pitaevskii ( gp ) equation governs the evolution of macroscopic wave function at ultra low temperatures @xcite . in particular , for cigar shaped becs , it has been shown that the gp equation can be reduced to the 1d variable coefficient nls equation @xcite @xmath0 where @xmath1 is the condensate wave function , @xmath2 is the dimensionless normalized time , @xmath3 is the dimensionless normalized coordinate in the axial direction , @xmath4 represents the effective scattering length and @xmath5 is the axial trap frequency .
a simple and straightforward way of exploring the localized / periodic structures of ( [ eq:1d - gp ] ) is by transforming it into a constant coefficient nls equation through a suitable transformation . from the known solutions of the latter equation
the solutions of the former equation can be identified . using this procedure a class of solutions , in particular various localized solutions , have been identified for the model ( [ eq:1d - gp ] ) , including soliton , breather , and rogue wave ( rw ) solutions . in order to appreciate the relevance of the above type of localized structures for eq .
( [ eq:1d - gp ] ) , we may first consider their existence in the case of constant coefficient nls equation , @xmath6 it is well known that the standard nls equation ( [ nls ] ) admits the following basic localized profiles and their higher - order versions @xcite .
[ [ i - envelope - soliton ] ] ( i ) envelope soliton : + + + + + + + + + + + + + + + + + + + + +    it is a solitary wave ( localized envelope along with a carrier wave ) that retains its characteristics ( amplitude / shape and velocity ) under collision with another soliton , except for a change in phase .
the intensity profile of the soliton is shown in fig . [ figa](a ) .
the typical form of the envelope soliton is @xmath7 where @xmath8 is a complex constant , @xmath9/2+\eta_{1i}^{(0)}$ ] and @xmath10 . here
@xmath11 and @xmath12 are constant parameters .
[ [ ii - breather ] ] ( ii ) breather : + + + + + + + + + + + + + +    it is a localized solution with temporally and/or spatially periodic structures having constant background exhibiting internal oscillations and bound states of nonlinear wave packets @xcite which is represented in fig .
[ figa](b ) .
its typical form reads @xcite    @xmath13   \label{nls2}\end{aligned}\ ] ]    where @xmath14 is the amplitude of the plane wave , @xmath15 , @xmath16 , @xmath17 , @xmath18 are constant parameters , @xmath19 , @xmath20 , @xmath21 , @xmath22 , @xmath23 , @xmath24 , @xmath25 , and @xmath26 is a complex constant .
this solution is also called a general breather ( gb ) since it is periodic both in space and in time .
two important special cases are the following : ( i ) when @xmath27 and @xmath28 , the gb solution corresponds to an akhmediev breather ( ab ) which is periodic in space and localized in time .
( ii ) if we take @xmath29 and @xmath30 , ( [ nls2 ] ) becomes a ma breather ( mb ) which is periodic in time only and localized in space .
[ [ iii - rogue - wave ] ] ( iii ) rogue wave : + + + + + + + + + + + + + + + + +    a further specialized structure , which is localized both in space and in time with a constant plane wave background , is the so called peregrine soliton or rw @xcite . it can be obtained by taking the limits @xmath31 and @xmath32 , and @xmath33 , where @xmath34 , @xmath35 and @xmath36 are constants .
( [ nls2 ] ) reduces to the following form , @xmath37 very often the above localized nonlinear wave is described as a wave that  appears from nowhere and disappears without a trace " @xcite as shown in fig . [ figa](c ) .
it was first observed in the area of oceanography @xcite and is traditionally defined as a wave whose wave - height ( the distance from trough to crest ) is more than twice the significant wave height ( swh ) . the latter is generally defined as the average wave height among one third of highest waves in a given time series @xcite .
further , it has also been explained that the above structure ( [ nls3 ] ) arises due to a modulation instability ( mi ) @xcite of the plane wave solution of the constant coefficient nls equation .
very recently higher - order rws ( hrws ) which correspond to the higher - order rational solutions of the nls equation ( [ nls ] ) have been deduced @xcite .
the explicit expressions for the second and third - order rws of the nls equation are given in appendix a. these hrws have higher amplitudes than the first - order rw .
the rws have also been observed experimentally in physical systems such as water wave tank @xcite , capillary waves @xcite and nonlinear optics @xcite .
several theoretical studies on the dynamics of rws in nonlinear fiber optics @xcite , plasma physics @xcite , laser - plasma interactions @xcite , and even econophysics @xcite , described by scalar nls equation , have been made in recent times .
now , as mentioned in the beginning , the dynamics of a cigar shaped bec at absolute zero temperature is usually described by the mean - field gp equation ( [ eq:1d - gp ] ) , which is a generalized form of the ubiquitous constant coefficient nls equation ( [ nls ] ) , for the wave function of the condensate . since the nls equation ( [ nls ] )
admits breather and rw solutions , it is natural to expect that rws and breathers may also be found in bec systems as well . in this context
, the rws can correspond to a sudden increase of peaks in the condensate clouds similar to the nature of high peaks in the open sea , while breathers are generalizations of the rws .
the formative mechanism for the matter rws in becs is the accumulation of energy and atoms towards its central part and their spreading out to a constant density background .
the formation of matter breather is the periodic exchange of atoms between the profile and the plane wave background . from an experimental point of view
, the existence of rw and breather structures in a bec system can be effectively controlled by tuning the nonlinear interaction between atoms by feshbach resonance technique @xcite and modulating the trapping frequency of the external potential .
it will therefore be of great interest to study the characteristics and/or controlling of structures of rws and breathers due to their localization both in space and in time in bec experiments .
past explorations of gp equation in becs have paved the way for important developments in manipulating coherent matter waves for application , including atom interferometry @xcite , coherent atom transport @xcite and quantum information processing or quantum computation @xcite .
therefore , it is of high significance to study the dynamics of rws and breather profiles of the gp equation ( [ eq:1d - gp ] ) .
however only few attempts have been made to identify and analyze the rws and breather solutions of ( [ eq:1d - gp ] ) @xcite . to the best of our knowledge
neither higher - order rw solutions ( with certain free parameters ) nor higher order breather solutions of ( [ eq:1d - gp ] ) have been taken up for study .
motivated by these observations , in this work , we construct the aforementioned localized and periodic solutions of ( [ eq:1d - gp ] ) . besides constructing these two families of solutions
we also investigate how to manipulate the rws and breathers through the effective scattering length and the strength of trap parameter .
having stated the motivation we now proceed to construct a transformation that transforms eq .
( [ eq:1d - gp ] ) to the standard nls equation ( [ nls ] ) . following the standard procedure
@xcite we find that the required similarity transformation should be of the form @xmath38 \notag \\   & \times \exp \left[{i\left(c_1r_0 ^ 2 rx-\frac{r_t}{r}x^2-\frac{1}{2}c_1 ^ 2r_0 ^ 4\int{r^2(t)}dt\right)}\right],\end{aligned}\ ] ] where    [ a4 ] @xmath39    and @xmath40 is the solution of the standard nls equation ( [ nls ] ) . in the above @xmath41
are arbitrary constants and the modulational functions @xmath4 and @xmath42 should satisfy the following condition @xmath43 which is a riccati type equation with dependent variable @xmath44 and independent variable @xmath2 .
regardless of what @xmath4 is , as long as the condition ( [ a5 ] ) is satisfied , the gp equation is integrable @xcite .
we also note here that the painlev singularity structure analysis performed on eq .
( [ eq:1d - gp ] ) confirms the same restriction ( [ a5 ] ) on the system coefficients @xcite .
we further note that the solution ( [ a7 ] ) provides us some flexibility to generate new structures related to the rws which may be useful for the bec experiments .    even though one can arbitrarily choose the functions @xmath4 and @xmath5 that satisfy the constraint ( [ a5 ] ) and generate the required solutions , in this paper we consider the trap frequency to be of the following three forms : ( i ) @xmath45 , ( ii ) @xmath46 $ ] and ( iii ) @xmath47 $ ] , where @xmath48 is a constant .
as shown in reference @xcite the effective scattering length for these three cases turn out to be ( i ) @xmath49 , ( ii ) @xmath50 and ( iii ) @xmath51 .
we consider the trap frequency to be in the above forms since it has been shown that they are valid forms in bec experiments @xcite .
we consider each one of the cases separately and substitute them in ( [ a7 ] ) along with the rw and breather solutions of ( [ nls ] ) .
we then analyze in detail how the nature of the rw and breather structures get modified by the above two functions @xmath5 and @xmath4 .
our analysis shows that the amplitude parameter @xmath52 plays a vital role in the formation of rw , and the trap frequency and the effective scattering length modify the structure of the rw and breather profiles , allowing one the possibility of manipulating the rws and breathers in specific ways .
the paper is organized as follows . in the following section ,
we construct rw solutions for time - independent and time - dependent traps and study their characteristics in detail . to the best of our knowledge , for the first time in the literature , we observe that the second- and third - order rws transform to first - order rw - like structures when we tune a parameter which appears in the harmonic traps ( time - independent and time - dependent traps ) . in sec .
iii , we also construct rw solutions with free parameters which allow us to split the symmetric form solution into a multi - peaked solution for ( [ eq:1d - gp ] ) and investigate how these rw structures get modified in the plane wave background by increasing the strength of the trap . in sec .
iv , we construct one - breather and two - breather solutions of ( [ eq:1d - gp ] ) and investigate their characteristics when we alter the trap parameter . finally , in sec .
v , we present a summary of the results and our conclusions .
to begin with , we consider the case in which the trap frequency is a constant , that is @xmath42 = constant = @xmath53 .
a time - independent trap frequency implies that the frequency does not change with time and space .
we then consider the trap frequency to be time - dependent and investigate the associated rw solutions .
substituting @xmath45 in the integrability condition @xmath54 , we find that the time - dependent interaction term should be of the form @xmath49 , where @xmath35 is an integration constant .        ) through split - step crank - nicolson method for the time - dependent nonlinearity coefficient @xmath49 and time - independent trap frequency @xmath45 .
the initial condition chosen corresponds to the analytic solution of fig .
[ fig8](d ) . ]    plugging this expression in ( [ a7 ] ) , we find @xmath55 where @xmath56\right\ } , \notag\end{aligned}\ ] ] and @xmath57 s , @xmath58 , are the first , second and third - order rw solutions of the nls equation ( [ nls ] ) whose explicit expressions are given in the appendix ( vide eqs .
( [ a8 ] ) , ( [ a11 ] ) and ( [ a12 ] ) ) . also @xmath59 and @xmath60 have the forms as given in eqs .
( [ a4 ] ) .                in fig .
[ fig8 ] , the first , second , and third columns represent the density profiles of the first , second , and third - order rws obtained from ( [ a15 ] ) . in this figure we present the formation of rws in cigar - shaped becs . in the context of becs , it is the fluctuation in the density of atoms , localized in space and time , which is what we observe as rws .
for instance consider the first - order rws depicted in figs .
[ fig8](a)-(d ) .
these may be interpreted as follows : atoms in the condensate , initially at a constant density background , suddenly accumulate to form a hump towards the center of the condensate at finite time while leaving voids in the density which appear as troughs in the rws , depending on the initial state .
the crest and troughs reach their extrema as time progresses , and then the condensate atoms spread out so as to recover the constant density background in finite time , thus revealing the unstable nature of rws . further , by tuning the amplitude parameter @xmath61 in ( [ a15 ] ) , we can visualize the formation and manipulation of rws .
for example , when we increase the amplitude parameter @xmath61 smoothly from @xmath62 to @xmath63 , we can observe in figs .
[ fig8](a)-(d ) the formation of more and more localized first - order rws with increasing amplitude . at @xmath64
one can visualize a large amplitude wave which is sufficiently localized both in space and in time which in turn confirms the formation of the first - order rw in bec .
the formation of second and third - order rws is also demonstrated in figs . [ fig8](c)-(h ) and [ fig8](i)-(l ) , respectively , for the same set of parameter values of @xmath61 . from these figures
we infer that the time - dependent nonlinear interaction between the atoms induces density fluctuations over the condensate which gets more and more localized both in space and in time as we increase the order of the rw .    in order to confirm the existence of the rws further
, we have also performed a direct numerical simulation of ( [ eq:1d - gp ] ) with the aid of the split - step crank - nicolson method using an initial wave function which is the same as the function ( [ a15 ] ) and with space step @xmath65 and time step @xmath66  @xcite . in fig .
[ num - rog ] , we have presented the computer generated density profile of the first - order rw and the corresponding contour plot with the parameters chosen as @xmath64 , @xmath67 , @xmath68 , and @xmath69 which are same as that of fig .
[ fig8](d ) .
one can easily observe a very good agreement between the numerical results and the analytical predictions for the emergence of the rws .
we have also verified numerically the presence of second- and third - order rws of ( [ eq:1d - gp ] ) as well , replicating figs .
[ fig8 ] .
next we demonstrate how these localized structures vary with respect to the trap parameter @xmath48 .
[ fig1 ] displays the first - order rw for the same nonlinearity management parameter @xmath49 and the external trap frequency @xmath45 .
the nature of the first - order rw for @xmath68 is depicted in fig .
[ fig1](a ) .
when we increase the strength of the trap parameter the density profiles corresponding to the first - order rw become more and more localized in time as shown in figs .
[ fig1](b ) and [ fig1](c ) , respectively .
the corresponding contour plots are given in figs .
[ fig1](d)-(f ) .
[ fig1a ] displays the density profiles of the second - order rws for the same nonlinearity management parameter as a function of @xmath48 . in fig .
[ fig1a](a ) we display the second - order rw for @xmath68 .
when the strength of the parameter @xmath48 is increased to @xmath63 , the wave subcrests start to stretch as shown in fig .
[ fig1a](b ) .
the wave subcrests become more and more localized in time when we increase the value of @xmath48 and finally the second - order rw attains a new structure as given in fig .
[ fig1a](c ) for @xmath70 .
the resultant structure looks almost like a first - order rw , see fig .
[ fig1](c ) . figs .
[ fig1a](d)-(f ) are the corresponding contour plots .
we also observe similar effects in the third - order rw case as well when we increase the value of @xmath48 .
[ fig3a ] demonstrates the changes in the third - order rw when we vary the interaction strength . at @xmath71 ,
the third - order rw transforms into the second - order rw - like structure as shown in fig .
[ fig3a](b ) .
when we increase the value of the parameter @xmath48 to @xmath72 , one obtains a first - order like rw which is displayed in fig .
[ fig3a](c ) .
these facts are also confirmed by the corresponding contour plots , figs .
[ fig3a](d)-(f ) .
the aforementioned results reveal that when we increase the strength of the trap parameter @xmath48 , the second and third - order rws become more localized in time and delocalized in space , approaching the structure of a first - order rw .
thus the robustness of the density profiles can be controlled by varying the strength of trap frequencies .
next we consider the time - dependent trap frequency in the form @xmath46 $ ] .                for this choice ,
the relation ( [ a5 ] ) fixes the interatomic interaction term to be of the form @xmath50 .
the first , second , and third - order rw solutions for this trap frequency and strength of interatomic interaction are found to be @xmath73where @xmath74 and    @xmath75}-c_1r_0 ^ 2\left[1+\tanh{\left(\frac{\beta_0 t}{2}\right)}\right]x + \frac{c_1 ^
2r_0 ^ 4\left(\beta_0 t+2\log\left[\cosh(\frac{\beta_0 t}{2})-\tanh(\frac{\beta_0 t}{2})\right]\right)}{\beta_0}\right]\right\}. \notag\end{aligned}\ ] ]    the qualitative nature of the first , second , and third - order rws for @xmath76 and @xmath77 $ ] turns out to be the same as in the previous case ( fig .
[ fig8](d ) , ( h ) and ( l ) ) when the amplitude parameter @xmath61 is varied and so we do not display the outcome here separately . on the other hand , we identify interesting structures while varying the parameter @xmath48 which is discussed in the following .    in fig .
[ fig4 ] , we depict the first - order rw for these choices of @xmath4 and @xmath5 . when @xmath78 the first - order rw is as shown in fig .
[ fig4](a ) ( see also fig .  [ fig4](d ) for the corresponding contour plot ) . by altering the value of the trap parameter @xmath48 to @xmath79 , the structure of the first - order rw
gets modified as shown in fig .
[ fig4](b)/fig .
[ fig4](e ) .
one can also see from fig .
[ fig4 ] that , as the trap parameter @xmath48 is increased , the rw gradually becomes more localized in time and the condensate atoms settle down to a slightly higher density background due to the attractive nature of the potential . when @xmath70 the modified structure of first - order rw is given in fig .
[ fig4](c)/fig .
[ fig4](f ) , where this feature is even more prominent .
the density profiles of the corresponding second - order rws are presented in fig .
when we tune the parameter @xmath48 from @xmath80 upwards the wave subcrests start to stretch . from the contour plots we can observe that the stretches occur on one side of the rw only .
when we increase the value of @xmath48 further the second - order rw gets modified to a first - order rw like structure which is demonstrated in fig .
[ fig5](c ) .
a similar transition has also been observed in the third - order rw case as well which is illustrated in fig .
the third - order rw acquires a new structure as shown in fig .
[ fig6](b)/fig .
[ fig6](e ) when we increase the value of the parameter @xmath48 from @xmath80 to @xmath79 . at @xmath81 , we observe that the third - order rw acquires a further modified structure which is displayed in fig .
[ fig6](c ) and fig .
[ fig6](f ) .
note the similarity in the central part with that of the first - order rw as given in fig .
[ fig4](c ) and fig .
[ fig4](f ) .      in the case of the third choice
we consider the time - dependent periodic trap frequency to be of the form @xmath47 $ ] so that the strength of the time - dependent periodic interatomic interaction turns out to be @xmath51 .
substituting these two expressions in ( [ a7 ] ) , we find @xmath82 where @xmath83\right\}.\notag\end{aligned}\ ] ] here @xmath57 , @xmath84 , are again the first , second and third - order rw solutions of the nls equation ( vide eqs .
( [ a8 ] ) , ( [ a11 ] ) , and ( [ a12 ] ) ) .        in fig .
[ fig7 ] , we present the density profiles of the first , second , and third - order rws ( first row ) and their corresponding contour plots ( second row ) for the strength of time - dependent interatomic interaction @xmath51 and time - dependent periodic trap frequency @xmath47 $ ] . fixing the value of @xmath48 at @xmath85 , from fig .
[ fig7 ] we see that the rw exists on a periodic background for the above forms of @xmath4 and @xmath86 .
further we have also verified that the behaviour of the rws as a function of @xmath48 follows the same qualitative picture discussed above for the other forms of @xmath4 and @xmath86 .
very recently it has been shown that one can generalize the expressions for the higher order rw solutions of the scalar nls equation ( [ nls ] ) given in appendix a further by introducing certain free parameters which allow one to split a symmetric form rw solution into a multi - peaked solution and that by varying these free parameters one can extract certain novel patterns of rws @xcite .
the introduction of free parameters decompose the higher - order rw solutions into @xmath87 first - order forms , where @xmath88 is the order of the rw , however maintaining the symmetry of the higher - order solutions even in their decomposed forms . further the free parameters are shown to determine the size and orientation of the first - order solutions @xcite .
triplets are symmetry preserving first - order structures revealing the fact that the ( higher- ) second - order rw solution is a family of three first - order rational solutions . it is evident that the existence of triplets in an ocean corresponds to three big waves on the water surface in a row or
three sister waves " @xcite .
inspired by such a possibility , in the following , we consider the second- and third - order rw solutions of ( [ eq:1d - gp ] ) with suitable free parameters and analyze the symmetrical structures that arise due to these free parameters when we vary the strength of nonlinearity and the trap parameter . to begin with ,
we confine our attention to the second - order rw solution . in this case
, we have the following modified expressions @xcite for @xmath89 , @xmath90 and @xmath91 in eq .
( [ a9 ] ) , that is @xmath92 , \notag \\   h_2 = & \ ,   24 \big\ { t\big [ 15 - 16 x^4 + 24 x^2 - 48 l x-8(1 - 4 x^2)t^2   \notag \\           & \ , -16 t^4 \big]+6m(1 - 4 t^2 + 4x^2)\big\ } , \notag   \\
d_2 = & \ ,   64 x^6 + 48 x^4 ( 4t^2 + 1)+12x^2(3 - 4t^2)^2 + 64 t^6\notag \\           & \ , + 432 t^4 + 396 t^2 + 9 + 48m\big[18m+t(9 + 4 t^2 \notag \\          & \ , -12x^2)\big]+48l\big[(18l+x(3 + 12t^2 -4x^2)\big ] .
\end{aligned}\ ] ] note that this rw solution now contains two free parameters , namely @xmath93 and @xmath94 .
the parameters @xmath93 and @xmath94 describe the relative positions of the first - order rws in the triplet . substituting the above expressions in ( [ a15 ] ) , ( [ a22 ] ) and ( [ a24 ] ) , for @xmath95 , we obtain the corresponding second - order rw solutions to the gp equation ( [ eq:1d - gp ] ) with the free parameters @xmath93 and @xmath94 included .        when @xmath96 , this solution coincides with the one given earlier ( vide eq .
( [ a11 ] ) ) which contains one largest crest and four subcrests with two deepest troughs ( fig . [ fig8](h ) ) . when @xmath93 and @xmath94 are not equal to @xmath97 , the second - order rw splits into three first - order rws .
these waves emerge in a triangular fashion ( a triplet pattern ) .
the three first - order rws form a triangular pattern with @xmath98 degrees of angular separation between them @xcite . in figs .
[ fig8b](a)-(c ) we display the triplet pattern for @xmath49 and @xmath45 when @xmath99 and @xmath100 . the triplet rw pattern for @xmath68 is shown in fig .  [ fig8b](a ) .
when we increase the parameter @xmath48 to @xmath101 we observe that the triplet pattern has started to collapse as shown in fig .
[ fig8b](b ) and a complete collapse is observed as in fig .
[ fig8b](c ) when we increase the @xmath48 value further to @xmath79 . figs .
[ fig8b](d)-(g ) represent the triplet pattern for @xmath76 and @xmath77 $ ] . the formation of triplet rws is shown in fig .  [ fig8b](d ) when @xmath68 . when we increase the value @xmath48 to @xmath101 one of the single rws in the triplet pattern vanishes which is illustrated in fig .
[ fig8b](e ) . by increasing
the parameter @xmath48 further we observe that two first - order rws are more localized in time as shown in fig .
[ fig8b](f ) . figs .
[ fig8b](g)-(i ) represent the triplet pattern for @xmath51 and @xmath47 $ ] .
the form of the triplet pattern for @xmath68 is displayed in fig .
[ fig8b](g ) . here
also when we increase the value of @xmath48 we observe the collapse of the triplet pattern in the periodic wave background ( vide fig .
[ fig8b](i ) ) .
we then move on to investigate the structure of the third - order rw solution with four free parameters , namely @xmath102 and @xmath103 .
the third - order rw solution with four free parameters is much lengthier than the one without free parameters and so we do not give the explicit expression here and analyze the results only graphically . here
also we analyze the solution with respect to the free parameters . when @xmath104 , we have the classical third - order rw solution which is shown in fig .  [ fig8](l ) .
it has one largest crest and six subcrests with two deepest troughs . for non - zero values of @xmath93 , @xmath94 , @xmath105 and @xmath103 , the third - order rw splits into six separated first - order rws .
when we increase the value of the free parameters , the six first - order rws take new positions .
the sextet pattern is displayed in figs .
[ fig8h](a)-(c ) for @xmath49 and @xmath45 when @xmath106 , @xmath107 , @xmath108 and @xmath109 .
for @xmath68 the set of six first - order rws is shown in fig .
[ fig8h](a ) .
when we increase the value of @xmath48 to @xmath110 , three peaks disappear and the remaining peaks start to bend as shown in fig .  [
fig8h](b ) .
when we increase the value of @xmath48 further the rws bend in the plane wave background as given in fig .
[ fig8h](c ) . figs .
[ fig8h](d)-(f ) represent six first - order rws for the time - dependent nonlinearity strength @xmath76 and the time - dependent external trap frequency @xmath77 $ ] .
the formation of six first - order rws at @xmath68 is shown in fig .
[ fig8h](d ) .
when we increase the strength of the parameter @xmath48 to @xmath110 , one of six first - order rws vanishes as seen in fig .
[ fig8h](e ) .
if we increase the value @xmath48 further three out of six peaks bend in the plane wave background and eventually collapse which is shown in fig .
[ fig8h](f ) . figs .
[ fig8h](g)-(i ) represent the sextet pattern for @xmath51 and @xmath47 $ ] .
the form of the six first - order rws for @xmath111 is displayed in fig .
[ fig8h](g ) and for further increase in @xmath48 , the modified structures in the periodic wave background are as shown in figs .
[ fig8h](h ) and [ fig8h](i ) .
in the previous two sections we have analyzed how the rw profiles get modified by the variations of the distributed coefficients present in the variable coefficient nls eq .
( [ eq:1d - gp ] ) . in this section
we analyze how the breather structures get modified in the condensates when we vary the strength of the external trap parameter .    to begin with
we consider the first - order breather solution of the nls equation ( [ nls ] ) , which is given in @xcite and is a special case of the gb solution ( [ nls2 ] ) , @xmath112 + 2 i \,f v \sinh[\alpha ( t - t_1)]}{2   \cosh[\alpha ( t - t_1 ) ] - 2v \cos [ f ( x - x_1 ) ] } \notag \\ & \ , -1\biggr ] \times \exp{(it ) } , \label{b1}\end{aligned}\]]where the parameters @xmath113 and @xmath114 are expressed in terms of a complex eigenvalue ( say @xmath115 ) , that is @xmath116 and @xmath117 , and @xmath118 and @xmath119 serve as coordinate shifts from the origin . the real part of the eigenvalue represents the angle that the one - dimensionally localized solutions form with the @xmath120 axis , and the imaginary part characterizes frequency of periodic modulation .
the parameter @xmath121 @xmath122 in ( [ b1 ] ) is the growth rate of modulation instability . substituting this breather solution of the nls equation into ( [ a15 ] ) , ( [ a22 ] ) and ( [ a24 ] ) , we study the underlying dynamics of ( [ eq:1d - gp ] ) .    , ( b ) mb for @xmath123 when @xmath49 and @xmath45 . ( c ) and ( d )
are their corresponding contour plots .
the other parameters are @xmath124 , @xmath125 , @xmath67 , and @xmath69 . ]    for illustration , we consider the case @xmath49 and @xmath45 , and plot the outcome in fig .  [ fig9 ]
the first row in fig .
[ fig9 ] represents an akhmediev breather ( ab ) solution for the eigenvalue @xmath126 and ma breather ( mb ) solution for @xmath123 and the second row represents their contour plots .
the ab and mb solutions are localized in time and space , respectively , as discussed in sec .
i. this is clearly demonstrated in the two columns of fig .
when we tune the parameter @xmath48 , a new structure against a breather background is obtained . to visualize this we fix the value of @xmath48 to be @xmath127 .
for this value a stretching occurs in space in the case of ab ( fig .
[ fig9a](a ) ) whereas in the case of mb ( fig .
[ fig9a](b ) ) the stretching occurs in time
[ fig9a](c)-(d ) illustrate the ab and mb profiles for the time - dependent nonlinearity coefficient @xmath76 and time - dependent trap frequency @xmath128 $ ] . here
also we tune the strength of the trap parameter @xmath48 to 0.8 and observe that stretching occurs over space in ab which is depicted in fig .
[ fig9a](c ) and the mb gets more localized and when @xmath129 the breather profile completely disappears as shown in fig .
[ fig9a](d ) .
our results reveal the fact that when we tune the parameter @xmath48 in the obtained breather solution , the breather gets a modified structure corresponding to a distortion of the breather profile .    ,
( b ) , ( d ) and ( f ) mb with @xmath123 for the three different forms of @xmath4 discussed in the text .
( a)-(d ) @xmath130 , ( e)-(f ) @xmath71 .
the other parameters are @xmath124 , @xmath67 and @xmath69 . ]
figs . [ fig9a](e ) and ( f ) respectively represent the ab and mb in the periodic background for @xmath51 and @xmath47 $ ] when @xmath71 .
next , we proceed to construct the two - breather solutions of ( [ eq:1d - gp ] ) and analyze how these solutions are distorted by the variations of modulation parameters .
the two - breather solution of the nls equation is given by @xcite @xmath131\exp{(it)},\end{aligned}\ ] ] where @xmath132 , @xmath133 , and @xmath134 are given by    @xmath135 ,   \end{aligned}\ ] ]    @xmath136 , \end{aligned}\ ] ]    @xmath137 , \end{aligned}\ ] ]    where the modulation frequencies , @xmath138 , @xmath139 , are described by the ( imaginary ) eigenvalues @xmath140 . in the above expressions , @xmath141 , @xmath142 , @xmath139 , represent the shifted point of origin , @xmath143 @xmath144 is the instability growth rate of each component and @xmath145 and @xmath146 are shifted variables .        with two purely imaginary eigenvalues , @xmath140 , @xmath139 ,
the solution ( [ 2b1 ] ) is capable of describing a variety of possible second - order breather structures .
the solution includes abs , mbs and the intersection of ab and mb solutions for certain combinations of eigenvalues .
for example , when the imaginary parts of both the eigenvalues @xmath147 , @xmath139 , lie between @xmath97 and @xmath148 , we obtain the abs . on the other hand when both of them are greater than one ( @xmath149 ) we obtain the mbs and in the mixed possibility case , that is one of the eigenvalues is less than one ( @xmath150 ) and the other eigenvalue ( @xmath151 ) is greater than one , we obtain the intersection of ab and mb solutions .    substituting the two - breather nls solution ( [ 2b1 ] ) in ( [ a15 ] ) , ( [ a22 ] ) and ( [ a24 ] )
we obtain the general two - breather solution of ( [ eq:1d - gp ] ) .
[ br2 ] displays the evolution of the two - breather solution of ( [ eq:1d - gp ] ) for @xmath49 and @xmath152 with imaginary eigenvalues . to obtain the abs from ( [ 2b1 ] ) we consider the situation where the magnitudes of both the eigenvalues @xmath153 and @xmath154 are less than @xmath148 ( @xmath155 and @xmath156 ) .
one ab developing with a time delay after another is shown in fig .
[ br2](a ) , while in fig .
[ br2](b ) we present the case where there is no such time delay . in fig .
[ br2](c ) we depict the case when one ab along with a rw coexist for the choice @xmath155 and @xmath157 .
when we take the eigenvalues @xmath153 and @xmath154 to be @xmath158 and @xmath159 , respectively , we obtain two mb solutions .
similarly the evolution of the two mbs with and without spatial delay is shown in figs .
[ br2](d ) and [ br2](e ) , respectively .
we also observe that the distance between the mbs increases when we set both the eigenvalues to be nearly equal , say for example @xmath160 and @xmath161 which is not displayed here .
when we take the eigenvalues as @xmath162 and @xmath163 , the ab intersects with the mb which is displayed in fig .
[ br2](f ) .
when we tune the strength of the trap parameter @xmath48 to @xmath164 , both the abs get stretched in the plane wave background which is demonstrated in fig .  [ br2a](a ) . when @xmath165 , in the intersection of ab and mb solutions we observe that the ab gets a bending structure while the mb fully disappears in the plane wave background which is shown in fig .
[ br2a](b ) . in fig .
[ br2a](c ) we note that both the mbs develop a bending structure in the plane wave background when @xmath166 .
figs . [ br2a](d)-(f ) display the evolution of the two - breather solution of ( [ eq:1d - gp ] ) for @xmath76 and @xmath77 $ ] with the imaginary eigenvalues .
when @xmath167 one of the abs gets stretched which is not shown here and on further increase of the value of @xmath48 to @xmath168 , one of the abs gets annihilated and the other ab bends in the plane wave background which is demonstrated in fig .
[ br2a](d ) . when @xmath169 the intersection of ab - mb structures is as shown in fig .
[ br2a](e ) . in fig .
[ br2a](f ) we observe that both the mbs get a bending structure in the plane wave background . figs .
[ br2a](g)-(i ) show the evolution of the two - breather solution of ( [ eq:1d - gp ] ) for @xmath51 and @xmath47 $ ] with the imaginary eigenvalues . here
also when we increase the value of @xmath48 we observe the collapse of the two - breather solution in the periodic wave background .
in this work , we have constructed higher order rw solutions with and without free parameters for the quasi one - dimensional gp equation with time - dependent interatomic interation and external trap through the similarity transformation technique . by mapping the variable coefficient nls equation onto the constant coefficient nls equation
we have derived these solutions .
we have shown that the mapping can be done when the external trap and the nonlinearly interatomic interaction of atoms satisfy a constraint . from the known higher order rw and breather solutions of the constant coefficient nls equation
, we have derived the solutions of ( [ eq:1d - gp ] ) . in our analysis
, we have considered the harmonic trap frequency in three different forms , namely ( i ) time - independent expulsive trap , ( ii ) time - dependent monotonous trap , and ( iii ) time - dependent periodic trap , and correspondingly fixed the effective scattering length .
we then studied the characteristics of the constructed rw solutions in detail .
we have observed that the second- and third - order rws transform to first - order rw - like structures when a parameter appearing in the harmonic trap ( time - independent and time - dependent traps ) is varied .
we have then analyzed the characteristics of triplet and sextet patterns of matter rws for ( [ eq:1d - gp ] ) .
we have also constructed one - breather and two breather solutions of ( [ eq:1d - gp ] ) .
we have investigated how these periodic localized waves change in the plane wave background when we tune the trap parameter in the obtained breather solutions .
our results may provide possibilities to manipulate rws experimentally in a bec system .
km thanks the university grants commission ( ugc - rfsms ) , government of india , for providing a research fellowship . the work of pm forms part of department of science and technology ( ref .
sr / s2/hep-03/2009 ) and council of scientific and industrial research ( ref .
03(1186)/10/emr - ii ) , government of india funded research projects . the work of ms forms part of a research project sponsored by nbhm , government of india .
the work forms part of an irhpa project and a ramanna fellowship project of ml , sponsored by the department of science and technology ( dst ) , government of india , who is also supported by a dae raja ramanna fellowship .      in the absence of the trap , and the nonlinearity strength
@xmath4 is equal to one , eq .
( [ eq:1d - gp ] ) reduces to the standard nls equation @xmath170 .
several localized and periodic structures of standard nls are documented in the literature @xcite .
equation @xmath170 admits @xmath171 order rw solution .
we present the rw solution of the nls equation in the following form @xcite , @xmath172\exp{(it)},\end{aligned}\ ] ] where @xmath173 and @xmath174 are polynomials in the variables @xmath120 and @xmath175 .    for the first - order @xmath176 rw solution @xmath177 , @xmath178 and @xmath179 .
from ( [ a6 ] ) we get @xmath180 .
for covenience we multiply this expression by @xmath181 $ ] and consider the solution in the form @xmath182}.\end{aligned}\ ] ] we use only this form of expression in our analysis .
this solution is same as the one given in eq .
( [ nls3 ] ) when @xmath183 , @xmath184 and @xmath185 . for the second - order @xmath186 rw solution ,
the function @xmath89 , @xmath90 , and @xmath91 turn out to be @xcite @xmath187 and @xmath188 then we get @xmath189\exp{(it)}.\end{aligned}\ ] ] the profiles of second - order rw and third - order rw of the constant coefficient nls equation are shown in fig .
[ nls - rog ] .        for the third - order @xmath190 rw solution
, we have @xmath191\exp{(it)},\end{aligned}\ ] ] where @xmath192 with @xmath193\notag \\ g_6 = &   \frac{2}{45}\left[73 + 14 ( 2 x)^2 + \frac{7}{3}(2 x)^4\right],\notag
\\ g_8 = &   \frac{1}{15 } ( 11 + ( 2 x)^2 ) , \;\;\;\ ;   g_{10 } = \frac{11}{675}\end{aligned}\ ] ] and @xmath194 with @xmath195 , \notag \\ h_2 = &   \frac{2}{3 } \left[-11 - 28 ( 2 x)^2 - 2 ( 2 x)^4 - \frac{28}{45 } ( 2 x)^6 + \frac{(2 x)^8}{45}\right],\notag \\ h_4 = &   \frac{4}{15 } \left[-107 + 19 ( 2 x)^2 - \frac{7}{3 } ( 2 x)^4 + \frac{(2 x)^6}{3 ^ 2}\right],\notag \\
h_6 = &   \frac{4}{45 } \left[-29 - 2 ( 2 x)^2 + \frac{(2 x)^4}{3}\right],\notag   \end{aligned}\ ] ] @xmath196 , \;\;\
; h_{10 } = \frac{2}{675}.\end{aligned}\ ] ] the denominator is represented by @xmath197 where @xmath198,\notag \end{aligned}\ ] ] @xmath199,\notag \end{aligned}\ ] ] @xmath200,\notag \\
d_8 = &   \frac{32}{15 } \left[\frac{83}{3 } + 2 ( 2 x)^2 + \frac{(2 x)^4}{3 ^ 2}\right),\notag \\ d_{10 } = &   \frac{2 ^ 8}{225 } \left[7 + \frac{(2 x)^2}{3}\right ] , \;\;\ ; d_{12 } = \frac{2 ^ 9}{2025}.\end{aligned}\ ] ]


the dynamical response of black holes to any external disturbance is dominated , at intermediate times , by a discrete set of damped oscillations called quasinormal modes ( qnms ) .
this result has been verified not only at linearized level , but also in fully numerical simulations for stellar gravitational collapse @xcite and collision of two black holes @xcite .
the possibility of observation in the forthcoming years of the radiation due to the qnms opens a new perspective for the gravitational wave astronomy . from these observations it will be possible to obtain direct evidence of black holes existence , as well as to estimate their parameters ( mass , electric charge , and angular momentum ) @xcite .
in addition to the astrophysical applications , an interpretation appeared some years ago of the qnms in terms of the ads / cft conjecture @xcite . according to this conjecture
, there is a correspondence relation ( or duality ) between string theory in an anti - de sitter ( ads ) spacetime and a conformal field theory ( cft ) on the boundary of this space @xcite .
it implies that a large black hole in the ads spacetime corresponds to an approximately thermal state in the cft .
perturbing the black hole is equivalent to perturbing the thermal state , and the perturbation damping time gives the time scale for the return to the thermal equilibrium .
therefore , by computing quasinormal ( qn ) frequencies in an ads space it is possible to obtain the thermalization time scale in the strongly coupled cft .
the early works studying perturbations of ads black holes considered the evolution of a conformally invariant scalar field in such background spacetimes @xcite .
after the advent of the ads / cft conjecture , qnms of spacetimes with a negative cosmological constant , such as kottler ( also called schwarzschild - ads ) @xcite , reissner - nordstrm - ads @xcite and kerr - ads @xcite , were computed for different fields , dimensions , and boundary conditions . within this program , other kinds of black hole solutions were also analyzed .
for instance , in the baados - teitelboim - zanelli solution @xcite , analytical expressions for the qn frequencies presented exact agreement with the poles of the retarded correlation function in the dual cft @xcite .
this result provided an important quantitative test of the ads / cft correspondence . as a complement to the aforementioned studies , analyses of the time evolution of some fields in asymptotically ads backgrounds
have been performed @xcite .
more recently , the conjecture relating the highly damped qnms to black hole area quantization has motivated a search for qn frequencies with large imaginary parts also in asymptotically ads spacetimes @xcite .
the studies of perturbations in ads spacetimes were also extended to black holes with nonspherical topology . in ref .
@xcite the evolution of a conformal scalar field in the family of topological backgrounds was investigated , which includes the plane - symmetric black hole as one of its members , while in ref .
@xcite the case of a nonminimally coupled scalar field in the same kind of spacetimes was considered . lately , the lowest qn frequencies associated with the scalar , electromagnetic , and gravitational perturbations of the plane - symmetric black hole background were computed @xcite .
the study of the scalar and electromagnetic fields was relatively complete .
however , important questions were left without answers in connection with the gravitational perturbations .
for instance , what are the wave equations governing general metric perturbations ? what is the physical interpretation of the modes with different wave numbers ?
what are the qn frequencies in the intermediate and small black hole regimes ? what is the relation between these frequencies and the wave number ?
why is the pure imaginary frequency not proportional to the horizon radius in the large black hole limit as it happens for ordinary modes ?
is there some relation between the latter and the algebraically special frequency ? and what is the behavior of the qnm spectrum in the regime of highly damped overtones ?
the objective of the present work is answering most of the previous questions as well as other ones that appeared during this work .
some of the conclusions drawn here are extended to the schwarzschild - ads black hole .
the plan of this paper is the following . in sec .
[ solucao ] , we write the equations defining the background spacetime in a newman - penrose ( np ) null tetrad frame @xcite . in sec .
[ perturb - np ] , we obtain schrdingerlike equations governing the axial and polar perturbations . in sec .
[ qnm ] , we define the qnm boundary conditions to be used and prove the black holes stability against perturbations with a non - negative potential .
a brief description of the horowitz - hubeny method ( sec . [ metodo ] ) follows , as well as subsequent exposition of the qn frequency values for different horizon radii and wave numbers ( sec .
[ resultados ] ) .
then , we study the large black hole limit of the axial qnms ( sec .
[ limite ] ) and the non - qn character of the algebraically special frequency ( sec .
[ algebrico ] ) . finally , we conclude by summarizing and commenting the main results ( sec .
[ conclusao ] ) .    throughout this work , we use the geometric system of units in which the speed of light and the newton s gravitational constant are set to unity .
the background spacetime , moreover , has a natural lengthscale given by the ads radius @xmath5 , where @xmath6 is the negative cosmological constant .
here we choose the length unit so that @xmath7 and measure everything in terms of the ads radius . in relation to notation , the asterisk ( @xmath8 ) denotes complex conjugation and the definition of the np quantities follows ref .
with the aim of investigating the evolution of linear gravitational perturbations in the background of a plane - symmetric ads black hole @xcite , we summarize here the basic properties of such a spacetime .
the metric may be written in the form @xmath9 where @xmath10 there is an event horizon at @xmath11 , and an essential singularity at @xmath12 .
depending on the ranges of coordinates @xmath13 and @xmath14 , the metric may represent a cylindrical ( @xmath15 , @xmath16 ) , toroidal ( @xmath17 ) or planar ( @xmath18 ) black hole spacetime .
for the torus , @xmath19 is related to the adm system mass @xcite ; for the cylinder , to the mass per length unit of @xmath14 line @xcite ; and for the plane , to the mass per area unit of the @xmath20 plane @xcite .    according to the petrov classification , the above metric is algebraically type d. it is then convenient to define the np quantities in terms of a kinnersley - like null frame @xcite . in this case
, the two real vectors of such a frame , @xmath21 and @xmath22 , are chosen to lie , respectively , on the ingoing and outgoing radial null geodesics and are defined by @xmath23 which are the double principal null directions of the weyl tensor . in order to complete the frame basis we have the vector @xmath24 and its complex conjugate @xmath25 .
both of these vectors , @xmath26 and @xmath27 , are orthogonal to @xmath28 and @xmath29 , and satisfy the normalization condition @xmath30 . with the chosen tetrad basis ,
the only nonvanishing spin coefficients are @xmath31 and the only nonvanishing weyl and ricci scalars are @xmath32
here the equations for the linear gravitational perturbations are obtained by using the np formalism , which allows us to investigate general properties of such perturbations without any further simplifying assumption . in this approach , changes in the metric coefficients are directly related to changes in the tetrad null vectors .
furthermore , perturbations in the vectors @xmath28 , @xmath29 and @xmath33 lead to first order alterations in weyl and ricci scalars , and also in the spin coefficients .
the full system of linearized perturbation equations is obtained by expanding the complete np equations on the plane - symmetric black hole background , and retaining only the first order terms in the perturbation functions .
generally , the linearized np equations form a large system of coupled equations whose solutions are difficult to be found . however , for petrov type d vacuum spacetimes , the problem reduces to solving a pair of decoupled equations for the gauge and tetrad invariants @xmath34 and @xmath35 @xcite . in the present case ,
these equations assume the simple forms @xmath36\psi_{0}=0 , \label{fund1}\\ [ ( \delta+3\gamma-\gamma^{\ast}+4\mu+\mu^{\ast})(d-\rho)- \delta^{\ast}\delta-3\psi_{2}]\psi_{4}=0 . \label{fund2}\end{gathered}\ ] ] the foregoing equations are already linearized since the weyl scalars @xmath34 and @xmath35 vanish in the equilibrium state , and then they are to be considered first order perturbations .
the other quantities take their unperturbed values given in eqs . - .
the background symmetries , together with the perturbation equations and , indicate that each one of the perturbation functions @xmath34 and @xmath35 can be written as a product of four different functions depending upon only one of the coordinates @xmath37 , @xmath38 , @xmath13 , and @xmath14 .
moreover , due to the wavelike character of the equations , the dependence of these functions on the coordinates @xmath37 , @xmath13 and @xmath14 is of the form @xmath39 $ ] , where @xmath40 is the frequency , and @xmath41 and @xmath42 are wave numbers which may assume discrete or continuous values , depending on whether @xmath13 and @xmath14 are compact or noncompact coordinates , respectively .
the radial functions , in turn , are conveniently written as @xmath43 where the subscripts @xmath44 are related to the spin weight and to the conformal weight of the weyl scalars @xmath34 and @xmath35 .
the resulting equations for @xmath45 take the following standard form : @xmath46 where we introduced the differential operators @xmath47 here @xmath48 is the tortoise coordinate , defined in such a way that @xmath49 , and @xmath50 with @xmath51 .    the chandrasekhar transformation theory @xcite can now be used to transform eqs . into one - dimensional schrdingerlike equations of the form @xmath52 this is done by means of the substitutions @xmath53 whose inverse transformation gets the form @xmath54 where @xmath55 in the foregoing equations , the potential functions @xmath56 are given by @xmath57 , \label{pot - polar}\ ] ] @xmath58 while the auxiliary functions @xmath59 are , respectively , @xmath60    the meaning of the signs @xmath61 comes from the metric approach to the perturbation problem , as it was performed in ref .
they are related to the parity property of the perturbed metric functions under the exchange @xmath62 ( or @xmath63 ) .
the variable @xmath64 represents the polar ( even ) perturbations , while @xmath65 represents the axial ( odd ) ones .    in a certain sense
, the wave equations can be thought of as a generalization to those found by cardoso and lemos @xcite . here
the perturbations depend on both coordinates @xmath13 and @xmath14 , while in the aforementioned work the analysis was restricted to @xmath14-independent metric variations .
this is equivalent to the vanishing of the wave number @xmath42 , resulting in @xmath66 .
however , the plane symmetry of the background spacetime renders possible to find new coordinates ( @xmath67 , @xmath68 ) in such a way that one of the wave numbers , @xmath69 or @xmath70 , is zero .
that is to say , it is always possible to find a coordinate system in which the perturbations depend only on two space coordinates ( @xmath38 and @xmath13 , or @xmath38 and @xmath14 ) .
another result derived from the above analysis is related to the starobinsky constant @xcite , @xmath71 , where @xmath72 are given by eq . .
the solutions of @xmath73 correspond to algebraically special frequencies in the sense of chandrasekhar @xcite . as we shall see in sec .
[ resultados ] , one of such solutions , given by @xmath74 appears in the numerical results in connection to the purely damped qnms .
the properties of the axial perturbations related to these algebraically special modes are analyzed in some detail in sec .
[ algebrico ] .
to solve the wave equations we must impose boundary conditions on the possible solutions .
first , however , consider the case of the asymptotically flat schwarzschild spacetime .
gravitational perturbations lead to schrdingerlike equations whose potentials vanish at infinity and at the event horizon .
hence , close to these extremities , the solutions are plane waves , @xmath75 , where the @xmath48 coordinate ranges from @xmath76 to @xmath77 .
qnms are then defined as solutions which are purely ingoing waves at the horizon , @xmath78 , and purely outgoing waves at infinity , @xmath79 .
these requirements are satisfied only for a discrete set of complex @xmath40 called the quasinormal frequencies .    in asymptotically ads spacetimes
, we keep the requirement of purely ingoing waves at the horizon . in the meantime , since now the tortoise coordinate @xmath48 takes a finite value at spatial infinity , the boundary condition at infinity must be changed .
this can be dirichlet , neumann , or robin boundary condition , depending on whether the field , its derivative , or a combination of both vanishes , respectively . here , as in most of the literature , we adopt dirichlet boundary conditions , since these lead to energy conservation in the field theory on ads space @xcite ( for alternative boundary conditions see ref .
@xcite ) .
having defined the qnms , we show below that the imaginary part of the qn frequencies is positive for non - negative potentials .
the demonstration follows the one carried out in the schwarzschild - ads case by cardoso and lemos @xcite .
initially , we exchange the variable @xmath80 by @xmath81 , so that eq . becomes @xmath82 where @xmath80 , @xmath83 and @xmath84 stand for both perturbation types , the polar @xmath85 and the axial ones @xmath86 .
near the horizon and at infinity , the solutions @xmath83 satisfying the qnm boundary conditions are given , respectively , by @xmath87 } , \label{assint1}\ ] ] @xmath88 where @xmath89 and @xmath90 are arbitrary constants , and @xmath91 and @xmath92 are definite constants whose values depend on the perturbation type .
now we multiply eq . by @xmath93 and integrate by parts using the asymptotic solutions and . after some algebra
we obtain @xmath94dr= \frac{|\omega|^{2}|\phi(r_{+})|^{2}}{\mbox{im}(\omega)}. \label{int - final}\ ] ] since @xmath95 assumes only positive values outside of the black hole , the sign of the quantity on the left hand side of eq .
depends crucially on the sign of the potential @xmath96 . for non - negative @xmath84
, it follows that the imaginary part of the qn frequencies is positive , @xmath97 .
moreover , it is seen from eq . that @xmath98 is positive in the region @xmath99 for any @xmath1 and @xmath0 , and
then we conclude that the plane - symmetric ads black holes are stable against polar perturbations .
the same is true for axial perturbations as far as the wave numbers are restricted by @xmath100 . on the other hand , for @xmath101 ,
the potential @xmath102 takes negative values and the theorem can not be applied .
another interesting property of frequencies associated to perturbations of a plane - symmetric black hole comes from the particular symmetry of such a background .
this symmetry allows us to rescale @xmath40 , @xmath1 and @xmath0 by a pure scale transformation of coordinates @xmath103 , @xmath104 , @xmath105 , and @xmath106 , where @xmath107 is a constant @xcite .
it implies that the frequency is a first degree homogeneous function of its parameters : @xmath108 .
in particular , we can take @xmath109 so that @xmath110 then , after computing the qn frequencies for a given @xmath1 , e. g. , for @xmath111 , and for different values of @xmath0 , we can obtain @xmath112 by using the previous equation , and the values of @xmath40 for a fixed wave number @xmath1 are naturally divided into three different regimes : large ( @xmath113 ) , intermediate ( @xmath114 ) and small ( @xmath115 ) black holes . in the large ( small ) black hole limit , @xmath1 ( @xmath0 ) is negligibly small and it is expected that @xmath116 ( @xmath117 ) .
we are now in a position to solve the perturbation equations and to find the values of the qn frequencies .
this task , however , is more easily done through numerical procedures .
therefore , it is appropriate summarizing here the method used for such a purpose .
in the last decades a series of methods has been developed to numerically calculate the qnms of black holes . here
we briefly describe one of these methods , used in ref .
@xcite , which is of particular interest for the present work since it is suited for asymptotically ads spacetimes . in order to apply such a method we rewrite here the perturbation equations in the appropriate form .    by introducing the independent variable @xmath118 , eq . becomes @xmath119 where @xmath120 and @xmath121 , with @xmath84 and @xmath122 now being considered as functions of @xmath123 . the foregoing equation has a regular singular point ( @xmath124 ) in the region of interest , @xmath125 .
then , by fuchs theorem , any solution of can be written as a frbenius series @xcite .
the behavior of the solutions near the horizon may be explored by writing @xmath126 , where @xmath127 is a constant to be determined . by substituting this expression into eq . and keeping terms only to the leading order , it follows that @xmath128 , which has the roots @xmath129 and @xmath130 .
these roots correspond exactly to the ingoing and outgoing waves at the horizon , respectively @xcite . since we want only ingoing modes , we take @xmath129 and look for solutions of the form @xmath131 substituting eq . into ,
we find expressions for the @xmath132 s in terms of the frequency @xmath40 ( and of the parameters @xmath1 and @xmath133 ) .
then , imposing the boundary condition at infinity , @xmath134 , we obtain @xmath135 thus , the problem of finding qn frequencies has been reduced to that of obtaining the roots of the infinite - order polynomial equation .
of course , we can not carry out numerically the full sum in expression .
in effect we truncate the series after a large number of terms and look for the zeros of this partial sum .
the precision of the results is then verified by computing the relative variation between the zeros as we perform higher partial sums .
we stop our search once we have a @xmath136 decimal digit precision . for each pair of @xmath1 and @xmath0
, we find a set of qn frequencies ( out of the infinity of possible frequencies ) .
these frequencies are labeled as usual with the principal quantum number @xmath137 and ordered in a set beginning with the roots with the lowest imaginary parts .
we use the horowitz - hubeny procedure described above to compute the gravitational qnms of plane - symmetric ads black holes . as it happens to other ads black holes , there are _ purely damped modes _ ,
whose frequencies are pure imaginary numbers , in addition to the _ ordinary modes _ , the frequencies of which also have nonzero real parts which are responsible for the oscillatory behavior of the perturbations .
as mentioned earlier , the numerical analysis is divided into three different regimes according to the ratio @xmath138 values . that is to say , we have ( i ) large black holes when the ratio @xmath138 is large compared to unity ; ( ii ) intermediate black holes if @xmath138 is of the order of unity ; and ( iii ) small black holes when @xmath138 is small compared to unity .
it is also worth mentioning here that some numerical analysis of the gravitational qnms was done in ref .
however , such a study was restricted mainly to the large black hole regime for a specific value of @xmath1 , without generalizations for a generic wave number . in
what follows , we perform a more detailed analysis , including small and intermediate black holes , and extend implicitly the numerical results to all values of @xmath1 .
we also search the qn frequency values in the limit of highly damped overtones .      to begin with
, we discuss the case of regular ( ordinary ) qnms for which the frequencies have nonvanishing real and imaginary parts , @xmath139 .
we list some numerical results of the lowest quasinormal frequencies for @xmath140 and some selected values of @xmath0 .
the computed quantities include both the axial and polar perturbations , and the three regimes defined above .
the data for large , intermediate , and small black holes are shown , respectively , in tables [ tlarge ] , [ tmedium ] , and [ tsmall ] .
.the lowest ordinary qn frequencies for @xmath140 and several values of @xmath0 in the large black hole regime . [ cols="^,^,^,^,^,^ " , ]
although numerical procedures need to be used to compute the qn frequencies , it is possible to study analytically the special case of axial perturbations of very large black holes .
here we show in some detail the process of finding an asymptotic expression for the frequency of the purely damped modes .
the analysis is performed in such a way that it holds for both the plane - symmetric and the schwarzschild - ads spacetimes .
our starting point is the horowitz - hubeny approach described in sec .
[ metodo ] . there
, the qn frequencies of axial perturbations are associated to solutions of the equation @xmath141 hence , we need to obtain an analytical expression for the sum of this series in order to find exact frequency values . such a task is very difficult to be carried out for a general @xmath133 , but it becomes feasible for the large black hole limit of purely damped modes .
the coefficients @xmath142 appearing in eq .
are the frbenius coefficients of the power series solution to either eq . written for the axial functions , or the corresponding equation in the schwarzschild - ads case ( see ref .
they are given in terms of the frequency , and of the other physical parameters of the problem , by the following recurrence formulas : @xmath143 .
\label{rec2}\end{gathered}\ ] ] the last equation is valid for all @xmath144 and we have shortened the equations by introducing the quantities @xmath145x_{+}^{2 } , \label{def1}\\ b_{m}=(m^{2}-2m-3)\left[1+\kappa x_{+}^{2}\right ] , \label{def2}\\ f_{m}=\left[m(3+\kappa x_{+}^{2})+2i\omega x_{+}\right]^{-1 } , \label{def3}\end{gathered}\ ] ] where @xmath146 and @xmath147 depending on whether it is considered a plane or spherical hole , respectively , and now the foregoing expressions hold for any @xmath148 .    in the large black hole regime , the wave number ( or the angular momentum ) of the perturbation is negligible and the frequency scales as @xmath149 . in terms of the variable @xmath133 , it implies that @xmath40 has a first order pole at @xmath150 .
therefore , @xmath151 has a laurent series about @xmath150 , which is valid over the range @xmath152 , for some @xmath153 @xcite .
such an expansion can be conveniently written as @xmath154 where @xmath155 are constants to be determined .
however , not all of the elements @xmath155 are independent , since part of them can be eliminated on the basis of the parity property of the function @xmath151 .
independently of what results from the sum on the left - hand side of eq . , the form of the quantities @xmath156 , @xmath157 , and @xmath158 ensures that the product @xmath159 is a function of @xmath160 .
this means that @xmath40 changes its sign under the transformation @xmath161 , being an odd function of @xmath162 .
accordingly , the series representing @xmath163 should be an even function of @xmath133 , which means that only the even powers of @xmath162 should survive , and then one should have @xmath164 for all @xmath165 .
this result allows us to write @xmath166 where @xmath167 and @xmath168 are constants to be determined .
the series representation of @xmath169 opens the possibility of taylor expanding the terms @xmath170 about @xmath171 . before doing this , we notice that the value of @xmath172 is arbitrary , since all the terms in the series of eq .
are proportional to this element and the full sum is equal to zero . for simplicity ,
we choose here @xmath173 .
then we substitute eq . into eq . and expand @xmath158 ( which depends on @xmath162 and also on @xmath40 ) about @xmath171 . when this expansion is substituted back into eqs . and
, we obtain a set of taylor series for the elements @xmath170 ( one of such series for each @xmath174 ) , which may be represented by @xmath175 where the elements @xmath176 are functions of @xmath177 , @xmath178 , @xmath179 , @xmath180 , ... , @xmath181 . substituting the previous expansions into eq . and collecting terms of the same power in @xmath14 , we find @xmath182 + \left[\eta_{1}^{1}(c_{0},c_{1})+\eta_{2}^{1}(c_{0},c_{1 } ) + ... \right]z+ ... =0 .
\label{colecao}\ ] ]    yet the left - hand side of eq . is a power series in the variable @xmath14 .
it implies that , in order for the equation being satisfied , all bracketed terms should vanish .
therefore , the zeroth order coefficient provides an equation for @xmath179 , the first order coefficient , together with the earlier calculated value of @xmath179 , provides an equation for @xmath180 , and so on . here
we begin with the computations that lead to @xmath179 .
the first two relevant terms for this task are @xmath183 and @xmath184}$ ] .
moreover , it can be seen from eq . that @xmath185 is identically zero , and as a consequence the term @xmath186 differs from @xmath187 only by a multiplicative factor .
this implies , in particular , that the element @xmath188 is related to @xmath189 by @xmath190 such a relation , together with the formulas and , implies that all the subsequent elements are also proportional to @xmath189 .
more specifically , the @xmath174th coefficient ( for @xmath191 ) can be written as @xmath192 , where @xmath193 is given by the recurrence formula @xmath194,\ ] ] with @xmath195 and @xmath196 .
taking into account the last results , the @xmath14-independent term of eq .
gives rise to the following equation for @xmath179 : @xmath197=0 .
\label{puredampedzeros}\ ] ]    clearly the above equation has the special solution @xmath198 , in addition to the roots that arise as zeros of the square - bracketed term .
the existence of this special solution explains why the purely damped modes do not scale with the horizon size , as it happens to ordinary modes in the large black hole limit .
in fact , the numerical computations show that purely damped modes are the only qnms which in first approximation are not proportional to @xmath0 .
this fact suggests they should be associated to the root @xmath198 , and then it follows from eq .
that for large horizon radii the corresponding frequencies go as @xmath199 . on the other hand , the numerical results also indicate that the ordinary axial qnms are proportional to @xmath0 ( in first approximation ) . from this one
concludes that the remaining roots that are obtained from the zeros of the term among brackets in eq .
correspond to the asymptotic values of the function @xmath200 for these ordinary qnms .
it is seen from the above analysis that , in order to calculate the leading term of the large horizon limit for the purely damped mode frequency , we must take @xmath198 and calculate the next parameter @xmath180 . in such a case , the coefficients @xmath201 and @xmath202 reduce to @xmath203/3 $ ] and @xmath204 , while the other elements @xmath205 ( for @xmath206 ) assume the form @xmath207 . substituting these results into the first order coefficient of we
obtain an equation that can be solved for @xmath180 , yielding @xmath208}{\left(4+\gamma \right)},\ ] ] where we have defined @xmath209    the problem now is summing the series that leads to @xmath210 .
fortunately , we succeeded in summing up this series ( see appendix b ) and the exact result is @xmath211 . with this
, the coefficient @xmath180 becomes @xmath212 where @xmath213 and @xmath214 are given by @xmath215    once we have obtained an expression for @xmath180 , the asymptotic formula of the pure imaginary qn frequency follows directly from eq . .
in fact , considering a plane - symmetric ads black hole , for which @xmath216 and @xmath217 , we have @xmath218 , and the interesting result follows : @xmath219 alternatively , for a schwarzschild - ads black hole one has @xmath220 and @xmath221 , which implies in @xmath222 and @xmath223 , leading to @xmath224 both of the last two formulas for the purely damped frequencies are exact results at the asymptotic limit of large black holes , in perfect agreement with the approximate fits found through the numerical calculations .
the numerical results of the qn frequencies revealed a close relation between the algebraically special modes and the purely damped qnms .
for a fixed wave number , the imaginary frequency @xmath225 approaches the values of the special frequency @xmath226 as the horizon radius goes to zero ( see the comments at the end of section [ perturb - np ] ) .
there also seems to exist a meeting point @xmath227 between the two frequencies for a special combination of parameters such that @xmath228 , as it can be seen from the numerical results shown in fig .
[ fpuro ] .
furthermore , the relation between the purely damped modes and the algebraically special frequency in the small horizon limit is not restricted to the plane - symmetric black hole case , but it also appears in axial perturbations in the schwarzschild - ads spacetime .
all these facts render it convenient to analyze here the qn character of the algebraically special perturbation modes .
the frequency function of these modes for the plane - symmetric black hole is given by eq .
, while the corresponding expression in the spherical case is found in ref .
@xcite . in a general form
, the frequency @xmath226 can be written as @xmath229 } { 6(1+\kappa x_{+}^{2})}x_{+}^{3},\ ] ] where again @xmath146 and @xmath147 , as it is considered a plane or spherical black hole , respectively .    to verify
if the algebraically special modes satisfy the qnm boundary conditions we need to obtain the solution @xmath230 associated with the algebraically special frequency @xmath226 . since
this function is a particular case of the general power series solution , we may use eqs . and to determine the coefficients @xmath132 for @xmath231 .
once again we may take @xmath173 , so that the coefficients @xmath232 and @xmath233 are given , respectively , by @xmath234 @xmath235 + 6i\omega ( 1+\kappa x_{+}^{2})}{4x_{+}(3+\kappa x_{+}^{2}+i\omega x_{+ } ) ( 3+\kappa x_{+}^{2}+2i\omega x_{+})}.\ ] ] it is seen from the last equation that @xmath233 vanishes for the frequency @xmath231 .
moreover , since all the other elements @xmath132 ( for @xmath206 ) differ from @xmath233 only by a multiplicative factor , we have @xmath236 .
the solution for @xmath230 then reduces to @xmath237x_{+}^{2}}\left(1-\frac{x}{x_{+}}\right ) .
\label{afunc}\ ] ]    we know from the horowitz - hubeny method that the function represents an ingoing wave at the event horizon . here
the problem is verifying if @xmath230 satisfies the dirichlet boundary condition at the spatial infinity . imposing the condition @xmath238
, we find an equation that is satisfied only if @xmath239 or @xmath240 , where @xmath214 and @xmath213 are given by eq . .
for a plane - symmetric ads black hole , this means we need a vanishing wave number in order to have a qn frequency equal to @xmath226 . in this case , however , the resulting perturbation functions would be independent of time , which could not represent gravitational waves .
similarly , for a schwarzschild - ads spacetime we need to have @xmath241 or @xmath242 .
the first root is a physically forbidden value , while the second one leads to a static perturbation , which only adds a small rotation to the black hole @xcite .
on the basis of these results , we conclude that @xmath226 can not be a qn frequency , because the associated wave function @xmath230 does not satisfy the dirichlet boundary condition at infinity .
although the numerical calculations show that in the small horizon limit @xmath225 is very close to @xmath226 , the above result proves that the frequency function associated with purely damped modes is not exactly equal to the algebraically special frequency . the solid line in fig .
[ fpuro ] does not coincide with the dots , which represent qnms , not even in the small black hole regime
. there is , however , a special point @xmath227 where @xmath226 seems to coincide with @xmath243 .
the answer to the question whether @xmath226 corresponds to a qn frequency at that point or not is given in the following .    in spite of the generality of the analysis performed above , it fails to test the qn character of @xmath226 in the case where @xmath244 .
this situation happens for a special family of parameters satisfying the condition @xmath245 in this case the roots of the frbenius index equation , which are @xmath129 and @xmath246 ( @xmath247 ) , differ by an integer number and the horowitz - hubeny method is not valid .
we can circumvent this problem by writing explicitly the axial version of eq . for @xmath231 and parameters satisfying eq . .
by doing this in such a way that the analysis holds for both the plane and the spherical ads black holes , it is found that @xmath248 \frac{d^{\,2}\phi_{a}^ { ( -)}}{dx^{2}}+ \left[3(1+\kappa x_{+}^{2})x + ( 3+\kappa x_{+}^{2})x_{+}\right ] \frac{d\phi_{a}^ { ( -)}}{dx}-3(1+\kappa x_{+}^{2 } ) \phi_{a}^ { ( -)}=0 .
\label{fund - new}\ ] ] here @xmath249 is a regular point of the above differential equation , and the solution for @xmath250 can be expressed as a taylor series , @xmath251 substituting this expression into eq . , we find as usual that all coefficients @xmath252 ( for @xmath253 ) can be expressed in terms of the first two , @xmath254 and @xmath255 .
the general solution of eq . then assumes the form @xmath256\\ & + d_{1}\left[(x - x_{+})+ \frac{(3 + 2\kappa x_{+}^{2})}{(3+\kappa x_{+}^{2 } ) } \sum_{m=2}^{\infty}\frac{\beta_{m}}{(-x_{+})^{m-1 } } ( x - x_{+})^{m}\right ] , \label{sol - geral } \end{split}\ ] ] where the parameters @xmath257 ( for all @xmath258 ) are given by the recurrence formula @xmath259 with @xmath260 and @xmath261 .    the first term among brackets in eq .
represents a perturbation function which is constant at the event horizon and corresponds to an ingoing wave , while the second term represents an outgoing wave , since it goes as @xmath262 in the limit @xmath263 . to satisfy the qnm boundary condition at the horizon , we must take @xmath264 and keep only the first term of the general solution .
in addition , without loss of generality we may choose @xmath265 .
finally , we must verify if the function @xmath230 satisfies the dirichlet boundary condition at the spatial infinity @xmath266 . by substituting @xmath267 , @xmath268 , and @xmath266 into eq .
it follows that @xmath269    for plane - symmetric ads black holes , our numerical calculations show that @xmath270 , independently of the horizon radius . on the other hand ,
we obtain different values for @xmath271 as we change the schwarzschild black hole size .
another way to see this variation is in terms of the angular momentum , since @xmath133 and @xmath4 are coupled through eq . .
[ falgeb ] , we plot the numerical values of @xmath272 ( dashed line ) and of @xmath271 for @xmath273 ( dots ) .
as the angular momentum goes to infinity , which means @xmath274 , the results found for the spherical case tend to the value of the plane black hole .
since @xmath275 is always different from zero , we conclude that the algebraically special modes are not qnms , not even for the special parameters that lead to @xmath228 . for the plane - symmetric black hole and for @xmath140 ,
this corresponds to the point @xmath227 in fig .
[ fpuro ] , where the lines representing the purely damped qn and the algebraically special frequencies tend to intercept each other . among other things
, this result implies that the qnm spectrum of some black holes does not present a particular pure imaginary frequency , for those specific modes whose wave numbers ( or angular momentum ) satisfy eq . .
for instance , an axial quadrupole ( @xmath276 ) perturbation is not able to excite a purely damped qnm in a schwarzschild - ads black hole with @xmath277 .
whether this conclusion has some physical implication to the ads / cft conjecture or not is still an open question .
we have studied the gravitational quasinormal modes of plane - symmetric anti - de sitter black holes .
these modes are responsible for the late time behavior of small gravitational disturbances .
the wave equations governing the axial and polar perturbations were derived by a procedure based on the newman - penrose formalism and on the chandrasekhar transformation theory .
this approach allowed us to treat equally the two spacelike coordinates @xmath13 and @xmath14 , as well as provided the algebraically special frequencies associated to the gravitational waves radiated by the black hole .    using the horowitz - hubeny method ,
we have computed the quasinormal frequencies for @xmath140 and several values of @xmath0 in the large , intermediate , and small black hole regimes .
the wave - number dependence of the modes was obtained on the basis of the homogeneity property of the function @xmath112 . for large black holes ,
the ordinary frequencies scale linearly with @xmath0 , and the pure imaginary frequency goes as @xmath199 . for small black holes , the fundamental modes present @xmath278 and @xmath279 , while @xmath225 approaches the algebraically special frequency @xmath280 . in the limit of highly damped overtones
, we found the qn frequencies are evenly spaced , with the spacing per unit of horizon radius being a constant which is independent of the wave - number value , perturbation parity , and black hole size .    in relation to stability , we have shown analytically that for a non - negative potential the imaginary part of the frequency is positive definite for every mode . from the numerical results it also follows that all qn frequencies have a positive imaginary part , including the cases where the potential takes negative values .
these results conjointly imply that plane - symmetric ads black holes are stable against gravitational perturbations .    to study the large black hole limit we have developed an analytical method based on the expansion of the frequency in powers of @xmath199 .
application of this method to axial perturbations showed that there is a special frequency for which the ratio @xmath281 goes to zero as @xmath282 .
the first nonvanishing term is @xmath283 for the plane hole and @xmath284 for the schwarzschild hole .
it is important to observe that these purely damped modes are particularly long living for large black holes , which is exactly the relevant regime for the ads / cft correspondence .
if the above asymptotic formulas could be derived on the cft side , it could provide an important quantitative test to the maldacena s conjecture .    as a last remark
, it should be noted that the discussion of the physical interpretation of the modes with different wave - number values remains .
preliminary results show that the zero wave - number axial perturbations yield only small rotations on the system , while the polar ones lead to an increase in the mass and the production of gravitational waves @xcite .
this issue , however , must be investigated in more detail and shall be considered in our future work .
we are indebted to v. cardoso for suggestions and for a critical reading of the manuscript .
we also thank j. p. s. lemos for stimulating conversations .
one of us ( a. s. m. ) thanks the conselho nacional de desenvolvimento cientfico e tecnolgico - cnpq , brazil , for a grant .
in this appendix we discuss briefly how to find the modes of gravitational perturbations in pure plane - symmetric ads spacetimes .
this is to be compared to the small black hole limit qnms , as calculated numerically in sec .
[ resultados ] .    the relevant wave equations for axial and polar perturbations may be obtained from eq . by assuming @xmath286 wherever it appears .
for instance , the potentials @xmath287 reduce to the same constant value @xmath288 the resulting equations are then @xmath289 whose solutions in terms of the coordinate @xmath290 are of the form @xmath291 where @xmath292 , @xmath293 , @xmath294 , and @xmath295 are integration constants .
the solution functions for @xmath296 are analytic in the whole range of the radial coordinate which is of interest for comparison to the small black hole qnms , @xmath297 .
hence , the only condition @xmath40 must satisfy is @xmath298 in order to avoid the wave growing exponentially in the limit @xmath299 .
however , additional restrictions to these modes arise by imposing the usual boundary conditions for qnms in ads spacetimes : no wave can come out from the black hole horizon , which is located at @xmath300 ; and no wave can propagate to spatial infinity , @xmath301 . by substituting these conditions into eqs .
it follows that @xmath302 and @xmath303 , which means waves whose frequencies satisfy the relation @xmath296 do vanish everywhere .
we are then left with the solutions for which @xmath304 only . in this case
there is in fact no wave in the radial direction , so that the solutions @xmath305 , for @xmath306 , represent waves traveling along the spatial directions orthogonal to the @xmath38 direction . in other words , there are no waves coming out from the horizon nor propagating to spacelike infinity .
considering these as the plane - symmetric pure ads modes , their frequencies , given by @xmath307 agree with the result found numerically for small black holes .
here we show how the series presented in eq . can be summed up exactly . since it is an alternating series , it is important not only to prove its convergence , but also to test it for absolute convergence . by mathematical induction , we find that the series can be written as a combination of seven terms @xmath308 where the elements @xmath309 ( @xmath310 , @xmath311 , ... , @xmath136 )
are given by @xmath312 then we are left with a sum of six series , and by proving the absolute convergence of each one of them , the absolute convergence of the whole series @xmath313 is also proved . in order to do that
, we calculate the limit @xmath314 and verify that it is smaller than unity for any @xmath315 , and then by applying the dalembert ratio test we conclude that every individual series is absolutely convergent .
this result implies that the final sum for @xmath210 is independent of the element ordering in the alternating series .
in particular , we can add the terms @xmath316 over all @xmath315 , for @xmath310 , @xmath311 , ... ,
@xmath136 , before performing the sum over @xmath174 . by doing that , we find the simple result @xmath317 for any value of @xmath174 , and hence the sum over @xmath174 gives zero .
therefore , the sum of the series ( [ comb1 ] ) is @xmath211 .