recently , there has been a lot of interest in cosmological solutions in the presence of a cosmological constant , when the latter is significant compared to the present total energy density of the universe .
indeed , the hubble diagram based on observations of type ia supernovae up to a redshift @xmath1 seems to imply that our universe is presently accelerating @xcite .
these data , when combined with the observed location of the first acoustic peak of the cmb temperature fluctuations , favor a spatially flat universe whose energy density is dominated by a `` cosmological constant''-like term .
the flatness of the universe is corroborated by the latest boomerang and maxima data @xcite , in accordance with the inflationary paradigm , though a marginally closed universe is still allowed by the position of the first acoustic ( doppler ) peak at @xmath2 .
a significant cosmological constant may help in resolving the dark matter problem
  for dustlike matter alone observations seem to imply @xmath3   and in reconciling flat cold dark matter ( cdm ) models with observations in the framework of @xmath4cdm models .
finally , a cosmological constant is an elegant way to allow a high hubble constant @xmath5 with @xmath6 and a sufficiently old universe @xmath7gyr @xcite ( see also , e.g. , @xcite for a recent comprehensive review and references therein ) .    therefore , this interpretation , if confirmed by future observations , constitutes a fundamental progress towards the solution of the dark matter problem and the formation of large - scale structure in the universe out of primordial fluctuations generated by some inflationary model .
that is certainly what makes it so appealing and gives it , maybe somehow prematurely , the status of new paradigm .
a striking consequence for our universe is then its present acceleration , for a large range of equations of state @xcite .
of course , from the point of view of particle physics , a pure cosmological constant of the order of magnitude @xmath8 , interpreted as the vacuum energy , is extremely problematic .
this is why attempts were made to find some alternative explanation to the origin of the acceleration under the form of some scalar field @xmath9 ( sometimes called quintessence @xcite , `` @xmath4''-field , etc . )
whose slowly varying energy density would mimic an effective cosmological constant .
this is very reminiscent of the mechanism producing the inflationary phase itself with the fundamental difference that this scalar field , which does not have to be _ a priori _ the inflaton , is accelerating the expansion today , therefore at a much lower energy scale .
this of course has problems of its own as this effective cosmological constant term started dominating the universe expansion only in the very recent past ( the so - called `` cosmic coincidence '' problem ) .
indeed , the energy density of the field @xmath9 must remain subdominant at very early stages and come to dominate in the recent past only .
hence , specific evolution properties are required to meet these constraints and were indeed shown to hold for particular potentials , partly alleviating the problem of the initial conditions . for inverse power - law potentials the energy density of the scalar field
was shown to decrease less rapidly than the background energy density so that it can be negligible in the early universe and still come to dominate in the recent past @xcite . for exponential potentials ,
@xcite the scalar field energy density has the very interesting behavior that it tends to a fixed fraction of the total energy density , these are the so - called `` tracker solutions '' .
hence a pure exponential potential is excluded if data confirm that the energy density of the scalar field is dominating today , as this fraction had to be small at the time of nucleosynthesis .
a slightly different potential is proposed in @xcite and a classification of the scaling behavior of the scalar field for various potentials has been given in @xcite .
hence , though a minimally coupled scalar field is an attractive possibility , some degree of fine tuning still remains in the parameters of the potential @xcite .
if one admits that it is some minimally coupled scalar field which plays the role of an effective cosmological constant while gravity is described by general relativity , the question immediately arises : what is the `` right '' potential @xmath10 of this scalar field ? in a recent work by starobinsky @xcite , the following `` phenomenological '' point of view was adopted : instead of looking for more or less well motivated models , like the interesting possibilities discussed above , it is perhaps more desirable to extract as much information as possible from the observations ( a similar approach can also be adopted to reconstruct the inflaton potential ) in order to reconstruct the scalar field potential , if the latter exists at all .
cosmological observations could then be used to constrain the particle physics model in which this scalar field is supposed to originate . in the context of general relativity plus a minimally coupled scalar field , it was shown that the reconstruction of @xmath10 can be implemented once the quantity @xmath11 , the luminosity distance as a function of redshift , is extracted from the observations @xcite , something that is expected in the near future .
can already be reconstructed from present experimental data , although not yet very accurately . ] the snap ( supernovae acceleration probe ) satellite will notably make measurements with an accuracy at the percent level up to @xmath12 .
of course , in this way only the recent past of our universe , up to redshifts @xmath13 ( for reference , we will push some of our simulations up to @xmath14 ) , is probed and so the reconstruction is made only for the corresponding part of the potential .
crucial information is therefore gained on the microscopic lagrangian of the theory through relatively `` low '' redshift cosmological observations .
a further step is to generalize the same mechanism in the framework of scalar - tensor theories of gravity , sometimes called `` generalized quintessence '' .
the usual minimally coupled models are certainly ruled out if , for example , it turns out that this component of the energy density obeys an equation of state @xmath15 with @xmath16 ( @xmath17 ) .
strangely enough , such an unexpected equation of state which in itself implies new physics , is in fair agreement with the observations @xcite .
also the inequality @xmath18 must hold for a minimally coupled scalar field , hence its violation would force us to consider more complicated theories , possibly scalar - tensor theories .
there are also strong theoretical motivations .
these theories , in which the scalar field participates in the gravitational interaction , are the most natural alternatives to general relativity ( gr ) .
indeed , scalar partners to the graviton generically arise in theoretical attempts at quantizing gravity or at unifying it with other interactions .
for instance , in superstrings theory , a dilaton is already present in the supermultiplet of the 10-dimensional graviton , and several other scalar fields ( called the moduli ) also appear when performing a kaluza - klein dimensional reduction to our usual spacetime .
moreover , contrary to other alternative theories of gravity , scalar - tensor theories respect most of gr s symmetries : conservation laws , constancy of ( non - gravitational ) constants , local lorentz invariance ( even if a subsystem is influenced by external masses ) , and they also have the capability of satisfying the weak equivalence principle ( universality of free fall of laboratory - size objects ) even for a strictly massless scalar field .
nevertheless , they can describe many possible deviations from gr , and their predictions have been thoroughly studied in various situations : solar - system experiments @xcite , binary - pulsar tests @xcite , gravitational - wave detection @xcite .
finally these scalar - tensor theories could play a crucial role in the very early universe , for example in the pre big bang inflationary model ( see e.g. @xcite ) .
thus , in this work we are investigating the possibility to have an accelerating universe in the context of scalar - tensor theories of gravity instead of pure gr .
this has indeed attracted a lot of interest recently and such cosmological models have been studied and possibly confronted with observations like cmb anisotropies , or the growth of energy density perturbations ( see for instance @xcite ) .
however , we emphasize once more that the central point of view adopted here , in analogy with starobinsky @xcite , is to constrain the model with the _ experimental _ knowledge of the hubble diagram up to @xmath13 .
this is precisely why use of the redshift @xmath19 as basic variable is crucial for our purpose : quantities like @xmath20 are directly observable , in contrast to , say , can be obtained from the knowledge of @xmath20 thanks to the relation @xmath21 $ ] , but the directly observable quantity is @xmath20 . ] @xmath22 or @xmath23 .
for instance , we have access to @xmath20 through the direct measurement of the luminosity distance in function of redshift @xmath11 .
in a recent letter @xcite , it was shown that the knowledge of both @xmath20 and @xmath24 is sufficient to reconstruct the full theory ( again , in the range probed by the data ) .
this means that we do _ not _ choose any specific theory _ a priori _ , but instead we reconstruct whatever theory possibly realized in nature .    as we will see , the knowledge of @xmath20 on its own , though insufficient in order to fully reconstruct a scalar - tensor theory unless one makes additional assumptions , turns out to be already very constraining when subclasses of models are considered .
this is particularly interesting because it means that cosmological observations at low redshifts implying an accelerated expansion might well give new constraints on scalar - tensor theories .
we will show that this is indeed the case .
throughout the paper , we use natural units for which @xmath25 , and the signature @xmath26 , together with the sign conventions of @xcite . in section
ii , we introduce the general formalism of scalar - tensor theories of gravity and their different parametrizations . in section iii , we briefly review the severe experimental restrictions imposed on these theories _ today_. in section iv , we consider frw universes in the framework of scalar - tensor gravity and we give the equations for the different parameterizations . in section
v , we review the full reconstruction problem . in section
vi , we give a detailed study of subclasses of models , which are investigated using the background equations . finally , in section vii , our results are summarized and discussed .
we are interested in a universe where gravity is described by a scalar - tensor theory , and we consider the action @xcite @xmath27\ . \label{s_jf}\ ] ] here , @xmath28 denotes the bare gravitational coupling constant ( which differs from the measured one , see eq .
( [ geff ] ) below ) , @xmath29 is the scalar curvature of @xmath30 , and @xmath31 its determinant . in ref .
@xcite , we used different conventions ( corresponding to the choice @xmath32 in the above action ) ; here , the quantity @xmath33 is dimensionless .
this factor @xmath33 needs to be positive for the gravitons to carry positive energy .
the action of matter @xmath34 is a functional of some matter fields @xmath35 and of the metric @xmath30 , but it does not involve the scalar field @xmath9 .
this ensures that the weak equivalence principle is exactly satisfied .
the dynamics of the real scalar field @xmath9 depends _ a priori _ on three functions : @xmath33 , @xmath36 , and the potential @xmath10 . however , one can always simplify @xmath36 by a redefinition of the scalar field , so that @xmath33 and @xmath36 can be reduced to only one unknown function .
two natural parametrizations are used in the literature : ( i )  the brans - dicke one , corresponding to @xmath37 and @xmath38 ; and ( ii )  the simple choice @xmath39 and @xmath33 arbitrary .
this second parametrization is however sometimes pathological .
[ the derivatives of @xmath9 can become imaginary in perfectly regular situations ; see the discussion about eq .
( [ positiveenergy ] ) below . ] in the following ,
we will write the field equations in terms of the two functions @xmath33 and @xmath36 , so that any particular choice can be recovered easily .    the variation of action ( [ s_jf ] ) gives straightforwardly    @xmath40    [ 2.2 ]    where @xmath41 is the trace of the matter energy - momentum tensor @xmath42 .
the scalar - field equation ( [ boxphi ] ) can of course be rewritten differently if one uses the trace of eq .
( [ einstein ] ) to replace the curvature scalar @xmath29 by its source , and one gets the brans - dicke - like equation @xmath43 where @xmath44 . [ in the brans - dicke representation where @xmath45 and @xmath46 , this factor @xmath47 reduces to the well - known expression @xmath48 . ] in the following , we will however use the form ( [ boxphi ] ) , which will simplify considerably our calculations .    the above equations are written in the so - called jordan frame ( jf ) . since in action
( [ s_jf ] ) , matter is universally coupled to @xmath30 , this `` jordan metric '' defines the lengths and times actually measured by laboratory rods and clocks ( which are made of matter ) .
all experimental data will thus have their usual interpretation in this frame .
in particular , the observed hubble parameter @xmath49 and the measured redshifts @xmath19 of distant objects are jordan - frame quantities .    however , it is usually much clearer to analyze the equations and the mathematical consistency of the solutions in the so - called einstein frame ( ef ) , defined by diagonalizing the kinetic terms of the graviton and the scalar field
. this is achieved thanks to a conformal transformation of the metric and a redefinition of the scalar field .
let us call @xmath50 and @xmath51 the new variables , and define    @xmath52    [ 2.4 ]    action ( [ s_jf ] ) then takes the form @xmath53\ , \label{s_ef}\ ] ] where @xmath54 is the determinant of @xmath50 , @xmath55 its inverse , and @xmath56 its scalar curvature . note that the first term looks like the action of general relativity , but that matter is now explicitly coupled to the scalar field @xmath51 through the conformal factor @xmath57 .
quantities referring to the einstein frame will always have an asterisk ( either in superscript or in subscript ) , e.g. @xmath58 and @xmath59 for the covariant derivative and the dalembertian with respect to the einstein metric .
the indices of einstein - frame tensors will also be lowered and raised with the einstein metric @xmath50 and its inverse @xmath55 .
the field equations deriving from action ( [ s_ef ] ) take the simple form    @xmath60    [ 2.6 ]    where @xmath61 is the coupling strength of the scalar field to matter sources @xcite , and @xmath62 is the trace of the matter energy - momentum tensor @xmath63 in einstein - frame units . from its definition
, one can deduce the relation @xmath64 with its jordan - frame counterpart .
let us underline that the cauchy problem is well posed in the einstein frame @xcite , because all the second - order derivatives of the fields are separated in the left - hand sides of eqs .
( [ 2.6 ] ) , whereas they are mixed in the jf equations ( [ 2.2 ] ) .
action ( [ s_ef ] ) also shows that the helicity-2 degree of freedom is described by the fluctuations of the einstein metric @xmath50 ( whose kinetic term is the standard einstein - hilbert one ) , and that the ef scalar @xmath51 is the true helicity-0 degree of freedom of the theory ( since its kinetic term has the standard form ) . on the other hand
, the fluctuations of the jordan metric @xmath30 actually describe a _ mixing _ of helicity-2 and helicity-0 excitations , and the jf scalar @xmath9 is related to the helicity-0 degree of freedom via the complicated relation ( [ varphi ] ) , because its kinetic term in action ( [ s_jf ] ) comes not only from the naive contribution @xmath65 but also from the cross term @xmath66 . in conclusion
, the theory can be mathematically well defined only if it is possible to write the ef action ( [ s_ef ] ) , notably with its negative sign for the scalar - field kinetic term ( so that @xmath51 carries positive energy ) .
if it happens that the transformation ( [ 2.4 ] ) is singular for particular values of @xmath9 , the consistency of the theory should be analyzed in the ef .
some singularities may be artifacts of the parametrization which is chosen to write action ( [ s_jf ] ) , and may not have any physical significance .
on the other hand , jordan - frame quantities may look sometimes regular while there is an actual singularity in the einstein frame ( a typical example is provided when @xmath33 vanishes ) . in this case , the solution should be considered as mathematically inconsistent . in the following
, we will see that the jf is better suited than the ef for our cosmological study , but we will always check the consistency of our results by finally translating them in terms of einstein - frame quantities .
the predictions of general relativity in weak - field conditions , and at present , are confirmed by solar - system experiments at the @xmath67 level @xcite .
one should therefore verify that the scalar - tensor models we are considering are presently close enough to einstein s theory .    if the scalar field is very massive
( say , if @xmath68 is large with respect to the inverse of the astronomical unit ) , its influence is exponentially small in solar - system experiments , even if it is strongly coupled to matter .
this situation corresponds to the particular scalar - tensor model considered in ref .
@xcite ( namely @xmath37 and @xmath69 in action ( [ s_jf ] ) , but assuming a large enough value for @xmath70 ) .
although this situation is phenomenologically acceptable , it remains somewhat problematic from a field theoretical viewpoint , since the massive scalar would _ a priori _ desintegrate into lighter ( matter ) particles .    on the contrary ,
if the scalar mass is small with respect to the inverse solar - system distances , it must be presently very weakly coupled to matter for the theory to be consistent with experimental data . at the first post - newtonian order ( @xmath71 with respect to the newtonian interaction )
, the deviations from general relativity can be parametrized by two real numbers , that eddington @xcite denoted as @xmath72 and @xmath73 . in the present framework
, they take the form @xcite
@xmath74    [ 3.1 ]    where the first expressions are given in terms of the einstein - frame notation ( [ s_ef])-([alpha ] ) , whereas the last ones correspond to the jordan - frame general representation ( [ s_jf ] ) . to simplify , the second expression of eq .
( [ betappn ] ) has been written in terms of the derivative of ( [ gammappn ] ) with respect to @xmath9 .
using the upper bounds on @xmath73 from solar - system measurements  @xcite , we thus get the constraint @xmath75 where an index @xmath76 means the present value of the corresponding quantity . on the other hand ,
the experimental bounds on @xmath72 can not be used to constrain the derivative @xmath77 appearing in eq .
( [ betappn ] ) , since it is multiplied by a factor @xmath78 consistent with 0 . because of nonperturbative strong - field effects , binary - pulsar tests are however directly sensitive to this derivative , _
i.e. _ , to the ratio @xmath79 . in a generic class of scalar - tensor models , refs .
@xcite have obtained the bound @xmath80    from action eq .
( [ s_jf ] ) , one can naively define newton s gravitational constant as the inverse factor of the curvature scalar @xmath29 : @xmath81 however , @xmath82 does not have the same physical meaning as newton s gravitational constant in gr .
indeed , the actual newtonian force measured ( in cavendish - type experiments ) between two close test masses @xmath83 and @xmath84 is of the form @xmath85 , where the effective gravitational constant reads @xcite @xmath86 the contribution @xmath87 is due to the exchange of a graviton between the two bodies , whereas @xmath88 comes from the exchange of a scalar particle between them .
of course , when the distance between the bodies becomes larger than the inverse mass of the scalar field , its influence becomes negligible and one gets @xmath89 .
note that as usual , the last expression in eq .
( [ geff ] ) , in terms of jordan - frame notation , is much more complicated than its einstein - frame counterpart . in the particular brans - dicke representation , @xmath45 and @xmath46
, it however reduces to the simpler ( and well - known ) form @xmath90 .
the experimental bound ( [ df0 ] ) shows that the present values of @xmath91 and @xmath82 differ by less than @xmath92 .
however , they can _ a priori _ differ significantly in the past .
it should be noted that the experimental limit on the time variation of the gravitational constant , @xmath93  @xcite , does _ not _ imply any constraint on @xmath94 .
indeed , @xmath91 can be almost constant even if @xmath95 ( or @xmath96 ) varies significantly .
a simple example is provided by barker s theory @xcite , in which @xmath97  : one gets @xmath98 , which is strictly constant independently of the time variations of @xmath99 .
nevertheless , as pointed out in @xcite , under reasonable cosmological assumptions , one can derive @xmath100 with @xmath101 accuracy up to redshifts @xmath102 .
the equations derived in this section generalize those of our previous paper @xcite in several ways .
first , we use the most general representation ( [ s_jf ] ) of the theory , instead of the simpler choice @xmath103 that was made in @xcite .
second , we take into account a possible spatial curvature of the universe , which will be an interesting possibility in our studies of sec .
third , we write the equations for an arbitrary pressure of the perfect fluid describing matter in the universe
. this will not be useful for our reconstruction program of the following sections , as matter can be assumed to be simply dustlike for the redshifts @xmath104 that we will consider , but these general equations may be interesting for further cosmological studies of earlier epochs of the universe . finally , we comment on the einstein - frame version of these equations , which are mathematically simpler , but actually more difficult to use for our purpose .
we consider a friedmann - robertson - walker ( frw ) universe whose background metric in the jordan frame is given by    @xmath105    [ 4.1 ]    where @xmath106 , @xmath76 , or @xmath107 for spatially open , flat , or closed universes respectively . the scalar field @xmath9 ( or @xmath51 , in the ef ) is also assumed to depend only on time .
since the relation between the ef and jf is given by @xmath108 , see eqs .
( [ 2.4 ] ) , our universe is still of the frw type in the ef , with @xmath109 and @xmath110 in the following , matter will be described by a perfect fluid , and we will write its energy - momentum tensor as @xmath111 where @xmath112 and @xmath113 are the spacetime components of the four - dimensional unit velocity of matter , in jf and ef units respectively . as we are interested in a frw background , the spatial components @xmath114 and @xmath115 ( @xmath116 ) all vanish . from ( [ t ] )
, we deduce the relation between the matter density and pressure in both frames : @xmath117    the background equations in the jf follow from ( [ einstein])([matter ] ) , and read    @xmath118    [ 4.5 ]    where @xmath119 , and a dot denotes differentiation with respect to the jordan - frame time @xmath120 . as usual , if @xmath121 , eq .
( [ dmatter ] ) is trivially integrated as @xmath122 ( and in particular @xmath123 for dustlike matter ) .
equation  ( [ ddotphi ] ) is actually a consequence of the other three , and we will not need it in the following .    since these equations correspond to the most general parametrization ( [ s_jf ] ) of scalar - tensor theories , many particular cases are easily recovered .
for instance , the case of a minimally coupled scalar field @xcite is obtained for constant values of @xmath96 and @xmath124 ( say , @xmath125 and @xmath126 ) , and the particular model considered in @xcite is recovered immediately for @xmath45 and @xmath127 .
the corresponding background equations in the ef are very similar to those in general relativity .
they follow from eq .
( [ einstein * ] ) , and read    @xmath128    [ 4.6 ]    where @xmath129 is the einstein - frame hubble parameter .
it is obvious from ( [ acc ] ) that a vanishing potential @xmath130 implies @xmath131 , so that the universe is decelerating in the einstein frame .
however , because of the relation @xmath132 , see eq .
( [ avsa * ] ) , the _ observed _ ( jordan - frame ) expansion rate @xmath133 may be positive even in this case , and we will see concrete examples in sec .
vi.a below .
this is an important point to remember : although we are looking for cosmological frw backgrounds whose expansion is accelerating , the sign of @xmath134 is _ a priori _ not fixed .
the scalar - field equation of motion in the ef follows from eq .
( [ boxvarphi ] ) , and reads @xmath135 it is also similar to the usual klein - gordon equation , with the notable difference of a source term on the right - hand side , with the coupling strength @xmath136 defined in eq .
( [ alpha ] ) above .
it is tempting to tackle our problem in the ef as the equations are simpler and we can rely on experience gained in general relativity .
however , a crucial difficulty that we encounter is that all physical quantities which appear in the ef background equations are not those that come from observations .
moreover , the behavior of matter in the ef is complicated by the relations ( [ rho*vsrho ] ) : instead of the simple power law @xmath137 for dustlike matter in the jf , one gets @xmath138 in the ef , where @xmath139 can have _ a priori _ any shape . to avoid these problems
, we will thus work in the jf , and show that the `` reconstruction '' program can equally well be implemented , like in general relativity , although it is mathematically very different .
we will nevertheless check at the end the consistency of the solutions obtained by translating them in terms of ef quantities .
we now consider the perturbations in the longitudinal gauge .
for this problem , we will restrict our discussion to the case of a spatially flat frw universe ( @xmath140 ) , and write the jf and ef metrics as    @xmath141    [ 4.8 ]    in the ef , the perturbation equations deriving from eq .
( [ einstein * ] ) are strictly the same as in general relativity plus a minimally coupled scalar field .
one thus finds notably @xmath142 . on the other hand ,
the equations for scalar - field and matter perturbations are modified by the matter - scalar coupling , proportional to @xmath136 in eqs .
( [ boxvarphi ] ) and ( [ matter * ] ) .    for our purpose
, it will be more useful to write the perturbation equations in the ( physical ) jf .
let us define the gauge invariant quantity introduced in @xcite : @xmath143 .
] @xmath144 where @xmath145 is the matter peculiar velocity potential ( such that @xmath146 is the perturbation of the four - dimensional unit velocity @xmath147 ) .
we now work in fourier space , and assume a spatial dependence @xmath148 , with @xmath149 .
the conservation equations of matter ( [ matter ] ) give    @xmath150    [ 4.10 ]    on the other hand , the einstein equations ( [ einstein ] ) give    @xmath151    [ 4.11 ]    note that @xmath152 in the jf , in contrast to the corresponding problem in general relativity or in the ef .
equation ( [ psi ] ) is actually an obvious consequence of the relation between @xmath50 and @xmath30 , eq .  (
[ g * ] ) , and of the fact that @xmath153 .
finally , eq .  ( [ boxphi ] ) yields the equation for the dilaton fluctuations @xmath154 : @xmath155\delta\phi = \nonumber\\ & & = \left[{k^2 \over a^2}(\phi-2\psi ) - 3(\ddot\psi+4h\dot\psi+h\dot\phi)\right]{1\over z}\,{df\over d\phi } + ( 3\dot\psi+\dot\phi)\dot\phi - 2\,{\phi\over z}\,{du\over d\phi}\ . \label{fluctphi}\end{aligned}\ ] ] in the particular representation @xmath103 used in ref .
@xcite , this equation reduces to the simpler form @xmath156\delta\phi = \nonumber\\ & = & \left[{k^2 \over a^2}(\phi-2\psi ) - 3(\ddot\psi+4h\dot\psi+h\dot\phi)\right]{df\over d\phi } + ( 3\dot\psi+\dot\phi)\dot\phi - 2\phi\,{du\over d\phi}\ . \label{fluctphiz1}\end{aligned}\ ] ]
the reconstruction of the potential @xmath10 was shown in @xcite to be possible in the framework of general relativity plus a minimally coupled scalar field , the @xmath4-field or quintessence , provided the hubble diagram ( and thus also @xmath20 ) can be extracted from the observations .
an essential difference arises when one deals with scalar - tensor theories : we have to reconstruct two unknown functions instead of one , hence we need to extract two quantities ( as functions of the redshift @xmath157 ) from the observations .
actually , in the minimally coupled case , the knowledge of the luminosity distance @xmath158 and of the clustering of matter @xmath159 , both in function of @xmath19 , provides two independent ways to reconstruct the scalar field potential @xcite . without any ambiguity in the minimally coupled case
, one needs to know both @xmath11 and the present energy density of dustlike matter @xmath160 , or both @xmath24 and the present value of the hubble constant @xmath5 . in our general scalar - tensor case ,
we need to know the two functions @xmath11 and @xmath24 , but no independent measurement of @xmath160 or @xmath5 is necessary . ] in our case , both quantities are necessary and the reconstruction itself is significantly more complicated .
the present section generalizes our previous results of ref .
@xcite not only by considering the most general parametrization ( [ s_jf ] ) of scalar - tensor theories and by taking into account the possible spatial curvature of the universe , but also by discussing particular cases that were excluded in this reference . from now on
, we will restrict our discussion to the case of a pressureless perfect fluid ( @xmath161 ) , because all matter in the universe will be assumed to be simply dustlike , of course besides that part needed to account for the present accelerated expansion ( _ i.e. _ , the scalar field in the present framework ) .
the first step of the reconstruction program is the same as in general relativity , since it is purely kinematical and does not depend on the field content of the theory : if the luminosity distance @xmath158 is experimentally determined as a function of the redshift @xmath19 , one can deduce the quantity @xmath20 from the relation @xmath162^{-1/2 } , \label{hvsdl}\ ] ] where the prime denotes the derivative with respect to @xmath19 .
the large square brackets contain a corrective factor involving the present energy contribution @xmath163 of the spatial curvature of the universe .
it was not written explicitly in refs .
@xcite , which focused their discussions on the flat - space case ( @xmath164 ) , but it is a straightforward consequence of eqs .
( 23)(25 ) of ref .
since present experimental data suggest that @xmath165 is small , the flat - space expression for @xmath166'$ ] is _ a priori _ a good approximation anyway .
note that even if one uses the exact equation ( [ hvsdl ] ) , it reduces to the flat - space expression for @xmath167 ( because @xmath168 ) , and therefore @xmath5 is always known without any ambiguity . to determine @xmath20 precisely at higher @xmath19 ,
one then needs to know both @xmath11 and @xmath169 .    by eliminating @xmath170 from the background equations ( [ h2 ] ) and ( [ doth ] ) , we then obtain the equation @xmath171 which , when rewritten in terms of the redshift @xmath19 , gives the fundamental equation @xmath172~f ' + \left[\frac{6}{(1+z)^2 } - \frac{2}{1+z}(\ln h)'-4\left({h_0\over h}\right)^2\omega_{\kappa,0}\right]~f = & & \nonumber\\ = \frac{2 u}{(1+z)^2 h^2 } + 3 ( 1+z ) \left({h_0\over h}\right)^2 f_0 \omega_{m,0}\ . & & \label{f}\end{aligned}\ ] ] as before , an index @xmath76 means the present value of the corresponding quantity , and we use again the notation @xmath173 . in this equation
, @xmath174 stands for the present energy density of dustlike matter relative to the critical density @xmath175 . to simplify
, this critical density is defined in terms of the present value of newton s gravitational constant ( [ g_n ] ) , @xmath176 , instead of the effective gravitational constant ( [ geff ] ) actually measured in cavendish - type experiments . indeed ,
solar - system experiments tell us that their present values differ by less than @xmath92 , as discussed in sec .
[ note in passing that by changing the value of @xmath28 , one can always set @xmath177 without loss of generality . ]    in conclusion , we are left with a non - homogeneous second order differential equation for the function @xmath178 , a situation very different from that prevailing in general relativity .
however , the right - hand side also depends on the unknown potential @xmath179 , so that this equation does not suffice to fully reconstruct the microscopic lagrangian of the theory . as we will show in sec .
vi below , it can nevertheless be used for a systematic study of several scalar - tensor models , provided one of the two unknown functions is given ( or a functional dependence between them is assumed ) .
this can be useful as we do not expect a simultaneous release of data yielding @xmath20 and @xmath24 .
we will see that such a study already yields powerful constraints on the family of theories which are viable .    on the other hand , if @xmath24 is also experimentally determined , and if we assume a spatially flat frw universe ( @xmath180 ) , we will see in the next subsection ( v.b ) that the value of @xmath160 as well as the function @xmath178 can be obtained independently of @xmath179 .
equation ( [ f ] ) then gives @xmath179 in an _ algebraic _ way from our knowledge of @xmath20 , @xmath178 and @xmath160 .
let us now assume that both @xmath178 and @xmath179 are known , either because one of them was given from theoretical naturalness assumptions , or because @xmath24 has been experimentally determined with sufficient accuracy .
we will also assume that both @xmath160 and @xmath169 are known .
it is then straightforward to reconstruct the various functions of @xmath9 entering the microscopic lagrangian ( [ s_jf ] ) . in the brans - dicke representation
, one has @xmath45 , therefore the knowledge of @xmath178 and @xmath179 suffices to reconstruct the potential @xmath10 in a parametric way .
however , to fully determine the theory , one also needs to know @xmath181 , or equivalently an equation giving the @xmath19-dependence of @xmath124 . on the other hand , in the simpler representation @xmath182 and
@xmath33 unknown , we need an equation giving the @xmath19-dependence of @xmath9 to reconstruct @xmath33 and @xmath10 parametrically .
these two cases , as well as any other possible parametrization of the theory , are solved thanks to eq .
( [ doth ] ) above , which reads in function of the redshift @xmath183~f ' + 2\left[\frac{(\ln h)'}{1+z } - \left({h_0\over h}\right)^2\omega_{\kappa,0}\right]~f \nonumber\\ & & - 3(1+z)\left({h_0\over h}\right)^2 f_0 \omega_{m,0}\ , \label{phi}\end{aligned}\ ] ] or equivalently @xmath184 in the @xmath182 representation , @xmath185 is thus obtained by a simple integration . in the brans - dicke representation , on the other hand
, @xmath186 is given by an _ algebraic _ equation in terms of @xmath20 , @xmath187 , and their derivatives .
it is rather obvious but anyway important to note that if the microscopic lagrangian ( [ s_jf ] ) can be reconstructed in the jf , it can also be directly obtained in the ef , eq .
( [ s_ef ] ) .
this allows us to check the mathematical consistency of the theory , and notably if the helicity-0 degree of freedom @xmath51 always carries positive energy .
one can also verify that the function @xmath188 defining the coupling of matter to the scalar field is well defined , and notably single valued .
finally , the second derivative of the potential @xmath130 also gives us the sign of the square of the scalar mass , and negative values would strongly indicate an instability of the model .
these important features can not easily be checked in the jf , because the sign of @xmath36 in eq .
( [ s_jf ] ) is not directly related to the positivity of the scalar - field energy ( see below ) , and also because the second derivative of @xmath10 does not give the precise value of its squared mass .
[ as shown by eq .
( [ v ] ) , the helicity-0 degree of freedom @xmath51 may have a mass , @xmath189 , even if @xmath10 is strictly constant , provided @xmath33 varies . ]
let us thus assume that @xmath20 , @xmath160 and @xmath169 are known , and that @xmath178 and @xmath179 were reconstructed as above .
equation ( [ a ] ) then gives @xmath190 , _ i.e. _ , the einstein - frame coupling factor @xmath95 as a function of the _ jordan - frame _ redshift @xmath19 ( which is the redshift we observe ) .
combining now eq .
( [ varphi ] ) with ( [ phi ] ) , we get    @xmath191{f'\over f } + { ( \ln h ) ' \over 1+z}\nonumber\\ & & -\left({h_0\over h}\right)^2\omega_{\kappa,0 } -{3\over 2}(1+z)\left({h_0\over h}\right)^2 { f_0\over f}\ , \omega_{m,0}\ , \label{varphiprime}\end{aligned}\ ] ]    [ 5.5 ]    or also , not eliminating the potential @xmath192 @xmath193 the ef scalar @xmath51 is thus also known as a function of the _ jordan - frame _ redshift @xmath19 ( up to an additive constant @xmath194 which can be chosen to vanish without loss of generality ) , and one can reconstruct @xmath188 in a parametric way .
similarly , the ef potential @xmath130 , eq .
( [ v ] ) , can be reconstructed from our knowledge of @xmath178 , @xmath179 and @xmath195 .
since @xmath51 describes the actual helicity-0 degree of freedom of the theory , this field must carry only positive energy excitations , and @xmath196 must be positive .
on the other hand , the tensor and scalar degrees of freedom are mixed in the jf , and the positivity of energy does not imply that @xmath197 should always be positive .
actually , eq .  ( [ positiveenergy ] ) shows that it can become negative when @xmath198 happens to be larger than @xmath199 , which can occur in perfectly regular situations .
[ we will see an explicit example in sec .
vi.a below .
] this underlines that the parametrization @xmath182 can sometimes be singular : the derivatives of @xmath9 may become purely imaginary although the scalar degree of freedom @xmath51 is well defined . on the other hand ,
the brans - dicke representation is well behaved ( @xmath200 remains always positive ) , and the positivity of energy simply implies the well - known inequality @xmath201 . actually , the particular value @xmath202 is also singular , as it corresponds to an infinite coupling strength @xmath203 between matter and the helicity-0 degree of freedom @xmath51 .
the domain for which the @xmath103 parametrization is pathological although the theory remains consistent simply corresponds to @xmath204 , or @xmath205 .      although the perturbations will not be used in sec .
vi below , we emphasize that the phenomenological reconstruction of the full microscopic lagrangian can be implemented without any ambiguity if fluctuations are taken into account . for completeness
, we review now this part of our program .
we assume that both @xmath20 and the matter density perturbation @xmath24 are experimentally determined with enough accuracy , and as in sec .
iv.b above , we focus our discussion on the case of a spatially flat frw universe ( @xmath180 ) .
we also assume that matter is dustlike ( @xmath206 ) , and the perturbation equations of sec .
iv.b are thus simplified .
in particular , eq .  ( [ conserv_z ] ) reduces to the mere identity @xmath207 .
we consider comoving wavelengths @xmath208 much shorter ( for recent times ) than the hubble radius @xmath209 , and also shorter than the inverse mass of the scalar field : @xmath210 two different reasonings can now be used to reach the same conclusions .
the first one , explained in ref .
@xcite , consists in taking the formal limit @xmath211 in the various perturbation equations .
then , the leading terms are either those containing @xmath159 or those multiplied by the large factor @xmath212 .
one also needs to consider only the growing adiabatic mode of eq .
( [ fluctphi ] ) , for which @xmath213 .
the other reasoning needs a simpler ( but _ a priori _ stronger ) hypothesis .
one assumes that the logarithmic time derivative of any quantity , say @xmath214 , is at most of order @xmath49  : @xmath215 suffices to derive straightforwardly all the following approximations .
note that both reasonings correspond in fact to the same physical situation of a weakly - coupled light scalar field . in the case of a strongly - coupled but very massive scalar ( see the second paragraph of sec .
iii ) , the equations can not be approximated as shown below , and the time evolution of density fluctuations does not follow the same law .
for instance , in the particular model considered in ref .
@xcite , one always finds a strong clustering of the scalar field at small scales .
indeed , this model corresponds to the choice @xmath45 and @xmath127 in action ( [ s_jf ] ) , and eq .  ( [ fluctphi ] ) can then be rewritten as @xmath216 .
therefore , even if the scalar field is very massive ( @xmath70 large ) , one finds that it is anyway strongly clustered for comoving wavelengths @xmath217 shorter than the inverse mass , _
i.e. _ , in the formal limit @xmath211 .
although this is _ a priori _ not forbidden by observations of gravitational clustering , since the inverse mass must be much smaller than the astronomical unit in this model , this is anyway an indication of its probable instability .
we will not consider such heavy scalar fields any longer in this paper , and we now come back to the class of weakly - coupled light - scalar models , which are the most natural alternatives to general relativity .    setting @xmath218 and making use of ( [ conserv_z ] )
, one can write ( [ conserv_0 ] ) as @xmath219 where the right - hand side is negligible with respect to each separate term of the left - hand side because of the above hypotheses .
note that ( [ ddotdeltam ] ) just reproduces the standard evolution equation for matter perturbations . using ( [ fluctphiz1 ] ) , we also arrive at @xmath220 where the second equality is a consequence of eq .
( [ psi ] ) . in the case of gr plus a minimally coupled scalar field ,
one finds that @xmath221 in the limit @xmath211 , so that the scalar field is not gravitationally clustered at small scales @xcite .
this is in agreement with the observational fact that the dark matter described by the @xmath4-term should remain unclustered up to comoving scales @xmath222 mpc ( where we recall that @xmath223 ) . on the other hand , in our scalar - tensor framework ,
( [ dphi ] ) shows that the scalar field _ is _ clustered at arbitrarily small scales , but only weakly because the derivative @xmath224 is experimentally known to be small [ see the solar - system constraint ( [ df0 ] ) , and the limit @xmath225 justified in @xcite for redshifts @xmath226 .
the class of models we are considering , involving a light scalar field weakly coupled to matter , is thus also in agreement with observations of gravitational clustering .    finally , still under the above hypotheses , eq .
( [ fluctmetr ] ) implies @xmath227 remembering the definition ( [ geff ] ) for @xmath91 , and using ( [ dphi ] ) above , eq .
( [ poi ] ) can be recast in a form which exhibits its physical content : @xmath228 poisson s equation is thus simply modified by the substitution of newton s constant @xmath229 by @xmath91 , the effective gravitational constant between two close test masses !
this conclusion was also reached in @xcite , but only for brans - dicke theory with a constant parameter @xmath230 , while we have derived it for an arbitrary ( light ) scalar - tensor theory . as discussed in sec .
iii above , expression ( [ geff ] ) is valid only if the distance between the test masses is negligible with respect to the inverse scalar mass .
the physical reason why this expression appears in poisson s equation ( [ poi2 ] ) is just that we are working in the short wavelength limit ( [ largek ] ) : the frequency of the waves we are considering is so large that the scalar field behaves as if it were massless .    combining ( [ ddotdeltam ] ) with ( [ poi2 ] )
, we now arrive at our final evolution equation for @xmath159  : @xmath231 in terms of the redshift @xmath19 , this reads @xmath232 provided we can extract from observation both physical quantities @xmath20 and @xmath24 with sufficient accuracy , the explicit reconstruction of the microscopic lagrangian is obtained in the following way . starting from ( [ delz ] ) and using the fact that _ today _ @xmath233 and @xmath234 differ by less than @xmath92 , eq .  ( [ delz ] ) evaluated at present gives us the cosmological parameter @xmath160 with the same accuracy . then , returning to eq .
( [ delz ] ) for arbitrary @xmath19 , we get @xmath235 , where @xmath236 is a known function of the observables @xmath20 , @xmath24 , and their derivatives . using now eq .  ( [ phi ] ) and expression ( [ geff ] ) for @xmath237 , we get a nonlinear second order differential equation for @xmath178 , which can be solved for given @xmath238 and @xmath239 [ one can always set @xmath177 without loss of generality , while @xmath239 is constrained by eq .
( [ df0 ] ) ] .
after we have found @xmath178 , we can plug it into ( [ f ] ) to determine @xmath179 in an algebraic way .
the final step is explained in the previous subsection , above eq .
( [ phi ] ) , for the various possible parametrizations of action ( [ s_jf ] ) : in the @xmath103 parametrization , @xmath185 is obtained by a simple integration of eq .
( [ phi ] ) , while in the brans - dicke parametrization ( @xmath37 ) , @xmath186 is given algebraically by the same eq .
( [ phi ] ) .
this enables us to reconstruct @xmath33 ( or @xmath240 ) and @xmath10 as functions of @xmath241 for that range corresponding to the data .    actually , for sufficiently low redshifts @xmath242 , eq .
( [ delz ] ) can be simplified without losing too much accuracy . indeed , as shown in ref .
@xcite , the square of the matter - scalar coupling strength @xmath243 , eq .
( [ alpha ] ) , is at most of order @xmath244 for such redshifts . moreover , under natural assumptions , much smaller values of @xmath245 are generically predicted in scalar - tensor theories @xcite . therefore , @xmath91 and @xmath82 differ by less than @xmath101 for redshifts
@xmath242 , and eq .  ( [ delz ] ) can be used to obtain @xmath246 with the same accuracy .
the interest of this simplification is that @xmath178 is now given by an _ algebraic _ equation . in the brans - dicke representation ,
_ all _ the steps of the reconstruction program are thus algebraic , eq .  ( [ f ] ) giving @xmath179 , and eq .
( [ phi ] ) giving @xmath186 .
the only non - algebraic step is the final parametric reconstruction of @xmath10 and @xmath240 .
let us end this section by a few comments on the observational accuracy which will be needed for this reconstruction program to be implemented .
first , eq .
( [ delz ] ) allows to reconstruct @xmath178 only if @xmath247 and @xmath248 are both determined with enough accuracy .
moreover , the second derivative of this reconstructed @xmath178 is needed in eq .
( [ f ] ) to obtain @xmath179 .
therefore , the actual reconstruction of the potential depends _ a priori _ on the fourth derivative of @xmath24 , so that extremely clean data seem to be necessary
. however , the situation is better than this naive derivative counting suggests .
indeed , the above estimates for @xmath245 show that @xmath178 does not vary much on the redshift interval @xmath249 .
therefore , the first two terms of eq .
( [ f ] ) , involving @xmath250 and @xmath251 , are expected to be negligible with respect to the third one involving @xmath96 .
a noisy experimental determination of @xmath252 and @xmath253 is thus not a serious difficulty for our reconstruction program . on the other hand ,
clean enough data are still needed to determine @xmath178 from eq .
( [ delz ] ) , using @xmath24 and its first two derivatives .
before such clean data are available , it will be sufficient to verify that eq .
( [ delz ] ) is consistent with a slowly varying @xmath178 . in the next section
, we will show that interesting theoretical constraints can anyway be obtained without knowing at all the density fluctuation @xmath24 , but using only the luminosity distance @xmath11 and consistency arguments within particular subclasses of scalar - tensor models .
@xcite , a fit of presently known supernovae events has been performed to obtain the luminosity distance @xmath11 up to redshifts @xmath1 , of course still with large uncertainties .
although this is not yet sufficient to constrain seriously scalar - tensor models , we can expect clean data on @xmath11 in the near future from additional supernovae events , and anyway earlier than for the density perturbations @xmath24 .
the snap satellite will in particular observe thousands supernovae events up to @xmath12 . in this section
, we will concentrate on the theoretical constraints that can be extracted from the knowledge of @xmath11 alone , and therefore of @xmath20 using eq .
( [ hvsdl ] ) .
we will thus only use the results of subsection v.a above . since the knowledge of this function does not suffice to fully reconstruct the microscopic lagrangian ( [ s_jf ] )
, we will need additional assumptions on one of the functions it involves , either @xmath96 ( or @xmath124 , depending on the parametrization ) or the potential @xmath192 .
one may also assume a functional relation between @xmath96 and @xmath192 ( for instance @xmath254 as in ref .
@xcite ) .    to emphasize as clearly as possible
what kind of constraints can be imposed on scalar - tensor theories , we shall consider the worst situation for them .
let us assume that the observed function @xmath20 will be exactly given by eq .
( [ h2 ] ) for @xmath140 , @xmath255 , and @xmath256  : @xmath257 of course , such an observation would _ a priori _ call for the following standard interpretation : gravity is correctly described by general relativity , and we live in a flat universe filled with dustlike matter and a cosmological constant , with corresponding present energy densities ( relative to the critical density ) @xmath160 and @xmath258 .
however , for our purpose , eq .
( [ h2gr ] ) should just be considered as kinematical .
it tells us how the universe expands with redshift @xmath19 , but we are free to assume that the dynamics of the expansion is governed by a scalar - tensor theory .
therefore , @xmath160 and @xmath258 are here mere parameters , whose names refer to their physical significance in the framework of gr . of course
, one should not forget that they do not have the same interpretation within scalar - tensor theories .    for our numerical applications
, we will further take the present estimates based on combined cmb fluctuations and supernovae observations ( they will be determined more accurately by future experiments ) : @xmath259 for these numerical values , ( [ h2gr ] ) is consistent with the presently available luminosity distance @xmath11 up to @xmath1 .
actually the best - fit universe , if we assume flatness , gives @xmath260 and @xmath261 .
we have chosen to work directly with the exact form ( [ h2gr ] ) , instead of the @xmath11 extracted from observation , in order to clarify the physical content of our results . indeed , the present observational estimates for @xmath11 are still too imprecise to constrain strongly the class of scalar - tensor theories we are considering . moreover , some of our results below depend crucially on the fact that @xmath20 keeps the form ( [ h2gr ] ) up to redshifts @xmath262 , which have not yet been reached experimentally . to relate our results to those obtained in @xcite using fitting functions or an expansion in powers of @xmath19 ,
one just needs to use eq .
( [ hvsdl ] ) : our exact expression ( [ h2gr ] ) for @xmath20 corresponds to some exact expression for @xmath11 .    to summarize , we are assuming in this section that future observations of the luminosity distance @xmath11 will provide a @xmath20 of the form ( [ h2gr ] ) with the numerical values ( [ omegavalues ] ) .
this implies notably that our universe is presently accelerating . on the other hand , we are _ not _ assuming that the correct theory of gravity is necessarily gr plus a cosmological constant .
the main question that we will address is therefore the following : would such an `` observed '' @xmath20 necessarily rule out the existence of a scalar partner to the graviton ? if not , would it be possible to reproduce ( [ h2gr ] ) within a more natural scalar - tensor theory , in which @xmath258 could be explained by a `` generalized quintessence '' mechanism ?
we will first analyze in subsection a the simplest subclass of scalar - tensor theories that we can consider , namely when @xmath263 in action ( [ s_jf ] ) .
since this is _ a priori _ the subclass which differs the most from gr plus a cosmological constant , this study will be rather detailed , and it will allow us to underline the mathematical and physical meaning of the constraints that are obtained .
subsection b will be again devoted to the case of a massless scalar field , but combined with a cosmological constant .
as its conclusions basically confirm those of subsection a , we will present them more concisely .
finally , we will briefly discuss in subsection c the cases where one imposes particular forms for the coupling function @xmath96 in action ( [ s_jf ] ) , and one reconstructs the potential @xmath192 from the background equations ( [ f])([phi ] ) .
the case of a given functional dependence between @xmath96 and @xmath192 will also be addressed .      since a cosmological constant
can be interpreted as a particular case of scalar - field potential , it is instructive to analyze whether an observed expansion like ( [ h2gr ] ) could be reproduced in a theory _ without _ any potential , and we now study eqs .
( [ f])([5.5 ] ) for @xmath264 .
this case can be analyzed using the second order differential equation ( [ f ] ) for @xmath96 , which simplifies significantly if one introduces a function @xmath214 such that @xmath265 [ as mentioned in sec .
v.a above , one can also set @xmath177 without loss of generality .
] then , using the assumed `` experimental '' expression ( [ h2gr ] ) for @xmath20 , and writing ( [ f ] ) in terms of @xmath266 , we get @xmath267 to avoid any confusion , let us recall that @xmath258 ( and the two occurrences of @xmath160 in the left - hand side ) comes from the `` observed '' cosmological function ( [ h2gr ] ) , notwithstanding the fact that there is _ no _ cosmological constant in the model we are considering .
the value @xmath160 appearing in the right - hand side stands for the present relative energy density of dustlike matter .
we assume that it takes the same numerical value ( [ omegavalues ] ) as in the `` observed '' @xmath20 ( [ h2gr ] ) .
equation ( [ eqf ] ) tells us how we should choose @xmath268 to mimic exactly this @xmath20 in the present potential - free theory . in other words ,
@xmath258 and @xmath160 are two numbers assumed to be given by experiment , and we wish to fit @xmath268 and @xmath169 to satisfy eq .
( [ eqf ] ) .    to integrate this second - order differential equation , we need two initial conditions for @xmath214 and its derivative .
the first one is an obvious consequence of eq .
( [ smallf ] ) taken at @xmath269 , and we simply get @xmath270
. the second one should be such that the solar - system bound ( [ df0 ] ) is satisfied .
for instance , if @xmath271 does not vanish , it is _ sufficient _ to impose @xmath272 , _ i.e. _ , @xmath273 using eq .
( [ smallf ] ) .
this corresponds to a scalar - tensor theory which has been attracted towards an extremum of @xmath96 during the cosmological expansion of the universe ( _ cf .
_ @xcite ) , so that it is presently strictly indistinguishable from general relativity in solar - system experiments .
[ the full allowed domain for @xmath274 will be explored below in a numerical way . ]
we consider first our potential - free model in a spatially flat frw universe ( @xmath164 ) .
then eq .
( [ eqf ] ) becomes a first - order differential equation for @xmath275 , and its integration yields @xmath276\ , , \label{fprime}\ ] ] where we have set @xmath277 , and where the final constant inside the square brackets has been chosen to impose @xmath273 ( _ i.e. _ , @xmath278 ) .
the function @xmath279 can be explicitly written in terms of generalized hypergeometric functions , but its complicated expression will not be useful for our purpose .
let us just quote the first order of its expansion in powers of @xmath280  : @xmath281 in conclusion , eq .
( [ eqf ] ) could be integrated analytically , in the particular case of a spatially flat universe .
this means that at least in the vicinity of @xmath269 , there _ a priori _ exists a potential - free scalar - tensor theory which exactly mimics general relativity plus a cosmological constant .
however , the theory is mathematically consistent only if @xmath178 remains strictly positive .
[ if @xmath96 vanishes , then the coupling function @xmath188 , eq .  ( [ a ] ) , between matter and the helicity-0 degree of freedom @xmath51 diverges , and if @xmath96 becomes negative , the graviton carries negative energy .
] let us thus compute the value @xmath282 for which @xmath283 , or @xmath284 , vanishes for the first time .
because of the complexity of the solution @xmath214 , we did not find a close analytical expression for @xmath282 , but its expansion in powers of @xmath280 can be obtained straightforwardly : @xmath285 { \omega_{m,0}^2\over\omega_{\lambda,0}^2 } + o\left({\omega_{m,0}^3\over\omega_{\lambda,0}^3}\right ) .
\label{zmax}\ ] ] numerically , for the values ( [ omegavalues ] ) of @xmath258 and @xmath160 , we find @xmath286 . in conclusion , this scalar - tensor model is able to mimic general relativity plus a cosmological constant , but only on the small interval @xmath287 . if future observations of type ia supernovae give a behavior of @xmath20 of the form ( [ h2gr ] ) on a larger interval , say up to @xmath1 , then the present scalar - tensor theory will be ruled out .
this example of a vanishing potential illustrates a conclusion that we will reobtain below for more general theories : the determination of the form of @xmath20 over some ( even rather small ) redshift interval is in fact more constraining than the precise value of the parameters @xmath288 themselves .
indeed , eq .
( [ zmax ] ) clearly shows that @xmath282 can not exceed 1 even in the presumably unrealistic case of @xmath289 .
[ a calculation using the exact expression for @xmath268 shows that @xmath290 would exceed 1 only for @xmath291 . ] note that all the results obtained are independent of the parameter @xmath5 .
one could try to increase @xmath282 by considering a spatially curved frw universe .
we did not solve eq .
( [ eqf ] ) in the most general case , but since we wish to compute the corrections to eq .
( [ zmax ] ) due to a small value of @xmath292 , it is sufficient to work at zeroth order in @xmath280 .
let us thus set @xmath293 in eq .
( [ eqf ] ) , which reduces to @xmath294 its solution is obviously a sine if @xmath295 ( _ i.e. _ , @xmath296 , closed universe ) , or a hyperbolic sine for @xmath297 ( _ i.e. _ , @xmath106 , open universe ) .
taking into account the initial conditions @xmath270 and @xmath298 , we thus get    @xmath299    [ fok0 ]    the first zero of @xmath300 is then reached either at @xmath301 or at @xmath302 . in both cases ,
the expansion in powers of @xmath303 gives @xmath304 .
working perturbatively , one can also compute the correction to this expression due to the nonzero value of @xmath160 , and one finds that @xmath282 is given by eq .
( [ zmax ] ) above plus the following correction : @xmath305 { \omega_{\kappa,0}\over\omega_{\lambda,0 } } + o\left({\omega_{\kappa,0}^2\over\omega_{\lambda,0}^2}\right)\ , .
\label{deltazmax}\ ] ] in conclusion , @xmath282 can be slightly enlarged if we consider our potential - free scalar - tensor theory in an open frw universe ( @xmath297 ) .
numerically , for the values ( [ omegavalues ] ) of @xmath258 and @xmath160 , we find @xmath306 .
since the latest experimental data on cmb temperature fluctuations already constrain @xmath165 to be small ( see the latest boomerang and maxima data ) , and actually an open universe is unlikely while a marginally closed universe is still acceptable , we thus recover the same qualitative conclusion as in the spatially flat case : it is possible to mimic general relativity plus a cosmological constant within a potential - free scalar - tensor theory only on a small redshift interval @xmath307 .
the above conclusions have been confirmed by numerical integrations of eqs .
( [ f])([5.5 ] ) , still assuming a hubble diagram consistent with ( [ h2gr ] ) . instead of considering only theories which are presently indistinguishable from general relativity ( @xmath278 )
, we imposed arbitrary initial conditions for @xmath250 , and computed the corresponding value of the present scalar - matter coupling strength @xmath308 , eq .
( [ alpha ] ) . in the case of a spatially flat frw universe
, we recovered that the solar - system bound ( [ df0 ] ) imposes the limit @xmath309 , consistently with the above analytical estimate ( [ zmax ] ) .
in other words , the constraint ( [ df0 ] ) is so tight that even taking the largest allowed value for @xmath310 does not change significantly @xmath282 .
figure 1 displays the reconstructed @xmath178 for this maximal @xmath310 , and one can note that its slope at @xmath167 is visually indistinguishable from the horizontal .
this figure also plots the einstein - frame scalar @xmath51 , eq .
( [ varphi ] ) , which is the actual helicity-0 degree of freedom of the theory .
notice that it diverges at @xmath282 , so that the theory loses its consistency beyond this value of the redshift .
curiously , we found that even if no experimental constraint like ( [ df0 ] ) is imposed on @xmath310 ( _ i.e. _ , even if we forget that solar - system experiments confirm very well general relativity ) , then the mathematical consistency of the theory anyway imposes @xmath311 .
in fact , eq .  ( [ f ] ) alone can be solved for arbitrary large values of @xmath19 , _
i.e. _ , there exist initial values of @xmath239 such that @xmath178 remains positive for any @xmath19 .
however , the values of @xmath239 needed to integrate eq .
( [ f ] ) beyond @xmath312 correspond to negative values of @xmath313 ( where @xmath314 denotes the present value of the brans - dicke parameter ) . in other words ,
the expression of @xmath196 given by eq .
( [ 5.5 ] ) would become negative around @xmath167 , and the helicity-0 degree of freedom would thus need to carry negative energy at least on a finite interval of @xmath19 , if one wished to integrate eqs .
( [ f])([5.5 ] ) beyond @xmath315 .
figure 2 displays the maximum redshift @xmath282 consistent with the positivity of energy of both the graviton and the scalar field , but for any value of the present matter - scalar coupling strength @xmath310 . as underlined above
, one finds that @xmath290 can never be larger than @xmath316 .
this figure also indicates the present solar system bound on @xmath310 , corresponding to @xmath309 as in fig .  1
. the limiting case of a vanishing @xmath310 , _
i.e. _ , of a scalar - tensor theory which is presently _ strictly _ indistinguishable from gr in the solar system , corresponds to @xmath286 , as was derived analytically in eq .
( [ zmax ] ) .
figure 2 also indicates the range of values for @xmath310 that are generically obtained in refs .
@xcite while studying the cosmological evolution of scalar - tensor theories at earlier epochs in the matter - dominated era : the theory is attracted towards a maximum of @xmath96 ( _ i.e. _ , a minimum of @xmath317 ) so that the present value of @xmath310 is expected to be extremely small .
finally , this figure also displays the maximum value of @xmath310 for which the parametrization @xmath318 of action ( [ s_jf ] ) has a meaning . beyond @xmath319 ( _
i.e. _ , for a brans - dicke parameter @xmath320 ) , one would get @xmath321 in this parametrization .
in other words , eqs .
( [ f])([phi ] ) can not be integrated consistently beyond @xmath322 if one sets @xmath39 , whereas the brans - dicke or the einstein - frame representations show that the theory can be mathematically consistent up to @xmath323 ( @xmath199 remains positive ) .
this underlines that the @xmath103 parametrization may be sometimes pathological .    our numerical integration of eq .
( [ varphiprime ] ) not only allowed us to check the positivity of the scalar field energy , but also to reconstruct parametrically the matter - scalar coupling function @xmath188 . since @xmath324 , eq .
( [ a ] ) , we know that @xmath325 is finite and strictly positive over the interval @xmath326 , but we also checked that it is single valued over this interval .
this means that if @xmath195 can take several times the same value for different @xmath19 , they must correspond also to the same value of @xmath325 . actually , since eq .
( [ varphiprime ] ) does not fix the sign of @xmath327 , one should keep in mind that @xmath51 can oscillate around a constant value @xmath328 . if the numerical integration confuses the two points @xmath329 , but if @xmath188 happens not to be symmetrical around @xmath328 , it may look like a bi - valued function .
when such a situation occurred in our programs , we always verified that a single - valued @xmath188 could be defined consistently by unfolding it around the oscillation points of @xmath51 .
figure 3 illustrates such a situation , for an intentionally unrealistic value of @xmath310 in order to clarify the plots . [ the value @xmath330 is inconsistent with the solar - system bound ( [ df0 ] ) , but it corresponds anyway to a mathematically consistent theory , although the @xmath103 parametrization can not be used in this case . ]
all the functions @xmath317 that we reconstructed have similar convex parabolic shapes .
this is consistent with the results of refs .
@xcite , showing that the scalar field is generically attracted towards a minimum of @xmath317 during the expansion of the universe . if we had found models such that the present epoch ( @xmath167 ) is close to a _
maximum _ of @xmath331 , this would have meant that the theory is unstable , and that we have extremely fine tuned it to be consistent with solar - system constraints . on the contrary , the convex functions @xmath317 that we obtained show that these scalar - tensor models are cosmologically stable , _
i.e. _ , that the tight bounds ( [ df0 ] ) are in fact natural consequences of the attractor mechanism described in @xcite .
we have checked that reducing the parameter @xmath258 allows us to extend the integration region in the past , consistently with the above analytical results .
for instance , when we vary @xmath258 , still satisfying @xmath332 and setting @xmath164 , we find that @xmath333 is required in order to integrate the equations up to a redshift @xmath334
. this would correspond to @xmath335 , _
i.e. _ , 100 times larger than present estimates .
we also added random noise to our assumed @xmath20 , eq .
( [ h2gr ] ) , and verified that the conclusions are not changed qualitatively provided @xmath20 is known over a wide enough redshift interval .
this means that the experimental determination of the luminosity distance @xmath11 needs not be very precise to be quite constraining , provided redshifts of order @xmath262 are probed . as an illustration ,
let us take the exact expression @xmath20 of eqs .
( [ h2gr])([omegavalues ] ) for discrete values of the redshift , say @xmath336 , and let us add or subtract randomly between @xmath76 and @xmath337 to the corresponding @xmath20 .
then , we may fit a polynomial through these `` noisy '' values of @xmath20 , and use our numerical programs to integrate the background equations ( [ f])([5.5 ] ) and reconstruct @xmath96 .
we found that there always exists a maximum redshift beyond which @xmath96 is negative ( and the theory thus inconsistent ) .
figure 4 displays the two extreme values of @xmath282 that we obtained with hundreds of such `` deformed '' @xmath20 : it is sometimes even smaller than for the `` exact '' @xmath20 of eq .
( [ h2gr ] ) , and sometimes larger but never greater that @xmath338 .
it should be noted that for the @xmath337 noise we chose , the @xmath20 of pure gr with a vanishing cosmological constant could have been obtained . in that case , a potential - free scalar - tensor model with @xmath339 would of course have fitted perfectly this @xmath20 up to @xmath340 .
the reason why we never managed to go beyond @xmath341 is that we considered _ random noise _ , instead of such a precise _ bias _ of our assumed function @xmath20 , eq .
( [ h2gr ] ) .
we are aware that our deformed functions of fig .  4
do not reproduce a realistic experimental noise .
however , they illustrate in a well defined way that an inaccurate determination of @xmath20 over a wide redshift interval is actually more constraining than a precise measurement over a small redshift interval only .
the conclusion of the present subsection is therefore that a scalar - tensor theory without potential can accommodate a hubble diagram consistent with ( [ h2gr ] ) , but only on a small redshift interval if @xmath258 is significant .
the experimental determination of the luminosity distance @xmath11 , either accurately for @xmath242 or even with large ( tens of percents ) uncertainties up to redshifts @xmath262 , severely constrains this subclass of theories .
future observations should thus be able to distinguish them from general relativity , and to confirm or rule them out without any ambiguity .
it is worth noting that such future determinations of @xmath11 would _ a priori _ be much more constraining than solar - system experiments and binary pulsars tests . indeed , although the precision of the latter is quite impressive ( see e.g. @xcite ) , they anyway probe only the first two derivatives of @xmath317 , eqs .
( [ df0])-([beta0 ] ) , whereas cosmological observations should give access to the full shape of this function .
let us also recall that the constraints we found crucially depend on the fact that the theory should contain only positive - energy excitations to be consistent , and notably that the function @xmath96 should remain always strictly positive .
we did not use any other cosmological observation , but obviously , once the microscopic lagrangian of a scalar - tensor theory has been reconstructed using @xmath11 , all its other cosmological predictions should also be checked .
for instance , a bound @xmath342 is given in ref .
@xcite for the value of the function @xmath96 at nucleosynthesis time .
are shown to be consistent with the observed abundances of light elements , provided @xmath343 is large enough , where @xmath188 is the matter - scalar coupling function ( [ a ] ) . ] if one assumes that @xmath178 is monotonic , the reconstructed function of fig .
1 would not be consistent with this nucleosynthesis bound beyond @xmath344
. this would be even more constraining than the bound @xmath345 we obtained just from mathematical consistency requirements .
alternatively , a reconstructed function @xmath178 like the one of fig .  1 would be consistent with the above nucleosynthesis bound only if it were non - monotonic beyond @xmath346 . although this would not be forbidden from a purely phenomenological point of view , this would be anyway unnatural , and more difficult to justify theoretically .      to confirm the results of the previous subsection ,
let us now consider the case of a massless scalar field together with a cosmological constant whose value _ differs _ from the one entering our assumed @xmath20 , eq .
( [ h2gr])-([omegavalues ] ) .
the question that we wish to address is the following : can _ part of _ the observed @xmath258 be due to the presence of a massless scalar field ?    to impose a cosmological constant in a scalar - tensor theory , one would naively choose a constant potential @xmath10 in action ( [ s_jf ] ) .
however , as shown by eq .
( [ v ] ) , the corresponding potential @xmath130 of the helicity-0 degree of freedom @xmath51 would not be constant in this case ( because @xmath33 is _ a priori _ varying ) , and its second derivative would give generically a nonvanishing scalar mass . to avoid any scalar self - interaction , and in particular to set its mass to @xmath76 ,
one needs in fact to impose @xmath347 in the einstein - frame action ( [ s_ef ] ) .
this defines a consistent cosmological `` constant '' in a massless scalar - tensor theory .
note that the corresponding jordan - frame potential @xmath10 is then proportional to @xmath348 , and therefore that it does not correspond to the usual notion of cosmological _ constant _ in action ( [ s_jf ] ) .
since our assumed `` observed '' @xmath20 involves a parameter denoted @xmath258 , eqs .
( [ h2gr])-([omegavalues ] ) , let us introduce a different notation for the contribution due to the constant potential @xmath349  : @xmath350 it is easily checked that for @xmath351 , the solution @xmath352 ( or @xmath353 ) is recovered , _
i.e. _ , a scalar field minimally coupled to gravity with a constant potential acting like a cosmological constant . indeed , in terms of the function @xmath268 defined in ( [ smallf ] ) , eq .
( [ f ] ) reads @xmath354 note that this is now a _
non_-linear equation in @xmath96 , contrary to eq .
( [ eqf ] ) above for the case of a vanishing potential .
if @xmath351 , one finds that @xmath355 is an obvious solution , _
i.e. _ , @xmath356 a constant scalar field @xmath9 ( or @xmath51 ) then satisfies eqs .
( [ phi])([varphiprimebis ] ) .
if we now consider a scalar - tensor theory for which @xmath357 differs from the `` observed '' @xmath358 , we find that like in the previous subsection , there exists a maximum redshift @xmath282 beyond which @xmath178 becomes negative , and therefore beyond which the theory loses its mathematical consistency .
figure 5 displays this maximum redshift as a function of @xmath357 .
we plot this figure for the initial condition @xmath278 ( _ i.e. _ , for a theory which is presently indistinguishable from gr in the solar system ) , but as before , we verified that the curve is almost identical if one takes the maximum value of @xmath359 consistent with the solar - system bound ( [ df0 ] ) .
we also assume @xmath164 ( spatially flat universe ) for this figure , as we know from the previous discussion that a value even as large as @xmath360 does not change qualitatively the results .    for @xmath361
, we recover the result @xmath362 derived above for a vanishing potential . when @xmath363 , @xmath282 becomes even smaller . as expected , this is worse than in the potential - free case . on the contrary ,
when @xmath357 is positive ( _ i.e. _ , when it contributes positively to part of the `` observed '' @xmath258 ) , the maximum redshift @xmath282 increases .
this is just due to the fact that our massless scalar field needs to mimic a smaller fraction of the `` observed '' @xmath258 , so that the theory can remain consistent over a wider redshift interval .
however , we find that @xmath282 is still smaller than @xmath364 for @xmath365 , and a @xmath20 of the form ( [ h2gr])-([omegavalues ] ) observed up to @xmath262 would thus suffice to rule out the model .
if such a @xmath20 could be confirmed up to @xmath14 , one would need @xmath366 for our massless scalar - tensor theory to fit it !
even so , the theory would anyway become pathological at slightly higher redshifts . in conclusion ,
a massless scalar _ can not _ account for a significant part of the observed cosmological constant if @xmath20 is experimentally found to be of the form ( [ h2gr ] ) over a wide redshift interval .
let us note finally that for @xmath367 , eq .  ( [ eqfvconst ] ) does admit strictly positive solutions for @xmath214 ( or @xmath96 ) up to arbitrarily large redshifts .
this ensures that the graviton energy is always positive .
however , it is now the scalar field which needs to carry negative energy . indeed , eq .
( [ varphiprimebis ] ) gives a negative value for @xmath199 , basically because of the presence of the large negative number @xmath368 in this equation .
[ in the @xmath318 parametrization , @xmath200 is obviously also negative , because of eq .
( [ positiveenergy ] ) , or directly from eq .
( [ phibis ] ) which also involves a @xmath368 term . ]
therefore , there is only one possibility for a consistent massless scalar - tensor theory to reproduce ( [ h2gr ] ) over a wide redshift interval : it must involve a cosmological constant , whose contribution @xmath357 is equal to ( or very slightly smaller than ) the parameter @xmath258 entering ( [ h2gr ] ) .
in other words , the theory should be extremely close to gr plus a cosmological constant , and the massless scalar field must have a negligible contribution .
this illustrates again the main conclusion of our paper : the experimental determination of the luminosity distance @xmath11 over a wide redshift interval , up to @xmath262 , will suffice to rule out ( or confirm ) the existence of a massless scalar partner to the graviton .      in the previous two subsections , the matter - scalar coupling function @xmath33 ( or @xmath188 ) was reconstructed from the assumed knowledge of @xmath20 , for theories whose potential @xmath10 ( or @xmath130 ) had a given form .
we now consider the inverse problem .
we still assume that future observations will provide a hubble diagram consistent with ( [ h2gr])-([omegavalues ] ) , but we wish now to reconstruct the scalar - field potential @xmath192 for given forms of the coupling function @xmath96 .
we first consider a generic two - parameter family of scalar - tensor theories , which has already been studied in great detail for solar - system , binary - pulsar and gravity - wave experiments @xcite , as well as for cosmology starting with the matter - dominated era @xcite and even back to nucleosynthesis @xcite .
its definition is simplified if we work in the einstein frame ( [ 2.4])-([s_ef ] ) . the matter - scalar coupling function is simply given by @xmath369 in which the present value of the scalar field , @xmath194 , may be chosen to vanish without loss of generality .
any analytical function @xmath317 may be expanded in such a way , but we here assume that no higher power of @xmath51 appears , _
i.e. _ , that @xmath370 is strictly parabolic : it depends only on the two parameters @xmath308 and @xmath371 .
the latter is a simplified notation for @xmath77 , and should not be confused with the post - newtonian parameter @xmath372 defined in ( [ betappn ] ) .
[ actually , this equation shows that @xmath373 . ]
solar - system experiments impose @xmath374 , eq .
( [ df0 ] ) , while binary pulsars give @xmath375 , eq .
( [ beta0 ] ) , for this class of theories .
we first study these models for the case of a spatially flat universe ( @xmath164 ) .    as shown by eq .
( [ kg ] ) , a constant scalar field @xmath376 may be a solution if @xmath377 ( so that @xmath378 vanishes too ) and if the potential @xmath130 is also constant .
our assumed @xmath20 , eqs .
( [ h2gr])-([omegavalues ] ) , can thus always be reproduced if the parameter @xmath308 vanishes identically , and the reconstructed potential merely reduces to the constant @xmath379 .
this corresponds simply to gr plus a cosmological constant , and the massless scalar degree of freedom @xmath51 remains unexcited , frozen at an extremum of the parabola ( [ generica ] ) .
actually , eq .  ( [ kg ] ) shows that this extremum corresponds to a stable situation only if it is a minimum , _
i.e. _ , if @xmath380 in ( [ generica ] ) .
this is consistent with the results of refs .
@xcite : if the theory involves a cosmological constant whose value equals the `` observed '' one in eqs .  ( [ h2gr])-([omegavalues ] ) , a massless scalar field is cosmologically attracted towards a minimum of the coupling function @xmath317 , and the present value of its slope , @xmath308 , is expected to be generically very small .    on the other hand ,
if @xmath308 is _ not _ assumed to vanish , say if its value is comparable to the solar - system bound ( [ df0 ] ) , then our reconstruction of the potential @xmath130 from eqs .  ( [ f])([5.5 ] ) leads to serious difficulties .
their nature depends on the magnitude of the curvature parameter @xmath371 of parabola ( [ generica ] ) .
string - inspired models @xcite suggest that @xmath371 may be as large as 10 , or even 40 . with such large values ( and assuming non - vanishing @xmath308 ) , our numerical integrations of eqs .
( [ f])([5.5 ] ) give concave potentials @xmath130 , unbounded from below .
this corresponds to unstable theories , and thereby to extremely fine - tuned initial conditions : changing slightly the derivative of the scalar field , @xmath327 , at high redshifts would _ a priori _ yield a totally different universe at present .
this result tells us that this kind of models can not be consistent over a wide redshift interval with the exact form of @xmath20 we chose in ( [ h2gr ] ) , unless the parameter @xmath308 is extremely small . actually , this is just another way to present the results of refs .
@xcite : since they predict that @xmath308 should be almost vanishing at present , assuming a significant non - zero value implies that the theory is unnatural .    to obtain convex - shaped potentials @xmath130 ( _ i.e. _ , stable theories ) while still assuming a non - vanishing @xmath308 , we typically need values of @xmath381
. however , the reconstructed potentials always exhibit sudden changes of their slope .
basically , they reproduce a cosmological constant over a finite interval around @xmath194 ( _ i.e. _ , around @xmath269 ) , and become rapidly divergent beyond a critical value of the scalar field ( depending on @xmath371 ) . therefore , as in subsection vi.b above
, we find that such scalar - tensor models can reproduce ( [ h2gr ] ) only if they involve a cosmological constant , whose energy contribution is close to the parameter @xmath258 entering @xmath20 , and if the scalar field has a negligible enough influence .
in other words , such models would not explain the small but nonzero value of the observed cosmological constant by a `` quintessence '' mechanism , and would not be more natural than merely assuming the existence of @xmath4 .
the above results are significantly changed if we take into account the possible spatial curvature of the universe .
indeed , smoother potentials @xmath130 are obtained for closed universes ( @xmath295 ) , and the present value of the cosmological constant thus becomes more `` natural '' .    to illustrate this feature ,
let us consider the case of a minimally coupled scalar field ( as in @xcite ) , corresponding to @xmath382 in eq .
( [ generica ] ) . for an open universe ( @xmath297 )
, we find from eq .
( [ 5.5 ] ) that the scalar field would need to carry negative energy to reproduce ( [ h2gr ] ) . on the other hand , for a closed universe ( @xmath295 )
, one can derive analytically the parametric form of the potential @xmath130 .
it can be expressed in terms of the hypergeometric function @xmath383 ( solution of the differential equation @xmath384{\cal f } ' - ab { \cal f } = 0 $ ] ) :    @xmath385    [ parametricv ]    where as before @xmath386 .
if @xmath165 is very small , we recover that @xmath130 exhibits a sudden change of slope , as was obtained above in the flat case .
this is illustrated by the left panel of fig .
6 . on the contrary , if @xmath165 is large enough , the same analytical expression ( [ parametricv ] ) gives nice regular potentials , like the one displayed in the right panel of fig .
this reconstructed @xmath130 , as well as those obtained numerically for weakly varying @xmath317 , eq .  ( [ generica ] ) , are natural in the sense that they can be approximated by the exponential of simple polynomials in @xmath51 . in that case , the observed value of the cosmological constant does not appear as a mere parameter introduced by hand in the lagrangian , but corresponds basically to the present value of @xmath387 . it should be noted that a value as large as @xmath388 is not excluded by the latest boomerang data , though it would be problematic in the framework of the inflationary paradigm .    in conclusion ,
the existence of non - singular solutions over a long period of time is again the constraining input .
a non - minimally coupled scalar field is essentially incompatible with ( [ h2gr ] ) over a wide redshift interval , unless the scalar field is frozen at a minimum of @xmath317 ( consistently with @xcite ) .
if future experiments provide a hubble diagram in accordance with ( [ h2gr ] ) and also give a very small value for @xmath169 , it will be possible to conclude that scalar - tensor theories ( either non - minimally or minimally coupled ) _ can not _ explain in a natural way the existence of a cosmological constant .
on the other hand , if the universe is closed and @xmath165 large enough , a `` quintessence '' mechanism in a scalar - tensor theory seems more natural than a mere cosmological constant .
the above conclusions can be confirmed by starting from a given @xmath178 ( or @xmath325 ) , rather than @xmath33 ( or @xmath188 ) .
we consider here `` scaling solutions '' , _ i.e. _ , we assume that these functions behave as some power of the scale factor @xmath389 .
one may for instance write @xmath390 , with @xmath391 . as before ,
our aim is to reconstruct a regular potential @xmath130 from the knowledge of @xmath20 , assumed to be of the form ( [ h2gr])-([omegavalues ] ) .
the strongest constraint on this class of theories is imposed by the solar - system bound ( [ df0 ] ) . indeed , using the definition ( [ alpha ] ) for @xmath136 , one can also write it as @xmath392 , and eq .
( [ 5.5 ] ) evaluated at @xmath167 then yields the following second - order equation for @xmath393  : @xmath394 note that this equation does not depend on the full form of eq .
( [ h2gr ] ) , but only on its first derivative at @xmath269 , _
i.e. _ , on the deceleration parameter @xmath395
. the constraints on @xmath393 derived below are thus valid as soon as @xmath396 is of order @xmath397 , consistently with the estimated value ( 6.2 ) for @xmath160 .    in the case of a spatially flat universe ( @xmath164 ) ,
( [ eq.p ] ) gives immediately @xmath398 , so that the solar - system bound ( [ df0 ] ) imposes @xmath399 .
therefore , the scalar field needs to be almost minimally coupled .
if @xmath393 vanishes identically , we recover as before the trivial solution of gr plus a cosmological constant , together with an unexcited minimally - coupled scalar field . on the other hand ,
if @xmath393 does not vanish , one finds that the scalar field needs to carry negative energy beyond @xmath400 . even without trying to reconstruct the potential @xmath130
, one can thus conclude that such scaling solutions would be ruled out by the observation of a @xmath20 of the form ( [ h2gr ] ) up to @xmath262 .
paradoxically , this result is valid even for an infinitesimal ( but nonzero ) value of @xmath393 .
indeed , there exists a discontinuity between the case of a strictly constant @xmath96 and that of a scaling solution @xmath401 . at first order in @xmath393 , and still assuming @xmath164 , one can write eq .
( [ 5.5 ] ) as @xmath402 this equation confirms that @xmath403 when @xmath404 , and therefore that the scalar field tends towards a constant in this limit .
however , it carries positive energy ( @xmath405 ) only if @xmath406 since the right - hand side is estimated to be @xmath407 , the maximum value of @xmath19 is obtained for @xmath408 , so that the large numerical factor coming from @xmath409 in the left - hand side is compensated by the small term inside the second parentheses .
working iteratively , this maximum redshift can be better approximated by @xmath410 , and the actual numerical resolution of equality ( [ psmallpositivee ] ) for the values ( [ omegavalues ] ) gives @xmath411 .
therefore , even if @xmath393 is vanishingly small , a scaling solution @xmath401 can not be consistent with ( [ h2gr ] ) beyond this maximum redshift .
this illustrates once more that the experimental determination of @xmath20 up to @xmath262 would be more constraining that solar - system experiments for this class of theories , provided one takes into account the requirement of positive energy .
let us underline that the above value for @xmath282 is valid for a monomial @xmath412 but not for more complicated polynomial expressions .
indeed , as shown for instance in sec .
vi.b above , there do exist scalar - tensor theories consistent with ( [ h2gr ] ) up to arbitrarily large redshifts , and they do not need to be _ strictly _ equivalent to gr plus a cosmological constant ( although they must be close enough to it ) . moreover , the above maximum redshift is a consequence of the _ exact _ form for @xmath20 we chose in eq .
( [ h2gr ] ) .
a slightly different function may of course allow a positive - energy scalar field up to much higher redshifts .
it suffices that the right - hand side of eq .
( [ varphiprime ] ) be strictly positive for @xmath413 , and the case of a closed universe discussed below provides an example , since the contribution @xmath414 is then positive in eq .
( [ varphiprime ] ) .
the case of a spatially open universe ( @xmath297 ) is forbidden by eq .
( [ eq.p ] ) , unless @xmath169 is smaller than @xmath415 .
such a situation would be indistinguishable from the spatially flat case .    in a spatially closed universe ( @xmath295 ) , @xmath393 is given by the positive root of the second - order equation eq .
( [ eq.p ] ) . remembering the solar - system bound ( [ df0 ] )
, one may consider the case @xmath416 , and one gets @xmath417 .
even if one considered values of @xmath169 as large as @xmath418 , this would limit @xmath393 to @xmath419 . therefore , in this case again , solar - system constraints impose that the scalar field should be almost minimally coupled , if one looks for such scaling solutions .
the difference with the spatially flat case is that eqs .
( [ f])([varphiprimebis ] ) can now be integrated for any redshift @xmath19 ( from future infinity , @xmath420 , to arbitrarily large @xmath19 ) .
since @xmath178 needs to be almost constant , we recover solution ( [ parametricv ] ) for the potential @xmath130 . as in sec .
vi.c.1 above , we can thus conclude that such models would be consistent with ( [ h2gr ] ) over a wide redshift interval only if they are ( almost ) minimally coupled , and they would provide a natural `` quintessence '' mechanism to explain the presently observed cosmological constant only if the universe is ( marginally ) closed .
let us end this paragraph by a remark concerning scaling solutions , for which the scalar - field energy density scales like a power of @xmath389 .
as mentioned in the introduction , they have attracted a lot of attention recently . for a minimally coupled field , the possible scaling behaviors and the corresponding potentials
can be classified @xcite . as for a non -
minimally coupled field , a subclass of theories was considered in @xcite , for which @xmath421 where @xmath422 and @xmath423 are constants . since , besides these two constants , there is only one unknown function of @xmath9 , the knowledge of @xmath20 suffices to reconstruct the full microscopic lagrangian from eqs .
( [ f])([5.5 ] ) above . however , the main conclusion of ref .
@xcite can be recovered from a simple argument , without any numerical integration .
indeed , it was shown in this reference that there exists a universal behavior of these theories , depending on @xmath423 but not on the precise shape of @xmath33 . as emphasized in @xcite , this result
was obtained in the strong coupling limit , corresponding formally to @xmath424 in action ( [ s_jf ] ) .
taking into account the assumed relation ( [ ufm ] ) , the class of theories under consideration is thus defined by @xmath425\ .
\label{s_jfz0}\ ] ] if we now introduce a new scalar variable @xmath426 , we notice that @xmath9 disappears totally from the action
. no physical result can thus depend on the precise form of @xmath33 , and we recover the conclusion of @xcite .
the constant @xmath422 may also be set to 1 by a change of length units , and this class of theories is thus parametrized by the single real number @xmath423 .
any physical prediction must therefore depend only on @xmath423 .
in this work we have investigated the constraints that arise from the _ experimental _ knowledge of the luminosity distance in function of the redshift up to @xmath262 , corresponding to @xmath20 given by ( [ h2gr])-([omegavalues ] ) .
in particular , our universe is then presently accelerating and we have studied the viability of subclasses of scalar - tensor theories of gravity .
we have shown that the subclass of models in which the scalar partner @xmath9 of the graviton has no potential at all , and which satisfy the present - day existing constraints , are inevitably ruled out if an expansion of the form of eq .
( [ h2gr ] ) holds even for a redshift interval as tiny as @xmath427 ( see the precise numbers in section vi ) .
we see that these theories become pathological in the form of a vanishing @xmath96 , already at such low redshifts for which @xmath20 will be experimentally accessible in the near future @xcite .
hence we show that a cosmological observation of the background evolution according to eq .
( [ h2gr ] ) in the `` recent '' epoch will be enough to rule out such models .
[ on the other hand , future observations might provide a @xmath20 which confirms the existence of a scalar partner to the graviton and rule out pure gr ! ] the main reason why we obtained so constraining results is that we took into account the mathematical consistency of the theory , _
i.e. _ , the fact that it should contain only positive - energy excitations to be well behaved .
this requirement severely restricts the class of viable models .
a non - flat universe can alleviate in some cases the tight constraints we found .
however , the latest cmb data released by boomerang and maxima @xcite favor a flat universe ( in accordance with the inflationary paradigm ) , and only a marginally closed universe is still allowed by the location of the first acoustic ( doppler ) peak at @xmath428 , while an open universe is more unlikely .
the most impressive conclusion is that future cosmological observations may prove to be more constraining for massless scalar - tensor theories than solar - system and binary - pulsar tests . indeed , even if the determination of the luminosity distance @xmath11 will not reach very quickly the impressive accuracy obtained in the solar system or with binary pulsars , it will nevertheless give access to the full coupling function @xmath33 in action ( [ s_jf ] ) , or @xmath188 in the einstein - frame rewriting ( [ s_ef ] ) , whereas only its first two derivatives are presently probed .
we thank our collaborators a.  starobinsky and b.  boisseau for discussions .
one of us ( gef ) also wishes to thank a.  riazuelo for informative discussions .
centre de physique thorique is unit propre de recherche 7061 .
s.  perlmutter , g.  aldering , m.  della valle _
_ , nature * 391 * , 51 ( 1998 ) ; s.  perlmutter , g.  aldering , g.  goldhaber _ et al .
_ , astroph . j. , * 517 * , 565 ( 1999 ) .
p.m.  garnavich , r.p .
kirshner , p.  challis _ et al .
_ , astrophys
. j. lett . * 493 * , l53 ( 1998 ) ; a.g .
riess , a.v .
philipenko , p.  challis _ et al .
. j. * 116 * , 1009 ( 1998 ) .
bernardis _ et al .
_ , nature * 404 * , 955 ( 2000 ) . s.  hanany _ et al .
_ , astro - ph/0005123 ; a.  balbi _ et al .
_ , astro - ph/0005124 .
bahcall , j.p .
ostriker , s.  perlmutter , and p.j .
steinhardt , science , * 284 * , 1481 ( 1999 ) .
v.  sahni and a.a .
starobinsky , ijmpd , to appear ( 2000 ) ; astro - ph/9904398 . m.  chevallier and d.  polarski , submitted , gr - qc/0009008 r.r .
caldwell , r.  dave , and p.j .
steinhardt , phys .
rev . lett . * 80 * , 1582 ( 1998 ) .
b.  ratra and p.j.e .
peebles , phys .
d * 37 * , 3406 ( 1988 ) .
ferreira and m.  joyce , phys .
d * 58 * , 023503 ( 1998 ) ; phys .
* 79 * , 4740 ( 1998 ) .
i.  zlatev , l.  wang and p.j .
steinhardt , phys .
lett . * 82 * , 896 ( 1999 ) .
a.d .  liddle and r.j .
scherrer , phys . rev .
d * 59 * , 0223509 , ( 1998 ) . c.  kolda and d.h .
lyth , phys .
b * 458 * , 197 ( 1999 ) .
starobinsky , jetp lett . * 68 * , 757 ( 1998 ) ; gravit . & cosmol . * 4 * ( suppl . ) , 88 ( 1998 ) .
d.  huterer and m.s .
turner , phys .
d * 60 * , 081301 ( 1999 ) ; t.  nakamura and t.  chiba , mon . not .
ast . soc .
* 306 * , 696 ( 1999 ) .
saini , s.  raychaudhury , v.  sahni , and a.a .
starobinsky , astro - ph/9910231 .
caldwell , astro - ph/9908168 .
will , _ theory and experiment in gravitational physics _ , ( cambridge university press , cambridge , england , 1993 ) .
t.  damour and g.  esposito - farse , class .
* 9 * , 2093 ( 1992 ) .
t.  damour and g.  esposito - farse , phys .
d * 53 * , 5541 ( 1996 ) . t.  damour and g.  esposito - farse , phys .
d * 54 * , 1474 ( 1996 ) .
.  will , phys .
d * 50 * , 6058 ( 1994 ) .
t.  damour and g.  esposito - farse , phys .
d * 58 * , 042001 ( 1998 ) .
m.  gasperini and g.  veneziano , astropart . phys .
* 1 * ( 1993 ) 317 ; mod . phys . lett . a*8 * ( 1993 ) 3701
. t.  damour and k.  nordtvedt , phys .
lett . * 70 * , 2217 ( 1993 ) ; phys .
d * 48 * , 3436 ( 1993 ) . d.i .
santiago , d.  kalligas , and r.v .
wagoner , phys .
d * 58 * , 124005 ( 1998 ) .
t.  damour and b.  pichon , phys .
d * 59 * , 123502 ( 1999 ) .
a.  navarro , a.  serna , and j.m .
alimi , phys .
d * 59 * , 124015 ( 1999 ) .
uzan , phys .
d * 59 * , 123510 ( 1999 ) .
t.  chiba , phys .
d * 60 * , 083508 ( 1999 ) .
x.  chen and m.  kamionkowski , phys .
d * 60 * , 104036 ( 1999 ) .
f.  perrotta , c.  baccigalupi , and s.  matarrese , astro - ph/9906066 .
holden and d.  wands , gr - qc/9908026 .
n.  bartolo and m.  pietroni , hep - ph/9908521 .
e.  gaztanaga and j.a .
lobo , astro - ph/0003129 .
b.  boisseau , g.  esposito - farse , d.  polarski , and a.a .
starobinsky , phys .
* 85 * ( 2000 ) 2236 .
misner , k.s .
thorne , j.a .
wheeler , _ gravitation _ ( san francisco , freeman , 1973 ) .
bergmann , int .
* 1 * , 25 ( 1968 ) ; k.  nordtvedt , astrophys . j. * 161 * , 1059 ( 1970 ) ; r.  wagoner , phys . rev .
d * 1 * , 3209 ( 1970 ) .
_ , bull .
phys .  soc .
, abstract # k  11.05 ( 1997 ) .
dickey _ et al .
_ , science * 265 * , 482 ( 1994 ) ; j.g .
williams , x.x .
newhall , and j.o .
dickey , phys .
d * 53 * , 6730 ( 1996 ) .
p.  fiziev , gr - qc/9911037 .
eddington , _ the mathematical theory of relativity _ , ( cambridge university press , cambridge , 1923 ) .
b.m .  barker , astrophys .
j. * 219 * , 5 ( 1978 ) .
bardeen , phys .
d * 22 * , 1882 ( 1980 ) .
l.  amendola , phys .
d * 60 * , 043501 ( 1999 ) .
j.  weller and a.  albrecht , astro - ph/0008314 ; i.  maor , r.  brustein and p.j .
steinhardt , astro - ph/0007297 ; t.  chiba and t.  nakamura , astro - ph/0008175 .
t.  damour and a.  polyakov , nucl .
* b423 * , 532 ( 1994 ) ; gen .
* 26 * , 1171 ( 1994 ) .