according to our current understanding @xcite , dark matter ( dm ) and dark energy ( de ) represent the two major components of the present universe .
surprisingly , it is found that the dm and de energy densities , @xmath0 and @xmath1 , are today roughly the same ( differing only by a factor of two ) , while their ratio has been varying by several orders of magnitude in the past history of the universe .
it seems quite natural , then , to explore the possibility of a dm  de interaction which could account for this coincidence .
this approach , however , is not free from problems if the de component is interpreted in terms of a dynamical quintessence @xcite scalar field .
indeed , such a scalar is constrained to be extremely light in order to fit the data , giving rise to unwanted long - range forces which may represent a severe threat to the equivalence principle .
in addition , couplings of the quintessence scalar with the gauge field strengths are potential sources of dangerous time variations of the fundamental constants .
a possible way of coping with a very light scalar while avoiding these shortcomings , is to choose to work in the framework of a scalar - tensor gravity ( st ) theory @xcite , where by construction matter has a purely metric coupling with gravity .
it has been shown @xcite that in this case , the scalar field can benefit from an attraction mechanism which , during the matter dominated era , makes st overlapping with standard general relativity ( gr ) . at the same time , st may possess a second `` attraction mechanism '' @xcite which will ensure the correct evolution of @xmath1 along a so - called ` tracking ' @xcite solution .    while st can very closely reproduce standard gr at the present time
, it may however lead to major differences in the past evolution of the universe , differences which may result in observable consequences for us today .
for example , it has been shown @xcite that st theories may have a profund impact on nucleosynthesis . at the same time , a curious fact has recently come to attention : a non  conventional dynamics of the quintessence scalar in the past history of the universe may remain ` hidden ' to the available cosmological observations , but manifest itself through the dm relic abundance @xcite .
it is then worth studying if st may provide a viable quintessence candidate and at the same time have an impact on the dm relic abundance .
we have considered this possibility and computed explicitly the differences from the standard case .
when considering st theories , the departures from standard cosmology are mainly due to the different expansion rate ( @xmath2 ) which they determine . such deviation from the usual expansion rate of gr bears potentially relevant consequences in all those phenomena which are closely dependent on the timing in which they occur .
the aim of the present work is to study the possible modifications of the expansion rate of the universe in st at the time of cold dark matter ( cdm ) freeze - out , focussing on weekly interacting massive particles ( wimps ) as the most natural candidates for cdm .
as it is well known , their present relic density depends on the precise moment they decouple and , in turn , on the precise moment the wimp annihilation rate equals the expansion rate of the universe .
we expect then that a variation of @xmath2 in the past may lead to measurable consequences on the wimp relic abundance .    in order to assess the allowed departure of @xmath2 from its standard value at wimps freeze - out
, we have to take into account the bounds imposed on st by phenomena at later epochs @xcite .
we already mentioned the crucial test of nucleosynthesis : passing the bbn exam will imply rather severe restrictions on the coupling of the scalar field of st to matter and also on its initial conditions at temperatures much higher than the wimp freeze - out .
coming to later epochs , we have to consider the photon decoupling and the consequent restrictions imposed by cmb data @xcite .
interestingly enough , we will show that the bbn `` filter '' on st is so efficient that it drastically limits any visible effect on the cmb spectrum ( in particular , shifts in the peaks positions are forced to be smaller than the experimental error ) .
more restrictive than the cmb probe turn out to be the gr tests , in particular after the tight bound on the post - newtonian parameter @xmath3 , recently provided by the cassini spacecraft @xcite .
this constraint becomes quite relevant when combined with that on the scalar equation of state of @xmath4 coming from sne ia data @xcite .
we will find that significant deviations from standard cosmology are possible only if the de equation of state differs appreciably from @xmath5 today .
in other words , if de is a cosmological constant , it is unlikely that future cosmological observations will be able to discriminate between st and gr .
the question we intend to explicitly tackle is the following : taking all the abovementioned restrictions ( bbn , cmb , gr tests ) into account , how much can the hubble parameter @xmath2 at the time of wimp freeze - out differ from its canonical value if st replaces gr ?
in other words , how much is the wimp relic density allowed to vary , if we consider st instead of gr ?    we find that in st theories the expansion rate of the universe at few gevs can profoudly differ from the usual value obtained in gr ( with variations up to five orders of magnitude ) and , yet , allow the correct light elements production at bbn .
this situation is perfectly analogous to the ` kination ' effect studied in @xcite , where a modification of @xmath2 at wimp freeze - out was induced by a short period of dominance of the scalar kinetic energy , although with some deal of fine - tuning . in the case considered here , the effect of st on @xmath2 depends on the strength of the scalar - matter coupling , however no particular fine - tuning is needed to pass the severe nucleosynthesis test even when large modifications of @xmath2 at freeze - out occur .
this means that the attraction of st towards gr proceeds very rapidly during the cooling of the universe from the few gevs of wimps freeze - out down to the mev range of nucleosynthesis .
the overlap of st with gr can subsequently be very efficient leading to st scenarios which can hardly be disentangled from ordinary gr in present tests at the post - newtonian level . the fact that st strongly affects the number of cdm particles may turn out to be the major signature of these theories .
> from the point of view of particle physics model building , these large variations in the wimps number density today is of utmost relevance
. particles which were not considered suitable to play a significant role in cdm scenarios can be rescued because of their enhanced number density . on the other hand , particles ( or regions of the parameter space for certain wimps candidates ) , which in usual gr scenarios constitute promising cdm candidates ,
would be excluded because their boosted number would overclose the universe .
these considerations become of particular interest if we focus on the case where the wimps correspond to the lightest supersymmetric particle .
a complete analysis of the cosmologically excluded and cosmologically interesting regions of the susy parameter spaces in different susy contexts , when st is considered , is presently in progress @xcite .
st theories represent a natural framework in which massless scalars may appear in the gravitational sector of the theory without being phenomenologically dangerous . in these theories
a metric coupling of matter with the scalar field is assumed , thus ensuring the equivalence principle and the constancy of all non - gravitational coupling constants @xcite . moreover , as discussed in @xcite , a large class of these models exhibit an attractor mechanism towards gr , that is , the expansion of the universe during the matter dominated era tends to drive the scalar fields toward a state where the theory becomes indistinguishable from gr .
st theories of gravity are defined by the action @xcite @xmath6 where & & s_g = d^4x .
[ jordan ] the matter fields @xmath7 are coupled only to the metric tensor @xmath8 and not to @xmath9 , _
i.e. _ @xmath10 $ ] .
@xmath11 is the ricci scalar constructed from the physical metric @xmath12 .
each st model is identified by the two functions @xmath13 and @xmath14 .
for instance , the well - known jordan - fierz - brans - dicke ( jfbd ) theory @xcite corresponds to @xmath15 ( constant ) and @xmath16 .
the matter energy - momentum tensor is conserved , masses and non - gravitational couplings are time independent , and in a locally inertial frame non gravitational physics laws take their usual form .
thus , the ` jordan ' frame variables @xmath12 and @xmath9 are also denoted as the ` physical ' ones in the literature . on the other hand , the equations of motion are rather cumbersome in this frame , as they mix spin-2 and spin-0 excitations . a more convenient formulation of the theory is obtained by defining two new gravitational field variables , @xmath17 and the dimensionless field @xmath18 , by means of the conformal transformation @xmath19 imposing the condition @xmath20 the gravitational action in the ` einstein frame ' reads @xmath21 , \ ] ] and matter couples to @xmath22 only through a purely metric coupling , @xmath23 \,\,\ , .\ ] ] in this frame masses and non - gravitational coupling constants are field - dependent , and the energy - momentum tensor of matter fields is not conserved separately , but only when summed with the scalar field one . on the other hand , the einstein frame planck mass @xmath24 is time - independent and the field equations have the simple form @xmath25 where @xmath26 \,,\ ] ] and @xmath27 is the matter energy - momentum tensor in the einstein frame .
the relevant point about the scalar field equation in ( [ eomgen ] ) is that its source is given by the trace of the matter energy - momentum tensor , @xmath28 , which implies the ( weak ) equivalence principle .
moreover , when @xmath29 the scalar field is decoupled from ordinary matter and the st theory is indistinguishable from ordinary gr .
we next consider an homogeneous cosmological space - time @xmath30 where the matter energy - momentum tensor admits the perfect - fluid representation @xmath31 with @xmath32 .
the friedmann - robertson - walker ( frw ) equations then take the form @xmath33 \\ \left ( \frac{\dot{a}}{a}\right ) ^{2}+\frac{k}{a^{2 } } & = & \frac{1}{3 \mp^2 } \left[\rho + \frac{\mp^2}{2}\dot\varphi^{2}+v(\varphi ) \right ] \\
\ddot{\varphi } + 3\frac{\dot{a}}{a}\dot\varphi & = &   -\frac{1}{m_{\ast}^2}\left [ \frac{\alpha ( \varphi ) } { \sqrt{2}}(\rho-3p)+   \frac{\partial v}{\partial \varphi }    \right ]   \ , \end{aligned}\ ] ] with the bianchi identity @xmath34    the physical proper time , scale factor , energy , and pressure , are related to their einstein frame counterparts by the relations @xmath35    defining new dimensionless variables @xmath36 and setting @xmath37 ( flat space ) the field equation of motion takes the more convenient form @xmath38 \varphi ^{\prime}= \nonumber \\ & & \quad\quad\quad\quad\quad - \sqrt{2}\ ; \alpha ( \varphi ) \;(1 - 3 w)-2 \;\lambda \;\frac{d\log v(\varphi ) } { d\varphi } , \label{eom}\end{aligned}\ ] ] where primes denote derivation with respect to @xmath39
. this will be our master equation .
the effect of the early presence of a scalar field on the physical processes will come through the jordan - frame hubble parameter @xmath40 : @xmath41 where @xmath42 is the einstein frame hubble parameter . in the flat
space case ( @xmath37 ) , eq  ( [ htilde ] ) finally gives : @xmath43 \,\,\,\ , .
\label{htilde2}\ ] ]
during radiation domination the scalar field energy density is suppressed , @xmath44 , and the first term in the rhs of eq .
( [ eom ] ) is proportional to & & 1 - 3w= = + & & = , [ rhs ] where the sum runs over all particles in thermal equilibrium , while @xmath45 is the contribution from the decoupled and pressureless matter abundance .    during radiation domination , @xmath46 , where @xmath47 is the jordan - frame temperature , and the contribution from a single particle in equilibrium gives @xmath48\,,\ ] ] with @xmath49 , @xmath50 the number of degrees of freedom of a , @xmath51 the number of relativistic degrees of freedom and @xmath52 \equiv \int_0^\infty dx \frac{x^2}{\varepsilon_a [ \exp(\varepsilon_a ) \pm 1]}\ , , \label{ffm}\ ] ] where @xmath53 and the minus ( plus ) sign in the denominator of the integrand holds for bosons ( fermions ) . in fig .
[ fm ] we plot @xmath54 $ ] . we see that it is different from zero only around @xmath55 , that is , for @xmath56 . for higher temperatures it is quadratically suppressed in @xmath57 , approaching the relativistic regime in which @xmath58 . for lower temperatures it is boltzmann- suppressed .
then , as emphasized in @xcite , the field @xmath59 evolves even during radiation domination , receiving a ` kick ' each time a particle in equilibrium becomes non - relativistic . the second term in eq .
( [ rhs ] ) is suppressed as @xmath60 , so it becomes relevant only as equivalence is approached .    in the following
, we will consider the evolution of the field from the freeze out of the dm particles down to today , so we will have to take into account all particles of masses between the freeze - out and matter - radiation equivalence .
= 3.2 in      during matter domination @xmath61 and , as long as the field energy density is subdominant ( @xmath62 ) , the rhs of the equation of motion ( [ eom ] ) is given by @xmath63 and the field evolution depends on the form of the coupling function @xmath64 .
as already mentioned , the jbd theory is given by a constant @xmath65 , and the value @xmath66 corresponds to gr .
a very attractive class of models is that in which the function @xmath64 has a zero with a positive slope , since this point , corresponding to gr , is an attractive fixed point for the field equation of motion @xcite .
it was emphasized in ref .
@xcite ( see also @xcite ) that the fixed point starts to be effective around matter - radiation equivalence , and that it governs the field evolution until recent epochs , when the quintessence potential becomes dominant .
if the latter has a run - away behavior , the same should be true for @xmath64 , so that the late - time behavior converges to gr .
the evolution of the field during the last redshifts depends on the nature of de .
we will consider two possibilities : a cosmological constant and a inverse - power law scalar potential for @xmath59 , which can be collectively represented by the potential @xmath67 @xmath68 corresponding to the cosmological constant .    in general , a cosmological constant in the einstein frame does not correspond to a cosmological constant in the jordan frame , as one can read from eq .
( [ transf ] ) .
however , present tests of gr ( see next section ) imply that at late times @xmath69 , so that the two frames are almost coincident and the expansion histories during the last few redshifts are practically indistinguishable .    for the purpose of this paper , that is the analysis of the impact of de on st cosmology , the situation in which the quintessence field is different from @xmath59 and decoupled from it would be basically the same as that with a cosmological constant , since in both cases the second term in the rhs of eq .
( [ eom ] ) vanishes .
the main feature of the potential in ( [ potential ] ) for @xmath70 is the existence of ` tracker ' solutions , which are attractors in field space @xcite . in the @xmath71 limit ,
the late - time behavior is completely determined by the two parameters @xmath72 and @xmath73 .
a non - vanishing @xmath65 would modify the behavior of the field today , hopefully keeping the desirable property of insensitivity to the initial conditions .    in this paper
, we will consider the following choice for @xmath74 , a()=1 + b e^- , [ aphi ] corresponding to ( ) = - , [ alphaphi ] which has a run - away behavior with positive slope , as required by the discussion at the end of the previous subsection .
the choice for the parameters @xmath75 and @xmath76 will be discussed in the following section .    in fig .
[ run ] we show the evolution of the background and field energy density .
we see that field energy densities corresponding to different initial conditions converge to the same solution , driven by @xmath77 .
notice also that the @xmath65-attractor becomes effective even before matter - radiation equivalence , due to the non - relativistic decoupling effect explained above .    in fig .
( [ levels ] ) we show the region of parameter plane @xmath78-@xmath73 giving @xmath79 , where @xmath80 and @xmath81 .
we see that in the st case , ( @xmath82 ) the region giving more negative values of the equation of state is somehow enlarged with respect to pure gr quintessence .
however , in the observationally allowed region for @xmath83 the influence of the @xmath75 parameter is negligible .
= 3.2 in    = 3.2 in
assuming @xmath84 in eq .
( [ htilde2 ] )  as we have checked numerically for the solutions relevant in this analysis  the jordan frame expansion rate during nucleosynthesis may be approximated as    ^2 a^2 ( ) .
the above expression should be compared to the gr one , in which the planck mass at nucleosynthesis was the same as today .
we obtain :    .|_nuc .
|_nuc = .
the change in the expansion rate is completely analogous to that obtained by adding extra neutrinos to the standard gr case . using the bound @xmath85 ( which is more restrictive than those obtained for instance in @xcite )
, we get < 1.08 .
[ bbn ]      at the post - newtonian level , the deviations from gr may be parametrized in terms of an effective field - dependent newtonian constant ) as well . ] @xmath86 and two dimensionless parameters @xmath87 and @xmath88 which , in the present theories turn out to be @xcite @xmath89 where @xmath90 .    a new constraint on the parameter
@xmath87 has been obtained recently using radio links with the cassini spacecraft @xcite , _
pn-1 = ( 2.12.3 ) 10 ^ -5 .
[ casbound ] present bounds on @xmath91 are @xmath92 and are less restrictive for our choice of @xmath77 , since @xmath93 .
the bound from the cassini spacecraft turns out to be quite strong when used in connection with the one on the equation of state @xmath94 from sne ia . in fig .
( [ betaw ] ) we show the excluded region in the @xmath76-@xmath94 plane implied by eq .
( [ casbound ] ) .
we see that an equation of state @xmath95 , as implied by sne ia data at @xmath96 c.l .
@xcite , requires either a large value for @xmath76 , or a very small @xmath75 . since the last case corresponds to an expansion history of the universe practically indistinguishable from gr , any non - standard behavior induced by the st theories in the past
should be likely accompanied by an equation of state different from @xmath5 today .
if de is a pure cosmological constant , then the bound from cassini implies @xmath97 ( making @xmath98 practically indistinguishable from one at least since bbn on ) , or unnaturally large values of @xmath76 .
= 3.2 in      the impact of a cosmological constant or quintessence on the cmb power spectrum has been extensively analyzed in refs .
@xcite . in the context of st theories ,
the problem has been studied in refs .
the main change with respect to a theory for de based on gr is due to a difference in the expansion rate , which affects the angular scale of the anisotropies .
the angle under which the first peak is seen goes as _
peak ~/l_peak ~v_s t_dec z_dec/ d(z_dec ) , where @xmath99 is the corresponding multipole , @xmath100 is the sound speed , @xmath101 and @xmath102 are the time and redshift of decoupling , and @xmath103 the distance to the last scattering surface .
the latter is given by d(z_dec ) = _ 1/(z_dec+1)^1 , [ dlast ] and is thus dominated by the behavior of @xmath2 close to the upper limit of integration , where @xmath104 is smaller .
for this reason , st theories passing the gr tests ( @xmath105 today , that is , @xmath106 ) imply a small deviation of the distance to the last scattering surface with respect to gr
.    on the other hand , the decoupling time might be significantly more perturbed .
it is given by an expression analogous to eq .
( [ dlast ] ) with the upper ( lower ) limit of integration replaced by @xmath107 ( @xmath108 ) . as a result ,
since the universe expanded faster than in gr at early times , we expect @xmath101 to be smaller , and the peak to move towards higher multiples .
we find , which is consistent with the numerical findings of ref .
@xcite .    due to the well known degeneracy of the cmb spectrum with respect to cosmological parameters , present measurements of the peak locations @xcite
do not translate straightforwardly into a bound on @xmath109 .
for instance , it is found that a shift in the peak multipole can be obtained also by varying the energy densities according to @xcite @xmath110 so that , in general , a full reanalysis of the cmb including the new parameter @xmath109 would be required .
however , in the present case we find that , once the bbn bound on @xmath111 has been imposed , the resulting values for @xmath109 are so close to unity as to give shifts in the peak multiples smaller than the experimental error .
thus , the cmb spectrum does not provide significant bounds to the present scenario .
having in mind all the bounds discussed in the previous section , we can now go on to compute the cosmological evolution of the scalar field and its impact on the dm relic abundance .    as a first step ,
we want to estimate if st can have a sizeable effect on the jordan - frame hubble parameter @xmath2 at the time of wimp decoupling , without violating any of the avaliable cosmological observations .
we will consider the function @xmath74 of eq .
( [ aphi ] ) , imposing on the parameters @xmath75 and @xmath76 the phenomenological constraints already discussed .
we will then compute the ratio @xmath112 at the decoupling time of a typical wimp of mass @xmath113  gev . in this way we will be able to get an estimate of the effect before going into further detail .
the tightest bound is that coming from eq .
( [ bbn ] ) .
it has an impact on both @xmath75 in eq .
( [ aphi ] ) and on the initial conditions of @xmath22 at temperatures higher than the wimp freeze - out .
indeed , since on the tracker solution the scalar field is @xmath114 today , it should have been @xmath115 at nucleosynthesis , otherwise it would not have reached the attractor in time @xcite .
this implies @xmath116 .
= 3.2 in    as already discussed , the equation for the dynamics of the scalar field @xmath59 is obtained by substituting the expression of eq .
( [ rhs ] ) in the rhs of eq .
( [ eom ] ) and choosing a coupling function @xmath64 as defined in eq .
( [ alphaphi ] ) . in the sum of eq .
( [ rhs ] ) only the terms corresponding to particles with @xmath117 have been considered , i.e. particles lighter than the critical temperature of the phase transition through which they acquire a mass ( see ref .
in particular we have taken into account the top quark , the @xmath118 , the @xmath119 , the bottom quark , the tau quark , the charmed quark , the pions , the muon , the electron and a wimp particle of mass @xmath113  gev .
numerical integration of the equation for @xmath59 has been carried out between approximately 500  gev and today .
we have then computed , through eq.([htilde2 ] ) , the modified expansion rate in the st theory at a temperature @xmath120  gev corresponding to a typical time of wimp decoupling and compared it to the expansion rate of the standard case at the same temperature .    in fig .
[ contorno ] we plot the ratio @xmath112 at @xmath120  gev as a function of the initial value of @xmath98 and initial ratio of the scalar to background energy density @xmath121 .
we have restricted the possible initial conditions to those regions of parameters values respecting the bbn bound of eq .
( [ bbn ] ) .
we see that we have been able to produce an enhancement of the expansion rate up to @xmath122 at the time of wimp decoupling .
it is then worth studying in more detail what happens to the wimp relic abundance .
let us now consider the calculation of the relic abundance of a dm wimp with mass @xmath123 and annihilation cross - section @xmath124 .
as already mentioned , laboratory clocks and rods measure the `` physical '' metric @xmath125 and so the standard laws of non - gravitational physics take their usual form in units of the interval @xmath126 . as outlined in ref.@xcite , the effect of the modified st gravity will enter the computation of particle physics processes ( like the wimp relic abundance ) through the `` physical '' expansion rate @xmath2 defined in eq .
( [ htilde ] ) .
we have therefore to implement the standard boltzmann equation with the modified physical hubble parameter @xmath2 : @xmath127 where @xmath128 , @xmath129 is the entropy density and @xmath130 is the wimp density per comoving volume .     a typical behaviour of the function @xmath74 defined in eq .  ( [ aphi ] ) , calculated for parameters @xmath131 and @xmath132 . ]
we have considered values of @xmath2 wich respect all the bounds discussed in section iv .
specifically , we have considered the function @xmath74 as given in eq .
( [ aphi ] ) with parameters @xmath131 and @xmath132 .
the function @xmath74 for this choice of parameters is plotted in fig .
[ fig : aphi ] , which shows that @xmath74 is very large at large temperatures , and then , at a temperature @xmath133 , sharply drops to values close to 1 before nucleosynthesis sets in .
a parametrization of the behaviour of @xmath74 for @xmath134 , that will be useful in the following discussion , is : @xmath135 where @xmath136 is the current temperature of the universe .
we have numerically checked that , in the regime we are considering , a good approximation to the physical hubble parameter is given by : @xmath137 the solution of the boltzmann equation is therefore formally the same as in the standard case , with the noticeable difference that now the hubble parameter gets an additional temperature dependence , given by the function @xmath74 .
this can be translated in a change in the effective number of degrees of freedom at temperature @xmath47 : @xmath138 an approximated solution of eq.([eq : boltzmann ] ) can be cast in a form analogous to the standard case : @xmath139 where @xmath140 and @xmath141 and @xmath142 are the wimp abundances per comovin volume today and at freeze  out , respectively .
numerical solution of the boltzmann equation eq .
( [ eq : boltzmann ] ) in a st cosmology for a toy  model of a dm wimp of mass @xmath143 gev and constant annihilation cross - section @xmath144 gev@xmath145 .
the temperature evolution of the wimp abundance @xmath146 clearly shows that freeze  out is anticipated , since the expansion rate of the universe is largely enhanced by the presence of the scalar field @xmath59 . at a value @xmath147
a re  annihilation phase occurs and @xmath146 drops to the present day value . ]
the freeze  out temperature is obtained by the following implicit equation : @xmath148 \label{eq : xf}\ ] ] where @xmath149 is the internal number of degrees of freedom of our wimp .
clearly , when @xmath150 we recover the standard case .
the relic abundance is then simply given by : @xmath151 where @xmath152 is the present entropy density and @xmath153 denotes the critical density .
a numerical solution of the boltzmann equation eq .
( [ eq : boltzmann ] ) in a st cosmology with the function @xmath154 given in fig .
[ fig : aphi ] is shown in fig .
[ fig : abundance ] for a toy  model of a dm wimp of mass @xmath143 gev and constant annihilation cross - section @xmath144 gev@xmath145 .
the temperature evolution of the wimp abundance @xmath146 clearly shows that freeze  out is anticipated , since the expansion rate of the universe is largely enhanced by the presence of the scalar field @xmath59 .
this effect is expected .
however , we note that a peculiar effect emerges : when the st theory approached gr ( a fact which is parametrized by @xmath155 at a temperature @xmath133 , which in our model is 0.1 gev ) , @xmath156 rapidly drops below the interaction rate @xmath157 establishing a short period during which the already frozen wimps are still abundant enough to start a sizeable re  annihilation .
this post - freeze  out `` re  annihilation phase '' has the effect of reducing the wimp abundance , which nevertheless remains much larger than in the standard case . for the specific case shown in fig .
[ fig : abundance ] the wimp relic abundance is @xmath158 for gr , while for a st cosmology becomes @xmath159 , with an increase of a factor of 44 .
the phenomenon of re  annihilation can be conveniently discussed in terms of the relation between the expansion rate of the universe @xmath156 and the wimp interaction rate @xmath160 . a numerical calculation of these two quantities is plotted in fig .
[ fig : rates ] as a function of the temperature .
the departure from equilibrium occurs earlier than in the gr case , because @xmath161 .
when decoupling is completed , the particles evolve with an approximately constant @xmath162 and @xmath163 , while the hubble rate evolves as @xmath164 , i.e. slower than in the standard case ( we have used here the approximate @xmath154 behavious of eq.([eq : aphi ] ) ) .    at the transition temperature
@xmath133 the hubble rate drops to its standard value @xmath165 and becomes smaller than the interaction rate : in this case the decoupled wimps start to annihilate again , for a short period . after this re  annihilation phase ,
the particles continue to evolve with an approximately constant abundace @xmath166 and @xmath157 recovers the behaviour @xmath167 , while @xmath168 as usual .
the expansion rate of the universe @xmath169 and the wimp interaction rate @xmath170 are plotted as a function of the temperature .
the re - annihilation effect discussed in the text is outlined .
the small drop in the rates at @xmath171 mev is due to the quark  hadron phase transition . ]
the ratio between the freeze  out values of @xmath172 in st cosmology and in gr as a function of the wimp mass .
the dashed , solid and dotted lines refer to @xmath173 gev@xmath145 , @xmath174 gev@xmath145 and @xmath175 gev@xmath145 , respectively . ]
we notice that a re  annihilation phase does not occur in the case of kination , i.e. in the case the energy density of the universe is dominated by the kinetic term of a scalar field @xcite . in this case the evolution of the expansion rate is @xmath176 during kination , and than evolves smootly into the standard behaviour @xmath177 . during kination
both @xmath2 and @xmath157 have the same @xmath47dependence and closely follow each other , until kination ends and the standard behaviour is recovered . re
annihilation is possibile if the phase during which the expansion rate has the transition toward its standard gr behaviour is faster than the post - freeze  out evolution of the interaction rate , i.e. faster than @xmath167 .    the change in the freeze  out temperature is shown in fig .
[ fig : freeze ] where we show the ratio between the freeze  out values of @xmath172 in st cosmology and in gr . the freeze  out temperature is anticipated about a factor of 2 , with a dependence also on the annihilation cross section , as is clear from eq .
( [ eq : xf ] ) : for very low values of @xmath178 the freeze - out temperature may be anticipated up to a factor of 5 .
for these low cross sections the relic abundance is anyway largely overabundant : we can therefore quantify the reduction in @xmath179 in a factor which ranges between 10% and 40% for wimps which can provide abundances in the cosmologically acceptable range .
increase in the wimp relic abundance in st cosmology with respect to the gr case .
the solid curve refers to an annihilation cross section constant in temperature , i.e. @xmath180 gev@xmath145 , while the dashed line stands for an annihilation cross section which evolves with temperature as @xmath181 gev@xmath182 . ]
the amount of increase in the relic abundance which is present in st cosmology is shown in fig .  [
fig : ratio ] .
the solid curve refers to an annihilation cross section constant in temperature , i.e. @xmath183 , while the dashed line stands for an annihilation cross section which evolves with temperature as : @xmath184 ( these two cases correspond to the two limiting situations of the usual non  relativistic expansion of the thermally averaged annihilation cross section : @xmath185 ) . in the case of @xmath186wave annihilation
the increase in relic abundance ranges from a factor of 10 up to a factor of 400 .
for a pure @xmath187 dependence , the enhancement can be as large as 3 orders of magnitude .
the behaviours shown in fig .
[ fig : ratio ] , which have been obtained by a numerical integration of the boltzmann equation eq .
( [ eq : boltzmann ] ) , can be understood by employing the approximate analytical solution ( [ eq : y0 ] ) .
in the case of @xmath183 , eq .  ( [ eq : y0 ] )
gives : @xmath188 in the standard gr case , and @xmath189 \nonumber \\ & + & { \cal g}\,m \left[g(x_\varphi)\frac{a}{x_\varphi}\right ] \label{eq : yphi}\end{aligned}\ ] ] in our st model where the @xmath154 function is given in eq .
( [ eq : aphi ] ) for @xmath134 and @xmath190 otherwise ( @xmath191 ) . for the sake of simplicity , in both solutions
we have dropped the term @xmath192 which adds a small correction , not relevant for the present approximate discussion . in both equations @xmath193 .
the ratio @xmath194 of the relic abundances is : @xmath195 where we have approximated @xmath196 and we have defined @xmath197 . by making explicit the mass dependencies we obtain : @xmath198 where the mass is expresses in gev , @xmath199 , @xmath200 , and the numerical values have been obtained for @xmath201 and @xmath202 ( since in our case @xmath133 is smaller than the quark  hadron phase transition which we have set at @xmath203 mev ) .
the analytic approximation of eq .  ( [ eq : ratio ] ) helps to explain the behaviour shown by the solid curve in fig .
[ fig : ratio ] , which has been obtained by numerical calculations which employ the exact form of the function @xmath74 . from eq .
( [ eq : ratio ] ) we can in fact derive that , for low masses , the ratio @xmath194 has the behaviour : @xmath204 which shows that in this mass regime @xmath194 grows almost linearly with the wimp mass @xmath123 , and it is larger for lower values of @xmath133 .
if we accept @xmath133 as low as the bbn scale , we can obtain a further increase in the relic abundance of a factor 100 on the top of the one showed in fig .
[ fig : ratioa ] for low values of @xmath123 .
when the wimp mass is very large , the ratio @xmath194 behaves as : @xmath205 with a slight drop with the mass .
the position of the maximum and the maximal value of @xmath194 are given by : @xmath206 and : @xmath207    these expressions show that the maximal effect is also obtained for the lowest values of @xmath133 ; in this case the position of @xmath208 is shifted toward lower masses . for @xmath133 at the bbn scale , the maximal increase in the wimp relic abundance is of the order of 3000 , instead of about 400 obtained for @xmath209 gev and shown in fig .
[ fig : ratioa ] .    an interesting property shown by eq .
( [ eq : ratio ] ) is that @xmath194 does not depend explicitely on the annihilation cross section @xmath183 , which drops out in the ratio .
an implicit dependence on the cross section is present through @xmath179 , as can be seen in eq .
( [ eq : xf ] ) .
this dependence however is only logarithmic and does not spoil the general behaviour of @xmath194 shown in fig .  [
fig : ratio ] .
this is shown in fig .
[ fig : ratioa ] : the largest difference occurs for very low annihilation cross sections , for which the deviation of @xmath179 is larger .
however , for cross sections of interest , i.e. cross section which provide relic abundances below the cosmologically acceptable upper bound , the values of @xmath194 are stable to a relatively good extent .
increase in the wimp relic abundance with an annihilation cross section @xmath183 , for different values of @xmath210 .
the dot  dashed , solid and dashed lines correspond to @xmath211 gev@xmath145 , @xmath174 gev@xmath145 and @xmath175 gev@xmath145 , respectively . ]
increase in the wimp relic abundance with an annihilation cross section @xmath212 , for different values of @xmath213 .
the dot  dashed , solid and dashed lines correspond to @xmath214 gev@xmath145 @xmath174 gev@xmath145 and @xmath215 gev@xmath145 , respectively . ]
the relic abundance in a st theory as a function of the wimp mass in the case of @xmath216 gev ( solid line ) and @xmath217 ( dashed line ) . the upper ( lower ) dotted lines corresponds to the gr case for @xmath218 gev and @xmath219 , respectively . ]
a similar analysis holds in the case of @xmath220 .
however , in this case the dependence of @xmath194 with @xmath179 is somehow stronger ( as obtained from the integration in eq .
( [ eq : y0 ] ) ) , and the effect of changing @xmath124 is slightely larger .
this effect can be seen in fig .
[ fig : ratiob ] , where @xmath194 is shown for different values of the parameter @xmath213 . notice that larger cross sections , which in the standard case provide lower values for the relic abundance , are the ones which get more enhanced in st cosmology .    finally ,
as an example we show in fig .
[ fig : omega ] the relic abundance as a function of the wimp mass in the case of @xmath221 gev and @xmath217 .
we see that , in this case , the relic abundance can be at the level required to explain the cdm content of the universe ( @xmath222 @xcite ) for a st theory , while it is underabundant in the standard case .
the models shown in fig .  [ fig : omega ] represent a case in which we can explain at the same time both the dm and de contents of the universe , and the interplay of the two component is crucial in determinig the right abundances of both dm and de .
an analysis of specific particle candidates of dm , in particular in supersymmetric models , will be examined elsewhere @xcite .
the idea of exploiting primordial ( ultralight ) scalars in order to shed some light on a dynamical interpretation of de is by now a widespread research topic in the literature . in this paper
we follow the promising proposal of considering the quintessence scalar as embedded in a scalar
tensor theory of gravity .
this approach is at variance with the usual interpretation of quintessence as a new light scalar whose interactions with matter are subject to the tight phenomenological constraints on the equivalence principle violation and time - variation of the fundamental coupling constants . identifying the quintessence field with the scalar component of a st theory , instead , does not pose any threat on the equivalence principle , since by construction matter has a purely metric coupling with gravity .
we focus on quintessence st models which possess a double `` attraction mechanism '' , one to gr and the other ensuring @xmath1 to follow a tracking solution .
these two simultaneous mechanisms act as a `` protection '' for the theory to prevent its fall into immediate troubles ( for instance , large departures from gr predictions ) .
nevertheless , we still obtain important phnomenological signatures which might disentangle this theory from gr or other alternative proposals .
the tests of our st scenario divide into two classes : deviations from gr and departures from standard cosmology , in particular concerning the expansion rate of the universe .
the latter effect has a big impact on the most distant epoch of the universe for which we have `` direct '' information , i.e. nucleosynthesis .
however , we pointed out that we can further extend the implications of a non - standard @xmath2 in the early universe to times prior to nucleosynthesis . sticking to the standard wimp picture of dm , one of the most relevant events
before bbn is the wimp decoupling which is expected to have occurred at a temperature of a few gevs .
our work shows that , despite the severe `` filters '' on st quintessence models which are provided by bbn , solar system tests of gravity and , to a lesser degree , by cmb , it is still possible to find remarkable enhancements on the expansion rate of the universe at wimp freeze - out , yielding to relic wimp abundances which can vary up to a few orders of magnitude with respect to the standard case .
in this paper we pointed out some general features of the new `` wimp story '' around its decoupling temperature in the presence of st quintessence .
in particular , we noticed that some unexpected effect can take place , such as a short phase of wimp `` re - annihilation '' when st approaches gr .
needless to say , such potentially ( very ) large deviations entail new prospects on the wimp characterization both for the choice of the cdm candidates and for their direct and indirect detection probes .
a thorough reconsideration of the `` traditional '' wimp identified with the lightest neutralino in susy extensions of the sm as well as the identification of other potentially viable cdm candidates in the st context is presently under way @xcite .
we acknowledge research grants funded by the italian ministero dellistruzione , delluniversit e della ricerca ( miur ) and by the istituto nazionale di fisica nucleare ( infn ) within the _ astroparticle physics project_. fr is partially supported by the university of padova fund for young researchers , research project n. @xmath223 .
s.  m.  carroll , phys .
* 81 * ( 1998 ) 3067 .
a.  masiero , m.  pietroni and f.  rosati , phys .
d * 61 * ( 2000 ) 023504 .
t.  damour , f.  piazza and g.  veneziano , phys .
d * 66 * ( 2002 ) 046007 ; t.  damour , f.  piazza and g.  veneziano , phys .
* 89 * ( 2002 ) 081601            a.  r.  liddle and r.  j.  scherrer , phys .
d * 59 * ( 1999 ) 023509 ; i.  zlatev , l.  m.  wang and p.  j.  steinhardt , phys .  rev .
* 82 * ( 1999 ) 896 ; p.  j.  steinhardt , l.  m.  wang and i.  zlatev , phys .  rev .
d * 59 * ( 1999 ) 123504 ; b.  ratra and p.  j.  e.  peebles , phys .  rev .
d * 37 * ( 1988 ) 3406 ; p.
j.  e.  peebles and b.  ratra , astrophys .
j.   * 325 * ( 1988 ) l17 .
a.  g.  riess _ et al . _
[ supernova search team collaboration ] , astron .
j.   * 116 * , 1009 ( 1998 ) ; s.  perlmutter _ et al .
_ [ supernova cosmology project collaboration ] , astrophys .
j.   * 517 * , 565 ( 1999 ) ; j.  l.  tonry _ et al .
_ , astrophys .
j.   * 594 * , 1 ( 2003 ) ; b.  j.  barris _ et al .
_ , astrophys .
j. * 602 * , 571 ( 2004 ) ; a.  g.  riess _ et al .
_ , astrophys .
j. * 607 * , 665 ( 2004 ) .
t.  damour , gr - qc/9606079 , lectures given at les houches 1992 , susy95 and corf 1995 . ;
g.  esposito - farese and d.  polarski , phys .
d * 63 * , 063504 ( 2001 ) ; b.  boisseau , g.  esposito - farese , d.  polarski and a.  a.  starobinsky , phys .
lett .   * 85 * , 2236 ( 2000 ) .
d.  huterer and m.  s.  turner , phys .
d * 64 * , 123527 ( 2001 ) ; s.  hannestad and e.  mortsell , phys .  rev .
d * 66 * , 063508 ( 2002 ) ; a.  balbi , c.  baccigalupi , f.  perrotta , s.  matarrese and n.  vittorio , astrophys .
j.   * 588 * , l5 ( 2003 ) .
x.  l.  chen and m.  kamionkowski , phys .
d * 60 * , 104036 ( 1999 ) ; f.  perrotta , c.  baccigalupi and s.  matarrese , phys .
d * 61 * , 023507 ( 2000 ) ; r.  nagata , t.  chiba and n.  sugiyama , phys .
d * 69 * , 083512 ( 2004 ) ; r.  nagata , t.  chiba and n.  sugiyama , phys .
d * 66 * , 103510 ( 2002 ) .