quark - antiquark annihilation subprocess of the hadron - hadron collision , mentioned in the abstract ( see @xcite and @xcite ) , may provide an interesting information about the quark dynamics inside the hadron .
the results of measuring of lepton pairs , produced in these collisions at high energies @xcite , were already used to enlarge our knowledge about the structure functions that describe the parton distributions in hadrons @xcite .    due to a limited size of a range of avaliable total energy panda
would be rather an intermediate than a high energy experiment . in this energy region
the important role plays the physics of bound states and resonances , which is strongly connected with the quark confinement problem .    at the same time
it is well known from the previous @xmath5 experiments , done at low energies ( like the @xmath5 annihilation experiments done at lear at cern ) , that many features of the resonance decays may be also quite well explained , basing on the simple phenomenological rules , like the so called quark line rule @xcite .
recently , the obelix collaboration ( cern ) has published the results of a detailed analysis of three reactions of the annihilation at rest with the three mesons production in the final state : @xmath6 @xcite , @xcite .
the analysis done in these articles also does not exclude the possibility that some features of the observed resonant structures , or enhancements , as well as the observation of some resonant structures suppression in the s- and p - waves of the @xmath5 system , which undergo the annihilation process at rest , may need the @xmath7 model application for their explanation .
from here we may conclude that if we shall consider a case of the antiproton scattering off the proton target in a case of the antiproton beam of 14 gev energy ( that corresponds to the center of mass energy of the @xmath5 system @xmath8 ) , then it would be quite natural to use a quark parton model to describe the @xmath9 process .
to this reason we have used here the pythia event generator @xcite , which is based on the ideas of the quark parton model and is well tested and widely used for the simulation of high energy hadron - hadron interactions .
some parameters of the generator were set here , following @xcite , to those values that allow fast and reliable simulation in the case of @xmath10=14 gev .
it is important to note that pythia simulation is based on the use of the amplitudes of the relativistic quantum field theory . in our case
it is the perturbative qcd amplitude of the @xmath0 subprocess , implemented into the pythia package .
this allows a proper account of the relativistic kinematics while simulation of the distributions of different physical variables , connected with the @xmath11- ( or @xmath12- ) pairs produced in event .
it is worth to mention also that in any case the dominant mechanism of the lepton - anti - lepton pair production would be one and the same in the most of another than perturbative qcd approaches or models .
really , in this special case of @xmath13 process the dominant contribution to the process amplitude in any model would come from one and the same amplitude of quark pair annihilation @xmath14 in the chosen range of beam energy value .
it should be mentioned that pythia generates the events in which the final state particles are produced in a result of a rather long chain of interactions .
this chain includes the primary partons , which are produced in a hard perturbative @xmath15 quark - parton fundamental subprocess and latter on do pass the stage of interaction between themselves ( which is described by the lund string fragmentation model @xcite - @xcite used in pythia ) .
after this the quarks , diquarks and strings , produced at this stage , do pass the hadronization phase .
the variety of the produced elementary particles then undergo a stage of decays in a case of the presence of the nonstable particles or resonances among them . at the end of the simulation there
are left only free particles in the final state . therefore , while using the pythia , we may get , as the result of simulation , only the distributions of free particles , unbiased by any detector effects like , for example , by a magnetic field .
the account of the influence of the detector effects would be a subject of the following papers .    in the section 2 , devoted to lepton distributions , the results of the previous note @xcite will be added by a number of new 3- and 2-dimensional plots that show different kind of correlations ( like energy - energy and angle - angle correlations ) between the physical variables of the leptons , produced in the quark level subprocess @xmath0 .
all plots of this note are collected in the appendix ( section 7:``figures '' ) .
let us mention that the main results of this note ( sections 2 , 3 and 4 ) are valid for a case of muon pair as well as for the electron - positron pair production .
to this reason in these three sections we shall speak about the distributions obtained for the lepton case because muon distributions are practically do not differ from the electron ones ( so , the word `` muon '' is equal to the `` electron '' in the sections 2,3 and 4 ) .
the main difference of @xmath12 case from that one of @xmath16 pair production would appear only while considering the `` background '' contribution from the fake muons ( see sections 5 and 6 ) .    the next important new addition , as comparing to @xcite , is presented in the section 3 , where the distributions of some global variables , that describe the produced muon or electron pair as a whole system , like for instance , the invariant mass of the signal lepton pair . among them of most interest is the total transverse momentum of a lepton pair , which value is connected with the intrinsic transverse velocity of a quark inside the proton .
recent results of the experiment e706 on the direct photon production , performed at fermilab @xcite , have shown that the averaged value of the quark transverse momentum in the proton may be of @xmath17 gev order .
two sources of a possible appearance of quark transverse momentum in the proton , i.e. the intrinsic fermi motion of quark inside the proton and the influence of an initial state radiation , are studied in section 3 .
it is pointed out that due to the fact that the transverse motion of quark in the nucleon does not depend on the value of the boost along the beam line , there is a unique possibility for panda to measure the value of quark @xmath2 , which is so important for a precise setting of the absolute scale of jet energy @xcite - @xcite and thus for a precise determination of particles mass ( like of top quark and may be of higgs boson and susy particles ) in the experimental measurements at tevatron and lhc .
let us noted that each generator has its own internal restrictions .
in present generation we had , due to the internal pythia cuts , a restriction on the value of the invariant mass , produced in a fundamental @xmath18 hard qcd subprocess . in practice ,
the value of the invariant mass of a muon ( electron ) pair , produced in a signal event , turned out to be restricted by the limit from below : @xmath19 gev .
the influence of this cut on the spectrum of the signal lepton pair energy will be discussed in the sections 3 and 4 .
the problem of fake muons background is discussed in the last section 5 .
we do not consider here a case of a contribution from the events , defined by another than @xmath20 subprocess , because , as far as we have an aim to study the muon pair production ( i.e. , a case of muon pair which has a common vertex ) the appearance of this pair from some vector mesons decays would not be so large .
nevertheless we plan to study it in the nearest future .
the main problem is caused by the presence of the fake muons from pions decays in the same @xmath5 process ( based on the annihilation @xmath20 subprocess ) .
it is clear that the appearance of this kind of fake electrons in the same signal annihilation @xmath21 process has a negligibly small probability as comparing with the case of fake muons from pion decays . ]    as the first step the number of events that would not contain at all the charged pions in the final state is estimated in section 5 .
the simulation has shown that the fraction of such events without fake decay muons may turn out to be not a small one .
the distribution of the number of fake muons versus their pt value , presented in section 5 , provides an information for the proper choice of an effective additional cut for their discrimination .
section 5 includes also the distributions of x-,y- and z - components of the position of the vertex of the fake muons production , which may be also useful for their separation . in section 6 the spectra of background leptons , mainly electrons , produced in other than @xmath0 qcd subprocesses ,
are presented .
we have done the generation of 100 000 `` @xmath9 '' events with @xmath18 qed @xmath23 subprocess . in what follows these events
should be called as `` signal events '' , and , respectively , the muons ( electrons ) that are produced in hard qed @xmath24 subprocess would be called as the `` signal '' leptons . those leptons that appear in events due to hadron ( mainly mesons ) decays will be called as `` decay '' leptons .
the distributions of the number of generated signal events versus the values of the energy @xmath25 of the signal leptons , as well as versus the modulus of the transverse momentum @xmath26 and of the polar angles @xmath27 , measured from the beam direction ( they are denoted in plots as @xmath28 ) are shown ( from top to bottom , respectively ) in the figure 1 of the section 7 ( appendix : `` figures '' ) , which contains all plots and figures .
the left hand column in the figure 1 is for @xmath29 distributions and the right one is for @xmath30 .
no difference is seen between @xmath30 and @xmath29 distributions .
one can see from the top raw of figure 1 that the energy of leptons may vary in the interval @xmath31 with a mean value @xmath32 gev and its spectrum has a peak at a rather small value of about @xmath33 = 0.4 gev .    the @xmath34 spectrum ( see the middle raw of figure 1 ) has an analogous peak position at @xmath35 gev . behind their peaks both @xmath36 and @xmath34 spectra falls rather steeply but the spectrum of @xmath34 is confined in a more narrow interval @xmath37 gev .
the number of events ( @xmath38 ) spectrum versus the polar angle @xmath39 ( see bottom raw ) has a peak around @xmath40 and the mean value @xmath41 .
we see that while the most of signal leptons are scattered into the forward direction ( @xmath42 ) still there is a small number of them which may fly even into the backward hemisphere ( @xmath43 ) .
it is useful to have the set of plots done separately for the signal leptons having the largest energy ( we shall call them as `` fast '' leptons in what follows ) in the lepton pair , produced in an event , and , correspondingly , for those , having a smaller energy ( `` slow '' lepton in a pair ) .
this set of plots for the signal leptons is given in figure 2 of the appendix , where the left hand column shows the distributions of `` slow '' leptons while the right one is for the `` fast '' leptons distribution .
one can see from figure 2 that the energy spectrum of the fast signal leptons ( top raw ) grows up rather fast from the point @xmath44 ( more than @xmath45 of fast leptons have @xmath46 ) up to the peak position at the point @xmath47 = @xmath48 and then practically vanishes at @xmath49 = @xmath50 .
in contrast to this picture the spectrum of the less energetic signal leptons ( left raw ) has a peak around the point @xmath51=0.4 gev ( where the spectrum of fast leptons only starts ) and then it vanishes at the point , which practically do corresponds to the mean value of the most energetic signal leptons , i.e. at @xmath52 = 3.8 gev .
one may see that the spectrum of slow leptons in a pair is a very different from that one of fast leptons .    the difference between the @xmath34 spectra ( see the middle raw of figure 2 ) of fast and slow leptons seams not to be so large .
it results only in about of 340 mev shift to the left of the peak position of slow leptons spectrum as well as in the corresponding end point shift to the left .
both of these spectra do demonstrate that the main part of slow as well as of fast leptons have @xmath53 gev .
this similarity of @xmath34 spectra of fast and slow leptons results ( due to a large difference of their energy spectra , i.e. of @xmath54 components ) in a large difference of their polar angle @xmath39 distributions ( see bottom raw in figure 2 ) . really , the @xmath55 spectrum is shifted to the right to the higher values , as comparing to that one of @xmath56 .
note also that their mean value @xmath57 = @xmath58 is a more than two times higher than the analogous mean value of the fast leptons : @xmath59 .
the most important difference that is seen from these angular distributions is that all fast leptons fly in the forward direction ( @xmath60 ) and their spectrum practically finishes at @xmath61 , while about of @xmath62 of slow leptons have @xmath63 .
it is worth to mention that about of @xmath64 of slow leptons may scatter into the back hemisphere .
the number ( or percentage ) of signal leptons , scattered into the forward or back hemispheres , may be easily found from the 3-dimensional `` angle - energy '' correlation plots , presented in figures 3,4 and 5 ( see section 7 , appendix ) done by help of root @xcite ) .
the left hand side of figure 3 contains the `` lego1 '' root plot with the distribution of the number of generated events versus : + @xmath65 the energy @xmath51 of the less energetic leptons in @xmath1 pair ( the energy axis is on the right hand side of this plot and its scale is given in gev units ) , as well as versus : + @xmath65 the polar angles @xmath55 ( the angle axis is on the left hand side of a plot and the scale on it is given in a degree measure units ) .
+ the right hand column contains the analogous `` angle- energy '' correlation plots with the distributions of the number of events versus the energy @xmath47 of the most energetic leptons in a @xmath1 pair , as well as versus their polar angles @xmath56 .
these plots show the correlation between the angle and the energy of slow and fast leptons in a @xmath1 pair in event , as well as the distribution of the number of events for the correlated values of these variables .
it becomes much more easy to see these numbers after making the projection of these surfaces onto a plane .
the projections of these 3-dimensional distributions onto the 2-dimensional @xmath66 -planes are shown in figures 4 and 5 .
the correspondence between the color contours in figure 4 ( surf2z plots ) and the number of events per unit of a bin square is defined by the vertical color scale , included into the right hand side of each plot . for a convenience of counting of the number of events in each bin square of the 2-dimensional @xmath39/@xmath36 planes we present in the figure
5 the `` text '' root analogs of plots from figure 4 .
the exact numbers in each square bin may be seen by a help of the zoom application to this place of .ps file plots . from these plots
one can see that the most of slow leptons that fly into the back hemisphere , i.e. @xmath67 , ( on total it is about about of 5.5% of events with @xmath1 pairs ) , do have the energy which is less than @xmath68 gev .
the number of events distribution versus the energies of leptons in a pair ( slow and fast leptons in a pair ) , i.e. that come from pions decays .
the `` energy - energy '' correlation , as well as versus the angles of both leptons in a pair , i.e. the `` angle - angle '' correlation , are shown in the 3-dimensional `` lego1 '' root plots in figure 6 of section 7 ( appendix ) .
the projections of these 3-dimensional plots onto the energy - energy @xmath69 plane ( left hand side ) , as well as onto the angle - angle @xmath70 plane ( right hand side ) , are shown in figure 7 ( as for the root `` surf2z '' color contour plots ) and in figure 8 ( as for the root `` text '' numerical plots ) .    from figures 7 and 8
we may conclude that among the signal events with @xmath1 pairs there are : + @xmath65 about of @xmath71 of them have lepton pairs with the polar angles belonging to the square ,     defined by the following bounds : +  @xmath72 ; + @xmath65 about of @xmath3 of events do have lepton pairs which angles are contained in a square : +  @xmath73 . + @xmath65 about of @xmath74 of events do have lepton pairs which angles fit into a more wide square : +  @xmath75 .
+ @xmath65 and there is about of @xmath76 of events , which angles are restricted by the relations : +   @xmath77 . +
it means that in case , when the lepton system would cover the angle region @xmath78 , then less than @xmath79 of events , which do contain the @xmath1 signal pairs , would be lost .
+ from the figure 2 and the energy - energy correlation plots shown in the figures 7 and 8 one may see that : + @xmath65 if the following loose cut from below would be applied to the energy of both ( fast     and slow ) leptons in a pair ( recall about the peak at @xmath80 ) +   @xmath81 , + then more than @xmath82 of simulated signal events with lepton pairs would be lost , + @xmath65 but in a case of setting of a more strict cut : +   @xmath83 , + then about of @xmath84 of simulated events with signal lepton pairs would be lost .
so far we considered the distributions of physical variables that describe the individual features of leptons belonging to a @xmath1 pair . here
we shall discuss a set of variables which characterize the lepton pair system ( and the quark - antiquark system ) as a whole .
figure 9 of section 7 ( appendix ) contains the plots that present the distributions of the number of generated events versus the most typical physical variables that may describe the initial and final states of two particle systems which appear in the annihilation subprocess @xmath23 .
let us firstly mention those of them which include different components of the total 4-momentum @xmath85 ( @xmath86 ) of a lepton pair :    @xmath87    the distribution of the total energy of a lepton system    @xmath88    is shown in the right hand plot of the middle raw of figure 9 of appendix ( and noted there as the `` e - l+l- distribution '' ) .
the longitudinal component distribution    @xmath89    ( noted there as `` pz - l+l- distribution '' ) is included into the left hand side plot of the same middle raw .
these two plots look very similar .
the explanation of this fact follows from the form of the plot for the transversal component    @xmath90    namely , of its modulus @xmath91 which is presented as @xmath92 in the left bottom corner of figure 9 . from here
one may see that the transversal component distribution is , really , much more narrower ( its mean value is about 1 gev ) than those of @xmath93 and @xmath94 distributions . thus , the main contribution to the value of the lepton pair energy comes from the longitudinal component @xmath93 .
it is worth to mention that both of @xmath93 and @xmath94 distributions do not start exactly from the zero values .
the starting point of the total energy @xmath95 spectrum is around 1 gev.- and @xmath96 - distributions that are shown in the upper raw of figure 9 .
see also the introduction and the discussion below .
note also that , according to the top raw plots of the figure 1 , the energy spectra of individual leptons in a pair may start , practically , from the zero values . ] it is important to mention that the spectrum of a lepton pair energy in @xmath97 collision has the mean value around 5 gev ( as well as the spectrum of the @xmath93 component ) .
it means that on average the lepton pair caries a bit more than 1/3 of the total energy ( in the rest frame of the proton ) of the @xmath97-system . from the same right hand plot of the middle raw of figure 9
it is clearly seen that in about a half of events the energy of a lepton pair may be higher than 5 gev and it may spread up to the value of 12 gev .
another important variable for describing the partonic subprocess @xmath23 is the modulus of the @xmath98 4-vector , which we shall denote as `` q '' and which is defined as follows :    @xmath99    it is clear from ( 1 ) that q has the meaning of the invariant mass of two signal leptons pair , i.e. @xmath100 the left hand side plot in the upper raw of figure 9 includes the distribution of the number of events versus the values of this variable ( it is denoted there as @xmath101 ) .
let us look now at the set of variables that describe the initial quark - antiquark state of the annihilation subprocess @xmath23 . in analogy with ( 1 ) one may consider the 4-momentum of @xmath102 system @xmath103 and construct with its help the variable of the invariant mass of quark - antiquark system @xmath104 ( in pythia it is named as `` @xmath105 '' ) .
thus ,    @xmath106    its distribution is shown in the right hand plot of the same first raw of figure 9 @xmath107 .
both of the considered here invariant mass distributions , i.e. @xmath108= @xmath105 and @xmath109 , look very similar . from the plots of figure 9 it is seen that the most of events are concentrated near the starting point @xmath110 gev and there are , practically , no events after the point @xmath111 = 2.5 gev .
so , the spectra of the invariant masses of produced quark - antiquark system as well as of lepton pair are rather narrow .
the appearance of the left boundary points @xmath112 = @xmath113 = 1 gev stems from the internal pythia restriction ( see the discussion above and the introduction ) for the lowest value of invariant mass of two body system produced in fundamental quark - parton level @xmath15 subprocess .
the distributions of bjorken x - variables are shown in figure 10 for up- and antiup- quarks ( top raw ) and for down- and antidown- quarks ( bottom raw ) correspondingly .
they represent the up- and antiup- quarks components of the proton structure function .
the number of entries in these plots reflects the valence quark flavor structure of the proton ( the analogous distributions and the number of entries for the antiquarks , given in figure 10 , show the absolute similarity of quark and antiquarks distributions , obtained from the generated pythia samples ) .
now let us consider , following @xcite-@xcite and in a complete analogy with the equations ( 4)and ( 5 ) , the vector sum of quark and anti - quark transverse momentum vectors @xmath114 and the modulus of this vector sum @xmath115 the last variable has a distribution ( see right hand side bottom corner of the figure 9 where it is denoted as  pt - qqbar-@xmath116 ) which shape practically coincides ( in accordance with the transverse momentum conservation law ) with the shape of the distribution of the number of events versus the lepton pair transverse momentum @xmath117 , shown in the left bottom corner .
the variable @xmath118 is of a special interest because it contains an important information about two physical features of quark dynamics inside the hadron .
really , in our case when a beam antiproton moves along the z - axis and scatters off the proton fixed target ( we neglect the effect of target drops transverse motion ) the only sources of transverse motion of quarks may be as follows :    \a ) quark internal fermi - motion ( with some transverse velocity ) inside a + proton , i.e. the so called `` @xmath119 '' ;    \b ) initial state radiation ( isr ) of gluons ( or photons ) from quarks before + hard quark - antiquark annihilation .    the importance of these two mentioned above effects was recently discussed in a number of papers devoted to the study of `` @xmath120 '' events at lhc @xcite , @xcite and tevatron @xcite .
a lot of attention was paid to this problem also in @xcite in connection to the interpretation of direct photons production study in the experiment e706 at fermilab @xcite .    in the absence of these two sources the initial state quark and antiquark ( from the colliding antiproton and proton , respectively ) should have zero transverse momenta . here
we have used the possibilities of pythia generator which allow to switch `` on / off '' each of `` a)''- and `` b)''- effects separately as well as to switch them `` off / on '' simultaneously .
let us remind ( see section 2 ) that both plots in the third raw of figure 9 are done for a case when `` a)-effect '' ( i.e. `` fermi motion '' or `` @xmath121-effect '' ) as well as `` b)-effect '' ( or isr ) were switched `` on''(i.e .
, we have chosen the parameters @xmath122).-effect '' ( see , for instance , @xcite ) . ] .
firstly we have checked that in a case when both sources a ) and b ) were switched `` off '' the corresponding plots for @xmath123 distribution `` ) and @xmath124 distribution '' ) really show zero in their distributions ( the initial state quark and antiquark have zero transverse momenta and thus they may produce a lepton - antilepton pair in which the leptons may have either zero or back - to - back transverse momenta ) .
therefore , any deviation from zero of the _ total transverse momentum of a lepton pair _ , produced in a hard quark - antiquark annihilation , may appear only due to the presence of one of the two mentioned above features `` a ) '' or  b) of quark dynamics , or as the simultaneous result of both of them .    to study the influence of isr - effect we have switched it `` off '' ( i.e. we set mstp(61)=0 ) while keeping the parameter of fermi motion , i.e. `` @xmath121-effect '' , to be switching `` on '' , ( i.e. mstp(91)=1 ) .
no difference ( in the shapes as well as in the mean values ) was found between the plots for @xmath125 and @xmath126 , obtained in this case ( when only the effect of fermi motion or `` @xmath121-effect '' was taken into account ) with the previously discussed plots shown in the lowest raw of figure 9 done with the account of both effects ( mstp(91)=mstp(61)=1 ) .
so , one may conclude that , within pythia simulation , switching `` on / off '' of the initial state radiation does not produce any effect on the total lepton pair transverse momentum at the chosen beam energy ( 14 gev ) .
this result may indicate ( still keeping in mind that this feature is based on the isr model used in pythia , i.e. , it has a model nature ) that the phase space limitations at the chosen value of beam energy @xmath127 gev may allow to measure at panda the fermi motion effect separately from the isr effect .
further study in this direction is necessary , of course .
thus , future measurements of the total transverse momentum @xmath117 of the lepton pairs at panda energies may be of a big interest for a proper treatment of many qcd results from lhc ( and tevatron ) as well as by itself .    to this reason we may suppose that the fraction of events that have the final state consisting of a @xmath128-pair and a pair of two charged pions ( see also the discussion in the following section ) may be the most suitable for the measuring of the primordial `` @xmath121-effect '' ( caused by the intrinsic fermi motion of quarks inside the proton ) by measuring of the total transverse momentum of the lepton - antilepton pair in the final state .
the signal events , defined by @xmath0 subprocess , should contain also some hadrons in the final state .
the pions , produced directly or in the decays cascades of other hadrons , may decay in the detector volume and thus serve as a main source of background muons that may fake the signal muons , produced in the annihilation subprocess .
figure 11 includes ( like in the figures 1 and 2 ) the distributions ( from top to bottom ) of the number of events versus the energy @xmath129= @xmath130 , the transverse momentum @xmath131 = @xmath132 and versus the polar angle @xmath133= @xmath134 of produced pions .
the right hand side bottom plot contains the distribution of the total number ( npi ) of charged @xmath135-mesons in the signal events .
one may see from this plot that there is a big number of signal events ( about @xmath136 ) which do not contain any charged pions at all ( their final states mostly include the nucleon - antinucleon pairs ) .
this result is rather interesting .
let us underline that pythia provides a good ( at least one of the bests if not the best one ) but still a model approximation to the hadronization effects that take part in a reality .
the right hand side plot of figure 11 indicates that it may happen so that about of @xmath3 of events with signal muon pairs may appear practically without any additional fake decay muons content . from the same right hand side plot
we also see that about of @xmath82 of events have only one charged pion in the final state ( appearing mostly from @xmath137 resonance decays ) and a bit less than @xmath82 of events may have two charged pions .
it demonstrates also that about of @xmath64 of events may have 3 charged pions and there may be about of @xmath138 of events with 4 final state charged pions .
figure 12 contains the plots with the fake muon distributions in signal muon events .
its left column includes , respectively , from top to bottom , the energy @xmath139 , transverse momentum @xmath140 and polar angle @xmath141 distributions for muons that appear from hadrons ( mostly charged @xmath135 -mesons ) decays .
the shapes of distributions follow the analogues spectra of parent pions , presented in figure 11 . by comparing of these plots with those from signal muon pairs ( see figure 1 ) , one may conclude that the mean value of the energy of signal muons @xmath142 = 2.6 gev corresponds to the point in the energy distribution of fake decay muons where the contribution of fake muons ( in the same signal events ) is already very low .
analogously the mean value of the pt- distribution of signal muons @xmath143 gev ( see figure 1 ) do corresponds to the point where the spectrum of @xmath144 practically vanishes .
the e- and pt- spectra of decay fake muons are rather similar to those of `` slow '' muons from figure 2 , but a bit softer than the last ones .
they are shifted by 300 mev ( i.e. by @xmath145 to the left ( see their mean values ) as comparing with those of signal `` slow '' muons .
analogously the fake decay muons polar angle spectrum is shifted also to the left by @xmath146 , i.e. by @xmath147 . still , due to a big similarity of `` decay '' and
signal muons the former may appear as a serious background .
therefore a reasonable cut upon the muons energy @xmath148 ( see discussion in section 2 ) as well as on the @xmath144 value may lead to an essential reduction of the decay mesons contribution and allow to keep the main part of signal events .
thus , a cut @xmath149 gev may allow to get rid of a half of decay mesons in the signal events .    another way to discriminate up to a some level
the signal muons from the `` decay '' ones may be connected with the use of the information about the position of the fake muon production vertex .
the right hand side column in figure 12 contains ( from top to bottom ) the vx- , vy- and vz- ( i.e. along the beam direction ) components of the 3-vector that defines the position of the fake muon production vertex , obtained , as everything at the pythia generator level of simulation , for a space that is free of the detector equipment .
the distances at these distribution plots are given in millimeters ( mm ) .
therefore , one may see that the tails of the vz - distribution may expand up to 100 meters ( @xmath150 of events have @xmath151 m ) , of events the decay muons may appear even at more larger than 100 m distances . ] while those that are along x- and y- axises practically do finish at the 20 meters distance . from these right hand plots it is clearly seen that in a case of the signal events that do include the `` decay '' fake muons ( @xmath152 of all events ) the account of the real size of the detector `` decay volume '' may strongly reduce the number of charged @xmath135-mesons , that potentially may produce muons in the decays ,    to this reason , one may expect that in the real experimental conditions the situation with the contribution of the additional `` decay '' fake muons in signal events may turn out to be even much more easier .
only @xmath82 of events have the position of vz - component ( as it may be seen from the right bottom plot of figure 12 ) within the panda detector volume , i.e. @xmath153 m. at the same time it is also seen that in about of @xmath64 of events the vz componet may be very close to zero .
the corresponding fake muons produced near the interaction point may pose the most difficult background .
thus the detailed geant simulation , which is in our nearest plans , should allow to get more definite predictions about the decay mesons contribution to the background in signal events .
the contribution of the background muons from the other type of events ( qcd mainly ) should be discussed in the following section .
in this section we consider the distributions of leptons that may be produced in other than qed @xmath0 subprocesses . for this aim
the generation of 100000 events was done with the use of more than 20 qcd subprocesses ( including the signal @xmath0 subprocess ) that exist in pythia .
the results have shown that the main contributions come from the events that include the following ( presented in the decreasing order of @xmath154 ) partonic subprocesses : +   @xmath65 @xmath155 +   ( gives @xmath156 of events and the corresponding value of cross section @xmath157 mb ) ; +   @xmath65 @xmath158 +   ( gives @xmath146 of events and the corresponding value of cross section @xmath159 mb ) ; +   @xmath65 @xmath160 +   ( gives @xmath161 of events and the corresponding value of cross section @xmath162 mb ) ; +   @xmath65 @xmath163 +   ( gives @xmath164 of events and the value of cross section @xmath165 mb ) ; + the signal @xmath166 processes , that includes the qed partonic subprocess : +   @xmath65 @xmath167 +   has the value of @xmath168 mb what explains its very small contribution into the total number of generated events : there is only one event of this kind has appeared among 100 000 of generated events .
the shapes of the distributions of muons , produced in the mentioned above generated background qcd events ( see left column of figure 13 in appendix ) are not so much differ ( what is natural ) from those of fake `` decay '' muons produced in the signal @xmath169 processes , shown in figure 12 .
the main difference appear in the number of entries : for the 100 000 of generated background events there are 3.3 times more entries for muons ( 335372 ) , i.e. there are , on average , about of 3 muons ( that mostly appear from pions decays ) per one qcd background event .
it means that the probability of appearing of a muon pair with the different signs of their charges ( like in signal events ) is a rather high one in qcd background events .
so , they may fake quite well the signal events .
the distribution plots for the production vertices of background muons are shown in the left column of figure 13 . from the left column bottom plot of figure 13
it is seen that in about of @xmath170 of events there are @xmath171 meters .
thus , in these events the potential parent hadrons definitely should not produce muons within the volume of panda detector . at the same time
one may see that for a lot of events the muon production vertices are rather close to the interaction point : in about of @xmath172 of events the values of vz are practicaly equal to zero and on total about of @xmath64 of events have @xmath173 metre and @xmath82 of events have @xmath153 metres .
so , these fake decay muons may be separated from the signal muons latter on at the level of full geant simualtion with account of real detector conditions if , firstly , the information from the tracker or latter on ( and what is more likely may happen ) the information from the calorimeter should allow to identify the muon parent particles as hadrons .
so , it seams that the separation of qed signal and qcd background events would be a rather hard problem for the muon channel and a more detailed geant simulation of muon processes with an account of the real detector conditions ( material volume and the magnetic field ) would be needed .
now let us turn to the situation with the background for the electron production processes .
the distributions for the electrons , produced in 100 000 of generated signal events based on the @xmath174 subprocess , are presented in the figure 14 of appendix . from the number of entries
one may see that , on average , the enters number of fake electrons is 105283 , i.e. close to muon case ( see figure 12 ) .
the analysis of the generated event listings has shown that about of @xmath156 of generated signal events have no other electrons than those produced in signal @xmath12 pairs .
in the most of other events the electrons are produced in muon decays .
some electrons do appear also as the dalitz pairs , i.e. from the process of neutral pion decays @xmath175 . these pions in their turn may appear from the decays of @xmath176 and @xmath177 mesons or from more heavy mesons or barions , produced as the resonance states according to the lund fragmentation model . due to the fact that the life time of these resonance states is rather short ,
this kind of electron pairs are produced practically at the interaction point , what is reflected in the right hand side column of figure 14 , that contains the values of the components of the vector of the production vertex position .
this fact shows some disadvantage of e+e- channel as comparing with muon one where the most of background muons were coming from the light charged pions decaying in flight at rather large distances from the interaction point .
but it does not mean that the vertex position information turns not be so efficient for separation of signal and background electron events as comparing to the background muons case .
really , as it is seen from the right hand side bottom plot for the vz - component of the vector of the fake electron vertex position , only less than @xmath178 of fake electrons are produced near the interaction point position @xmath179 .
about of @xmath180 of events have the value @xmath181 metres ( `` overflow = 1.001e+05 '' ) , i.e. most likely are produced from the decays of fast muons ( less energetic muons may give about of @xmath164 of fake electrons inside the @xmath153 m volume of panda detector , as it is seen from the `` integral '' value in this right corner plot ) .
let us tern now to the case of electrons produced in the non - signal events .
the corresponding distributions , obtained from the same sample of 100 000 generated qcd background events , are shown in figure 15 of appendix . from the number of entries
one may see that , on average , the fraction of fake electrons in qcd background events ( 5487 ) is about of 60 times less than the analogous muons fraction ( 335372 ) ( see figure 13 ) .
this is an obvious advantage of the electron channel . from the left hand side column , that contains the values of the components of background electron vertex production position
, we see that the qcd background electrons are also produced mostly at the origin , i.e. from dalitz decays of neutral mesons like pion , @xmath176 , @xmath177 and others .
the number of `` overflow '' events ( 443 ) may be connected with the fast muon decays out of the detector size , ( i.e. with @xmath181 m ) .
now , taking into account that , as it became clear from the pythia event listings , in qcd events the decay electrons , that are produced within the panda volume @xmath153 m , do appear mostly in dalitz pairs , we may estimate the fraction of qcd background events , that really may pose the background @xmath182-pairs , as of about of @xmath183 order . as it is seen from the background electron energy- and pt - distributions , shown in the left hand side of figure 15 , these spectra fall much more steeply ( as pion energy in dalitz process have to be shared between 3 particles ) than the analogous background muon spectra of figures 12 and 13 .
to this reason the energy cut @xmath184 gev or the pt cut @xmath185 gev may kill a lot of background events .
nevertheless , these cuts may not be sufficient for the compensation of five orders difference in the rates of signal and background events for their clear separation .    a possible way to solve this task ( taking into account
the discussed above suppression of electron background in qcd events down to only @xmath183 content in generated events ) may be connected with the detector ( a tracker and an electromagnetic calorimeter ) possibilities of the identification of those electron pairs that may be associated with their parent decayed mesons by the reconstruction of the invariant mass of the system composed of the electron pair and of one of the registered photons , which may appear from the dalitz decay process @xmath175 .
this way , combined with the imposing of the cut from bellow on electron energy and its transverse momentum , may be even much more efficient in a case if the electromagnetic calorimeter possibilities would also allow to reconstruct the masses of the neutral pions from the other photon pair signals in event .
so , one may conclude from here that a more detailed study , based on full geant reconstruction with account of the detector conditions , would be very important for further study of electron pair production channel as well as for the muon one .
the sample of @xmath13 events , simulated by help of pythia event generator with the use of the hard partonic @xmath0 subprocess , is used to make a set of 50 plots which give a reasonable , but approximate , physical picture of the lepton pair production in a case when the antiproton beam of intermediate energy @xmath127 gev would scatter off the target proton .
the distributions of the most of essential kinematic variables of individual leptons , produced in the parton subprocess @xmath0 , are described in the section 2 and are presented in the figures 1 and 2 of the section 7 ( appendix ) .
they show that the energy and the angle spectra of the fast ( most energetic ) leptons in a pair are very different from those of slow leptons : the mean values @xmath186 gev are about three times higher than those of slow ones @xmath187 gev .
the simulation has shown also a tendency , which may be rather a general one : the fast leptons do fly predominantly at smaller angles @xmath188 as comparing to the angles of the slow ones @xmath189= @xmath190 .
it is worth to mention that about of @xmath79 of events may have slow leptons in @xmath1 pairs that may scatter into the back hemisphere , i.e. , @xmath67 .
the energy - energy and angle - angle correlations among a slow and a fast lepton in one and the same lepton pair in event are also described in section 2 and are presented in the plots of figures 3 - 8 in section 7 together with the corresponding distributions of the number of events versus the corresponding lepton energies and angles .
it should be noted that the present simulation , based on pythia , does not take into account , of course , the detector effects ( like the magnetic field , for example ) .
hence , the obtained plots , which describe the unbiased distributiond of free produced particles , may be useful as for the preliminary estimations mainly .
nevertheless , there is a big probability that the main features of the process would be similar to those presented at these plots .
thus , from the figures 2 and 8 it is clear that one may expect that the cut @xmath191 gev , which is a very reasonable due to the effect of the detector material , may lead to a loss of @xmath82 of events with signal lepton pairs , mainly due to the presence of the peak in the slow lepton energy spectrum in this region .
the data from the generated samples are planned to be used for a further geant simulation of these processes within the detector conditions .
the study of the 4-momentum components variables of the lepton pair system as a whole , done in the section 3 , has shown that the spectrum of the invariant mass of the lepton pair , which starts ( due to the already discussed in introduction restriction onto the value of the lepton pair invariant mass ) at the value of @xmath192 gev , decreases rather fast and do vanishes at @xmath193 gev .
at the same time the lepton pairs total energy @xmath95 spectrum starts at around 1 gev and spreads up to the value of 12 gev . in section 3
it was shown also that about a half of events do have the energy of a lepton pair higher than 5 gev .
so , in about a half of events the lepton pairs are rather energetic and they take away of more than 1/3 of the total energy of the produced @xmath1 system .    in section 4
it was pointed out that the total transverse momentum @xmath118 of the annihilating @xmath194-pair may be caused , in principle , by two sources
. it may appear due to the intrinsic transversal fermi motion of quarks inside the nucleon or in a result of the initial state radiation of gluons .
the results of the simulation have shown that , at least within the framework of pythia simulation , the initial state radiation ( isr ) gives a negligible contribution at the chosen value of energy .
therefore , it was argued that the study of the spectrum of the total transverse momentum of @xmath195 pair ( which coincides with the spectrum of the total transverse momentum @xmath118 of the annihilating @xmath194-pair due to the 4-momentum conservation law ) may allow to measure the intrinsic fermi motion of quarks inside the nucleon .
the results of such a measurement may be very useful for an interpretation of many qcd processes at future experiments at lhc as well as of the data collected at tevatron .
all what was said above and discussed in the sections 2 - 4 may be to a good approximation applied both to the case of @xmath11 as well as to the case of @xmath12 pair production in the final state    an important question of a background is considered in sections 5 and 6 . in section 5
we have concentrated on the additional muons that may appear from meson decays in the same signal process , defined by the mentioned above quark - antiquark annihilation to muon pair .
the problem is that the spectra of  slow  muons , which are present in the signal muon pair , are very similar to those of fake `` decay '' muons .
the most important observation , done in the section 5 , is that pythia predicts that there would be about of @xmath136 of events that would not contain at all the charged pions in the final state and thus they would be free of fake muons .
it is also shown there that the information about the coordinate position of the muon production vertex , as well as the imposing of the cut on the values of @xmath196 or of @xmath197 , may be useful for the discrimination of the fake muons .    in section 6
the study of background lepton events is presented .
it has been shown that in both muon and electron cases the separation of signal and muon background should not be a simple task and it needs more efforts .
it was pointed out that the electron pair production channel has an evident advantage as been compared with the muon one .
it is argued that the detailed geant simulation with an account of the detector design ( and based on the simulated and described here pythia event sample ) may be very useful .
the results of this simulation should be a subject of our following publications .
the authors are grateful to g.d.alexeev for the suggestion of this topic for study , the interest to the work and multiple stimulating discussions of the questions concerned .    50    v.a .
matveev , r.m .
muradian , a.n .
tavkhelidze , jinr p2 - 4543 , jinr , dubna , 1969 ; slac - trans-0098 , jinr r2 - 4543 , jun 1969 ; 27p
.    s.d .
drell , t.m .
yan , slac - pub-0755 , jun 1970 , 12p .
; phys.rev.lett . * 25*(1970)316 - 320 , 1970 .
cern ua1 collaboration , c.  albajar et al .
lett . , * b209 * ( 1988 ) 397 ; fnal e772 collaboration , p.l .
gaughey at al .
phys . rev . *
d50 * ( 1994 ) 3038 ;    a.d .
martin _ et al .
j. * c4 * ( 1998 ) 463 .
s.okubo , phys.lett .
b 5 , 165 ( 1963 ) ; g.zweig , cern report 8419/th 412 ( 1964 ) ; i.izuka , prog . theor .
phys ( suppl 37 ) 38 21 ( 1966 ) ;    m.bargiotti et al . ,
epjdirect a 2,1 ( 2002 ) ;    m.bargiotti et al . ,
j. c 35 , 177 ( 2004 ) ;    a.n.skachkova , n.b.skachkov hep - ph/0412279 , 2004 .
t.sjostrand , computer phys . commun.*39 * ( 1986 ) 347 , t.sjostrand and m.bengtsson , computer phys . commun.*43 * ( 1987 ) 367 .
r.  brun and f. rademakers , root - an object oriented data analysis framework , proceedings aihenp 96 workshop , lausanne , sep .
1996 , nucl .
inst . & meth .
* a389 * ( 1997 ) 81 .
see also http://root.cern.ch/.    fermilab e706 collab .
, l. apanasevich et al .
phys.rev.d70:092009 , 2004 , hep - ex/0407011 ; phys.rev.lett.81:2642-2645,1998 , hep - ex/9711017 .
skachkov , v.f .
konoplyanikov d.v .
bandourin , third annual rdms cms collaboration meeting .
cms - document , 1997168 .
cern , december 16 - 17 , 1997 , p.139 - 153 .
bandurin , v.f .
konoplyanikov , n.b .
skachkov , hep - ex/0207028 ;    d.v .
bandurin , n.b .
skachkov , d0-note-3948,fnal , mar .
2002 , hep - ex/0203003 , 2002 ; phys.part.nucl.35:66-106 , 2004 ( fiz.elem.chast.atom.yadra 35:113 - 177,2004 ) , hep - ex/0304010 ; l. apanasevich et al .
, phys.rev.d59:074007,1999 ; hep - ph/9808467 ; + j. huston ; int.j.mod.phys.a16s1a:205-208,2001 ;    b. andersson , g. gustafson and c. peterson , z. phys . * c1 * ( 1979 ) 105 ;    b. andersson , g. gustafson , z. phys . * c3 * ( 1980 ) 22 ;    b. andersson , g. gustafson and t. sjostrand , z. phys .
* c6 * ( 1980 ) 235 ; z. phys . *
c12 * ( 1982 ) 49 ;    b. andersson , g. gustafson and t. sjostrand , phys . lett . * b94 * ( 1980 ) 211 ;    b. andersson , g. gustafson , i. holgersson and o. mansson , nucl . phys .
* b178 * ( 1981 ) 242 ;    b. andersson , g. gustafson , g.ingelman and t. sjostrand , z. phys .
* c9 * ( 1981 ) 233 ;    b. andersson , g. gustafson and t. sjostrand , nucl .
* b197 * ( 1982 ) 45 ;    b. andersson , g. gustafson , g. ingelman and t. sjostrand , + phys .
* 97 * ( 1983 ) 31 ;
-5 mm   ( top raw ) , the modulus of the transverse momentum ( in gev ) @xmath26 ( middle raw ) and the polar angle @xmath27 ( bottom raw ) distributions . left column is for @xmath29 and right one for @xmath30.,title="fig:",width=680,height=755 ]    -5 mm   ( left column ) and the smaller energy ( `` slow '' leptons ) @xmath51 ( right column ) .
top raw includes their energies , in middle raw are @xmath34 and in bottom raw are @xmath198.,title="fig:",width=680,height=755 ]      -5 mm /@xmath36 plane . in each plot
the right hand vertical scale shows the one - to - one correspondence between the color and the number of events per square bin ( see fig.3 `` lego '' and fig.5 `` text '' plots).,title="fig:",width=680,height=755 ]        -5 mm /@xmath47
( left side ) and @xmath199/@xmath200 ( right side ) correlations but as been projected onto @xmath51/@xmath47 and @xmath199/@xmath200 planes .
the right hand vertical scale in each plot shows the one - to - one correspondence between the color and the number of events per square bin ( see fig.6 `` lego '' and fig.8 `` text '' plots).,title="fig:",width=680,height=755 ]      -5 mm   ( left column ) and of quark - antiquark system @xmath108 ( right column ) .
the middle raw contains the distributions of pz- and energy- components of @xmath201 system . in the bottom raw the total transversal momenta of lepton system ( left corner ) and of of @xmath194 system ( right side )
are shown .
, title="fig:",width=680,height=755 ]      -5 mm = @xmath130 , the transverse momentum @xmath202= @xmath203 and versus the polar angle @xmath133= @xmath204 of produced pions .
the right hand plot shows the distribution of the total number ( npi ) of charged @xmath135-mesons in the signal events.,title="fig:",width=680,height=755 ]    -5 mm , the transverse momentum @xmath205 and the polar angle @xmath206 of decay muons , produced in hadron decays in signal events that include the @xmath20 subprocess . the right column ( from top to bottom ) shows the distributions ( in mm ) of @xmath207- , @xmath208- and @xmath209- components of the position vector of a decay electron production vertex.,title="fig:",width=680,height=755 ]    -5 mm , the transverse momentum @xmath205 and the polar angle @xmath206 of muons , produced in qcd background events . the right column ( from top to bottom ) shows the distributions ( in mm ) of the @xmath207- , @xmath208- and @xmath209- components of the position vector of a decay electron production vertex.,title="fig:",width=680,height=755 ]    -5 mm , the transverse momentum @xmath210 and the polar angle @xmath211 of electrons , produced in hadron decays in signal events that include the @xmath212 subprocess . the right column ( from top to bottom ) shows the distributions ( in mm ) of @xmath207- , @xmath208- and @xmath209- components of the position vector of a decay electron production vertex.,title="fig:",width=680,height=755 ]    -5 mm , the transverse momentum @xmath210 and the polar angle @xmath211 of electrons produced in qcd background events . the right column ( from top to bottom ) shows the distributions ( in mm ) of @xmath207- , @xmath208- and @xmath209- components of the position vector of a decay electron production vertex.,title="fig:",width=680,height=755 ]


the present universe shows a variety of structures .
how such a structure is formed in the evolution of the universe .
one of the most plausible explanations is that the nonlinear dynamics of a self - gravitating system provides such a scale - free structure during the evolution of the universe . in order to clarify whether such a dynamics really gives an appropriate observed feature , we have so far three approaches : @xmath0-body simulation , the eulerian perturbation approach and the lagrangian one .
although the final answer for a structure formation would be obtained by the @xmath0-body simulation , it may be difficult to obtain enough resolution to discuss a fine structure of the universe such as a scaling property .
as for a perturbation approach , however , it is just an approximation and will break down in a nonlinear regime , although the lagrangian approach would be better if we are interested in density perturbations .
this is just because a density fluctuation @xmath1 and a peculiar velocity @xmath2 are perturbed quantities in the eulerian approach@xcite , while a displacement of particles from uniform distribution is assumed to be small in the lagrangian approach@xcite .
its first order solution is the so - called zeldovich approximation ( za ) @xcite .
the lagrangian approach is confirmed to be better than the eulerian approach by comparison of exact solutions in several cases@xcite .
therefore , we will adopt the lagrangian approximation in this paper and discuss about how to improve it .
in the standard approach of newtonian cosmology , the global cosmological parameters such as hubble expansion rate and mean density are given first by a solution of the einstein equations , i.e. the so - called friedmann - robertson - walker ( frw ) universe , which is an isotropic and homogeneous spacetime . according to observation
, however , a local structure in the universe is definitely not homogeneous and isotropic@xcite . in the standard approach ,
the density averaged over the whole space ( or a horizon scale ) is assumed to be the energy density of the frw spacetime .
however , here the problems of how to average inhomogeneous matter fluid and how to define an averaged isotropic and homogeneous spacetime are arisen .
this problem has been discussed in the framework of newtonian cosmology and general relativity by many authors @xcite .
the most important point to discuss is how inhomogeneities affect the expansion law of the universe and whether such a backreaction affects the evolution of the density perturbations .
this is our subject to discuss in this paper .
although many authors proposed different methods of averaging procedure in general relativity @xcite , a serious problem remains .
different gauge choices make the problem more complicated .
since a relativistic effect may not be so important , in this paper , in order to avoid such a difficulty , we discuss only an averaging procedure in the newtonian framework . proposing a averaging procedure , which is defined by spatial average of physical quantities , buchert and ehlers lead the averaged raychaudhuri s equation@xcite . the equation describes how the averaged expansion rate of domain with a finite volume evolves .
this equation has an additional term , which we call a backreaction term of inhomogeneities on averaged expansion . then , buchert , kerscher , and sicka estimated the backreaction term using the conventional lagrangian perturbation approach@xcite .
they first consider density perturbations in einstein - de sitter ( e - ds ) universe , then calculated the backreaction term and solved the averaged raychaudhuri s equation .
they showed difference between a cosmological parameter such as hubble expansion rate in their averaged model and that of the e - ds model .
although they included a backreaction term to estimate the averaged variables in a local domain , they used the perturbed quantities from the e - ds universe .
in other words , they have not take into account a backreaction on the evolution of density perturbations .    here
we improve their approach , i.e. we include a backreaction effect to averaged expansion rate and solve the averaged raychaudhuri s equation ( the generalized friedmann s equation ) with evolution equation of perturbations in the averaged domain . in an averaged domain ,
the averaged density is either higher or lower than that of the e - ds universe .
this difference will change the evolution of density perturbations .
in fact , if the domain is overdense , growth rate of perturbations behaves as that in the closed universe , while , if it is underdense , it is just like a solution in the open universe .    this paper is organized as follows : in @xmath3 , we shortly derive the generalized friedmann equation , following ehlers and buchert , and estimate a backreaction term .
we present our formalism in @xmath4 .
we show some relation between two lagrangian descriptions ( the conventional one and ours ) and give our initial setting in  .
4 . in  5 , we analyze a simple example , i.e . a plane - symmetric one - dimensional model , to show the probability distribution of the observed quantities such as hubble parameter , or density fluctuation , or peculiar velocity .
the conclusion and discussion follow in @xmath5 .
in the newtonian cosmology , the expansion of a domain is influenced by inhomogeneity inside the domain .
such an effect may be evaluated by spatial integration of field variables in the lagrangian domain , which evolves with matter fluid . hence in this paper
we study fields averaged over a simply - connected spatial lagrangian domain @xmath6 at time @xmath7 , which evolved out of the initial domain @xmath8 at time @xmath9 .
the locally averaged scale factor @xmath10 , depending on the content , shape and location of the domain @xmath6 , is defined by the volume of domain @xmath11 and its initial volume @xmath12 as @xmath13    we define a spatial averaging for any rank tensor field @xmath14 by the volume integral normalized by the volume of the domain as @xmath15 the average density is then given by @xmath16 where the total mass @xmath17 for a domain @xmath6 is conserved . using this averaging
, we can derive the generalized friedmann equation@xcite , which we will shortly derive here .
we start from the hydrodynamic equations for a self - gravitating perfect fluid ; @xmath18 although it is easy to extend the case with pressure of fluid , we consider just dust matter and ignore a pressure term .
to discuss inhomogeneity , a spatial derivative of field will play a crucial role .
in particular , it turns out that the spatial derivative of velocity * u * , which is divided into three variables as @xmath19 is very important , where @xmath20 the magnitudes of shear @xmath21 and of rotation @xmath22 are defined by @xmath23 using these variables , we find the basic equations for density , shear and rotation from eqs .
( [ eqn : conti - eq ] ) , ( [ eqn : euler - eq ] ) as @xmath24    the volume of a spatial domain @xmath25 described by its initial domain @xmath8 through a transformation from the eulerian coordinates @xmath26 to the lagrangian one @xmath27 as @xmath28 where @xmath29 is the jacobian of the transformation @xmath30 the change rate of the volume of a domain is then given as @xmath31 where @xmath32 denotes the lagrangian time derivative .
since @xmath33 is the averaged volume expansion rate of the domain @xmath6 , it may be natural to define the effective hubble expansion rate of the domain @xmath6 by @xmath34 . from the definition of the volume @xmath35
, we find @xmath36 .    the lagrangian time derivative does not commute with spatial averaging . for an arbitrary tensor field @xmath37 ,
ehlers and buchert showed the commutation rule : @xmath38    applying this commutation rule to the basic equations , we find @xmath39 replacing @xmath40 in eq .
( [ eqn : theta_evolution ] ) with @xmath41 , we find the averaged raychaudhuri s equation @xmath42 where @xmath43 describes the backreaction due to inhomogeneity of the universe . for @xmath44 case , the equation becomes the conventional friedmann equation with a scale factor @xmath45 for homogeneous and isotropic universe .
this is regarded as the generalized friedmann equation .      in order to evaluate the back reaction term @xmath46
, it may be convenient to introduce three principal scalar invariants as follows : for 2-rank tensor field @xmath47   in cartesian coordinates , those are defined by @xmath48 in particular , for the velocity gradient of matter fluid  @xmath49  , we find @xmath50 where we have used the following relations @xmath51    using those invariants , we can describe @xmath46 as @xmath52    in the following discussion , we will use the lagrangian description .
then , it will be convenient to rewrite eq .
( [ eqn : qdineuler ] ) by the lagrangian variables .
the lagrangian domain is just the initial domain @xmath8 .
the spatial average of a tensor field @xmath53 is described by the lagrangian coordinates as @xmath54 where the lagrangian average @xmath55 is defined by @xmath56 using this relation , the backreaction term is described by the volume average in the lagrangian domain @xmath8 as @xmath57    so far , we have not made any approximation .
however , when we evaluate those invariants explicitly , we need a further ansatz . here
we adopt a perturbation method . in the conventional perturbation approach ,
we expand the variables around the background friedmann universe .
in fact , buchert , kerscher , and sicka evaluated the backreaction term @xmath46 by this perturbation method@xcite , which we will show in the next section .
we have also another possibility to divide a perturbed part from a uniform background part as follows : since we average the variables in some domain and take into account its backreaction effect , we can adopt the averaged variables as unperturbed parts . this way to extract
the perturbed parts seems to be more consistent and may provide some difference from the previous approach by buchert et al .
we will discuss the detail next .
the conventional lagrangian approximation ( for example , za @xcite ) was constructed in a homogeneous and isotropic universe . however , when a nonlinear structure is formed , it is not so clear what is the background universe and how to divide the perturbed part from the unperturbed one .
although the simplest way is the conventional approach , as we discussed in previous section , a local expansion rate of the universe in some finite domain will be affected by its inhomogeneity
. then the growth rate of density fluctuations may also be affected by the inhomogeneity .
we shall derive the consistent equations for the lagrangian perturbations with the backreaction term @xmath46 .
first , we introduce the  comoving " eulerian coordinates @xmath58 and peculiar velocity @xmath59 as @xmath60 in the conventional lagrangian perturbation approach , we set @xmath61 where @xmath62 is just a scale factor of the background hubble flow of the whole universe and satisfies the friedmann equation .
however , because of the above reason , the expansion of a domain @xmath6 may be described by @xmath10 .
hence , we have not fixed a scale factor @xmath63 yet .
it will be determined later by adding a further ansatz .    with the relations @xmath64
the basic equations ( [ eqn : conti - eq])-([eqn : poisson - eq2 ] ) with @xmath65 are rewritten as @xmath66 where the density @xmath67 is divided into an unperturbed uniform part @xmath68 , which is defined by the mass conservation @xmath69= constant , and its fluctuation @xmath1 defined by @xmath70 @xmath71 is the gravitational potential defined by @xmath72 .    from eqs .
( [ eqn : euler1 ] ) and ( [ eqn : poisson1 ] ) , we obtain the following equation : @xmath73 if we specify @xmath63 , this equation provides the equation for inhomogeneous part subtracted a background uniform expansion . for example , in the conventional lagrangian approximation , @xmath63 is taken to be a scale factor of the frw universe @xmath62 , which satisfies the friedmann equation @xmath74 where @xmath75 is the mean energy density of the universe .
on the other hand , if @xmath63 is chosen to be @xmath10 , which evolution is given by the averaged raychaudhuri s equation ( [ eqn : g - friedmann ] ) , the perturbation equations are modified because of the backreaction effect @xmath46 as @xmath76 where the variables with subscript @xmath6 are defined those with the scale factor @xmath10 and @xmath77 .    the rotation of ( [ eqn : euler1 ] ) yields @xmath78    now we will describe those equations in terms of the lagrangian coordinates . using the lagrangian time derivative along the fluid flow @xmath79 eqs . ( [ eqn : evolution2 ] ) and ( [ eqn : rot ] ) are rewritten as @xmath80 where we have introduced @xmath81 , which is zero for the conventional lagrangian approach without backreaction but is chosen to be @xmath46 when we include the backreaction due to inhomogeneity .
note that there are still some space to estimate this backreaction term depending on a choice of the background scale factor @xmath63 .
the lagrangian coordinates @xmath82 , which follows the background flow , is defined by the initial values of the comoving eulerian coordinates @xmath58 , and the lagrangian perturbations are described as @xmath83 where @xmath84 denotes the displacement from the uniform distribution assuming that the scale factor is given by @xmath63 . the continuity equation ( [ eqn : conti1 ] ) or equivalently the mass conservation yields @xmath85 where @xmath86 is the jacobian of the coordinate transformation @xmath87 . with @xmath69= constant
, we have @xmath88 or equivalently for density contrast @xmath89    now all physical quantities are found to be written in terms of @xmath84 and it remains only to find solutions for @xmath84 .
we obtain the following equations for @xmath84 from equations ( [ eqn : longitudinal ] ) and ( [ eqn : transverse ] ) : @xmath90 @xmath91    in what follows , the subscript @xmath92 will be dropped except for @xmath63 and @xmath93 , but we may put it when we have to distinguish them . in order to write eqs .
( [ eqn : poissonl ] ) and ( [ eqn : curlfreel ] ) in terms of @xmath94 derivative , we use the relation between @xmath95 and @xmath96 , which is given by the definition as @xmath97 using this equation iteratively , we have @xmath98    this gives the differential equations in terms of @xmath27 formally . in order to write the basic equations explicitly , however , we need perturbative approach .
we then expand the deviation vector @xmath99 as @xmath100 .
superscripts @xmath101 and @xmath102 denote first and second order quantities in a perturbative expansion with respect to magnitude of small displacement @xmath103 from uniform distribution , respectively . from eq .
( [ qd ] ) , the leading terms of the backreaction is second order of @xmath103 .
then we have to derive the perturbation equations until at least second order . from eq .
( [ eqn : transverse ] ) , we obtain the first and second order equations as follows : @xmath104    next we consider equation ( [ eqn : longitudinal ] ) .
the jacobian @xmath105 is expanded as @xmath106 + \cdots \ , .
\label{eq : jacobian}\ ] ] thus we obtain the first order equation , @xmath107 and the second order one , @xmath108 in order to solve the above perturbation equations , we decompose @xmath109 and @xmath110 into the longitudinal and the transverse parts in the form @xmath111 where @xmath112 and @xmath113 are first - order and second - order scalar functions , respectively . as for their physical meanings , the first - order longitudinal and transverse parts are related to linear density and vortical perturbations , respectively . for the second - order level , however , such a simple interpretation of the perturbation modes does not hold any more .    the first - order perturbation equations ( [ eqn : first1 ] ) and ( [ eqn : firstt ] ) become @xmath114 the second - order perturbation equations ( [ eqn : second1 ] ) and ( [ eqn : secondt ] ) are also divided into the longitudinal and the transverse parts as @xmath115 @xmath116_i=4 \pi g \rho_b \varepsilon_{ijk } \mbox{\boldmath $ \psi$}^t_{l|j } \psi_{|kl } \label{eqn : second_trans}\;.\ ] ]    in what follows , we will fix @xmath63 to be @xmath10 .
note that if we set @xmath117 , the above equations are equivalent to those by buchert et al . since the displacement vector @xmath99 denotes a deviation from homogeneous distribution in the local domain @xmath6 ,
the average density @xmath118 , which is written by @xmath119 in what follows , is not equal to that in the whole universe .
the density fluctuation described by @xmath99 should be defined as @xmath120 here we use the subscript @xmath6 to show what is the background mean density .    to solve the above equations , we carry out further procedure .
since the first order equations ( [ eqn : first2 ] ) are the same as those of za except for a scale factor , we can easily solve them .
the perturbation variables are separated into time and spatial functions as @xmath121 the evolution equation for @xmath122 is @xmath123 while @xmath124 is determined from the initial data .    as for the second order equations , we need a little more complicated procedure . in appendix
a , we rewrite the backreaction term @xmath46 in context of the lagrangian approximation and present its explicit form up to the second order perturbations ( @xmath112,@xmath125,@xmath113 , and @xmath126 ) . since the transverse modes ( @xmath125 and @xmath126 ) may not be so important
, we shall ignore them in what follows .
then the leading order of the backreaction term @xmath46 becomes very simple as @xmath127 where @xmath128 ( see eq .
( [ eqn : la_qd_2nd ] ) . using ( [ eqn:1st - psi ] )
, we find @xmath129 which depends only on time @xmath7 .
( [ eqn : second_long ] ) for @xmath130 is now @xmath131 since the first term of r.h.s .
( [ eqn : second_long ] ) @xmath132 depends on spatial coordinate @xmath27 as well as @xmath7 , we have to divide a solution into two parts : one is inhomogeneous term and the other is homogeneous one .
first , in order to find a homogeneous term , we commute the time and spacial derivatives in eq .
( [ eqn : second - long ] ) , and then integrate it over @xmath8 , finding @xmath133 which is an ordinary differential equation for @xmath134 . by introduction of new variable @xmath135 by
@xmath136 we can eliminate the backreaction term @xmath46 .
the equation for @xmath135 is now separable by setting @xmath137 as @xmath138 where @xmath53 is a separation constant .
finally we obtain a set of basic equations ( [ eqn : g - friedmann]),([eqn : bd ] ) , ( [ eqn : average - zetad ] ) , and ( [ eqn : cd ] ) with the definition ( [ eqn : backreaction ] ) for dynamical variables @xmath10 , @xmath122 , @xmath134 , and @xmath139 . the spatial inhomogeneities
@xmath140 and @xmath141 are determined by initial distributions and eq .
( [ eqn:2nd_init_condition ] ) .
in order to solve our basic equations for lagrangian perturbations , we have to set up our initial data .
since the initial fluctuation is given by deviation from uniform distribution of the whole universe , which is described by the eulerian coordinates , we have to construct our initial data from those data by an appropriate transformation .
furthermore , since our equations are valid only in each domain @xmath6 , if we wish to analyze some statistical properties of our results , we have to go back to the whole universe . here
we will first discuss such a transformation , and then set up our initial data .
suppose that we have some inhomogeneous distribution in the eulerian coordinates @xmath26 .
we then find @xmath62 and @xmath10 by integrating a volume over the whole universe and over a domain @xmath6 , respectively . here
we set these initial values equal to unity .
these scale factors fix @xmath142 ( then @xmath143 and @xmath144 ) , and @xmath145 ( then @xmath146 and @xmath147 ) .
we also find @xmath75 and @xmath148 by averaging the density @xmath149 over the whole universe and over the domain @xmath6 , respectively .
the density perturbation in our basic equations is deviation from a density averaged in a domain @xmath6 .
hence a perturbation in the conventional lagrangian approximation does not satisfy the condition ( [ eqn : local - delta1 ] ) .
we have to reconstruct initial conditions in our local domain from fluctuations given in the whole universe . in order to set up initial conditions
, we may need only the relation between the eulerian coordinates @xmath26 and our lagrangian coordinate @xmath146 .
however , we also need the relation between two lagrangian coordinates , @xmath146 and @xmath143 from the following reason . in an inhomogeneous universe
, we do not know in which domain we are living .
then , when we analyze our results , we need a statistical analysis in the whole universe to know how the observed values are plausible .
in particular , one of the most interesting observed values is a peculiar velocity , for which we have to introduce a hubble flow @xmath62 and the comoving eulerian coordinate @xmath142 . as a result
, we need a relation between @xmath143 and @xmath146 as well .
note that although both lagrangian coordinates follow matter fluid , those scales are different . because we have mass conservation for infinitely small domain ( @xmath150 , or equivalently ( @xmath151 , i.e. @xmath152 assuming @xmath153 we find @xmath154^{1/3 }   \;. \label{eqn : qd - qh1}\ ] ] since the density perturbation in the conventional lagrangian approximation is given by @xmath155 we find @xmath156 from the definition of @xmath148 .
this gives the relation between @xmath157 and @xmath158 as @xmath159 we then find the relation between @xmath160 and @xmath119 as follows : from @xmath161 with eq.([eqn : qd - qh2 ] ) , we find @xmath162 } -1 \ ; , \label{eqn : convert - delta}\end{aligned}\ ] ] we also obtain a peculiar velocity against the hubble flow as @xmath163 we will use this definition in analysis of a peculiar velocity .
now we can set up our initial condition .
as for the density perturbations , we find from eq .
( [ eqn : convert - delta ] ) as @xmath164 } -1 \;. \label{eqn : initial_delta}\end{aligned}\ ] ] from this equation with definitions ( [ eq : delta ] ) , ( [ eq : jacobian ] ) , ( [ psi_def ] ) , and ( [ eqn:1st - psi ] ) , we find the initial condition for @xmath165 as @xmath166 where we set @xmath167 .    as for the scale factor @xmath10 of the domain @xmath6 , its initial value
is set to be unity ( eq.([eqn : initial_a ] ) ) , but its time derivative is non - trivial .
in fact , from ( [ eqn : local - h ] ) , we find @xmath168 as @xmath169 \;.
\label{eqn : ad - init}\ ] ]
here , we study a simple model to show new aspect in our approach .
we assume a plane - symmetric 1-dimensional model .
the lagrangian perturbation is given by @xmath170 in the conventional lagrangian approximation , za gives an exact solution in a plane - symmetric case .
however , it does not take into account a backreaction effect of inhomogeneity on the hubble expansion .
since our approach includes the backreaction effect , we will analyze this simple one - dimensional model and compare our results with those by za .
we also look at a difference of the backreaction term estimated by the conventional lagrangian approach ( buchert et al ) .
when we consider our backreaction effect , we have a set of equations : @xmath171    since we analyze the present model numerically , we have to introduce the size of `` whole universe '' , which is @xmath172 .
we the introduce the scale of averaging domain @xmath173 .
as for initial conditions , we adopt a power law spectrum with the index @xmath174 : @xmath175 we also introduce a cutoff at small scale , which wave number is @xmath176 , where @xmath177 .
we set that the initial time is @xmath178 .
then the amplitude of fluctuation is chosen so that the first shell - crossing occurs at @xmath179 .
the number of grids is @xmath180 and we use a periodic boundary condition .
we take an ensemble average over 500 samples , which initial conditions are given by random gaussian . in our approach , since we do not know where we are living , we study its statistical properties .
in particular , we will see the scale dependence of the averaged variables and the probability distribution of the hubble parameter , deceleration parameter , and pair - wise velocity .
first we analyze the expansion rate of local domain .
if we fix the hubble parameter @xmath181 by local observation in a domain @xmath25 , the most probable value of @xmath181 is given by the averaged expansion rate of the domain , which is @xmath182 .
several authors so far discussed such a local measurement of the hubble parameter @xcite . in particular ,
shi and turner @xcite estimated a possible value of the hubble constant measured locally and discussed a deviation from the global value , using linear perturbation theory with the cdm model .
they found that for small samples of objects that only extend to 10,000 km @xmath183 , the variance can reach 4% , while for large samples of objects to 40,000 km @xmath183 , the variance is about 1 - 2 % .
we solve eq .
( [ eqn : g - friedmann ] ) for @xmath10 with a backreaction due to inhomogeneity .
if we are living in an underdense region on average , the expansion rate will be faster than the hubble one for the whole universe . while , if we stay in an overdense region , the rate will be slower than the global hubble one .
fig.[fig : pdf - h0 ] shows the pdf of a local hubble parameter in our calculation .
if our domain is small , the deviation from @xmath181 gets large .
for example , the dispersion of the hubble parameter is about 1.2 % for the @xmath184-grid domain , while 0.66 % for the@xmath185-grid domain .
the dispersion of our model is consistent with the result by shi and turner .
next we show the pdf of density fluctuations . in the eulerian linear approximation ,
if initial data is given by random gaussian distribution , the pdf of density fluctuations will remain its gaussian form during evolution . on the other hand , in the lagrangian approximation
, there appears a nonlinear effect .
in fact , kofman et al shows that the pdf approaches to a log - normal function rather than a gaussian function in the cases of the lagrangian approximation and n - body simulation@xcite .
padmanabhan and subramanian also discussed the pdf with the za and found a non - gaussian distribution@xcite .    here
we analyze the pdf of density fluctuations using our approximation .
the results are shown in fig.[fig : pdf - delta ] . from comparison with the result of za , the void region ( i.e. an underdense region ; @xmath186 )
is found in higher probability in our approximation .
especially , if the size of a domain is smaller , the difference gets larger . on the other hand , the probability to find an overdense region ( @xmath187 ) decreases in our approximation
the reason is very simple : during evolution , an overdense region shrinks and a nonlinear structure is formed as the zeldovich s pancake . on the other hand ,
an underdense region expands .
therefore , although the initial volumes of overdense underdense regions are the same , the volume of the latter gets larger than that of the former in a nonlinear stage .
in addition to this lagrangian nonlinear effect , we take into account a backreaction effect .
this effect enhances expansion of an underdense region and contraction of an overdense region . as a result
, the above difference between za and our approximation appears .      in the conventional lagrangian approximation ,
if an initial condition is given by a random gaussian distribution , the pdf of a peculiar velocity also remains its gaussian form@xcite . because a peculiar velocity @xmath188 in the conventional lagrangian approximation ( or the za ) is given by @xmath189 the spectrum of a peculiar velocity is proportional to that of a density fluctuation as @xmath190 .
however , in our case , a peculiar velocity is given by ( [ eqn : convert - v ] ) and the growing factor @xmath122 is different in each domain , the pdf could deviate from a gaussian distribution .
however , from our numerical analysis ( see fig .
[ fig : pdf - v ] ) , the pdf of a peculiar velocity seems to still be a gaussian
. this may be because we need only a linear term of the lagrangian perturbations in the present analysis .
the pdf of a radial pairwise peculiar velocity is known to show an exponential form from analysis of n - body simulation and the za @xcite . for this variable , even if initial data is given by a random gaussian distribution , the pdf approaches an exponential form as the universe evolves@xcite .
the origin of this result could be understood by nonlinearity of gravity . in the case of one - dimensional plane - symmetric system in the conventional lagrangian approximation ( in fact
, the za is exact ) , we do not find any non - gaussian structure .
however , even in the za , the pdf shows non - gaussian behavior in the case of three - dimensional case@xcite .
if we take into account our backreaction effect , does non - gaussian behavior emerge even in one - dimensional case ?
the pairwise peculiar velocity is defined as follows : @xmath191 where @xmath192 and @xmath193 represent components parallel and perpendicular to @xmath194 , respectively . in a plane - symmetric case
, @xmath195 only appears .
hereafter we write this by @xmath196 .
if matter distribution is clustering , @xmath196 is expected to be negative .
giving initial data by random gaussian , the pdf of a pairwise peculiar velocity is gaussian at initial time . during evolution ,
the pdf will deviate from gaussian .
in fact , from fig.[fig : pdf - pv ] , which shows the pdf of a peculiar velocity in nonlinear regime , we find that it is not gaussian and approaches an exponential form in a small velocity region . as a reference , we show that the pdf for the za , which shows a gaussian form . the reason may be understood as follows : in a plane - symmetric model , gravitational potential is proportional to a distance between two sheets ( @xmath197 ) . then
, even if two sheets approaches very closely , a gravitational force does not become strong but keep constant .
when we take into account a backreaction effect , however , a gravitational force will be strong in a clustering region , because the expansion rate of a local domain is slow down .
the strengthening of a gravitational force in a cluster region may make a deviation of the pdf of pairwise peculiar velocity from its gaussian form .
note that in the 3d system , which shows non - gaussian pdf , a gravitational force increases when two particles approach .
another interesting observable variable is a deceleration parameter .
the recent observation of type ia supernova may suggest an acceleration of the universe@xcite .
although this result may naively suggest an existence of dark energy such as a cosmological constant @xmath199 , we could find some effective model without dark energy which explain the observation .
then we shall estimate a deceleration parameter averaged in a local domain here .
we define local deceleration parameter @xmath157 as @xmath200 which can be evaluated by eq .
( [ eqn : g - friedmann ] ) and @xmath201 .
buchert et al@xcite showed the evolution of deceleration parameter for the za .
they picked up overdense and underdense regions of three-@xmath21 fluctuations and found that @xmath157 for an underdense region could be a present day value , which is smaller by more than 200% than that of the background e - ds universe , although such a region is still decelerating .
although our approach includes a backreaction consistently , our analysis shows that a deviation of @xmath157 does not get so large .
the difference of @xmath157 from the za is very little even just before the shell crossing .
we show the time evolution of @xmath157 for a plane - symmetric 1-dimensional model in fig .
[ fig : qd - evolv ] . even if a domain is extremely underdense , the domain is decelerating .
this may be because our approach is still perturbative .
we will discuss it further in the next section .
we propose new lagrangian perturbation theory with a backreaction effect by inhomogeneity of density perturbations and present a set of basic equations .
the inhomogeneity affects the expansion rate in a local domain and its own growing rate . in a one - dimensional plane - symmetric model ,
we have numerically analyzed our approach , and calculated the growing rate density perturbations and the pdf of several observed variables .
we set our initial conditions as random gaussian distribution . from our analysis
, we show that the expansion rate of an overdense region is faster than that of the whole universe as expected .
we also show that the local hubble parameter may deviate from the global one by about 1.2 % for a @xmath184-grid domain
. it may be too small to distinguish its effect in the present observations@xcite , but it could become important in fully non - linear stage , which we can not describe in the present approach .
the pdf of density is slightly different from that of the za . in our model ,
an underdense region expands faster than that in the e - ds model and then its volume gets larger . hence , a probability for a negative @xmath1 region increases as seen in the pdf . as for a peculiar velocity ,
even if we take into account a backreaction , its pdf is still gaussian .
the pdf of pairwise peculiar velocity , however , shows an effective difference from the conventional lagrangian approach . in one - dimensional plane symmetric case ,
the pdf in the conventional lagrangian approximation ( the za ) is gaussian , but ours is not but approaches an exponential form in a small relative - velocity region , which agree with the n - body simulation and the 3d lagrangian approximation@xcite .
finally , we mention about recent observation about cosmological parameters . according to the observation of type ia supernova
, the expansion of the universe seems to accelerate@xcite . combining observation of the cosmic microwave background radiation ( cmbr ) ,
the result suggests existence of dark energy such as a cosmological constant @xmath199@xcite .
however , this produces another difficulty , that is the so - called cosmological constant problem . to avoid such a difficulty , if we could explain the observation without cosmological constant , it would be more natural .
recently , tomita discussed such possibility assuming we are in a large local void@xcite .
globally the universe is flat ( eds universe ) , but we are sitting near the center of a local void , which existence is observationally confirmed .
then he calculated the luminosity distance , finding that the observation can be explain by such a model .
we then wonder whether we could have the similar explanation if we are living in an effective void , which is a domain with an averaged energy density below the critical value . since we have to treat a strongly nonlinear structure to explain the observation , this is beyond our present approach .
however there is some indication . if we analyze the time dependence of the backreaction term @xmath46 , which evolves as @xmath202 in a linear perturbation level .
this time dependence is the same as a perfect fluid with the equation of state @xmath203 , which could be dark energy .
if we could explain the above observation without introduction of any strange matter but just by inhomogeneity of density distribution in the universe , we will have a natural understanding of the universe .
this is under investigation .
in addition , there are two further approaches which are related to the present work .
takada and futamase@xcite proposed that they divided lagrangian perturbation to large - scale and small - scale perturbations , then discussed interaction between those scales .
taruya and soda@xcite discussed dynamics of averaged variables in the case of a spherical infall model taking into account a backreaction effect .
since those are interesting approaches , it may be useful to use their approaches to discuss the present subject and compare those results with ours in future
.    we would like to thank t. buchert , m. morikawa , m. morita , and y. sota for useful discussions and comments .
this work was supported partially by the waseda university grant for special research projects .
j. e. peebles , _ the large scale structure of the universe _ ( princeton university press , princeton , new jersey , 1980 ) .    ya .
b. zeldovich , astron . astrophys . * 5 * , 84 ( 1970 ) .
t. buchert , astron . astrophys . * 223 * , 9 ( 1989 ) .
t. buchert , mon . not .
* 254 * , 729 ( 1992 ) .
p. catelan , mon . not .
. soc . * 276 * , 115 ( 1995 ) .
v. sahni and p. coles , phys . rep . *
262 * , 1 ( 1995 ) .
d. munshi , v. sahni , and a. a. starobinsky , astrophys .
j. * 436 * , 517 ( 1994 ) .
v. sahni and s. shandarin , mon . not .
. soc . * 282 * , 641 ( 1996 ) .
a. yoshisato , t. matsubara , and m. morikawa , astrophys .
j. * 498 * , 49 ( 1998 ) .    m. j. geller and j. p.
huchra , science * 246 * , 897 ( 1989 ) .
s. bildhauer and t. futamase , gen .
rel . grav .
* 23 * , 1251 ( 1991 ) .
t. futamase , phys .
d * 53 * , 681 ( 1996 )    h. russ , m. h. soffel , m. kasai , and g. brner , phys .
d * 56 * , 2044 ( 1997 ) .
w. r. storger , a. helmi , and d. f. torres , gr - qc/9904020 .
t. buchert , gen .
. grav . * 32 * , 105 ( 2000 ) .
t. buchert and j. ehlers , astron .
astrophys . * 320 * , 1 ( 1997 ) .
j. ehlers and t. buchert , gen .
* 29 * , 733 ( 1997 ) .
t. buchert , m. kerscher , and c. sicka , phys . rev .
d * 62 * , 043525 ( 2000 ) .
e. l. turner , r. cen , and j. p. ostriker , astron .
j. * 103 * , 1427 ( 1992 ) .
y. suto , t. suginohara , and y. inagaki , prog .
. phys . * 93 * , 839 ( 1995 ) .
t. t. nakamura and y. suto , astrophys .
j. * 447 * , l65 ( 1995 ) .
wu , z. deng , z. zou , l .- z .
fang , and b. qin , astrophys .
j. * 448 * , l65 ( 1995 ) .
x. shi and m. turner , astrophys .
j. * 493 * , 519 ( 1998 ) .
l. kofman , e. berschinger , j. m. gelb , a. nusser , and a. dekel , astrophys .
j. * 420 * , 44 ( 1994 ) .
t. padmanabhan and k. subramanian , astrophys .
j. * 410 * , 482 ( 1993 ) .
w. h. zurek , p. j. quinn , j. k. salmon , and m. s. warren , astrophys .
j. * 431 * , 559 ( 1994 ) .    n. seto and j. yokoyama , astrophys .
j. * 492 * , 421 ( 1998 ) .
w. l. freedman _
et al . _ , astrophys j. * 553 * , 47 ( 2001 ) .
s. perlmutter _
et al . _ ,
j. * 517 * , 565 ( 1999 ) .
a. h. jaffe _
et al . _ ,
lett . * 86 * , 3475 ( 2001 ) .
k. tomita , prog .
. phys . * 105 * , 419 ( 2001 ) .
k. tomita , mon.not.roy.astron.soc .
( to be published ) , [ astro - ph/0011484 ] .    k. tomita , astro - ph/0104141 .
m. takada and t. futamase , gen .
. grav . * 31 * , 461 ( 1999 ) .
a. taruya and j. soda , mon . not .
. soc . * 317 * , 873 ( 2000 ) .
using the lagrangian approximation , we estimate a backreaction term @xmath46 . using the lagrangian perturbation @xmath99 , we shall rewrite the backreaction term . for convenience ,
first we define the functional determinant of three functions @xmath204 : @xmath205 the invariants of the velocity gradient @xmath206 are given in terms of functional determinants as @xmath207 introducing the r.h.s .
of the lagrangian coordinates @xmath208 and perturbations @xmath209 by @xmath210 , the r.h.s . of these equations
are rewritten as @xmath211 \nonumber \\      & & + a^3\left[{\cal j}(\dot{s}_i , q_j , q_k)+2{\cal j }      ( \dot{s}_i , s_j , q_k ) + { \cal j}(\dot{s}_i , s_j , s_k)\right ] \ ; , \label{eqn : func_det1}\\ { \cal j}(\dot{r}_i , \dot{r}_j , r_k ) & = & a\dot{a}^2\left[{\cal j}(q_i , q_j , q_k)+3 { \cal j}(s_i , q_j , q_k)+3{\cal j}(s_i , s_j , q_k)+{\cal j}(s_i , s_j , s_k)\right ] \nonumber \\ & & + 2a^2\dot a\left[{\cal j}(\dot{s}_i , q_j , q_k)+2{\cal j}(\dot{s}_i , s_j , q_k ) + { \cal j}(\dot{s}_i , s_j , s_k)\right ] \nonumber \\ & & + a^3\left[{\cal j}(\dot{s}_i , \dot{s}_j , q_k)+{\cal j } ( \dot{s}_i , \dot{s}_j , s_k ) \right ] \;. \label{eqn : func_det2}\end{aligned}\ ] ] introducing a simplified description follows : @xmath212 we find that most terms in ( [ eq : a2])-([eq : a4 ] ) with ( [ eqn : func_det1 ] ) , ( [ eqn : func_det2 ] ) are given by the following quantities ; @xmath213 @xmath214 and @xmath215 are written as @xmath216 but the last term @xmath217 can not be written by any of @xmath218 , @xmath219 , @xmath220 and those time derivatives .
if we take into account the perturbations up to the second order , we find @xmath223 where we use the variables in eqs .
( [ psi_def ] ) and ( [ zeta_def ] ) .
note that the second order perturbation @xmath224 do not appear in eq .
( [ eqn : la_qd_2nd ] ) .
buchert et al estimated the backreaction term @xmath46 assuming za , i.e. @xmath117 and the lagrangian perturbation @xmath99 is separable as @xmath225 and found the similar form to ( [ eqn : la_qd]).@xcite