the past few years have witnessed significant progress in the building of realistic inflation models within string theory @xcite .
the most popular approach has used the interaction potential between branes and antibranes as the source of energy driving inflation ; hence the interbrane separation plays the role of the inflaton .
early attempts suffered from the assumption that moduli could be fixed by fiat , without specifying the stabilization mechanism .
however , explicit examples have shown that properly fixing the moduli can interfere with the flatness which is required for the inflaton potential @xcite ( ) .
flux compactifications in type iib string theory have provided a robust way of stabilizing the dilaton and complex structure moduli of the calabi - yau manifold @xcite , with the added benefit of providing warped throats @xcite which can be useful for tuning the scale of inflation or of the standard model .
nonperturbative effects have been plausibly invoked to stabilize the remaining khler moduli @xcite .
this last step , unfortunately , generically spoils whatever flatness of the inflaton potential which was accomplished by warping , and thus requires additional fine tuning of parameters .
one can introduce extra parameters associated with corrections to the superpotential , which must be tuned to a part in 100 in order to obtain the minimum amount of 60 e - foldings of inflation @xcite .
alternatively , it is possible to tune parameters already present within the model but at the level of a part in 1000 @xcite .
the fine - tuning problem of the  model has been considered by several authors @xcite .
shift symmetry of the lagrangian has been proposed as a possible solution @xcite .
it has also been pointed out that if lyman @xmath1 constraints are viewed with skepticism , the degree of tuning required is considerably relaxed @xcite . of course , this observation also applies to other models of inflation .    in this paper
we reexamine the fine - tuning issue , starting within the context of the original  potential , without introducing additional corrections .
a closer consideration of the potential between the brane and antibrane reveals the importance of a parameter ( which we call @xmath0 ) which was set to a certain value in .
it is the distance which the brane will travel before annihilating with the antibrane in the warped throat , assuming that the brane starts somewhere closer to the the interior of the calabi - yau at an initial position @xmath2 .
the other important parameters governing the shape of the inflaton potential are the warp factor of the throat , which we denote as @xmath3 , and the khler modulus at its stabilized value , @xmath4 .
the warp factor is tunable through the ratio of ns - ns and rr fluxes @xmath5 : @xmath6 @xcite .
the vev of the khler modulus @xmath4 is determined by the nonperturbative superpotential @xmath7 invoked by ref .
@xcite , where @xmath8 in the case of gluino condensation in an @xmath9 gauge theory sector which can naturally arise in these compactifications .
if the parameter @xmath0 is set to zero as in , then indeed it is necessary to introduce additional sources of dependence on the inflaton into the potential , in order to cancel the unwanted contribution to the inflaton mass coming from the khler modulus .
however , in this work we point out that for a range of values @xmath10 of the combination @xmath11 , it is possible to achieve 60 or more e - foldings of inflation without any additional corrections to the potential
. fortunately this restriction on @xmath12 applies only to the simplest  model .
we will show that the restriction disappears when one considers more general examples in which corrections to the khler and superpotentials are taken into account .
another crucial parameter which we consider is the number of mobile mobiles branes , @xmath13 , generalizing the simplest case of only a single brane playing the role of the inflaton .
we refer to this possibility as multibrane inflation . within the aforementioned range of @xmath12 values
, we observe a very interesting phenomenon , whereby starting with a sufficiently large number @xmath13 dynamically leads to a sufficiently flat inflaton potential . at large @xmath13 ,
the branes start out being confined in a metastable minimum of the potential .
but as successive branes tunnel out of this minimum , it becomes shallower , until finally all the remaining branes can roll together into the throat .
this is illustrated in figure [ fig1 ] .
the potential is often very flat when the branes start rolling , leading to a long period of inflation .
this is a novel , dynamical mechanism for flattening the inflaton potential , which appears to be quite specific to brane - antibrane inflation .
before looking at the dynamics of inflation with stacks of branes and antibranes , we consider under what conditions multiple branes are compatible with the stabilization of the extra dimensions , _
i.e. , _ the khler modulus , @xmath4 , in supergravity language . for each extra brane
we have to introduce a corresponding antibrane such that the tadpole cancellation condition is statisfied . adding antibranes lifts the local minimum of the potential by an additional amount , given by the sum of their tensions , reduced by the warp factor at their location in the bottom of the klebanov - strassler throat .
there is a critical number of antibranes that can be added such that the potential still has a local minimum along the @xmath4 direction and therefore the khler modulus remains stable .
the stronger the warping , the larger the number of antibranes that can be added . in our analysis
we will add antibranes such that the khler modulus remains stabilized and heavy .
thus it does not act like a second inflaton field , since it changes very little during the slow - rolling of the brane - antibrane separation @xmath14 .
not all models are suitable for constructing such a multi - brane setup .
we will give here two examples , one which works and one which does not .
the original kklt @xcite setup does not accommodate a large number of mobile branes ; just one additional antibrane is sufficient to lift the vacuum from anti - desitter ( ads ) to minkowski ( or desitter , ds ) .
the barrier that prevents tha khler modulus from rolling has a height roughly equal to the depth of the ads minimum .
addind a second brane lifts the potential by an amount that is roughly equal to the height of the barrier , which destroys the local minimum . in figure [ fig1n ]
we take the parameters of the potential and the antibrane tension as given by @xcite .
one sees that adding a second brane results in a very shallow minimum and three branes destabilizes the khler modulus .
the potential of the khler modulus is given by v = aae^-a ( 13aa e^-a + w_0 + a e^-a)+nd^3 with @xmath15 , @xmath16 , @xmath17 @xmath18 .
-0.25 cm     and @xmath19 antibranes , showing there can be at most 2 antibranes without destabilizing the khler modulus .
[ fig1n ] , scaledwidth=45.0% ]    the situation is dramatically improved if one uses a racetrack superpotential as in ref .
the minimum is already minkowski and adding antibranes in a highly warped throat lifts the vacuum to ds without destabilizing the khler modulus .
the potential is : & & v= e^-2(a+b)6 ^ 2(b b e^a+ a a e^b ) + & & + where we choose the same parameters as in @xcite , @xmath20 and @xmath21 . with these values
the critical number of branes turns out to be @xmath22 , illustrated in fig .
[ fig2n ]    , scaledwidth=45.0% ]    in the remaining analysis , we will not be concerned with the particular model of khler modulus stabilization , beyond the assumption that it admits the required number of additional branes .
it will also be found that a large number is not required , so the racetrack model mentioned here is sufficient to establish the existence of a viable model .
our starting point for the dynamics of inflation is the lagrangian for the mobile brane position @xmath14 in the  framework . to simplify the analysis
, we will assume that all moduli ( notably the khler modulus @xmath4 ) are heavy except for the inflaton , and that the 4d effective cosmological constant vanishes at the end of inflation . the lagrangian can be parametrized as [ lag ]
l = 6r^2 ^ 2 - tr^2(1 + ( -_0)^4)^-1 where @xmath23 @xcite , @xmath24 is related to the tension @xmath25 of the 3-brane or antibrane by @xmath26 @xcite , @xmath27 is the warp factor , @xmath14 is the 3-brane position , and @xmath0 is the location of the antibrane . to keep our work self - contained ,
this expression is derived in appendix a , in just the same way as in except for the omission of two approximations : ( 1 ) we do not assume that @xmath28 , and ( 2 ) we do not taylor - expand the term @xmath29 .
let us first comment on point ( 1 ) above .
in  and certain other investigations of brane - antibrane inflation , it was assumed that inflation takes place sufficiently far from the antibrane to justify the approximation @xmath30 . on the other hand ,
the @xmath0 parameter was kept in @xcite and @xcite .
physically , it may be important to allow for nonzero @xmath0 because of the possibility that the @xmath14-dependence coming from the @xmath31 factor in the potential will tend to stabilize the brane at some point other than the bottom of the warped throat , in the absence of the brane - antibrane interaction .
an example where this is particularly clear is that of two symmetrically placed warped throats , each with its own antibrane @xcite . in that case
the potential ( translated to our notation ) was taken to be v = tr^2(2 - ( -_0)^4 - ( + _ 0)^4 ) from symmetry considerations it is obvious in this model that there exists a metastable or unstable equilibrium point exactly midway between the two throats , at @xmath32 , which can not coincide with the bottom of either throat .
if we simply remove one of the throats , the result corresponds to our starting point ( [ lag ] ) .
the presence or absence of the @xmath0 parameter explains differences in the amount of required fine - tuning of the potential found in different investigations . to elucidate this ,
recall that the @xmath14-dependence in the volume modulus @xmath33 is the origin of the @xmath34 problem pointed out in the mass of the inflaton is of order @xmath35 in planck units , ruining the slow roll condition .
it was noticed by @xcite that it is possible to tune parameters of a model with a potential similar to that in ( [ lag ] ) so that the slow - roll conditions are satisfied , even without adding any extra @xmath14-dependent corrections to the superpotential , by varying only the parameters which already appear .
however , the degree of tuning needed was more severe than in where corrections to the superpotential were invoked .
one way of understanding this is to realize that if @xmath36 and only corrections of order @xmath37 are added , then the potential automatically has @xmath38 , so that one of the slow roll parameters @xmath39 is small without any tuning .
then only the @xmath34 parameter ( @xmath40 ) needs to be tuned . in our potential
however , both slow - roll parameters must be tuned , for generic parameters in the potential . in @xcite
it was found that 60 e - foldings of inflation could be obtained by tuning the warp factor to one part in 2000 , whereas  only needed one part in 100 tuning .
although this sounds contrary to the spirit of the present paper , it was due to the fact that no systematic study of the effect of varying @xmath0 was done in @xcite ; instead only one fixed value was considered .
we note that the behavior of the potential as @xmath41 is classically correct , since by derivation it is of the form @xmath42 where @xmath43 satisfies @xmath44 in the 6d calabi - yau manifold .
although the full brane position must be specified by 3 complex @xmath14 coordinates , we are keeping the angular ones fixed and following only the radial coordinate for the brane - antibrane separation .
we also comment on our choice of avoiding the second approximation ( point ( 2 ) above ) , by not taylor expanding @xmath45 .
we choose not to do so because there is no advantage in making this expansion
it only leads to an unphysical divergence in the potential energy as the brane - antibrane distance shrinks .
the unexpanded form has good behavior as @xmath46 , including the vanishing of the potential at the minimum .
we will be interested below in the generalization to more than one brane rolling .
the lagrangian for @xmath13 branes and @xmath13 antibranes is [ lagm ]
l = 6r^2_i=1^n_i^2 - ntr^2(1 + _ i=1^n(_i-_0)^4)^-1 where now the volume modulus is r = 2- _ i=1^n_i^2 later we will specialize to the situation where all the branes are coincident , @xmath47 . then the potential takes the form @xmath48 .
the optimal potential for achieving slow roll inflation is one where both the first and second derivatives vanish simultaneously at some critical point @xmath49 .
although this requires fine tuning , we will discuss a mechanism where the tuning can automatically occur through a dynamical process . to characterize the tuning , we introduce new parameters a & = & n_0 ^ 4 , b^2 = 2_0 ^ 6 , x = _ 0 such that the potential is proportional to [ hatv ] v(x ) = 1(b^2/ a - x^2)^2 ( 1 + a(x-1)^4)^-1 the conditions for the potential to have a flat point , as defined above , are [ crit ] a = ( 6x-1)(x-1)^4 , 5x^2(x-1)^4 = b^2    these equations can be solved analytically , as the solution of a cubic equation , but the expression is not very enlightening .
instead we give an approximate description , by first considering @xmath12 as a free parameter , and determining as a function of @xmath12 the critical values @xmath50 and @xmath51 which satisfy ( [ crit ] ) . for consistency
, we should also require that @xmath52 ; otherwise the potential blows up at @xmath53 , before the brane reaches the antibrane .
the range of values where @xmath54 ( obviously also required for a physically sensible solution ) and @xmath55 turns out to be rather narrow , [ range ] 0.259 < b < 0.286 , and in this region @xmath50 is to a good approximation a linear function of @xmath12 : [ ac ] a_c 3b - 0.77 one can also approximate the critical value @xmath51 by a quadratic form , @xmath56 .
we will show later that the narrowness of the interval ( [ range ] ) is an accident , in the sense that physically well - motivated perturbations to the potential can remove this restriction
. therefore we will ultimately not regard it as a fine tuning .
however the special value @xmath50 for a given value of @xmath12 clearly does represent a fine tuning of the parameters .
it is this condition which will be addressed by a dynamical mechanism .
namely , by starting with a large enough number @xmath13 of branes , one can initially have @xmath57 ( but still @xmath58 ) .
the potential always has a metastable minimum from which branes can tunnel in this case .
after each tunneling event , the value of @xmath59 diminishes by a discrete amount . at some critical number of branes , @xmath60 ,
the curvature changes so that the branes roll instead of being trapped .
it becomes a quantitative question , depending on the values of @xmath4 and @xmath27 , whether the potential is sufficiently flat to give enough inflation .
we will show that in fact it is almost always flat enough to give at least 60 e - foldings of inflation .    to find out how likely it is that adequate inflation will occur ,
we have explored the parameter space @xmath61 , subject to the aforementioned restrictions and the assumption that @xmath13 starts out being larger than the critical value , where all the branes start to roll . for each value of @xmath62 , the quantities @xmath0 and @xmath60 are determined .
the fields start from rest at the value where the potential has a local minimum with @xmath63 branes , on the assumption that the nontunneling branes do not move during the tunneling event .
the equations of motion are = -3h_i - v_i , = h for the canonical momentum _
i = 12r^2_i and number of e - foldings @xmath64 , where the hubble parameter is given by @xmath65 . since all the fields move together , we need consider only one equation for the inflatons , which we do numerically .
the result is that , for @xmath12 satisfying ( [ range ] ) , there is a high probability of getting enough inflation , for random values of @xmath66 selected from within a reasonable range , @xmath67 , @xmath68 .
scanning uniformly over this region , we find that @xmath69 of tries result in more than 55 e - foldings of inflation when @xmath70 , and an even greater percentage occurs for other values of @xmath12 .
the results are illustrated for @xmath70 in figure [ fig2 ] , showing the number of e - foldings , @xmath71 ( omitting values where @xmath72 for clarity ) versus @xmath73 and @xmath74 .    in figure [ fig2 ] it can be seen that @xmath71 is constant along lines of constant @xmath75 .
this can be understood by rewriting the potential in terms of the canonically normalized field @xmath76 : v = tna^3 ( 1 - 16 ^ 2)^-2 ( 1 - 36(n / a)^3(-_0)^4)^-1 where @xmath77
. thus for fixed @xmath12 , the shape of the potential , power spectrum , and duration of inflation depend only upon combination of parameters in @xmath78 , and not on the orthogonal combination @xmath79 .
we can therefore restrict our exploration of parameter space to a one - dimensional scan along the direction of @xmath78 .
moreover , we can show that the number of branes @xmath13 depends in a very simple way on the value of @xmath78 : @xmath80 .
this follows from eliminating @xmath0 from the definitions of @xmath59 and @xmath12 .
therefore the critical number of branes is given by [ nc ] n_c = a_c b^-4/3 ( 2)^2/3 ^ 1/3 and for a fixed value of @xmath12 , we can treat the parameters @xmath78 and @xmath60 interchangeably .
of course ( [ nc ] ) ca nt be satisfied exactly , because @xmath60 has to be an integer .
the accidental discrepancies between the optimal value of @xmath60 and the closest integer below this value will affect how long inflation actually lasts .
= 0.5    in figure [ fig3 ] we show the correlation between number of e - foldings and the critical number of branes , scanning uniformly in the interval @xmath81 of parameter space .
although the two quantities are generally correlated , there are order - of - magnitude variations in @xmath71 over narrow ranges of @xmath60 .
figure [ fig3 ] demonstrates that both features can be matched fairly well by an empirical relationship .
let @xmath82 be the discrepancy between the ideal value of @xmath60 and the closest integer below this value .
then we find that [ emp ] n_e 10 n_c the factor @xmath83 , is understandable since if @xmath84 , the potential is tuned to be perfectly flat at some position @xmath85 .
= 0.5    the factor @xmath86 in ( [ emp ] ) can be partly understood in terms of assisted inflation @xcite . with @xmath13 fields rolling simultaneously
, the kinetic term is multiplied by @xmath13 relative to one field .
we should rescale the fields by @xmath87 , and the slow - roll parameters by @xmath88 .
one therefore expects an enhancement in the length of inflation by some power of @xmath60 .
it would of course be satisfying to understand why that power is @xmath89 in the present model .
we have not yet succeeded in analytically deriving this result .
however it is not surprising that @xmath71 fails to be enhanced by a full factor of @xmath60 as would be the case for ordinary assisted inflation , because the potential contains @xmath90 in the khler modulus .
only if the potential was independent of @xmath13 would we be sure that @xmath91 . in that case , transforming to the canonically normalized field would cause @xmath92 and the slow roll parameters would go like @xmath93 .
the surprise is that there is any assisted inflation effect at all .
the reason for it is that at the critical value of @xmath13 , the potential is flatter than it would generically be for an arbitrary function of @xmath90 .
the above results were for @xmath70 .
allowing @xmath12 to vary , ones finds that smaller @xmath12 tends to give more inflation with a smaller number of branes .
this dependence is shown in figure [ fig4 ] .
the main observation , though , is that the great majority of all cases yield more than @xmath94 e - foldings of inflation .
in addition to having a long enough period of inflation , we must obtain a spectrum of density perturbations in agreement with observational constraints from the cosmic microwave background . to compute the power spectrum , we use the following expression @xcite which applies to a set of noncanonically normalized fields with lagrangian @xmath95 : p(k ) = v75 ^ 2 m_p^2 n_e^i n_e^j g^ij = n_c v h^275 ^ 2 m_p^2 r6 ^ 2 these quantities
are evaluated at horizon crossing of the relevant wave number , when @xmath96 . to normalize the spectrum
, we evaluate @xmath97 for modes which crossed the horizon 55 e - foldings before inflation ended , corresponding to a scale of inflation around @xmath98 gev
. the consistency of this assumption will be verified .
the cobe normalization implies that @xmath99 for these modes .
this allows us to normalize the scale of the potential by fixing the value of the parameter @xmath24 which is proportional to the unwarped brane tension .
the results are expressed in terms of the scale of inflation relative to the planck scale . in figure [ fig5 ]
we plot @xmath100 as a function of the number of branes .
( in the previous section it was explained that , just like @xmath71 , spectral properties depend only on the same combination of @xmath4 and @xmath27 which determine @xmath60 . ) from the figure one sees that the inflationary scale ranges between @xmath101 and @xmath102 , roughly consistent with our assumption of @xmath103 e - foldings of inflation .
the string scale , which determines the brane tension , must be in about the same range .
if the number of branes is small enough , the inflationary scale can be sufficiently high for the tensor contribution to the cosmic microwave background to be observable in future experiments .
= 0.5    the scalar spectral index is given by @xmath104 .
we have numerically evaluated it over the relevant range of parameters . like the other quantities ,
it is most conveniently displayed as a function of the number of branes .
the deviation @xmath105 at 55 e - foldings before the end of inflation is shown as a function of @xmath60 in figure [ fig6 ] .
there is a wide range of possible values , @xmath106 .
= 0.5    it is also interesting to consider the running of the spectral index , @xmath107 , for which wmap has given a hint of a larger - than - expected negative value , @xmath108 @xcite .
although the error bars are large enough to be consistent with no running , the suggestion of large running has attracted attention .
generically , within the slow - roll formalism , one expects @xmath109 , which is usually smaller than the measured central value .
our numerical survey shows that it is possible to have rare instances of running as negative as @xmath110 , even though inflation lasts long enough .
however , as shown in figure [ fig7 ] , it is still true that @xmath109 .
the rare cases of large running correspond to occasional extreme values of the spectral index .
we have checked that the spectral index and running are in good agreement with the predictions of slow roll inflation , @xmath111 , @xmath112 , where @xmath113 , @xmath114 and @xmath115 , in units where @xmath116 .
the derivatives must be taken with respect to the canonically normalized fields , @xmath117 . in all cases , @xmath118 , whereas @xmath119 and @xmath120 .
the running is completely dominated by @xmath121 .
= 0.5    it is interesting to compare the predicted values of the running and the spectral index with experimental constraints from combined wmap and sloan digital sky survey lyman @xmath1 data @xcite the allowed region is superimposed on the predictions in figure [ fig8 ] . although most of the predicted points are allowed , we note that the few examples with large running have too large a deviation of @xmath122 from unity . thus if large running should be experimentally confirmed in the future , it will rule out the models under consideration .
the previous analysis shows that within a narrow range of the parameter @xmath123 , a long period of inflation will always ensue , as long as the initial number of mobile branes exceeded the critical number @xmath60 .
but since @xmath12 looks fine - tuned , this is not a complete solution to the problem of naturalness for the inflaton potential .
could this be an accident of the particular form of the simplest potential ? to answer this , we explore the effect of perturbing the model by corrections which are expected to be present , but which we have so far ignored for the sake of simplicity .    at least two kinds of modifications to the basic  model are generically expected : corrections to the superpotential and to the khler potential .
a simple example of the former was considered in @xcite , where the superpotential for the khler modulus @xmath124 and inflaton was taken to be w = w_0 + g(t)(1 + ^2 ) leading to a contribution to the potential of the form [ deltav ] v = + ^2 + ^4 ( 2- ^2)^2 with @xmath4-dependent coefficients @xmath125 .
this should be added to the brane - antibrane potential which we have been using so far . in the case of @xmath13 branes , we take @xmath126 . to simplify the analysis , we consider examples where @xmath125 are chosen in such a way that the full potential continues to have a global minimum at @xmath127 , where @xmath128 is assumed to vanish : v ~n we have included the factors of @xmath13 appropriate to the multibrane case , when all branes are coincident .
the second kind of correction involves the khler potential , which for @xmath13 branes was taken to be @xmath129 @xcite .
this was always understood to be an approximate expression , good when @xmath130 , but subject to higher order corrections in @xmath131 , for example @xmath132 .
the denominator in @xmath133 would then be replaced by ( 2- n^2)^2 ( 2- n(^2 + d^4))^2 as shown in appendix b , the khler corrections also give contributions to the numerator of the f - term potential , whose taylor series starts with a term of order @xmath134 . for small @xmath14 , we can absorb this leading term into the @xmath135 coefficient of ( [ deltav ] ) .
the kinetic term @xmath136 also gets @xmath137 corrections , which can be inferred from eq.([kmetric ] ) , but these are negligible whenever the branes are slowly rolling .    by the above reasoning , the two kinds of corrections enter into the rescaled potential ( [ hatv ] ) as [ newhatv ] v ( c(x^2 - 1)^2 + ( 1 + a(x-1)^4)^-1 ) by experimenting with this generalized potential , one finds that the new coefficient @xmath138 allows for the removal of the lower bound on @xmath12 in ( [ range ] ) .
moreover , negative values of @xmath139 can remove the upper bound .
it is easy to find values of @xmath138 and @xmath139 for which the potential has the desired behavior as a function of @xmath13 , as in figure [ fig1 ] , regardless of the value of @xmath12 .
for example @xmath140 , @xmath141 yields the potentials shown in figure [ fig9 ] , at @xmath142 and @xmath143 .
it is important that the new parameters @xmath138 and @xmath139 are independent of the number of branes .
thus @xmath59 continues to be the only parameter which changes due to a brane tunneling event .
= 0.5    as in the previous section , we scanned over parameters for which @xmath78 varies , hence the number of branes , while keeping @xmath144 fixed , for the two examples with @xmath142 and @xmath143 .
figure [ fig10 ] shows the number of e - foldings versus number of branes .
the general growth of @xmath145 is evident , as it was for the simple model of the previous section .
= 0.5    figure [ fig11 ] shows the correlation of @xmath146 with @xmath122 .
it shows that the small-@xmath12 model tends to have a red spectrum , with @xmath147 , and negligible running .
the large-@xmath12 model also has @xmath147 for most realizations , but for rare cases with @xmath148 , it is possible to get large deviations from a flat spectrum with @xmath149 .
these are accompanied by sizable running , but as in figure [ fig8 ] , such points occur outside of the experimentally allowed region .
the overall conclusion is that the generalized models have similar qualitative behavior to the simplest class which was studied in sections iii - iv .
in this work we have presented a very general means whereby a long period of inflation , with a suitably flat spectrum of density perturbations , can be automatically obtained after tunneling of branes , from a metastable position in a calabi - yau manifold , into a warped throat containing antibranes .
it depends only on the fact that the brane - antibrane potential generically has such a metastable minimum when the number of branes is sufficiently large , and the possibility of starting at a place in the string theory landscape where the number of mobile branes exceeds the threshold for having a minimum .
once a certain critical number of branes @xmath60 is reached , the minimum disappears and the branes start rolling slowly to the true minimum .
they automatically start from the same position , which is guaranteed to be close to the flattened region of the potential .    in this scenario , an arbitrarily long time will be spent in a metastable de sitter phase resembling old inflation , before the slow rolling period ( similarly to ref .
each tunneling event results in a small bubble of open universe , within which another period of quasi - de sitter expansion will occur .
the process repeats itself until the critical number of branes is reached . for us , these long periods of old inflation are irrelevant  we have focused entirely on the final phase of conventional inflation associated with the smooth motion of the branes until they annihilate with antibranes in the throat , hopefully leading to the reheating of the universe @xcite , and perhaps the production of cosmic superstrings @xcite .
we numerically scanned over large ranges of the logarithm of relevant combinations of parameters which govern the shape of the inflaton potential , and found that large fractions ( @xmath150 ) of the trials resulted in inflation with no need for tuning , whenever the flattening mechanism could occur .
this seems to be a quite promising approach for ameliorating the need for fine - tuning within the framework of the string - theory landscape @xcite .    for convenience of numerical integration , we chose a form of the potential which has good behavior as the branes reach the bottom of the throat .
this form , which differs from the commonly used taylor - expanded form , is in fact the result of a direct perturbative computation , and does not require any special assumptions .
however , the perturbative derivation does assume that @xmath151 in order to be valid .
one might therefore wonder how sensitive our results are to the behavior of the potential when this approximation is breaking down .
we have investigated what happens if the taylor expansion @xmath152 is used in place of @xmath153 , and we also checked different ways of regulating the singularity , like @xmath154 .
for all these alternatives , we found that although the specific parameter values where a flat potential occurs may be somewhat sensitive to the behavior of the potential near @xmath0 , the mechanism of multibrane inflation nevertheless works in the same way as we have presented .
we therefore believe it to be a robust result , which is qualitatively quite different from field theory models of inflation .
although the dynamical flattening of the potential requires a rather special range of parameters in the simplest  model , we showed that this restriction can be lifted once corrections to the superpotential and the khler potential are taken into account .
we exhibited an example of such corrections that works , without systematically exploring the parameter space of such corrections .
subjectively , it seems that the parameter space ( referred to as @xmath155 in section v ) where multibrane inflation works will be large , given that it was easy to find a successful example , but a more complete investigation could be made .
moreover , such corrections should in principle be computed from a fully string theoretic starting point , which might restrict the allowed ranges of @xmath155 .
another issue which deserves further study is the precise relationship between the radial coordinate @xmath156 of the brane within the warped throat , and the field @xmath14 which is used to describe the inflaton position in the supergravity description . the exact correspondence between these fields was not needed in because the main dependence on @xmath14 was taken to be through the khler modulus factor ( and superpotential corrections ) , rather than a combination of these with the brane - antibrane potential .
we have simply identified the two fields in our analysis , but the actual relation could be more complicated , possibly introducing additional parameters into the lagrangian .
we thank c.  burgess and f.  quevedo for collaboration during the early stages of this work .
we also thank h.  firouzjahi , s.  kachru , r.  kallosh , l.  kofman , a.  linde , i.  neupane , e.  silverstein , l.  susskind , and h.  tye for helpful comments and discussions .
jc and hs are supported by grants from nserc of canada and fqrnt of qubec .
we consider the background of @xmath157 space with the @xmath158 metric : ds_5 ^ 2 = ( -dt^2+d*x*^2 ) + dr^2 where @xmath159 is the radius of the @xmath158 space created as the near - horizon geometry of a stack of @xmath160 @xmath161 branes .
there is also an @xmath162 background field created by the stack of branes : a_x^0 , x^1 , x^2 , x^3 = _
r0()^4 if one now writes the action for a probe brane is this background , one finds : s_d3 & = & -t_3d^4 + & + & t_3d^4a_x^0 , x^1 , x^2 , x^3 we will use the convention that @xmath163 , but the rest of the coordinates will describe the location of the brane , including fluctuations ( local bendings ) , @xmath164 .
we make here a further assumption , namely that only its radial location varies .    _
^ax^m_^bx^ng_nm = _ ^ax^0_^bx^0g_00 + _ ^ax^i_^bx^ig_ii + _ ^ar_^brg_rr the metric will look like : (    cccc g_00 + g_rr()^2 & g_rr & g_rr & g_rr + g_rr & g_xx + g_rr()^2 & g_rr & g_rr + g_rr & g_rr & g_yy + g_rr()^2 & g_rr + g_rr & g_rr & g_rr & g_zz + g_rr()^2    ) = (    cccc g_00 & 0 & 0 & 0 + 0 & g_xx & 0 & 0 + 0 & 0 & g_yy & 0 + 0 & 0 & 0 & g_zz    ) + ( g_rr ) we can apply the usual formula : ( a+b ) = a(1+a^-1b ) a ( 1 + tr ( a^-1b ) ) writing the metric at the location of the brane as : @xmath165 we obtain : ( _ ^ax^m_^bx^ng_nm)= -()^8(1-()^4 f^ )    the factor @xmath166 in the bracket comes from the product g^g_rr = ( ) ^4f^ using these results in the @xmath167 action we obtain : [ dbi+cs ] s_d3 & = & -t_3d^4()^4 + & & t_3d^4 ( ) ( ) ^4_t^1 ^ 2 ^ 3 expanding the square root term we obtain for a @xmath161 brane ( sign @xmath168 ) : s_d3 = t_3d^4 f^ while for a @xmath169 ( sign @xmath170 ) : [ d3_bar_action ] s_3 & = & t_3d^4 f^ + & - & 2t_3d^4()^4 the extra piece in the case of the @xmath169 comes from the non - bps nature of the antibrane , since the graviton exchange and the rr exchange forces now add instead of canceling ( for 3-branes the dilaton is constant and contributes no force ) .
the potential energy for the @xmath169 in the gravitational + rr field will make it favorable for the brane to move towards @xmath171 , that is , the bottom of the @xmath158 throat . for the setup in fig.([kklmmt_setup ] ) we have to add the piece corresponding to the graviton + rr exchange between the @xmath161 and the @xmath169 .    ,
scaledwidth=40.0% ]    however in this setup we have to include the effect of the @xmath161 in addition to that of the stack
. the derivation of the background perturbed by the presence of the @xmath161 placed away from the stack can be found in appendix b of .
it starts by writing the metric in the form : ds_10 ^ 2 & = & h^-(-dt^2+d*x*^2)+ h^(dr^2+_abdy^ady^b ) + h&= & placing the extra @xmath161 brane at @xmath172 and solving for the perturbed harmonic function h(r)= + h(r ) which satisfies the laplace equation in the space transverse to the brane , which is spanned by @xmath156 and @xmath173 : ^2_(6)h(r ) = c^(6 ) ( - _ 1 ) which amounts to solving : _
( - _ 1 ) which give @xmath174 .
a similar equation is satisfied by the unperturbed @xmath175 .
the constants @xmath176 and @xmath177 depend on the total `` charge '' , more precisely the total tension of the brane or stack of branes , so they are related by : c = and since the solution for the unperturbed @xmath178 is already known , we can immediately find the solution for @xmath179 h(r ) = = we can use the perturbed metric to calculate the action for the @xmath169 using the new function @xmath180 .
now modify the setup and place at @xmath181 a stack of @xmath182 @xmath169 s with @xmath183 and an equal number of @xmath161 s placed at @xmath184 , as shown in figure [ modified_setup ] .
we can simply generalize the analysis above to come to the result : h(r ) = + _ j=1^m there is no interaction between the @xmath161 s or between the @xmath169 s in the second stack since they are mutually bps .
we can now go back to eq .
( [ dbi+cs ] ) and replace the expression @xmath185 by @xmath186 .
the relevant term is the second term from eq .
( [ d3_bar_action ] ) , which in the limit @xmath187 becomes :    2t_3d^4 = 2t_3d^4 d^4\{2t_3(1-_j=1^m ) }    where @xmath181 is the location of the stack of antibranes .
the above result is for one single @xmath169 in the stack .
the final form of the effective action should be multiplied by @xmath182 , the number of antibranes in the stack .
, scaledwidth=40.0% ]
the khler potential is assumed to have the following form : k=-3(r)= -3(t+-f(_i,_i ) ) and we will also assume that @xmath188 decomposes into a sum of functions , each one depending on one field @xmath189 : [ kpot ] k=-3(t+-_if_i(_i,_i ) ) the simplest form of the functions @xmath190 is f_i(_i,_i ) = _ i_i but we can consider more complicated expressions , for example : f_i(_i,_i ) = _ i_i - ( _ i_i)^2    allowing for the more complicated expression for @xmath191 changes the structure of the singularities of the khler potential and of the f - term . in the general case
the khler metric derived from the khler potential of eq .
( [ kpot ] ) has the form :              the f - term is calculated using the usual formula : v_f = e^k(k^bd_bw-3|w|^2 ) where = -3(r ) the covariant derivatives are expressed as : d_tw & = & _ tw+w_tk = _ tw - w + d__iw & = & _ _
iw+w__ik = w + & = & _ + _ k = _ - + & = & _ + _ k =    using the above expression for the inverse khler metric one obtains the following terms : k^td_tw & = & ( _ - ) ( _ tw - w ) ( r+ +  + ) + k^_id__iw & = & ( _ -)w ( +  + ) + k^td_tw & = & ( _ tw - w ) ( +  + ) + k^_jd__jw & = & w ( +  + )          g.  r.  dvali and s.  h.  h.  tye , `` brane inflation , '' phys .  lett .
b * 450 * , 72 ( 1999 ) [ arxiv : hep - ph/9812483 ] .
c.  p.  burgess , m.  majumdar , d.  nolte , f.  quevedo , g.  rajesh and r.  j.  zhang , `` the inflationary brane - antibrane universe , '' jhep * 0107 * , 047 ( 2001 ) [ arxiv : hep - th/0105204 ] .
j.  garcia - bellido , r.  rabadan and f.  zamora , `` inflationary scenarios from branes at angles , '' jhep * 0201 * , 036 ( 2002 ) [ arxiv : hep - th/0112147 ] .
c.  p.  burgess , p.  martineau , f.  quevedo , g.  rajesh and r.  j.  zhang , `` brane antibrane inflation in orbifold and orientifold models , '' jhep * 0203 * , 052 ( 2002 ) [ arxiv : hep - th/0111025 ] .
r.  blumenhagen , b.  kors , d.  lust and t.  ott , `` hybrid inflation in intersecting brane worlds , '' nucl .
b * 641 * , 235 ( 2002 ) [ arxiv : hep - th/0202124 ] .
n.  jones , h.  stoica and s.  h.  h.  tye , `` brane interaction as the origin of inflation , '' jhep * 0207 * , 051 ( 2002 ) [ arxiv : hep - th/0203163 ] .
k.  dasgupta , c.  herdeiro , s.  hirano and r.  kallosh , `` d3/d7 inflationary model and m - theory , '' phys .
d * 65 * , 126002 ( 2002 ) [ arxiv : hep - th/0203019 ] .
o.  dewolfe , s.  kachru and h.  l.  verlinde , `` the giant inflaton , '' jhep * 0405 * , 017 ( 2004 ) [ arxiv : hep - th/0403123 ] . j.  j.  blanco - pillado _ et al . _ ,
`` racetrack inflation , '' jhep * 0411 * , 063 ( 2004 ) [ arxiv : hep - th/0406230 ] .
e.  silverstein and d.  tong , `` scalar speed limits and cosmology : acceleration from d - cceleration , '' phys .
d * 70 * , 103505 ( 2004 ) [ arxiv : hep - th/0310221 ] .
m.  alishahiha , e.  silverstein and d.  tong , `` dbi in the sky , '' phys .
d * 70 * , 123505 ( 2004 ) [ arxiv : hep - th/0404084 ] .
x.  g.  chen , `` multi - throat brane inflation , '' phys .
d * 71 * , 063506 ( 2005 ) [ arxiv : hep - th/0408084 ] ; `` inflation from warped space , '' arxiv : hep - th/0501184 .
k.  becker , m.  becker and a.  krause , `` m - theory inflation from multi m5-brane dynamics , '' nucl .
b * 715 * , 349 ( 2005 ) [ arxiv : hep - th/0501130 ] .
s.  kachru , r.  kallosh , a.  linde , j.  maldacena , l.  mcallister and s.  p.  trivedi , `` towards inflation in string theory , '' jcap * 0310 * , 013 ( 2003 ) [ arxiv : hep - th/0308055 ] .
s.  b.  giddings , s.  kachru and j.  polchinski , `` hierarchies from fluxes in string compactifications , '' phys .
d * 66 * , 106006 ( 2002 ) [ arxiv : hep - th/0105097 ] . i.  r.  klebanov and m.  j.  strassler , `` supergravity and a confining gauge theory : duality cascades and @xmath194sb - resolution of naked singularities , '' jhep * 0008 * , 052 ( 2000 ) [ arxiv : hep - th/0007191 ] .
s.  kachru , r.  kallosh , a.  linde and s.  p.  trivedi , `` de sitter vacua in string theory , '' phys .
d * 68 * , 046005 ( 2003 ) [ arxiv : hep - th/0301240 ] .
de alwis , `` effective potentials for light moduli , '' arxiv : hep - th/0506266 .
c.  p.  burgess , j.  m.  cline , h.  stoica and f.  quevedo , `` inflation in realistic d - brane models , '' jhep * 0409 * , 033 ( 2004 ) [ arxiv : hep - th/0403119 ] .
m.  berg , m.  haack and b.  kors , `` loop corrections to volume moduli and inflation in string theory , '' phys
.  rev .
d * 71 * , 026005 ( 2005 ) [ arxiv : hep - th/0404087 ] .
j.  p.  hsu , r.  kallosh and s.  prokushkin , `` on brane inflation with volume stabilization , '' jcap * 0312 * , 009 ( 2003 ) [ arxiv : hep - th/0311077 ] .
j.  p.  hsu and r.  kallosh , `` volume stabilization and the origin of the inflaton shift symmetry in string theory , '' jhep * 0404 * , 042 ( 2004 ) [ arxiv : hep - th/0402047 ] .
h.  firouzjahi and s.  h.  h.  tye , `` closer towards inflation in string theory , '' phys .
b * 584 * , 147 ( 2004 ) [ arxiv : hep - th/0312020 ] .
h.  firouzjahi and s.  h.  tye , `` brane inflation and cosmic string tension in superstring theory , '' jcap * 0503 * , 009 ( 2005 ) [ arxiv : hep - th/0501099 ] .
r.  kallosh and a.  linde , `` landscape , the scale of susy breaking , and inflation , '' jhep * 0412 * , 004 ( 2004 ) [ arxiv : hep - th/0411011 ] .
o.  dewolfe and s.  b.  giddings , `` scales and hierarchies in warped compactifications and brane worlds , '' phys .  rev .
d * 67 * , 066008 ( 2003 ) [ arxiv : hep - th/0208123 ] .
s.  b.  giddings and a.  maharana , `` dynamics of warped compactifications and the shape of the warped landscape , '' arxiv : hep - th/0507158 .
n.  iizuka and s.  p.  trivedi , phys .
d * 70 * , 043519 ( 2004 ) [ arxiv : hep - th/0403203 ] .
a.  r.  liddle , a.  mazumdar and f.  e.  schunck , `` assisted inflation , '' phys .  rev .
d * 58 * , 061301 ( 1998 ) [ arxiv : astro - ph/9804177 ] .
a.  mazumdar , s.  panda and a.  perez - lorenzana , `` assisted inflation via tachyon condensation , '' nucl .
b * 614 * , 101 ( 2001 ) [ arxiv : hep - ph/0107058 ] .
y.  s.  piao , r.  g.  cai , x.  m.  zhang and y.  z.  zhang , `` assisted tachyonic inflation , '' phys .
d * 66 * , 121301 ( 2002 ) [ arxiv : hep - ph/0207143 ] .
m.  sasaki and e.  d.  stewart , `` a general analytic formula for the spectral index of the density perturbations produced during inflation , '' prog .
* 95 * , 71 ( 1996 ) [ arxiv : astro - ph/9507001 ] .
h.  v.  peiris _ et al .
_ , `` first year wilkinson microwave anisotropy probe ( wmap ) observations : implications for inflation , '' astrophys .
j.  suppl .
* 148 * , 213 ( 2003 ) [ arxiv : astro - ph/0302225 ] .
u.  seljak _ et al .
_ , `` cosmological parameter analysis including sdss ly - alpha forest and galaxy bias : constraints on the primordial spectrum of fluctuations , neutrino mass , and dark energy , '' phys .  rev .
d * 71 * , 103515 ( 2005 ) [ arxiv : astro - ph/0407372 ] .
n.  barnaby and j.  m.  cline , `` creating the universe from brane - antibrane annihilation , '' phys .
d * 70 * , 023506 ( 2004 ) [ arxiv : hep - th/0403223 ] .
n.  barnaby , c.  p.  burgess and j.  m.  cline , `` warped reheating in brane - antibrane inflation , '' jcap * 0504 * , 007 ( 2005 ) [ arxiv : hep - th/0412040 ] .
l.  kofman and p.
yi , `` reheating the universe after string theory inflation , '' arxiv : hep - th/0507257 .
s.  sarangi and s.  h.  h.  tye , `` cosmic string production towards the end of brane inflation , '' phys .
b * 536 * , 185 ( 2002 ) [ arxiv : hep - th/0204074 ] .
n.  t.  jones , h.  stoica and s.  h.  h.  tye , `` the production , spectrum and evolution of cosmic strings in brane inflation , '' phys .  lett .
b * 563 * , 6 ( 2003 ) [ arxiv : hep - th/0303269 ] .
l.  pogosian , s.  h.  h.  tye , i.  wasserman and m.  wyman , `` observational constraints on cosmic string production during brane inflation , '' phys .
d * 68 * , 023506 ( 2003 ) [ arxiv : hep - th/0304188 ] .
e.  j.  copeland , r.  c.  myers and j.  polchinski , `` cosmic f- and d - strings , '' jhep * 0406 * , 013 ( 2004 ) [ arxiv : hep - th/0312067 ] . t.  damour and a.  vilenkin , `` gravitational radiation from cosmic ( super)strings : bursts , stochastic background , and observational windows , '' phys .  rev .
d * 71 * , 063510 ( 2005 ) [ arxiv : hep - th/0410222 ] .
n.  barnaby , a.  berndsen , j.  m.  cline and h.  stoica , `` overproduction of cosmic superstrings , '' jhep * 0506 * , 075 ( 2005 ) [ arxiv : hep - th/0412095 ] .
l.  susskind , `` the anthropic landscape of string theory , '' arxiv : hep - th/0302219 .
t.  banks , m.  dine and e.  gorbatov , `` is there a string theory landscape ?
, '' jhep * 0408 * , 058 ( 2004 ) [ arxiv : hep - th/0309170 ] .


dense plasmas have recently gained growing interest due to their relevance for the interior of giant planets as well as for laser interaction with matter and inertial confinement fusion scenarios .
examples of recent experimental studies include the ultrafast thermalization of laser plasmas @xcite or free electron laser excited plasmas @xcite , inertial confinement fusion experiments at the national ignition facility @xcite and magnetized z - pinch experiments at sandia @xcite .
questions of fundamental theoretical importance are the conductivity and heat conduction , the energy loss of energetic particles ( stopping power ) in such a plasma , e.g. @xcite or the temperature equilibration of the electronic and ionic components @xcite .    despite recent advances in modeling and computer simulations
a fully selfconsistent treatment of these , in general , highly nonequilibrium electron - ion plasmas has not been possible so far due to the requirement of the simultaneous account of electronic quantum and spin effects together with the ( possibly ) strong ionic correlations .
the main problem here are the vastly different time scales of electrons and ions resulting from their different masses .
a possible solution of this dilemma is a multi - scale approach that has been proposed by ludwig _
it takes advantage of the weak electron - ion coupling that allows for a linear response treatment of the electrons .
this idea has been used by graziani _
_ to decouple the electron kinetic equation using an stls ( singwi - tosi - land - sjlander ) scheme @xcite or a recently derived extension @xcite .
the key of this multiscale approach is to absorb the fast electron kinetics into an effective screened potential @xmath0 of the heavy ions with charge @xmath2 where the screening is provided by the electrons via a proper dielectric function @xmath3 , e.g. @xcite @xmath4 taken in the static limit .
we note that also the effect of streaming electrons maybe important in warm dense matter
. then the frequency argument becomes @xmath5 which leads to wake effects , that are well investigated theoretically and experimentally for dusty plasmas , e.g. @xcite and quantum plasmas , e.g. ref .
@xcite and references therein . however , this is beyond the scope of this paper .
a similar but even simpler approach that is based on a quantum hydrodynamic model ( qhd ) and replaces the linear response potential of the ions by a simplified expression that is derived from linearized qhd ( lqhd ) , @xmath6 . the first expression for @xmath7
was derived by shukla and eliasson ( se ) who predicted an attractive interaction between ions , even in the absence of streaming , @xmath8 @xcite .
comparisons with density functional theory revealed that this is incorrect @xcite .
this also underlined the limitations of qhd models to weakly correlated dense plasmas @xcite .
a recent overview and more references can be found in ref .
@xcite .
recently , modified expressions for @xmath7 were derived by akbari
moghanjoughi ( am ) @xcite and stanton and murillo ( sm ) @xcite .
the latter is particularly interesting because it is also applicable to finite temperatures , in contrast to the se and am potentials that neglect thermal effects .
the question arises whether the am and sm potentials are more accurate than the se potential and whether they allow for an extension of the applicability limits of qhd .
it is the goal of the present paper to perform such an analysis . to quantify the accuracy of these three potentials we use , as a basis , the standard static potentials of weakly correlated electrons that can be derived from quantum kinetic theory @xcite by linearization in the external perturbation
this directly leads to the static potential ( [ pot_stat ] ) involving the random phase approximation ( rpa ) for the dielectric function and the mermin dielectric function , respectively .
both include kinetic effects , finite temperature and , in the latter case , also collisions .
being rooted in a kinetic approach , the corresponding static potentials are , by construction , more accurate than the potentials @xmath7 and , thus , allow for a rigorous test of the validity of the latter .
the main conclusions of our analysis are the following : 1 . )
the se potential is qualitatively different from the kinetic results . 2 .
) the sm potential exhibits very good agreement with the rpa .
the am potential is in agreement with the sm potential , at @xmath9 , differing only in the notation and 4 . )
the sm potential is in good agreement with the rpa result even at elevated temperature .
this confirms the correctness of the recently reported @xcite correction factor @xmath1 that has to be included in front of the bohm term of the qhd equations for fermions .
the paper is organized as follows : we first recall the rpa and mermin dielectric functions in sec .  [
s : mermin ] and discuss their basic properties such as the reproduction of friedel oscillations .
after this we recall the se , am and sm potentials in sec .
[ s : potential ] and present numerical results that compare them to the potential @xmath0 computed from the rpa and mermin dielectric functions .
the paper concludes with a brief discussion of the physical origin of the factor @xmath1 .
the main steps of its derivation are outlined in an appendix .
inset shows an enlargement of the potential around its zero values to resolve the friedel oscillations . ]    the mean - field result for the electron dielectric function is given by the random phase ( rpa or hartree or quantum vlasov ) approximation , @xmath10 , where correlation effects are neglected , and damping of collective oscillations is entirely due to landau damping .
temperature effects are straightforwardly included by using the corresponding fermi function .
for an overview on the final expression , see the appendix of ref .
@xcite . the effective ion potential ( [ pot_stat ] ) computed from the rpa dielectric function in the zero temperature limit
is shown in figs .
[ fig : rpa_friedel ] and [ fig : rpa_friedel2 ] for different densities .
the zero temperature rpa potential reproduces the friedel oscillations , as is illustrated in fig .
[ fig : rpa_friedel2 ] by comparison to the large @xmath11 asymptotics , @xmath12 , ( red dashed line ) .
while the oscillations are very shallow and hardly of relevance for the thermodynamics and transport properties in warm dense matter , the reproduction of the correct periodicity is a useful consistency and accuracy test of the fourier transformation in eq .
( [ pot_stat ] ) .    for later reference we also provide the static long wavelength limit of the effective ion potential which is nothing but the yukawa potential ( [ pot_stat ] ) @xmath13 with the familiar inverse yukawa screening length , @xmath14 , that interpolates between the debye and thomas - fermi expressions , in the non - degenerate and zero temperature limits , respectively , @xmath15 here
@xmath16 is the thomas - fermi wave number , and @xmath17 is the fermi integral of order @xmath18 . note that the potential ( [ yukawa ] ) depends  via @xmath19on density and temperature .    in the warm dense matter regime
we expect the electrons to be weakly to moderately coupled so that correlation effects are , in general , not negligible .
the simplest quantum dielectric function ( df ) which takes collisions into account in a conserving fashion via a relaxation time approximation is the mermin df @xcite @xmath20}{\omega+i\nu[\epsilon_{\rm rpa}(k,\omega+i\nu)-1]/[\epsilon_{\rm rpa}(k,0)-1]},\end{gathered}\ ] ] which involves the finite temperature rpa df and the electron collision frequency @xmath21 . in ref .
@xcite various choices for the collision frequency were studied , so here we will just use a typical set of values to highlight the main trends .
( top ) and @xmath22 ( bottom ) and @xmath23 .
the zoom into the small potential values confirms the correct oscillation frequency at large distances that is expected for friedel oscillations , as seen from the comparison to the ( red ) dashed curve .
here @xmath24 corresponds to the distance indicated by a vertical dashed line .
the grey box shows that the oscillation period of the full potential agrees with the known asymptotics.,title="fig : " ] +   ( top ) and @xmath22 ( bottom ) and @xmath23 .
the zoom into the small potential values confirms the correct oscillation frequency at large distances that is expected for friedel oscillations , as seen from the comparison to the ( red ) dashed curve . here
@xmath24 corresponds to the distance indicated by a vertical dashed line .
the grey box shows that the oscillation period of the full potential agrees with the known asymptotics.,title="fig : " ]
the shukla - eliasson potential has been derived from linearized qhd , below we reproduce its most recent variant @xcite .
it involves the electron plasma frequency @xmath25 , [ @xmath26 is the effective electron mass ] and the electron fermi speed @xmath27 , where @xmath28 is proportional to the wigner - seitz radius .
the term @xmath29^{1/2}$ ] can be regarded as ( optional ) correction that is supposed to account for exchange - correlation effects [ it was proposed in ref .
@xcite based on the exchange - correlation potential of density functional theory ] . since the other potentials do not include such a correction and it does not qualitatively alter the potential , in the comparisons below we will use @xmath30 .    the se potential , for @xmath31 ,
is given by @xcite : @xmath32 e^{-k'_+ r } , \label{eq : sep}\end{aligned}\ ] ] with the parameter @xmath33 .
the wave vectors are given by  @xcite @xmath34 , where @xmath35 is similar ( but not identical ) to the inverse thomas - fermi screening length @xmath19 , eq .
( [ ksc ] ) .
the se potential for the second parameter range , @xmath36 , is given by : @xmath37.\end{aligned}\ ] ] with the definitions @xmath38 and @xmath39 .
note that the shukla - eliasson potential is a zero - temperature approximation , and the potential depends only on density .
numerical results are included in fig .
[ fig : all_vs_rpa0 ] .
here we reproduce the electrostatic potential derived by akbari - moghanjoughi ( we use the most recent corrected version , ref .
this potential has the same mathematical form as the se potential @xmath40 , for @xmath41 : @xmath42 \ , e^{-k'_+ r}. \label{eq : amp>}\end{aligned}\ ] ] the difference lies in the definition of the parameter @xmath43 and of the wave numbers  @xcite @xmath44 where @xmath45 is the inverse thomas - fermi screening length .
these formulas are connected to the se versions by @xmath46 and an additional factor of @xmath1 in the definition of @xmath47 .
the potential for @xmath48 is not given in ref .
@xcite . assuming the same relation to the se
result it has the form @xmath49 , \label{eq : amp<}\end{aligned}\ ] ] with @xmath50 and @xmath51 .
as the se potential , this is a zero - temperature result which , therefore , depends only on the density . numerical results are included in fig .
[ fig : all_vs_rpa0 ] .      here
we reproduce the screened potential of stanton and murillo following the latest corrected version @xcite . again
, there are two cases . for @xmath52 the potential is given by @xcite @xmath53 \ , e^{-k'_+ r},\end{aligned}\ ] ] with @xmath54 , @xmath55 , @xmath56 , and @xmath57 .
@xmath58 denotes the fermi integral and @xmath59 its derivative with respect to @xmath60 .
@xmath61 is determined by the normalization , @xmath62 with the inverse temperature @xmath63 .
the inverse thomas - fermi screening length for finite temperatures is given as @xmath64 .    a second version of the potential follows for the case @xmath65 : @xmath66,\end{aligned}\ ] ] where @xmath67 and @xmath68 .
we note that the sm expressions are identical to the am version in the limit @xmath9 , if @xmath69 is rescaled by a factor of @xmath70 , @xmath71 .
numerical results are included in figs .
[ fig : all_vs_rpa0 ] and [ fig : ms_rpa_t ] .    ]
potentials to the zero - temperature full rpa potential ( full black line ) for three densities .
top : @xmath72 , center : @xmath73 , bottom : @xmath74 .
the sm potential shows increasing deviations from the rpa potential @xmath0 for decreasing density , cf .
top and middle figures . at high density
( bottom figure ) the sm potential exhibits very good agreement with the rpa potential . the se potential exhibits large deviations from the rpa at all densities .
insets show the behavior for large distances and confirm that all models ( except for the rpa ) are unable to describe the friedel oscillations correctly.,title="fig : " ] + ] potentials to the zero - temperature full rpa potential ( full black line ) for three densities .
top : @xmath72 , center : @xmath73 , bottom : @xmath74 .
the sm potential shows increasing deviations from the rpa potential @xmath0 for decreasing density , cf .
top and middle figures . at high density ( bottom figure ) the sm potential exhibits very good agreement with the rpa potential .
the se potential exhibits large deviations from the rpa at all densities .
insets show the behavior for large distances and confirm that all models ( except for the rpa ) are unable to describe the friedel oscillations correctly.,title="fig : " ] + ] potentials to the zero - temperature full rpa potential ( full black line ) for three densities .
top : @xmath72 , center : @xmath73 , bottom : @xmath74 .
the sm potential shows increasing deviations from the rpa potential @xmath0 for decreasing density , cf .
top and middle figures . at high density
( bottom figure ) the sm potential exhibits very good agreement with the rpa potential .
the se potential exhibits large deviations from the rpa at all densities .
insets show the behavior for large distances and confirm that all models ( except for the rpa ) are unable to describe the friedel oscillations correctly.,title="fig : " ]      we now present a numerical comparison of the three lqhd potentials to the rpa potential ( [ pot_stat ] ) . in the following ,
all results are given in hartree atomic units .
since the am potential conincides with the sm potential taken at @xmath9 we will , in the following , not distinguish both potentials .
we use the shortcut `` sm '' for both cases .      to begin the analysis , we recall that the se , sm [ and am ] ion potentials are all derived by applying the long - wavelength limit to the full potential . correspondingly , in fourier space one should expect accurate results only for small wave numbers , @xmath75 , where @xmath76 is the fermi wavenumber .    for a comprehensive comparison we consider a broad density range , covering 6 orders of magnitude with the three values of the brueckner parameter , @xmath77 ( similar values were studied in ref .
@xcite ) . the first observation , cf .
[ fig : all_vs_rpa0 ] , is that all four potentials show the same overall behavior with the deviations growing when the density is reduced .
the next observation is that the sm potential conincides with the rpa potential at large distances : for @xmath78 , for the @xmath79 and @xmath22 , respectively .
however , a closer inspection of the large distance behavior shows that substantial deviations remain : none of the three approximate potentials reproduces the friedel oscillations of the rpa potential , see also figs .
[ fig : rpa_friedel ] and [ fig : rpa_friedel2 ] .
nevertheless , the sm potential provides a good fit to the rpa potential if one discards the ( small ) fourier component with @xmath80 .
consider now the case @xmath73 , cf .
middle part of fig .
[ fig : all_vs_rpa0 ] .
here one clearly recognizes the difference between the sm ( and am ) potentials , compared to the se potential .
while the se potential exhibits strong deviations from the rpa potential , the sm potential is very close to the rpa result . at the highest density , @xmath74 , bottom figure , the deviations between sm and se are very small . but this is not surprising as , in this case , the bohm term in the qhd equations is negligible , so the influence of the different prefactors of this term ( @xmath81 , for se , versus @xmath1 , for sm ) is insignificant . finally , at the lowest density ,
@xmath72 , top figure , all three potentials show substantial deviations from the rpa for small distances below @xmath82 , although the deviations of the se potential are significantly larger , in particular for larger distances
. however , at these low densities , the plasma is strongly correlated and the rpa potential itself is not applicable .
thus , from the zero - temperature behavior we conclude that the sm ( and am ) potential shows a substantial improvement over the se potential , confirming the correctnesss and importance of the prefactor @xmath1 in front of the bohm term @xcite .
in fact , we will show in sec .
[ s : dis ] that this coefficient is not a free parameter but , for the description of the long - wavelength properties of the plasma ( such as the statically screened potential ) the value @xmath1 follows rigorously .
quantum hydrodynamics , as used in the plasma physics community , is a zero - temperature theory , so the results of se @xcite for the effective ion potential and the recent correction by akbari - moghanjoughi @xcite are restricted to the ground state . at the same time , for relevant applications to dense plasmas , warm dense matter or laser plasmas , finite temperature effects are usually non - negligible . it was shown in refs .
@xcite that the qhd equations are directly linked to the thomas - fermi theory which allows for a straightforward incorporation of finite - temperature effects .
this link was exploited in ref .
@xcite to derive a statically screened potential for arbitrary temperatures .
it is now interesting to compare the corresponding finite - temperature sm potential to the rpa potential at the same temperature .
this is done in fig .
[ fig : ms_rpa_t ] for different values of the degeneracy parameter @xmath83 at constant density .
the figure shows that , for high density , @xmath84 , the agreement is excellent , for all temperatures .
this is not surprising since , at these densities gradient corrections ( the bohm term of qhd ) are small , as noted above . here
the temperature dependence is governed by the one of the ideal fermi gas ( fermi pressure ) .
more interesting is the comparison at lower density , @xmath73 , bottom part of fig .
[ fig : ms_rpa_t ] where correlation effects ( and the bohm contribution ) are significant . here
the agreement between sm and rpa potentials is slightly worse and improves with increasing temperature .
this is due to smoothening of the friedel oscillations with increased temperature .
overall , we conclude that the agreement between the sm potential and the rpa , even at finite temperature , is very impressive .
this confirms again that the sm potential correctly captures the long - wavelength properties of the static rpa potential ( [ pot_stat ] ) .    ,
bottom : @xmath73 . as it is seen
the sm potential accurately describes the screening effect at high temperatures and densities.,title="fig : " ] , bottom : @xmath73 .
as it is seen the sm potential accurately describes the screening effect at high temperatures and densities.,title="fig : " ]     and two collision frequencies ( in units of the plasma frequency ) .
the behavior is representative for other densities and temperatures and indicates that the electron collisions do not affect the static screened potential , within the relaxation time approximation ( mermin model ) . ]    finally , let us analyze the effect of collisions ( correlations ) on the dielectric function .
it is well known that these effects become increasingly important outside the weak - coupling limit , i.e. for @xmath85 .
yet it is not clear , a priori , how important these effects are for the static ion potential .
we , therefore , computed the static potential ( [ pot_stat ] ) with the mermin dielectric function for a broad range of temperatures , densities and collision frequencies .
a typical result is shown in fig .
[ fig : ms_mermin_t ] where the rpa potential is compared to the mermin dielectric function with the collision frequency @xmath86 , at zero temperature .
the differences are extremely small , and the same behavior is observed at other densities and elevated temperatures .
the explanation is , of course , that correlation effects are particularly relevant at short distances where the coulomb potential is large anyway .
we note that the situation is very different in streaming plasmas that are out of equilibrium . here
scattering effects have an important impact on the effective ion potential @xcite .
we have presented an analysis of three recently proposed static ion potentials in quantum plasmas and compared them to the static limit of the rpa potential ( [ pot_stat ] ) .
the first of these potentials ( se ) was the one derived by shukla and eliasson @xcite who predicted an attractive minimum of the potential [ cf .
middle part of fig .
[ fig : all_vs_rpa0 ] ] that would be responsible for ordering of ions in dense plasmas .
our analysis reveals that no such minimum exists in the rpa potential , so this minimum should be regarded an artefact of the used lqhd model @xcite .    the comparison with the two other potentials , the ones of stanton and murillo ( sm ) @xcite and of akbari - moghanjoughi ( am ) @xcite which do not exhibit a comparable minimum reveals the origin of this discrepancy ( see below ) .
the sm ( and am ) potential exhibits very good agreement with the rpa potential at @xmath9 indicating that it correctly captures the long - wavelength properties of the rpa , in contrast to the se potential .
since the only difference between the se potential ( [ eq : sep ] ) and the sm ( and am ) potential ( [ eq : sm > ] ) is the different prefactor of the bohm term in the underlying qhd equations , the good accuracy of the latter indicates that the choice of the prefactor @xmath1 is the correct one .
in fact , this question is easily answered by recalling the one - to - one correspondence between qhd and thomas - fermi theory with gradient corrections @xcite . in thomas - fermi theory and ,
more generally , in density functional theory , it is well known that an expansion of the ( non - interacting ) kinetic energy ( together with exchange contributions ) in a power series in terms of density gradients converges to the full hartree - fock energy @xcite . in the present case , we start from the statically screened potential ( [ pot_stat ] ) in rpa , so we should expect agreement from a gradient expansion of the thomas - fermi kinetic energy @xmath87 $ ] without exchange and correlation corrections . since the bohm term of qhd corresponds just to the first non - vanishing gradient correction , @xmath88 , ( quadratic in density gradients ,
see appendix ) , the pre - factors of both are linked inseparably via the equation @xcite @xmath89}{\delta n } = \gamma\ , \frac{\hbar^2}{8m}\left ( \left| \frac{\nabla n}{n } \right|^2 - 2 \frac{\nabla^2 n}{n } \right )    = \gamma \,v^{\rm mh}_b[n ] , \label{eq : grad = vb1}\ ] ] where @xmath90 is the bohm potential for quantum systems with bose statistics . for the case of quantum plasmas ( fermions ) the same potential ( i.e. @xmath91 )
was postulated by manfredi and haas @xcite .
however , the correct value ( at @xmath9 ) is @xmath92 which is clearly confirmed by the comparison of the se and sm potentials to the rpa in the present work .
let us briefly discuss the origin and validity of the prefactor @xmath1 .
this , coefficient was first obtained by kompaneets @xcite and kirzhnitz @xcite , almost 60 years ago , who used a matsubara green functions method to obtain the gradient correction to the thomas - fermi energy functional . in the same year
golden @xcite derived this factor by expanding the density matrix of a many - electron system .
later hohenberg and kohn @xcite presented a systematic derivation in terms of a gradient expansion to all orders , on the basis of the electron polarization function for zero temperature .
mermin @xcite extended the hohenberg and kohn method to finite temperature but did not give a derivation of the gradient correction at finite temperature .
this was done by krizhnitz _
_ @xcite and by perrot @xcite who minimized the grand canonical potential using the known form of the free energy of the electrons . at last ,
the feynman path - integral method was used by yang @xcite to obtain the single - particle green function and to prove the factor @xmath1 . over the recent decades
occasionally additional derivations have appeared . in particular , the recent analysis of akbari - moghanjoughi @xcite essentially reproduces earlier derivations based on the calculation of the long - wavelength limit of the inverse of the lindhard dielectric function , @xmath93 at @xmath9 , whereas the analysis of stanton and murillo @xcite directly starts from perrot s result @xcite .
we also note that there exists a different approach to the gradient corrections in thomas - fermi theory where the prefactor of the respective term is used as a fit parameter to reproduce improved simulations ( that go beyond the rpa ) , e.g. @xcite .
similarly , for the case of electrons bound in atoms , the prefactor can be optimized to better match the hartree - fock ground state energy , yielding @xmath94 @xcite .
however , in the context of the screened potential of free charged particles in a plasma one has to use the static long - wavelength limit of the inverse dielectric function which is exactly given by the rpa , so there exists no freedom of choice , see sec .
[ ss : p - expansion ] .
a key input for the gradient correction of the kinetic energy is the wavenumber expansion of the inverse of the static density response function , i.e. of the longitudinal polarization function @xmath95 .
the result is ( for the derivation , see the appendix ) @xmath96 this is a completely general result where the whole microphysics is contained in the value of the expansion coefficients . from this
we can directly obtain systematic approximations for the statically screened ion potential ( [ pot_stat ] ) without recourse to the connection between thomas - fermi theory and qhd , eq .  ( [ eq : grad = vb1 ] ) . indeed , what is needed to evaluate ( [ pot_stat ] ) is the inverse of the dielectric function which is related to the polarization function by @xcite @xmath97 from this and eq .
( [ eq : pinv_exp ] ) we obtain the desired wavenumber expansion of the inverse static dielectric function [ and of the screened potential ( [ pot_stat ] ) ] @xmath98    using , as the lowest order approximation , the rpa polarization , we obtain for the first five non - vanishing coeffcients ( see appendix ) @xmath99 where the @xmath100 involve fermi integrals of different orders [ see eq .
( [ eq : bi ] ) of the appendix ) ] .
note that the finite temperature results for @xmath101 and @xmath102 are given here for the first time .
this result allows one to systematically derive approximations for the screened potential by considering the long - wavelength limit , @xmath103 of the inverse dielectric function ( [ eq : eps_inv_exp ] ) .
* lowest order result : * the lowest order is given by neglecting all @xmath104-dependent terms in the parantheses of eq .
( [ eq : eps_inv_exp ] ) , @xmath105 , and we obtain @xmath106 inserting this into eq .
( [ pot_stat ] ) , immediately yields the yukawa potential ( [ yukawa ] ) .
* second order result : * the next order is obtained by retaining the terms quadratic in @xmath104 , @xmath107 the result for @xmath108 is ( see appendix ) @xmath109 in the ground state , @xmath9 , this coefficient becomes @xmath110 .
equation ( [ eq : eps_2 ] ) is the approximation used by stanton and murillo @xcite and ( for @xmath9 ) by akbari
moghanjoughi ( am ) @xcite .
the approximation of shukla and eliasson is obtained for @xmath9 by using @xmath111 without the prefactor @xmath112 .
this coefficient yields the ratio of the corresponding results for fermions and bosons ( or spinless particles ) @xcite . from the temperature dependence of @xmath108
we can extract the temperature dependence of the coefficient @xmath113 which is displayed in fig .
[ fig : alpha_3 ] .
@xmath113 increases monotonically from @xmath1 to @xmath114 @xcite , for temperatures large compared to @xmath115 .
evidently this factor is crucial for the correct treatment of fermions within thomas - fermi theory or in quantum hydrodynamics .
in front of the bohm potential @xmath116 , cf .
( [ eq : grad = vb1 ] ) , for fermions as a function of temperature . in the zero - temperature limit
, @xmath113 approaches @xmath1 , and it increases monotonically to @xmath114 @xcite .
in contrast , for bosons or spinless particles this factor equals @xmath81 ( dashed line ) , as in the case of a single quantum particle @xcite . ]    to compare the first and second approximations to the static inverse dielectric function , eqs .
( [ eq : eps_0 ] ) and ( [ eq : eps_2 ] ) , we plot again the associated static potentials computed according to eq .
( [ pot_stat ] ) . in the top two parts of fig .
[ fig : all_vs_rpa0_tf ] we show again the zero temperature results of se and sm and , in addition the thomas - fermi screended yukawa result , eq .
( [ yukawa ] ) .
thus , we are able to assess the quality of the different orders of the wavenumber expansion of @xmath117 , in comparison to the full rpa result . for @xmath73
we see that the yukawa potential ( zeroth order of the expansion ) provides a good approximation , but the next correction , proportional to @xmath118 , that is contained in eq .
( [ eq : eps_2 ] ) and , thus , in the sm potential constitutes an improvement .
in contrast , the se potential which contains the @xmath118 correction with a nine times larger prefactor exhibits a poor performance and is even significantly less accurate than the zeroth order correction alone .
the same trend is observed at higher density ( top figure ) where the zeroth and first order are both indistinguishable from the full result whereas the se potential is substantially less accurate .    the bottom part of fig .
[ fig : all_vs_rpa0_tf ] shows the performance of the zeroth and second order potentials [ tf and sm , eqs .
( [ eq : eps_0 ] ) and ( [ eq : eps_2 ] ) , respectively ] at finite temperature .
while at low temperature , @xmath119 , the second order provides a significant impovement , for higher temperature , @xmath120 , the expansion exhibits sign alternating convergence and both approximations are of comparable quality with a slightly better performance of the sm potential .
in other words , for @xmath120 , the standard static yukawa potential provides a fairly accurate description of the static limit of the rpa .     but with the yukawa potential , eq .
( [ yukawa ] ) , `` tf '' , added for comparison of the ground state behavior .
bottom : temperature dependence of the potentials .
same as bottom part of fig .
[ fig : ms_rpa_t ] with the yukawa potential ( [ yukawa ] ) added.,title="fig : " ] +   but with the yukawa potential , eq .
( [ yukawa ] ) , `` tf '' , added for comparison of the ground state behavior .
bottom : temperature dependence of the potentials .
same as bottom part of fig .
[ fig : ms_rpa_t ] with the yukawa potential ( [ yukawa ] ) added.,title="fig : " ] +   but with the yukawa potential , eq .
( [ yukawa ] ) , `` tf '' , added for comparison of the ground state behavior .
bottom : temperature dependence of the potentials .
same as bottom part of fig .
[ fig : ms_rpa_t ] with the yukawa potential ( [ yukawa ] ) added.,title="fig : " ]    * higher order results .
correlations : * it poses no principal problem to extend the above results for the inverse dielectric function to higher powers of @xmath104 . with the equations ( [ eq : a4][eq : a8 ] )
the next three higher orders are available , for arbitrary temperature .
this should improve the behavior of the static potential also at shorter distances .
however , for the short - distance behavior also other effects will be important , most noticeably exchanges and correlations .
one way to include correlations is to use static local field corrections , as proposed in ref .
a more systematic and powerful approach is to use many - body theory which provides dynamic correlation corrections , e.g. refs .
@xcite and @xcite and the cited references .
the conclusion from the above comparison of approximate static ion potentials to the rpa result is that , in fact , there is no freedom in the choice of the prefactor of the bohm term for quantum hydrodynamics for fermions .
the qhd equations follow from the thomas - fermi theory by performing a systematic gradient expansion of the kinetic energy or , equivalently , a long - wavelength expansion of the inverse polarization function @xmath121 , eq .
( [ eq : pinv_exp ] ) .
this expansion directly leads to the coefficient @xmath1 , at @xmath9 , in front of the bohm term compared to the cases of bosons or of a single - particle , whereas for higher temperatures the coefficient increases monotonically to @xmath122 ( in 1d and 2d different coefficients apply @xcite ) .
this coefficient is a consequence of the pauli principle and has important implications for a hydrodynamic modeling of dense quantum plasmas .
for example , we have shown that a neglect of this coefficient gives rise to a statically screened ion potential that is far less accurate than the standard thomas - fermi ( yukawa ) potential , rather than an improvement .
it also indicates that earlier derivations of qhd equations are invalid for quantum plasmas and earlier qhd results that were using the prefactor @xmath81 do not describe fermions , in general , nor dense plasmas , in particular .
this work has been supported by the deutsche forschungsgemeinschaft via grant sfb - tr 24 , project a9 and the ministry of education and science of kazakhstan .
here we present a systematic derivation of the gradient correction to the thomas fermi energy using mermin s approach for finite temperatures .
this appears to be the most systematic and general approach and allows one to derive higher order corrections as well .
following hohenberg and kohn @xcite , and mermin @xcite , the electron energy is a functional of the density : @xmath123=\int \ !
v({\mathbf{r}})n({\mathbf{r}})\ , \mathrm{d}{\mathbf{r}}+\frac{1}{2}\int \ !
\frac{n({\mathbf{r}})n({\mathbf{r}}^\prime)}{\mid { \mathbf{r}}-{\mathbf{r}}^\prime\mid } \ , \mathrm{d}{\mathbf{r}}\mathrm{d}{\mathbf{r}}^\prime+ t[n],\ ]
] where the first term is the energy due to external field and the second is the interaction energy of the electrons in mean field approximation .
$ ] is the kinetic energy ( we neglect additional exchange - correlation terms ) which is well known in the local approximation ( thomas - fermi theory ) , @xmath124    = \int d{\bf r } \ , t[n({\bf r})]\ , n({\bf r}),\ ] ] where the kinetic energy density is related to the ( local ) fermi energy , @xmath125=\frac{3}{5}e_f[n]$ ] . to obtain corrections to this local approximation
, we expand the kinetic energy in terms of the density perturbation , @xmath126 , @xmath127=t_0[n_0]+\int \ !
k({\mathbf{r}}-{\mathbf{r}}^\prime)\tilde{n}({\mathbf{r}})\tilde{n}({\mathbf{r^{\prime}}})\ , \mathrm{d}{\mathbf{r}}\mathrm{d}{\mathbf{r}}^\prime+\\ + \int \ !
l({\mathbf{r}},{\mathbf{r}}^\prime,{\mathbf{r}}^{\prime \prime})\tilde{n}({\mathbf{r}})\tilde{n}({\mathbf{r^{\prime}}})\tilde{n}({\mathbf{r^{\prime \prime}}})\ , \mathrm{d}{\mathbf{r}}\mathrm{d}{\mathbf{r}}^\prime\mathrm{d}{\mathbf{r}}^{\prime \prime}+ .... \end{gathered}\ ] ] we now transform the second term to an expansion in powers of the density gradient ( gradient corrections ) , by fourier expanding the kernel ( @xmath128 is the volume ) , @xmath129 from eqs .
( [ gn ] ) and ( [ k ] ) we have , using the convolution theorem , @xmath130=t_0[n_0]+\frac{1}{\omega } \sum_{{\mathbf{k } } } { \tilde k}({\mathbf{k}})\tilde{n}({\mathbf{k}})^2 + \dots.\ ] ] we are now looking for the long - wavelength limit and expand @xmath131 since @xmath132 is real , eq .
( [ k ] ) leads to @xmath133 , therefore , the expansion ( [ k3 ] ) contains only even powers of @xmath134 . from this , we obtain for @xmath135 , using ( [ k ] ) , @xmath136 substituting ( [ k4 ] ) into ( [ gn ] ) one obtains @xmath137 & = t_0[n_0 ] + \sum_{l=1 }   t_{2l}[n ] , \\
t_{2l}[n ] & =   \int \mathrm{d}{\mathbf{r}}\,a_{2l}\left[n\right ] \mid \nabla^l n({\mathbf{r}})\mid^{2 } ,   \label{term_2l } \ ] ] where we took into account that @xmath138 and identifiy @xmath139 .
thus we have rewritten the kinetic energy in terms of a local term plus gradient corrections .
it was shown in ref .
@xcite that this energy functional leads directly , via functional differentiation with respect to the density profile , to quantum hydrodynamic equations with the local term yielding the fermi energy and the first gradient term giving rise to a bohm - type potential , @xmath89}{\delta n } = \gamma\ , \frac{\hbar^2}{8m}\left ( \left| \frac{\nabla n}{n } \right|^2 - 2 \frac{\nabla^2 n}{n } \right )    = \gamma \,v^{\rm mh}_b[n ] , \label{eq : grad = vb}\ ] ] where @xmath90 is the potential postulated by manfredi and haas @xcite .
what remains is to establish the relation between the coefficient @xmath113 and the coeffient @xmath140 of the expansion ( [ term_2l ] ) . in the following we prove that , for zero temperature , to lowest order in the correlation energy , @xmath141 which means that @xmath142 .      to evaluate the coefficient @xmath139 we recall that the energy associated with a charged particle density fluctuation can be expressed , in linear response , in terms of an effective potential ( external plus induced potential ) @xmath143=t_0[n_0]+\frac{1}{2\omega } \sum_{{\mathbf{k}}}\tilde{n}({\mathbf{k } } ) \tilde{u}_{\rm eff}({\mathbf{k}}).\ ] ] recalling the definition of the longitudinal polarization function , @xmath144 , eqs .
( [ gn2 ] ) and ( [ gn3 ] ) allow us to identify @xmath145 and to make use of known results for @xmath146 below .
the lowest order many - body approximation for @xmath146 is the random phase approximation ( rpa ) which reads , for an arbitrary temperature @xcite @xmath147,\ ] ] where @xmath148 , @xmath149 , @xmath150 , @xmath151 , @xmath152 , and @xmath153 where @xmath154 , and @xmath155 is the chemical potential . to obtain the long - wavelength limit of @xmath121 , eq .
( [ k3 ] ) , we now expand the inverse of the static ( i.e. @xmath156 ) polarization function in the limit @xmath157 :    @xmath158 h_1(\eta ) & = \frac{\sqrt{\theta}}{2 } \ , i_{-1/2}(\eta ) , \quad h_2(\eta)=\frac{1}{2\sqrt{\theta } } \,i_{-3/2}(\eta).\end{aligned}\ ] ]    the first term in ( [ k5 ] ) yields the coefficient @xmath159 in eq .
( [ k3 ] ) which does not contribute to the kinetic energy since the integral of @xmath160 vanishes whereas the second term in ( [ k5 ] ) yields the coefficient @xmath140 : @xmath161,\theta\right ) & = -\frac{\pi e^2}{2k_f^2\chi_0 ^ 2 h_1(\eta ) } , \\
a_2\left(\left[n\right],\theta\right ) & = -\frac{\hbar^2 h_2(\eta)}{72m_e n h_1 ^ 2(\eta)}.\end{aligned}\ ] ] in the ground state , @xmath162 , we have the limits @xmath163 and @xmath164 , so @xmath165 , and ( [ a2 ] ) gives @xmath166 meaning that the coefficient in front of the bohm potential in eq .
( [ eq : grad = vb ] ) becomes @xmath142 .
the temperature dependence of this coefficient was displayed in fig .
[ fig : alpha_3 ] .    finally , substituting @xmath140 from ( [ a2 ] ) to the formula ( [ gn4 ] ) , for the kinetic energy at finite temperature we have @xcite : @xmath167=t_0[n_0]+ \frac{1}{9 } \frac{\hbar^2}{8m_e } 4\beta^{3/2 } \int \ !
e_{f}^{3/2}(n)\frac{i_{-1/2}^{\prime}(\eta)}{i_{-1/2}^2(\eta ) } \frac{\mid \nabla n\mid^{2}}{n } \ , \mathrm{d}{\mathbf{r}},\ ] ] were we eliminated the fermi integral of order @xmath168 using the formula @xmath169 .
now consider the higher order ( @xmath171 ) terms of the gradient expansion ( [ gn4 ] ) . in the long - wavelength limit ,
@xmath172 the kernel @xmath121 can be written as : @xmath173 where @xmath174,\theta\right)=\frac{\theta^{-i}}{2i+1 } \frac{i_{-i-1/2}(\eta)}{i_{-1/2}(\eta)}. \ ] ] expanding ( [ k6 ] ) it is straightforward to obtain the higher ( even ) order terms of the gradient correction : @xmath175=\frac{\hbar^2 \pi^2}{16 m_e }   \int \ ! \frac{b_1 ^
2-b_2}{\theta\left[n\right]^{1/2}k_f^5i_{-1/2}(\eta)}\mid \nabla^2 n\mid^{2 } \ , \mathrm{d}{\mathbf{r } } , \nonumber\ ] ] @xmath176=\frac{\hbar^2 \pi^2}{64 m_e }   \int \ !
\frac{b_1 ^ 3 - 2b_1b_2+b_3}{\theta\left[n\right]^{1/2}k_f^7i_{-1/2}(\eta)}\mid \nabla^3 n\mid^{2 } \ ,
\mathrm{d}{\mathbf{r } } , \nonumber\ ] ] @xmath177=\frac{\hbar^2 \pi^2}{256 m_e}\times \\   \int \ !
\frac{b_1 ^ 4 - 3b_1 ^ 2b_2+b_2 ^ 2 + 2b_1b_3-b_4}{\theta\left[n\right]^{1/2}k_f^9i_{-1/2}(\eta)}\mid \nabla^4 n\mid^{2 } \ , \mathrm{d}{\mathbf{r}}. \nonumber   \end{gathered}\ ] ] in the zero - temperature limit , @xmath178 and @xmath179 were obtained by hodges @xcite and murphy @xcite , respectively .
for finite temperature these terms and , additionally , @xmath180 are given here for the first time .
here we do not discuss the convergence of this gradient series .
obviously , before using higher order terms one should verify that all terms of the same order arising from three - particle correlations [ third term in ( [ gn ] ) ] are included as well .
further , it is well possible that there exist comparable terms arising from correlation effects beyond the rpa , e.g. @xcite .
finally we note that these results apply only to three - dimensional plasmas .
results for 1d and 2d systems were discussed in ref .
extensions of these results to finite temperatures are given in ref .
10 t. kluge , c. gutt , l. huang , j. metzkes , u. schramm , m. bussmann , and t.e cowan , _ using xfels for probing of complex interaction dynamics of ultra - intense lasers with solid matter_. phys .
plasmas * 21 * , 033110 ( 2014 )
.        m.r .
_ , _ experimental demonstration of fusion - relevant conditions in magnetized liner inertial fusion _ , phys . rev . lett . * 113 * , 155003 ( 2014 ) ; p.f .
et al . _ ,
. lett . * 113 * , 155004 ( 2014 ) .
h. khlert , g. kalman , and m. bonitz , _ dynamics of strongly correlated confined plasmas _ ,
e * 90 * , 011101(r ) ( 2014 ) ; _ linear fluid theory for weakly inhomogeneous plasmas with strong correlations _
plasma phys .
* 55 * , 352 ( 2015 ) .
p.  k. shukla , and b. eliasson , _ novel attractive force between ions in quantum plasmas _ ,
108 * , 165007 ( 2012 ) ; _ erratum _ , _ ibid _ : * 108 * , 219902 ( 2012 ) ; _ erratum _ , _ ibid _ : * 109 * , 019901 ( 2012 ) .
m. bonitz , e. pehlke , and t. schoof,_comment on  discussion on novel attractive force between ions in quantum plasmas failure of simulations based on a density functional approach " _ , phys .
scr . * 88 * , 057001 ( 2013 )    govind s. krishnaswami , rajaram nityananda , abhijit sen , anantanarayanan thyagaraja , _ a critique of recent semi - classical spin - half quantum plasma theories _ , contrib .
plasma phys .
( 2014 ) , doi : 10.1002/ctpp.201400043 , arxiv:1407.6865 sergey v. vladimirov , _ plasma and potential distributions of moving objects in classical and quantum plasmas _ , _ plasma phys .
control . fusion _ * 53 * , 074005 ( 2011 ) .
n. crouseilles , p .- a .
hervieux , g. manfredi , _ quantum hydrodynamic models for nonlinear electron dynamics in thin metal films _ ,
b * 78 * , 155412 ( 2008 ) . in the original formulas @xmath181
are denoted by @xmath182 and @xmath183 , respectively .
a. moldabekov , p. ludwig , j.p .
joost , m. bonitz , and t. ramazanov , _ dynamical screening and wake effects in classical , quantum , and ultrarelativistic plasmas , _ , contrib .
plasma phys .
* 55 * , 186 ( 2015 ) .
a. s. kompaneets , e. s. pavlovskii , zh .
31 * , 427 ( 1956 ) .
d. a. kirzhnitz , _ quantum corrections to the thomas - fermi equation _ , zh . eksp .
. fiz . * 32 * , 115 ( 1957)(sov
jetp , * 5 , 64 * ) .
sidney golden , _ statistical theory of many - electron systems .
general considerations pertaining to the thomas - fermi theory _ , phys .
rev . * 105 * , 604 ( 1957 ) .
p. hohenberg , w. kohn , _ inhomogeneous electron gas _ , phys .
rev . * 136 * , b864 ( 1964 ) .
n. david mermin , _ thermal properties of the inhomogeneous electron gas _ , phys . rev . * 137 * , a1441 ( 1965 ) .
kirzhnitz , yu.e .
lozovik , and g. v. shpatakovskaya , _ statistical model of matter _ , usp .
nauk * 117 * , 3 - 47 ( 1975 ) [ sov .
, * 18 * , no . 9 , ( 1976 ) ] f. perrot , _ gradient correction to the statistical electronic free energy at nonzero temperatures : application to equation - of - state calculations _ , phys . rev .
a , * 20 * , 586 ( 1979 ) .
weito yang , _ gradient correction in thomas - fermi theory _
a , * 34 * , 4575 ( 1986 ) . g. manfredi , and f. haas , _ self - consistent fluid model for a quantum electron gas _ , phys .
b * 64 * , 075316 ( 2001 ) .
nestor r. arista and werner brandt , _ dielectric response of quantum plasmas in thermal equilibrium _ ,
a. , * 29 * , 1471 ( 1984 ) .
h.  xu , and j .-
hansen , _ plasma - insulator transition of spin - polarized hydrogen _ ,
e * 60 * , r9(r ) ( 1999 ) .
k. yonei , and y. tomishima , j. phys .
japan * 20 * , 1051 ( 1965 ) .
gross , and r.m .
dreizler , phys .
a * 20 * , 1789 ( 1979 ) . c. h. hodges , _ quantum corrections to the thomas - fermi approximation  the kirzhnits method _ , can
. j. phys .
* 51 * , 1428 ( 1973 ) .
d. r. murphy , _ sixth - order term of the gradient expansion of the kinetic - energy density functional _ , phys .
a. , * 24 * , 1682 ( 1984 ) .
kwong , and m. bonitz , _ real - time kadanoff - baym approach to plasma oscillations in a correlated electron gas _ , _ phys .
lett . _ * 84 * , 1768 ( 2000 ) .
a. moldabekov , t.  s. ramazanov , and m. bonitz , to be published .