a variety of nonequilibrium phenomena such as open flows in hydrodynamics , thermal convection in pure fluids and binary mixtures , processes in lasers and nucleation during the first - order phase transitions , can be modelled by the complex ginzburg - landau equations .
( for review and reference , see e.g. @xcite . ) of primary importance are pulse - like solitary wave solutions , which represent localised structures widely observed in nonequilibrium systems . it is well known that in the cubic ginzburg - landau equation , @xmath2 solitary waves are unstable for all positive @xmath0 , @xmath3 and real @xmath4 . in the @xmath5 case
, however , they can be stabilised by the addition of quintic terms : @xmath6 with positive @xmath7 @xcite .
this example suggests that solitary pulses can be stable in a more general class of ginzburg - landau equations where the zero solution undergoes a subcritical ( rather than supercritical ) bifurcation to a flat nonzero solution @xcite . in eq.([glquintic ] ) , the terms with @xmath3 and @xmath0 account for linear homogeneous and nonhomogeneous losses , respectively , while the terms with @xmath4 and @xmath7 describe the cubic gain and quintic dissipation .    the present work deals with another equation of the ginzburg - landau type exhibiting the subcritical bifurcation , viz .
the _ parametrically driven _
ginzburg - landau : @xmath8 here , as in ( [ glcubic])-([glquintic ] ) , the positive @xmath0 and @xmath3 are the homogeneous loss and diffusion coefficients , respectively .
the term with an asterisk ( indicating the complex conjugation ) represents the parametric driver .
( the driver s amplitude @xmath9 can be chosen positive without loss of generality . )
when @xmath10 , eq.([eq : ddcgle_bis ] ) gives the parametrically driven damped nonlinear schrdinger equation .
this special case has been studied extensively ; in particular , stable solitary waves @xcite and their stable complexes @xcite were found , and their bifurcations and supercritical dynamics analysed @xcite .
the objective of the present work is to advance beyond the nonlinear schrdinger limit .
we will still keep @xmath11 but allow for a nonzero diffusion coefficient , @xmath0 .
as we will see , even such a minimal generalisation gives rise to a new phenomenology of localised solutions which includes the multistability of pulses and pulse - front transitions .
like the ginzburg - landau equations with intrinsic gain , the parametrically driven equation ( [ eq : ddcgle_bis ] ) arises in a wide range of physical applications .
these include nonlinear faraday resonance in vertically vibrated layers of water @xcite and nonlinear lattices @xcite ; commensurate - incommensurate transitions in convective systems @xcite ; waves in nematic and cholesteric liquid crystals in rotating magnetic fields @xcite ; magnetisation waves in easy - plane ferromagnet exposed to microwave fields @xcite ; domain walls in the easy - axis ferromagnet near the curie temperature @xcite and in the easy - plane magnet in the stationary magnetic field @xcite ; nonlinear fiber lines with phase - sensitive amplification and mean - field models of degenerate optical parametric oscillators under continuous - wave pumping @xcite ; pulsed optical parametric oscillators with spectral filtering and lasers with intracavity parametric amplification @xcite . in most cases the models considered in literature are either purely diffusive ( @xmath12 ) @xcite or purely dispersive ( @xmath1 ) @xcite .
however there are situations where it is crucial that @xmath0 be finite but nonzero .
one physical phenomenon to which both diffusion and dispersion make essential contributions , is faraday resonance in strongly dissipative media such as viscous fluids and granular materials @xcite .
another situation where @xmath13 , corresponds to nondegenerate optical parametric oscillators @xcite and fiber - optic telecommunication links ; in these contexts , the term @xmath14 represents spectral filtering in the spatial or temporal domain , respectively . in fact the applicability of eq.([eq : ddcgle_bis ] ) with nonzero @xmath0 and @xmath4 may happen to be even wider ; it is commonly held that this equation provides a phenomenological description to a broad range of pattern - forming systems @xcite .    to study the solitary wave phenomenology introduced by taking the ( @xmath15)-term into account , we will use the diffusion - free limit ( @xmath1 ) as a starting point , and perform the numerical continuation of analytical solutions available in that case , to nonzero @xmath0 .
the stability of solutions obtained in this way will also be studied numerically .
we will show that  switching on " the diffusion gives rise to a sequence of stable multihumped pulses occurring in the vicinity of a certain particular value of the diffusion coefficient , @xmath16 .
the closer @xmath0 is to @xmath17 , the greater is the number of multihump solutions coexisting at this @xmath0 .
the solutions with more than 5 or 6 humps describe a flat  plateau " ( where @xmath18 is equal to the stable flat nonzero solution ) sandwiched between two fronts of opposite polarity .
the paper is organised as follows .
section [ flat_backgrounds ] deals mainly with spatially homogeneous solutions .
in particular , we show that there is a stable flat nonzero solution for sufficiently large @xmath0 , and this uniform solution can serve as a background for solitary waves . in section [ continua ]
we use perturbation - type arguments to demonstrate that both @xmath19 and @xmath20 solitons are continuable in @xmath0 and in section [ sec : adiabatic ] we construct the solutions with nonzero @xmath0 in the adiabatic approximation . the actual continuation is carried out numerically in section [ sec : continuation ] where we also examine the stability of the continued solutions . some additional insight into the bifurcation of stationary pulses
is gained in section [ h_vs_c ] .
finally , section [ sec : conclusion ] summarises results of this work .
we close this introduction by mentioning a recent publication @xcite which was devoted to the study of the single - humped solution of eq.([eq : ddcgle_bis ] ) , by means of an averaging technique and direct numerical simulations .
neither the multistability of pulses nor the pulse - front transitions were dealt with in ref.@xcite .
the transformation @xmath21 casts eq.([eq : ddcgle_bis ] ) in a ` standard ' autonomous form @xmath22 where we set @xmath11 and rescaled @xmath23 so that @xmath24 .
this is the representation that we are going to work with in what follows .
the equation ( [ eq : ddcgle ] ) has three time - independent spatially uniform solutions , or ` flat backgrounds ' , for short .
( we do not distinguish between solutions different only in the overall sign here . )
one flat solution is @xmath25 ; it will be central for this work where we are focussing on solutions decaying to zero at infinities .
the other two flat solutions are given by @xmath26 where @xmath27 note that these solutions do not depend on @xmath0 .
one can easily check that the zero solution is stable as long as @xmath28 , irrespectively of the value of @xmath29 .
the analysis of the flat nonzero solutions is somewhat more laborious .    linearising eq.([eq : ddcgle ] ) about @xmath30 , and assuming a perturbation @xmath31 $ ] yields the dispersion relation @xmath32 where @xmath33 ( here @xmath34 stands for @xmath35 or @xmath36 , depending on which solution we are linearizing about . )
the flat solution is stable iff @xmath37 , i.e. , when @xmath38 , for all real @xmath39
. the latter condition amounts to an inequality @xmath40 where @xmath41 stands for @xmath42 .
let , first , @xmath43 . since @xmath44 , the inequality ( [ inequalityfork ] ) does not hold for @xmath45 ; hence the ` low ' background @xmath46 is always unstable .
let now @xmath47 , the ` high ' background . the quadratic expression in ( [ inequalityfork ] ) will be positive for all @xmath48 if either its two roots @xmath49 and @xmath50 are both real negative or complex .
( note that we can not have two real roots of opposite signs as the quadratic s constant term , @xmath51 , is always positive for @xmath52 . )
whether the roots are real or complex is determined by the discriminant of the quadratic ( [ inequalityfork ] ) which can be written as @xmath53(c - c_+)(c - c_-).\ ] ] here we have introduced @xmath54    we need to consider two cases .
assume , first , that @xmath55 .
in this case the discriminant is negative ( and hence , the roots @xmath56 are complex ) provided @xmath0 lies in the interval @xmath57 . on the other hand , the roots @xmath56 are real negative in this case if @xmath58 is @xmath59 and the coefficient in front of the middle term in ( [ inequalityfork ] ) is positive : @xmath60 since in the case at hand we have @xmath61 , the union of the above two ` stable ' regions , @xmath57 and @xmath62 , with the endpoints included , is simply @xmath63 .
the second case is defined by inequality @xmath64 .
here we have @xmath65 and the quadratic ( [ inequalityfork ] ) can not have negative real roots as the region ( [ c0 ] ) does not overlap with the region where @xmath66 .
however , it can have complex roots  provided @xmath67 .
thus we arrive at a simple stability criterion for the flat nonzero solution @xmath68 , valid for all @xmath9 and @xmath3 : @xmath69 with @xmath70 as in ( [ eq : cpmroots ] ) .
note that the stability threshold @xmath71 is always strictly positive .
this implies , in particular , that the solution @xmath68 is always unstable in the case @xmath1 .
( this fact has already been mentioned in the literature @xcite . ) therefore the ` focusing ' , or ` attractive ' , damped - driven schrdinger equation ( eq.([eq : ddcgle ] ) with @xmath1 ) does not have stable flat backgrounds except the trivial one , @xmath72 .
the analysis of localised solutions over flat _ nonzero _
backgrounds becomes meaningful only within the full ginzburg - landau equation , i.e. eq.([eq : ddcgle ] ) with positive @xmath0 .
no less important is the question of when a flat solution can serve as an asymptotic value to a localised waveform  in other words , when is spatial decay to the flat solution possible . to find the corresponding criterion
, we set @xmath73 in equation ( [ eq : dispersioneqn ] ) .
this results in a quadratic equation @xmath74 where @xmath75 . the spatial decay to a flat background is possible unless both roots of eq.([eq : squad ] ) , @xmath49 and @xmath50 , are positive .    in the case of the @xmath46 background ,
the discriminant of equation ( [ eq : squad ] ) is positive while the constant term is negative , whence @xmath76 and @xmath77 .
consequently , the decay to @xmath46 may occur for all values of @xmath9 , @xmath78 and @xmath0 .
( this fact is of little importance , however , since we have just shown that this background is always unstable . )    in the case of the @xmath68 solution , we _ may _ have two positive roots  provided @xmath79 and @xmath80 , with @xmath81 as in ( [ c0 ] ) .
following the steps in the previous subsection , one can readily show that this situation occurs only if @xmath82 , with @xmath70 as in ( [ eq : cpmroots ] ) .
thus the _ stable _ constant solutions , defined by inequality ( [ flat_stab ] ) , can always serve as backgrounds to fronts and pulses .
we will come across localised solutions over nonvanishing backgrounds in section [ sec : continuation ] below .    finally , to examine the case of the zero background , we set both @xmath83 and @xmath73 in equation ( [ eq : dispersioneqn ] ) .
this yields @xmath84 where , as before , @xmath75 .
since the middle term in ( [ zero_decay ] ) has a positive coefficient , we have @xmath85 , which means that the situation where both @xmath49 and @xmath50 are positive , can never occur in this case .
thus the decay to the zero background is possible for all @xmath9 , @xmath3 and @xmath0 .
our last result in this section concerns the range of driving strengths which can support localised solutions in the presence of dissipation .
multiplying eq.([eq : ddcgle ] ) by @xmath86 and subtracting the complex conjugate of the resulting equation gives what would be a local conservation law of the number of particles if @xmath0 , @xmath9 and @xmath3 were equal to zero : @xmath87 -2i\gamma \left| \psi \right| ^{2 } + ic\left ( \psi _ { xx}\psi ^{\ast } + \psi _ { xx}^{\ast } \psi
\right ) \ , .
\label{eq : ada - eq1x}\end{aligned}\ ] ] assuming that @xmath18 and @xmath88 as @xmath89 and integrating ( [ eq : ada - eq1x ] ) over the real line , we get @xmath90 dx-2c \int|\psi_x|^2 dx , \label{decay}\ ] ] where we have denoted @xmath91 the phase of the complex field @xmath18 : @xmath92 .
let @xmath93 ( and remember that @xmath94 ) . in this case
it follows from ( [ decay ] ) that the time derivative of @xmath95 remains negative at all times .
hence as @xmath96 , @xmath97 for all @xmath98 .
no stationary , time - periodic , quasiperiodic or chaotic solutions , decaying to zero as @xmath89 , can arise if @xmath93 .
we will make use of this simple criterion in what follows .
in the limit @xmath1 , eq.([eq : ddcgle ] ) becomes the parametrically driven damped nonlinear schrdinger ( nls ) equation , which has exact time - independent solitary wave solutions @xcite of the form @xmath99 with @xmath100 and @xmath101 as in ( [ plusminus ] ) .
the solution @xmath102 is unstable for all @xmath9 and @xmath78 , while @xmath103 is stable in a certain parameter region @xcite .
the purpose of this section is to show that the solitary pulse solutions @xmath104 of the nls equation persist for @xmath105 .
we restrict ourselves to stationary solutions ( @xmath106 ) .
writing @xmath107 with @xmath108 a constant phase to be chosen later , eq.([eq : ddcgle ] ) becomes an ordinary differential equation for @xmath109 : @xmath110 to find out whether solutions available at @xmath1 can be continued to nonzero @xmath0 , we expand @xmath111 in power series @xmath112 , substitute into eq.([eq : eqnforphi ] ) and match like powers of @xmath0 .
it is convenient to choose @xmath108 to be @xmath113 in the case of the soliton @xmath103 , and @xmath114 for @xmath102 ; this choice makes @xmath115 real . matching terms linear in @xmath0 and decomposing @xmath116 into its real and imaginary part ( @xmath117 ) , yields an equation @xmath118 where the primes stand for the derivatives in @xmath98 and the operators @xmath119 are defined by    @xmath120    ( the subscripts @xmath121 and @xmath122 indicate whether we are checking the continuability of @xmath103 or @xmath102 . )    according to fredholm s alternative , eq.([eq : eqnuv ] ) has bounded solutions if and only if its right - hand side is orthogonal to the kernel of the adjoint operator @xmath123 . transforming to @xmath124 and @xmath125 in the case of the @xmath103 and @xmath102 solution , respectively , the equation for the zero modes spanning @xmath126 becomes @xmath127 where @xmath128 , @xmath129 , @xmath130 and we have introduced the standard pschl - teller operators with familiar spectral properties : @xmath131    the only discrete eigenvalue of @xmath132 is @xmath133 , while the continuous spectrum occupies the semiaxis @xmath134 .
therefore , the operator @xmath135 is invertible as long as @xmath136 and @xmath137 . assuming that the zero background is stable ( i.e. , assuming that @xmath138 ) , the quantity @xmath139 is indeed less than one but greater than zero . on the other hand ,
the @xmath140 is negative .
( an exception is the point @xmath141 , where we have @xmath142 . )
thus we conclude that @xmath143 is invertible except in the special case @xmath141 .
( however , even in this special case the operator @xmath144 is invertible on _ odd _ functions because the eigenfunction @xmath145 associated with the zero eigenvalue is even . )
consequently , @xmath146 can not be equal to zero for @xmath147  otherwise the bottom equation in ( [ eq : ker - l - dag ] ) would imply that @xmath148 , too .
fortunately , the operator @xmath149 does have a zero eigenvalue , and therefore @xmath146 can be a nonzero multiple of the corresponding eigenfunction ( which is @xmath150 . )
thus in the case @xmath151 the kernel of @xmath152 is spanned by just one zero mode , namely @xmath153 on the other hand , when @xmath141 , the dimension of the kernel space is two .
firstly , the zero mode ( [ eq : zerom - hneqg ] ) persists for @xmath154 as the operator @xmath155 is defined on odd functions .
secondly , there is another zero mode given by @xmath156    it is obvious that the vector - function ( [ eq : zerom - hneqg ] ) , whose both components are odd functions of @xmath157 , is orthogonal to the right - hand side of eq.([eq : eqnuv ] )  which is even .
hence eq.([eq : eqnuv ] ) is solvable for @xmath158 .
it is also easy to check that the mode ( [ eq : zerom ] ) is _ not _ orthogonal to the right - hand side of ( [ eq : eqnuv ] ) and so the solvability condition is _ not _
satisfied for @xmath141 .
thus , having started from the two nonlinear schrdinger solitons , we constructed two one - pulse solutions of the ginzburg - landau equation to within the first order in the small parameter @xmath0 : @xmath159 .
consequently , we expect to be able to continue the schrdinger solitons @xmath160 into the region @xmath161 ( provided that @xmath158 . )
this expectation is born out by results displayed in section [ sec : continuation ] below , which present an outcome of the _ numerical _ continuation of @xmath104 in @xmath0 .
before attempting the full - scale numerical continuation , it is instructive to construct approximate solutions . our approximation will be valid for small @xmath0 and exploit the fact that when @xmath1 , the stationary pulse - like solutions of equation ( [ eq : ddcgle ] ) have the form @xmath162 where @xmath163 and @xmath164 are constants defined by eqs.([plusminus ] ) .
for @xmath0 small nonzero , approximate solutions can be sought for by assuming that @xmath18 retains the form ( [ eq : ada - soln ] ) , but @xmath165 and @xmath166 become slowly varying functions of @xmath23 .    to obtain an expression for @xmath167 ( the overdot stands for the derivative in @xmath23 )
, we substitute the ansatz ( [ eq : ada - soln ] ) into ( [ eq : ada - eq1x ] ) . integrating over @xmath98 and using the boundary conditions @xmath168 as @xmath169 produces an evolution equation for the pulse s amplitude : @xmath170 where @xmath171 .
an equation for the pulse s phase arises by multiplying ( [ eq : ddcgle ] ) by @xmath172 , adding with its complex conjugate , substituting ( [ eq : ada - soln ] ) for @xmath18 and integrating over @xmath98 : @xmath173    fixed points of the system ( [ eq : ada - eq1 ] ) and ( [ eq : ada - eq2 ] ) correspond to stationary solutions of ( [ eq : ddcgle ] ) .
these can be easily found explicitly : eliminating @xmath174 from the stationary system @xmath175 produces a quadratic equation @xmath176 with roots @xmath177 the corresponding @xmath178 are defined by their sine and cosine in ( [ eq : stat - eq2 ] ) .
it is not difficult to determine when the roots ( [ eq : ada - soln - for - a ] ) are real and positive .
we assume that the zero background solution of eq.([eq : ddcgle ] ) is stable , i.e. the constant term in ( [ eq : stat - eqnfora ] ) is positive . hence if the roots are real , they are of the same sign , and this sign is opposite to that of the middle term in ( [ eq : stat - eqnfora ] )
. therefore we have two positive roots provided the discriminant @xmath179 is nonnegative , and , at the same time , the inequality @xmath180 holds true . in ( [ eq : ada - realaa ] ) we have introduced @xmath181 where the @xmath121 corresponds to @xmath182 and @xmath122 to @xmath183 .
it is straightforward to verify that for small @xmath9 , @xmath184 , we have @xmath185 , while for larger @xmath9 , @xmath186 , we have @xmath187 .
( here we have made use of the inequality @xmath188 ; as proved in section [ sec : threshold ] , no stationary localised solutions exist for @xmath189 . ) therefore , the region of @xmath0 values where the roots of ( [ eq : stat - eqnfora ] ) are positive , is given by the inequality @xmath190  for all @xmath9 .
thus we conclude that the adiabatic equations have two stationary points for @xmath191 , and none for @xmath192 .
we complete the adiabatic analysis by classifying their stability and bifurcation .
linearising equations ( [ eq : ada - eq2 ] ) about the stationary points and assuming small perturbations of the form @xmath193 and @xmath194 , yields a characteristic equation @xmath195=0.\ ] ] since the coefficient in front of the middle term in this equation is positive , either its two roots are complex with negative real parts , or we have two real roots , of which one is negative .
the case when the second root is positive ( and hence the fixed point is unstable ) occurs if the constant term is negative .
conversely , if the constant term is nonnegative : @xmath196 \ge 0 , \label{stabby}\ ] ] the fixed point is stable . simplifying eq.([stabby ] )
, we get the stability condition in the form @xmath197 comparing this to the expressions for @xmath198 in ( [ eq : ada - soln - for - a ] ) , we conclude that for all @xmath0 , @xmath9 and @xmath3 , the points @xmath199 and @xmath200 are a stable node and saddle , respectively .
the two fixed points come into being through a saddle - node bifurcation which occurs as the diffusion coefficient @xmath0 is decreased past @xmath201 for the fixed @xmath9 and @xmath3 or , alternatively , as the driving amplitude @xmath9 is increased for the fixed dissipation coefficients @xmath0 and @xmath3 .
one can easily find the bifurcation value of @xmath9 , at which two complex roots of the quadratic equation ( [ eq : stat - eqnfora ] ) converge on the positive real axis .
( here we are assuming that @xmath202 is smaller than @xmath203 . )
writing the discriminant ( [ eq : ada - realaa ] ) as @xmath204,\ ] ] the threshold driving strength is given by @xmath205 where the subscript  ad " serves to remind that eq.([threshold ] ) was obtained in the adiabatic approximation .
it is important to emphasise that this formula is valid only for small @xmath202 .
equivalently , the expression ( [ eq : ada - soln - for - c ] ) for the turning point @xmath206 is valid only for @xmath9 close to @xmath3 .
in this section we describe the bifurcation diagram obtained by the numerical continuation of the solitons @xmath207 in the parameter @xmath0 .
the diagram is presented in fig.[biffy ] and displays the sobolev norm of the solution , @xmath208 as a function of @xmath0 .      for stationary solutions
, @xmath209 , equation ( [ eq : ddcgle ] ) reduces to an ordinary differential equation @xmath210 for the numerical continuation of solutions to eq.([stat_gle ] ) we utilised the auto97 software package @xcite .
the infinite line was approximated by a finite interval @xmath211 , with @xmath212 , and the boundary conditions @xmath213 .
the tolerance of the computation ( the grid norm of the difference between the right and left - hand side of ( [ stat_gle ] ) ) was set to be @xmath214 .
the solid and broken curves in fig.[biffy ] represent stable and unstable branches , respectively .
the stability was examined by linearising eq.([eq : ddcgle ] ) about the corresponding stationary solution . choosing the small perturbation in the form
@xmath215 e^{\lambda t}$ ] , with real @xmath216 and @xmath217 , we arrive at an eigenvalue problem @xmath218 where    @xmath219    and @xmath220 the @xmath221 and @xmath222 are the real and imaginary parts of the solution : @xmath223 .
we solved the eigenvalue problem by expanding @xmath224 and @xmath225 over 500 fourier modes in the interval @xmath226 .
the continuation in @xmath0 was performed for fixed values of @xmath9 and @xmath3 .
we selected @xmath227 and @xmath228 ; for these @xmath9 and @xmath3 the nonlinear schrdinger soliton @xmath19 is stable @xcite . before proceeding to the @xmath19 soliton
, however , we briefly deal with the @xmath20-case . in agreement with predictions of section [ continua ] , the soliton @xmath20 was found to persist both for @xmath229 and @xmath94 .
we were in fact able to continue it indefinitely without encountering any obstacles in both directions .
as @xmath230 , the width of the localised solution of ( [ stat_gle ] ) grows and the solution tends to @xmath231 , where @xmath232 stands for a pulse - shaped solution of equation @xmath233 in a similar way , as @xmath234 , the solution tends to @xmath235 . when continued to @xmath94
, the pulse remained unstable for all @xmath0 , with a single positive eigenvalue in its spectrum .
( as for the negative-@xmath0 region , all solutions there are _ a priori _ unstable against continuous spectrum perturbations with arbitrarily large @xmath236 . )    the continuation of the @xmath19 soliton proved to be more rewarding from the stability viewpoint .
the corresponding bifurcation diagram is displayed in fig.[biffy ] .
as with the soliton @xmath20 , the @xmath19 persists both for @xmath237 and @xmath238  in agreement with section [ continua ] .
continuing into the @xmath94 region , we have found that the solution gradually changes its shape , with the hump in the imaginary part splitting into two ( see fig.[conti](a ) ) . at @xmath239
the branch turns back ( fig.[biffy ] ) .
note that the turning point occurs for @xmath0 much smaller than @xmath240 , the saddle - node bifurcation point predicted by the adiabatic analysis .
what is more important , the @xmath19 solution turns not into the @xmath20 ( as the adiabatic approach predicted ) but into some other solution which has two well - separated humps in the imaginary part .
the discrepancy is not surprising as the adiabatic approximation is valid only for small @xmath0 , where the shape of the solution is still well reproduced by the single - hump constant - phase trial function ( [ eq : ada - soln ] ) . in the entire stretch between @xmath1 and the turning point ,
the solution remains stable .
continuing this branch further , additional humps are added as the solution passes a sequence of turning points .
each pass of a turning point results in the creation of a new hump in the middle of the solitary wave .
as we move along the branch , the distance between successive turning points ( the difference between the corresponding values of @xmath0 ) becomes smaller and the new humps come with smaller amplitudes . as a result ,
a long plateau is formed which keeps on expanding as we continue the branch ( fig.[conti](c ) ) .
the broadening plateau accounts for the vertical segment of the curve in fig.[biffy ] .
if the bifurcation parameter @xmath0 is seen as a function of the sobolev norm @xmath241 ( that is , if we turn fig.[biffy ] by @xmath242 ) , the curve @xmath243 has the form of a decaying oscillation .
the stability of the solution alternates at each successive turning point .
these changes are due to a single real eigenvalue which moves back and forth through the origin on the real line .
( at the turning points the eigenvalue is right at the origin , of course . )
the lengths of the incursions this eigenvalue makes into the positive and negative real lines decrease with each new turning point until the eigenvalue becomes indistinguishable from zero .
therefore the branch becomes ( neutrally ) stable sufficiently ` high up ' in @xmath241 in fig.[biffy ] ( i.e. for sufficiently long plateaus . )    despite the fact that negative values of the diffusivity @xmath0 are not physically meaningful , we did continue to @xmath229  in the hope that the resulting branch would reach a turning point and then return to the positive @xmath0-semiaxis ( which it did indeed do ) . the resulting branch is also shown in fig.[biffy ] .
as we move into the region @xmath229 , the @xmath19 solution gradually develops into a three - humped state and when the curve returns to @xmath1 , we have a complex of the @xmath244 type , with a large separation between the individual solitons .
( this complex of three solitons of the parametrically driven nonlinear schrdinger equation was previously found in ref.@xcite . ) crossing into the @xmath94 region , the central @xmath19 soliton in the complex transforms as if it did not have the @xmath20 solitons attached to its flanks . as a result ,
the @xmath94-portion of the corresponding @xmath245 curve has virtually the same shape as the curve resulting from the continuation of @xmath19 to the region @xmath94 ; the only difference is that the curve emanating out of @xmath244 is shifted upwards relative to the curve emanating out of @xmath19 .
similarly to the continuation of @xmath19 , the continuation of @xmath244 goes via a series of turning points , with each pass of the turning point resulting in the creation of another hump in the middle of the central region which becomes a long plateau .
the lateral @xmath20 solitons are not affected by this process .
the linearised spectrum is the union of the spectrum of the long pulse described in the previous paragraph and spectra of two @xmath20 solitons .
in particular , it includes two positive real eigenvalues contributed by the @xmath20 s and so the entire branch resulting from the continuation of @xmath244 is unstable .
we disregard it in what follows .
the plateau arising in the final stage of continuation of the two branches shown in fig.[biffy ] , is nothing but an interval on the @xmath98-axis where @xmath18 equals @xmath68 , the flat nonzero solution given by eq.([flat_back ] ) .
the corresponding value of @xmath0 , @xmath246 , falls within the region @xmath63 where the background @xmath68 is stable . here
@xmath247 is given by eq.([eq : cpmroots ] ) ; in particular , @xmath248 . in the language of phase transitions , the long pulse shown in fig.[conti](c ) can be seen as a  bubble " of one stable phase in another one .
the long pulse can also be interpreted as a bound state of two fronts interpolating between different stable backgrounds , @xmath72 and @xmath52 .
this intuitively - appealing idea can be put on the quantitative footing by considering the spatial decay of perturbations to the flat background @xmath68 .
indeed , consider the turning point separating the second stable branch from the first unstable branch in fig.[biffy ] .
( this is a point with coordinates @xmath249 and @xmath250 .
note that we are only considering branches obtained by the continuation of the soliton @xmath19 to positive @xmath0 .
branches obtained by the continuation of @xmath20 to negative @xmath0 first , are disregarded here . )
it is at this point that the modulus squared of @xmath251 becomes double - humped ; before that , that is on the branch that leads to this point ( the first unstable branch ) , the function @xmath252 remained single - humped despite the double - humped imaginary part . letting @xmath228 , @xmath253 and @xmath249 we check that @xmath64 and hence , according to section [ spatially_uniform ] ,
the quadratic ( [ inequalityfork ] ) , ( [ eq : squad ] ) has two complex roots , @xmath254 and @xmath255 .
therefore the wavenumber @xmath39 is complex as well : @xmath256 . solving equation ( [ eq : squad ] ) for @xmath41 we obtain , subsequently , @xmath257 ( we have chosen positive values for @xmath258 and @xmath259 here . )    for @xmath0 near the turning point in question
, the plateau has not yet formed between the two humps and so they can be crudely thought of as two overlapping nls - like solitons with oscillations on their adjacent ( i.e. partner - facing ) tails .
the bound state arises when one solitary wave is trapped in the potential well formed by the oscillatory tail of its partner . making use of the potential of interaction of two attractive nls solitons @xcite , @xmath260 ( where @xmath261 stands for the distance between the two humps ) , and taking into account that @xmath262 , we obtain a rough estimate for the separation : @xmath263 .
this formula could also be derived from purely qualitative considerations ; all one needs to notice is that the absolute value of the field @xmath18 obtained from the linear superposition of two in - phase solitons is maximised ( and so the energy is minimised ) by placing one soliton at the first maximum of @xmath264 on the tail of its partner .
substituting @xmath258 from ( [ kk ] ) , the approximate formula gives @xmath265 .
this is in qualitative agreement with the numerically found value @xmath266 .    moving further along the ( second stable ) branch
, the plateau appears and the two humps can no longer be approximated by the nls - like solitons .
this makes the above estimate invalid .
the solution can still be regarded as a bound state of two fronts but this time , in order to calculate the characteristic width of the pulse one would need to know the full profile of the front .
we conjecture that for a given @xmath9 and @xmath3 , a free - standing stationary front exists just for a single value of @xmath0 , namely @xmath267 .
( we are planning to verify this conjecture in future publications . ) on the other hand , stable and unstable bound states of fronts exist in a finite interval of @xmath0 values containing @xmath267 as an internal point .
it is fitting to note here that similar pulse - to - front transformations occur also in the other system featuring the subcritical bifurcation , viz . the cubic - quintic ginzburg - landau equation with internal gain @xcite .
as we mentioned in the previous section , there is a certain discrepancy between the adiabatic analysis and numerical continuation .
numerically , the @xmath268 soliton was found to be continuable all the way to @xmath269 whereas the adiabatic approach predicted the existence of a turning point at @xmath270 , where the @xmath268 should have merged with the @xmath271 branch . as for the @xmath19 solution
, we found that it turns into a pulse with the double - humped imaginary part ( and not into the @xmath268 branch as suggested by the adiabatic approximation . ) in order to shed some light on the possible source of the discrepancy we performed the numerical continuation of the pulse @xmath268 in @xmath9 , for several fixed values of @xmath0 . here , by the @xmath268 we mean the pulse solution which results from the continuation of the schrdinger @xmath268 soliton to positive @xmath0 , for some fixed large value of @xmath9 ( in our case for @xmath272 . )
having obtained this starting - point solution for several values of @xmath0 , we then continued it in @xmath9 , from @xmath272 to smaller @xmath9 .
the stability of the arising solutions was examined by computing eigenvalues of the operator ( [ ev_problem ] ) at sample values of @xmath9 .    in each case
considered , the @xmath268 branch was found to turn into the @xmath271 solution as @xmath9 reached the threshold value @xmath273 .
( that is , for @xmath274 there are two branches of solutions whereas for @xmath275 , there is none ; see fig.[conti_in_h](a ) ) .
the entire @xmath20 branch is unstable ; the single positive real eigenvalue moves to the negative semiaxis as the branch is continued past the turning point .
continuing the arising @xmath271 branch to larger @xmath9 , we reach another turning point at @xmath276 , where the @xmath271 solution transforms into a pulse with the double - humped imaginary part .
the values @xmath277 and @xmath278 are shown in fig.[conti_in_h](b ) , as functions of @xmath0 ( for the fixed @xmath279 ) .
as @xmath280 , the difference @xmath281 decreases but remains nonzero .
we verified this by computing @xmath277 and @xmath278 for equation ( [ c_infty ] ) which pertains to @xmath12 . in the same plot
we display the function ( [ threshold ] ) which gives the adiabatic approximation to the curve @xmath282 .
( note that for small @xmath0 , there is a good agreement between numerical and approximate values but as @xmath0 grows , the two curves diverge . )    continuing the @xmath271 branch past the second turning point , the solution adds another hump in the middle of the pulse , turns back again , adds another one , and so on .
( see fig.[conti_in_h](a ) ) .
a long plateau developes in the middle of the pulse , just like when it was continued in @xmath0 .
similarly to the @xmath0-continuation , the turning points on the @xmath9-axis separate regions of stability from regions of instability , with the instability being caused by a single positive real eigenvalue ( which moves through @xmath283 at the turning points . )    from figure [ conti_in_h](b ) it is clear why the saddle - node bifurcation point where the @xmath271 and @xmath268 solutions would merge , did not appear in fig.[biffy ] .
the reason is that fig.[biffy ] was plotted for a relatively large value of @xmath9 @xmath284 whereas according to fig.[conti_in_h](b ) , a horisontal line @xmath285 with @xmath286 can have no intersections with the saddle - node curve @xmath287 .
( here all the numbers are for @xmath279 . ) the same fig.[conti_in_h](b ) explains what semed to be a discrepancy between the adiabatic analysis and the numerical continuation of the soliton @xmath20 in @xmath0 . the numerical result that seemed to contradict
the adiabatics was that for @xmath228 , @xmath20 could be continued without bounds .
it is now obvious from fig.[conti_in_h](b ) that the unbounded continuation is only possible for @xmath9 greater than @xmath288 .
continuing the @xmath268 soliton to positive @xmath0 for @xmath9 _ smaller _ than @xmath288 , the branch turns back ( already as the @xmath271 pulse ) after hitting the lower solid curve in fig.[conti_in_h](b ) .
therefore , the pattern arising for @xmath9 close to @xmath3 actually _ is _ in qualitative agreement with the adiabatic analysis , which was expected to be valid precisely for small @xmath0 or , equivalently , for small @xmath289 differences .
in this paper we studied a cubic complex ginzburg - landau equation in which linear losses and diffusion are compensated by the linear parametric drive .
the nonlinear term in the equation was taken to be purely conservative .
there are three stationary homogeneous solutions to equation ( [ eq : ddcgle ] ) , and we have shown that the @xmath72 solution is stable if @xmath290 , as long as @xmath29 .
this stability condition coincides with the corresponding condition for the nonlinear schrdinger case ( i.e. for @xmath1 ) . on the other hand ,
the stability properties of the _ nonzero _ homogeneous solutions are not the same as in the @xmath1 case . indeed , unlike for @xmath1
, there is a stable flat _
solution @xmath52 for sufficiently large diffusion coefficients , @xmath291 .    having established the persistence of the nls solitons @xmath19 and @xmath20 for small nonzero @xmath0 , we continued them in @xmath0 numerically .
the continuation of the soliton @xmath19 yields a sequence of pulse - like solutions , separated by turning points , with increasing number of humps .
the stability of these solutions changes at each turning point , so that stable multihump solutions coexist with unstable ones .
it is fitting to note here that the multistability of multipulse solutions is _ not _ observed in the schrdinger limit where only the two - soliton complex was found to be stable @xcite . as @xmath292 , where @xmath293 , the solution takes the form of a long plateau ( an interval of the stable background @xmath68 ) sandwiched between two fronts .
we also performed the continuation in @xmath9 , for the fixed @xmath0 . for each @xmath94 ,
two localised solutions are born in a saddle - node bifurcation as @xmath9 exceeds a threshold value .
( we obtained an analytic formula for the threshold in the adiabatic approximation ; it is in agreement with the numerical results for small @xmath0 . )
the subsequent continuation gives rise to a sequence of coexisting stable multihump solutions culminating in a bound state of two widely separated fronts .
solitary pulses in the form of long shelves ( plateaus ) can allow easy experimental observation in physical systems described by our model .
in particular , they may be employed as a natural basis for the non - return - to - zero ( nrz ) format of the data transmission in optical telecommunications . in the nrz format , the 1 and 0 bits are coded , respectively , by sending or withholding the signal within a standard time slot .
a string of several 1 s looks as a long uniform pulse of an essentially arbitrary length .
the stability of such pulses is crucial to maintain the fixed shape of the long array of 1 s , and to prevent the _ inter - symbol interference _ ,
i.e. , the blurring of empty intervals between such strings , which represent ( strings of ) 0 s ( see , e.g. ref . @xcite and references therein . ) in the case of lasers , which can also be described by the present model @xcite , the possibility of the generation of stable long pulses of arbitrary duration , i.e. , an effective _ tunability _ of the output , is an essential advantage too .
we thank nora alexeeva for her advice on numerics .
the first author s ( i.b.s ) work was supported by the nrf of south africa under grant no.2053723 , by the johnson bequest fund and the urc of the university of cape town . the second author ( s.c . ) was supported by the nrf and the harry crossley foundation . b.a.m .
appreciates hospitality of the department of physics at the universitt erlangen - nrnberg , and a partial support from the european science foundation , through the `` @xmath294-shift '' programme .
n. akhmediev and a. ankiewicz .
nonlinear pulses and beams . chapman and hall , 1997 ; n. akhmediev and a. ankiewicz .
solitons of the complex ginzburg - landau equation . in : s. trillo , w. torruellas ( eds . ) spatial solitons .
springer series in optical sciences , vol .
82 , 2001      proceedings of the workshop on the complex ginzburg - landau equation : theoretical analysis and experimental applications in the dynamics of extended systems . edited by m. br and a. torcini .
physica d * 174 * , 1 - 220 ( 2003 )          i.v .
barashenkov and e.v .
zemlyanaya , phys .
83 , 2568 ( 1999 ) m. bondila , i.v .
barashenkov and m.m .
bogdan , physica d * 87 * , 314 ( 1995 ) ; s. longhi , g. steinmeyer , and w.s .
wong , j. opt .
* b 14 * , 2167 ( 1997 ) ; n.v .
alexeeva , i.v .
barashenkov and d.e .
pelinovsky , nonlinearity * 12 * , 103 ( 1999 ) ; v.s . shchesnovich , i.v .
barashenkov , physica d * 164 * , 83 ( 2002 )          b. denardo , b. galvin , a. greenfield , a. larraza , s. putterman , and w. wright , phys .
lett . * 68 * , 1731 ( 1992 ) ; s .- y . lou and g. huang , mod .
b * 9 * , 1231 ( 1995 ) ; g. huang , s .- y .
lou , and m. velarde , int .
j. bif . and chaos * 6 * , 1775 ( 1996 ) ; j. yu , phys .
lett . a * 259 * , 366 ( 1999 ) ; n.v .
alexeeva , i.v .
barashenkov and g.p .
tsironis , phys .
* 84 * , 3053 ( 2000 ) ; w. chen , b. hu , and h. zhang , phys .
b * 65 * , 13402 ( 2002 )            a. mecozzi , w.l .
kath , p. kumar , and c.g .
goedde , opt . lett . * 19 * , 2050 ( 1994 ) ; s. longhi , opt . lett .
* 20 * , 695 ( 1995 ) ; k. staliunas , j. mod
* 42 * , 1261 ( 1995 ) ; s. longhi , j. mod . opt . * 43 * , 1089 ( 1996 ) ; g .-
oppo , a.j .
scroggie , and w.j .
firth , phys .
e * 63 * , 66209 ( 2001 )          p. coullet , j. lega , b. houchmanzadeh , and j. lajzerowicz , phys .
lett . * 63 * , 1352 ( 1990 ) ; p. coullet and k. emilsson , physica * d 61 * , 119 ( 1992 ) ; b.a .
malomed and a.a .
nepomnyashchy , europhys .
* 27 * , 649 ( 1994 ) ; c. elphick , a. hagberg , b.a .
malomed and e. meron , phys .
a * 230 * , 33 ( 1997 ) ; c. utzny , w. zimmermann , and m. br , europhys . lett . * 57 * , 113 ( 2002 ) ; g.j .
de valcrcel , i. prez - arjona , and e. roldn , phys .
lett . * 89 * , 164101 ( 2002 )


the application of the general relativity on the description of astrophysical objects has been matter of interest in the theoretical physics .
exact solutions to the einstein field equations with spherical symmetry are of great relevance if systems like stars , planets or spherical clusters are considered .
now , in order to model systems like galaxies , accretion disks or rotating black holes , axially symmetric solutions must be considered . for the case of galaxies ,
the first disks models in general relativity were proposed by bonnor and sackfield in 1968 @xcite and the second one by morgan and morgan in 1969 @xcite .
these solutions increased the interest in this kind of solutions , and later , many other solution were proposed : solutions for static sources are presented in @xcite , while disks with a rotational motion ( stationary ) are presented in @xcite .    the study of the gravitational field produced by an astrophysical source , from the point of view of general relativity , has been a difficult task due to the nonlinear nature of the einstein equations . the situation is not better when we try to find solutions that describe the gravitational field of a system consisting of two or more parts .
examples of this kind of systems are the galaxies because they have ( in general ) a disk , a galactic halo , a bulb , and in many cases , a central black hole .
the case of a disk immersed in a galactic halo is interesting because the halo is the largest and most massive part of a galaxy .
this problem has already been treated by vogt and letelier in @xcite and by gutirrez , gonzlez and quevedo in @xcite . in the first paper
the authors apply the displacement , cut and reflect method to some known disk - like solutions in order to generate disks with haloes , and in the second paper the authors present a new family of disk - like solutions with halo for a conformastatic spacetime in einstein - maxwell gravity .
in this paper we present two families of axially symmetric thin disks of dust immersed in a spheroidal halo .
the disks are obtained from solutions of the einstein equations for an axially symmetric conformastatic spacetime in which the metric tensor is characterized only by one metric function . by introducing a finite discontinuity on the first derivatives of the metric tensor ,
solutions with a singularity of the delta function type are obtained , so describing thin disks .
the nonzero components of the energy - momentum tensor , both for the disk and the halo , are obtained from the einstein equations . in this way ,
the energy densities and pressures of the sources are determined . by imposing the fulfillment of all the energy conditions we obtain a constraint over the solutions , in such a way that the metric function can be properly expressed in terms of a solution of the laplace equation . by using the solution to laplace equation in cylindrical coordinates
we find infinite disks and by using the solution to laplace equation in oblate spheroidal coordinates we find finite disks . in both cases
we obtain particular solutions with energy densities and pressures well behaved everywhere .
we also show that the masses of de disks and the haloes are finite .
finally we solve the geodesic equation for circular orbits in the plane of the disk .
some preliminaries results of this work have been presented in @xcite .
the structure of this article is the follow : in section [ sec : sols ] we solve the einstein equations for a conformastatic metric in order to find the energy density and stress of the halo in terms of solutions to the laplace equation . in section [ sec : emtensor ] we apply the thin shell formalism to find a disk - like source that satisfy the energy conditions ( physically well behaved ) . in section [ sec : mass ]
we introduce the gravitational mass , both for the disk and for the halo , and we study the dynamics of a particle moving in circular orbits in the plane of the disk in order to obtain the rotational curves , which describe the radial dependency of the speed of a particle that moves in circular paths around the center of the galaxy and in the plane of the disk .
we then present , in section [ sec : kuzmin ] , a particular family of relativistic solutions of infinite thin disks with halo and we analyze the behavior of the physical properties and , in section [ sec : morgan ] , we present a second family of solutions of finite thin disks with halo and discuss the physical features of this solutions .
finally , in section [ sec : concluding ] , we present a brief discussion of our main results .
in order to formulate the einstein field equations for an axially symmetric static spacetime with an infinitesimally thin discoidal source which is immersed in a fluid or matter halo , we first introduce the cylindrical coordinates @xmath1 .
now , we consider the spacetime as described by the conformastatic metric @xcite @xmath2 where we demand @xmath3 to be a continuos function of the @xmath4 and @xmath0 coordinates . the einstein field equations for this metric yield the energy - momentum tensor for the halo fluid @xmath5 , whose non - null components are    [ e0116 ] @xmath6    now , to give a physical interpretation to the energy - momentum tensor of the halo , we write @xmath5 in a locally minkowskian reference frame whose space - like vectors are oriented in the main stresses directions .
this reference frame is described by the orthonormal basis ,    [ e0116 ] @xmath7    where @xmath8 . by writing @xmath5 in this reference frame , we find the energy density @xmath9 the main stresses of the halo , which are given by @xmath10 and the isotropic pressure @xmath11 which is defined as the average of the main stresses .
the physical quantities we have already obtained depend on the unknown metric function @xmath3 , so we need an additional equation to get @xmath3 . on the other hand ,
if we demand @xmath12 and @xmath13 to satisfy the energy conditions @xcite , then we get the inequality @xmath14 these conditions are important to yield physically well behaved solutions , in other words , to fulfill the minimum conditions for gravitational systems .
one way to satisfy ( 7 ) is considering @xmath3 functions that are solutions to the equation @xmath15 with @xmath16 . this particular choice to fulfill the energy conditions leads us to the state equation @xmath17 for the halo . here
we can see that the pressure and the energy density are linearly related by a constant @xmath18 that depends on @xmath19 value .
it means that a specific kind of fluid in the halo depends on a particular choice of @xmath19 ; for example , if @xmath19 tends to infinity , then @xmath13 tends to zero and the fluid in the halo would be made of dust .
if we choose @xmath19 equal to 1 , then the state equation would describe a fluid made of radiation .
henceforth , we will not analyze the behavior of @xmath13 , because it is related @xmath12 via the state equation .
now , in order to get solutions for the @xmath3 function , we rewrite equation ( [ eq8 ] ) as @xmath20 in such a way that solutions for @xmath3 can be obtained form solutions of the laplace equation .
furthermore , as the spacetime needs to be asymptotically flat , then @xmath21 has to be @xmath22 at infinity , and we can write @xcite @xmath23 where @xmath24 es any solution of the laplace equation that vanishes at infinity .
this relation allows us to write @xmath12 in terms of @xmath24 as @xmath25 it is clear that we can obtain different solutions depending on the choice of @xmath24 ; nevertheless , in the next section we will see that @xmath24 must satisfy some conditions in order to have a disk solution .
in the preceding section we computed the energy - momentum tensor for the fluid in the halo ; now , we want to compute the one of the disk . to accomplish this objective we use the distributional approach @xcite and the thin shells formalism @xcite to describe the disk - like source .
mathematically the disk is described by an hypersurface @xmath26 whose function is @xmath27 and its normal vector is @xmath28 .
this hypersurface separates the spacetime in two regions : @xmath29 ( above @xmath26 ) and @xmath30 ( bellow @xmath26 ) . by using distributional functions ,
the energy - momentum tensor can be written as @xmath31 where @xmath32 is the heaviside function , @xmath33 is the dirac function , @xmath34 is the energy - momentum tensor in @xmath35 , and @xmath36 is the energy - momentum tensor in @xmath26 .    now , to describe properly a surface layer it is necesary to satisfy the condition of continuity of the metric across the plane @xmath37 , which is mathematically expressed via the jump function as @xmath38=g_{\alpha\beta}|_{_{z=0^{+}}}-g_{\alpha\beta}|_{_{z=0^{-}}}=0 .
\label{eq14}\ ] ] this condition can be set in terms of the solution to laplace equation through ( [ eq1 ] ) and ( [ eq10 ] ) as @xmath39 we also need to ensure that @xmath40 is discontinuous in the plane @xmath37 , i.e. @xmath41=2g_{\alpha\beta , z}|_{_{z=0^{+}}}\neq 0 .
\label{eq16b}\ ] ] with the metric tensor defined in ( [ eq1 ] ) and the relation ( [ eq10 ] ) it is easy to show that the last condition is satisfied if we choose a function @xmath24 such that @xmath42    on the other hand , we can compute @xmath36 through the expression @xmath43 where @xmath44 is the part of the ricci tensor associated to @xmath26 . here
it is important to note that @xmath45 by virtue of the condition ( [ eq16b ] ) .
computing @xmath36 we find that its only non - null component is @xmath46 which is evaluated in @xmath47 .
now , the surface energy - momentum tensor of the disk is computed from @xmath36 through the integral @xmath48 where @xmath49 .
this integral yield the non - null component of @xmath50 , which we write in the tetrad ( [ e0116 ] ) to interpret it physically as the surface energy density of the disk - like source , @xmath51 this solution can be rewritten in terms of @xmath24 as @xmath52 which is evaluated in @xmath47    as in the previous section , @xmath53 needs also fulfill the energy conditions , but in this case it is easier than before because our disk is made of dust ( @xmath54 ) , i.e. the particles in the fluid are collisionless ; so all the energy conditions are satisfy if @xmath55 here we do not consider the case @xmath56 because it implies @xmath57 and @xmath58 .
once the energy - momentum tensor has been specified , it is possible to compute the total mass of the system . to do it
, we use the komar formulae @xcite , which is written as @xmath59 where @xmath60 is the normal vector to the space - like hypersurface , @xmath61 is the time - like killing vector , @xmath62 is the determinant of the induced metric and @xmath63 is the volume element . in our case , @xmath64 , @xmath65 , and @xmath66 .    by using the non null components of @xmath5 and the relation ( [ eq10 ] ) in the komar formulae
we can compute the mass of the halo as @xmath67 where @xmath68 is the usual euclidean volume element of the space .
on the other hand , if we take @xmath69 in the komar formulae we find that the mass of the disk is @xmath70 finially , the total mass of the system is given by @xmath71 .
now , we can also study the movement of a particle in a spacetime described by the line element ( [ eq1 ] ) .
the particle is interacting with a gravitational field produced by a disk - like distribution of energy ( [ eq20s ] ) surrounded by a halo whose energy distribution is described by ( [ eq10u ] ) .
now , a very simple case is the one of a particle that moves in a circular orbit in the plane of the disk ( @xmath37 ) , with velocity @xmath72 tangential to its path .
this is interesting , from the physical point of view , because we have observational information for @xmath73 from spiral galaxies @xcite .
the relation between the tangential velocity and the distant from the center of the galaxy ( @xmath4 ) is described by the rotational curves .    for a circular path in the plane of the disk @xmath37 and @xmath4
is constant ; therefore , the velocity is described only by two components , @xmath74 where @xmath75 is the angular velocity . by imposing the normalization condition
@xmath76 , we find that @xmath77 where @xmath78 is the circular velocity or tangential velocity .
on the other hand , we solve the geodesic equation , @xmath79 to obtain an expression for @xmath72 in terms of the function @xmath24 .
as we saw , the spacetime has time symmetry and axial symmetry , so @xmath80 and @xmath81 are not going to change in @xmath82 , so , the geodesic equation reduces to @xmath83 , and it yield the expression for the tangential velocity , @xmath84 from which we obtain the rotational curves .
the energy density of the halo and the disk , the stress of the halo , the gravitational mass and the circular velocity are the physical properties that describe the disk - halo system .
all these properties are written in terms of the solution to laplace equation @xmath24 . nevertheless , as we saw in sec .
[ sec : emtensor ] , @xmath85 needs to satisfy the conditions ( [ eq15 ] ) and ( [ eq17u ] ) in order to define a surface energy density in the plane @xmath37 .    to find a particular family of models we choose @xmath24 to be the solution to laplace equation in spherical coordinates @xmath86 .
now , if we use the usual coordinates transformations @xmath87 and @xmath88 , then the condition ( [ eq17u ] ) will not be satisfied because @xmath57 . to solve this problem we use
the displacement , cut and reflect method used by kuzmin ( 1956 ) and toomre ( 1963 ) @xcite to generate disk - like sources in the frame of newtonian gravity .
the method is based on the transformation @xmath89 where @xmath90 is a constant . with this
, the function @xmath24 takes the form @xmath91 where @xmath92 , @xmath93 , @xmath94 are arbitrary constants and @xmath95 are the legendre polynomials .
the solution to the laplace equation presented in ( [ ukuz ] ) allows us to write the energy density of the halo as    @xmath96^{2 } + \left[\sum_{l=0}^{n}\left(a_{l}(l+1)p_{l+1}(\frac{|z|+a}{r})\right)/r^{l+2}\right]^{2}}{\left[1+\sum_{l=0}^{n}a_{l}p_{l}(\frac{a}{r})/r^{l+1}\right]^{2(k+1)/k } } \right\}. \label{rhokuzmin}\ ] ]    in a similar way , we can write the energy density of the disk as @xmath97/r_{0}^{l+2}}{\left[1+\sum_{l=0}^{n}a_{l}p_{l}(a / r_{0})/r_{0}^{l+1}\right]^{(k+1)/k } } , \label{sigmakuzmin}\ ] ] where @xmath98 .
so , from the behavior of @xmath95 and the fact that @xmath99 vanishes at infinity , we can see that @xmath100 and @xmath101 are always positive , have their maximum value at the center of the system ( @xmath102 , @xmath37 ) and both vanish at infinity .
now , we are going to analyze the behavior of @xmath103 and @xmath104 , and to do so we will consider only the models with @xmath105 .
the solutions that describe the energy density of the halo for the first three models are given by    @xmath106^{2}+[\tilde{z}(\tilde{a}_{0}\tilde{r}^{2}+3\tilde{a}_{1}\tilde{z})-\tilde{a}_{1}\tilde{r}^{2}]^2}{\tilde{r}^{2(2k-3)/k}[\tilde{r}^{3}+\tilde{r}^{2}\tilde{a}_{0}+\tilde{a}_{1}\tilde{z}]^{2(1+k)/k } } , \label{r1}\\ \nonumber\\ \tilde{\rho}_{2}&=&\frac{\tilde{r}^{2(5 - 2k)/k}\left[\tilde{a}_{0}\tilde{r}\tilde{r}^{4}+(3\tilde{a}_{1}\tilde{z}+\tilde{a}_{2})\tilde{r}\tilde{r}^{2}+\frac{5}{2}\tilde{a}_{2}(2\tilde{z}^{2}-\tilde{r}^{2})\right]^{2}}{\left[\tilde{r}^{5}+\tilde{a}_{0}\tilde{r}^{4}+\tilde{a}_{1}\tilde{r}^{2}\tilde{z}+\frac{1}{2}\tilde{a}_{2}(\tilde{z}^{2}-\tilde{r}^{2})\right]^{2(1+k)/k } } \nonumber \\ \nonumber \\ & & + \frac{\tilde{r}^{2(5 - 2k)/k}\left[(\tilde{a}_{0}\tilde{z}-\tilde{a}_{1})\tilde{r}^{4}+3\tilde{z}(\tilde{a}_{1}\tilde{z}-\tilde{a}_{2})\tilde{r}^{2}+\frac{5}{2}\tilde{a}_{2}\tilde{z}(2\tilde{z}^{2}-\tilde{r}^{2})\right]^{2}}{\left[\tilde{r}^{5}+\tilde{a}_{0}\tilde{r}^{4}+\tilde{a}_{1}\tilde{r}^{2}\tilde{z}+\frac{1}{2}\tilde{a}_{2}(\tilde{z}^{2}-\tilde{r}^{2})\right]^{2(1+k)/k } } , \label{r2}\end{aligned}\ ] ]     + where @xmath107 , @xmath108 , @xmath109 , @xmath110 , @xmath111 , @xmath112 , @xmath113 , and @xmath114\rho_{n}$ ] .    in fig .
[ denhalokuz ] we plot ( [ r0]-[r2 ] ) to analyze the behavior of the energy density of the halo .
this figure allows us to see that the most of the energy in the halo is located at the center of the system , i.e. in the region near to @xmath115 and @xmath116 .
the density presents a maximum at the center and , it tends to zero at infinity with a rate that depends on the constants values in each particular model
.    on the other hand , the energy density of the disk for the models @xmath105 is given by    @xmath117\tilde{r}_{0}^{\frac{3 - 2k}{k}}}{\left[\tilde{r}_{0}^{3}+\tilde{a}_{0}\tilde{r}_{0}^{2}+\tilde{a}_{1}\right]^{\frac{1+k}{k}}},\label{s1}\\ \nonumber \\ \tilde{\sigma}_{2}&=&\frac{\left[(\tilde{a}_{0}-\tilde{a}_{1})\tilde{r}_{0}^{4}+3(\tilde{a}_{1}-\tilde{a}_{2})\tilde{r}_{0}^{2}+\frac{5}{2}\tilde{a}_{2}(2-\tilde{r}^2)\right]\tilde{r}_{0}^{\frac{2(5 - 2k)}{k}}}{\left[2\tilde{r}_{0}^{5}+2\tilde{a}_{0}\tilde{r}_{0}^{4}+2\tilde{a}_{1}\tilde{r}_{0}^{2}+\tilde{a}_{2}(2-\tilde{r}^2)\right]^{\frac{1+k}{k}}},\label{s2}\end{aligned}\ ] ]    where @xmath118 , @xmath110 , @xmath111 , @xmath112 , @xmath113 , and @xmath119 .
we plot , in fig .
[ densidiskuz ] , the energy density @xmath120 for the first models with @xmath121 . in each graph
we plot some curves for different values of the constants @xmath122 and @xmath19 in order to analyze different behaviors .
hereafter , we introduce a convenient notation in which each curve is labeled by the set of numbers ( @xmath123 ) .
it is important to clarify that the constants are chosen so that the energy density is positive ; this assure the condition ( [ eq21 ] ) .
the surface energy density profiles , presented in fig .
[ densidiskuz ] , go to zero at infinity from a maximum value located at the center of the disk .
we can change the maximum value of @xmath124 and its decrease rate by changing the values of the constants . in the three models ,
the most of the energy is concentrated at the central region of the disk , but as we can see , when we increase @xmath125 , the energy is less distributed in the space , i.e. the energy density goes to zero faster .
we can compute the total mass for the disk - halo system , which is defined through ( [ mh ] ) and ( [ md ] ) , for any model characterized by the integer value of @xmath125 . by doing the integrals for @xmath126
we find that , @xmath127.\label{mh0}\end{aligned}\ ] ] now , since @xmath104 and @xmath103 vanish at infinity , it is important to show that @xmath128 will converge for all @xmath125 .
to do so , we start by putting ( [ md ] ) in the form @xmath129 where @xmath98 and @xmath130 with ( [ ukuz ] ) it is easy to show that @xmath131 and by using the limit comparison test @xcite we can demonstrate that @xmath132 will be finite if @xmath133 converges .
now , since @xmath134 has a finite value , then the mass of the disk for every model will be finite too .
to guarantee the convergence of @xmath135 we proceed in a similar manner : the integral ( [ mh ] ) , with @xmath136 can be rewritten as @xmath137 where @xmath138 and @xmath139 .
it is easy to show that @xmath140 so by virtue of the limit comparison test , if @xmath141 converge , then @xmath142 will converge too and from ( [ mhproof ] ) , @xmath135 and @xmath143 will be finite .
since @xmath144 has a finite value , then the mass of the halo for all models converge .
@xmath145    finally , we present the circular velocity for the models we are considering . to compute the expressions for this physical quantity we use ( [ ukuz ] ) and ( [ eq33 ] ) to obtain    @xmath146}{k[m^{2}+m\tilde{a}_{1}+m^{3}\tilde{a}_{0}+\frac{1}{2}\tilde{a}_{2}(3-m)]m^{3}-\tilde{r}^{2}[3\tilde{a}_{2}+3\tilde{a}_{1}m^{2}+\tilde{a}_{0}m^{4}+\frac{3}{2}\tilde{a}_{2}(3-m)m ] } , \label{v2}\end{aligned}\ ] ]     + where @xmath147 , @xmath110 , @xmath111 , @xmath112 , and @xmath113 .
now , we plot these expressions to see clearly the dependence of @xmath148 with the coordinate @xmath149 for a time - like particle in a circular orbit .
the resulting plots are presented in fig .
[ circvelkuz ] . here
we choose the same constants values that we use to perform fig .
[ densidiskuz ] in order to compare the rotational curve with its corresponding surface energy density profile .    in the plot for @xmath126 model
we see that the velocity increases very fast until a maximum value from which we distinguish clearly two kind of behavior : in the first one , the velocity decrease with @xmath149 , and in the second one it remains approximately constant . if we compare with the corresponding plot in fig .
[ densidiskuz ] , we see that the second kind of behavior is associated with energy density profiles whose maximum is higher than the cases where the velocity decrease after reaching its maximum .    in the plot for @xmath150 ,
all the curves show that the velocity remains constant after reaching a maximum , but , by comparing with fig .
[ densidiskuz ] , we see that if the maximum is higher , the energy density reaches a lower value in @xmath115 . in the plot for @xmath151
we obtain different kind of behavior .
of particular interest are those with two peaks because actually , the observational data shows that the velocity is not exactly constant after the first maximum , it has fluctuations qualitatively similar to those of our plots . finally , comparing with the plot for @xmath152 in fig .
[ densidiskuz ] , we see that , when the rotational curves present two peaks , the corresponding energy density profiles have two peaks as well .    @xmath153    @xmath154
in order to obtain a new family of finite thin disks with halo , we use oblate spheroidal coordinates ( @xmath155 ) to solve the laplace equation , because they adapt to the geometry of the source and introduce , in a natural way , a finite radius for the disk .
the transformation relations with the cylindrical coordinates are    @xmath156    where @xmath157 and @xmath158 .
the coordinates of the disk of radius @xmath159 are @xmath160 and @xmath161 , and if we cross the disk , @xmath162 changes sign .    in these coordinates ,
the solution to the laplace equation takes the form @xmath163 where @xmath164 are the legendre polynomials of order @xmath165 , @xmath166 , with @xmath167 the legendre function of second kind of imaginary argument @xcite , and @xmath168 are constants choosen to satisfy both the boundary condition ( [ eq14 ] ) and the condition for the allowed values for the speed of a particle in a circular trajectory , @xmath169 . according to ( [ eq15 ] ) and ( [ eq17u ] ) , to define properly a disk - like source ,
we want the metric tensor to be continuous and its first @xmath0-derivate to be discontinuous across the disk ; therefore , since @xmath170 is continuous across the disk and @xmath162 is discontinuous , we only need to consider even legendre polynomials .
therefore , @xmath24 in ( [ uobla ] ) is formed only by the even terms of the general solution to laplace equation .    to analyze the behavior of the physical quantities we have to consider particular solutions by restricting the sum in ( [ uobla ] ) .
one way to do it is to consider @xmath24 to be the gravitational newtonian potential of a finite disk , where the constants @xmath168 are found by imposing the orthogonality property of the legendre polynomials to give @xcite @xmath171 where , @xmath172 , \label{eq35}\ ] ] if , @xmath173 , and @xmath174 if @xmath175 , for @xmath176 . in the above equation @xmath177 and @xmath159
are the total mass and the radius of the disk respectively .
depending of the value of @xmath178 we will have a particular solution because the constants @xmath168 are restricted by @xmath179 ; so we are going to choose only the three first models ( @xmath180 ) to see the behavior of the energy densities of the disk and the halo , their masses , the state equation of the halo and the rotational curves .
the energy density of the halo for the first three models are computed by using ( [ uobla ] ) in ( [ eq10u ] ) in such a way that @xmath181 , \label{eq35}\ ] ] where    [ u ] @xmath182 , \\ u_{2 } & = - { \tilde m } \left [ \operatorname{acot}\xi + \frac{10 a}{7 } ( 3\eta^{2 } - 1 ) + b ( 35 \eta^{4 } - 30 \eta^{2 } + 3 ) \right ] ,   \\ u_{3 } & = - { \tilde m } \left [ \operatorname{acot}\xi + \frac{10 a}{6 } ( 3 \eta^{2 } - 1 )   + \frac{21 b}{11 } ( 35 \eta^{4 } - 30 \eta^{2 } + 3 )   + c ( 231 \eta^{6 } - 315 \eta^{4 } + 105 \eta^{2 } - 5 ) \right],\end{aligned}\ ] ]    with @xmath183 and    @xmath184 , \\ b & = \frac{3}{448 } \left [ ( 35 \xi^{4 } + 30 \xi^{2 } + 3 ) \operatorname{acot}\xi - 35 \xi^{3 } - \frac{55}{3 } \xi \right ] , \\ c & = \frac{5}{8448 } \left [ ( 231 \xi^{6 } + 315 \xi^{4 } + 105 \xi^{2 } + 5 ) \operatorname{acot}\xi   - 231 \xi^{5 } - 238 \xi^{3 } - \frac{231}{5 } \xi \right].\end{aligned}\ ] ]    note that @xmath185 and @xmath162 are related with the cylindrical coordinates by virtue of the transformation relations ( [ robla ] ) and ( [ zobla ] ) , so we can compute the first @xmath4-derivate of the above expressions ,    [ ur ] @xmath186 , \\
u_{2,r } = & \frac{15\tilde{m}\tilde{r}}{16d(1+\xi^{2 } ) } \left\ { \xi-13\eta^{2}\xi+(3 - 15\eta^{2})\xi^{3}+(1+\xi^{2 } ) [ 1 - 3\xi^{2}+3\eta^{2}(1 + 5\xi^{2 } ) ] \operatorname{acot}\xi\right\ } , \\
u_{3,r } = & \frac{35\tilde{m}\tilde{r}}{128d(1+\xi^{2})}\left\ { \xi \left [ 3 - 4 \xi^{2 } - 15 \xi^{4 } + \eta^{2 } ( 6 + 200 \xi^{2 } + 210 \xi^{4 } ) - \eta^{4}(113 + 420 \xi^{2 } + 315 \xi^{4 } ) \right ] \right . \notag \\ & \left .
+ 3 ( 1+\xi^{2 } ) \left [ 1 - 2 \xi^{2 } + 5 \xi^{4 } + \eta^{2}(2 - \xi^{2}-70\xi^{4})+5\eta^{4}(1 + 14\xi^{2}+21\xi^{4 } ) \right ] \operatorname{acot}\xi \right\},\end{aligned}\ ] ]    and also , the first @xmath0-derivative , which are    [ uz ] @xmath187 , \\
u_{2,z } = & \frac{5\tilde{m}\eta}{4d } \left\ { -9 \xi^{2 } + \eta^{2 } ( 4 + 15 \xi^{2})+ 3 \xi \left [ 1 + 3 \xi^{2 } - \eta^{2 } ( 3 + 5 \xi^{2 } ) \right ] \operatorname{acot}\xi \right\ } , \\
u_{3,z } = & -\frac{7\tilde{m}\eta}{64d } \left\ { 50 \eta^{2 } \xi^{2}(11 + 21 \xi^{2 } ) - 15 ( \xi^{2 } + 15 \xi^{4 } ) - \eta^{4 } ( 64 + 735 \xi^{2 } + 945 \xi^{4 } ) \right . \notag \\ & \left .
+ 15 \xi \left [ -1 + 6 \xi^{2 } + 15 \xi^{4 } - 2 \eta^{2}(3 + 30 \xi^{2 } + 35 \xi^{4 } ) + \eta^{4}(15 + 70 \xi^{2 } + 63 \xi^{4 } ) \right ] \operatorname{acot}\xi \right\}.\end{aligned}\ ] ]    @xmath188    in fig .
[ densihalokal ] , we show the plots of the rescaled energy density @xmath189 as a function of @xmath190 and @xmath191 .
we can observe two kinds of behavior , both of them present in all the particular solutions discussed in this paper , that can be obtained by choosing properly the constant values @xmath19
, @xmath192 .
the first behavior is the one presented in the plots for @xmath193 , where the energy density reaches its maximum at the center of the halo ( @xmath194 ) , and then goes to zero at infinity .
it is clear , in this case , that the energy is more concentrated in the region near the disk .
the second behavior is shown in the plot associated to the @xmath195 model .
here we see an unusual distribution of energy ( if we think , only , in the visible matter of the halo ) , where the maximum of density is located at @xmath196 and @xmath197 , i.e. near the edge of disk , and goes to zero at infinity .
we can also get analytic expressions for the energy density of the disk for the particular models considered in this section by using ( [ uobla ] ) in ( [ eq20s ] ) , so that    [ s ] @xmath198^{\frac{1+k}{k}}},\label{s1}\\   \tilde{\sigma}_{2}&=&\frac{20(1-\tilde{r}^{2})^{3/2}}{\left[1+\frac{15\pi}{128}\tilde{m}a_{2}(\tilde{r})\right]^{\frac{1+k}{k}}},\label{s2}\\   \tilde{\sigma}_{3}&=&\frac{28(1-\tilde{r}^{2})^{5/2}}{\left[1+\frac{35\pi}{512}\tilde{m}a_{3}(\tilde{r})\right]^{\frac{1+k}{k}}},\label{s3}\end{aligned}\ ] ]     + where @xmath199 , @xmath190 , @xmath200 , and    @xmath201    in fig .
[ dendiskal ] , we plot the energy density @xmath120 of the disk as a function of the normalized radius for the models @xmath180 . in the first plot
we compare the energy density for the constants values @xmath202 and @xmath203 , and for different values of @xmath178 . in this plot
we see that @xmath120 has a maximum at the center of the system @xmath204 , and goes to zero at the edge ( @xmath205 ) .
the spatial distribution of the energy depends on the particular value of @xmath178 : for @xmath195 , we see that the energy density remains approximately constant , and starts to fall from @xmath206 ; but if we increase the value of @xmath178 , then the energy density will be more concentrated at the center , and its maximum value will be higher .
the remaining plots show how , by changing @xmath19 and @xmath192 , we can get different kinds of behavior in each particular model .
it is clear that we can manipulate the maximum value of the energy density and its rate of decline by varying the constants , but in all cases the behavior is very similar and the energy is concentrated in the center of the disk .    now , with ( [ eq8 ] ) we can write the gravitational mass ( [ mh ] ) equivalently as @xmath207 where @xmath208 is the area of a closed surface which contains a mass @xmath177 . we choose @xmath208 to be the surface of an ellipsoid determined by @xmath185 constant .
accordingly , we can define the mass of the disk throw the integral @xmath209 and the total mass of the system by @xmath210 so that , the mass of the halo is @xmath211 .
+    @xmath212    despite the fact that in this model the halo has infinite extension , it is possible to compute the total mass of the system through the integral given on ( [ mtobla ] ) , and show that @xmath128 converges . to do it ,
we first take advantage of the fact that at infinity , @xmath24 tends to zero .
so , by using ( [ uobla ] ) , we find that @xmath213 now , with the recurrence relation , @xmath214 we can transform ( [ eq25 - 1 ] ) on @xmath215 as @xmath216 becomes zero at infinity , and @xmath217 goes to @xmath218 if @xmath219 for all @xmath220 , then @xmath221 , and the total mass of the disk - halo system becomes @xmath222 which depends on the radius of the disk , the first constant @xmath223 and the particular kind of fluid in the halo ( determined by the value of @xmath19 ) . with this
, we proof that the total mass is finite , even though the halo has infinite extension .
this proof is very general , so @xmath128 is the same for all particular models .
now , for the first model ( @xmath195 ) the integral ( [ mdobla ] ) yield to @xmath224 , \label{eq36}\ ] ] where @xmath225 ; so the mass of the relativistic disk is related to the newtonian mass .
the analytic expression for the circular velocity of a particle on the plane @xmath37 has two parts : one part correspond to the particle moving on the disk ( @xmath226 ) , for which we need to do @xmath160 on ( [ eq33 ] ) ; and the other part correspond to the velocity of a particle moving outside the disk ( @xmath227 ) , for which we have to do @xmath228 on the expression for @xmath229 .
both parts must coincide exactly at @xmath205 .
now , by using the expressions for @xmath230 and @xmath231 in ( [ eq33 ] ) , and doing the appropriate evaluation in each region we get the solutions for the circular velocity for the models we are considering , which are given by    @xmath232 - 6 \tilde{m } \pi \tilde{r}^{2 } } ,   & 0
\leq \tilde{r } \leq 1 , \\               \\
\dfrac{6 \tilde{m } ( \tilde{r}^{2 } \operatorname{acot}j - j ) } { 6 \tilde{m } j + k ( 4 + 3 \tilde{m } j ) - 3 \tilde{m } [ 2 \tilde{r}^{2 } + k ( \tilde{r}^{2 } - 2 ) ] \operatorname{acot}j } , & \tilde{r } \geq 1 , \end{array } \right .
\\ \nonumber \\ v_{c_{2}}^{2 } & = \left\ { \begin{array}{ll } \dfrac{- 60 \tilde{m } \pi \tilde{r}^{2 } l_{1}}{60 \tilde{m } \pi \tilde{r}^{2}l_{1 } + k ( 128 + 15 \tilde{m } \pi l_{2 } ) } , & 0 \leq \tilde{r } \leq 1 , \\
\\ \dfrac{60 \tilde{m } [ \tilde{r}^{2}l_{1 } \operatorname{acot}j - j l_{3 } ] } { 60\tilde{m } j l_{3 } + k [ 45\tilde{m}(\tilde{r}^{2 } - 2 ) j - 64 ] - 15 \tilde{m } ( 4 \tilde{r}^{2 } l_{1 } + k l_{2 } ) \operatorname{acot}j } , & \tilde{r}\geq 1 , \end{array } \right .
\\ \nonumber \\ v_{c_{3}}^{2 } & = \left\ { \begin{array}{ll } \dfrac{-210 \tilde{m } \pi \tilde{r}^{2 } n_{1}}{210 \tilde{m } \pi \tilde{r}^{2}n_{1 } + k ( 35\tilde{m}\pi n_{2 } - 512 ) } , &   0\leq \tilde{r}\leq 1 , \\               \\
\dfrac{210 \tilde{m } [ 3 \tilde{r}^{2 } n_{1 } \operatorname{acot}j - j n_{3 } ] } { 210 \tilde{m } j n_{3 } + k ( 768 + 35 \tilde{m } j n_{4 } ) - 105 \tilde{m } ( 6 \tilde{r}^{2}n_{1 } + k n_{2 } ) \operatorname{acot}j } , &   \tilde{r}\geq 1 , \end{array } \right.\end{aligned}\ ] ]    where @xmath233 , @xmath234 , @xmath235 , @xmath236 , @xmath237 , @xmath238 , @xmath239 and @xmath240 .    in fig .
[ velcirckal ] , we show the profiles of the circular velocity as a function of @xmath149 for some constants values @xmath241 .
all the constants values , in this paper , are chosen in such a way that the physical quantities are well behaved ( @xmath242 and @xmath243 ) . in the first plot , in fig .
[ velcirckal ] , we compare the rotational curve in the models @xmath180 for @xmath19 and @xmath192 fixed . here
we can see that if @xmath178 is higher , then the peak would be less acute , the union between the two curves mentioned in the last paragraph would be smoother and the maximum value of @xmath72 would be attained closer to the center of the disk .
+ in all particular solutions , the behavior of @xmath72 is the same : the velocity starts from zero and reaches a maximum that depends on the particular constants values , then the velocity falls to a value that remains about constant ; this is very important because many real rotational curves from spiral galaxies @xcite has the same behavior ( qualitatively speaking ) . an specific shape for the rotational curve can be obtained by choosing properly the values for @xmath244 and @xmath192 , e.g , in the plots for @xmath180 , the curves associated to the values @xmath245 and @xmath246 show that we can also get circular velocities that fall after reaching the maximum
.    @xmath247
we obtained two sets of new exact solutions to the einstein equations for an axially symmetric static spacetime from the conformastatic metric .
these solutions describe the physical properties of an infinite disk of dust , i.e. without stresses , which is surrounded by a fluid or halo with an equation of state @xmath248^{-1}\rho$ ] .
this equation implies different kind of fluids : radiation if @xmath249 , dust if @xmath19 tends to infinity or another kind of fluid according to the value of @xmath19 .
the disk - like solutions obtained here are physically well behaved because they satisfy the energy conditions required for any source of gravitational field , and because they are asymptotically flat .
the asymptotic behavior was ensured by writing an appropriate relation between the metric function @xmath3 and the auxiliary function @xmath24 .    according to the solutions , we got for the halo a fluid whose main stresses are different from zero ( @xmath250 ) and they are in the same direction of the main axes . on the contrary ,
the solutions obtained for the disk describe a fluid with null stress . in the kuzmin - toomre disks with halo we obtained that the most of the energy density was located at the center of the system
; so it is possible to see a well - defined central region with a maximum at ( @xmath115 and @xmath116 ) ; nevertheless , as far as @xmath125 increases , the amount of constants increases too , so we have more free parameters to adjust the energy density profile to a particular behavior . for instance , we presented in fig .
[ densidiskuz ] one case in which the maximum energy density is not at the center of the system .
we found that the morgan - morgan relativistic thin disks models presented in section [ sec : morgan ] does not have the same freedom as the kuzmin - toomre ones . in these models
we only have two parameters to adjust : @xmath19 and @xmath192 , so by varing these parameters we only change the maximum of @xmath53 , @xmath12 and @xmath72 , and their corresponding rates of decrease .
we computed the total mass for the @xmath126 model of the kuzmin - toomre relativistic disks with halo , turning out that @xmath251 , and we proved that the mass of the system converges for all @xmath125 . on the other hand , in the morgan - morgan relativistic disks solutions we computed the total gravitational mass for all the family of solutions which is @xmath252 .
finally , we presented the rotational curves for some particular models in both families of solutions .
we saw that , for some values of the constants , it is possible to get a behavior in which the square of the circular velocity increases from @xmath115 until a radius where @xmath72 stays approximately constant . nevertheless , we got some rotational curves with two peaks in the model @xmath253 of the kuzmin - toomre relativistic solutions .
this is a very interesting behavior because the observed rotational curves in spiral galaxies exhibit these kind of features @xcite .
gag was supported in part by vie - uis , under grants number 1347 and 1838 , and colciencias , colombia , under grant number 8840 .
omp wants to thank the financial support from the program _
jvenes investigadores _ of colciencias , colombia .