Stability of cosmic structures in scalar-tensor theories of gravity

In the present work we study a concrete model of Scalar--Tensor theory of gravity characterized by two free parameters, and compare its predictions against observational data and constraints coming from supernovae, solar system tests as well as the stability of cosmic structures. First an exact analytical solution at the background level is obtained. Then using that solution the expression for the turnaround radius is computed. Finally we show graphically how current data and limits put bounds on the parameters of the model at hand.


Introduction
In the end of the 90's the most dramatic discoveries in Particle Physics and Cosmology were from the one hand the neutrino oscillations and on the other hand the current acceleration of the Universe [1,2]. Nowadays many well established observational data from Astrophysics and Cosmology show that we live in a spatially flat Universe that expands in an accelerating rate [3]. Dark energy, the fluid component that dominates the evolution of the Universe and drives the current cosmic acceleration, is one of the biggest challenges of modern cosmology, as its nature and origin still remains a mystery. The ΛCDM model with a constant equation of state w = −1 is the most economical one in excellent agreement with current data. However, given the cosmological constant problems other alternatives with an evolving equation of state have been studied in the literature over the a e-mail: grigorios.panotopoulos@tecnico.ulisboa.pt b e-mail: arrincon@uc.cl years. In general all dark energy models fall into two broad classes, namely from the one hand dynamical dark energy models in which one has to introduce a new dynamical field assuming Einstein's General Relativity (GR) [4][5][6][7], and from the other hand geometrical dark energy models in which one assumes an alternative theory of gravity that modifies GR at cosmological scales. In the latter category we find the well-known examples of f (R) theories of gravity [8][9][10][11], the Dvali-Gabadadze-Porrati brane model [12] or Scalar-Tensor theories of gravity (ST), with the Brans-Dicke [13] model being the archetypical one and recently the scale dependent approach previously applied to certain black holes problems [14][15][16][17][18].
Until a few years ago the observational data used to constrain dark energy models were mainly the temperature anisotropies of the cosmic microwave background, galaxy surveys and supernovae data. However, recently a new potentially local check was proposed in [19] based on two facts, namely a) the motion of a test particle depends on the interplay between the initial momentum of the Big-Bang, the attractive nature of gravity and the repulsive nature of dark energy, and b) for a given mass of a spherical structure there is a maximum radius, called the turnaround point R T , beyond which a test particle cannot stay bound due to the antigravity effect of dark energy. This is very similar to what happens in neutron stars where the mass-to-radius relation depends crucially on the poorly known equation-ofstate [20], and observed pulsars with a mass at two solar masses rule out equations-of-state that predict a lower higher value for the star mass [21]. Then in [22] the authors considered dark energy models with a constant equation-of-state w, and soon after that the idea was further pursued in subsequent works applied to generic dark anergy models [23], Brans-Dicke theory [24] and DGP brane model [25]. Unfortunately, in novel cosmologies characterized by non-standard Friedmann-like equations it is highly non-trivial to see the implications of the prediction of the models regarding the maximum turnaround point, although a formula for R T may exist.
ST theories of gravity are straight-forward generalizations and in fact the simplest extension of GR. Given that the stability of cosmic structures based on the maximum turnaround point has not adequately analyzed yet, in the present work we ask ourselves the question what the stability of cosmic structures together with other observational data and limits can tell us about cosmologies based on ST theories of gravity. The goal of this work is two-fold. First we present an exact analytical cosmological solution of a ST theory of gravity, which is always desirable, and then we compute the turnaround radius of the model. Our work is organized as follows: After this introduction we present the model, the cosmological equations at the background level as well as the exact analytical solution in section two. In the third section we compute the expression for the turnaround radius, and in section 4 we use current data to constrain the parameters of the model. Finally we conclude our work in the fifth section. We use natural units such that c = 8πG =h = 1 and metric signature (−, +, +, +).

The model, the cosmological equations and the exact solution
We start by defining the model where S m is the action of matter fields, g is the determinant of the metric g µν , R is the Ricci scalar, φ is the scalar field and V (φ) is the scalar self-interaction potential. The dimensionless function F (φ) describes the variation of the effective gravitational constant. This is a generalization of quintessence models and it is characterized not only by the scalar potential but also by F (φ). Considering a flat FRW ansatz for the metric with t being the cosmic time and a(t) being the scale factor, and assuming no interaction between the scalar field and the non-relativistic matter with pressure p = 0 and energy density ρ, one obtains the following cosmological equation for the background [26,27] 3F where H =ȧ/a is the Hubble parameter, the dot denotes differentiation with respect to cosmic time, while the , φ denotes differentiation with respect to the scalar field. Clearly, when F (φ) = 1 we recover the standard equations valid in GR. Note that there are four equations in total, but only three of them are independent.
In addition, power-law solutions are very common in cosmological models based on GR in various contexts.
Besides the trivial examples of the radiation and the matter dominated era, one can mention the well-known cases of the power-law inflation [28] as well as the study of cosmological scaling solutions [29]. Thus, we seek power-law solution of the form where the subindex 0 denotes present values and a 0 is defined at the present time as unity, and p > 1 corresponds to accelerating solutions. Plugging everything into the equations one can check that all of them are satisfied provided that and it is easy to verify that combining any two of the last three equation we obtain the third. Therefore, we can choose p, λ to be the free parameters of the model, while φ 0 , V 0 are determined by the previous expressions. This is our first main result in this work. We remark in passing that exact analytical solutions have been obtained in [30], but without matter.
Next the behaviour of the set of fields {φ, F, V } as functions of red-shift z = −1 + 1/a is investigated. We combine the aforementioned fields to obtain the dimensionless functions involved, namelỹ is dimensionless by definition. Note that in our scalartensor model, the potential V (z) does not depend on the free parameter p. We plotφ(z), F (z) andṼ (z) as functions of red-shift for different values of p. The three quantities are shown in Figures 1, 2 and 3 respectively.

The maximum turnaround radius in the ST model
To study cosmic structures we need to study the evolution of the metric scalar perturbations Ψ (η, x i ), Φ(η, x i ) defined by [31][32][33] with dη = dt/a being the conformal time. We also need the conservation equation for the peculiar velocity δu i [31][32][33] where H is the conformal Hubble parameter. Following [25] we consider a shell of backreactionless cold dark matter fluid moving just outside the structure. The physical spatial coordinate corresponding to the cold dark matter perturbation is We can now obtain the velocity as well as the acceleration of this element as follows: First the velocity is computed to be while taking the derivative once more we obtain the acceleration In a non-standard cosmology the Poisson equation for sub-horizon scales becomes [25] where δρ is the perturbation of the matter energy density, and G eff is the effective Newton's constant which is different than G N . In the last step we approximate the whole structure as a point mass located at r = 0, and the source in the Poisson equation reads δρ = Mδ 3 (r) (25) with M being the mass of the structure. Then the Poisson equation becomes and therefore the solution reads The maximum turnaround radius by definition is computed by requiring that the acceleration vanishes at that point. Thus we obtain [25] HR which implies We see that the maximum turnaround point depends on the interplay between the background evolutionḢ and the effective Newton's constant G eff . In Scalar-Tensor theories of gravity the effective Newton's constant is given by [26,27] assuming that F F 2 ,φ (which in the end we check that it indeed holds), while constraints from solar system tests require that [26] Given the exact solution we obtained in the previous section we can now compute bothḢ and G eff , and we finally obtain for R T the expression where we have evaluated all the time-dependent quantities at today since the cosmic structures we consider here are nearby structures and thus they correspond to z 0. This is our second main result in the present article.
We recall at this point that in [22,23] it was found that for the ΛCDM model, where w = −1, the turnaround point is given by Therefore to compare with the ΛCDM model we write the previous formula equivalently as follows Therefore we see that the ST cosmological model studied here agrees with the ΛCDM model when the ratio x = p/λ ∼ 1.

Comparison of the model with data
Finally, in this section we briefly compare the ST model considered here against observational data from a) supernovae data, b) solar system tests, and c) stability of cosmic structures based on the maximum turnaround radius obtained in the previous section.

Supernovae data
The Hubble parameter as a function of the red-shift z = −1 + 1/a is computed to be while the luminosity distance is given by [34] d L (z) = (1 + z) Finally the supernovae distance modulus µ = m − M , where M is the absolute and m the apparent magnitude, is given by [35,36] µ(z) = 25 + 5 log 10 d L (z) Mpc (37) In Fig. 4 we show the distance modulus as a function of the red-shift both for ΛCDM and for the ST model studied here for p = 1.25. Observational data from the Union 2 compilation [37] are shown too.

Solar system tests
As already mentioned, constraints from solar system tests require that [26] Given that F (φ) = λ(φ/φ 0 ) 2 and using eq. (15) we obtain the following expression for ω −1 where Ω m,0 is today's value of the normalized density of matter. Fig. 5 shows ω −1 0 as a function of λ for Ω m,0 = 0.27 and p = 1.25. It is easy to check that the constraint from solar system tests requires that λ < 0.013.

Stability of cosmic structures
We recall that in [22,23] it was shown that in dark energy models with a constant equation-of-state parameter w in GR, the stability of cosmic structures requires that w > −2.3. In Fig. 6 we show the prediction for the maximum turnaround radius a) for the ΛCDM model

Conclusions
In this article we have analyzed a concrete model of Scalar-Tensor theory of gravity and we have obtained an exact power-law analytical solution (with matter included). Given that solution the expansion history as well as the effective Newton's constant can be computed explicitly, and the maximum turnaround radius can be computed in terms of the two free parameters of the model. Finally we have used several current observational data and constraints coming from supernovae, solar system tests and stability of cosmic structures to put bounds on the parameters of the model.