Matter Instability in Theories with Modified Induced Gravity

Matter stability is a necessary condition to have a cosmologically viable model. Modified gravity in the spirit of f R theories suffers from matter instability in some subdomains of the model parameter space. It has been shown recently that the late-time cosmic speedup can be explained through an f R -modified induced gravity program. In this paper, we study the issue of matter instability in a braneworld setup with modified induced gravity.


Introduction
There are many lines of astronomical evidence supporting the idea that our universe is currently undergoing a speedup expansion 1-4 .Several approaches are proposed in order to explain the origin of this novel phenomenon.These approaches can be classified in two main categories: models based on the notion of dark energy which modify the matter sector of the gravitational field equations and those models that modify the geometric part of the field equations are generally dubbed as dark geometry in the literature 5-13 .From a relatively different viewpoint but in the spirit of dark geometry proposal , the braneworld model proposed by Dvali, Gabadadze, and Porrati DGP 14 explains the late-time cosmic speedup phase in its self-accelerating branch without recourse to dark energy 15, 16 .However, existence of ghost instabilities in this branch of the solutions makes its unfavorable in some senses 17, 18 .Fortunately, it has been revealed recently that the normal, ghost-free DGP branch has the potential to explain late-time cosmic speedup if we incorporate possible modification of the induced gravity on the brane in the spirit of f R theories 19-24 .This extension can be considered as a manifestation of the scalar-tensor gravity on the brane.Some features of this extension are studied recently 25-29 .Modified gravity in the spirit of f R theories have the capability to provide a unified gravitational alternative to dark energy and inflation 30-37 .A number of viable modified f R gravities are proposed in recent years see 38, 39 and references therein .The cosmological viability of these theories is a necessary condition, and in this respect, there are important criteria for viability such as the fulfillment of the solar system tests.Among these requirements, one of the most important ones is related to the so-called matter instability 40-57 in f R gravity.Matter instability is related to the fact that spherical body solution in general relativity may not be the solution in modified theory in general.This instability may appear when the energy density or the curvature is large compared with the average one in the universe, as is the case inside of a star 58 .In a simple term, matter instability means that the curvature inside a matter sphere becomes very large, leading to a very strong gravitational field.It was indicated that such matter instability may be dangerous in the relativistic star formation processes 59-61 due to the appearance of the corresponding singularity.In this respect, for a model to be cosmologically viable, it is necessary to have matter stability in the model.For a detailed study of the issue of matter instability in f R theories, see 40-58 .Since the f R -modified induced gravity brane f R gravity has the capability to bring the normal DGP solutions to be self-accelerating, it is desirable to see whether this model is cosmologically viable from matter instability viewpoint.So, this letter is devoted to the issue of matter instability in a brane f R gravity.

The Basic Equations
Modified gravity in the form of f R theories is derived by generalization of the Einstein-Hilbert action so that R the Ricci scalar is replaced by a generic function f R in the action where L m is the matter Lagrangian and κ 2 8πG.Varying this action with respect to the metric gives where T m μν diag ρ, −p, −p, −p is the stress-energy tensor for standard matter, which is assumed to be a perfect fluid and by definition f ≡ df/dR.Also, T f μν is the stress-energy tensor attributed to the curvature defined as follows: By substituting a flat FRW metric into the field equations, one achieves the analogue of the Friedmann equations as follows 30-37 : where a dot marks the differentiation with respect to the cosmic time.In the next step, following 28, 29 we suppose that the induced gravity on the DGP brane is modified in the spirit of f R gravity.The action of this DGP-inspired f R gravity is given by where g AB is a five-dimensional bulk metric with Ricci scalar R, while q ab is an induced metric on the brane with induced Ricci scalar R. The Friedmann equation in the normal branch of this scenario is written as 28, 29 3f where r c G 5 /G 4 κ 2 5 /2κ 2 is the DGP crossover scale which has the dimension of length and marks the IR infrared behavior of the DGP model.The Raychaudhuri equation is written as follows: To achieve this equation, we have used the continuity equation for ρ f as where the energy density and pressure of the curvature fluid are defined as follows:

2.9
After presentation of the necessary field equations, in what follows we study the issue of matter instability and cosmological viability in this setup.

The Issue of Matter Instability
In 4 dimensions, the conditions f R > 0 and f R > 0 are necessary conditions for f R theories to be free from ghosts and other instabilities 40-44, 62 .In our brane f R scenario, in addition to f R , there is another piece of information in the action 2.5 the DGP character of the model which should be taken into account when discussing the issue of instabilities.To study possible instabilities in this setup, we proceed as follows: variation of the action 2.5 with respect to the metric yields the induced modified Einstein equations on the brane where E αβ which we neglect it in our forthcoming arguments is the projection of the bulk Weyl tensor on the brane where n M is a unit vector normal to the brane, and S αβ as the quadratic energy-momentum correction into the Einstein field equations is defined as follows: τ αβ as the effective energy-momentum tensor localized on the brane is defined as 25-27 Following 40-44 , the trace of 3.1 , which can be interpreted as the equation of motion for f R , is obtained as

3.5
We parameterize the deviation from Einstein gravity as where is a small parameter with the dimension of an inverse-squared length, and ϕ is arranged to be dimensionless.Typically, ≈ H 2 0 10 −66 eV 2 see, for instance, 42 .By evaluating f as 3.5 can be rewritten as follows:

3.8
This equation to first order of gives 3.9 Now we consider a small region of spacetime in the weak-field regime in which curvature and the metric can locally be approximated by g ab η ab h ab , 3.10 respectively, where R 1 is curvature perturbation, and η ab is the Minkowski metric.In this case, the metric can be approximately taken as a flat one, so ∂ 2 t − ∇ 2 and ∇ α R ∇ α R Ṙ2 − ∇R 2 .Now, 3.9 up to first order of R 1 gives we set for simplicity κ 2 1

6 Advances in High Energy Physics
This relation can be recast in the following suitable form:

3.12
The coefficient of R 1 in the fourth term on the left-hand side is the square of an effective mass defined as

3.13
This quantity is dominated by the term 1/ ϕ due to the extremely small value of needed for these theories.It is therefore obvious that the theory will be stable i.e., m 2 eff > 0 if ϕ f < 0, while an instability arises if this effective mass squared is negative, that is, if ϕ f > 0. Based on this fact and as an example, the f R R γR −n function with γ as a positive quantity in the spirit of normal DGP braneworld scenario suffers from a matter instability for n > 0 and n < −1.For this kind of f R function, γ plays the same role as in 3.6 and is supposed to be positive note that is a small parameter with dimension of an inverse length squared , so ϕ R R −n .The condition for matter stability ϕ f < 0 leads to −1 < n < 0.
As another important point, we focus on the stability of the de Sitter accelerating solution under small homogeneous perturbations in the normal branch see 28, 62, 63 for a similar argument .It is useful to rewrite the Friedmann equation corresponding to a de Sitter brane with Hubble rate H 0 in a form that exhibits the effect of an extra dimension on a 4D f R model as follows 28, 62, 63 :

3.14
Subscript 0 stands for quantities evaluated in the de Sitter space time.We also note that the de Sitter brane is described by the scalar factor a t a 0 exp ξt which leads to R 0 12ξ 2 .H 2 4 is defined as Advances in High Energy Physics 7 Therefore, the presence of the extra dimension implies a shift on the Hubble rate.One can perturb the Hubble parameter up to the first order as δH H t −H 0 .Also by a perturbed Friedmann equation based on 2.6 , one can achieve an evolution equation for δH as 64 The stability condition for the de Sitter solution in the DGP normal branch is positivity of the effective mass squared, M eff 2 > 0. Now M eff 2 can be written as the sum of three In the 4D version of the f R gravity, this summation reduces to 64 In the braneworld version, we except the crossover distance to affect the effective mass.In this respect, m 2 back is a purely background effect due to the shift on the Hubble parameter with respect to the standard 4D case, while m 2 pert is a purely perturbative extradimensional effect 28, 63 .In our setup, these quantities are defined as follows:

3.18
The de Sitter brane is close to the standard 4D regime as long as H 2 0 ∼ H 2 4 , which leads to The assumption that we are slightly perturbing the Hilbert-Einstein action of the brane, that is, f 0 ∼ R 0 > 0, and also the positivity of the effective gravitational constant on the brane at late time, that is, f 0 > 0, implies that f 0 f 0 1 for the last inequality.Now the stability of the de Sitter accelerating solution under small homogeneous perturbations in the normal DGP branch of the model since m 2 back > 0 and m 2 pert < 0 is guaranteed if m 2 eff > m 2   4   which leads to the following condition for more details, see 28, 63, 64 : Based on this condition, f R R γR −n function in the spirit of the normal DGP braneworld exhibits a de Sitter stability if the following condition is satisfied: < nγ 6 4n .

3.21
Figure 1 shows the behavior of M ≡ m 2 eff − m 2 4 versus n and γ.The de Sitter phase is stable in this setup if n −4.12.This is shown with more resolution in Figure 2. Note that the matter stability in this induced gravity braneworld scenario occurs in those values of n that the corresponding de Sitter phase is not stable.

Summary and Conclusion
Matter stability is a necessary condition for cosmological viability of a gravitational theory.Recently, it has been shown that the normal, non-self-accelerating branch of the DGP cosmological solutions self-accelerates if the induced gravity on the brane is modified in the spirit of f R gravity.In this letter, we have studied the issue of matter stability in an induced gravity, brane-f R scenario.We obtained the condition for matter stability in this setup via a perturbative scheme, and we applied our condition for an specific model of the type f R R γR −n .For this type of modified induced gravity, the matter is stabilized on the brane for −1 < n < 0. We have also studied the stability of the de Sitter phase for this type of modified induced gravity.For this type of the modified induced gravity, the de Sitter phase is stable for n − 4.12.Albeit for these values of n, matter is not stable.So, these types of modified induced gravity are not suitable candidates for late-time cosmological evolution.We note however that other types of modified induced gravity such as f R R n exp η/R with  η a constant have simultaneous matter stability and stable de Sitter phase in some subspaces of the model parameter space see 65 .

Figure 2 :
Figure 2: Stability of the de Sitter phase.The de Sitter phase is stable for n −4.12.