# Gravity with Perturbative Constraints: Dark Energy Without New Degrees of Freedom

###### Abstract

Major observational efforts in the coming decade are designed to probe the equation of state of dark energy. Measuring a deviation of the equation-of-state parameter from -1 would indicate a dark energy that cannot be represented solely by a cosmological constant. While it is commonly assumed that any implied modification to the CDM model amounts to the addition of new dynamical fields, we propose here a framework for investigating whether or not such new fields are required when cosmological observations are combined with a set of minimal assumptions about the nature of gravitational physics. In our approach, we treat the additional degrees of freedom as perturbatively constrained and calculate a number of observable quantities, such as the Hubble expansion rate and the cosmic acceleration, for a homogeneous Universe. We show that current observations place our Universe within the perturbative validity of our framework and allow for the presence of non-dynamical gravitational degrees of freedom at cosmological scales.

## I Introduction

Cosmological observations suggest the existence of dark matter. There is, also, by now ample evidence that the Universe at late times is accelerating. Do either of these facts require the introduction of new, dynamical degrees of freedom in the gravitational and matter sectors?

In the case of dark matter, the analysis of “bullet cluster” objects such as 1E 0657-56 Clowe et al. (2006) and MACSJ0025.4-1222 Bradač et al. (2008) comes as conclusive evidence for the existence of new degrees of freedom independent of the visible baryonic matter. While the exact nature of these new degrees of freedom are still under debate, the simplest choice – a single, non-relativistic, dark matter fluid – is the center of the cosmological standard model.

Even if we choose to account for the galactic rotation curves with a modification to gravity, it is difficult to do it without the addition of extra degrees of freedom. A consequence of the Lovelock-Grigore theorem Lovelock (1970); Grigore (1992) is that most local, classical, and covariant modifications to the gravitational action lead to new fields Skordis (2008). Alternatives to dark matter that have these properties – such as TeVeS Bekenstein (2004) and STVG Moffat (2006) – themselves have new fields as a consequence of the Lovelock-Grigore theorem.

Today, cold dark matter with a cosmological constant (CDM) satisfies all current observational constraints Komatsu et al. (2008). However, a number of theoretical issues Weinberg (1989) have prompted observational campaigns to search for deviations from the predicted CDM behavior Albrecht et al. (2006). These observations will measure the equation-of-state parameter of the dark energy and search for deviations from the value that corresponds to a cosmological constant. Should their results force us to abandon CDM, we will face a problem for dark energy similar to that we faced for dark matter: do such observations require new dynamical degrees of freedom? So far, the two major approaches to dark energy modeling that predict values of , quintessence Armendariz-Picon et al. (2000) and gravity Carroll et al. (2004), incorporate new dynamical degrees of freedom, which can be cast, in both cases, in the form of a single scalar field.

A number of attempts have been made for modifying General Relativity in a way that generates cosmologies with but without introducing new dynamical fields. Modified source gravity Carroll et al. (2006), based on insights from the Palatini formulation of the Einstein equations Flanagan (2004) forces the (Einstein-frame) scalar field to be non-dynamical by erasing the kinetic term. Cuscuton cosmology Afshordi et al. (2007a, b) also modifies the kinetic term; it appears to have new dynamical degrees of freedom but these are frozen out at the perturbative level. Non-local Barvinsky (2003); Soussa and Woodard (2003) or holographic Larsen et al. (2002) theories may also avoid the introduction of new degrees of freedom. Albeit ambitious, these attempts can face serious problems Iglesias et al. (2007).

In the opposite extreme, because of the Lovelock-Grigore theorem, the theoretical pressures on purely phenomenological models that deviate from CDM but that do not introduce new dynamical degrees of freedom are strong. The fact that observations favoring dark energy are made on relativistic, Hubble scales makes it even harder to avoid the consequence of the Lovelock-Grigore theorem, while remaining internally consistent. Requiring general covariance from the start, for example, prevents one from even phrasing many minimal theories the way MOND did on galactic scales Milgrom (1983). In general, while many theories similar or equivalent to Brans-Dicke Brans and Dicke (1961) may, in the limit of small post-Newtonian corrections, appear only to modify the General Relativistic equations of motion, the new scalars such theories introduce appear as truly independent degrees of freedom in the relativistic regime – for example, on scales approaching the horizon size.

Approaches, however, that tie together two particular sets of observations – homogeneous expansion and linear structure formation – in a phenomenological but self-consistent fashion have also been recently explored. For example, the parametrized post-Friedmann framework Hu and Sawicki (2007) makes minimal assumptions about the underlying physics, requiring only causality, metric structure, and energy-momentum conservation. The Cardassian model Freese and Lewis (2002) introduces a power-law term in the Friedmann equation, which leads to accelerated expansion even in a matter dominated universe. The approach of Ref. Skordis (2008) is similar in spirit, and explicitly shows how a subset of modifications will introduce new degrees of freedom. Although these approaches are general and simple to use in connection to observational data, they can not, by their very nature, generate predictions for other astrophysical settings – the solar system, for example, or nearby (non-cosmological) event horizons. As a result, they can not guarantee that dynamical degrees of freedom are not required for consistency in these other situations.

We follow here a different approach intermediate to the two outlined above, which we call Gravity with Perturbative Constraints (see also Ref. DeDeo and Psaltis (2008) for an initial discussion of cosmology with perturbative constraints). In contrast to the fundamental approach, we do not propose a basic theory but rather remain largely neutral on the details of how a particular set of observations are generated. In contrast to the phenomenological approach, our models have a wide range of applicability and one can, for example, ask whether the cosmological predictions that do not add extra degrees of freedom are consistent with observations in the solar system or of compact objects.

Our aim is to cover as large a class as possible of modifications to gravity that do not require new degrees of freedom. Within our framework, it is possible in principle to incorporate all such modifications that are covariant, metric (in the sense of obeying the Einstein equivalence principle), and conserve energy-momentum. The price we will pay for such generality is that we must remain perturbatively close to General Relativity. While we discuss in detail the exact conditions for perturbative validity below, we caution here that this valid- ity does not require weak fields (as defined with reference to either the Minkowski or homogenous FRW metrics as the background).

Our method of generating such models is known in the literature as that of perturbative constraints Simon (1990); it is an alternative way of associating, in a physically unambiguous fashion, a general action of the gravitational field with a field equation that is of second order. Developed to deal with approximate local actions derived from fundamentally non-local theories, it has a wide range of applicability and does not rely on the existence of an underlying non-locality.

The method of perturbative constraints we use here allows us to retain many of the important features of the action of the gravitational field. In particular, it allows us to retain the translational symmetries that, under Noether’s principle, lead to energy-momentum conservation, as well as enforcing diffeomorphism invariance – the invariance of physical law under changes of coordinates. It does so without introducing new degrees of freedom, while at the same time producing behavior that differs from the General Relativistic case, and allowing for a consistent set of predictions across a wide variety of spacetimes.

## Ii Gravity with Perturbative Constraints

In this paper, we specialize to the case of a flat, homogenous FRW metric,

(1) |

It is important to emphasize here, however, that because we will be modifying gravity in ways that affect cosmological expansion, standard relationships between, e.g., the energy density of various components and the rate of expansion that hold in a flat Universe within General Relativity will not necessarily hold here.

We start with an action for the gravitational field that is a general function of the Ricci scalar curvature. The field equation for this general theory is

(2) |

where . For a flat FRW spacetime, the Ricci scalar curvature is

(3) |

The above field equation obeys the Bianchi identity and hence , which results in the usual conservation law

(4) |

Here and are the energy density and pressure, respectively of the constituent of the Universe. In this paper, we will consider Lagrangian densities of the form

(5) |

where n and are free parameters and the last term in the action is considered only at the perturbative level. The case , was first discussed in Ref. DeDeo and Psaltis (2008).

Henceforth we shall adopt geometerized units where . Following standard procedure, we evaluate the component of the field equation as well as its trace and obtain

(6) |

and

(7) |

In this last equation, and denote the total density and pressure of all constituents of the Universe.

When the function is non linear in the Ricci scalar curvature, Eq. 6 is not a first-order ordinary differential equation in time, as is the standard Friedmann equation. Indeed, as discussed earlier, a non-linear Lagrangian introduces additional degrees of freedom to the gravitational field. The key feature of gravity with perturbative constraints is a consistent way of treating these additional degrees of freedom as fixed and not allowing them to evolve independently Woodard (1989); Simon (1990); DeDeo and Psaltis (2008).

In practice, we combine Eq. 6 with the conservation law of Eq. 4, and obtain a single differential equation for the scale factor that we formally denote by

(8) |

Solving this equation with , i.e.,

(9) |

and requiring the scale factor to be unity at the present epoch, we obtain the General Relativistic solution, which we denote by . We then use this solution to evaluate the derivatives that correspond to all the new degrees of freedom in Eq. (8), i.e.,

(10) |

This last equation is again an ordinary differential equation that is of first order in time and we solve it requiring the scale factor to be unity in the present epoch.

Inspection of the field equations or of the Lagrangian action shows that we can expand the final solution for the scale factor in orders of the parameter , where is the unperturbed Ricci scalar curvature of the Universe in the present epoch. For the case with , we can write the solution for the scale factor in the form

(11) | |||||

Note that, although the modified Friedmann equation includes only terms of order , its solution has higher-order terms because of the non-linearity of the differential equation. Consistency requires that terms of order and higher are negligible and can, therefore, be dropped.

When discussing our results, we will occasionally express them in terms of the contributions of matter and of the cosmological constant to the energy density of the Universe, which are given by

(12) |

and

(13) |

respectively, where . Finally, in order to compare our results to the usual parametrization of the dark energy equation of state, we will also write the Friedmann equation in the XCDM form Peebles and Ratra (2003) as

(14) |

and the equation for acceleration as

(15) |

where is the energy density of an equivalent “dark energy” component with an equation of state . Note that this term includes contributions from both the cosmological constant and by the perturbative term that modifies the Einstein equations.

### ii.1 Case : ,

This situation corresponds to a Universe with zero cosmological constant and . Since we are interested in deviations in the matter dominated regime, we will only consider matter with . In this case, the component of the field equation and its trace become

(16) |

and

(17) |

respectively.

When , the equation for the scale factor describes the familiar matter-dominated Einstein cosmology, for which

(18) |

As a result

(19) |

and

(20) |

Now to order Eq. (16) gives

(21) | |||||

Solving for using Eq. (19) and (20) the contribution to order is

(22) |

We then perform the same expansion up to order . The contribution is

(23) |

As a boundary condition, we require that at the present epoch, which we denote by , the scale factor is unity. Thus, as we evaluate to higher order in , we get corrections to the value of from General Relativity, i.e.,

(24) |

As discussed earlier, it is practically impossible to check the convergence of our method without knowing the form of the Lagrangian action to higher order in the parameter . We can provide, however, a necessary condition for convergence in calculating an observable quantity, by requiring that the correction to the value of this observable calculated to order is not larger than the correction calculated to order .

As an example, we calculate the value of the Hubble constant at the present epoch, , as

(25) |

Similarly we evaluate the deceleration parameter at the present epoch as

(26) |

Requiring the correction to order in the expression for the Hubble constant to be smaller than the correction to order results in the bound , whereas the same requirement for the deceleration parameter leads to . Both of these constraints at the present epoch are very similar and require to have very small values, which cannot lead to an accelerated expansion within the limits of perturbative validity.

An additional problem of using a cosmology with and to account for the observations can be seen by writing the Friedmann equation to order (as in DeDeo and Psaltis (2008)), i.e.,

(27) |

Transforming this equation in the XCDM form requires that

(28) |

Clearly the matter-dominated perturbative gravity theory to order behaves as a dark energy fluid with which is inconsistent with the combined WMAP and BAO data Komatsu et al. (2008).

### ii.2 Case : ,

We have found that, while the case with and may generate accelerating solutions, these solutions lie outside the range where the perturbative expansion to the action is valid. Intuitively, this comes from the fact that by the time acceleration has begun, the matter density is low enough that the zeroth-order relation, , implies that the contribution of the additional term, , is too large to be treated as a small perturbation to the standard gravitational action. This is not the case, however, in a Universe with a non-zero cosmological constant, as we will discuss below.

In the case of a Universe with a non-zero cosmological constant, the component of the field equation and the trace are

(29) |

and

(30) |

respectively.

We will again seek a perturbative solution to this equation of the form (11). When , the Friedmann equation is the same as for CDM, i.e.,

(31) |

The time at the present epoch, i.e., the one for which the CDM solution gives a scale factor of unity, is

(32) |

At orders and , the solutions of the equation for the scale factor, and , respectively, as well as the time that corresponds to the present epoch can be found in closed form but are too long to be displayed here (we present the solution to order in Appendix A). Instead, we write the Friedmann equation to order

(33) | |||||

and the acceleration equation to the same order

(34) | |||||

These two equations allow us to calculate to order two observable quantities at the present epoch, as a function of the parameters and : the Hubble constant, or equivalently the parameter , and the deceleration parameter . We show the result of this calculation in Fig. 1 and 2 respectively, together with the regions where the calculation of each observable is not perturbatively valid.

In both figures, the horizontal axis, i.e., , corresponds to a General Relativistic Universe with a cosmological constant. The WMAP results, interpreted within General Relativity, correspond to a Universe with Komatsu et al. (2008). Keeping this value constant and moving towards positive values of , the predicted value of increases, because the perturbative term in the modified Friedman equation contributes negatively to the rate of expansion of the universe, and the deceleration parameter decreases for the same reason. The opposite is true when the parameter becomes increasingly negative.

The hatched areas in both figures show the regions of the parameter space where the corrections to each observable calculated to order are not negligible and hence the solution is no longer of the same order as the field equation. This is the criterion we have discussed earlier in determining the perturbative validity of our calculations. The particular shapes of the hatched regions are determined by the dependence of the perturbative terms in Eqs. 33 and 34, both of which can be negative, positive, or even zero for different values of the ratio . In both figures, however, the excluded regions lie far from the parameters of the General Relativistic Universe that are consistent with the WMAP data, suggesting that small potential deviations from the General Relativistic predictions can be modeled successfuly within our framework.

Finally, we investigate the parameter at the present epoch for an equivalent XCDM cosmology that would result in the same cosmological expansion rate and acceleration as that predicted by the last two equations. We calculate the equivalent equation-of-state parameter by comparing the last two equations to those of the XCDM framework. The result is, to order ,

(35) | |||||

In Fig. 3, we show contours of constant values of the parameter . The current observed limits on the equation-of-state parameter Komatsu et al. (2008) suggest that deviations from the CDM model that do not introduce new degrees of freedom fall within the perturbative limits of our approach and cannot be ruled out a priori.

### ii.3 Case : ,

Finally, we examine a more general class of theories with that are nevertheless matter dominated. We focus here on the case in order to isolate the expected dependence of our previous results on the value of the parameter .

In this case, the field equation becomes

(36) |

and its solution can be expanded again in terms of the perturbative parameter . The constrain on this parameter for perturbative validity becomes increasingly less restrictive as the exponent decreases.

As with Case 1, we can write the Friedmann equation to order and obtain

(37) |

At the same order the acceleration equation is given by

(38) |

If we write the Friedmann equations in the flat XCDM form, Eq. (14), in the present epoch, we obtain for the equivalent equation-of-state parameter,

(39) |

As a result, increasing the value of the exponent makes the equation-of-state parameter more negative.

## Iii Conclusions

We have presented a formalism – Gravity with Perturbative Constraints – that describes deviations from General Relativistic predictions at cosmological scales without introducing new degrees of freedom. Our formalism preserves diffeomorphism invariance, the Einstein equivalence principle, and energy-momentum conservation. Because it is based on modifications of the Lagrangian of the gravitational field, it allows us to link a wide variety of observations – from compact objects and solar system phenomena, to non-linear clustering and linear structure formation – in order to test for the consistency of assuming no new gravitational degrees of freedom.

Our method is quite different from the Palatini formulation of 1/R gravity discussed by, for example, Ref. Vollick (2003), and expanded on as “modified source gravity” in Ref. Carroll et al. (2006). Since both our methods and those related to Palatini formulations, however, eliminate the extra degree of freedom associated with modifying General Relativity it is worth comparing some of their features. In particular, because we are able to treat our system perturbatively – and, explicitly, because our predictions for densities do not demand a particular behavior when – we do not run into the issues found by Ref. Flanagan (2004), where the microscopic discreteness expected for ordinary baryons and dark matter particles no longers “averages out” correctly to return the homogenous Friedmann equations.

Put explicitly, Ref. Flanagan (2004) noted first that the microscopic discreteness of matter did not strongly affect the averaged metric in General Relativity: i.e., when taking the matter field to be composed of discrete particles (small but still larger than their Schwartzschild radius) the gravitational fields relative to the background are still small and can be consistently be replaced by plus a small, linear perturbation. This behavior was then contrasted with the Palatini formulation, where the field, a function of the matter density that appears as a modification to the Einstein equation, has strong variations from point to point since the behavior of is strongly non-linear in the range regardless of the length scale of the perturbation.

In the case of gravity with perturbative constraints, however, describing the behavior of the equations of motion when may involve increasingly higher order terms in the approximate equations of motion, and – likely for many, though perhaps not all, models where the first terms in the equations of motion are generated by cases discussed in this paper – eventually goes outside the radius of convergence of the series. That our theory can not in all cases make predictions for atomic-scale physics is a feature of the perturbative approach; in contrast to Palatini , the behavior at decouples from these more difficult questions and no particular atomic-scale behavior is required by us for consistency with larger scales.

In this first approach to the problem, we chose to look for modifications of the Einstein-Hilbert action that involve the addition of terms proportional to an inverse power of the Ricci scalar curvature. We were lead to this choice by our goal to describe potential deviations at cosmological scales, without affecting significantly the behavior of gravity in the solar system. One might object that our particular choice of parametrization – a two dimensional space, (, ) – too-drastically narrows down the possible models considered. In a sense, this is a problem with any attempt to cover a space of functions with a finite number of parameters.

Our approach, however, follows the general direction of testing for modifications of gravity in the solar system and in other astrophysical systems by using the Parametrized Post-Newtonian (PPN) framework Will (1971). In a similar fashion to our work, the PPN framework provided a way to link together different observations in various astrophysical settings in a consistent and physically rigorous way. While the PPN model itself required a certain degree of arbitrariness – the functional forms of the potentials – it turned out that the original choices of these functional forms Will and Nordtvedt (1972) were sufficient to cover nearly every theory proposed in the following thirty-seven years. Only in a few cases have failures of PPN been due to an overly restrictive choice of parameters (see, e.g. Alexander and Yunes (2007)). The usefulness of the framework we propose here will depend, of course, on its ability to capture the behavior of plausible gravitational theories that go beyond General Relativity.

The particular set of terms that we included in the Lagrangian of the gravitational field allow us also to study in a more systematic way other phenomenological cosmologies that produce accelerating expansion. For example, the set of parameters discussed in Sec. II.3, i.e., and , provide a homogenous cosmology identical to the Cardassian model of Ref. Freese and Lewis (2002). Within our formalism, however, we can also make predictions for the formation of structure in this cosmology, by taking explicitly into account the modifications to the Poisson equation.

Our formalism is also capable of describing the cosmological expansion histories that are generated by the phenomenological XCDM model Turner and White (1997). This is demonstrated in Fig. 3, for a flat universe, where the values of the equivalent equation-of-state parameter are shown as a function of the perturbative parameter , the matter density, and the value of the cosmological constant. In a similar manner, our framework can also describe the cosmological expansion generated by the model suggested by the Dark Energy Task Force Albrecht et al. (2006), in which the equation-of-state parameter is allowed to change in time as , since the values of the equivalent parameter calculated here depend indeed on cosmic time (cf. Eq. 35). In contrast to the phenomenological model, however, our framework offers additional insight as to how the resulting gravity modifications agree with structure formation observations, or with “fifth-force” constraints from other arenas.

We will study the behavior of the cosmological expansion predicted by gravity with perturbative constraints for a Universe with a finite spatial curvature as well as compare these predictions directly to observations in forthcoming papers.

###### Acknowledgements.

A. C. thanks the members of the Kavli Institute for Cosmological Physics for their hospitality. S. D. thanks the Perimeter Institute (Canada) for hospitality, and Wayne Hu, Justin Khoury and members of the Perimeter Institute cosmology group for helpful discussions. We thank Éanna Flanagan for helpful comments. We also thank Feryal Özel for carefully reading the manuscript. This work was supported in part by an NSF CAREER award at the University of Arizona.## Iv Appendix: Perturbative solution to 4th order for

As discussed in the main text, perturbations to CDM at order has closed-form solution. In particular:

## References

- Clowe et al. (2006) D. Clowe, M. Bradač, A. H. Gonzalez, M. Markevitch, S. W. Randall, C. Jones, and D. Zaritsky, ApJ 648, L109 (2006).
- Bradač et al. (2008) M. Bradač, S. W. Allen, T. Treu, H. Ebeling, R. Massey, R. G. Morris, A. von der Linden, and D. Applegate, astro-ph (2008), eprint arXiv:0806.2320.
- Lovelock (1970) D. Lovelock, Arch. Rat. Mech. Anal. 36, 293 (1970).
- Grigore (1992) D. R. Grigore, Class. Quant. Grav. 9, 1555 (1992).
- Skordis (2008) C. Skordis, gr-qc (2008), eprint arXiv:0806.1238v1.
- Bekenstein (2004) J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004).
- Moffat (2006) J. W. Moffat, JCAP 03, 004 (2006).
- Komatsu et al. (2008) E. Komatsu, J. Dunkley, M. R. Nolta, C. L. Bennett, B. Gold, G. Hinshaw, N. Jarosik, D. Larson, M. Limon, L. Page, et al., astro-ph (2008), eprint arXiv:0803.0547.
- Weinberg (1989) S. Weinberg, Rev. Mod. Phys. 61, 1 (1989).
- Albrecht et al. (2006) A. Albrecht, G. Bernstein, R. Cahn, W. L. Freedman, J. Hewitt, W. Hu, J. Huth, M. Kamionkowski, E. W. Kolb, L. Knox, et al. (2006), eprint astro-ph/0609591v1.
- Armendariz-Picon et al. (2000) C. Armendariz-Picon, V. Mukhanov, and P. J. Steinhardt, Phys. Rev. Lett. 85, 4438 (2000).
- Carroll et al. (2004) S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, Phys. Rev. D 70, 043528 (2004).
- Carroll et al. (2006) S. M. Carroll, I. Sawicki, A. Silvestri, and M. Trodden, New J. Phys. 8, 323 (2006).
- Flanagan (2004) É. É. Flanagan, Phys. Rev. Lett. 92, 071101 (2004).
- Afshordi et al. (2007a) N. Afshordi, D. J. H. Chung, and G. Geshnizjani, Phys. Rev. D 75, 083513 (2007a).
- Afshordi et al. (2007b) N. Afshordi, D. J. H. Chung, M. Doran, and G. Geshnizjani, Phys. Rev. D 75, 123509 (2007b).
- Barvinsky (2003) A. O. Barvinsky, Phys. Lett. B 572, 109 (2003).
- Soussa and Woodard (2003) M. E. Soussa and R. P. Woodard, Class. Quant. Grav. 20, 2737 (2003).
- Larsen et al. (2002) F. Larsen, J. P. van der Schaar, and R. G. Leigh, Journal of High Energy Physics 04, 047 (2002).
- Iglesias et al. (2007) A. Iglesias, N. Kaloper, A. Padilla, and M. Park, Phys. Rev. D 76, 104001 (2007).
- Milgrom (1983) M. Milgrom, ApJ 270, 365 (1983).
- Hu and Sawicki (2007) W. Hu and I. Sawicki, Phys. Rev. D 76, 104043 (2007).
- Freese and Lewis (2002) K. Freese and M. Lewis, Physics Letters B 540, 1 (2002), eprint arXiv:astro-ph/0201229.
- DeDeo and Psaltis (2008) S. DeDeo and D. Psaltis, Phys. Rev. D 78, 064013 (2008), eprint arXiv:0712.3939.
- Simon (1990) J. Z. Simon, Phys. Rev. D 41, 3720 (1990).
- Woodard (1989) D. A. E. . R. P. Woodard, Nucl. Phys. B 325, 389 (1989).
- Peebles and Ratra (2003) P. J. Peebles and B. Ratra, Reviews of Modern Physics 75, 559 (2003), eprint arXiv:astro-ph/0207347.
- Will (1971) C. M. Will, ApJ. 163, 611 (1971).
- Will and Nordtvedt (1972) C. M. Will and K. Nordtvedt, ApJ 177, 757 (1972).
- Alexander and Yunes (2007) S. Alexander and N. Yunes, Phys. Rev. Lett. 99, 241101 (2007).
- Turner and White (1997) M. S. Turner and M. White, Phys. Rev. D 56, 4439 (1997), eprint arXiv:astro-ph/9701138.
- Vollick (2003) D. S.Vollick, Phys. Rev. D 68, 063510 (2003), eprint arXiv:astro-ph/0306630.
- Brans and Dicke (1961) C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961).