modified gravity ( for a review , see e.g.  @xcite ) constitutes an interesting dynamical alternative to @xmath0cdm cosmology in that it is also able to describe with success the current acceleration in the expansion of our universe , the present dark energy epoch .
the specific class of modified @xmath1 gravities ( for a review , see e.g.  @xcite ) has undergone many studies which suggest that this family of gravitational alternatives for dark energy is able to pass the constringent solar system tests . the investigation of cosmic acceleration as well as the study of the cosmological properties of @xmath1 models has been done in refs .
the possibility of a natural unification of early - time inflation with late - time acceleration becomes a realistic and quite natural possibility in such models , as is demonstrated e.g. in ref .
@xcite .
recently , the importance of modified gravity models of this kind has been reassessed with the appearance of the so - called ` viable ' @xmath1 models @xcite .
those are theories which satisfy both the cosmological as well as the local gravity constraints , which had caused in the past a number of serious problems to some of the first - generation theories , that had to be considered now as only approximate descriptions from more realistic theories .
the final aim of all those phenomenological models is to describe a segment as large as possible of the entire history of our universe , as well as to recover all local predictions of einstein s gravity that have been already verified experimentally , with very good accuracy , at the solar system scale .
it is remarkable that , as was demonstrated in refs .
@xcite , some of these realistic models lead to a natural unification of the early - time inflation epoch with the late - time acceleration stage .
let us recall that , in general ( see e.g. @xcite , for a review ) , the total action for the modified @xmath1 gravitational models can be written as [ xxx7 ] s = d^4 x + s_(m ) . here
@xmath2 is a suitable function of the scalar curvature @xmath3 , which defines the modified gravitational part of the model .
the general equation of motion in @xmath4 gravity with matter is given by [ xxx22 ] g _
f(r ) - r _ f(r ) - g _ f(r ) + _ _ f(r ) = - t_(m )  , where @xmath5 is the matter energy - momentum tensor and @xmath6 is the derivative of @xmath1 with respect to its argument @xmath3 . for a generic @xmath1 model
is not easy to find exact static solutions .
however , if one impose some restrictions , one can proceed along the following lines .
first , we may require the existence of solutions with _ constant _ scalar curvature @xmath7 , and we arrive at f(r_0)r_=g_. [ geneqrc ] taking the trace , we have the condition 2f(r_0)=r_0f(r_0 ) [ b ] and this means that the solutions are einstein s spaces , namely they have to satisfy the equation r_=g_= g _ , [ e ] @xmath8 being a solution of eq .
( [ b ] ) .
this gives rise to an effective cosmological constant , namely _ eff==. [ l ] the purpose of our work will be to study the appearance of multiply de sitter space solutions in several realistic models of modified gravity .
the occurrence of multiply de sitter solutions plays a fundamental role in modified gravity because it permits to describe the inflation stage as well as current @xmath0cdm cosmology in terms of this theory alone without any need for either fine - tuning of a cosmological constant nor of introducing extra scalar fields . in summary , this is a minimal and at the same time very powerful approach which circumvents some of the hardest problems of present day physics .
one may argue that these theories are equivalent to introducing extra scalar fields , but this equivalence has been proven to hold at the classical level only , not at the quantum one .
in addition , the cosmological interpretation of modified gravity solutions is different from that of scalar field cosmology .
the paper is organized as follows . in the next section
we discuss a viable modified gravity model @xcite in an attempt to study de sitter solutions there .
it is shown that , for some values of the parameters , it is possible to find several de sitter spaces ( as a rule , numerically , for small values of the curvature exponent ) ; one of them may serve for the inflationary stage , while the other one can be used for the description of the dark energy stage .
the evolution of the effective equation of state parameter is investigated numerically .
3 is devoted to the analysis of the same question in a slightly generalized model which is known to describe the unification of the early - time inflation with the late - time acceleration epochs @xcite . the same numerical investigation is carried out , and a multiply de sitter universe solution is constructed . in section 4 , the corresponding problem is investigated for a viable model of tangential modified gravity , proposed in ref .
@xcite . having in mind the possibility to use a sds universe for the description of cosmic acceleration , those solutions for modified gravity are investigated in sect .
comparison of the free energies of the de sitter and the sds solutions , for the viable model of the second section , is then made .
finally , some outlook is given , together with the conclusions , at the end of the paper .
in this section we will study de sitter solutions in realistic modified gravity models using a dynamical system approach .
we shall start from the following form for the initial action @xmath9 ,   \label{1.1}\ ] ] where @xmath10 is some function of the geometrical invariants . in order to make use of the method of dynamical systems , we work in the metric corresponding to a spatially - flat frw universe , namely @xmath11 from here
, the frw equation can be written in the following way ( for more details see ref .
@xcite ) @xmath12 where @xmath13 the `` dot '' over the symbol means derivative with respect to time @xmath14 , while @xmath15 and @xmath16 represent the derivative of @xmath10 with respect to @xmath17 and @xmath18 respectively . in what follows
we will mainly concentrate on the simplest case @xmath19 .
let us consider the following choice for the function @xmath20 @xcite , which represents a very interesting subclass of viable modified gravitational models @xmath21 where the constant @xmath22 has dimension of mass , while @xmath23 , @xmath24 , @xmath25 are some positive dimensionless constants .
note that here @xmath26 $ ] , @xmath27 , @xmath28 . after making variation ,
it is chosen @xmath29 .
introducing a nonzero constant @xmath25 is here equivalent to introducing a shift in the effective cosmological constant , and for this reason we will assume @xmath30 in our further computations . finally , we can write eq .
( [ 1.3 ] ) for the function in ( [ 1.5 ] ) as @xmath31 ,   \end{array }   \label{1.6}\ ] ] where @xmath32 .
this equation can be rewritten as a dynamical system , namely @xmath33 it is easy to see that the critical points of this system are the de sitter points ( @xmath34 ) . to investigate the nature of these points we need to determine them explicitly , a non - trivial problem in the general case .
it is however easy to obtain from ( [ 1.6 ] ) the equation satisfied by the critical points @xmath35 : @xmath36 the same result can be obtained directly starting from equation ( [ b ] ) of the previous section .
it is convenient to introduce the notation @xmath37 , for further investigation of this equation .
so , finally , we have @xmath38 first of all we find one ( trivial ) root of this equation @xmath39 ( it corresponds to @xmath40 ) , what allows us to reduce the order of the equation .
( note also that , if @xmath41 , eq .
( [ 1.9 ] ) takes a more complicate form and then it does not have the trivial root @xmath39 . ) but nevertheless the equation still is of @xmath42-order and in the interesting case ( @xmath43 ) ; this is too high and the roots can not be found algebraically .
we need a specific discussion of this problem .
it is clear that if there are too many ( say 10 or 20 ) de sitter points in the theory , it looks like a classical analogue ( at a reduced scale , of course ) @xcite of the string landscape vacuum structure in which case it will be far from trivial to obtain the standard cosmology .
nevertheless , something can be done even in this case , by comparing the energies of the corresponding de sitter solutions , which should in fact differ , as is discussed in sect .  v. using descartes rule of signs
, we find that eq .
( [ 1.9 ] ) can have 2 or 0 positive roots .
more detailed information can be obtained by using sturm s theorem .
unfortunately it is not possible to apply it in the general case ( for arbitrary @xmath44 ) , so we need to investigate separately each of the different cases , @xmath45 , and see which one is more interesting from a cosmological viewpoint .
we have therefore undertaken here a systematic analysis of all possible cases in this region of values of @xmath44 .
computations are rather involved ( specially for larger @xmath44 ) , but the final results are not so difficult to describe analytically .
we have found , for any of the values of @xmath44 considered , in the range above , that the number of roots depends on the parameter @xmath46 only , and that there is an @xmath47 such that , for @xmath48 , eq .  (
[ 1.9 ] ) has no positive roots , while for @xmath49 , eq .  ( [ 1.9 ] ) has two positive roots .
we here now enumerate our results systematically : @xmath50 ; @xmath51 ; @xmath52 ; @xmath53 ; @xmath54 ; @xmath55 ; @xmath56 ; @xmath57 . note also that , near the critical point @xmath47 , both solution are very close to each other when they exist ( a real value of @xmath49 ) , and they are complex conjugate to each other for @xmath48 .
this means that , in a situation where these two roots differ substantially , it must necessarily be @xmath58 .
the ds - points described above can be used for the construction of inflationary or late - time acceleration behavior ( depending on the value of the scale factor @xmath59 ) . to this aim
, the corresponding ds point must be unstable in the inflationary stage but can be either unstable or stable for the late - time acceleration one .
however , this is in fact not a strictly necessary condition , since even for stable ds inflation , the exit from it can be achieved by a coupling with matter , through the effect of small non - local term or by some other mechanism .
unfortunately , an all - round investigation of stability of the ds point turns out to be very difficult in the general case .
for the recent analysis of critical points in more general modified gravity theories depending on all geometrical invariant see ref .
@xcite .    here
we can carry out the analysis on the stability of de sitter points for the model we are dealing with in this section , only for some chosen numerical set of parameters .
such computations show that one of the ds points is very likely to be always unstable ( with a smaller value of @xmath60 ) and that another one is very likely to be always stable ( with larger value of @xmath60 ) .
this means that there is actually an easy way for the inflationary stage construction ( if the initial conditions lie sufficiently close to the unstable ds point , see fig .  [ fig.1 ] ) .
however , only within the model under consideration it might be a problem with late - time acceleration , originated from the stable point .
the only possibility for late time acceleration is to choose very special parameters , in which the stable point is situated very far away from the unstable one .
of course , the problem disappears when one describes only late - time acceleration within such a model , or when one takes into account other terms like the local ones .
note however that specific values of the parameters in which both points are unstable can be chosen too , which give rise to a sound theory which provides a unified description of the inflationary and late - time acceleration epochs .    ) near the equilibrium points .
stars denote the beginning of the trajectories . dashed and dot
dashed lines lead to bottom.,scaledwidth=50.0% ]     near an unstable point for the dashed and dot - dashed lines of fig .
[ fig.1 ] .
the time origin corresponds to the stars from the previous graph.,scaledwidth=50.0% ]    let us now consider specific numerical results .
first , note that a full numerical investigation starting from the inflation ( or from our ds point ) and extending to the late - time acceleration epoch and going through a frw - like stage is not feasible because of the presence of numerical instabilities that prevent such possibility .
for this reason , we will investigate the regions near the equilibrium points only , and especially near the unstable one . in fig .
[ fig.1 ] the trajectories in the phase plane ( @xmath61 ) are presented for the set of parameters : @xmath62 , @xmath63 , @xmath64 , @xmath65 .
such small values of @xmath23 and @xmath24 are motivated in @xcite , but the picture is typical in any case .
the position of the ds point is specified .
we may see one stable point ( with bigger value of @xmath60 ) and another , unstable one .
the stars in the graph correspond to the beginning of the trajectories .
we see clearly that trajectories can pass sufficiently close to the ds point , and that they escape from it , lead by the two attractors .
one of them ( the upper one ) leads to a singularity , while another one ( heading towards the bottom ) , leads to a solution going somewhere in the frw - like region ( this is difficult to analyze exactly , due to numerical instabilities ) . in fig .
[ fig.2 ] the evolution of @xmath66 is depicted for two trajectories : the upper one corresponds to the dashed line of fig .
[ fig.1 ] , the bottom one to the dot - dashed line of fig .
[ fig.1 ] ( both of them lead to the bottom in fig .
[ fig.1 ] ) .
we see that there actually are trajectories which start from phantom - like solutions and tend to a normal evolution ( when this is possible in the model ) . of course , these are qualitative considerations only , owing to the fact that no matter is taken into account .
most probably , inclusion of matter will change the situation .
however , it is just remarkable that both the early - time inflation epoch as well as the late - time accelerating one could be obtained in a unified and natural way in such a model , owing to the presence of several ds points .
a slightly different choice for the function @xmath2 is motivated by the realistic and viable model which was proposed for the unification of early - time inflation and late - time acceleration in ref .
@xcite , namely @xmath67}.   \label{2.1}\ ] ] here we imply that @xmath68 , @xmath69 and @xmath8 is positive , what seems reasonable .
now , our frw - like equation ( [ 1.3 ] ) reads @xmath70 ,   \end{array }   \label{2.2}\ ] ] where @xmath32 . rewriting this equation under the form of a dynamical system ( [ 1.7 ] )
we find , by introducing the new variable @xmath71 , that the equilibrium points of this system ( which are ds points indeed ) are obtained from the following equation @xmath72x^m+2\kappa^2r_0f_0mx^{m-1}+\\ \\
( f_0+f_1r_0^m)^2x+r_0(f_0+f_1r_0^m ) ( f_0+f_1r_0^m-4\kappa^2r_0^{m-1})=0 .
\end{array }   \label{2.3}\ ] ] it is clear that , in the general case ( even for fixed @xmath73 ) , solving this equation will not be an easy thing . we can make use of the following hint .
let us require that @xmath74 be a solution of eq .
( [ 2.3 ] ) .
this is possible only if there is some relation among the constants of our theory . as we can easily see from ( [ 2.3 ] ) , it must be that @xmath75 .
moreover , since @xmath76 , we need that @xmath77 . substituting this condition into ( [ 2.3 ] ) we find an essential simplification of the former equation , namely @xmath78 where the root @xmath74 is already excluded .
we now discuss the physical meaning of finding a solution @xmath74 . as a rule
this will imply that @xmath8 is the value of the scalar curvature at present time , but in principle @xmath8 is just a parameter of the theory , which could be given any likely value .
we can , for instance , fix @xmath8 to have the ds point value which is situated exactly at present time , @xmath79 , or at @xmath80 years .
and there is also the possibility to explain late time acceleration if this point is stable .
now let us consider other possible roots of eq .
( [ 2.4 ] ) . using sturm s theorem as in the previous section ,
we find several positive roots ( which correspond to ds points in the past ) and a number of negative roots ( needless to say , only values which are @xmath81 have a physical meaning ) , corresponding to ds points in the future .
unfortunately , computations are harder than in the previous case and our calculations have been performed for three specific values of the parameter @xmath82 : 3 , 5 and 7 , only . in any case , the results are very similar for the three situations , what gives us hope that for larger values of @xmath82 result will be also similar .
we may predict the upper limit of the ds points in the future only , which is 2 , but in any case this point is not much interesting .    concerning the ds points in the past
, the picture turns out to be much like the one in the previous case : the number of ds points depends only on a dimensionless parameter for all investigated values of @xmath82 .
that is , @xmath83 . and
there are two values , @xmath84 and @xmath85 , which are different for each @xmath82 , so for @xmath86 and @xmath87 there are two ds points and for @xmath88 there is no ds point .
we have not been able to obtain exact analytical expressions for @xmath84 and @xmath85 and thus give here numerical results only . for @xmath89 @xmath90 , @xmath91 ( as positive roots of @xmath92 ) . for @xmath93 @xmath94 , @xmath95 ( as positive roots of @xmath96 ) . for @xmath97 @xmath98 , @xmath99 ( as positive roots of @xmath100 ) .
as we can see , the interval where there are no ds points shrinks when the value of the parameter @xmath82 is increased .
if we assume that for larger @xmath82 the number of ds points depends on the parameter @xmath101 only  as clearly happens in the investigated cases  we do find @xmath84 and @xmath85 for any value of @xmath82 .
such investigations show that @xmath84 and @xmath85 slowly change when @xmath82 increases and that the interval without ds points is present even for very big value of the parameter @xmath82 .
for example , for @xmath102 we have @xmath103 and @xmath104 .
a simple numerical investigation of stability for the existing ds points shows that one of them ( the one with the smaller value of @xmath60 ) is always unstable , but the other one ( with the larger value of @xmath60 ) can either be stable or unstable . from our numerical results ,
the appearance of unstable ds points pairs relative to the interval @xmath86 seems most probable .
note also that the condition @xmath76 means that it actually must be @xmath105 and that in this case it is most likely the one stable and one unstable point show up .
now let us consider the evolution equation ( [ 2.2 ] ) from a different point of view .
as we already noted , @xmath8 is a parameter of the theory which corresponds to the value of the scalar curvature at some epoch .
this means that during the normal evolution of our universe , from large to little ( or zero ) curvature , it becomes @xmath7 at some moment .
but , as we can see from eq .
( [ 2.2 ] ) , this means that that at this point the coefficient of the higher derivative term ( @xmath106 ) is equal to zero .
this is a well - known mathematical problem , which needs special investigation . as we know from mathematics there are two possibilities : the solution of the perturbed equation ( with a higher derivative term ) may tend to the solution of the degenerate equation ( without higher derivatives ) , and then the coefficient of the higher derivative term may tend to zero or not tend to the solution of the degenerate equation .
a special investigation of this problem shows that solutions of eq .
( [ 2.2 ] ) tend to solutions of the degenerate equation , which is @xmath107 when @xmath108 but not so when @xmath109 .
this means that the point @xmath8 is reachable during the evolution from large @xmath3 to zero , but it is not the final point of the evolution , because there is an instability in the future . on the other hand
, we have @xmath7 at this point , and therefore @xmath110 . substituting this into ( [ 2.5 ] ) we find a relation among the parameters of the theory which were introduced before , @xmath75 .
this means that , generically , if one wants the point @xmath8 to be reachable during evolution , a relation of this sort among the parameters of the theory must be fulfilled .
that is , we have two independent parameters only ( or even just one if we consider that @xmath8 is strictly related to our epoch ) . note also that it is impossible to use standard numerical methods near this point for solving eq .
( [ 2.2 ] ) because all available methods need solving the equation with respect to the highest derivative , and it turns out that , near this point , the numerical solution is unstable .    )
near equilibrium points .
stars denote beginning of trajectories . dashed and dot dashed lines lead to bottom.,scaledwidth=50.0% ]     near unstable point for dashed and dot - dashed lines of fig.[fig.1 ] .
zero in time corresponds to the stars from the previous graph.,scaledwidth=50.0% ]    now let us consider numerical results near the stable points .
since the value of @xmath8 is very small in our epoch ( @xmath111 @xmath112 ) , it is impossible to use real values in computations . thus , we have produced a qualitative analysis using the following set of parameters : @xmath113 , @xmath114 , @xmath115 , @xmath97 .
as we can see , @xmath116 , and there must be two ds points corresponding to positive @xmath117 .
the evolution lines near these points are represented in fig .
[ fig.3 ] , where we can trace one stable and one unstable point , stars denoting the beginning of trajectories , as in the previous case . actually , there are three ds points for these chosen parameters
. one of them , @xmath118 , which corresponds to @xmath74 , is degenerate and situated out of the graph .
the time evolution of @xmath119 is represented in fig .
[ fig.4 ] by the dashed and dot - dashed lines .
the picture is very similar to the one in the previous section , and the same comments also apply here .
thus , we have here shown that the unification of early - time inflation with late - time acceleration is in principle possible , and even very likely , due to the appearance of several ds points in the evolution of the universe .
for a more realistic study , the presence of ( minimally and non - minimally coupled ) matter should be taken into account with care .
as a third example , now we discuss a model proposed in @xcite .
such model is defined by means of the function [ tan1 ] f(r)=r - f(r)=r - a , where @xmath120 , @xmath121 and @xmath8 are arbitrary parameters .
one immediately sees that @xmath122 , as required and , moreover , that [ tan3 ] _
rf(r)=2_eff a. if @xmath123 in the present universe then @xmath124 plays the role of the effective cosmological constant .
we also observe that the derivative [ tan4 ] f(r)=1- has a minimum when @xmath7 , which reads [ tan5 ] f(r_0)=1- . in order to avoid antigravity
one needs to require [ tan6 ] 0<f(r)<f(r_0)= 1- . the model given by eq .
( [ tan1 ] ) is able to describe late acceleration , since , in general , de sitter critical points exist .
they are the solutions of ( [ b ] ) .
then , we set k(r)=2f(r)-rf(r ) , [ rkr ] and study ( numerically ) the zeros of this transcendental function ( see fig .  [
p : tanh1 ] ) .
we see that for suitable choices of the parameters there are one or two de sitter critical points .     for @xmath125 and different choices of the parameter @xmath121
: @xmath126 , @xmath127 ( picture on the left ) ; @xmath128 , @xmath129 ( picture on the right).,title="fig:",scaledwidth=40.0% ]   for @xmath125 and different choices of the parameter @xmath121 : @xmath126 , @xmath127 ( picture on the left ) ; @xmath128 , @xmath129 ( picture on the right).,title="fig:",scaledwidth=40.0% ]    now , we impose @xmath8 to be a de sitter critical point for the model in eq .
( [ tan1 ] ) .
this means that ( [ rkr ] ) has to be satisfied for @xmath7 and such a condition fixes one parameter , say @xmath120 .
then , we get a= , and introducing , for convenience , the dimensionless variables x=,b_0=r_0b , we obtain f(r_0)=1- . _ eff= , k(r)=- . in order to have @xmath130 and , at the same time , to avoid antigravity
, the parameter @xmath131 has to be negative and in the range @xmath132 . by varying @xmath131 in that range
, the value of @xmath124 can acquire any desired value .
moreover the model has always two de sitter critical points with constant curvatures @xmath133 and @xmath134 .
one can see that @xmath135 , while @xmath136 . in fig .
( [ p : tanh2 ] ) we have plotted @xmath137 for some values of @xmath131 .    , with @xmath8 as a fixed de sitter solution , for different choices of the parameter @xmath131 : @xmath138 , @xmath139 ( picture on the left ) ; @xmath140 , @xmath141 ( picture on the right).,title="fig:",scaledwidth=40.0% ] , with @xmath8 as a fixed de sitter solution , for different choices of the parameter @xmath131 : @xmath138 , @xmath139 ( picture on the left ) ; @xmath140 , @xmath141 ( picture on the right).,title="fig:",scaledwidth=40.0% ]    thus , our study shows the existence of multiply - de sitter universes in more complicated , tangential models .
one can prove that other viable @xmath2 gravities introduced in ref .
@xcite naturally lead to de sitter solutions too , as demonstrated in ref .
in this section we will discuss spherically symmetric exact solutions for the modified gravities above .
we mainly concentrate on multiply black hole solutions and on their related entropies and free energies .
we have seen in the introduction that , requiring constant curvature solutions , one may consider the simplified form of the equations of motion given in ( [ e ] ) . as a result
, @xmath1 modified gravity models admit the following class of general static neutral black hole solutions  in four dimensions and with a non vanishing cosmological constant . in order to describe them ,
recall the metric ds^2=-v(r)dt^2++r^2d_k^2,v(r)>0 , [ ds ] where @xmath142 and the possible horizon manifolds are @xmath143 , the two dimensional sphere , @xmath144 , the two dimensional torus , and @xmath145 , the two dimensional compact riemann surface .
the scalar curvature for the _ ansatz _ ( [ ds ] ) reads r=- [ r^2v+4rv+2v-2k ] and this means that every constant curvature solution with @xmath7 has to satisfy the equation r^2v+4rv+2v-2k =- r^2 r_0 .
the general solution of this differential equation depends on two integration constants , @xmath121 and @xmath146 , and reads v(r)=+k -- r^2 .
the ricci tensor and the scalar curvature have to satisfy the equations
r_=g_,r = r_0 , [ e1 ] where @xmath8 is implicitly given by ( [ b ] ) , that is r_0= , now , it is easy to verify that eqs .
( [ e1 ] ) are satisfied only if @xmath147 , while @xmath146 is an arbitrary parameter , which is usually assumed to be non - negative and related to the mass @xmath148 of the black hole by @xmath149 .
the special case @xmath150 is also admissible .
then , we have v(r)=k -- r^2,c0,v(r)>0 .
[ vr ] in principle there are physical solutions of the latter equation for @xmath151 and @xmath152 , which give rise , respectively , to a torus topology and a hyperbolic topology for the horizon manifold ( the so called topological black holes @xcite ) . here
we are mainly interested in the usual spherical symmetric horizons and so we only consider in detail the case @xmath153 .
as it is well known , in the special case @xmath154 , @xmath155 in ( [ vr ] ) is always positive when @xmath156 , and this corresponds to the anti de sitter ( ads ) solution . on the contrary , when @xmath157 , @xmath155 is positive , for @xmath158 , and this corresponds to the de sitter solution . if @xmath153 and @xmath159 , one has black hole solutions but only if c^2r_0 - 032mg1 .
[ bkc ] we see that ( [ bkc ] ) is always satisfied if @xmath156 .
this corresponds to the schwarzschild - anti - de sitter ( sads ) black hole . in this case
@xmath160 , @xmath161 being the positive root of @xmath162 ( horizon radius ) . if @xmath157 , there are solutions only if @xmath163 . in this case
@xmath164 , @xmath161 and @xmath165 being the positive roots of @xmath162 ( horizon and cosmological radius respectively ) .
the extremal case @xmath166 is also admissible ( nariai solution ) , but the thermodynamics of such a black hole have to be discussed separately @xcite .    at this point , we provide a brief discussion regarding the thermodynamical properties of the above black hole solutions . if one make use of the noether charge method for evaluating the entropy associated with the black hole solutions with constant curvature @xcite in modified @xmath1 gravity models , one has @xcite
s= f(r_h ) .
[ ve ] where the factor @xmath1 has to be evaluated on the horizon , with area @xmath167 . as a consequence ,
one obtains a modification of the `` area law '' .
several examples have been discussed in @xcite . in the case of a constant curvature solution one
has simply @xmath168 , @xmath8 being the solution of eq .
( [ b ] ) , which has been investigated in previous sections .
with regard to this , we have found sufficient conditions to have two de sitter solutions .
thus , we may investigate their thermodynamical behavior evaluating the associated free energy .
the free energy @xmath169 is related to the canonical partition function @xmath170 by = - . on the other hand ,
a semiclassical approximation gives ze^-i_e , where @xmath171 is the euclidean classical action associated with the de sitter solution .
a direct calculation leads to i_e=- and from ( [ b ] ) and ( [ ve ] ) it directly follows that i_e =- s_h=- , [ fe ] @xmath172 and @xmath173 being , respectively , the entropy and the inverse temperature of the black hole .    for a generic schwarzschild - de sitter ( sds ) solution , eq .
( [ fe ] ) reads = - .
[ fe1 ] the temperature @xmath174 is related to the horizon radius @xmath161 by t_h=== |-| @xmath161 being a positive solution of the algebraic equation r_h^3-r_h+=0 .
[ ho ] as a result , we finally have = -|r_h-3mg| .
[ fe2 ]    in the pure de sitter case @xmath175 , @xmath176 , and so one has
[ fds ] for sds one needs to consider separately the two admissible cases @xmath166 and @xmath163 . here
we discuss the second one only , that is , the proper sds solution where @xmath177 . in this case , eq .  ( [ ho ] ) has one negative root and two distinct positive roots , @xmath161 , the event horizon and @xmath165 , the cosmological horizon , with @xmath178 .
the positive roots can be written in the form r_c=,12<_c<1 , r_h= , 0<_h<12 , where _ c=12and _ h=12(3-_c ) . in this way
, we get _ c = f_ds , [ fc ] _ h = f_ds , [ fh ] @xmath179 being the expression ( [ fds ] ) .
it has to be noted that , in the latter equations , the factors in front of @xmath179 are always smaller than 1 and thus the energies @xmath180 and @xmath181 are always smaller than @xmath179 , independently of the value of the mass .
moreover , the factor @xmath182 is a monotone decreasing function of @xmath183 , which is equal to 1 for @xmath184 and to 0 for @xmath166 , while the factor @xmath185 is exactly equal to 0 for @xmath186 and quite small otherwise .
its maximum value , a little bit more than @xmath187 , is reached for @xmath188 .
this means that the minimum value for the free energy @xmath181 of a sds black hole is obtained for @xmath189 .
now , we shall study some models and explicitly compute the corresponding free energy .
first of all , as a trivial example , we consider the @xmath190 model described by @xmath191 .
this has a stable de sitter solution , with @xmath192 .
then , we immediately have _ ds=- .
the corresponding sds cosmological and black - hole free energies are plotted in fig .
[ p : lacdm ] , as a function of the mass @xmath148 ( in energy units @xmath193 ) .
and @xmath181 sds free energies for the @xmath0cdm model.,scaledwidth=40.0% ]    as a second example we consider the viable model in ( [ 1.5 ] ) , with @xmath63 , @xmath194 , @xmath195 , @xmath30 , as discussed in sect .
ii . with this choice of parameters , eq .  ( [ b ] ) has three real solutions for @xmath196 ( minkowski ) , @xmath197 and @xmath198 .
the last one ( @xmath134 ) is stable , while the other two are unstable . for the free energies of the two de sitter solutions we get ^(1)_ds=-23 -,^(2)_ds=-23 - .
[ minimods ] we see that ^(2)_ds < f^(1)_ds,~3.636 .
independently of the parameter @xmath22 and so the stable de sitter solution with curvature @xmath134 is favorite , from the energetic point of view , with respect to the one with smaller curvature @xmath133 .
the sds free energies are plotted in fig .
[ p : fcfh ] , as functions of the mass @xmath148 ( in energy units @xmath199 , for @xmath200 .
however , it has to be stressed that , for the sds corresponding to the de sitter solution with curvature @xmath134 , the mass needs to be restricted to @xmath201 and thus , for a fixed mass @xmath148 , we can compare the corresponding free energies in this range only .
for both the cosmological and black - hole horizons there are critical values of the black hole mass @xmath148 , say @xmath202 , for which @xmath203 .
then , we can distinguish three different regions : + ( a ) @xmath204 : in this first case @xmath205 ; + ( b ) @xmath206 : in this second case @xmath207 , but @xmath208 ; + ( c ) finally , for @xmath209 , we have @xmath210 .    as we have seen from general considerations , the free energies corresponding to the cosmological horizons are monotonous functions of the mass , while the ones corresponding to the black - hole horizons are convex functions , which reach the minimum at @xmath211 . for our two cases , @xmath133 and @xmath134
, we get m_1~~ , ^(1)_h , min~- , m_2~~ , ^(2)_h , min~- .
[ minimum ] then we see that , for the model in ( [ 1.5 ] ) we are dealing with , the configuration with the minimum free energy corresponds to sds with mass @xmath212 and curvature @xmath213 .
the associated free energy is given in ( [ minimum ] ) .
we also observe that this is smaller than the free energy of the pure de sitter configuration as it is clear from ( [ minimods ] )    thus , it is demonstrated here that multiple de sitter solutions can appear also under the form of sds solutions .
in other words , the number of multiply solutions becomes significantly bigger . with the appearance of both ds and sds universe solutions
one can suggest various scenarios for the universe evolution .
for instance , we can conjecture that a ( pre-)inflationary universe is described by some sds spacetime . as
time proceeds , this universe decays and enters into the well known radiation / matter dominance phase . in its further evolution , the universe transits to the ds ( or almost ds ) era , by stability and least - energy principle considerations .
the future universe may again appear as an sds spacetime .
): cosmological free energies on the left and black - hole free energies on the right.,title="fig:",scaledwidth=40.0% ] ): cosmological free energies on the left and black - hole free energies on the right.,title="fig:",scaledwidth=40.0% ]
.    in summary
, we have investigated in this paper several viable models of modified gravity which satisfy both the constraints of local as well as cosmological tests . by means of a numerical study
, it is demonstrated that some versions of highly non - linear models exhibit multiply de sitter universe solutions which often appear in pairs , being one of them stable and the other unstable .
the numerical evolution of the effective equation of state parameter is presented too . as a result
, these models can be considered as natural candidates for the unification of early - time inflation with late - time acceleration through ds critical points .
moreover , based on the de sitter solutions , multiply sds solutions can also be constructed .
further , we have investigated the thermodynamic properties of these sds universes : their corresponding entropies and free energies have been calculated and compared .
sds universe might also appear at the ( pre-)inflationary stage .
owing to the highly non - linear structure of the theories under discussion , the ds universes had to be constructed with numerical tools mainly .
moreover , in order to simplify the problem , at this first stage of the investigation we did not to consider matter contributions .
it is clear that , in the next step , we must necessarily include matter and reconsider the problem in its presence .
this has the potential to lead to a sufficiently realistic quantitative description of the universe expansion history , in which modified gravity would be responsible for both acceleration stages : the inflation epoch and the dark energy one .
this quantitative analysis will be presented elsewhere .
* acknowledgements . *
this paper is an outcome of the collaboration program infn ( italy ) and dgicyt ( spain ) .
it has been also supported in part by mec ( spain ) , projects fis2006 - 02842 and pie2007 - 50i023 , by agaur ( generalitat de catalunya ) , contract 2005sgr-00790 and grant 2008be1 - 00180 , and by rfbr , grant 06 - 01 - 00609 ( russia ) .
the work of p.t . was partially supported by rfbr , grant 08 - 02 - 00923 , and with the scientific school grant 4899.2008.2 of the russian ministry of science and technology .
pt thanks alexey toporensky for some useful discussions .
the research of ee is partly based on work done while on leave at the department of physics and astronomy , dartmouth college , 6127 wilder laboratory , hanover , nh 03755 , usa .
s. nojiri and s.d .
odintsov , arxiv:0807.0685[hep - th ] ; t. sotiriou and v. faraoni , arxiv:0805.1726[gr - qc ] ; s. capozziello and m. francaviglia , gen.relat.grav . *
40 * , 357 ( 2008 ) ; f. lobo , arxiv:0807.1640[gr - qc ] .
s. capozziello , int.j.mod.phys .
d * 11 * , 483 ( 2002 ) ; m.  ruggiero and l.  iorio , arxiv : gr - qc/0607093 ; a.  cruz - dombriz and a.  dobado , arxiv : gr - qc/0607118 ; n.  poplawski , arxiv : gr - qc/0610133 ; m.  c.  b.  abdalla et al , class .
quant .  grav .
* 22 * , l35 ( 2005 ) [ arxiv : hep - th/0409177 ] ; a.  brookfield , c.  van  de  bruck and l.  hall , arxiv : hep - th/0608015 ; h. oyaizu , m. lima and w. hu , arxiv:0807.2462[astro - ph ] ; b.  li , k.  chan and m.  chu , arxiv : astro - ph/0610794 ; x.  jin , d.  liu and x.  li , arxiv : astro - ph/0610854 ; t.  sotiriou and s.  liberati , arxiv : gr - qc/0604006 ; t.  sotiriou , arxiv:0805.1160[gr - qc ] ; i.  navarro and k.  van  acoleyen , arxiv : gr - qc/0611127 ; a.  bustelo and d.  barraco , arxiv : gr - qc/0611149 ; f.  briscese , e.  elizalde , s.  nojiri and s.  d.  odintsov , phys
.  lett .
b * 646 * , 105 ( 2007 ) [ arxiv : hep - th/0612220 ] ; b.  li and j.  barrow , arxiv : gr - qc/0701111 ; v.  faraoni , arxiv : gr - qc/0703044 ; arxiv:0706.1223[gr - qc ] ; arxiv:0805.1249[gr - qc ] ; arxiv:0810.2602[gr - qc ] ; o.  bertolami , c.  boehmer , t.  harko and f.  lobo , arxiv:0704.1733[gr - qc ] ; o.  bertolami and j.  paramos , arxiv:0805.1241[gr - qc ] ; s.  carloni , a.  troisi and p.  dunsby , arxiv:0706.0452[gr - qc ] ; arxiv:0707.0106[gr - qc ] ; arxiv:0809.3673[astro - ph ] ; j.  santos , j.  alcaniz , f.  carvalho et al , arxiv:0708.0411[astro - ph ] ; arxiv:0807.2443[astro - ph ] ; h.  mohseni  sadjadi , arxiv:0709.2435[gr - qc ] ; n.  lanahan - tremblay and v.  faraoni , arxiv:0709.4414[gr - qc ] ; n.  goheer , j.  leach and p.  dunsby , arxiv:0710.0814[gr - qc ] ; arxiv:0809.5247[gr - qc ] ; o.  lecian and g.  montani , arxiv:0807.4428[gr - qc ] ; t. harko , arxiv:0810.0742[gr - qc ] ; b. vakili , arxiv:0809.4591[gr - qc ] ; s. srivastava , arxiv:0809.1950[gr - qc ] ; d. saez , arxiv:0809.1311[hep - th ] ; k. bamba et al , arxiv:0803.3384[hep - th ] , phys.rev . d*77 * , 123532 ( 2008 ) ; arxiv:0801.0954[astro - ph ] , jcap * 0804*,024 ( 2008 ) ; arxiv:0807.2575[hep - th ] .
p.  mahato , annales de la fondation louis de broglie 32 ( 2007 ) 297 - 310;arxiv:0710.3821 [ astro - ph ] .
s.  nojiri and s.  odintsov , gen .
grav .  * 36 * , 1765 ( 2004 ) [ arxiv : hep - th/0308176 ] ; phys .
b * 576 * , 5 ( 2003 ) [ arxiv : hep - th/0307071 ] ; phys .
* b599 * , 137 ( 2004 ) [ arxiv : astro - ph/0403622 ] ; arxiv:0801.4843[astro - ph ] ; g.  cognola , e.  elizalde , s.  nojiri , s.  d.  odintsov and s.  zerbini , phys .  rev .  * d73 * , 084007 ( 2006 ) , [ arxiv : hep - th/0601008 ] ; s.  capozziello et al , arxiv : astro - ph/0501426 ; arxiv : hep - th/0512118 ; g.  allemandi , a.  borowiec , m.  francaviglia and s.  d.  odintsov , arxiv : gr - qc/0504057 ; g.  allemandi , m.  francaviglia , m.  ruggiero and a.  tartaglia , arxiv : gr - qc/0506123 ; t.  multamaki and i.  vilja , arxiv : astro - ph/0612775 ; arxiv:0709.3422[astro - ph ] ; j.  a.  r.  cembranos , phys .  rev .
d * 73 * , 064029 ( 2006 ) [ arxiv : gr - qc/0507039 ] ; t.  koivisto and h.  kurki - suonio , arxiv : astro - ph/0509422 ; t.  clifton and j.  barrow , arxiv : gr - qc/0509059 ; i.  brevik , arxiv : gr - qc/0601100 ; t.  faulkner , m.  tegmark , e.  bunn and y.  mao , arxiv : astro - ph/0612569 ; s.  capozziello et al , arxiv:0804.2163[gr - qc ] ; arxiv:0712.2980[gr - qc ] ; arxiv:0810.2549[gr - qc ] ; s.  nojiri , s.  d.  odintsov and p.  tretyakov , arxiv:0704.2520[hep - th ] ; m.  movahed , s.  baghram and s.  rahvar , arxiv:0705.0889[astro - ph ] ; c.  boehmer , t.  harko and f.  lobo , arxiv:0710.0966[gr - qc ] ; j.  miritzis , arxiv:0708.1396[gr - qc ] ; d.  vollick , arxiv:0710.1859[gr - qc ] ; j.  evans , l.  hall and p.  caillol , arxiv:0711.3695[astro - ph ] ; t. clifton , arxiv:0807.4682[gr - qc ] ; h. stefancic , arxiv:0807.4428[gr - qc ] ; a. dev et al , arxiv:0807.3445[hep - th ] ; t. kobayashi and k. maeda , arxiv:0807.2503[astro - ph ] ;    s.  nojiri , s.  d.  odintsov , phys .
d * 74 * , ( 2006 ) 086005 [ arxiv : hep - th/0608008 ] ; arxiv : hep - th/0611071 ; s.  capozziello , s.  nojiri , s.  d.  odintsov and a.  troisi , phys .
b * 639 * , 135 ( 2006 ) [ arxiv : astro - ph/0604431 ] .
s.  tsujikawa , k.  uddin and r.  tavakol , arxiv:0712.0082 ; s.  fay , s.  nesseris and l.  perivolaropoulos , arxiv : gr - qc/0703006 .
g.  cognola , m.  gastaldi and s.  zerbini , arxiv : gr - qc/0701138 ; j.l .
cortes and j. indurain , arxiv:0805.3481[astro - ph ] ; s. tsujikawa and t. tatekawa , arxiv:0804.4343[astro - ph ] ; n. deruelle , m. sasaki and y. sendouda , phys.rev .
d * 77*,124024 ( 2008 ) ; k. kainulainen and d. sunhede , arxiv:0803.0867[gr - qc ] ; g. cognola and s. zerbini , j.phys .
a * 39 * , 6245 ( 2006 ) .
s.  a.  appleby and r.  a.  battye , arxiv:0705.3199[astro - ph ] ; l.  pogosian and a.  silvestri , arxiv:0709.0296[astro - ph ] ; s.  tsujikawa , arxiv:0709.1391[astro - ph ] ; s.  capozziello and s.  tsujikawa , arxiv:0712.2268[gr - qc ] .
i.brevik , s. nojiri , s.d .
odintsov and l. vanzo , hep - th 0401073 , phys.rev .
d * 70 * , 043520 ( 2004 ) ; g. cognola , s. nojiri , s.d .
odintsov and s. zerbini , hep - th/0501096 , jcap * 0502 * , 010 ( 2005 ) ; f.  briscese and e.  elizalde , phys .
d * 77 * , 044009 ( 2008 ) [ arxiv:0708.0432 [ hep - th ] ] .


physics of graphene perception begins when one compares landau levels in two dimensional schrdinger and dirac theories .
such spectacular phenomenon as unconventional quantum hall effect @xcite is caused by the anomaly of the lowest landau level ( lll ) @xcite which for dirac fermions in graphene is field independent and can accommodate only half the usual number of the states from the conduction band and takes the other half from the valence band .
the easiest way to accomplish this peculiar feature of the lll is to solve a pair of the dirac equations that describe excitations near two inequivalent @xmath0 points of graphene s brillouin zone .
normally this is done in a constant homogeneous magnetic field , although this property of the lll for dirac fermions is topologically protected for inhomogeneous field configurations and in the presence of ripples @xcite .
the simplest inhomogeneous field configuration which contains nontrivial aharonov - bohm physics can be created by an infinitesimally thin solenoid . in practice
such magnetic field configuration may be obtained when a type - ii superconductor is placed on top of graphene or semicondunducting heterojunction hosting a 2d electron gas ( 2deg ) with quadratic dispersion .
while graphene devices still have to be fabricated , devices like this with a superconducting film grown on top of a semiconducting heterojunction ( such as gaas / algaas ) hosting a 2deg have in fact been fabricated twenty years ago @xcite and theoretically well studied ( see e.g. refs .  ) .
theoretically a problem of the dirac fermions in the field of aharonov - bohm flux was encountered in the context of cosmic strings by gerbert and jackiw @xcite .
while for the solutions of the dirac equation with nonzero angular momentum the square integrability requirement specifies which of the two independent solutions should be taken , they noticed this is not the case for the zero angular momentum .
for the zero momentum there is an ambiguity as both solutions are square integrable , but divergent as @xmath1 at the origin , where @xmath2 is the space coordinate .
the ambiguity of the solution selection is caused by the singular nature of the infinitesimally thin solenoid vector potential at the origin .
this problem has initiated a vast theoretical literature which addresses interesting aspects related to the rigorous treatment of the solutions of the dirac equation with the aharonov - bohm potential ( see e.g. refs .   for a review ) . in the condensed matter context the dirac fermions in the field of solenoid emerged during the study of the dirac - bogolyubov - de gennes quasiparticles in the vortex state of @xmath3-wave superconductors
@xcite ( see also @xcite for a review ) . due to
the divergence of the zero mode solutions theory predicts a formation of nonzero local density of states ( ldos ) near the vortex center .
however , this theoretical prediction based on the dirac nature of quasiparticles in @xmath3-wave superconductors does not agree with the results of scanning tunneling spectroscopy ( sts ) measurements @xcite in high - temperature superconductors .
finally we mention a related problem of the description of topological defects in graphene based on the dirac equation with a pointlike pseudomagnetic vortex which has also been studied intensively , see e.g. refs .  .
the purpose of the present paper is to study the electronic excitations in graphene in the field of the aharonov - bohm flux and compare them with the corresponding results for 2deg with a quadratic dispersion .
we rely on the existing studies of the dirac fermions in the aharonov - bohm potential , but focus on the specific features of graphene such as the presence of two inequivalent @xmath0 points which implies that one should consider the solutions for both inequivalent irreducible representations of the dirac @xmath4 matrices . also to avoid unnecessary formal complications we consider the physical regularization of the problem modeling a finite radius flux tube created by the abrikosov vortex .
we utilize the simplest case of magnetic field concentrated in a thin cylindrical shell of small but finite radius @xmath5 when @xmath6 @xcite .
in contrast to high - temperature superconductors dirac description of the quasiparticles in graphene is proven valid under the different conditions . in particular , sts measurements of graphene flakes on graphite
@xcite exhibit the structural and electronic properties expected of pristine graphene such as the development of a single sequence of pronounced landau level peaks corresponding to massless dirac fermions in a homogeneous magnetic field .
we propose to perform sts measurements for graphene penetrated by vortices from a type - ii superconductor , because the dirac theory predicts rather peculiar behavior of ldos not expected for the 2deg with a quadratic dispersion of carriers .
the paper is organized as follows . in sec .
[ sec : model ] we introduce the model hamiltonians and discuss the regularization of the aharonov - bohm potential used in this work .
[ sec : schrodinger ] is devoted to the nonrelativistic case , and the relativistic case is discussed in detail in sec .
[ sec : dirac ] .
the structure of both sections is the same : we consider the solution of the corresponding schrdinger or dirac equation , construct the green s function ( gf ) with coinciding arguments , obtain the dos and study the behavior of the ldos . in sec .
[ sec : concl ] our final results are summarized .
we consider nonrelativistic and relativistic hamiltonians .
the 2d nonrelativistic ( schrdinger ) hamiltonian has the standard form @xmath7 where @xmath8 , @xmath9 with the vector potential @xmath10 , planck s constant @xmath11 and the velocity of light @xmath12 describes a spinless particle with a mass @xmath13 and charge @xmath14 .
the dirac quasiparticle in graphene is described by the hamiltonian @xmath15 where the matrices @xmath16 and @xmath17 are defined in terms of the pauli matrices as @xmath18 here @xmath19 labels two unitary inequivalent representations of @xmath4 gamma matrices , so that one considers a pair of dirac equations corresponding to two inequivalent @xmath20 points of graphene s brillouin zone . in eq .
( [ hamilton - rel ] ) @xmath21 is the fermi velocity and @xmath22 is the dirac mass ( gap ) , which is introduced in the hamiltonian for generality .
note that we consider the simplest case when the gap has the same sign for @xmath19 [ see ref .   for a discussion of more general cases ] . while tight binding calculations show that the quasiparticle excitations in graphene have a linear dispersion at low energies @xcite and are described by the massless dirac equation with @xmath23 @xcite , recent sts measurements revealed a mass gap near the dirac point in a single layer graphene sample suspended above a graphite substrate @xcite .
since this gap and its origin are intensively studied both theoretically and experimentally in the last few years , here we consider a generic case with a finite value of @xmath22 .
the vector potential of a vortex at the origin directed in the @xmath24 direction is @xmath25 where @xmath26 is the flux of the vortex expressed via magnetic flux quantum of the electron @xmath27 with @xmath28 , . ] where the value @xmath29 corresponds to the flux created by the abrikosov vortex .
the magnetic field is then @xmath30    the essential difference between the schrdinger @xmath31 and dirac @xmath32 equations in this case can be seen if one squares the latter : @xmath33   \right ) \psi = ( e^2 - \delta^2 ) \psi,\ ] ] where the commutator @xmath34 = -\frac{e}{\hbar c } b_z(\mathbf{r})\ ] ] which introduces pseudo - zeeman term which is related to the sublattice rather than the spin degree of freedom .
it should be mentioned that in the case of graphene the components of the spinor @xmath35 are associated with a sublattice rather than a spin degree of freedom .
since the hamiltonian ( [ hamilton - rel ] ) originates from a nonrelativistic many - body theory , the zeeman interaction term has to be explicitly added to this hamiltonian .
this resembles the case of the nonrelativistic hamiltonian ( [ hamilton - nonrel ] ) which becomes pauli one when the interaction between the magnetic moment of the spin and an external magnetic field is added . in the present paper
we do not include the spin degree of freedom neither in ( [ hamilton - nonrel ] ) nor in ( [ hamilton - rel ] ) .
( [ dirac - squared ] ) and ( [ commutator ] ) identify the origin of complications @xcite in the problem with a singular vortex ( [ a - singular ] ) when a singularity in @xmath36 occurs at a singular point of the differential equation ( [ dirac - squared ] ) . to avoid these complications one can consider a vortex with a finite radius flux tube @xcite , i.e. with the magnetic field and vector potential written in cylindric coordinates @xmath37 : @xmath38 where @xmath39 is a profile function with a compact support satisfying the normalization @xmath40 and connected to the profile function @xmath41 by the relation @xmath42 .    the simplest choice of the field distribution @xmath39 which regularizes the problem with the solutions solely expressed in terms of bessel functions is a magnetic field concentrated on the surface of the cylinder of radius @xmath5 , @xmath43 .
then , the corresponding profile function @xmath44 in the limit @xmath45 we recover the aharonov - bohm potential ( [ a - singular ] ) but avoiding formal complications . as shown in ref .
@xcite , there is no dependence on the detailed form of @xmath39 in the limit @xmath6 provided that @xmath46 . in the present paper
we restrict ourselves by considering the profile function ( [ step - profile ] ) .
in this section we recapitulate the results of refs .   for nonrelativistic case .
they are important not only for comparison with relativistic case , but also because the relativistic result is constructed using the nonrelativistic one .      in the limit @xmath47
the admissible solution of the schrdinger equation is always a regular solution which in polar coordinates @xmath48 takes the form @xmath49 where @xmath50 is the bessel function with the wave vector @xmath51 which is related to the quasiparticle energy @xmath52 via @xmath53 .    the eigenfunction expansion for the retarded schrdinger gf reads @xmath54 or after substituting the wave function ( [ solution - nonrel ] ) it becomes @xmath55 where @xmath56 . since an analytic continuation of the gf ( [ nonrel - gf ] ) on the imaginary axis in the complex momentum plane , @xmath57 is free of singularities , it is convenient to work with the corresponding gf @xmath58 where @xmath59 in eq .
( [ nonrel - gf - continued ] ) we already set two arguments coinciding , @xmath60 , because in the present work we consider the dos only . as we will see below , the function @xmath61 is also used in the representation of the dos for the dirac fermions . after the calculation of the gf ( [ g - continued ] ) is done ,
the ldos per spin projection can be found by returning back to the real momentum axis @xmath62 the gf @xmath63 was calculated in ref .   using the contour integration technique . a weak point of this calculation
was discussed in @xcite , where the same method was applied to obtain the @xmath64-dependent contribution to the gf , @xmath65 with coinciding arguments when the approach of @xcite is valid . referring to the derivation of @xcite , here we simply start from the corresponding expression for @xmath66 obtained in ref .
@xmath67 notice that eq .
( [ delta - g - marino ] ) coincides with the corresponding expression from @xcite up to a coefficient . changing the variables to @xmath68
, @xmath69 one can obtain from ( [ delta - g - marino ] ) the final expression @xmath70}{\cosh x } e^{-2 \mathcal{q } r \cosh x\cosh y}\ ] ] which we will use in what follows .
it turns out that the investigation of the full dos is simpler than the analysis of the ldos .
thus in the next sec .
[ sec : dos - nonrel ] we firstly consider the dos and return to the ldos ( [ ldos - def ] ) below in sec .
[ sec : ldos - nonrel ] .
the full dos per spin projection is obtained from the ldos ( [ ldos - def ] ) by integrating over the space coordinates @xmath71 since we have the integral representation ( [ delta - g ] ) for @xmath72 , it is straightforward to calculate directly the perturbation of dos , @xmath73 induced by the aharonov - bohm potential . here @xmath74 is a free dos of 2deg per unit area and @xmath75 is the 2d volume ( area ) of the system .
firstly integrating @xmath72 over space , one obtains @xcite @xmath76 then integrating over @xmath77 and @xmath78 we obtain that @xmath79 returning to the real @xmath80 axis we reproduce the usual aharonov - bohm depletion of the dos with respect to the free dos @xmath81 @xcite : @xmath82 writing eq .
( [ norel - dos - depletion ] ) we have also included a case of the opposite field direction .
now we come back to the ldos ( [ ldos - def ] ) . as in the case of the full dos , it is convenient to consider the excess ldos , @xmath83
. then the value @xmath84 can be obtained by calculating the function @xmath85 given by eq .
( [ delta - g ] ) and substituting the result to eq .
( [ nonrel - gf - continued ] ) .
the analytic continuation @xmath86 described by eq .
( [ ldos - def ] ) has to be done at the very last step of the calculation .
thus our purpose is to derive a simple representation for the function @xmath87 .
first , one can rewrite it in the form @xmath88 where we used that @xmath89 . differentiating ( [ delta - g ] )
we get the integrand of the last expression @xmath90\int_0^\infty dy   \cosh y e^{-2\mathcal{q }
r\cosh x\cosh y}.\ ] ] using the integral representation for the modified bessel function @xmath91 @xcite @xmath92 and the formula ( 2.16.13.2 ) from @xcite @xmath93 we come to the equation @xmath94 integrating the last expression with mathematica , we get @xmath95 - \ln \frac{\mathcal{q } r}{2 } \\ &   - \frac{\mathcal{q}^2 r^2}{4 \eta(1-\eta)}\,_3f_4(1,1,3/2;2,2,2-\eta , 1+\eta;\mathcal{q}^2 r^2 )   + \frac{1}{2}\psi(1-\eta)+ \frac{1}{2}\psi(\eta ) , \end{split}\ ] ] where @xmath96 is the generalized hypergeometric function and @xmath97 is the logarithmic derivative of the gamma function @xmath98 .
after the analytic continuation @xmath99 is made only two terms in the square brackets ( with the hypergeometric function itself remaining real ) and logarithmic term contribute in @xmath100 , so that @xmath101 -1\}\ ] ] with @xmath102 note that the last term @xmath103 in eq
.  ( [ ldos - function ] ) arises from the logarithmic term of ( [ hypergeometric - integral ] ) .
in the important limiting cases the last expression is greatly simplified .
for example , in the limit @xmath104 one obtains @xmath105 also in the physically important case @xmath29 the ldos is expressed in terms of the sine integral , @xmath106 as follows @xmath107.\ ] ] using the asymptotic of the integral sine @xcite @xmath108 for @xmath109 from ( [ ldos - si ] ) we recover the previous expression ( [ ldos - large ] ) valid for @xmath29 and @xmath104 . in the opposite limit @xmath110 for @xmath111 , we see that @xmath112 $ ] .    in fig .
[ fig:1 ] we show the dependence ( [ ldos - si ] ) of the induced ldos @xmath113 on the distance from the center of the vortex @xmath2 .    we observe that in the case of nonrelativistic 2deg the presence of the vortex induced the depletion of the ldos for small @xmath114 .
the function @xmath115 crosses zero near @xmath116 and for @xmath117 the function it reaches the maximal value @xmath118 . in fig .
we model a situation when the stm tip is positioned at some distance from the center of vortex and a bias voltage is swept to explore the energy dependence of the ldos . to take into account the presence of a finite carrier density in 2deg
, we introduce a finite fermi energy @xmath119 , so that the ldos at zero energy , @xmath120 , corresponds its value at the fermi level , i.e. @xmath121 . to choose the appropriate units we set the distance scale @xmath122 to be the order of the lattice constant .
then the energy scale , @xmath123 is the order of the bandwidth .
the dimensionless variable @xmath124 can now be rewritten as follows @xmath125 .
the dependence of @xmath126 is shown in fig .
[ fig:2 ] for three values of @xmath127 : for @xmath128  solid ( blue ) curve , for @xmath129  long - dashed ( red ) curve and for @xmath130  short - dashed ( black ) curve .
the chemical potential is taken @xmath131 .
we observe that only for the smallest ratio @xmath132 the values of the ldos are significantly depleted below the free ldos @xmath133 . as we saw
, the depletion of the ldos occurs for @xmath134 . since the presence of the fermi surface makes the value of @xmath80 large , so that the region of small @xmath135 is accessible only for @xmath136 .
indeed , we observe that only for the smallest ratio @xmath132 the values of the ldos are depleted to a half of the value of free ldos @xmath133 .
since the realistic values of the vortex core size are at least of the order of magnitude larger than @xmath122 , this implies that the region of a significant depletion of the ldos is not accessible experimentally .
still , due to the slow decay of @xmath137 even for @xmath138 the amplitude of @xmath139 oscillations is of order of @xmath140 , so that this behavior can be probably observed experimentally .
in this section we examine the density of states for dirac particles in the potential of the infinitesimally thin solenoid . to avoid formal complications we consider the physical regularization of the problem with the magnetic field concentrated in a thin cylindrical shell of small but finite radius @xmath5 and
take the limit @xmath45 at the end of the calculation .
all answers are presented in the form convenient for comparison with the case of the schrdinger equation .
the dirac equation ( [ dirac - eq ] ) in the field of the regularized vortex ( [ profile ] ) , ( [ step - profile ] ) ( see also ( [ a - regular ] ) ) is solved in appendix  [ sec : a ] . in our consideration
we follow refs .  .
the profile ( [ a - regular ] ) implies that for @xmath141 the particle obeys the dirac equation ( [ system : r < r ] ) for a free particle , while for @xmath142 the particle moves in the field of aharonov - bohm vortex ( [ system : r > r ] ) .
accordingly , for @xmath141 the squared dirac equation ( [ dirac - squared ] ) is equivalent to free schrdinger equations for the components of the spinor @xmath143 notice that in the appendix  [ sec : a ] the definition ( [ psi - def - appendix ] ) for @xmath144 explicitly includes the factor @xmath145 .
for @xmath142 the components @xmath35 satisfy the schrdinger equations with aharonov - bohm potential and the commutator ( [ commutator ] ) is singular at @xmath146 : @xmath147=-\frac{\eta}{r}\delta(r - r).\ ] ] the solution of the problem can found by matching the solutions obtained in the domains @xmath148 and @xmath142 [ see e.g. ( [ solution - psi1 ] ) ] .
the radial components of the spinor @xmath149 have to be continuous : @xmath150 and the singularity of the commutator ( [ commutator - cylinder ] ) is taken into account by a condition on the derivatives @xmath151 we stress that in contrast to the dirac equation case , for the nonrelativistic case when the solution ( [ solution - nonrel ] ) is obtained using the same regularization procedure , both the wave function and its derivative should be continuous .
this as we saw from eqs .
( [ dirac - squared ] ) , ( [ commutator ] ) and ( [ components - cartesian ] ) is related to the pseudo - zeeman term .
after the limit @xmath45 is taken we obtain for the case @xmath152 the following solutions : @xmath153 for positive value of the energy @xmath154 and @xmath155 ; @xmath156 for negative value of energy @xmath157 and @xmath155 . here the upper and lower sign of the second spinor s component corresponds to @xmath158 and @xmath159 solutions , respectively .
the @xmath160 solution turns out to be special , while the upper component is regular at @xmath161 , the lower component diverges as @xmath162 : @xmath163 for @xmath164 and @xmath165 for @xmath166 .
one can check this property by using the matching conditions for @xmath167 component and then finding @xmath168 from eq .
( [ system : r > r - for2 ] ) . on the other hand
, one can discover this singularity in the limit @xmath6 by directly analyzing the matching conditions for @xmath168 and then finding nonsingular @xmath167 from eq .
( [ system : r > r - for1 ] ) . comparing eqs .
( [ m=0.e>0 ] ) and ( [ m=0.e<0 ] ) we observe that a singular at @xmath161 zero - mode solution @xmath169 is hole - like , because a singular electron - like solution vanishes due to the @xmath170 factor .
concluding the discussion of the solutions for the case @xmath152 we stress that for the opposite field direction it is the upper component of the spinor which is singular in the solution equivalent to eqs .
( [ m=0.e>0 ] ) and ( [ m=0.e<0 ] ) .
moreover , the nonvanishing zero mode is now electron - like . as it was firstly noticed in ref .
, this behavior under the change in the field direction breaks the symmetry under @xmath171 ( see also ref .   for a discussion ) .
the set of solutions for the case @xmath172 is the following : @xmath173 for @xmath164 and @xmath174 for @xmath166 .
the prescription for the upper and lower sign of the second spinor s component is the same as for @xmath152 .
the @xmath160 solutions in the case @xmath172 are the following : @xmath175 for @xmath164 and @xmath176 for @xmath166 .
we observe that in this case it is the upper component of the spinor which is singular at @xmath161 and nonvanishing is the electron - like zero mode .
now we are at the position to construct the gf for the dirac fermions using the presented above solutions . however , before going to this , we should stress that as shown in@xcite the field configuration with @xmath177 confined to the surface of a cylinder of radius @xmath5 is not essential for the main result . in simple words , an
more general and not singular at the origin potential can be considered as a set of concentric shells . or putting more formally
, a profile function @xmath39 should satisfy the condition given below eq .
( [ step - profile ] ) which excludes delta - funtion at the origin .    the eigenfunction expansion for the retarded dirac gf s now includes both positive and negative energy solutions @xmath178 where depending on the sign of @xmath179 the sum is taken over either @xmath152 or @xmath172 solutions .
accordingly , for the diagonal matric elements of @xmath180 , for instance , in the @xmath152 case we obtain @xmath181 and @xmath182 , \end{split}\ ] ] where @xmath183 .
then the gf with coinciding arguments , @xmath60 acquires the following form @xmath184 \end{split}\ ] ] for @xmath152 , and @xmath185 , \\
g^{\mathrm{d}}_{\eta22}(\mathbf{r},\mathbf{r},e+i0;\zeta=-1 ) & = \frac{e-\delta}{2\pi(\hbar v_f)^2}g_\eta(\mathbf{r},q ) , \end{split}\ ] ] for @xmath172 . here
the function @xmath186 is related to the function @xmath61 defined in eq .
( [ g - continued ] ) by analytic continuation @xmath187 . note that while the function @xmath61 is identical for both nonrelativistic and relativistic cases ,
the real momentum @xmath80 function @xmath186 has a different analytical properties in these cases reflecting the fact that in contrast to the nonrelativistic case the relativistic spectrum contains positive and negative energy branches . similarly , the function @xmath188.\ ] ] is the analytic continuation to the imaginary axis of the function @xmath189 .
the @xmath190 contribution to @xmath191 originates from the zero mode solutions of the dirac equation . using the integral [
( 2.12.32.12 ) from@xcite ] @xmath192 where @xmath193 is the modified bessel function and @xmath194 is the macdonald function , we obtain the following simple result for @xmath190 @xmath195      due to the matrix structure of the gf and the presence of the valley degree of freedom , @xmath196 , in contrast to the nonrelativistic expression ( [ dos - nonrel ] ) , the full dos in the relativistic case involves not only the integration over area , but also summing over the diagonal components of the gf and @xmath20 valleys . for a better understanding of the final result for the full dos ,
it is instructive to consider separate expressions for the dos corresponding separate @xmath20 points @xmath197.\ ] ] recall that for @xmath198 the free dos of the dirac quasiparticles per spin and one valley is equal to @xmath199 .    using eqs .
( [ dirac - g - zeta=1 ] ) and ( [ dirac - g - zeta=-1 ] ) one can express the dos ( [ dirac - dos - zeta ] ) via the functions @xmath200 and @xmath201 defined by eqs .
( [ g - continued ] ) and ( [ f - continued ] ) , respectively , as follows @xmath202 . \end{split}\ ] ] the integral of used in the nonrelativistic case function @xmath200 , or to be precise of the function @xmath87 , over the space coordinates is calculated in eq .
( [ g - integral ] ) .
the corresponding space integration of the function @xmath201 given by eq .
( [ f - macdonald ] ) produces the result @xmath203 notice that since @xmath204 for @xmath198 , there is no need to introduce a function @xmath205 . having these results for the spatial integration we can calculate the imaginary part . since both eqs .
( [ g - integral ] ) and ( [ f - integral ] ) depend on @xmath206 , the imaginary part is obtained by using a simple prescription @xmath207/(\hbar v_f)^2 $ ] .
this gives the final expression for the perturbed dos , @xmath208 by the aharonov - bohm vortex @xmath209 +    ( [ dirac - dos - final ] ) which is @xmath210 originates from @xmath211 part of eq .
( [ dirac - dos - g&f ] ) and resembles the nonrelativistic result ( [ norel - dos - depletion ] ) . for @xmath23 it turns out to be twice larger than ( [ norel - dos - depletion ] ) simply because in eq .
( [ dirac - dos - final ] ) we summed over the diagonal components of the gf @xmath180 which are related to the two sublattices of graphene .
the last term of ( [ dirac - dos - final ] ) originates from @xmath190 part of eq .
( [ dirac - dos - g&f ] ) and thus is related to the zero mode solution of the dirac equation .
the mentioned above fact that the singular component of the zero mode solution is hole - like for @xmath152 [ see eq .
( [ m=0.e<0 ] ) ] and it is electron - like for @xmath172 [ see eq .
( [ m=0.e>0-negative - zeta ] ) ] finds its reflection in the asymmetric form of the last term of eq .
( [ dirac - dos - final ] ) which also corresponds to the holes ( electrons ) for @xmath212 ( @xmath172 ) . in accordance with the behavior of the zero mode solutions described in sec .
[ sec : dirac - solutions ] , when the direction of the field is reversed the expressions @xmath213 and @xmath214 are interchanged .
the full excess dos @xmath215\ ] ] is obviously symmetric in energy .
in contrast to the nonrelativistic case , in the dirac case the aharonov - bohm vortex induces the excess of the states which is related to the presence of the last term of eq .
( [ dirac - dos - final ] ) and caused by the zero modes .
we note that in ref .   the corresponding term of eq .
( [ dirac - dos - final ] ) has a wrong sign .
the positiveness of @xmath216 can also be understood by the following simple argument . for @xmath217 the free dos @xmath218 .
therefore , since the dos has to be positive , the value @xmath216 should also be positive .
now we investigate the ldos for the dirac case . while it was useful to consider each valley separately , especially
because field theoretical studies of the problem often involve only one unitary inequivalent representation of @xmath219 gamma matrices [ see e.g.@xcite ] , the ldos measurement picks up both valleys together . on the other hand , ldos distinguishes sublattices .
thus we consider separately the ldos for @xmath220 and @xmath177 sublattices which are defined as follows @xmath221 , \\
n_\eta^{\mathrm{d } ( b)}(\mathbf{r},e)=-\frac{1}{\pi } \mbox{im } \left [ g_{\eta 22}(\mathbf{r},\mathbf{r},e+i0;\zeta=1 ) + g_{\eta22}(\mathbf{r},\mathbf{r},e+i0;\zeta=-1 ) \right ] .
\end{split}\ ] ] again we consider the perturbation of the ldos induced by the aharonov - bohm potential , @xmath222 .
using eqs .
( [ dirac - g - zeta=1 ] ) and ( [ dirac - g - zeta=-1 ] ) we rewrite the ldos s @xmath223 in terms of the functions @xmath224 and @xmath201 as follows @xmath225.\ ] ] here the upper ( lower ) sign corresponds to @xmath220 ( @xmath177 ) sublattice . to obtain the final expression for ldos
we should make the analytic continuation @xmath226 .
for the function @xmath227 given by eq .
( [ hypergeometric - integral ] ) one obtains @xmath228 - 1\right\},\ ] ] where the function @xmath229 is defined in eq .
( [ f - eta ] ) .
one can easily see that because in the nonrelativistic case the analytic continuation to the real momentum @xmath99 is different from the analytic continuation @xmath226 in the relativistic case , the function ( [ g - real - rel ] ) differs from the function ( [ ldos - function ] ) . using the relationships @xcite @xmath230 between the macdonald function of the imaginary argument , the hankel function of the second kind @xmath231 and the bessel function of the first @xmath232 and second @xmath233 kinds , we obtain that @xmath234,\ ] ] where we used the property @xmath235 .
accordingly , the analytic continuation of the function @xmath236 from eq .
( [ f - macdonald ] ) takes the form @xmath237 - 2\cos\pi\eta j_\eta(qr)y_\eta(qr)\}.\ ] ] substituting ( [ g - real - rel ] ) and ( [ f - real - rel ] ) in eq .
( [ ldos - ab - f&g ] ) we arrive at the final main result @xmath238 - 1\\ + & \frac{\sin^2(\pi\eta)}{2 } [ y^2_\eta(qr)- j^2_\eta(qr)]- \frac{\sin ( 2 \pi\eta)}{2 } j_\eta(qr)y_\eta(qr ) \biggr\ } , \qquad q = \frac{\sqrt{e^2 - \delta^2}}{\hbar v_f } , \end{split}\ ] ] where @xmath239 is free dos of the dirac quasiparticles per spin and one sublattice ( or valley ) for @xmath23 .
the first part of eq .
( [ ldos - ab - final ] ) which includes @xmath240 and @xmath241 is identical to the nonrelativistic expression ( [ ldos - function ] ) , while the second part of eq .
( [ ldos - ab - final ] ) with bessel functions originates from the zero mode contribution .    in the limit
@xmath104 the last expression acquires a simple form @xmath242 comparing eqs .
( [ ldos - ab - large - qr ] ) and ( [ ldos - large ] ) we conclude that in the dirac case the impact of the vortex is more localized than in the nonrelativistic case .    as we saw in sec .
[ sec : ldos - nonrel ] in the physically important case @xmath29 the expression ( [ ldos - function ] ) is significantly simplified to the result ( [ ldos - si ] ) .
the same remains true for eq .
( [ ldos - ab - final ] ) , because half - integer bessel functions are expressed in terms of the elementary functions and we obtain that @xmath243.\ ] ] from eq .
( [ ldos - ab - si - cos ] ) we immediately observe the main difference between the relativistic and nonrelativistic cases .
the presence of zero modes causes a positive divergence of the ldos , @xmath244 for @xmath114 near the center of the vortex .
integrally this results in the excess of the states in the full dos ( [ dirac - dos - full ] ) . using the asymptotic expansion of @xmath245 given below eq .
( [ ldos - si ] ) we recover the previous expression ( [ ldos - ab - large - qr ] ) valid for @xmath29 and @xmath104 .    in fig .
[ fig:3 ] we show the dependence ( [ ldos - ab - si - cos ] ) of the induced ldos @xmath246 on the distance from the center of the vortex @xmath2 for @xmath23 .
we observe the features expected from the analytic expressions such as the excess of the ldos for small @xmath114 and faster than in 2deg decay of @xmath247 for @xmath104 .    in fig .
[ fig:4 ] we consider a situation similar to fig .  [ fig:2 ] .
we fix the distance at @xmath248 and plot the energy dependence of the relativistic ldos ( [ ldos - ab - si - cos ] ) .
we consider the most interesting case of undoped graphene with zero carrier density .
in contrast to the 2deg , in graphene is easily tuned to this regime .
again , we introduce the distance scale @xmath122 of the order of the lattice constant .
then for the dirac case the energy scale is @xmath249 , and accordingly the dimensionless variable is @xmath250 .    to make sublattices inequivalent we introduce a finite gap @xmath251 which introduces the asymmetry between the the ldos on @xmath220 and @xmath177 sublattices .
the ldos @xmath252 shown as a solid ( blue ) curve has a sharp peak near @xmath253 and practically no peak at @xmath253 , while the ldos @xmath254 shown as dashed ( red ) curve has a peak near @xmath255 and no peak at @xmath256 .
when @xmath257 is increasing the ldos very quickly reduces to its value for the dirac system in the absence of the vortex .
we note that this behavior of the ldos is seen for the large value of @xmath138 , i.e. far from the center of the vortex , indicates that it should be possible to observe this excess of the ldos in experiments on graphene .    finishing our discussion of the dirac case ,
we mention a recent work @xcite , where the effect of the vacuum polarization in the field of an infinitesimally thin solenoid at the distances much larger than the radius of solenoid was studied . constructing the gf the authors neglected the delta - function ( [ discontinuity ] ) motivated by the fact that they are interested in the regime @xmath258 .
one of their interesting conclusions is that for @xmath29 the induced current is zero .
the main motivation of this work was to address a natural question whether one can distinguish graphene from 2deg measuring the ldos near the abrikosov vortex penetrating them . just by comparing figs .
[ fig:2 ] and [ fig:4 ] we find a positive answer . indeed , by putting the stm tip close to the vortex , for example , at the distance , @xmath259 , where @xmath122 is the lattice constant , we observe that in 2deg the value of the ldos will be close to its background value @xmath260 in the free system , while in undoped graphene one should observe a strong ldos enhancement near the fermi level . or to be more precise , if there is no gap @xmath22 in the quasiparticle spectrum , the peak of the ldos should be observed at the fermi level , or if @xmath261 the peaks should be observed at the energies @xmath262 and the sign of the peak energy should depend on the sublattice .
this peaked behavior of the ldos with @xmath244 for @xmath114 reflects the specific feature of the dirac fermions such as the presence of the divergent as @xmath1 at the origin zero mode solution of the dirac equation .
thus the observation of this feature in sts measurements would contribute to the expanding list of the experimental manifestations of the dirac fermions in graphene .
our second conclusion is that while in the nonrelativistic case the presence of the aharonov - bohm vortex leads to a depletion of the full dos , in the dirac case the full dos is enhanced .
this result can likely be checked by analyzing the sts maps .
we remind at the end that in this paper we did not consider the effect of disorder in the presence of magnetic field which will definitely affect the behavior of the ldos in magnetic field @xcite and should be taken into account in the analysis of the sts data .
we thank e.v .  gorbar and v.p .
gusynin for many stimulating discussions .
s.g.s . thanks c.g .
beneventano and e.m .
santangelo for an interesting discussion and bringing ref .   to his attention .
this work was supported by the scopes grant no .
iz73z0`_`128026 of swiss nsf , by the grant simtech no .
246937 of the european fp7 program and by the program of fundamental research of the physics and astronomy division of the national academy of sciences of ukraine .
was also supported by the ukrainian state foundation for fundamental research under grant no .
a positive energy solution of the time - dependent dirac equation has a form @xmath263 with a two - component spinor @xmath35 satisfying the time - independent dirac equation ( [ dirac - eq ] ) : @xmath264\psi({\mathbf r})=0.\ ] ] where @xmath19 distinguishes two unitary inequivalent representations of @xmath265 gamma matrices .
it is convenient to denote the components of two - component spinor @xmath35 @xmath266 with the factor @xmath145 explicitly included in the definition of the lower component
. then we can rewrite dirac equation ( [ dirac - eq - append ] ) in the components as follows : @xmath267 since we consider a cylindrically symmetric configuration of the field with a vector potential @xmath268 , the system ( [ components - cartesian ] ) has to be rewritten in the polar coordinates @xmath269 : @xmath270 as discussed in sec .
[ sec : model ] , to analyze the problem with a singular at @xmath161 aharonov - bohm potential ( [ a - singular ] ) one has to do a self - adjoint extension of the dirac operator , see e.g. refs .  . to avoid this complication we consider a regularized field configuration ( [ profile ] ) suggested in refs .   with the profile function ( [ step - profile ] ) , so that @xmath271 from now on , we consider the specific case @xmath212 and seek for a solution of eq .
( [ components - polar ] ) in the following form @xmath272 for @xmath273 we obtain a system of the radial equations for a free dirac particle            [ 2nd - order : r < r ] @xmath277\psi_1(r)=0 , \label{psi1:r < r } \\ & \frac{d^2}{dr^2}\psi_2(r)+\frac{1}{r}\frac{d}{dr}\psi_2(r)- \left[\frac{m^2}{r^2}-\frac{e^2-\delta^2}{(\hbar v_f)^2}\right]\psi_2(r)=0 ,
\label{psi2:r < r}\end{aligned}\ ] ]      [ 2nd - order : r > r ] @xmath278\psi_1(r)=0 , \label{psi1:r > r } \\
\frac{d^2}{dr^2}\psi_2(r)+\frac{1}{r}\frac{d}{dr}\psi_2(r)- \left[\frac{(m+\eta)^2}{r^2}-\frac{e^2-\delta^2}{(\hbar v_f)^2}\right]\psi_2(r)=0 \label{psi2:r > r}\end{aligned}\ ] ]    for @xmath142 . the solutions of eqs .
( [ 2nd - order : r < r ] ) and ( [ 2nd - order : r > r ] ) are expressed in terms of the usual bessel functions .
for example , for the solutions of the equations ( [ psi1:r < r ] ) and ( [ psi1:r > r ] ) for the component @xmath167 are given by      where @xmath280 , @xmath281 , and @xmath282 are constants and the @xmath283 s are the bessel functions .
the solution ( [ solution - psi1:r < r ] ) is standard due to normalizability and the absence of delta function at @xmath161 .
the coefficients @xmath280 , @xmath281 , and @xmath282 to be found from the matching conditions ( [ continuity ] ) and ( [ discontinuity ] ) . to find the second component , @xmath168 one can substitute the result for @xmath167 in eq .
( [ system : r > r - for2 ] ) . or , equivalently
, one can start from eqs .
( [ psi2:r < r ] ) and ( [ psi2:r > r ] ) for the component @xmath168 which also have a solution in the form ( [ solution - psi1 ] ) , find the corresponding constants from the matching conditions and then use eq .
( [ system : r > r - for1 ] ) to obtain @xmath167 .
finally , the overall factor before the solution is determined by the normalization condition @xmath284