Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: when a commensurate one - dimensional charge density wave ( cdw ) is doped , the resulting charges may be viewed as defects in the cdw . if the interactions between them are sufficiently strong compared to their effective kinetic energy , then they will form an ordered lattice ( which can be stabilized by a crossover to 3d order if there are many @xmath2 chains coupled together ) of defects and the cdw will become incommensurate with the lattice . this can be more energetically favorable than simply doping holes into the rigid band formed in the presence of cdw order at a fixed wavevector . the reason is that , in the latter case , we must pay the single - particle energy gap of the cdw . in a two - dimensional @xmath0-density wave ( ddw ) state , however , there are nodes in the order parameter , so it is relatively painless for doped holes to simply go to the nodes , which expand into fermi pockets ( in quasi-@xmath2 systems , however , this might not be possible , so doping must lead to incommensurability @xcite ) . however , this can not continue indefinitely if the ddw state is stabilized by approximate nesting , since the fermi surface ( fs ) would eventually move away from the nesting wavevector . thus , incommensurate ddw order is a strong possibility , at least over some range of hole dopings . in this paper , we explore this possibility . near half - filling , the hubbard model , the @xmath3 model , and generalizations of these appear to have many phases which are nearby in energy @xcite . thus , small perturbations can strongly influence the competition between them . consequently , it should not be too surprising that experiments on the cuprate superconductors have also uncovered evidence for a cornucopia of phases , particularly on the underdoped side of the phase diagram , which appear in particular materials and for certain values of the doping level , temperature , magnetic field , etc . @xcite . depending on the interfacial energies between these phases , one way in which their competition can be resolved is through the formation of stripes @xcite or other inhomogeneous patterns of microscale phase - separation . neutron scattering and scanning tunneling microscopy experiments point towards this possibility@xcite . furthermore , as the cuprates are doped , their fermi surfaces evolve away from the commensurate nesting wavevector @xmath4 , so we would expect translational symmetry - breaking order parameters to occur at incommensurate wavevectors . hence , it would seem important for any theory of the pseudogap to incorporate the tendency towards incommenuration , which we take as further impetus to study incommensurate ( ic ) order parameters related to ddw order . in this paper we consider various ways in which @xmath0-density wave ( ddw ) order can go incommensurate . before turning to this discussion , we review the reasons for expecting ddw order to occur in the pseudogap regime of the cuprates . it has been suggested@xcite that competition between ddw order and @xmath5-wave superconductivity ( dsc ) can explain many of the anomalous properties of the pseudogap@xcite . the ddw order parameter is given by @xmath6 with @xmath7 and a _ commensurate _ ( c ) ordering wavevector @xmath8 . the particle - hole condensate in breaks translation by one lattice spacing , parity , and time - reversal but is invariant under the combination of any two of the above symmetries . in real space it is represented by currents which alternate from one plaquette to the next . one motivation for considering ddw order is to explain @xcite the abrupt depletion of the superfluid density below the critical @xmath9 in bi2212 and similar dopings in other materials @xcite . because the order has a @xmath5-wave momentum dependence , its onset would explain arpes measurements above @xmath10 in the pseudogap @xcite . furthermore , since ddw is a spin - singlet condensate , it would also explain the suppression of the spin susceptibility in nmr measurements and opening of a spin gap in inelastic neutron scattering @xcite . also , the depletion of the c - axis conductivity @xcite , whose contribution mostly comes from the antinodal regions , would also be naturally explained by the development of ddw order , which gaps these regions . direct attempts to measure ddw order through neutron scattering have neither ruled it out @xcite nor unambiguously verified its presence , but there are intriguing suggestions that it may be present@xcite . thus , it is natural to ask whether ddw order is compatible with incommensurability . ( for a discussion of the experimental signatures of ddw order , the reader is referred to the literature @xcite . the case of ddw in @xmath11 is discussed in@xcite . for numerical evidence of the existence of ddw in other systems , one is referred to@xcite . ) the plan of our paper is as follows . in section [ susc ] we develop some preliminary intuition about the energetics of ic density wave order by analyzing the ddw susceptibility of a fermi liquid ( both the nested and non - nested cases ) . in section [ c - ic ] we review some of the theory of commensurate ( c ) ddw order and show various ways to make it incommensurate . we also suggest the kind of microscopic physics that could lead to c and ic ddw , and in addition discuss possible experimental signatures of ic ddw . finally , in section [ num ] , we compute the mean - field phase diagrams of an extended hubbard model . our main results are that ( a ) ic ddw develops for a wide range of parameters , but ( b ) the ic wavevectors can be very small . their relevance to the cuprates may , at most , be for a rather narrow set of dopings . we also offer a physical explanation of our results , which sheds some light on the nature of c - ic ddw transition in extended hubbard models . we note that similar incommensurate order has been found in refs . the easiest way to see that ic ddw phases could arise in effective models ( which , in this paper , will , in turn , emerge from microscopic extended hubbard models with correlated hopping terms@xcite ; see section [ micro ] for details ) is to consider the bare susceptibility as a function of the incommensurability wavevector @xmath12 , where @xmath13 : @xmath14 where @xmath15 is defined after eq . [ cop ] . here , for simplicity , ddw is considered without the presence of any other order . in eq . [ suscep ] , @xmath16 is the fermi occupation number , and the energy dispersion contains nearest , @xmath17 , and next - nearest neighbor , @xmath18 , hopping parameters , @xmath19 . the ddw equivalent of the stoner criterion is satisfied when the right - hand - side of this equation approaches the inverse of the pertinent coupling @xmath20 . this observation allows us to interpret the positions of the peaks of @xmath21 as the wavevectors at which ddw order is likely to occur as we lower the temperature . a similar analysis was done by schulz@xcite for the case of af and ic sdw order in the hubbard model close to half - filling , where perfect nesting at half - filling was asumed . in fact ,     is at nonzero @xmath12 . here , @xmath22 is measured along the zone diagonal in units of @xmath23 , where @xmath24 is the lattice spacing.,height=326 ]    his results can be recovered by letting @xmath25 in if we set @xmath26 . the support of @xmath27 in the integrand in is at @xmath28 and symmetry - related antinodal points of the fs since this is where the ddw order parameter is large . however , since these points are also the location of van hove singularities , the regions around them also dominate the af susceptibility . thus , at least in the @xmath26 case , one expects that whenever an ic sdw instability occurs , an ic instability for ddw should also occur . indeed , for both orders , the peak of @xmath21 occurs at @xmath29 for low enough temperatures and sufficiently away from half - filling , e.g. if @xmath30 , as shown in fig . [ suscfig ] . note that the same analysis repeated with any other density wave or spin - density wave order would yield an ic instability as long as the dominant channels of the modulation coincide with fermi surface points which give large contributions to the susceptibility . for instance , a @xmath31-density wave@xcite with wavevector @xmath32 has a peak in @xmath33 for @xmath34 away from half filling at @xmath35 , in the case @xmath26 . a natural question to ask is : to what extent does the ic ddw instability mimic the ic sdw one ? to answer this , one needs to consider specific models as we do in [ micro ] . however , two differences can be anticipated even beforehand . first , the ddw order parameter has nodes , while the af order parameter does not . second , in a more realistic model , one must consider competition with dsc , which also could change features of the transition to ic state . a better understanding of ic ddw order parameter is necessary to better understand the the situation away from half - filling . given that the susceptibility points towards an ic instability , the next issue is what ic ddw order would look like . before addressing this question , let us first summarize some features of the c - case.@xcite in the c case , a general singlet particle - hole condensate can be written as ( neglecting the @xmath36 contribution ) : @xmath37 here , @xmath38 is the magnitude of the order parameter . the wavevector @xmath39 is a reciprocal lattice vector . the spin index @xmath40 included above will be omitted in the future for simplicity unless clarity requires otherwise . in , the form factor @xmath15 , which can transform non - trivially under the point group of the 2-d lattice , is related to the angular momentum of the particle - hole condensate . for ddw order , @xmath41 which corresponds to @xmath42 and @xmath0 symmetry . a @xmath31-density - wave with @xmath43 has @xmath44 , which corresponds to @xmath45 .    in real space , becomes @xmath46 where the vertices @xmath47 of the orders discussed above are given by @xmath48 note that the factors of @xmath49 in @xmath50 and @xmath51 signify that the corresponding phases break time - reversal , e.g. current flows along the bonds . positive / negative signs in represent current going in / out of the vertex . translation invariance , by default , is broken by any density wave order since there is a preferred vector @xmath52 , but in the c case , invariance under a _ subset _ of the lattice group is retained . ddw order , for example , is invariant under translation by linear combinations of even reciprocal lattice vectors along each bond direction . the current along the bonds is given by : @xmath53 here @xmath54 is along either the @xmath55- or @xmath56-direction . because of the equivalence of @xmath52 and @xmath57 in the c case , the current is simply @xmath58 . in the ddw case this results in an alternating plaquette current / checkerboard magnetic field patterns.@xcite      to obtain an ic version of ddw order , one first notes that ddw is chiral and hence breaks the ising symmetry which reverses all the clockwise plaquette currents into counter - clockwise and vice versa . therefore one can construct anti - phase domain walls as in fig.[fig : domains ] . a simple caricature of a bond - oriented domain wall is given in fig.[fig : domains]a . its real space representation is @xmath59 where @xmath60 is the step function , so @xmath61 with @xmath62 . at the midpoint of the domain wall , the currents in the @xmath63-direction vanish while the currents in the @xmath64-direction are halved . the diagonal domain wall depicted in figure [ fig : domains]b has a similar representation . note that because of the halving of the @xmath64-direction currents at the midpoint of the bond - domain wall in , we need to include the term : @xmath65 where the new vertices @xmath66 are given by : @xmath67 such vertices seem necessary for ic transitions of current density waves . one could imagine situations in which the transition to a disordered state is driven by thermal proliferation of such vertices , similar to earlier analysis of the 6-vertex model@xcite . although the analysis of such a transition is an interesting topic on its own , we will not pursue it in this paper . we are more interested in the energetics of the ic phase , so a long - wavelength mean - field treatment is sufficient .    to analyze the c - ic transition , one then considers an array of alternating domain walls with average separation @xmath68 . then , in the single harmonic approximation one retains only the @xmath69 term in and identifies @xmath70 with the incommensurability wavevector . the new short - range vertex contributions are also neglected in the lowest harmonic approximation . the neglect of higher harmonics should not change the transition temperature@xcite or the value of the critical doping for the onset of ic order since @xmath38 is small and the corrections are of higher order in @xmath38 . an incommensurate instability in a clean system occurs when @xmath68 becomes finite . ( in a disordered system , it occurs when when @xmath71 , where @xmath72 is a characteristic length scale for the disorder . )    within the single harmonic approximation , an array of anti - phase domain walls with separation @xmath73 , looks simpler in momentum space : @xmath74 note that in the limit @xmath75 , the above reduces to . one could have guessed the result by writing @xmath76 in the commensurate case , and then let @xmath52 go incommensurate . this suggests other ways to get ic order . in particular the oddness of ddw order with respect to rotations by @xmath77 and transpositions about , say , the @xmath4 direction lead to the following example of a ddw _ checkerboard pattern _ : @xmath78 where the operation @xmath79 can be either a transposition , or a rotaion by @xmath77 . the checkerboard pattern in corresponds to simply superimposing two domain walls rotated by @xmath77 with respect to each other .    in a similar spirit , one could also investigate the possibility of _ non - topological _ domain walls@xcite , e.g. ic domains as the above , superimposed with an uniform order , so the overall order parameter does not change chirality . however , as we will see in sec . [ num ] , much of the physics of the c - ic transition will already be transparent from the simpler choices and of order parameter . finally , note that within the single harmonic approximation , and are current - conserving only to lowest order in @xmath22 . for the purpose of clarifying the energetics of the ic order close to the c - ic transition , or in the cases where the magnitude of the incommensurability is small , this lowest order approximation is good enough . our numerical results , sec . [ num ] , _ a posteriori _ show that we are in such regime .      we will consider the above order parameters in the context of the following hamiltonian @xcite : @xmath80 in this formula , @xmath81 is hopping with @xmath82 for nearest neighbors , @xmath83 for next nearest neighbors and @xmath84 otherwise . @xmath85 is a correlated hopping term which simultaneously hops an electron from site @xmath86 to site @xmath49 and hops an electron from @xmath87 into the vacated site @xmath86 . the on - site and nearest - neighbor repulsions are , respectively , @xmath88 and @xmath89 . the indices @xmath90 signify lattice sites and @xmath91 the spin . we will consider this hamiltonian at the mean - field level . at this level , the energetics of ddw order is the same as for a ` bcs - reduced ' hamiltonian for ddw order of the form : @xmath92 with @xmath93 . a similar reduced hamiltonian for @xmath5-wave superconductivity has @xmath94 . note that the processes responsible for both dsc and ddw are essentially kinetic . a hamiltonian of the form ( [ micro - ham ] ) has one important virtue : it has a @xmath5-wave superconducting ground state over a range of dopings ( and an antiferromagnetic ground state at half - filling ) . since this is the _ sine qua non _ for any description of the cuprates , we believe that this is a good starting point . one may object that correlated hopping terms have to be put in by hand in a microscopic hamiltonian . however , they naturally arise even from the one - band hubbard model away from half - filling , where their contribution can be computed in @xmath95 perturbation theory : @xmath96 at dopings @xmath97 , for @xmath98 , we would have @xmath99 . thus , not far away from half - filling , correlated hopping terms are certainly not negligible . let us generalize the mean field description of c ddw to the ic case . the path we will take is similar to that in sec . [ icop ] , where we wrote down the _ commensurate _ order parameter without making use of the @xmath100 symmetry , and as a result the generalization to incommensurate @xmath52 was straightforward . we thus obtain a functional @xmath101 , equation , whose minimization yields self - consistently both the order parameters @xmath102 and the deviation from commensurability @xmath103 . the only subtlety , as we discuss at the end of this section , is that although @xmath104 reduces to the mean - field free energy in the c phase , in the ic phase it is _ not _ the free energy even within the single harmonic approximation . instead , minimization of @xmath104 is equivalent to a variational search of the ground state . our derivations generalize the treatment of nayak and pivovarov @xcite of competition of c ddw with dsc and possibly af . the c ddw mean field hamiltonian we start with is : @xmath105 where @xmath106 . note that in the @xmath107 symmetry for _ incommensurate _ @xmath52 is not equivalent to the hermiticity of the hamiltonian due to the @xmath108-dependance of @xmath109 ( e.g. @xmath110 for ic * q * ) . this is unlike the case of af order , where the interaction is simply @xmath111 and we have @xmath112 .    in order to check with earlier results on the c - ic transition of af order in the _ weak coupling _ hubbard model@xcite ( see sec . [ susc ] ) , we will also include a term like in our full reduced hamiltonian . of course , in the strong coupling , @xmath113 , limit , one is not justified in treating af as a fermi surface instability , especially close to half filling . a @xmath3 description is more appropriate in this case . however , as mentioned in sec . [ micro ] , both ddw and dsc are essentially kinetic in origin within the extended hubbard model , so even in the strong onsite repulsion limit , there is a good justification to use mean field theory to describe their competition away from half filling .    in the c case , the standard dsc term in our reduced hamiltonian is : @xmath114 where @xmath115 . using a nambu basis @xmath116 , the full hamiltonian @xmath117 can be rewritten as @xmath118 where : @xmath119 with @xmath120 . the action derived from the above hamiltonian yields , upon integrating the fermions out , the free energy in the ordered state : @xmath121 with @xmath122 the energy eigenvalues of , and @xmath123 if we assume @xmath100 symmetry , @xmath124 would be given by : @xmath125 ^ 2}\ ] ] where @xmath126 and @xmath127 . however , in the ic case , @xmath100 need not hold , in which case the above would generalize to :    @xmath128\delta^2_\text{dsc}+ f^4_-({\bf k})\delta^4_\text{dsc}. \end{split}\ ] ]    the effect of a nonzero incommensurability @xmath12 can be seen not only in the appearance of the terms @xmath129 and @xmath130 , but also in the dsc - and doping - dependent function @xmath131 . our final step before letting @xmath52 be incommensurate is to pass the integration in from rbz to bz : @xmath132 once we convert the rbz integrals to bz integrals as in , we numerically minimize ( see sec . [ num ] ) the resulting ` free energy ' with respect to both the order parameters @xmath102 and the deviation from commensurability @xmath12 . ev and @xmath133ev are large ( compare to figs . [ fig : ddwpd3 ] and [ fig : ddwpd4 ] ) @xmath134ev . the next nearest neighbor hopping , @xmath135ev , is small compared to typical values fitted to arpes@xcite , but ic order was nevertheless robust away from half - filling . inset : ic wavevector @xmath22 at @xmath136ev as a function of doping . wavevectors are in units of @xmath137 , where @xmath24 is the lattice spacing.,height=288 ]    it must be emphasized that the expression to be minimized in sec . [ num ] is _ not _ the true free energy of the system if @xmath52 is incommensurate . one way to see this is to note that the eigenvalues are not invariant under @xmath138 if @xmath139 and therefore they can not be interpreted as band energies . in a c - ic transition , the excitations split into a generally uncountably many bands even within the single harmonic approximation@xcite . therefore , the problem of finding the mean - field free energy is notoriously difficult even in the single harmonic limit . the calculation we are doing , on the other hand , is essentially a mean - field variational search for the ground state , and not for the excitations . a similar calculation was done by schulz@xcite for the case of the c - ic transition of the weak - coupling hubbard model close to half - filling . by continuity of the true free energy@xcite , one expects that the resulting variational estimates for @xmath140 would be accurate only close to the c - ic transition , where @xmath22 is small . ev , @xmath141ev are ` large ' ( note temperature scale ) similar to fig . [ fig : ddwpd1 ] , with @xmath142ev . however this time we taka a more realistic @xmath143ev . again ic order is robust at dopings relevant for the pseudogap . inset : ic wavevector @xmath144 at @xmath136ev as a function of doping . [ fig : qdw4 - 2tp012 ] for a more detailed view of @xmath144 ( and also fig . [ fig : qdw4 - 2tp006 ] for the case of @xmath145ev ) .,height=288 ]    ev , @xmath146ev , @xmath147ev . the next nearest neighbor hopping , @xmath148ev . note the appearance of dsc proximate to ` magic ' dopings . inset : the ic wavevector @xmath149 at @xmath150ev . see fig.[fig : qdw25 - 1tp012 ] for comparison between @xmath149 and @xmath151 . , height=288 ]    to obtain the phase diagrams for the models in , for various fixed dopings @xmath63 and temperatures @xmath152 , we minimize @xmath153 with respect to both the order parameters @xmath102 ( as was done previously in the c - case@xcite ) and the ic wavevector magnitude @xmath22 . each choice of ic order parameter ( e.g. bond or diagonal ) corresponds to a separate minimization procedure . to distinguish which ic order dominates at a particular point in @xmath154 space , we compared the free energies of the orders considered , and chose the ic phase of minimal energy . note that , as given in , the free energy is a function of @xmath155 , so to obtain the dependance on the doping @xmath63 , we first perform a legendre transformation and minimize instead @xmath156 where the chemical potential is obtained from @xmath157 with @xmath158 . our approach is slightly different than that previously used in the c - case@xcite . there , the competition between af , ddw and dsc was analyzed , and some sensitivity on the details of the microscopics was found . phase diagrams resembling that of the cuprates were obtained . the general trend was that any phase diagram which included @xmath5-wave superconductivity also had ddw order competing with it . by changing the bare coupling constants , rather different phase diagrams could be obtained and these generally did not have superconductivity .    in our treatment here , we do not consider the competition with af order , which is not a fermi surface instability anyway in the context of the cuprate superconductors . in , we still include af , but only in order to check with previous results@xcite of ic af order away from half filling in the _ weak coupling _ @xmath159 limit . on the other hand , the af order appearing at half filling in the cuprates is a strong coupling phenomenon . a reduced hamiltonian is the wrong starting point at half filling in the limit @xmath113 . therefore , we only consider competition between ddw and dsc and look for ic order in this case . all the phase diagrams constructed this way are robust in the sense that ddw and dsc are found to coexist at mean field level for almost all reasonable values of the coupling constants . sample phase diagrams are shown in figs . [ fig : ddwpd1]-[fig : ddwpd4 ] . in all of them , ic order is present at dopings @xmath160 where ddw and dsc compete . in the first two of these figures , we vary @xmath18 at relatively high @xmath161 and @xmath162 , while in the last two we explore the phase diagram for more realistic ddw and dsc couplings ( compare the temperature scales to those in the cuprates ) . note the complexity of figs . [ fig : ddwpd3 ] , figs . [ fig : ddwpd4 ] compared to figs . [ fig : ddwpd1 ] , figs . [ fig : ddwpd2 ] . at realistic couplings , the ic wavevector profile is complicated , and furthermore dsc order appears in more than one doping region within the ddw dome . in fact , the appearance of incommensurability seems correlated with the appearance of dsc . a more systematic analysis to be discussed below reveals that generally the types of c - ic transitions that were produced in our phase diagrams fall into two categories : ( a ) those which correspond to motion of the bare fermi surface away from commensurate nesting ( shaded regions ) ;     except @xmath163ev . there is a distinct ic response ( inset ) at @xmath164 and @xmath165 . [ fig : qdw2 - 1tp012 ] for comparison between @xmath151 and @xmath149.,height=288 ]     where ddw and dsc compete . note that nonzero @xmath22 is found above as well as below the dsc critical temperature.,height=288 ] \(b ) those which appear near certain ` magic ' filling fractions . case ( a ) is similar to the familiar c - ic mechanism discussed by schulz@xcite , while case ( b ) appears to be a consequence of negative effective stiffness at special fractions , the exact nature of which at present is not clear . in both cases @xmath22 is independent of the onset of dsc , fig . [ fig : qoft ] , but because incommensurability relieves the ddw - dsc competition , the onset of dsc correlates to the presence of ic ddw . let us first discuss ( a ) . it is well known@xcite that in the weak - coupling hubbard model at @xmath26 , the af instability at half - filling is driven ic for any non - zero doping ( see sec . [ suscep ] ) . the source of the c - ic transition can be traced to the approximate nesting of the bare fermi surface ( fs ) at ic wavevectors away from half - filling . close to half - filling , this picture would produce a linear deviation from commensurability , @xmath166 , fig . [ fig : afpd1 ] . the incommensurability ( deviation from @xmath4 ) is horizontal or vertical , not diagonal , as we would expect because the bare susceptibility is more divergent due to a greater overlap between nested portions of the fermi surface@xcite . a similar scenario happens in the case of ddw order . however , two notable differences from the af case are present . first , ddw is less sensitive to @xmath167 , since it is driven mainly by antinodal regions of high density of states ( d.o.s . ) . therefore , while there is no af instability without nesting ( at half - filling af fermi surface instability is only possible if @xmath168 ) , ddw is still possible in this regime . this suggests that whenever the antinodal regions of the bare fs evolve away from the anti - nodal regons ( which happen to also be regions with a high density of states as a result of the presence of van hove singularities ) , a c - ic transition is likely to occur .    . the results of schulz@xcite were recovered . here , we used @xmath169ev with @xmath134ev . bond - oriented ic order was found to be energetically favorable , as a direct computation@xcite in the weak repulsion regime would yield . _ inset _ : as expected , the ic wavevector at @xmath170 is linear in @xmath63 . , height=288 ]    the second difference from the af case is that ddw has nodes . this implies that the chemical potential in the ordered state need not be higher than the maximal ddw gap . as a result , ddw order can stay commensurate for wider ranges of doping . this is indeed what we find . for example , for all @xmath171 , we obtain @xmath172 near zero temperature . in general , we find that the doping values @xmath173 for which the above type of c - ic ddw is possible are unrealistically high in the context of the cuprates , if one considers more realistic bare fs , @xmath174 . furthermore , in the ic ddw case , there need not be a simple relation @xmath166 between the doping and doping and deviation from commensurability . let us now discuss the c - ic transitions that are present in our phase diagrams proximate to ` magic ' filling fractions @xmath175 . although our understanding as to why they arise is incomplete , we would like to point out several interesting features about our results . first , to our knowledge , there is nothing special about the bare fs close to these filling fractions , likewise for the fs of the commensurate ordered state . the c - ic transition is not driven by nesting , since upon variation of @xmath18 we see similar signatures proximate to the ` magic ' dopings . second , the magnitude of the ic wavevector is generally small ( corresponding to periodicity of more than @xmath176 lattice spacings ) and shows no clear dependance on @xmath177 , figs . [ fig : qdw4 - 2tp006 ] -[fig : qdw2 - 1tp012 ] . third , we checked that similar behavior is present deep into the af ordered state in the total free energy . [ fig : qaf004t0tp004 ] . unlike the ddw case , however , the ic lattice spacing in the af case depends on @xmath177 : @xmath22 follows an increasing staircase pattern at ` magic ' dopings . note that our results for af order seem to contradict schulz s results@xcite which predict a linear dependence of @xmath149 . the naive contradiction is easily resolved by realizing that schulz s computation@xcite of @xmath149 relies on the stoner criterion , which is strictly valid at @xmath10 . indeed , while deep into the ordered state the ic wavevector is stepwise increasing at ` magic ' fractions , at higher temperatures we regain the linearity of @xmath178 , fig . [ fig : qaf004t035tp004 ] . moreover , the fact that both the bare and the gapped fs in the c case are generally not nested at the the ` magic ' fractions suggests that the transition is not driven by a divergent susceptibility , but is first order . indeed , at low temperatures , the chemical potential near the special dopings is non - monotonic for both the af and the ddw+dsc diagrams . similar behavior of @xmath151 was found by nayak and pivovarov in the context of c ddw competing with dsc@xcite . a simple view of the energetics is given by the following landau - ginsburg free energy of two nodeless competing orders ( the generalization to order parameters with nodes is straightforward ) : @xmath179 in this free energy , @xmath180 plays the role of the ddw order parameter and @xmath181 the role of superconductivity . when @xmath182 becomes negative , @xmath180 orders at non - zero wavevector @xmath183 . this occurs at the dopings indicated in figures [ fig : ddwpd3 ] and [ fig : ddwpd4 ] . when @xmath181 orders , this further favors incommensurability if @xmath184 since the competition between the two orders is partially alleviated . this alleviation is a possible explanation of why the onset of dsc at low @xmath152 is correlated to the presence of ic - ddw .     and @xmath151 ( inset ) at @xmath136 for @xmath185ev , @xmath186ev . note that @xmath187ev which is smaller than the value taken in fig . [ fig : ddwpd2].,height=288 ]     except this time @xmath143ev as in fig . [ fig : ddwpd2].,height=288 ]     and @xmath149 ( inset ) at @xmath150ev corresponding to the parameters in fig . [ fig : ddwpd3 ] . although here the ic response is not at ` magic ' fractions , note that in fig . [ fig : ddwpd3 ] there is dsc both at @xmath164 and @xmath165 . , height=288 ]     ( and at the same temperature ) except this time there is a clear c - ic transition at ` magic ' fractions ( inset ) . at the same time @xmath151 is non - monotonic at these dopings . the coupling constants @xmath188ev and @xmath189ev are the same as in fig . [ fig : ddwpd4].,height=288 ]     ev ) staircase ic spectrum of af around ` magic ' filling fractions . the parameters we have taken here are @xmath190 ev , @xmath142 ev . inset : @xmath151 is non - monotonic at both jumps of @xmath149 , but is smoothly decreasing whenever @xmath22 is pinned.,height=288 ]    the same phenomenon can be viewed from the opposite perspective : when ddw order becomes incommensurate , there are regions in which it is suppressed ( as in the domain walls of fig . 2 ) . superconductivity can occur more easily there . hence , incommensurability alleviates the competition between dsc and ddw order at the mean - field level , thereby enhancing superconductivity . of course , ic order appears to suppress superconductivity in the cuprates , but this could be due to effects beyond a simple mean - field theory .    when the superconducting order parameter becomes strong enough , it may become energetically favorable for phase separation to occur . we can guess when this occurs from the usual maxwell construction whenever the mean field chemical potential @xmath151 does not vary monotonically . this appears to occur at a few isolated ` magic ' doping values . however , in a real system , proximity to a second order critical point@xcite , coulomb repulsion@xcite , or disorder could drive the system in a state of mesoscopic phase separation . moreover , the lattice pins any ic density wave order , so once the system is ic in order to minimize the free energy , the deviation from commensurability @xmath22 would be pinned at a finite value . finally , it is important to remember that disorder could have important effects . it would not only disorder the array of domain walls , but could lead to droplets of one phase in regions of another .     regains its linearity around @xmath191 as predicted by schulz@xcite . here @xmath192ev and all other parameters are as in fig . [ fig : qaf004t0tp004 ] . inset : @xmath151 is monotonic throughout , suggesting that at higher temperatures , the c - ic transition is driven purely by ic nesting of the fs away from half - filling.,height=288 ]    we stress that while the general trend towards ic order may be clear , many of the quantitative aspects are not clear , in particular the precise copings at which ddw order becomes incommensurate . it is not clear to us why such c - ic transitions happen near ` magic ' filling fractions . the fact that similar behavior occurs in the case of af order suggests that the tendency of phase separation at ` magic ' filling fractions is a general trend for a variety of density wave order . because the af only case was explored upon taking the @xmath193 limit in the full free energy , it is not clear yet whether the tendency of phase separation at ` magic ' filling fractions is purely a feature of density waves competing with other order parameters , or of the lattice . it might be worth reconsidering the af case in the presence of other order parameters ( such as ddw , dsc ) . what makes our results surprising is that unlike in the frustrated phase separation picture for stripe formation , there is no coulomb repulsion present . hence , there is no obvious connection between phase separation and incommensurability . ( please note that while ic ddw order can induce charge density wave order at twice the wavevector , we have not included this effect , which we expect to be small , in our calculation . ) therefore , while coulomb repulsion may play a role in determining the cdw spacing at @xmath164 in @xmath194 , in our calculation , the phase separation at ` magic ' fractions can only arise as interplay between the lattice , competing order physics and energetics . what exactly is the combination of the three deserves some future research as the analog of the fulde - ferrel - larkin - ovchinnikov ( fflo ) state of dsc in a magnetic field parallel to the plane@xcite . in the absence of next - nearest neighbor hopping , the two mean fields are related via a particle - hole transformation in the same way _ commensurate _ ddw is related to dsc at half - filling . however , while in the latter case one can rotate c ddw to dsc@xcite , in the case of finite @xmath155 , the particle - hole transformation can not interchange fflo and ic ddw . while the mean - field hamiltonians of ic ddw and the fflo are related by a particle - hole transformation , the spectra are not  in the ic ddw state there is a splitting to a generally uncountably many bands ( if @xmath195 is irrational as we discussed in sec . [ num ] ) since the order parameter connects @xmath108 and @xmath196 while in the case of dsc there is no such splitting . we are indebted to the authors of @xcite for pointing out their work to us . ] . according to the discussion following it seems likely that within the dsc state a competing ddw order parameter is ic as it develops . hence , we briefly describe the experimental signatures that would distinguish between c and ic ddw order . originally , c ddw was also named ` hidden order ' in the context of the pseudogap , because most of the experimental signatures associated with it are indirect . the magnetic field , created by the alternating currents around neighboring plaquettes , is too small to be detected by neutron scattering ( @xmath197 g for @xmath198 mev ) . at the same time , because of the symmetry of the c ddw magnetic field , there is no signature in nmr . the charge of the c ddw state is uniform so stm would not show signatures , unless they are associated with interference of nodal quasiparticles due to scattering off impurities . it has been argued recently@xcite that such interference disperses weakly enough to account for the pseudogap d.o.s . incommensurate modulations observed by stm . for a further discussion of the signatures of c ddw , the reader is referred to the literature @xcite . an ic ddw is more easily detectable directly . the local electronic density of states would be modulated due to the ic order parameter , and hence could be observable by stm . secondly , the ic modulations of the staggered currents would produce nmr line splitting of both the cu and o atoms when the applied magnetic fields are _ perpendicular _ to the a - b planes . this is in sharp contrast with ic sdw , where there is an nmr line - splitting due to spatially modulated local magnetization , which is directed _ along _ the a - b plane . let us estimate the modulated magnetic field strength arising from ic ddw ( the reader is referred to@xcite where similar considerations are applied to the case of orbital antiferromagnetism proposed to exist in @xmath199 ) . suppose we have an array of domain walls , fig . [ fig : domains ] , along the y - axis . deep between the domains the field is zero as in the case of c ddw . along a domain at @xmath200 , where the field is maximal , again most of the fields cancel except those created by the vertical alternating currents of strength @xmath201 at @xmath202 and those of strength @xmath203 at @xmath204 . equivalently this would be the field of a single vertical array of alternating currents of strength @xmath203 at @xmath202 : @xmath205 where @xmath206/[\sqrt{((i-1)a + y_0)^2 + x^2_0}]$ ] and @xmath207 for o , cu atoms respectively . clearly at @xmath208 , cu , atoms the field is zero by symmetry . for o atoms the sum can be performed yielding @xmath209 . in the dilute domain limit , this field is independent of the distance between the domains .    for a sinusoidal current modulation , @xmath210 would depend on the modulation length @xmath211 . if @xmath212 then a good estimate of @xmath210 is given by with @xmath201 replaced by @xmath213 where @xmath177 is the periodicity of the modulation ( in units of the lattice spacing ) . therefore , for large @xmath177 we obtain @xmath214 . incommensurate order is an important possibility when competing orders are considered . at a first - order transition the system is phase separated and hence necessarily inhomogeneous . if the competition between the orders is weak , coexistence is possible , but modulations of the competing orders are likely to be induced . the effects of the lattice , disorder , lomg - range coulomb repulsion@xcite and proximity to criticality@xcite is to generally stabilize incommensurability on mesoscopic length scales . this is unlike the case of , say , the liquid to vapor phase transition where the competition results in phase separation on macroscopic length scales .    in the case of ddw and dsc , the competition is certainly weak within the framework of the extended hubbard models we have considered . because the processes stabilizing both orders are of similar nature and are at the same time smaller than all of onsite repulsion , exchange and hopping , it does not cost much energy to convert ddw order into dsc and vice versa . a simple resolution of the competition is for ddw order to become incommensurate . we find that this occurs for some doping levels , but the incommensuration tends to be relatively small . such ic ddw could , in principle , be observable at low temperatures by nmr since it creates an inhomogeneous magnetic field at o atoms _ perpendicular _ to the ab - plane . most of nmr experiments test for much stronger magnetic fields _ parallel _ to the ab - plane on the cu atoms , and at higher temperatures such ic ddw phase would be smeared by thermal fluctuations , which may exlpain why such ic order , if it exists , it hasnt been seen by experiment . we would like to thank s. brown , s. chakravarty , r. jamei , s. kivelson , and e. pivovarov for helpful discussions . special thanks to e. pivovarov for letting i. d. generalize some of his computer code . this work has been supported by the nsf under grant no . dmr-0411800 . 99 u. schollwck , s. chakravarty , j. o. fjaerestad , j. b. marston , and m. troyer , phys . lett . * 90 * , 186401 ( 2003 ) ; e. dagotto , rev . phys . * 66 * , 763 ( 1994 ) ; c. s. hellberg and e. manousakis , phys . lett . * 83 * , 132 ( 1999 ) ; s. sorella , g.b . martins , f. becca , c. gazza , l. capriotti , a. parola , e. dagotto , phys . lett . * 88 * , 117002 ( 2002 ) . s.  r. white and d.  j. scalapino , phys . b * 61 * , 6320 ( 2000 ) ; d. f. b. ten haaf , p. w. brouwer , p. j. h. denteneer , and j. m. j. w. o. putikka , m. u. luchini , and r. r. p. singh , phys . lett . * 81 * , 2966 ( 1998 ) ; l.  p. pryadko , s.  a. kivelson , o. zachar , phys . rev . lett . * 92 * , 067002 ( 2004 ) ; e. dagotto and t. m. rice , science 271 , 618 - 623 ( 1996 ) ; g. roux , s. r. white , s. capponi , a. laeuchli , d. poilblanc , phys . b 72 , 014523 ( 2005 ) . experiments in the cuprates have been interpreted in terms of competing orders discussed in refs . ; and also in c.  m. varma , phys . b 55 , 14554 ( 1997 ) ; s.  c. zhang , science 275 , 1089 ( 1997 ) . e. w. carlson , v. j. emery , s. a. kivelson , and d. orgad , in _ the physics of conventional and unconventional superconductors _ , k. h. bennemann and j. b. ketterson ( springer - verlag 2003 ) , cond - mat/0206217 ; a.  h. castro neto and c. morais smith , in _ strong interactions in low dimensions _ , ed . d. baeriswyl and l. degiorgi ( kluwer acad . pub . 2004 ) , cond - mat/0304094 ; v.  j. emery , s.  a. kivelson , and j.  m. tranquada , proc . usa * 96 * , 8814 ( 1999 ) ; m. vojta and s. sachdev , phys . lett . * 83 * , 3916 - 3919 ( 1999 ) ; j. zaanen , nature 404 , 714 ( 2000 ) ; and references therein . j.  m. tranquada , physica b 241 - 243 , 745 ( 1998 ) ; in _ neutron scattering in layered copper - oxide superconductors _ , ed.a . furrer ( kluwer , dordrecht , the netherlands ) , p. 225 . j. e. hoffman , k. mcelroy , d .- h . lee , k. m lang , h. eisaki , s. uchida , and j. c. davis , science * 297 * , 1148 ( 2002 ) ; k. mcelroy , r. w. simmonds , j. e. hoffman , d .- h . lee , j. orenstein , h. eisaki , s. uchida , j. c. davis , nature * 422 * , 592 ( 2003 ) ; c. howald , h. eisaki , n. kaneko , m. greven , and a. kapitulnik , phys . b * 67 * , 014533 ( 2003 ) . s. chakravarty , r. b. laughlin , d. k. morr , and c. nayak , phys . rev . * b 63 * , 094503 ( 2001 ) . c. nayak , phys b 62 * , 4880 ( 2000 ) . t. timusk and b. statt , rep . progr . phys . * 62 * , 61 - 122 ( 1999 ) . s. tewari _ et al . _ , , 224516 ( 2001 ) . j.l . tallon and j.w . loram , physica c * 349 * , 53 ( 2001 ) ; j.w.loram _ et al . _ , j. phys . solids @xmath215,2091 ( 1998 ) ; j.l.tallon _ et al . b @xmath216 , 531 ( 1999 ) . y. j. uemura _ et al . _ , lett , @xmath217 , 2317 ( 1989 ) ; c. panagopoulos _ et al . b * 60 * , 14617 ( 1999 ) ; c. panagopoulos _ et al . _ , physica c , * 341 - 348 * , 843 ( 2000 ) . s. chakravarty , c. nayak , s. tewari , , 100504/1 - 4 ( 2003 ) . homes _ et al . _ , , 1645 ( 1993 ) . c. stock _ et al _ , , 024505 ( 2002 ) . h. a. mook _ et al . _ , , 144513 ( 2002 ) . s. chakravarty , h. y. kee , and c. nayak , _ intl . j. phys . _ * 15 * , 2901 ( 2001 ) . s. chakravarty , c. nayak , s. tewari , x. yang , , 277003/1 - 4 ( 2002 ) . et al . _ , nature * 424 * , 912 ( 2003 ) . s. tewari _ et al . _ , , 014514 ( 2004 ) . _ , , 137002 ( 2004 ) . c. wu , j. zaanen , s - c zhang , , 247007 ( 2005 ) . u. schollwoeck _ et al . _ , 186401 ; j. o. fjaerestad , j. b. marston , u. schollwoeck ann . ( n.y . ) * 321 * , 894 ( 2006 ) . ; s. capponi , c. wu , s - c . zhang , 220505(r ) ( 2004 ) . e. cappelluti and r. zeyher , phys . rev . * b 59 * , 6475 - 6486 ( 1999 ) . kee and y. b. kim , phys . rev . * b 66 * , 012505 ( 2002 ) . c. nayak , e. pivovarov , , 064508 ( 2002 ) . h. j. schulz , j. de physique * 50 * , 2833 ( 1989 ) . s. chakravarty , phys . rev.*b 66 * , 224505 ( 2002 ) . a. kotani , j. phys . jpn , * 42 * , 408 ( 1977 ) . l. pryadko , _ et . al . _ , phys . rev . * b 60 * , 7541 ( 1999 ) . a. damascelli , z. hussain , z - x shen , rev . phys . * 75 * , 473 ( 2003 ) . s. kivelson , _ et . al . _ , rev . phys . * 75 * , 1201 ( 2003 ) . r. jamei , s. kivelson , b. spivak , cond - mat/0408066 . c. bena _ et . _ , phys . rev . * b 69 * , 134517 ( 2004 ) ; c. bena _ et . _ cond - mat/0311299 ; a. ghosal , a. kopp , s. chakravarty , cond - mat/0412241 . v. tripathi , p. chandra , p. coleman , j. phys . : condens . matter 17 , 5285 ( 2005 ) .
Abstract:

Answer: we consider various incommensurate ( ic ) order parameters for electrons on a square lattice which reduce to @xmath0-density wave ( ddw ) order when the ordering wavevector @xmath1 . we describe the associated charge and current distributions and their experimental signatures . such orders can arise at the mean - field level in extended hubbard models . we compare the phase diagrams of these models with experiments in the underdoped cuprates , where ( 1 ) ddw order is a possible explanation of the pseudogap , and ( 2 ) there are experimental indications of incommensurability . we find various types of ic ddw and discuss their possible relevance to the physics of the cuprates . our main finding is that ic ddw order is generally accompanied by superconducting order , but the magnitude of the ic wavevector can be small . a comparison with the analogous af - icsdw transition is given .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: 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 .
Abstract:

Answer: the effective potential @xmath0 of a classical ion in a weakly correlated quantum plasma in thermodynamic equilibrium at finite temperature is well described by the rpa screened coulomb potential . additionally , collision effects can be included via a relaxation time ansatz ( mermin dielectric function ) . these potentials are used to study the quality of various statically screened potentials that were recently proposed by shukla and eliasson ( se ) [ phys . rev . lett . * 108 * , 165007 ( 2012 ) ] , akbari  moghanjoughi ( am ) [ phys . plasmas * 22 * , 022103 ( 2015 ) ] and stanton and murillo ( sm ) [ phys . rev . e * 91 * , 033104 ( 2015 ) ] starting from quantum hydrodynamic theory ( qhd ) . our analysis reveals that the se potential is qualitatively different from the full potential , whereas the sm potential ( at any temperature ) and the am potential ( at zero temperature ) are significantly more accurate . this confirms the correctness of the recently derived [ michta _ et al . _ , contrib . plasma phys . * 55 * , ( 2015 ) ] pre - factor @xmath1 in front of the bohm term of qhd for fermions .