# Interacting Dark Energy in cosmology : A Dynamical System analysis

###### Abstract

The present work deals with an interacting dark energy model in the framework of cosmology. A cosmologically viable form of is chosen ( is the torsion scalar in teleparallelism ) in the background of flat homogeneous and isotropic Friedmann-Robertson-Walker ( FRW ) space-time model of the universe. The matter content of the universe is chosen as dust interacting with minimally coupled scalar field. The evolution equations are reduced to an autonomous system of ordinary differential equations by suitable transformation of variables. The nature of critical points are analyzed by evaluating the eigenvalues of the linearized Jacobi matrix and stable attractors are examined from the point of view of cosmology. Finally, both classical and quantum stability of the model have been discussed.

Keywords: cosmology, Dark energy, Dynamical system, phase space.

PACS Numbers: 98.80.-k ,95.36.+x, 04.20.-q, 04.20.Cv.

## I Introduction

The prediction through the SNeIa Hubble diagram [1,2] that the universe is at present going through on accelerating phase, has been confirmed by wide range of data, from more recent SNeIa data to BAOs and CMBR anisotropies [3-7]. However, the overwhelming abundance of observational evidences for cosmic speed up does not match with standard cosmology in the frame work of General Relativity (GR). To resolve this paradox, an unexpected new ingredient in the form of a negative pressured component has been introduced but it poses difficult questions on its nature and introduces further problems hard to be solved. So it is naturally speculated that these observational evidences may be the first signal of a breakdown of our understanding of the laws of gravity on cosmological scale. As a result several modified gravity theories has been developed of which -gravity theory[8,9] gets much attention. In this theory the scalar curvature is replaced by a suitably chosen function in the gravity Lagrangian.

An equivalent formulation is represented by teleparallelism where instead of curvature, torsion is responsible for the gravitational interaction [10-12]. This model was first proposed by Einstein for unifying electromagnetism and gravity on Weitzenböck non-Riemannian manifold. The Weitzenböck connection replaces the Levi-Civita one on the underlying Riemann-Cartan space-time. Thus gravitational interaction is not purely geometrical rather, the torsion acts as a force, allowing the interpretation of gravity as a gauge theory of the translation group[13]. Although, conceptually, GR and teleparallel gravity are distinct but they yield equivalent dynamics at least at the classical level.

In analogy with -gravity, a generalization to teleparallel gravity is obtained by replacing ( the torsion scalar ) with a generic function . A particular important consequence is the breakdown of the equivalence with the classical GR and hence predicts different dynamics [14-16]. Further, modified teleparallel gravity preserves the field equations to be still second order in the field derivatives in contrast to fourth order equations in -gravity. However, this modified theory suffers from the lack of Local Lorentz Invariance (LLI), hence all the sixteen components of the vierbien are independent and one can not fix six of them by a gauge choice[17].

For the past few years, models based on DE interacting with dark matter ( DM ) or any other exotic matter components have gained great impetus. Such interacting DE models can successfully explain numerous cosmological puzzles namely, phantom crossing, cosmic coincidence and cosmic age problem [18-23].

The present paper is devoted to the study of dynamics of interacting dark energy in cosmology. The interaction between DE and DM could be a major issue to be confronted in studying the physics of DE. However, due to the nature of these two components remaining unknown, it will not be possible to derive the precise form of the interaction from first principles. One has to assume a specific coupling from the outset or determine it from phenomenological requirements. Further, in the framework of field theory it is natural to consider the inevitable interaction between the dark components. An appropriate interaction between DE and DM can provide a mechanism to alleviate the coincidence problem. Moreover, complementary observational signatures of the interaction between DE and DM have been obtained from the cosmic expansion history by using WMAP, SNIa, BAO and SDSS data as well as the growth of cosmic structure. Interestingly, it was disclosed that the late integrated Sachs-Wolfe effect has the unique ability to provide insight into the coupling between the dark sectors. Further, in view of the continuity equations, the interaction between DE and DM must be a function of the energy densities multiplied by a quantity with units of the inverse of time which has the natural choice as Hubble parameter.

Usually one assumes a phenomenological form of interaction between matter and dark energy which are the dominant components of the cosmic composition [24-28] today. In fact, the interaction term gives the rate of exchange of energy density in the dark sector.

The motivation for choosing such a complex system is to test of such a weird model could explain the overall cosmological evolution. Due to complicated field equations, we study the dynamics in the phase space associated to this scenario around both hyperbolic and non-hyperbolic critical points. Considering first order perturbations near the critical points, we examine the nature of critical points by considering the eigenvalues of first order perturbed matrix. Previously, similar dynamical system analysis has been done by us in Dirac-Born-infeld gravity theory [29] and in Brane Scenario [30]. In DBI model, critical points are evaluated, analyzed and their stability have been discussed. In brane gravity model which is essentially on higher dimensional theory, in addition to the above analysis the classical stability of the model itself has been examined. In the present gravity model, choosing the function ’’ appropriately the evolution equations are converted to an autonomous system by suitable choice of the auxiliary variables. In addition to the above course of analysis in the present work we shall also examine the stability of the system both classically as well as quantum mechanically.

We follow the plan : In section II, we present the basic theory of gravity. In section III, we introduce the basic equations and formation of dynamical system. While the critical points of the dynamical system and phase space analysis to the corresponding critical points are shown in section IV. In section V and VI, we work out the stability analysis of the critical points. The cosmological implications at the critical points are discussed in section VII. We conclude and summarize in section VIII. Throughout the paper we use natural units ( ).

## Ii - gravity theory

A possible way to modify gravity beyond general relativity is to consider the Weitzenböck connection, which has no curvature but torsion, rather than the curvature defined by the Levi-Civita connection. Such an approach is termed as ’teleparallelism’[12,31-33], which was also considered by Einstein himself[10,11]. However, in the context of the present accelerated expansion of the universe, the teleparallel Lagrangian density described by the torsion scalar has been extended to a function of ( ) [34,35], equivalent to the concept of gravity.

In this new modified gravity theory ( called -gravity ), the gravity is no longer caused by curved space-time but torsion and moreover, the field equations are only second order unlike the fourth order equations in the theory.

In gravity theory the action is written as

(1) |

Where is the torsion scalar, is a differentiable function of the torsion, corresponds to the matter Lagrangian, and . The torsion scalar is defined as

(2) |

with

(3) |

(4) |

(5) |

In teleparallelism, orthogonal tetrad components are used as dynamical objects and are cosidered as an orthonormal basis for the tangent space at each point of the manifold. So we have

. Here each vector can be described by its components , =0,1,2,3 ( =0,1,2,3) in a coordinate basis i.e. . Here capital letters refer to the tangent space while greek indices label coordinates on the manifold. The metric tensor is obtained from the dual vierbein as

(6) |

In the above is the curvature less Weitzenböck connection and it encompasses all the information about the gravitational field. The contorsion tensor gives the difference between Weitzenböck and Levi-Civita connections.

It is straight forward to show that this equation of motion reduces to Einstein gravity when . Indeed, it is the equivalency between the teleparallel theory and Einstein gravity [12]. The theory has been found to address the issue of cosmic acceleration in the early and late evolution of universe [36] but this crucially depends on the choice of . For instance, exponential or power-law form of can not lead to phantom crossing [37]. But subsequently, it has been shown that crossing of the phantom divide line for the effective equation of state of two specific models in gravity is possible[38] and the best fit results suggest that the observations favour a crossing of the phantom barrier. Then reconstruction of models has been reported in [39,40] while a detailed cosmological analysis is performed in [41,42] and thermodynamics of cosmology [43] including the generalized 2nd law of thermodynamics has been investigated.

Further, it should be noted that a constant acts like as a Cosmological Constant while linear in (i.e ) is simply a redefinition of Newton’s Constant . A desirable choice of is such that at high redshift general relativity should hold so as to agree with primordial nucleosynthesis and cosmic microwave background constraints i.e at early times ( ) [44]. On the other hand, in the asymptotic future there will be de Sitter state. A simplest choice of is the polynomial form [44]

(7) |

where is arbitrary constant.

In particular for the model gives the same expansion history as DGP gravity [44,45] and hence gravity can be considered to be related to higher dimension theories. Note that the restriction gives a viable model compared to current observed data set and also puts the rescale factor to Newton’s constant to be small. The effective DE equation of state varies from in the past to in the future.

Throughout the work we consider the flat homogeneous and isotropic FRW universe with the metric

(8) |

where t is cosmic time. For this metric

where is cosmological scale factor. By combining with (5), (3) and (4) one obtains

where is the Hubble parameter and dot stands for derivative with respect to t.

## Iii Basic Equations in cosmology and formation of Dynamical system :

In gravity theory, we consider flat, homogeneous and isotropic FRW space-time as the model of our universe and the matter is chosen as dark matter in the form of dust(having energy density ) interacting with the dark energy which is chosen as a minimally coupled scalar field having self interacting potential . So the energy density and pressure corresponding to this scalar field are

(9) |

Thus the modified Friedmann equations in - gravity are[14,15]

(10) |

(11) |

The energy balance equations for the individual dark components are

(12) |

and

(13) |

Where the interaction term corresponds to energy exchange between dark energy and dark matter. The positivity of indicates a transfer of energy from dark energy to dark matter. This is required to alleviate the coincidence problem and is compatible with the 2nd. law of thermodynamics. For the time being is unspecified, only it is assumed that does not change sign during the cosmic evolution. Due to the unknown nature of the two dark components (DM and DE), the precise form of the interaction can not be determined from the outset or from phenomenological requirements. However, it is speculated that the interaction term may lead to a major issue to be confronted in studying the physics of DE. Also, from the view point of field theory, interaction between the dark components appears naturally. Moreover, an appropriate interaction between the DE and DM can provide a mechanism to alleviate the coincidence problem. Further, in view of the continuity equations ( i.e. equations (12) and (13) ) the interaction between DE and DM should be a function of energy densities multiplied by a factor having dimension inverse of time and Hubble parameter is a natural choice for it. Thus phenomenologically, can be chosen as (i) , (ii) , (iii) or more generally (iv) . In the present work, for simplicity is chosen as , where the coupling parameter is assumed to be small.

Now, using (9) in the continuity equation (13) the evolution of the scalar field is given by

(14) |

(the modified Klein-Gordon equations)

Due to complicated form of the evolution equations, it is not possible to have analytic solution so to have a qualitative idea about the cosmological behavior we shall put the evolution equations into an autonomous dynamical system. For this, we introduce the new variables

(15) |

which are normalized over Hubble scale. is the density parameter for DM. As a result, the evolution equations reduce to the following autonomous system of ordinary differential equations ( after some algebra )

(16) |

where we choose , a constant such that and the independent variable is chosen as , which is called the e-folding parameter. The system of equations in (16) is analyzed by first equating them to zero to obtain the critical points. Next we perturb equations up to first order about the critical points and check their stability.

In deriving the above autonomous system of ODE we choose (i.e in equation (7) ) and the interaction is chosen in the form . As mentioned earlier, the positive coupling is very reassuring in view of coincidence problem. However, models with negative coupling parameter (indicates decay of DM into DE ) allow for the possibility that there is no DE field in the very early universe and that DE ’Condenses’ as result of the slow decay of DM [46]. Also, it has been shown [47] that the coupling parameter is weakly constrained to negative values by Planck measurements. Further, the negative coupling can not be counted to resolve the tension between the Planck and HST measurements of the Hubble parameter [47]. Although, the negative coupling does not help to alleviate the coincidence problem, it appears in the observed data fittings that models with negative coupling shows most significant departure from zero coupling.

Further, using the normalized variables in the first modified Friedmann equation, we obtain the density parameter for the dark matter as

(17) |

Due to the energy condition , so for fix , lies on the circle . In fact, the phase space of the autonomous system (16) forms a paraboloid ( ) bounded by and . Hence the phase space is finite .

The cosmological parameters related to the scalar field namely the equation of state parameter and the density parameter can be expressed by newly defined variables as

(18) |

and

(19) |

since . Also, the deceleration parameter has the explicit form :

(20) |

and the equation of state parameter of the equivalent geometric matter is

(21) |

## Iv Critical points and Phase -Space Analysis :

We analyze the stability of the corresponding dynamical system about the critical points. We shall plot the phase and evolutionary diagrams accordingly. For this reason we must find the critical points of the system of equations (16) forming an autonomous dynamical system and then we linearize the system near the critical points. The system of equations (16) has the following six critical points for the positive coupling parameter i.e. decaying of DE into matter
:

I. Critical points: ,

II. Critical Points :

III. Critical points :

where , .

On the other hand, for negative coupling parameter we have the following two critical points

IV. Critical points : where .

These critical points and the relevant physical parameters at those points are shown in table I and table II.

1 | 0 | 0 | 1 | 1 | 1 | 2 | 0 | |

-1 | 0 | 0 | 1 | 1 | 1 | 2 | 0 | |

0 | 1 | |||||||

0 | 1 | |||||||

0 | 1 | |||||||

0 | 1 |

From table I, we see that the critical points and always exist ( for all ) while and exist only for . Also, the four critical points , , and represent only the DE components ( DM is absent ). The critical points and exist for and correspond a combination of DM and DE with the ratio of two energy densities
.
For all critical points DE behave like a perfect fluid.
Further, in table II we see that both the critical points and are combination of DM and DE with the ratio of density parameter
. The equivalent geometric matter behaves as dust for the critical points and , it behaves as quintessence matter for the critical points and while the geometric matter acts as phantom fluid for the critical points and .

## V Stability Analysis

We shall now discuss the stability of the critical points ( presented in Table I and the table II ) of the autonomous system (16), considering first order perturbations near the critical points. To examine the nature of critical points, one has to study the eigenvalues of the first order perturbation matrix which have been presented in Table III and table IV below :

0 | |||

0 |

0 | |||

0 |

For classical stability () | For quantum stability | |
---|---|---|

( , ) | ||

and | ||

From table III, for , the critical point is non-hyperbolic in nature otherwise hyperbolic. On the other hand, for , the point is non-hyperbolic in nature otherwise hyperbolic.

The critical point is saddle node for ( see figure 2. )
and is saddle node for ( see figure 1. )
while is unstable node for ( fig: 1. )
and is unstable node for ( fig: 2. )
There are no acceleration phase of the universe near the critical points and . Also, these two critical points represent solutions without matter part of the universe. They are only the kinetic energy dominated solutions.

On the other hand, the critical points and are same in all respect and are hyperbolic in nature. These two critical points are stable solutions if . Also, they behave like attractor solutions for ( i.e. ). But when the two critical points and behave like saddle node ( these are shown in the phase portrait ). There exists an accelerating phase of the universe near these critical points and for .

and are also same in all respect. They are combination of both DM and DE. These solutions are dominated by DE components only for . and correspond to solutions in accelerating or decelerating phase of universe according as .

Now, in the table IV, for negative coupling parameter, we get two critical points namely and which are non-hyperbolic in nature. The solutions are saddle like ( because one eigenvalue is always positive in both of the cases ) due to instability in the eigen direction associated with positive eigenvalue and the stability of an eigen direction associated to a negative eigenvalue. These solutions are dominated by both DM and DE components. There is no accelerating phase near these critical points.

Thus, our main results can be summarized as follows ( see tables I, III, II and IV ):

For all and , the kinetic energy dominated solutions ( points and in table I ) are always unstable points in phase space. In particular, and represent source or repealer ( unstable nodes ) in phase space for ) and are always decelerating () (see figures 3, 4 and 7 ).

For , ( i.e. ); the scalar field dominated solutions ( points and in table I ) are the late time attractor ( see figures 3 , 4 and 7 ), while for the solutions are saddle node ( see figures 5 and 6 ). Also, for , there exists an accelerating phase of the universe near the critical points and .

For all and two critical points and ( see table I ) are non-hyperbolic and are combination of DE and DM components. These are the solutions accelerating or decelerating completely depend on the coupling parameter . By the linear stability theory, it is not easy to get the stability criteria of these critical points.

For , the two non-hyperbolic critical points and and ( see in table II ) can be obtained which are saddle points and are both combination of DE and DM. They are always decelerating ( , see table II ).

In the above figures 1-7 the variables and are chosen in the range [-2,2] with grid points at a difference of . The model parameters and are chosen as and . To examine, whether the present model depends continuously on data we have drawn figures 8 and 9. In fig. 8 we have changed the grid difference to 0.52 while in fig. 9 we have changed the model parameters as and . It is found that both the figures 8 and 9 do not differ significantly from figure 7. So we may conclude that the model depends continuously on data.

local stability | classical stability | quantum stability | ||||

1 | 0 | 0 | unstable | stable | stable if | |

-1 | 0 | 0 | unstable | stable | stable if | |

0 | stable (sink) if | |||||

and | stable if | stable if and | ||||

0 | stable (sink) if | |||||

and | stable if | stable if and | ||||

linear stability fails | stable if | stable if and | ||||

linear stability fails | stable if | stable if and | ||||

0 | unstable (saddle) | unstable | unstable | |||

0 | unstable (saddle) | unstable | unstable |

## Vi Equilibrium points and Stability criteria

In the present three dimensional autonomous system the local stability criteria of an equilibrium point is characterized by the eigenvalues of the perturbation matrix ( presented in table III and table IV for different values of coupling parameter ) and then discussion about local stability is presented . We shall now investigate the classical as well as quantum stability of the model.

In cosmological perturbation, sound speed ( ) has a crucial role in characterizing classical stability. In fact, appears as a coefficient of the term ( is the comoving momentum and ’’ is the usual scale factor ) and classical fluctuations may be considered to be stable when is positive. On the other hand, for quantum instabilities at UV scale we decompose the scalar field into a homogeneous part ( ) and a fluctuation as

(22) |

Then by expanding the pressure , up to second order in , the Hamiltonian for the fluctuations takes the form [48,49]

(23) |

where suffix stands for differentiation with respect to the corresponding variable[50].

For positive definiteness of the Hamiltonian we must have

(24) |

where the first two inequalities are related to quantum stability.

Further, should be less than unity. Because otherwise it is possible to send signals along space-like world lines and this would open a Pandora’s box of classical time travel Paradoxes and this also violates the unitarity of quantum theory. Further, if then from the effective field theory the fluctuations propagate faster than the speed of light and therefore there will be violation of Causality.
In the present cosmological scenario we have

(25) |

so for classical stability ( assuming )

(26) |

while for quantum stability the restrictions are

where which is positive, H is the Hubble parameter and the values of , and are given in the system of equations (16).

We have shown both classical and quantum stability criteria of the model in table V. We shall now discuss about the criteria for the model stability at the equilibrium points ( presented in table I and table II ) when x,y and take the corresponding values at the equilibrium points. From the tables I and III we see that the equilibrium points and are not locally stable( saddle node or unstable node ), they behave like source for and from the above model stability analysis they correspond to classical stability only but conditional quantum stable i.e corresponds to quantum stability if while for , the critical point corresponds to quantum stability. As the equilibrium points and are same in all respect so they correspond to classical stability of the model if while for quantum stability the restrictions are and ( see table VI ). From local stability analysis, we see that the critical points and ( table I ) are stable if . From the above model stability analysis we see that the two equilibrium points and are classical stable if where as they correspond to quantum stability if and ( presented in table VI ). Finally, from the tables II and IV, we see that two equilibrium points and are not locally stable (saddle) and from the above analysis we see that they are not classical as well as not quantum stable. The corresponding condition for stability at each equilibrium point are presented in table VI.

## Vii Equilibrium points and Cosmological implications

From the above phase space analysis of the model we shall now discuss about the cosmological behavior of the model at the equilibrium points. From table I, we find that and indicate the universe is completely dominated by kinetic energy of the scalar field and late time acceleration is not possible. They correspond to flat non-accelerating unstable universe without matter part and so they are not of much interest in the present context.

On the other hand, the points and are interesting from the cosmological point of view (for ). There exists an accelerated expansion of the universe near and which correspond to DE model in the quintessence era and the critical points are stable as well as the model is locally stable. It should be mentioned that these critical points correspond to cosmological solutions which describe the recently observed late time acceleration of the universe ( details are shown in figures ). Thus, our model predicts the correct evolutionary scheme for DE density in the regime of this quintessence scalar field. However, for the restriction the points and are stable from local analysis and at these points the potential energy dominates over the kinetic part and the scalar field behaves as exotic fluid. The equilibrium points and are the potential energy dominated non hyperbolic solutions so we cannot investigate the local stability criteria of the system by the linear stability theory. They are the combination of both DM and DE. Accelerating or decelerating phase of the universe near the points and are fully depended on the coupling parameter of the interaction term of DM and DE. The equilibrium points and are same in all respect for the negative coupling parameter. These are the kinetic energy dominated solutions of the scalar field and there exists an non-accelerating universe near and .

## Viii Discussion and Concluding Remarks

In this paper we have performed a dynamical system analysis of a complicated model based on an alternative theory of gravity- the so called gravity. The matter is chosen as dark species - dark matter and dark energy which are represented by dust and a scalar field with an exponential potential respectively. We have investigated the stability and phase space description of interacting dark energy in cosmology by introducing a simple interaction between dark energy and the matter content of the universe. We have written the general dynamical system equations from basic evolution equations by the suitable variables which are normalized over Hubble scale and have found two sets of critical points for positive and negative coupling parameter ( see table I and II ). In the context of gravity the negative coupling parameter shows that the energy is being transferred from matter to dark energy. We get a finite phase space of the system of evolution equations ( autonomous system ) and phase space forms a paraboloid bounded by and . From the phase space analysis and stability analysis of the critical points ( presented in tables I and II ) we find the stable and unstable solutions ( see figures 1-7 ) for the different values of parameters involved. Also the accelerating or decelerating behaviour of the universe near both the hyperbolic and non-hyperbolic critical points have been given in the last section. It is worth noting that the recently observed accelerated expansion of the universe near the equilibrium points and can be seen from the autonomous system (16) in the quitessence scalar field. These support the recent evolutionary scheme of the universe.

Further, we have investigated the classical as well as quantum stability of the model. Note that these two types of stability are not interrelated because the stability of the critical point is related to the perturbations , and ( the corresponding variables of critical point ). On the other hand, the classical stability of the model is connected to the perturbations (and depends on the conditions ) while the quantum stability is related to the perturbations and the conditions take the form of inequalities (24). Thus the critical points can be classified into three categories namely

(i) unstable points at which the model is stable .

(ii) stable points at which model is unstable and

(iii) stable points with stable ( both classical and quantum ) model.

From the table I, the equilibrium points and are interesting in the present context, they purely describe the scalar field dominated solutions in the model and describe the late time acceleration. So, we conclude that although the equilibrium points , correspond to stable model configuration but locally they are not stable equilibrium points. However, it should be noted that the dimensionless variables introduced in equation (15) are not suitable for a complete analysis. Though, the analysis of the static solutions that correspond to equilibrium points at infinity could be analyzed by Poincare projection technique, but for bouncing scenario the expansion rate may pass through zero making the state space non-compact. Therefore, for future work, it is desirable to choose appropriate compact variables to care of expansions, collapses, static solutions and bounces as well.

###### Acknowledgements.

One of the authors (S C) is thankful to UGC-DRS programme , in the Department of the Mathematics, J.U . S C is also thankful to IUCAA , Pune for research facilities at Library.## References

- (1) A.G. Riess et al., AJ, 116, 1009 (1998)
- (2) S. Perlmutter et al., ApJ, 517, 565 (1999)
- (3) D.N. Spergel et al., ApJS, 148, 175 (2003)
- (4) S. Cole et al., MNRAS, 362, 505 (2005)
- (5) M. Kowalski et al., ApJ, 686, 749 (2008)
- (6) D.N. Spergel et al., ApJS, 170, 377 (2007)
- (7) E. Komatsu, K.M. Smith, J. Dunkley, C.L. Bennett, B.Gold, et al., ApJS 192, 18 (2011)
- (8) S. Capozziello ; S. Carloni; A. Troisi Rec. Res. Develop. Astron. Astrophys. 1, 625 (2003)
- (9) S.M. Carroll, V. Duvvuri, M. Trodden, M. S. Turner, Phys. Rev. D 70, 043528 (2004)
- (10) A. Einstein, Sitzungsber. Preuss. Akad. Wiss. Phys. Math. KI. 217, (1928); ibid., 401, (1930)
- (11) A. Einstein,Math. Annal. 102, 685 (1930)
- (12) K. Hayashi, T. Shirafuji, Phys. Rev. D 19, 3524 (1979) [Addendum - ibid., 24, 3312 (1982)]
- (13) H. I. Arcos, J. G. Pereira, Int. J. Mod. Phys. D 13, 2193 (2004)
- (14) R. Ferraro, F. Fiorini, Phys. Rev. D 75, 084031 (2007)
- (15) R. Ferraro, F. Fiorini, Phys. Rev. D 78, 124019 (2008)
- (16) F. Fiorini, R. Ferraro, Int. J. Mod. Phys. A 24, 1686 (2009)
- (17) B. Li, T.P. Sotiriou, J.D. Barrow, Phys. Rev. D 83, 064035 (2011)
- (18) M. Jamil, M. A. Rashid, Eur. Phys. J. C. 60, 141 (2009)
- (19) M. Jamil, M. A. Rashid, Eur. Phys. J. C. 58, 111 (2008)
- (20) M. Jamil, M. A. Rashid, Eur. Phys. J. C. 56, 429 (2008)
- (21) M. Jamil, F. Rahaman, Eur. Phys. J. C. 64, 97 (2009)
- (22) M. Jamil, Int. J. Theor. Phys. 49, 144 (2010)
- (23) M. Jamil, M. Raja, U. Debnath, Astrophys. Space Sci. 337, 799 (2012)
- (24) H. M. Sadjadi, M. Alimohammadi, Phys. Rev. D 74, 103007 (2006)
- (25) T. Clifton, J. D. Barrow, Phys. Rev. D 73, 104022 (2006)
- (26) G. M. Phys. Rev. D 68, 123507 (2003)
- (27) M. R. Setare, Phys. Lett. B 648, 329 (2007)
- (28) M. R. Setare, Int. J. Mod. Phys. D 18, 419 (2009)
- (29) N. Mahata, S. Chakraborty; arXiv:1501.04441, (2015)
- (30) S. Kr. Biswas, S. Chakraborty; Gen.Relativ.Gravit.47 ,22 (2015)
- (31) F. W. Hehl, P. Von Der Heyde, G. D. Kerlick and J. M. Nester, Rev. Mod. Phys. 48, 393 (1976)
- (32) E. E. Flanagan and E. Rosenthal, Phys. Rev. D 75, 124016 (2007)
- (33) J. Garecki, arXiv:1010.2654 [gr-qc].
- (34) G. R. Bengochea and R. Ferraro, Phys. Rev. D 79, 124019 (2009)
- (35) E. V. Linder, Phys. Rev. D 81, 127301 (2010) [Erratum-ibid. D 82, 109902 (2010)].
- (36) D. Liu, P.Wu, H. Yu, arXiv: 1203. 2016[gr-qc]
- (37) P. Wu, H. Yu, Phys. Lett. B 693, 415 (2010)
- (38) P. Wu, H. Yu, Eur. Phys. J. C 71, 1552 (2011)
- (39) M. R. Setare, M. J. S. Houndjo, arXiv: 1203. 1315v1[gr-qc]
- (40) M. H. Daouda, M. E. Rodrigues,M. J. S. Houndjo, Eur. Phys. J. C 72, 1893 (2012)
- (41) V. F. Cardone, N. Radicella and S. Camera Phys. Rev. D 85, 124007 (2012)
- (42) S. Camera, V. F. Cardone and N. Radicella Phys. Rev. D 89, 083520 (2014)
- (43) K. Bamba, M. Jamil, D. Momeni, R. Myrzakulov, arXiv: 1202. 6114v1[physics.gen-ph]
- (44) G.Dvali, G.Gabadadze and M. Porrati Phys. Lett. B 485, 208 (2000)
- (45) C.Deffayet, G.Dvali and G.Gabadadze Phys. Rev. D 65, 044023 (2002)
- (46) A.A. Costa, X. D. Xu, B. Wang, E. G. M. Ferreira and E Abdalla Phys. Rev. D 89, 103531 (2014)
- (47) C.G. Bohmer, G. Caldera-Cabral, R. Lazkoz and R. Maartens; Phys. Rev. D 78, 023505 (2008)
- (48) F. Piazza, S. Tsujikawa ; JCAP 0407, 004 (2004)
- (49) N. Mahata, S. Chakraborty ; Gen. Rel. Grav. 46 , 1721 (2014)
- (50) Usually , the fluctuations of a time - varying scalar obey Lorentz - violating dispersion relations. However , in cosmological context , Lorentz invariance is always violated considering preferred ( CMB -, comoving observers-) frame.