the accelerated expansion of the universe has now been confirmed by several independent observations including those of high redshift type ia supernovae , and the cosmic microwave background ( cmb ) combined with the large scale structure of the universe @xcite .
another way of presenting this kinematic property of the universe is to postulate the existence of a new entity  dark energy ( de ) .
this latter statement is dynamical in nature and therefore requires some assumptions to be made about the form of gravitational field equations governing the evolution of the ( observed part of the ) universe .
although observationally well established , no single theoretical model provides an entirely compelling framework within which cosmic acceleration or de can be understood .
indeed , the very many models of de existing in the literature illustrate that its nature is still very much an enigma . at present , all existing observational data are in agreement with the simplest possibility of de being a cosmological constant @xmath2 with @xmath3gev@xmath4 ( inside @xmath5 error bars in the worst case ) .
is used throughout the paper . ]
this case is internally self - consistent and non - contradictory .
the extreme smallness of the cosmological constant expressed in either planck , or even atomic units means only that its origin is _ not _ related to strong , electromagnetic and weak interactions ( in particular , to the problem of the energy density of their vacuum fluctuations ) . where @xmath6 is some characteristic neutrino rest - mass ( the lightest one ? ) suggests that vacuum energy of the interaction responsible for non - zero neutrino rest - masses may be relevant for a non - zero @xmath2 . ]
although in this case de reduces to only a single fundamental constant we still have no derivation from any underlying quantum field theory for its small value .    within this context it is but natural that other possibilities admitting a ( slightly ) variable dark energy have also been actively studied by the scientific community in recent years .
moreover , it is interesting that properties of the currently observed ` late de ' are _ qualitatively _ similar to those of an ` early de ' which is believed to have given rise to accelerated expansion ( inflation ) in the early universe .
however in the case of the latter , there are sufficient grounds to support the view that ` early de ' was unstable and , thus , more complicated than a cosmological constant .
so , it is natural to conjecture by analogy that the same might also be true of ` late de ' .
de models proposed to account for the present cosmic acceleration include : +
( i ) quiessence with @xmath7 , the @xmath2 ( @xmath8 ) is a special member of this class .
+ ( ii ) quintessence models which are inspired by the simplest class of inflationary models of the early universe and employ a scalar field rolling down a potential @xmath9 to achieve late - time acceleration .
quintessence potentials with @xmath10 have the attractive property that approaches a common evolutionary ` tracker path ' from a wide range of initial conditions .
+ ( iii ) the chaplygin gas model ( cg ) has the equation of state @xmath11 and evolves as @xmath12 where @xmath13 is the redshift , @xmath14 .
it therefore behaves like dark matter at early times ( @xmath15 ) and like the at late times .
cg appears to be the simplest model attempting to unify de and non - baryonic cold dark matter .
+ ( iv ) ` phantom ' de ( @xmath16 ) .
+ ( v ) oscillating de . + ( vi ) models with interactions between de and dark matter .
+ ( vii ) scalar - tensor de models .
+ ( viii ) modified gravity de models in which the gravitational lagrangian is changed from @xmath17 to @xmath18 where @xmath17 is the scalar curvature and @xmath19 is an arbitrary function .
+ ( ix ) dark energy driven by quantum effects . +
( x ) higher dimensional ` braneworld ' models in which acceleration is caused by the leakage of gravity into extra dimensions .
+ ( xi ) holographic dark energy , etc .
+ see the reviews @xcite for an exhaustive list of models and references .
however , none of these models leads to the reduction of the number of fundamental constants ( parameters in the microscopic lagrangian ) as compared to standard @xmath2cdm .
in other words , at the current state - of - the - art , de requires at least one new parameter whose value is set from observations .
models with variable de can be broadly divided into two main classes :    , in these models de is the energy density of some new , very weakly interacting physical field .
( otherwise dubbed modified gravity models ) . in these models
the gravity equations do not coincide with those of einsteinian general relativity .
however , it is usually possible to re - write the new equations in the conventional _ einsteinian _ form by transferring all additional terms from the l.h.s . into the r.h.s . of the einstein equations and referring to them as an effective energy - momentum tensor of de ( see sec . 2 below ) .
( precisely this happened to the cosmological constant which originally appeared in the l.h.s .
of einstein s field equations but is now felt by many to constitute an effective matter term such as vacuum energy , etc . )
another category , which is even more important from the observational point of view , arises in response to the question whether or not the description of de requires a new field degree of freedom (= a new kind of matter ) .
if the answer is in the affirmative then de may be considered as being ` induced ' by other kinds of matter .
all physical de models and many geometrical ones belong to this category but there do exist geometrical de models which do not ( for instance the @xmath18 model with the palatini variation of its action ) .    faced with the increasing proliferation of de models each with its own physical motivations and assumptions , a concerned cosmologist can proceed in either of two ways :    \(i ) test _ every single _ model against observations .
\(ii ) try and ascertain properties of dark energy in a _ model independent manner_.    in this article , we proceed along route ( ii ) and attempt to review both the successes as well as difficulties faced by methods attempting to reconstruct the properties of dark energy directly from observations in a model independent manner .
before attempting to determine its properties , we first need to provide a definition of _ dark energy_. a traditional approach is to use the einstein form of the gravitational field equations [ einst ] r _ - 12g_r = 8 g ( _ a t_^(a ) + t_^de ) as providing a _ definition _ of the effective energy - momentum tensor @xmath20 of de . here , the summation over @xmath21 in the r.h.s .
includes all types of matter known from laboratory experiments ( protons , neutrons , photons , neutrinos , etc . ) as well as non - relativistic non - baryonic cold dark matter ( whose energy - momentum tensor is dust - like in the first approximation @xmath22 ) .
@xmath23 is the present value of newton s gravitational constant .
this definition has a number of advantages : + ( i ) it is simple , well - defined and self - consistent ; + ( ii ) it treats physical and geometrical de on an equal footing ; + ( iii ) in the absence of direct physical interaction between de with known forms of matter or with cold dark matter , the de energy - momentum tensor is conserved : @xmath24 .
+ one should stress here that using ( [ einst ] ) we automatically ascribe terms describing ( possible ) gravitational interactions between de and non - relativistic matter , as well as the matter energy - momentum tensor multiplied by a change in the effective gravitational constant , to the de energy - momentum tensor .
the latter possibility arises , for instance , in scalar - tensor de models .    when applied to a homogeneous and isotropic friedmann - robertson - walker ( frw ) cosmological model , eqs .
( [ einst ] ) reduce to two algebraically independent equations : h^2 & = & ( _ a _ a + _ de ) -  , [ eq : acc ] + & = & -(_a ( _ a + 3p_a ) + _ de + 3p_de ) [ eq : acc1 ] where @xmath25 is the frw scale factor and the hubble parameter @xmath26 .    at late times when radiation may be neglected , one gets [ eq : energy ] _ de = ( 1 - _ m ) where we have omitted the contribution from the curvature term in ( [ eq : acc ] ) for simplicity .
@xmath27 is the total density of non - relativistic matter in terms of its critical value .
similarly , the expression for the deceleration parameter q - a / ah^2 =   x - 1 ,   x = 1 + z [ eq : decel1 ] ( prime implies differentiation with respect to @xmath28 ) takes the form : p_de = ( q - )  .
[ eq : pressure ]    dividing ( [ eq : pressure ] ) by ( [ eq : energy ] ) we get the following expression for the _ effective _ equation of state ( eos ) of dark energy ( @xmath29 ) : [ eq : state ] w(x ) = 2 q(x ) - 1 3 ( 1 - _ m(x ) ) . for physical de ,
the eos makes physical sense .
however , this is not so for geometrical de for which the acceleration of the universe is caused by the fact that the field equations describing gravity are not einsteinian .
as a specific example of geometrical de consider the widely studied dgp braneworld @xcite for which [ eq : dgp ] h = +  , where @xmath30 is a new length scale and @xmath31 and @xmath32 refer respectively to the four and five dimensional planck mass . the acceleration of the universe in this model arises not because of the presence of de but due to the fact that gravity becomes five dimensional on length scales @xmath33 .
the contrast between ( [ eq : acc ] ) and ( [ eq : dgp ] ) makes it abundantly clear that for models such as dgp the eos in ( [ eq : state ] ) is an _ effective _ quantity ( @xmath34 ) , which may still be useful for descriptive purposes but which no longer represents any fundamental physical property of an accelerating universe .
indeed , instances are known when @xmath35 even when matter itself satisfies the weak energy condition @xmath36 @xcite .
it may be instructive to note that , for geometrical de , @xmath1 may show pathological behaviour in certain cases , such as the presence of poles at finite values of redshift , @xmath37 , even though the underlying cosmological model is completely well behaved ( see for instance @xcite ) .
for such models , the deceleration parameter @xmath38 and other geometrical parameters prove to be more robust for determining de properties than the eos .
we shall return to these important issues in section [ sec : statefinder ] .
observational tests of de rely on an accurate measurement of at least one of the following quantities :    1 .   the luminosity distance [ eq : lumdis ] = _ 0^z  .
the angular size distance [ eq : angdis ] d_a(z ) = _ 0^z  .
the coordinate distance [ eq : coordis ] r(z ) = _ 0^z  .    here
, we have assumed that the universe is spatially flat for simplicity .
note that the distance duality relation @xmath39 which follows from ( [ eq : lumdis ] ) and ( [ eq : angdis ] ) is valid only for metric theories of gravity @xcite .
its violation ( if observed ) could , therefore , be used to probe alternative theories of gravity .
existing data , however , appear to support it ( within observational errors ) @xcite .    in all of the above expressions the value of the hubble parameter
can be ` reconstructed ' through a relation similar to the one given below for the luminosity distance @xcite : [ eq : h ] h(z)=^-1  .
differentiating a second time allows one to reconstruct the equation of state of de ( [ eq : state ] ) .
equations ( [ eq : state ] ) and ( [ eq : h ] ) immediately inform us that @xmath1 will be a noisier quantity than @xmath0 since two successive differentiations are needed for the reconstruction @xmath40 while a single differentiation suffices for @xmath41 .
this has led several people to suggest @xmath0 ( or @xmath42 ) as being better suited for providing a model independent description of properties of de .
another important difference between @xmath0 defined in ( [ eq : h ] ) and @xmath1 in ( [ eq : state ] ) is that the former is independent of the value of the matter density parameter @xmath27 while the latter is not . as a result
, uncertainties in the current value of @xmath27 affect the reconstruction of the eos far more profoundly than they do @xmath0 .
we shall return to this issue in section [ sec : obstacle ] .    knowing @xmath0 ( either through ( [ eq : h ] ) or using corresponding relations for other observational tests ) , allows us to extend cosmological reconstruction to other important physical properties of the universe including :    * its age [ eq : age ] t(z ) = _ z^ . *
the deceleration parameter ( [ eq : decel1 ] ) and the equation of state ( [ eq : state ] ) ; * the electron - scattering optical depth to a redshift @xmath43 ( z_reion ) = c_0^z_reion , [ eq : reion ] where @xmath44 is the electron density and @xmath45 is the thomson cross - section describing scattering between electrons and cmb photons . * the product @xmath46 , which plays a key role in the alcock  paczynski anisotropy test @xcite . * the product @xmath47 , which is used in the volume - redshift test @xcite . * the parameter @xmath48 associated with the determination of the baryon acoustic peak @xcite a =  _ 0^z_1  ^2/3 =  0.469 0.017  , where @xmath49 and @xmath50 is the redshift at which the acoustic scale has been measured in the redshift sample . * the ` shift ' parameter @xmath17 associated with the cmb @xcite + r = _ 0^z_lss = 1.7 0.03 .    for quintessence , it is additionally possible to reconstruct its potential , since the einstein equations @xmath51 can be rewritten as ( see , for instance , @xcite ) : @xmath52 integrating ( [ eqn : phidot ] ) , one determines @xmath53 to within an additive constant .
the inversion @xmath54 followed by substitution into ( [ eqn : vzed ] ) allows us to reconstruct @xmath9 , since @xmath0 and its first derivative can be determined from observations using ( [ eq : h ] ) .
the presence of @xmath55 in ( [ eqn : vzed ] ) and ( [ eqn : phidot ] ) implies that the value of the matter density must be known rather precisely for an accurate reconstruction of @xmath9 .
in several important cases ( [ eqn : vzed ] ) and ( [ eqn : phidot ] ) have led to a closed form expression for @xmath9 , for instance :    * de with a constant equation of state @xmath56 is described by @xcite + v()=3h_0 ^ 2(1-w)(1-_m0)^1/|w| 16g_m0^ ^-2(|w| ( -_0+_1 ) )  , [ quequi ] where = ,  _
0=(t_0 ) ,  _ 1=  . consequently
, a universe consisting of such a scalar field will have expansion properties which closely mimic a different kind of de , e.g. , a tangled network of cosmic strings ( @xmath57 ) or domain walls ( @xmath58 ) .
* the chaplygin gas model which unifies dark matter and de and has @xmath59 can be described by the scalar field potential @xcite ( in the absence of additional cold dark matter and neglecting baryons and photons ) : + v ( ) = ( ( 2 ) + )  .
note that the behaviour of the chaplygin gas may also be modelled completely differently using a scalar field with the born - infeld kinetic term @xcite ( a specific type of k - essence ) .
this illustrates once more that , even for physical de , the equation of state @xmath1 does not uniquely define an underlying field - theoretical model is quite different from the quintessence lagrangian @xmath60 . ]
( see also @xcite ) .
from the fact that the left - hand side of eq . (
[ eqn : phidot ] ) is always non - negative , follows an important restriction on the expansion law for the universe which must be satisfied before attempts to reconstruct the potential are made , namely 3h_0 ^ 2(1+z)^2  .
[ ineq ] equation ( [ ineq ] ) is simply a restatement of the weak energy condition @xmath61 . integrating ( [ ineq ] ) ,
we get the relation h^2(z)h_0 ^ 2(1+(1+z)^3 ) [ int - ineq ] which is easier to verify observationally .    in the case of physical de ,
with its implied minimal coupling to gravity , the expansion history @xmath0 can also be reconstructed from the growth rate of inhomogeneous density perturbations in the non - relativistic matter component on scales significantly less than the hubble radius @xmath62 , ( provided perturbations in the de component can be neglected and the effective gravitational constant does not change with time ) @xcite .
indeed , the linearized perturbation equation and other quantities referring to the evolution of a frw background that , in turn , can be used to determine @xmath0 @xcite . ] + 2h - 4g_m= 0
[ eq : recon3 ] can easily be inverted , with the result : h(z)= h(0 ) - 3_0^z dz^. [ eq : recon7 ] furthermore , an interesting relationship follows between the current value of the matter density @xmath55 and @xmath63 : = ^2(0)(3_0^dz)^-1 , [ eq : recon6 ] which could provide a consistency check on direct observational determinations of @xmath55 . at sufficiently low redshifts
@xmath64 , the advent of deep redshift surveys probing large scale structure may help in reconstructing @xmath0 using ( [ eq : recon7 ] ) . although it is unlikely that the value of @xmath63 for @xmath65 will be reliably known in the near future , the fact that the universe is expected to become matter dominated at fairly low redshifts , @xmath66 at @xmath67 , allows us to use the spatially flat matter dominated solution @xmath68 to extrapolate to higher redshifts and thereby evaluate ( [ eq : recon6 ] ) .    as a result , for physical de
we now have two independent methods to reconstruct the same quantity @xmath69 : ( i ) through quantities like @xmath70 referring to an unperturbed frw background , and ( ii ) through @xmath63 describing the growth of perturbations on small scales .
it therefore follows that one should be able to reconstruct the density perturbation @xmath63 from observations of the luminosity distance .
consider the quantity @xmath71 which can also be written as the difference of conformal times : @xmath72 .
for our purpose it will be useful to invert the quantity @xmath73 ( determined from observations ) to @xmath74 .
next , we again assume that de does not cluster and rewrite ( [ eq : recon3 ] ) in the following integral form : ( e(z ) ) = ( 0 ) + (0)_0^e ( 1+z(e))de + _ 0^e ( 1+z(e_1))de_1_0^e_1 ( e_2)de_2 . [ eq : recondelta ] here @xmath75 since @xmath76 for @xmath77 .
clearly ( [ eq : recondelta ] ) can be solved iteratively if the value of @xmath78 is known .
the value of @xmath79 does not really matter since we are only interested in the ratio @xmath80 .
the requirement @xmath81 leads to the integral condition : (0)= - ( 3/2)_m_0^e_lss(e)de  , and , as in the case of ( [ eq : recon6 ] ) , one needs to postulate some reasonable interpolating behaviour for @xmath82 at large redshifts for which no supernova data exist at present . since @xmath83 , the value of this quantity should be possible to determine from the cmb .
both ( [ eq : recon7 ] ) and ( [ eq : recondelta ] ) can be used as consistency checks for physical de which , by assumption , is minimally coupled to gravity . for geometrical de on the other hand ,
the linearized perturbation equation is usually modified from its conventional form ( [ eq : recon3 ] ) , see for instance @xcite for scalar - tensor gravity or @xcite for the dgp model .
is expected to be different from the reconstructed expression ( [ eq : recondelta ] ) provides an interesting possibility to distinguish between physical and geometrical models of de . ] as a result , de reconstruction becomes more complicated . in particular , a generic scalar - tensor model of de
depends on two arbitrary functions @xmath9 and @xmath84  the scalar field potential and its coupling to gravity respectively . in this case one needs information from _ both _ @xmath70 and @xmath63 in order to reconstruct the de model unambiguously @xcite .    to summarize ,
we have shown in this section that the reconstruction of the expansion rate @xmath0 and other de properties from observational data is well defined and unambiguous from the mathematical point of view . in practice , however , the situation is much more difficult since all observational functions such as @xmath70 and @xmath63 are noisy and known only at discrete values of the redshift @xmath85 ( associated with the redshifts of @xmath86 supernovae in the case of @xmath70 ) .
thus , it is impossible to directly differentiate them with respect to redshift as formulae ( [ eq : h ] ) and ( [ eq : recon7 ] ) require .
add to this the fact that the dispersion in the luminosity distance is not expected to get significantly better than @xcite @xmath87 , and one is confronted with a noisy quantity @xmath88 sampled at a set of discrete intervals @xmath89
. therefore , to convert from @xmath88 to the function @xmath0 defined at _ all _ redshift values within the interval ( say ) @xmath90 using ( [ eq : h ] ) clearly requires a crucial additional step : some sort of _ smoothing _ procedure .
this is usually accomplished using either parametric or non - parametric reconstruction which we turn to next .
this approach is based on the assumption that the quantities @xmath91 vary ` sufficiently slowly ' with redshift and can therefore be approximated by a fitting formula ( ansatz ) which relies on a small number of free parameters @xmath92 .
the ansatz for @xmath93 is compared against observations , and the values of the free parameters @xmath94 are determined using a minimization procedure ( usually maximum likelyhood ) .
quite clearly a successful fitting function should , in principle , be able to faithfully reproduce the properties of an entire class of de models .
different parameterizations have been suggested for : @xmath95 @xcite , @xmath0 @xcite , @xmath1 @xcite and @xmath96 @xcite .
we attempt to summarize some of these approaches below .    1 .
* fitting functions to the luminosity distance @xmath95 * + an approximation to a function can easily be generated by expanding it in a taylor series about some redshift @xmath97 .
applied to the luminosity distance ( with @xmath98 ) , this method gives the following fitting function @xcite = _
i=1^n a_iz^i .
[ eq : poly ] unfortunately , it appears that in order to accurately determine quantities of interest such as @xmath0 and @xmath1 one must take a large number of terms in ( [ eq : poly ] ) which greatly increases the errors of reconstruction @xcite and reduces the reliability of this ansatz . + a more versatile pad - type ansatz
was suggested in @xcite = 2  ,   x = 1+z  , [ eqn : star ] which is able to _ exactly _ reproduce the results both for cdm ( @xmath99 ) and the steady state model ( @xmath100 ) . +
another accurate fit with a greater number of free parameters is @xcite & & = ( 1)-(y )  , + & & ( y)=2 ^ -1/8  ,   y = 1/  .
[ eq : fit_chiba ] 2 .
* fitting functions to the de density * + the dark energy density can be written as a truncated taylor series polynomial in @xmath101 , @xmath102 .
this leads to the following ansatz for the hubble parameter @xcite [ eq : poly_fit ] h(x ) = h_0x^3 + a_1 + a_2x + a_3 x^2^ , [ eq : taylor ] which , when substituted in the expression for the luminosity distance ( [ eq : lumdis ] ) , yields = _ 1 ^ 1+z .
[ eq : taylor1 ] this fitting function gives exact results for the cosmological constant @xmath8 ( @xmath103 ) as well as for quiessence with @xmath104 ( @xmath105 ) and @xmath106 ( @xmath107 ) .
the presence of the term @xmath108 in ( [ eq : poly_fit ] ) ensures that the ansatz correctly reproduces the matter dominated epoch at early times ( @xmath67 ) . for quintessence models as well as the chaplygin gas , the luminosity distance @xmath70 in ( [ eq : taylor1 ] )
can be determined to an accuracy of better than @xmath109 if @xmath110 @xcite .
* fitting functions to the equation of state * + in this approach one assumes that the de equation of state @xmath1 is an unknown variable whose behaviour is ` guessed ' by means of a suitable fitting function @xmath111 . since h^2(z ) & = & h_0 ^ 2 ( 1+z)^3 + ^2  , +   + & = & ( 1- )  , [ eq : hubble_recon ] the values of @xmath94 can be determined by comparing the luminosity distance ( [ eq : lumdis ] ) or @xmath112 , against observations .
+ several possible fits for @xmath1 have been suggested in the literature .
perhaps the simplest is the taylor expansion @xcite w(z ) = _
i=0^n w_iz^i , [ eq : wpol ] which , for @xmath113 , gives results significantly better than the taylor expansion for the luminosity distance ( [ eq : poly ] ) .
the simple two parameter representation @xmath114 is however of limited utility since it is only valid for @xmath115 .
a considerably more versatile four parameter ansatz has been suggested in @xcite . +
a popular two parameter fit was suggested by chevallier and polarski @xcite and by linder @xcite w(a ) = w_0+w_1(1-a ) = w_0 + w_1  , [ eq : cpl ] the associated luminosity distance can be obtained by substitution into ( [ eq : hubble_recon ] ) and ( [ eq : lumdis ] ) . a more general form for this fit is w(a ) = w_p + ( a_p - a)w_a  , where @xmath116 is the value of the scale factor at the ` sweet spot ' where the equation of state @xmath117 is most tightly constrained .
the value of @xmath116 usually depends upon the data set being used @xcite .
@xmath118 \epsfxsize=2.7 in \epsffile{h_exp_w_g.epsi } & \epsfxsize=3.3 in \epsffile{huterer3.ps } \\
\end{array}$ ]    as mentioned earlier , for an ansatz to be regarded as being successful it should embrace within its fold the behaviour of a reasonably wide class of de models .
most of the ansatz s discussed above successfully accommodate de whose equation of state evolves moderately with redshift .
it is quite clear that these simple fits can not be used to rule out ( using observations ) models with rapidly evolving @xmath1 .
one reason for this is that applying an ansatz such as ( [ eq : taylor ] ) or ( [ eq : cpl ] ) to sne data is equivalent to smoothing the evolution of the universe over a redshift interval @xmath119 where @xmath86 is the number of free parameters in the ansatz ( the concrete form of an ansatz defines the form of a smoothing redshift function ) .
clearly , an implicit smoothing such as this will cause rapid transitions in dark energy to disappear not because of disagreement with data but simply because of the manner in which the data have been ` massaged ' @xcite .
( an example of the disastrous results of applying the prior @xmath120 constant to models with an evolving eos is discussed in section [ sec : obstacle ] . )
models with a fast phase transition in dark energy ( @xmath121 over a narrow range of redshift @xmath122 ) have been discussed in @xcite . to accommodate such models
the following ansatz was suggested @xcite w(z)=w_i+ , [ eq : bass ] where @xmath123 is the initial equation of state at high redshifts , @xmath124 is a transition redshift at which the equation of state falls to @xmath125 and @xmath126 describes the rate of change of @xmath1 .
substitution into ( [ eq : hubble_recon ] ) and ( [ eq : lumdis ] ) gives the luminosity distance .
finally , models with oscillating de have also been discussed in the literature @xcite . to accommodate an oscillating eos the following ansatz has been suggested @xcite w ( ) = w_0 + w_1where @xmath127 are free parameters whose values must be obtained by fitting to observations .
it is encouraging to note that , in spite of some ambiguity in the form of the different fits , when applied to the same supernova data set most of them give consistent results in the range @xmath128 ( where there is sufficient data ) @xcite .
one should also note that although increasing the number of parameters usually increases the accuracy of reconstruction of the ` best fit ' , this is often accompanied by severe degeneracies which limit the utility of introducing a large number of free parameters .
in addition , extra free parameters are quite severely penalized by information criteria such the akaike information criterion @xcite and the bayesian information criterion @xcite , see also @xcite .
however , there exists a subtlety when applying these information criteria to fitting functions of @xmath129 etc . as opposed to concrete theoretical models of de . in the case of the former
, one should keep in mind the possibility that a fit with a larger number of parameters might in fact be describing the behaviour of a fundamental de model ( not yet known ) containing a smaller number of truly free parameters .
the reverse is also true , a primitive fit ( hence not penalizable ) may , if the information criteria are applied , detract our attention from a more complicated but also more fundamental explanation of a phenomenon .
let us illustrate this with an example .
the theory of gravitational instability informs us that the linear density contrast @xmath130 and the peculiar velocity field @xmath131 are related as @xcite = - [ eq : velocity ] where the function @xmath132 can be approximated as @xmath133 .
an exact calculation however reveals @xmath134 to be a more complicated function involving elliptical integrals , etc .
therefore a naive application of information criteria to ( [ eq : velocity ] ) would needlessly penalize the latter , which is correct but complicated , in favour of the simpler @xmath133 .
another popular approach the so - called principal component analysis is based on expanding @xmath1 in terms of a basis of orthogonal functions .
the form of these functions depends upon the kind of data used and its constraining capabilities @xcite .
we end this section by noting that the reconstruction approach has also been extended to scalar - tensor gravity @xcite , string inspired cosmology with the gauss - bonnet term coupled to a scalar field @xcite and dark energy with non - canonical kinetic energy terms @xcite .
other applications of cosmological reconstruction may be found in @xcite .
non - parametric smoothing usually involves directly smoothing either @xmath95 , or some other quantity appropriately binned in redshift space with some characteristic smoothing scale .
different ways of implementing this approach have been discussed in @xcite .
all reconstruction methods must deal with the fact that the data sample is usually sparse and the coverage of redshift space uneven .
consider for instance the comoving coordinate distance to a single supernova within a larger sample r_i = r(z_i ) + n_i  . here
, @xmath135 is the noise term @xmath136 , and @xmath137 is related to the distance modulus @xmath138 of the supernova by @xcite @xmath139
. direct differentiation of @xmath140 to give @xmath141 will obviously result in a very noisy quantity .
one way to tackle this is to fit @xmath142 piece - wise using a set of basis functions . following this approach , daly and djorgovsky @xcite reconstructed the deceleration parameter from the dimensionless coordinate distance @xmath143 by : ( i ) fitting @xmath144 locally within a small redshift bin by means of a second order polynomial requiring that at least 10 data points lie within each bin ; ( ii ) determining the first and second derivatives of @xmath144 from the fit coefficients ( for each bin ) and reconstructing the deceleration parameter from @xmath144 by means of the relation -q(z ) = 1 + ( 1+z )  . applying their method to a set of radio galaxies in addition to type ia sne , daly and djorgovsky @xcite found that although the cosmological constant provided a very good fit to the data , modest evolution in the dark energy density was also perfectly acceptable .
the transition from deceleration to acceleration occured at @xmath145 , with a best fit value of @xmath146 .
these results are in broad agreement with those obtained using parametric approaches @xcite .
a different approach to non - parametric reconstruction is discussed by shafieloo @xcite who generalize a smoothing ansatz widely used in the analysis of large scale structure . according to this method
, a smoothed quantity @xmath147 is constructed from a fluctuating ` raw ' quantity @xmath148 using a low pass filter @xmath19 having a smoothing scale @xmath126 d^s(*x * , ) = d(*x * ) f(|*x*-*x*| ; )  d*x *  . in large scale structure studies , @xmath149 is the density field @xcite whereas for cosmological reconstruction @xmath149 could be either of @xmath150 .
commonly used filters include : ( i ) the ` top - hat ' filter , which has a sharp cutoff @xmath151 where @xmath152 is the heaviside step function ( @xmath153 for @xmath154 , @xmath155 for @xmath156 ) and ( ii ) the gaussian filter f_g  .
when applied to snap - quality data using a gaussian filter , this method reconstructs the hubble parameter to an accuracy of @xmath157 within the redshift interval @xmath158 . the look - back time t(z)=t(z=0)-t(z)=h_0
^ -1_0^z is reconstructed to an even better accuracy of @xmath159 at @xmath160 .    despite the considerable success of these approaches , neither they , nor the parametric methods discussed earlier provide independent measurements of either @xmath0 or @xmath38 within a given redshift range .
one approach towards an uncorrelated determination of cosmological quantities is to bin the data in redshift bins of size @xmath126 with a weight ( or window ) function assigned to each bin @xcite .
each window function vanishes outside its redshift bin and the two supernovae falling on either side of a bin boundary are discarded .
this method ensures that , since two adjacent window functions have no common supernovae , the measurement of quantities within these bins will have uncorrelated error bars .
this method is quite successful and complements other reconstruction approaches in many respects .
one weakness , which practitioners of this approach have to guard against is that , since some information at the bin boundary is lost , discontinuities in the fit values at bin edges may be present @xcite .
these , in turn , could lead to unphysically large derivatives of derived quantities such as @xmath0 and @xmath1 .
results obtained by applying a variant of this method to recent cmb+sne+lss data are shown in figure [ fig : wangm ] and appear to agree with the results of the parametric approach in figure [ fig : eos ] .
other promising methods for obtaining uncorrelated estimates of cosmological evolution are discussed in @xcite .    in all of the above approaches , it is important to choose the value of the smoothing scale @xmath126 ( or bin size )
optimally , so that ( i ) a sufficient number of data points is accommodated within each bin to reduce the effects of shot noise , ( ii ) the bin size is not too large to cause excessive smoothing .
the following useful formula giving the relative error bars on @xmath0 comes to our aid @xcite [ eq : tegmark ] , where n is the total number of supernovae ( assuming an approximately uniform distribution ) and @xmath161 is the noise of the data .
clearly , a reduction in the value of @xmath126 by 3 should be compensated for by having 27 many times more supernovae if one wants to keep @xmath162 unchanged .
the situation is far worse for @xmath1 which has a @xmath163 scaling , since the data are being differentiated twice to get @xmath1 .
with respect to cosmological reconstruction , considerable caution must be exercised when setting priors on the values of cosmological parameters .
this is true both for the density parameter @xmath27 whose value is currently known to about @xmath164 accuracy and the de equation of state .
maor @xcite gave a particularly insightful example of the dangers of incorrect reconstruction by assuming a fiducial de model with @xmath165 and an evolving eos @xmath166 .
in their reconstruction , the results of which are shown in figure [ fig : maor ] , it was assumed that @xmath167 was a constant and , additionally , that @xmath168 .
these two ( incorrect ) priors led to a gross underestimation of @xmath167 and an overestimation of @xmath27 as shown in figure [ fig : maor ] .
the above example shows that the eos of dark energy can be badly reconstructed if we assume a misleading prior for @xmath169 .
since @xmath1 depends upon the value of @xmath27 through ( [ eq : state ] ) , it follows that an incorrect assumption about the value of the matter density has the potential to affect cosmological reconstruction quite significantly .
this is demonstrated in figure [ fig : arman ] in which the de eos is reconstructed for a fiducial @xmath2cdm model with @xmath170 and @xmath165 .
cosmological reconstruction is based on the ansatz ( [ eq : poly_fit ] ) in which the _ incorrect _ value @xmath171 is assumed .
this results in the reconstructed eos evolving with redshift when , in fact , no such evolution is present in the fiducial model !
as we have seen , the eos has certain blemishes  it is obtained from @xmath95 after differentiating _ twice _ and it is sensitive to the _ precise _ value of @xmath55 .
both these features hamper its reconstruction which has led several authors to propose the hubble parameter or the de density as being more appropriate for cosmological reconstruction @xcite .
carrying these ideas further , the possibility of extracting information about the equation of state from the reconstructed hubble parameter by constructing a _ weighted average _ of the equation of state was explored in @xcite .
this quantity , dubbed the @xmath169-probe , gleans information about the equation of state from the _ first derivative _ of the luminosity distance .
therefore , it is less noisy and better determined than @xmath1 .
the weighted average of the equation of state is defined as @xcite [ eq : weight1 ] 1+|w= ( 1+w(z ) ) .
an important feature of the @xmath169-probe ( @xmath172 ) is that it can be expressed in terms of the difference in dark energy density over a given redshift range [ eq : weight2 ] 1+|w(z_1,z_2 ) & = &   + & & ( )  .
+ here @xmath130 denotes the total change of a variable between integration limits and @xmath173 ( @xmath174 ) .
thus @xmath172 is easy to determine if the hubble parameter has been accurately reconstructed .
an important property of @xmath172 is that it is less sensitive to uncertainties in the value of @xmath55 than @xmath1 .
this is demonstrated in figure [ fig : margin ] where the results for @xmath172 are shown after marginalizing over the matter density .
remarkably , the value of @xmath172 for the fiducial @xmath2cdm model remains close to @xmath175 , while @xmath172 for the evolving de model shows a clear signature of evolution .
thus , small uncertainties in the value of the matter density do not appear to adversely affect the accuracy of the reconstructed @xmath169-probe .
furthermore , several excellent methods for determining @xmath176 and @xmath0 have been suggested in the literature @xcite , any of which could be used to determine @xmath172 using ( [ eq : weight2 ] ) .
cosmological observations made during the past two decades have brought about both qualitative and quantitative changes in our perception of the universe .
observational support has come for many theories ( such as inflation ) which earlier were regarded as being at the level of hypothesis . on the other hand , fundamentally new properties of our universe , such its current state of acceleration ,
have also been unravelled .
the remarkable qualitative similarity between current acceleration ( fueled by de ) and inflationary acceleration ( supported by what may be called primordial dark energy ) may help us understand in which direction cosmology may be headed .
the discovery by cobe of an approximately scale - invariant spectrum of primordial fluctuations ( @xmath177 ) was regarded by many as providing tacit support for the inflationary scenario .
however , virtually all models of inflation predict small departures from the scale invariance @xmath178 , where @xmath179 is a model dependent quantity .
strict scale invariance ( @xmath180 ) arises for a specific 2-parameter family of inflaton potentials only @xcite , which reduces to @xmath181 in the slow - roll approximation . in this context
, the recent 3-year data release of the wmap experiment suggests @xmath182 , which is in excellent agreement with predictions made by the simplest inflationary models such as chaotic inflation with @xmath183 .    the present situation concerning de resembles , in some respects , the status of inflation just after the release of the cobe data in 1992 .
the cosmological constant is in excellent agreement with data and departures ( if any ) from @xmath184 are believed to be quite small . nevertheless , since at present no fundamental theory predicting the ( small ) value of the @xmath2-term exists , it is of utmost importance that departures from @xmath184 be probed to accuracies of at least @xmath109 by future generations of experiments . since virtually all de models ( other than @xmath2 )
have either @xmath185 or @xmath186 , the need for probing both the equation of state as well as its first derivative become crucial if any deep insight is to be gained into the nature of de .
it is useful to recall that several de models which agree well with the current data arise because of modifications to the gravitational sector of the theory @xcite . for these _ geometrical de _
models , the eos no longer plays the role of a fundamental physical quantity and it would be very useful if we could supplement it with a diagnostic which could unambiguously probe the properties of all classes of de models .
the @xmath169-probe discussed earlier , provides one such method , since it is based on the expansion history @xmath187 . since the expansion factor @xmath25 is an essential feature of all metric theories of gravity , it is worthwhile trying to explore the properties of de by considering the generic form [ eq : taylor0 ] a(t ) = a(t_0 ) + a_0 ( t - t_0 ) + ( t - t_0)^2 + ( t - t_0)^3 + ...  .
cosmic acceleration appears to be a fairly recent phenomenon @xcite , so we can confine our attention to small values of @xmath188 in ( [ eq : taylor0 ] ) .
the second and third terms in the rhs of ( [ eq : taylor0 ] ) have played a crucial role in the development of cosmology .
indeed , for a long time the large errors in @xmath189 and @xmath190 impaired a precise picture of cosmic expansion and only recently has the sign of @xmath191 been determined to sufficient accuracy for us to make the statement that the universe is accelerating .
links the ricci scalar @xmath17 to @xmath191 .
therefore , an accelerating ( @xmath192 ) spatially flat universe corresponds to @xmath193 . ] keeping in mind the significant progress expected in observational cosmology over the next decade , we feel the time to be ripe to supplement @xmath194 and @xmath191 with @xmath195 , which is the next logical step in the hierarchy of cosmological parameters .
it is remarkable that @xmath196 for the spatially flat @xmath2cdm model which is the ` standard ' cosmological model at present , as well as for the cdm model .
supplementing @xmath197 ( the first statefinder ) with the second statefinder @xmath198 permits us to break this degeneracy and to characterize different de models in a very informative manner .
was also dubbed the cosmic jerk in @xcite . ]
the statefinder pair @xmath199 is a ` geometrical ' diagnostic since it depends upon @xmath25 and hence upon the space - time geometry .
an important property of the statefinder is that spatially flat lcdm corresponds to the fixed point _
lcdm = 1,0  .
[ eq : fixedpt ] so the basic question of whether de is the @xmath2-term or ` something else ' can be rephrased into whether the equality ( [ eq : fixedpt ] ) is satisfied for de . indeed , the departure of a de model from @xmath200 provides a good way of establishing the ` distance ' of this model from lcdm @xcite .
the fact that different classes of de models show distinctly different behaviour when plotted in the @xmath199-plane adds to the practical utility of this diagnostic .
see figure [ fig : statefinder ] and @xcite for some applications of the statefinder .
finally , the pair @xmath199 can also be expressed in terms of @xmath201 quite simply [ eq : statefinder1 ] r & = & 1 + ( 1+w ) - , + s & = & 1+w- .
it is instructive to note that the second statefinder can be rewritten as @xcite s=(+ p)p   [ eq : statefinder2 ] where @xmath202 is the _ total _ pressure including that of de . from ( [ eq : statefinder2 ] ) we find that @xmath203 is very sensitive to the epoch when the total pressure in the universe vanishes . for @xmath2cdm ,
this takes place at the redshift @xmath204 , when the positive radiation pressure cancels the negative pressure in @xmath2 .
consequently , @xmath205 as @xmath206 . in general , the redshift @xmath207 at which the total pressure in the universe becomes zero is quite sensitive to the nature of de .
-18pt    lll level & geometrical parameter & physical parameter +   + 1 & @xmath208 &  @xmath209 , + & &  @xmath210 +   +   + 2 & @xmath211 &  @xmath212 + & @xmath213 &  @xmath214 +   +   + 3 & @xmath215 &  @xmath216 ,  @xmath217 +   + & @xmath218 & +   +    [ table : dark ]    finally , returning to our analogy between an inflaton and dark energy , an interesting correspondence exists between the @xmath219 pair describing de and the _ slow roll parameters _
@xmath220 , @xmath221 , describing inflation q + 1 & = & + r - 1 & = & ( 2 - 3 )  .
we end this section by noting that our current description of cosmology relies on parameters which are either _ geometrical _ or _ physical _ in nature . in the table above , we have divided cosmological parameters according to their _ level _ which is related to the number of differentiations of the expansion factor needed to construct that parameter .
note that all relative energy densities in the table @xmath222 are defined using the present value @xmath223 of the newton gravitational constant in agreement with the definition ( [ einst ] ) .
each higher level requires an additional differentiation of observational data @xmath224 and , therefore , demands a higher level of accuracy for the latter . for current data
, we broadly have about 10% accuracy for quantities belonging to the first level , 50% for the second level , while determination of third level parameters lies in the future . by smoothing ,
we lower the level of a parameter by unity .    relations between geometrical and physical parameters are the following :    w(z ) = 2 q(z ) - 1 3 ( 1 - _ m(z ) )  ; + _ v(z ) = - _
m0(1+z)^3  ; + _ t(z ) = - _
mo(1+z)^3  ; + _ ( z ) = ( r - 3q - 4 + _ m0 ( 1+z)^3 )  .
[ relation ]
the nature of dark energy is clearly one of the outstanding physical and cosmological puzzles of this century . although many distinct theoretical models have been advanced to explain the cosmic acceleration , an alternative approach to study and understand de is to determine de properties from observational data in a model independent fashion , thus , _ reconstructing _ de from observations . in this review ,
we have attempted to briefly summarize some important approaches of cosmological reconstruction .
these approaches are briefly classifiable into parametric and non - parametric methods .
both methods appear to give consistent results when applied to current data .
future directions for model independent reconstruction of dark energy have also been briefly discussed . in the first approximation , within current observational errors
, the reconstructed de behaves like a cosmological constant . but
future , much more precise data , may well show some deviations from this behaviour .
however small , these could have a profound influence on the whole of physics .
_ acknowledgments _ : we are grateful to dragan huterer , irith maor , yun wang and jun - qing xia for agreeing to provide figures for this review .
vs acknowledges a useful conversation with ranjeev misra . as was partially supported by the russian foundation for basic research , grant 05 - 02 - 17450 , and by the research program
 astronomy " of the russian academy of sciences .
he also thanks the centre emile borel , institut henri poincare , paris for hospitality in the period when this paper was finished .
a. g. riess , 1998 , , 1009 [ astro - ph/9805201 ] ; s. j. perlmutter , 1999 , , 565 [ astro - ph/9812133 ] ; r. a. knop , 2003 , , 102(k ) [ astro - ph/0309368 ] .
a. g. riess , 2004 , , 665(r ) [ astro - ph/0402512 ] .
p. astier , 2005 , astron .
astroph . * 447 * , 31 [ astro - ph/0510447 ] .
s  m. carroll , 2001 , living rev.rel . * 4 * , 1 [ astro - ph/0004075 ] ; p.
j  e.  peebles and b.  ratra , 2003 , rev.mod.phys . *
75 * , 559 [ astro - ph/0207347 ] ; t. padmanabhan , 2003 , phys . rep . * 380 * , 235 [ hep - th/0212290 ] ; v. sahni , 2004 , lect .
notes phys .
* 653 * , 141 [ astro - ph/0403324 ] ; 2005 , astro - ph/0502032 ; s.  nojiri and s.  d.  odintsov , 2006 , hep - th/0601213 ; r.  p.  woodard , 2006 , astro - ph/0601672 ; e. j. copeland , m. sami and s. tsujikawa , 2006 , hep - th/0603057 ; s. bludman , 2006 , astro - ph/0605198 .                                                          a. lue , r. scoccimarro and g.  d.  starkman , 2004 , , 124015 [ astro - ph/0401515 ] ; k.  koyama and r.  maartens , 2006 , jcap * 0601 * , 016 [ astro - ph/0511634 ] ; g.  deffayet , g.  gabadadze and a.  inglesias , 2006 , jcap * 0608 * , 012 [ hep - th/0607099 ] .
v. sahni and l.  m. wang , 2000 , 103517 [ astro - ph/9910097 ] ; s. dodelson , m. kaplinghat and e. stewart , 2000 , , 5276 [ astro - ph/0002360 ] ; c.  rubano , p.  scudellaro , e.  piedipalumbo and s.  capozziello , 2003 , , 123501 [ astro - ph/0311535 ] ; g.  barenboim , o.  mena , c.  quigg , 2005 , , 063533 [ astro - ph/0412010 ] ; s.  nojiri and s.  d.  odintsov , 2006 , , 139 [ hep - th/0603062 ] .
j. m. f. maia and j. a. s. lima , 2002 , , 083513 [ astro - ph/0112091 ] ; m. szydlowski and w. czaja , 2004 , , 083507 [ astro - ph/0309191 ] ; 2004 , , 083518 [ gr - qc/0305033 ] ; n.  j.  nunes and j.  e.  lidsey , 2004 , , 123511 [ astro - ph/0310882 ] ; u.  alam , v.  sahni , t.  d.  saini and a.  a.  starobinsky , 2004 , astro - ph/0406672 ; v.  b.  johri and p.  k.  rath , 2005 , astro - ph/0510017 ; m. sahln , a. r. liddle and d. parkinson , 2005 , , 083511 ; s. tsujikawa , 2005 , , 083512 ; d. polarski and a. ranquet , 2005 , , 1 ; y. wang and k. freese , 2006 , , 449 ; f.  simpson and s.  bridle , 2006 , , 083001 ; a. a. sen , 2006 , jcap * 0603 * , 010 [ astro - ph/0512406 ] ; x. zhang , 2006 , astro - ph/0604484 ; astro - ph/0609699 .                                                            v. gorini , a. kamenshchik and u. moschella , 2003 , , 063509 [ astro - ph/0209395 ] ; w. zimdahl and d. pavon , 2004 , gen .
. grav . * 36 * , 1483 [ gr - qc/0311067 ] ; a. k. d. evans , i. k. wehus , o. gron and o. elgaroy , 2005 , astron . astrophys . * 430 * , 399 [ astro - ph/0406407 ] ; u.  debnath , a.  banerjee and s.  chakraborty , 2004 , class . quant . grav . * 21 * , 5609 [ gr - qc/0411015 ] ; x. zhang , 2005 , , 1 [ astro - ph/0503075 ] ; x. zhang , 2005 , int.j.mod.phys . *
d14 * , 1597 , [ astro - ph/0504586 ] ; x. zhang , feng - quan wu , j. zhang , 2006 , jcap * 0601 * , 003 [ astro - ph/0411221 ] ; p. wu and h. yu , 2005 , int
d * 14 * , 1873 [ gr - qc/05090365 ] ; e.  e.  o.  ishida , 2006 , astro - ph/0609614 ; d. rapetti , s. w. allen , m. a. amin , r. d. blandford , astro - ph/0605683 .