A&A 454, 707-714 (2006)
M. Amarzguioui1 - Ø. Elgarøy1,2 - D. F. Mota1 - T. Multamäki2,3
1 - Institute of Theoretical Astrophysics, University of Oslo, Box 1029, 0315 Oslo, Norway
2 - NORDITA, Blegdamsvej 17, 2100 Copenhagen, Denmark
3 - Department of Physics, University of Turku, 20014 Turku, Finland
Received 9 February 2006 / Accepted 30 March 2006
We investigate f(R) theories of gravity within the Palatini approach and show how one can determine the expansion history, H(a), for an arbitrary choice of f(R). As an example, we consider cosmological constraints on such theories arising from the supernova type Ia, large-scale structure formation, and cosmic microwave background observations. We find that the best fit to the data is a non-null leading order correction to the Einstein gravity. However, the current data exhibits no significant trend toward such corrections compared to the concordance CDM model. Our results show that the oft-considered 1/R models are not compatible with the data. The results demonstrate that background expansion alone can act as a good discriminator between modified gravity models when multiple data sets are used.
Key words: gravitation - cosmological parameters - cosmology: observations - cosmology: theory - large-scale structure of Universe
The combination of Einstein's General Relativity (GR) and ordinary matter, as described by the standard model of particle physics (Eidelman et al. 2004), cannot explain the current cosmological observations. Key observations confronting the matter-only universe are the luminosity-redshift relationship from observations of supernovae of type Ia (SNIa, Riess et al. 2004), the matter power spectrum of large-scale structures as inferred from galaxy redshift surveys such as the Sloan Digital Sky Survey (SDSS, Tegmark et al. 2004), and the 2dF Galaxy Redshift Survey (2dFGRS, Colless et al. 2001), and from the anisotropies in the cosmic microwave background radiation (CMBR, Spergel et al. 2003). Two exotic components are required in the matter-energy budget of the Universe to account for the results from all of these cosmological probes within GR. These two components are dark matter, a collisionless and pressureless fluid that contributes about of the universe energy budget, and a negative pressure fluid called dark energy. Currently the dark energy component dominates the energy density of the universe, causing accelerating expansion.
Despite the very good agreement between the so-called concordance model and the astrophysical data, the nature of dark matter and dark energy is one of the greatest mysteries of modern cosmology. In fact, none of the dark matter candidates from high energy physics beyond the standard model (Amendola 2000; Ellis 2000; Brookfield et al. 2005; Bertone et al. 2005) have ever been observed. One should bear in mind that the existence of the dark sector is only inferred from the motion of ordinary matter in a gravitational field, so one can ask whether the necessity of including dark matter and dark energy in the energy budget might not be a sign of our lack of understanding of gravitational physics (Lue 2003; Lue et al. 2004b,a). A natural alternative to adding new exotic fluids is to modify gravitational physics. Several ways of modifying gravity have been proposed to dispense with dark matter (Sanders & McGaugh 2002; Skordis et al. 2006; Moffat 2004; Bekenstein 2005; Sellwood & Kosowsky 2001; Milgrom 1994) or dark energy (Nojiri & Odintsov 2003; Meng & Wang 2003; Arkani-Hamed et al. 2002; Bento et al. 2003; Dvali & Turner 2003; Shao et al. 2005; Gong et al. 2004; Lima 2004; Dvali et al. 2000; Deffayet 2001; Nojiri & Odintsov 2004b,a; Abdalla et al. 2005; Bento et al. 2002; Ahmed et al. 2004; Nojiri & Odintsov 2005; Dolgov & Kawasaki 2003; Freese & Lewis 2002). They can account for the observations with some degree of success (Multamaki et al. 2003; Bento et al. 2003; Elgaroy & Multamaki 2005; Bento et al. 2002; Deffayet et al. 2002; Amarzguioui et al. 2005; Multamaki & Elgaroy 2004; Lue 2003; Lue et al. 2004b,a; Carroll et al. 2005; Multamaki et al. 2004; Koivisto et al. 2005).
In this article we investigate a family of alternative models to the dark energy paradigm based on a generalization of the Einstein-Hilbert Lagrangian. These models are called nonlinear theories of gravity (Magnano & Sokolowski 1994; Allemandi et al. 2006) or f(R) theories (Cognola et al. 2005; Ezawa et al. 2003; Bronnikov & Chernakova 2005; Rippl et al. 1996; Barrow & Ottewill 1983; Schmidt 1998; Nunez & Solganik 2004; Barraco et al. 2000), since the scalar curvature R in the Einstein-Hilbert Lagrangian is replaced by a general function f(R). The main motivation for this generalization is that higher-order terms in curvature invariants (such as R2, , , etc.) have to be added to the effective Lagrangian of the gravitational field when quantum corrections are considered (Birrell & Davies 1982; Buchbinder et al. 1992; Vilkovisky 1992; Gasperini & Veneziano 1992). Furthermore, there is no a priori reason to restrict oneself to the simple Einstein-Hilbert action when a more general formulation is allowed. Higher-order terms in the gravitational action also have interesting consequences in cosmology, like natural inflationary behavior at early times (Barrow & Cotsakis 1991; Starobinsky 1980; La & Steinhardt 1989; Barrow & Cotsakis 1988), and late time acceleration of the Universe (Meng & Wang 2003,2004b; Capozziello et al. 2005; Carroll et al. 2005). Consequently, several authors have investigated whether such theories are indeed compatible with current cosmological observations, big bang nucleosynthesis, and solar system constraints in its weak field limit (Clifton & Barrow 2005; Olmo 2005b; Capozziello 2002; Carloni et al. 2005; Barraco et al. 2002; Mena et al. 2005; Olmo 2005a; Quandt & Schmidt 1991; Sotiriou 2005a; Barrow & Clifton 2006; Sotiriou 2005b; Allemandi et al. 2005; Dominguez & Barraco 2004; Allemandi et al. 2004a,b; Hwang & Noh 2001). Most of these investigations have been model dependent and the conclusions are somewhat contradictory (Flanagan 2004; Vollick 2004).
When dealing with f(R) theories of gravity, the choice of the independent fields to vary in the action is a fundamental issue. The so-called Palatini approach considers the metric and the connection to be independent of each other, and the resulting field equations are in general different from those one gets from varying only the metric in the so-called metric approach. The two approaches lead to the same equations only if f(R) is linear in R. The correct choice of approach to derive the field equations is still a hot topic of research. Initially, f(R)theories were investigated in the metric approach. However, since this method leads to fourth-order equations and the Palatini approach leads to second-order equations, the latter is appealing because of its simplicity. Moreover, the equations resulting from the metric approach seem to have instability problems in many interesting cases (Dolgov & Kawasaki 2003; Chiba 2003) that the Palatini approach does not have. However, recent work has cast doubts on these instabilities (Cembranos 2005). In the present paper we concentrate on the Palatini approach.
The aim of this article is to use current cosmological data to consider possible deviations from GR by combining a number of different cosmological observations. The data used are the latest Supernovae Ia gold set (Riess et al. 2004), the CMBR shift parameter (Bond et al. 1997), the baryon oscillation length scale (Eisenstein et al. 2005), and the linear growth factor at the 2dFGRS effective redshift (Hawkins et al. 2003; Wang & Tegmark 2004). As far as we are aware, this is the first time one uses all of the main cosmological data sets in order to constrain these models.
Newtonian and Solar system constraints within the Palatini formulation have been analyzed by several authors, but there doesn't seem to be any consensus on this topic at present time. To illustrate this we list briefly the conclusion reached by a few authors: Sotiriou (2005a) claims that "any reasonable f(R) model will give the correct Newtonian limit''. Meng & Wang (2004a) claim that all models with inverse powers of R give a correct Newtonian limit; i.e. there are no Newtonian constraints on such models. Dominguez & Barraco (2004), however, claim that the Newtonian limit places very strong constraints on the form of f(R). Using a scalar-tensor representation of the f(R) model, Olmo (2005b) also finds that there are Newtonian constraints on the Lagrangian.
The structure of the paper is as follows. In Sect. 2 we deduce and summarize the basic equations and properties of general f(R)gravities in the Palatini approach. In Sect. 3, we investigate observational constraints based on the evolution of the background of the Universe. In particular, we consider fits to the SNIa and CMB shift parameter from the WMAP data. In Sect. 4, we analyze the evolution of linear perturbations and probe large-scale structure formation in these models using the linear growth factor derived from the 2dFGRS data. Finally, Sect. 5 contains a summary of our work and our conclusions.
The action that defines f(R) gravity theories in the Palatini formalism is
Varying the above action with respect to
the metric, we obtain the generalized Einstein equations
The generalized Ricci tensor can now be written as
Since we are interested in cosmological solutions, we consider the
spatially flat FRW metric,
The generalized Friedmann equation can be derived straightforwardly using the generalized
Ricci tensor (7). This gives
If the equation of state of the fluid, ,
one can use the continuity equation of the fluid together with Eq. (9)
In the case of a constant equation of state,
and hence the generalized Friedmann equation
can be written as
In the special case w=0, applicable when the universe is matter dominated,
these equations reduce to
To investigate to what extent observations allow deviations from
general relativity (where f(R)=R), we consider
the following gravity Lagrangian:
Defining and choosing units so that H0=1, we can solve for R0 from Eq. (17). Consistency then requires that substituting the obtained value of R0 into Eq. (12) must give H0=1. Hence, given and , is fixed. As an example, consider the case , which corresponds to the CDM model. From Eq. (12), we have and from Eq. (17), , so is fixed.
Armed with the modified Friedmann equation, we can now consider the constraints arising from cosmological observations. In this section we consider quantities related to the background expansion of the Universe: the SNIa luminosity distance-redshift relationship, the CMBR shift parameter, and the baryon oscillation length scale.
The CMBR shift parameter (Bond et al. 1997; Odman et al. 2003; Melchiorri et al. 2003) in a spatially
flat universe is given by
Since we do not have an explicit expression for the
Hubble parameter in terms of the redshift, it is useful to rewrite
the shift parameter in terms of the curvature scalar, R:
To incorporate measurements from SNIa, it is useful to rewrite the expression for the
luminosity distance as
|Figure 1: The 68, 95, and confidence contours arising from fitting the CMBR shift parameter. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
|Figure 2: The 68, 95, and confidence contours arising from fitting the SN Ia data. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
With the added information from the CMBR in the form of the shift parameter,
the situation improves as one can see from Fig. 3.
Still quite a large degeneracy persists on the level, but
on the level, the model is quite well-constrained and centered
around the concordance CDM model.
|Figure 3: The combined 68, 95, and confidence contours arising from fitting the SN Ia and the CMBR shift parameter data. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
The baryon oscillations in the galaxy power spectrum are imprints from acoustic oscillations prior to recombination, which are also responsible for the the acoustic peaks seen in the CMBR temperature power spectrum. The physical length scale associated with the oscillations is set by the sound horizon at recombination, which can be estimated from the CMBR data (Spergel et al. 2003). Measuring the apparent size of the oscillations in a galaxy survey allows one to measure the angular diameter distance at the survey redshift. Together with the angular size of the CMB sound horizon, the baryon oscillation size is a powerful probe of the properties and evolution of the universe.
The imprint of the primordial baryon-photon acoustic oscillations
in the matter power spectrum therefore provides us with a "standard ruler''
via the dimensionless quantity A (Linder 2005,2003; Eisenstein & Hu 1998; Hu & Sugiyama 1996; Eisenstein & White 2004):
Recently the acoustic signature associated with the baryonic oscillations
has been identified at low redshifts in the
distributions of luminous red galaxies in the Sloan Digital Sky Survey
(Eisenstein et al. 2005), with a value of
|Figure 4: The 68, 95, and confidence contours arising from fitting the length scale associated with the baryon oscillations. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
Using the observed baryon oscillation-length scale, one can constrain
the cosmological model. The confidence contours for the modified
gravity model we are concerned with in this paper,
are shown in Fig. 4.
Once again there is a large degeneracy in the
similar to the case when we fit SNIa (Fig. 2) and the CMBR
shift parameter (Fig. 1). Such degeneracies are
strongly restricted, however, when one combines all of the data sets
in one single plot (see Fig. 5).
|Figure 5: The 68, 95, and confidence contours arising from fitting the combined data from the SNIa, CMBR shift parameter and the length scale associated to the baryon oscillations. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
Using the baryon oscillation-length scale and combining it with the SNIa and the CMBR shift parameter data hence imposes strong constraints on possible deviations from Einstein's General Relativity. This also demonstrates how combining the current data, one can efficiently study and constrain cosmological models from the background expansion only. The best-fit model to the three data sets is and , which is slightly different from the CDM model; but the CDM model is well within the confidence contour.
So far we have only considered observables related to the background evolution. In order to get further information, it is useful to go beyond these "zeroth order'' tests and consider perturbations within f(R) models. Knowledge of the evolution of perturbations allows one to confront models with large-scale structure observations from galaxy surveys. Several authors have investigated cosmological perturbations in generalized gravity theories in the metric approach (Hwang 1991b,a). However within the Palatini formalism, the first steps were made only very recently (Koivisto & Kurki-Suonio 2005). Here we follow a spherical collapse formalism (Multamaki et al. 2003; Lue et al. 2004a) where, by requiring that a general gravity theory respects the Jebsen-Birkhoff theorem, one derives the modified gravitational force law necessary to describe the evolution of the density perturbations at astrophysical scales.
Although the spherical collapse model has been set and used for a long time (Peebles 1993; Padmanabhan 1993) and in many different contexts (Mota & Barrow 2004b; Lue et al. 2004a; Mota & Barrow 2004a; Mota & van de Bruck 2004; Clifton et al. 2005), it has not been applied previously to nonlinear gravity theories. Lue et al. (2004a) start by assuming a generalization of the Jebsen-Birkhoff theorem: for any test particle outside a spherically symmetric matter source, the metric observed by that test particle is equivalent to that of a point source of the same mass located at the center of the sphere. With this one assumption, one can deduce the Schwarzschild-like metric of the new hypothetical gravity theory.
Armed with the Schwarzschild-like metric one can then investigate the evolution of spherical matter overdensities and compare it to the latest large-scale structure data. The idea is to consider a uniform sphere of dust. Imagine that the evolution inside the sphere is exactly cosmological, while outside the sphere is empty space, whose metric (given the Jebsen-Birkhoff theorem) is Schwarzschild-like (as defined by the metric Eq. (2.2) in Lue et al. 2004a). The mass of the matter source (as determined by the form of the metric at short distances) is unchanged throughout its time-evolution. The surface of the spherical mass therefore charts out the metric throughout space as the sphere expands with time, as long as we demand that the cosmological metric just inside the surface of the sphere smoothly matches the Schwarzschild solution just outside. In order to see how the metric depends on the mass of the central source, we just take a sphere of dust of a different initial size and watch its surface chart out a new metric. The procedure for determining the metric from the cosmological evolution is described in Lue et al. (2004a), so we refer the reader to this article for further information and details.
An open question is the validity of the Jebsen-Birkhoff theorem in f(R) theories. Although one cannot explicitly show that the only possible solution to the field equations, when a spherically symmetric ansatz is inserted into them, is the Schwarzschild metric, there are strong indications that this is indeed the case. In fact, this was shown for the case of f(R)=R+R2, where (Ramaswamy & Yasskin 1979). And it was generalized later for any type of invariant of the form R2, even in the case of a non-null torsion (Neville 1980). The authors also claimed that similar results would most probably be valid even for the case of higher-order curvature invariants, such as R3, R4, etc. However, there is no mathematical proof of this as yet, even though there are several studies and proofs for other complex cases such as multidimensional gravities, Einstein-Yang-Mills systems, and conformally transformed metrics (Brodbeck & Straumann 1993; Bronnikov & Melnikov 1995).
We want to follow the evolution of a lump of matter throughout the history
of our Universe. We start by considering a Universe with a
In order to investigate the evolution of density perturbations in this
scenario, we use the spherical collapse model (Peebles 1993).
Consider a top-hat overdensity
of a spherical distribution of
dust with mass M and radius r defined by
The large-scale structure information we choose to use here
is the linear growth rate
measured by the
2dFGRS (Hawkins et al. 2003; Knop et al. 2003; Verde et al. 2002), where
We compare the theoretical value we get for the linear growth
rate of our model with the value measured by the 2dFGRS at its
The constraints arising from the
linear growth rate are plotted in Fig. 6.
Combining this with all the other constraints leads to the
confidence contours shown in Fig. 7.
|Figure 6: The 68, 95, and confidence contours arising from fitting the linear growth factor using the spherical collapse approach. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
|Figure 7: The 68, 95, and confidence contours arising from fitting to the Combined SNIa, shift parameter, baryon oscillations and the linear growth data sets. The parameter values corresponding to the concordance CDM model ( or ) are marked with a cross. The gray area represents a section of the parameter space that is not allowed.|
|Open with DEXTER|
We have investigated observational constraints on f(R) theories within the Palatini formalism. In order to relate these theories to observations, we show how one can determine the expansion history for a given f(R). In a matter dominated universe in particular, determining H(a)is straightforward as expressed by Eqs. (14) and (15).
We have investigated the possible form of the leading correction to standard GR, parameterizing the Lagrangian for gravity as , and used a combination of data sets to determine the allowed ranges of and . This is by no means an exhaustive study. Other interesting forms of f(R) that are definitely worth studying include and perhaps especially f(R)=R-c1/R+c2 R3 (Sotiriou 2005b), but our main purpose here has been to set up the formalism and demonstrate the effectiveness of combining the current data sets.
Using a combination of data sets that probe the background evolution, we find that the current data efficiently constrains the allowed parameter space of the leading correction to GR in the Palatini approach. The best-fit models to the individual data sets are for the SNIa, for the CMBR shift parameter, and for the baryon oscillations. The best fit to the combination of these data sets is , but the CDM concordance model is well within the contour. Note, however, that the commonly considered 1/R model is strongly disfavored by the data.
In order to bring in additional information from the current galaxy surveys, we also considered the growth of structures in these models of modified gravity. By assuming that the new gravitational physics obeys a limited version of the Jebsen-Birkhoff theorem, we can describe the evolution of overdensities in f(R) gravity theories at sub-horizon scales. We find the best-fit model to the linear growth factor alone to be , but the allowed parameter range is degenerate and does not improve constraints derived from the background evolution. To fully utilize the information available from the galaxy survey in the form of the large-scale matter power spectrum, a more detailed analysis is needed along the lines presented in Koivisto & Kurki-Suonio (2005).
In summary, modified gravities provide us with an interesting alternative to the cosmological concordance model with a dominant dark energy component. Modern cosmological data can efficiently constrain such models. These data indicate that currently there is no compelling evidence for non-standard gravity.
We would like to thank T. Koivisto for useful discussion. M.A., Ø.E., and D.F.M. acknowledge support from the Research Council of Norway through project numbers 159637/V30 and 162830/V00. T.M. is supported by the Academy of Finland through project number 108658.