Instruction: Write an abstract for the following acadamic paper.

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

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

Answer: viable models of modified gravity which satisfy both local as well as cosmological tests are investigated . it is demonstrated that some versions of such highly non - linear models exhibit multiply de sitter universe solutions , which often appear in pairs , being one of them stable and the other unstable . it is explicitly shown that , for some values of the parameters , it is possible to find several de sitter spaces ( as a rule , numerically ) ; one of them may serve for the inflationary stage , while the other can be used for the description of the dark energy epoch . the numerical evolution of the effective equation of state parameter is also presented , showing that these models can be considered as natural candidates for the unification of early - time inflation with late - time acceleration through ds critical points . moreover , based on the de sitter solutions , multiply sds universes are constructed which might also appear at the ( pre-)inflationary stage . their thermodynamics are studied and free energies are compared . = 5000


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the precision test of the standard model ( sm ) at the lep gave a possibility not only to determine all the parameters and particle masses at the level of radiative corrections but also afforded an opportunity for searching for signals of new heavy particles beyond the energy scale of it . on the base of the lep2 experiments the low bounds on parameters of various models extending the sm have been estimated and the scale of new physics was obtained @xcite . although no new particles were discovered , a general believe is that the energy scale of new physics to be of order 1 tev , that may serve as a guide for experiments at the tevatron and lhc . in this situation , any information about new heavy particles obtained on the base of the present day data is desirable and important . numerous extended models include the @xmath0 gauge boson  massive neutral vector particle associated with the extra @xmath1 subgroup of an underlying group . searching for this particle as a virtual state is widely discussed in the literature ( see for references @xcite ) . in the content of searching for @xmath0 at the lhc and the ilc an essential information and prospects for future investigations are given in lectures @xcite . such aspects as the mass of @xmath0 , couplings to the sm particles , @xmath2 mixing and its influence in various processes and particles parameters , distinctions between different models are discussed in details . we shall turn to these papers in what follows . as concerned a searching for @xmath0 in the lep experiments and the experiments at tevatron @xcite , it was carried out mainly in a model - dependent way . a wide class of popular models has been investigated and low bounds on the mass @xmath3 were estimated ( see , @xcite ) . as it is occurred , the low masses are varying in a wide energy interval 400 - 1800 gev dependently on a specific model . these bounds are a little bit different in the lep and tevatron experiments . in this situation a model - independent analysis is of interest .    in the papers @xcite of the present authors a new approach for the model - independent search for @xmath0-boson was proposed which , in contrast to other model - independent searches , gives a possibility to pick out uniquely this virtual state and determine its characteristics . the corresponding observables have also been introduced and applied to analyze the lep2 experiment data . our consideration is based on two constituents : 1 ) the relations between the effective low - energy couplings derived from the renormalization group ( rg ) equation for fermion scattering amplitudes . we called them the rg relations . due to these relations , a number of unknown @xmath0 parameters entering the amplitudes of different scattering processes considerably decreases . 2 ) when these relations are accounted for , some kinematics properties of the amplitude become uniquely correlated with this virtual state and the @xmath0 signals exhibit themselves . the rg relations allow to introduce observables correlated uniquely with the @xmath0-boson . comparing the mean values of the observables with the necessary specific values , one could arrive at a conclusion about the @xmath0 existence . the confidence level ( cl ) of these values has been estimated and adduced in addition . without taking into consideration the rg relations the determination of @xmath0-boson requires a supplementary specification due to a larger number of different couplings contributing to the observables . similar situation takes place in the `` helicity model fits '' of lep collaborations @xcite when different virtual states contribute to each of the specific models ( aa , vv , and so on ) . therefore these fits had the goal to discover any signals of new physics independently of the particular states which may cause deviations from the sm . note that the lep collaborations saw no indications of new contact four fermion interactions in these fits .    in refs . @xcite the one - parametric observables were introduced and the signals ( hints in fact ) of the @xmath0 have been determined at the 1@xmath4 cl in the @xmath5 process , and at the 2@xmath4 cl in the bhabha process . the @xmath0 mass was estimated to be 11.2 tev . an increase in statistics could make these signals more pronounced and there is a good chance to discover this particle at the lhc .    in ref . @xcite the updated results of the one - parameter fit and the complete many - parametric fit of the lep2 data were performed with the goal to estimate a possible signal of the @xmath0-boson with accounting for the final data of the lep collaborations delphi and opal @xcite . usually , in a many - parametric fit the uncertainty of the result increases drastically because of extra parameters . on the contrary , in our approach due to the rg relations between the low - energy couplings there are only 2 - 3 independent parameters for the investigated leptonic scattering processes . as it was showed in ref . @xcite , an inevitable increase of confidence areas in the many - parametric space was compensated due to accounting for all accessible experimental information . therefore , the uncertainty of the many - parametric fit was estimated as the comparable with previous one - parametric fits in refs . @xcite . in this approach the combined data fit for all lepton processes is also possible .    from the results obtained on the searching for abelian @xmath0 within the lep experiment data set we conclude that it is insufficient for convincing discovery of this particle as the virtual state . in this situation it is reasonable to analyze the data by using the neural network approach which is able to make a realistic prognoses for the parameters of interest . this investigation was done within the two parametric global fit of the lep2 data on the bhabha scattering process . as the result of all these considerations we derive at the 2@xmath4 cl the characteristics of the @xmath0 ( the vector @xmath6 and axial - vector @xmath7 couplings of the @xmath0 with sm leptons and the @xmath8 mixing ) . the @xmath0 mass is also estimated . due to the universality of the @xmath7 we also derived the model independent estimate of the @xmath0 axial - vector couplings to quarks , @xmath9 . note that the hints for the @xmath0 have been determined in all the processes considered that increases the reliability of the signal . these results may serve as a good input into the future lhc and ilc experiments and used in various aspects . to underline the importance of them we mention that there are many tools at the lhc for the identification of @xmath0 . but many of them are only applicable if @xmath0 is relatively light . the knowledge of the @xmath0 couplings to sm fermions also have important consequences . the paper is organized as follows . in sect . 2 we give a necessary information about the description of @xmath0 at low energies . in sects . 3 - 5 we discuss the origin of the rg relations , their explicit forms for the case of the heavy @xmath0 and consequences of the relations for scattering processes investigated . in sect . 6 the cross sections and the observables to pick out uniquely the virtual @xmath0 in the @xmath10 processes are given . the fits of data are described and discussed . then in sect . 7 the same is present for the bhabha process @xmath11 . the one parametric and two parametric fits are discussed . in sect . 8 the analysis of this process is carried out by using the neuron network approach . the criteria for training the network are introduced which guarantee the 2@xmath4 cl deviations of the data from the model containing the sm with extra @xmath12 . the obtained parameters of the @xmath0 practically coincide with that of derived in the one parameter analysis . in this way we determine the characteristics of the @xmath0 coming from the lep experiments . in sect . 9 we discuss the role of the present model - independent analysis for the lhc experiments . the discussion and comparison with results of other approaches are given in sect . 10 . in the appendix we describe the two - mass - scale yukawa model and analyze in detail how the decoupling of the loop contributions due to heavy virtual states is realized when the mixing of fields is taken into consideration . this point is an essential element of the approach developed . let us adduce a necessary information about the abelian @xmath0-boson . this particle is predicted by a number of grand unification models . among them the @xmath13 and @xmath14 based models @xcite ( for instance , lr , @xmath15 and so on ) are often discussed in the literature . in all the models , the abelian @xmath0-boson is described by a low - energy @xmath16 gauge subgroup originated in some symmetry breaking pattern .    at low energies , the @xmath0-boson can manifest itself by means of the couplings to the sm fermions and scalars as a virtual intermediate state . moreover , the @xmath17-boson couplings are also modified due to a @xmath17@xmath0 mixing . in principle , arbitrary effective @xmath0 interactions to the sm fields could be considered at low energies . however , the couplings of non - renormalizable types have to be suppressed by heavy mass scales because of decoupling . therefore , significant signals beyond the sm can be inspired by the couplings of renormalizable types . such couplings can be derived by adding new @xmath16-terms to the electroweak covariant derivatives @xmath18 in the lagrangian @xcite ( review , @xcite ) @xmath19 where summation over all the sm left - handed fermion doublets , leptons and quarks , @xmath20 , and the right - handed singlets , @xmath21 , is understood . @xmath22 denotes the charge of @xmath23 in positron charge units , @xmath24 , and @xmath25 for leptons and 1/3 for quarks .    for general purposes we derive the rg relations for the @xmath0 beyond the sm with two light higgs doublets ( thdm ) @xcite . @xmath0 interactions with the scalar doublets can be parametrized in a model - independent way as follows , @xmath26 in these formulas , @xmath27 are the charges associated with the @xmath28 and the @xmath0 gauge groups , respectively , @xmath29 are the pauli matrices , @xmath30 is the generator corresponding to the gauge group of the @xmath0 boson , and @xmath31 is the @xmath32 hypercharge .    the yukawa lagrangian can be written in the form @xmath33 \right.\nonumber\\&&\left . + g_{f_{u , i } } [ \bar{f}_l \phi_i^c ( f_u)_r + ( \bar{f}_u)_r \phi^{c+}_i f_l]\right),\end{aligned}\ ] ] where @xmath34 is the charge conjugated scalar doublet .    the lagrangian ( [ lscal ] ) leads to the @xmath17@xmath0 mixing . the @xmath17@xmath0 mixing angle @xmath35 is determined by the coupling @xmath36 as follows @xmath37 where @xmath38 is the sm weinberg angle , and @xmath39 is the electromagnetic fine structure constant . although the mixing angle is a small quantity of order @xmath40 , it contributes to the @xmath17-boson exchange amplitude and can not be neglected at the lep energies .    in what follows we will also use the @xmath0 couplings to the vector and axial - vector fermion currents defined as @xmath41 the lagrangian ( [ lf ] ) leads to the following interactions between the fermions and the @xmath17 and @xmath0 mass eigenstates : @xmath42f , \nonumber\\ { \cal l}_{z'\bar{f}f}&=&\frac{1}{2 } z'_\mu\bar{f}\gamma^\mu\left [ ( v_f+\gamma^5 a_f)\cos\theta_0 -\right.\nonumber\\&&\quad\left . -(v^\mathrm{sm}_{fz}+\gamma^5 a^\mathrm{sm}_{fz})\sin\theta_0\right]f,\end{aligned}\ ] ] where @xmath23 is an arbitrary sm fermion state ; @xmath43 , @xmath44 are the sm couplings of the @xmath17-boson .    since the @xmath0 couplings enter the cross - section together with the inverse @xmath0 mass , it is convenient to introduce the dimensionless couplings @xmath45 which can be constrained by experiments . low energy parameters @xmath46 , @xmath47 , @xmath48 , @xmath49 must be fitted in experiments . in most investigations they were considered as independent ones . in a particular model , the couplings @xmath50 , @xmath51 take some specific values . in case when the model is unknown , these parameters remain potentially arbitrary numbers . however , this is not the case if one assumes that the underlying extended model is a renormalizable one . in the papers @xcite it was shown that these parameters are correlated due to renormalizability . we called them the rg relations . since this notion is a key - point of our consideration , we discuss it in detail . what is rg relation ? generally speaking , this is a correlation between low energy parameters of interactions of a heavy new particle with known light particles of the sm following from the requirement that full unknown yet theory extending sm is to be renormalizable . strictly speaking , rg relations are the consequence of two constituencies :    1 . rg equation for a scattering amplitude ; 2 .   decoupling theorem . the latter one describes the modification of both the rg operator @xmath52 and an amplitude at the energy threshold @xmath53 of new physics . here , @xmath54- and @xmath55-functions correspond to all the charges @xmath56 and fields and masses @xmath57 of the underlying theory . the rg equation for a scattering amplitude @xmath23 reads , @xmath58 where @xmath23 accounts for as intermediate states either the light or heavy virtual particles of the full theory . the standard usage of the rg equation is to improve the amplitude by solving this equation for the operator @xmath59 calculated in a given order of perturbation theory . however , to search for heavy virtual particles , we will use eq . ( [ rge ] ) in another way . first we note that for any renormalizable theory , the rg equation is just identity , if @xmath23 and @xmath59 are calculated in a given order of loop expansion . in this case ( [ rge ] ) expresses the well known fact that the structure of the divergent term coincides with the structure of the corresponding term in a tree - level lagrangian . for example , in massless qed , the tree - level plus one - loop one - particle - irreducible vertex function describing scattering of electron in an external electromagnetic field @xmath60 , @xmath61 , is        if we calculate the rg operator in one - loop order @xmath62 where @xmath63 , @xmath64 are the beta - function and the anomalous dimensions of electromagnetic and electron fields , respectively , and apply it to @xmath65 , we obtain @xmath66 then , accounting for the values of @xmath67 and the factor @xmath68 in @xmath69 , we observe that the first and the last terms in the r.h.s . cancel . since @xmath70-dependent term in @xmath71 is @xmath72 , we see that eq.([rge1 ] ) is identity in the order @xmath73 . next important point is that in a theory with different mass scales the decoupling of heavy - loop contributions at the threshold of heavy masses , @xmath53 , results in the following property : the running of all functions is regulated by the loops of light particles . therefore , the @xmath74 and @xmath75 functions at low energies are determined by the sm particles , only . this fact is the consequence of the decoupling theorem @xcite . the decoupling results in the redefinition of parameters at the scale @xmath53 and removing heavy - particle loop contributions from rg equation @xcite : @xmath76 where @xmath77 and @xmath78 denote the parameters of the sm . they are calculated assuming that no heavy particles are excited inside loops . the matching between both sets of parameters @xmath77 , @xmath78 and @xmath56 , @xmath57 is chosen at the normalization point @xmath79 , @xmath80    the differential operator @xmath81 in the rg equation is in fact unique ; the apparently different @xmath81 in both theories are the same !    note that if a theory with different mass scales is specified one can freely replace the parameters @xmath77 , @xmath78 and @xmath56 , @xmath57 by each other @xcite . an example of the derivation and the main features of the rg relations are shown in the appendix for a simple model with different mass scales .    if underlying theory is not specified , the set of @xmath56 , @xmath57 is unknown . the low energy theory consists of the sm plus the effective lagrangian generated by the interactions of light particles with virtual heavy particle states . the low energy parameters @xmath82 of these interactions are arbitrary numbers which must be constrained by experiments . by calculating the rg operator @xmath81 and the scattering amplitudes of light particles in this ` external field ' in a chosen order of loop expansion , it is possible to obtain the model - independent correlations between @xmath82 . these are just the re relations . let us derive the correlations between @xmath50 , @xmath83 appearing due to renormalizability of the underlying theory containing @xmath0 .    in our case , the rg invariance of the vertex leads to the equation @xmath84 where @xmath85 and @xmath86 are computed with taking into account the loops of light particles . now we derive the rg relations following from the one - loop consideration . in accordance with the previous sections , the one - loop rg equation for the vertex function reads @xmath87 where @xmath88 and @xmath89 are the tree - level and one - loop contributions to the fermion-@xmath0 vertex . @xmath90 is the one - loop level part of the rg operator . vertex at the one - loop level . ]        to calculate these functions , only the divergent parts of the one - loop vertices are to be calculated . the corresponding diagrams are shown in fig . [ f-2 ] . the fermion anomalous dimensions @xmath91 can be calculated by using the diagrams in fig . [ f-3 ] . then , eq.([rgez ] ) leads to algebraic equations for the parameters @xmath50 , @xmath92 and @xmath49 which have two sets of solutions @xcite : @xmath93 and @xmath94 here @xmath23 and @xmath95 are the partners of the @xmath96 fermion doublet ( @xmath97 and @xmath98 ) , @xmath99 is the third component of weak isospin . the first of these relations describes the @xmath0 boson analogous to the third component of the @xmath96 gauge field . the couplings to the right - handed singlet are absent . the second relation corresponds to the abelian @xmath0 . in this case the sm lagrangian appears to be invariant with respect to the @xmath16 group associated with the @xmath0 . the last relation in eq.([rgr2 ] ) ensures the @xmath100 eq.([lyukawa ] ) is to be invariant with respect to the @xmath16 transformations . introducing the @xmath0 couplings to the vector and axial - vector fermion currents ( [ av ] ) , the last line in eq . ( [ rgr2 ] ) yields @xmath101 the couplings of the abelian @xmath0 to the axial - vector fermion current have a universal absolute value proportional to the @xmath0 coupling to the scalar doublet . these relations are model independent . in particular , they hold in all the known models containing the abelian @xmath0 . the most discussed models are derived from the @xmath102 group ( the so called lr , @xmath103-@xmath104 models ) . the tree - level @xmath0 couplings to the sm fermions in the models are shown in table 1 . [ zpcoup ]    .the @xmath105 couplings to the sm fermions in the most discussed @xmath102-based models . [ cols="<,^,^,^,^ " , ]     the standard @xmath106-fit gives the following constraints for the coupling @xmath107 at the 68% cl : @xmath108 hence it is seen that the most precise data of delphi and opal collaborations give no signal of the chiral @xmath0 at the 1@xmath4 cl . the combined value also shows no signal at the 1@xmath4 cl . from the combined fit the 95% cl bound on the value of @xmath107 can be derived , @xmath109 . supposing the @xmath0 coupling constant @xmath110 to be of the order of the electroweak one , @xmath111 , the corresponding @xmath0 mass has to be larger than @xmath112 tev .      a complete two parametric fit of experimental data based directly on the differential cross - sections has been carried out in ref . @xcite . due to only two independent couplings this fit is efficient . in the fitting the available final data for the differential cross - sections of the bhabha process were used . the data set consists of 299 bins including the data of aleph at 130 - 183 gev , delphi at 189 - 207 gev , l3 at 183 - 189 gev , and opal at 130 - 207 gev @xcite . the fitting procedure is similar to that of discussed above for the abelian @xmath0 . the results can be summarized as follows . the parameter space of the chiral @xmath0 is the plane ( @xmath113 , @xmath114 ) . the minimum of the @xmath106-function , @xmath115 , is reached at zero value of @xmath113 ( @xmath116 ) and almost independent of the value of @xmath114 ( the maximal - likelihood values of the couplings ) . the 95% cl area ( @xmath117 ) is shown in figure 3 . [ fig - ch - many ]    as one can see , the zero point , @xmath118 ( the absence of the chiral @xmath0 boson ) is inside the confidence area . the value of @xmath106 in this point ( 238.62 ) is indistinguishable from the @xmath119 . in other words , the set of experimental data can not determine the signal of the chiral @xmath0-boson . as also is seen , the value of @xmath113 is constrained as @xmath120 at the 95% cl . this upper bound is in an agreement with the corresponding result of the one - parameter fit ( @xmath121 ) . thus , the @xmath0 mass has to be larger than @xmath122 tev , if the @xmath0 coupling constant @xmath110 is again supposed to be of the order of the electroweak one , @xmath111 . the fit of the differential cross - sections leads to a better accuracy for @xmath113 than the fit of the integrated cross - sections based on the same data . without accounting for the model - independent relations between the @xmath0 couplings it is impossible to obtain such results in this section we discuss all the assumptions giving a possibility to pick out the @xmath0 signal and determine its characteristics in a model independent way . we also note the role of the present results for the future lhc and ilc experiments . as it was already stressed , in searching for this particle at the lep and tevatron a model dependent analysis was mainly used . the motivation for this was the different number of chiral fermions involved in different models ( see , for example , @xcite ) . in this approach , the low bounds on @xmath3 have been estimated and a smallness of the @xmath2 mixing was also observed . on the contrary , in our model independent approach the rg relations between the parameters of the effective low energy lagrangians have been accounted for that gave a possibility to determine not only the bounds but also the mass and other parameters of the @xmath0 . first we note all the assumptions used in our considerations . we analyzed the four - fermion scattering amplitudes of order @xmath123 generated by the @xmath0 virtual states . the vertices linear in @xmath0 were included into the effective low - energy lagrangian . we also impose a number of natural conditions . the interactions of a renormalizable type are dominant at low energies @xmath124 . the non - renormalizable interactions generated at high energies due to radiation corrections are suppressed by the inverse heavy mass @xmath125 and neglected . we also assumed that the @xmath126 gauge group of the sm is a subgroup of the gut group . as a consequence , all structure constants connecting two sm gauge bosons with @xmath0 are to be zero . hence , the interactions of gauge fields of the types @xmath127 , and other are absent at the tree level . our effective lagrangian is also consistent with the absence of the tree - level flavor - changing neutral currents ( fcncs ) in the fermion sector . the renormalizable interactions of fermions and scalars are described by the yukawa lagrangian ( [ lyukawa ] ) that accounts for different possibilities of the yukawa sector without the tree - level fcncs . these assumptions are quite general and satisfied in a wide class of @xmath13 inspired models .    within these constraints for the low energy effective lagrangian the rg correlations have been derived . correspondingly , the model independent estimates of the mass @xmath3 and other parameters are regulated by the noted requirements . therefore , the extended underlying model has also to accept them .    in this regard , let us discuss the role of the obtained estimates for the lhc . as it is well known ( see , for example , @xcite ) , there are many tools at the lhc for @xmath0 identification . but many of them are only applicable if this particle is relatively light . our results are in favor to this case . next important point is the determination of @xmath0 couplings to the various sm fermions . as we have shown , the axial - vector couplings of the @xmath0 to the sm fermions are universal and proportional to its coupling to the higgs field . hence we have obtained an estimate of the @xmath128 couplings for both leptons and quarks . this is an essential input because experimental analyses for the lhc have mainly concentrated on being able to distinguish models and not on actual couplings . the vector coupling @xmath129 was also estimated that , in particular , may help to distinguish the decay of the @xmath0 resonance state to @xmath130 pairs . since the couplings @xmath131 and @xmath129 were estimated there is a possibility to distinguish this process from the decay of the @xmath132 system . in the literature @xcite on searching for the @xmath0 it is also mentioned that the determination of the @xmath0 couplings to fermions could be fulfilled channel by channel , @xmath133 ,  . in almost all these considerations the relations between the parameters have not been taken into account . but this is very essential for treating of experimental data and introducing relevant observables to measure .    other parameter is the @xmath2 mixing , which is responsible for different decay processes and the effective interaction vertices generated at the lhc @xcite . it is also determined by the axial - vector coupling ( see eq . ( [ grgav ] ) ) and estimated in a model independent way . remind that in our analysis ( in contrast to the approaches of the lep collaborations ) the mixing was systematically taken into account . its value is of the same order of magnitude as the parameters that were fitted in experiments . note also that the existence of other heavy particles with masses @xmath134 does not influence the rg relations which are the consequences of the necessary condition for renormalizability . in fact , this condition ( the structure of a divergence generated by radiation corrections coincides with that of the tree - level vertex ) holds for each renormalizable type interaction . an important role of the model independent results for searching for @xmath0 at the lhc and ilc consists , in particular , in possibility to determine the particle as a virtual state due to a large amount of relevant events . we mentioned already that , in principle , lep experiments were able to determine it if the statistics was sufficiently large . experiments at the ilc will increase numerously the data set of interest . in fact , the observables , introduced in sects . 6 and 7 for picking out uniquely @xmath135 and @xmath136 couplings in the leptonic scattering process , are also effective at energies @xmath137 gev and could be applied in future experiments at ilc . other model independent methods of searching for the @xmath0 as a resonance state are proposed in the literature ( see refs . we do not discuss them here because they take no relations between the parameters into consideration . besides , the main goal of the present paper is to adduce model independent information about the @xmath0 followed from experiments at low energies . different aspects of @xmath0 physics at the lhc are out of the scope of it . in this section we collect in a convenient form all the results obtained and make a comparison with other investigations on searching for @xmath0 at low energies . in fact , this is a large area to discuss . references to numerous results obtained in either model dependent or model independent approaches can be found in the surveys @xcite . further subdivision can be done into the considerations accounting for any type correlations between the parameters of the low energy effective interactions and that of assuming complete independence of them . because of a large amount of fitting parameters the latter are less predictable . now , for a convenience of readers let us present the results of fits of the @xmath0 parameters in terms of the popular notations @xcite . the lagrangian reads @xmath138f , \nonumber\\ { \cal l}_{z'\bar{f}f}&=&\frac{1}{2 } z'_\mu\bar{f}\gamma^\mu\left [ ( v'_f-\gamma^5 a'_f)\right]f,\end{aligned}\ ] ] with the sm values of the @xmath17 couplings @xmath139 where @xmath68 is the positron charge , @xmath22 is the fermion charge in the units of @xmath68 , @xmath140 for the neutrinos and @xmath141-type quarks , and @xmath142 for the charged leptons and @xmath143-type quarks . @xmath144 & @xmath145 , @xmath146 & @xmath147 , @xmath146 & @xmath148 , @xmath144 +   + @xmath149 & @xmath150 & - & @xmath151 & 0.437 +   + @xmath149 & - & @xmath152 & - & - + @xmath153 & @xmath154 & - & @xmath155 & 1.278 + @xmath156 & @xmath157 & - & @xmath158 & 0.464 +   + @xmath149 , @xmath159 & - & @xmath160 & - & - +     @xmath144 & @xmath145 , @xmath146 & @xmath147 , @xmath146 & @xmath148 , @xmath144 +   + @xmath149 & 68% & @xmath161 & - & @xmath162 & @xmath163 +   + @xmath149 & 95% & - & @xmath164 & - & - + @xmath153 & 95% & @xmath165 & - & @xmath166 & @xmath167 + @xmath153 , & & & & & + @xmath168 & 95% & @xmath169 & - & @xmath170 & @xmath171 +   + @xmath149 , & & & & & + @xmath153 , & 95% & @xmath172 & @xmath173 & @xmath174 & @xmath175 + @xmath168 & & & & & + @xmath149 , & & & & & + @xmath159 & 68% & @xmath176 & @xmath177 & @xmath178 & @xmath179 +    the results of the fits of the @xmath0 couplings to the sm leptons obtained from the analysis of the lep experiments are adduced in the tables [ fitmlv]-[fitcli ] and fig . [ fig : fitfigs ] . remind that due to the universality of the axial - vector coupling @xmath180 the same estimates take also place for quarks . first of all , one parameter fits of lep experiments as well as the many - parameter fit for the @xmath181 backward bins show the hints of the @xmath0 boson at the 1 - 2@xmath4 cl . due to this fact , the fits allow to determine the maximum likelihood values of @xmath0 parameters . in spite of uncertainties , these values can be used as a guiding line for the estimation of possible @xmath0 effects in the lhc experiments . the maximum likelihood values are given in table [ fitmlv ] . as is seen , different fits and processes lead to the comparable values of the @xmath0 parameters .    in table [ fitcli ] we present the confidence intervals for the fitted parameters . with this table one is able to estimate the uncertainty of the @xmath0 couplings as well as the lower bounds on the parameters . +   +    in fig . [ fig : fitfigs ] the maximum likelihood values and the cl intervals are drawn for the different values of the @xmath0 mass . all the plots exploit the same color scheme . the values excluded at 95% cl by the many - parameter fit of all the lep2 leptonic processes @xmath182 are shown in gray . the 95% confidence intervals from the one - parameter fit of lep2 @xmath183 are drawn in pink with the maximum likelihood values as the dashed red line . the corresponding results with taking into account the @xmath184 process only are shown in yellow with solid red line for the maximum likelihood values . the maximum likelihood values from the lep1 experiments are represented as dotted blue line . the 95% confidence interval from the one - parameter fit of the lep2 bhabha scattering is shown as the blue crosshatched region with the maximum likelihood values as the dashed blue line . finally , the 68% confidence interval and the maximum likelihood values from the many - parameter fit of the backward bins of lep2 bhabha scattering are shown in green . now we compare the above results with the ones of other fits accounting for the @xmath0 presence . as it was mentioned in introduction , lep collaborations have determined the model dependent low bounds on the @xmath0 mass which vary in a wide energy interval dependently on a model . the same has also been done for tevatron experiments . the modern low bound is @xmath185 gev . it is also well known that though almost all the present day data are described by the sm @xcite , the overall fit to the standard model is not very good . in ref . @xcite it was showed that the large difference in @xmath186 from the forward - backward asymmetry @xmath187 of the bottom quarks and the measurements from the slac sld experiment can be explained for physically reasonable higgs boson mass if one allows for one or more extra @xmath1 fields , that is @xmath0 . a specific model to describe @xmath0 physics of interest was proposed which introduces two type couplings to the hyper charge @xmath188 and to the baryon - minus - lepton number @xmath189 . within this model by using a number of precision measurements from lep1 , lep2 , sld and tevatron experiments the parameters @xmath190 and @xmath191 of the model were fitted . the presence of @xmath0 was not excluded at 68% cl . the value of @xmath190 was estimated to be of the same order of magnitude as in our analysis and is comparable with values of other parameters detected in the lep experiments . the erroneous claim that @xmath190 is two order less then the value derived from our table 4 is , probably , a consequence of some missed factors . the upper limit on the mass was also obtained @xmath192 tev . these two analyzes are different but complementary . a common feature of them is an accounting for the @xmath0 gauge boson as a necessary element of the data fits . the results are in accordance at 68 - 95% cl with the existence of the not heavy @xmath0 which has a good chance to be discovered at tevatron and/or lhc . in this appendix we are going to investigate the yukawa model with a heavy scalar field @xmath103 and a light scalar field @xmath193 @xcite . the goals of our investigation are two fold : 1 ) to derive the one - loop rg relation for the four - fermion scattering amplitude in the decoupling region and 2 ) to find out the possibility of reducing this relation in the equation for vertex describing the scattering of light particles on the external field when the mixing between heavy and light virtual states takes place . consider the four - fermion scattering @xmath196 . the @xmath197-matrix element at the one - loop level is given by @xmath198,\nonumber \\ s_\mathrm{1pr}&=&\sum\limits_{{\phi}_{1},{\phi}_{2}=\varphi , \chi}g_{{\phi}_{1 } } g_{{\phi}_{2}}\left ( \frac{{\delta}_{{\phi}_{1}{\phi}_{2}}}{s - m^2_{{\phi}_{1 } } } + \frac{\pi_{\phi_1\phi_2}(s)}{(s- m^2_{\phi_1})(s- m^2_{\phi_2})}\right )   \nonumber\\ & & \times { \bar\psi}\left(p_{4}\right)\left[1+\gamma\left ( p_{3 } , -p_{4}- p_{3}\right )   \right ] \psi\left(p_{3}\right)\nonumber \\ & & \times { \bar\psi}\left ( p_{1}\right ) \left[1+\gamma\left ( p_{2 } , -p_{1}- p_{2}\right )   \right ] \psi\left ( p_{2}\right),\end{aligned}\ ] ] where @xmath199 , @xmath200 is the contribution from the one - particle reducible diagrams shown in figs . [ fig : tree]-[fig : loop ] , and @xmath201 is the contribution from box diagrams . the one - loop polarization operator of scalar fields @xmath202 and the one - loop vertex function @xmath65 are related to the green functions as @xmath203 where @xmath204 is the spinor propagator in the momentum representation . using the dimensional regularization ( the dimension of the momentum space is @xmath206 ) and the @xmath207 renormalization scheme one can compute the renormalization constants @xmath208 , \nonumber\\ \delta\lambda^2&=&\frac{1}{4\pi^2\varepsilon } \left[\left(g_\chi^2 + 6\xi\right)\lambda^2 - 6g_\chi^2m^2-\rho m^2\right ] , \nonumber\\ z_g^{-1}&=&\left[1-\frac{3}{16\pi^2\varepsilon}\left(g_\varphi^2+g_\chi^2\right)\right ] { \left(z_\phi^{1/2}\right)}^\mathrm{t}.\end{aligned}\ ] ]    from eq . ( [ app5 ] ) we obtain the appropriate @xmath74 and @xmath75 functions at the one - loop level : @xmath209 then , the @xmath197-matrix element can be expressed in terms of the renormalized quantities . the contribution from the one - particle reducible diagrams becomes @xmath210\psi(p_3 ) \nonumber\\&&\times \bar\psi(p_1)\left[1+\gamma^\mathrm{(fin)}\left(p_2,-p_1-p_2\right)\right]\psi(p_2),\end{aligned}\ ] ] where the functions @xmath211 and @xmath212 are the expressions @xmath213 and @xmath65 without the terms proportional to @xmath214 . since the quantity @xmath201 is finite , the renormalization leaves it without changes .    introducing the rg operator at the one - loop level @xcite @xmath215 we determine that the following relation holds for the @xmath197-matrix element @xmath216 where the @xmath217 and the @xmath218 are the contributions to the @xmath200 at the tree level and at the one - loop level , respectively : @xmath219 @xmath220.\end{aligned}\ ] ]    the first term in eq . ( [ app11 ] ) is originated from the one - loop correction to the fermion - scalar vertex . the rest terms are connected with the polarization operator of scalars . the third term describes the one - loop mixing between the scalar fields . it is canceled in the rg relation ( [ app9 ] ) by the mass - dependent terms in the @xmath74 functions produced by the non - diagonal elements in @xmath221 . ( [ app9 ] ) is the consequence of the renormalizability of the model . it insures the leading logarithm terms of the one - loop @xmath197-matrix element to reproduce the appropriate tree - level structure . in contrast to the familiar treatment we are not going to improve scattering amplitudes by solving eq . ( [ app9 ] ) . we will use it as an algebraic identity implemented in the renormalizable theory . naturally if one knows the explicit couplings expressed in terms of the basic set of parameters of the model , this rg relation is trivially fulfilled . but the situation changes when the couplings are represented by unknown arbitrary parameters . in this case the rg relations are the algebraic equations dependent on these parameters and appropriate @xmath74 and @xmath75 functions . in the presence of a symmetry the number of @xmath74 and @xmath75 functions is less than the number of rg relations . so , one has non trivial system of equations relating the unknown couplings . for example , such a scenario is realized for the gauge coupling . although the considered simple model has no gauge couplings , we are able to demonstrate the general procedure of deriving the rg relations .      at energies @xmath222 the heavy scalar field @xmath103 is decoupled . this means , that the four - fermion scattering amplitude is described by the model with no heavy field @xmath103 plus terms of the order @xmath223 . at the tree level , this is the obvious consequence of the expansion of the heavy scalar propagator @xmath224,\ ] ] which is resulted in the effective contact four - fermion interaction in eq . ( [ app10 ] ) @xmath225 so , the tree level contribution to the scattering amplitude becomes @xmath226 \bar\psi\psi\times\bar\psi\psi,\ ] ] and the lowest order effects of the heavy scalar in the decoupling region are described by the parameter @xmath227 , only . decoupling of heavy particles is present also at the level of radiative corrections . the radiative corrections are generally described by various loop integrals in the momentum space ( the passarino  veltman functions ) . considering a passarino  veltman function with at least one heavy mass @xmath228 inside loop in the low - energy limit , one can see the following asymptotic behavior : the function splits into 1 ) possible energy - independent divergent part ( including also @xmath229 ) and 2 ) energy - dependent finite part which can be expanded by inverse powers of @xmath53 and vanishes at the small energies . the important property is that the @xmath229-term in the divergent part reproduces the logarithm of the cut - off scale . so , such a potentially large term has to be automatically absorbed by the renormalization at low energies and leads to no observable effects . however , if the renormalization is actually performed at high energies ( as in the @xmath207 renormalization scheme ) the potentially large @xmath229-terms should be re - summed manually by the redefinition of the physical couplings and masses at the scale @xmath53 . what is the form of the rg relations in the limit of large @xmath53 ? the method of constructing the rg equation in the decoupling region was proposed in @xcite . it introduces the redefinition of the parameters of the model allowing to remove all the heavy particle loop contributions to eq . ( [ app11 ] ) . let us define a new set of fields , charges and masses @xmath230 , @xmath231 , @xmath232 , @xmath233 , @xmath234 , @xmath235 @xmath236 where dots stand for the higher powers of @xmath229 responsible for the decoupling at higher loop orders . the differential operator ( [ app8 ] ) ban be rewritten in terms of these new low - energy parameters : @xmath237 where @xmath238 and @xmath239 functions are obtained from the one - loop relations ( [ app6 ] ) and ( [ app15 ] ) @xmath240 hence , one immediately notices that @xmath238 and @xmath239 functions contain only the light particle loop contributions , and all the heavy particle loop terms are completely removed from them . the @xmath197-matrix element expressed in terms of new parameters satisfies the following rg relation @xmath241 @xmath242 @xmath243 where @xmath244 is the redefined effective four - fermion coupling . as one can see , eq . ( [ app20 ] ) includes all the terms of eq . ( [ app11 ] ) except for the heavy particle loop contributions . it depends on the low energy quantities @xmath230 , @xmath231 , @xmath245 , @xmath246 , @xmath234 , @xmath235 . the first and the second terms in eq . ( [ app20 ] ) are just the one - loop amplitude calculated within the model with no heavy particles . the third and the fourth terms describe the light particle loop correction to the effective four - fermion coupling and the mixing of heavy and light virtual fields . due to the mixing term it is impossible to split the rg relation ( [ app18 ] ) for the @xmath197-matrix element into the one for vertices . hence , we are not able to consider eq . ( [ app18 ] ) in the framework of the scattering of light particles on an external field induced by the heavy virtual scalar . but this is an important step in deriving the rg relation for el parameters . fortunately , there is a simple procedure allowing to avoid the mixing in eq . ( [ app20 ] ) . the way is to incorporate the diagonalization of the leading logarithm terms of the scalar polarization operator into the redefinition of the @xmath247 , @xmath248 , @xmath231 , @xmath232 : @xmath249 \left(\zeta^{-1/2}\right)^\mathrm{t } \left(\begin{array}{c}\tilde{g}_\varphi \\ \tilde{g}_\chi \end{array}\right ) , \\ \nonumber\zeta^{1/2}&= & 1-\frac{\tilde{g}_\varphi\tilde{g}_\chi}{8\pi^2\left(\tilde\lambda^2-\tilde{m}^2\right ) } \log\frac{\kappa^2}{\tilde\lambda^2 } % \nonumber\\&&\times \left(\begin{array}{cc}0 & \tilde\lambda^2 - 6\tilde{m}^2\\-\tilde{m}^2 - 6\tilde{m}^2 & 0\end{array}\right).\end{aligned}\ ] ] the appropriate @xmath238 functions @xmath250 contain no terms connected with mixing between light and heavy scalars . so , the fourth term in eq . ( [ app20 ] ) is removed , and the rg relation for the @xmath197-matrix element becomes          at @xmath254 eq . ( [ app23 ] ) is just the rg identity for the scattering amplitude calculated in the absence of the heavy particles . the terms of order @xmath245 describe the rg relation for the effective low - energy four - fermion interaction in the decoupling region . the last one can be reduced in the rg relation for the vertex describing the scattering of the light particle ( fermion ) on the external field @xmath255 substituting the virtual heavy scalar : @xmath256 where @xmath257    eqs . ( [ app23])-([app27 ] ) is the main result of our investigation . one can derive them with only the knowledge about the low - energy couplings of heavy particle ( [ app13 ] ) and the lagrangian of the model with no heavy particles . one also has to ignore all the heavy particle loop contributions to the rg relation and the one - loop mixing between the heavy and the light fields . ( [ app23])-([app27 ] ) depend on the effective low - energy parameters , only . but as the difference between the original set of parameters and the low - energy one is of one - loop order , one may freely substitute them in eqs . ( [ app23])-([app26 ] ) . it is also possible to reduce the rg relation for scattering amplitudes to the one for vertex describing the scattering of light particles on the ` external ' field induced by the heavy virtual particle . in fact , this result is independent on the specific features of the considered model . 99 the lep collaborations aleph , delphi , l3 , opal , and the lep electroweak working group , hep - ex/0612034 . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j. * c33 * ( 2004 ) 173 [ hep - ex/0309053 ] ; eur . j. * c6 * ( 1999 ) 1 ; k.  ackerstaff _ et al . _ [ opal collaboration ] , eur . phys . j. * c2 * ( 1998 ) 441 . j.  abdallah _ et al . _ [ delphi collaboration ] , eur . phys . j. * c45 * ( 2006 ) 589 [ hep - ex/0512012 ] . a. leike , phys . rep . * 317 * ( 1999 ) 143 . p. langacker , arxiv:0801.1345 [ hep - ph ] . t. rizzo , hep - ph/0610104 . a. ferroglia , a. lorca , and j.j . van der bij , annalen phys * 16 * ( 2007 ) 563 - 578 [ hep - ph/0611174 ] . a. gulov and v. skalozub , eur . j. * c17 * ( 2000 ) 685 . a. gulov and v. skalozub , phys . rev . * d61 * ( 2000 ) 055007 . v. demchik , a. gulov , v. skalozub , and a. tishchenko , phys . * 67 * ( 2004 ) 1312 [ yad . fiz . * 67 * ( 2004 ) 1335 ] . a. gulov and v. skalozub , phys . rev . * d70 * ( 2004 ) 115010 . a. gulov and v. skalozub , phys . * d76 * ( 2007 ) 075008 . j.hewett and t.rizzo , phys . rep . * 183 * ( 1989 ) 193 . m. cvetic and b.w . lynn , phys . rev . * d35 * ( 1987 ) 51 . g.  degrassi and a.  sirlin , phys . rev . * d40 * ( 1989 ) 3066 . t. appelquist and j. carazzone , phys . rev . * d11 * ( 1975 ) 2856 . j. collins , f. wilczek , and a. zee , phys . rev . * d18 * ( 1978 ) 242 . m. bando , t. kugo , n. maekawa , and h. nakano , progress of theor . * 90 * ( 1993 ) 405 . t. rizzo , phys . * d55 * ( 1997 ) 5483 . r.  barate _ [ aleph collaboration ] , eur . j. * c12 * ( 2000 ) 183 [ hep - ex/9904011 ] . m.  acciarri _ et al . _ [ l3 collaboration ] , phys . * b479 * ( 2000 ) 101 [ hep - ex/0002034 ] . a. babich , g. della ricca , j. holt , p. osland , a. pankov , and n. paver , eur . j. * c29 * ( 2003 ) 103 . w. eadie , d. dryard , f. james , m. roos , and b. sadoulet , _ statistical methods in experimental physics _ , amsterdam , north - holland , 1971 . aleph collaboration , delphi collaboration , l3 collaboration , opal collaboration , sld collaboration , lep electroweak working group , sld electroweak group , and sld heavy flavour group , phys . * 427 * ( 2006 ) 257 [ hep - ex/0509008 ] . v. abazov _ [ d0 collaboration ] , phys . * b517 * ( 2001 ) 282 . a. buryk and v. skalozub , arxiv:0802.1486 [ hep - ph ] . schwind / mlpfit.html a. gulov and v. skalozub , phys . * 70 * ( 2007 ) 1100 - 1106 [ hep - ph/0510354 ] . m. dittmar , a .- s . nicollerat , and a. djouadi , phys . b583 * ( 2004 ) 111 - 120 [ hep - ph/0307020 ] . c. coriano , a. faraggi , and m. guzzi , phys . * d78 * ( 2008 ) 015012 [ arxiv:0802.1792 ] . f. petriello and s. quackenbush , phys . rev . * d77 * ( 2008 ) 115004 [ arxiv:0801.4389 ] . a. gulov and v. skalozub , phys . atom . nucl . * 63 * ( 2000 ) 139 - 143 [ yad . fiz . * 63 * ( 2000 ) 152 - 157 ] .
Abstract:

Answer: an approach to the model - independent searching for the @xmath0 gauge boson as a virtual state in scattering processes is developed . it accounts for as a basic requirement the renormalizability of underlying unspecified in other respects model . this results in a set of relations between low energy couplings of @xmath0 to fermions that reduces in an essential way the number of parameters to be fitted in experiments . on this ground the observables which uniquely pick out the @xmath0 boson in leptonic processes are introduced and the data of lep experiments analyzed . the @xmath0 couplings to leptons and quarks are estimated at 95% confidence level . these estimates may serve as a guide for experiments at the tevatron and/or lhc . a comparison with other approaches and results is given .