Nonadiabatic perturbations in Ricci dark energy model
Abstract
We show that the nonadiabatic perturbations between Ricci dark energy and matter can grow both on superhorizon and subhorizon scales, and these nonadiabatic perturbations on subhorizon scales can lead to instability in this dark energy model. The rapidly growing nonadiabatic modes on subhorizon scales always occur when the equation of state parameter of dark energy starts to drop towards near the end of matter era, except that the parameter of Ricci dark energy equals to . In the case where , the rapidly growing nonadiabatic modes disappear when the perturbations in dark energy and matter are adiabatic initially. However, an adiabaticity between dark energy and matter perturbations at early time implies a nonadiabaticity between matter and radiation, this can influence the ordinary SachsWolfe (OSW) effect. Since the amount of Ricci dark energy is not small during matter domination, the integrated SachsWolfe (ISW) effect is greatly modified by density perturbations of dark energy, leading to a wrong shape of CMB power spectrum. The instability in Ricci dark energy is difficult to be alleviated if the effects of coupling between baryon and photon on dark energy perturbations are included.
Keywords: Cosmological Perturbation Theory, Dark Energy,. CMBR
1 Introduction
The present acceleration of the universe, suggested by various observations [1, 2, 3],is one of the most important puzzles in cosmology. Possible explanations for this puzzle can be based on the modification of theory of gravity or introducing a new form of energy whose pressure is negative called dark energy [4]  [8]. Many models of dark energy have been proposed in literature from various motivations. The interesting one is motivated by the holographic principle [9].
The holographic principle implies a relation between the short distance (UV) cutoff and the long distance (IR) cutoff in quantum field theory [9]. Thus the amplitude of vacuum energy, which may play the role of dark energy, is controled by the long distance cutoff whose length scale should be comparable with one of the cosmological length scales. Choosing the Hubble radius as an IR cutoff, the present amplitude of vacuum energy is in agreement with observations [9]. Nevertheless, the vacuum energy evolves like matter, so that it cannot play the role of dark energy. The equation of state parameter of vacuum energy can be negative if particle or future event horizons are used as the IR cutoff [10, 11], but for the choice of particle horizon it cannot be negative enough to drive the present acceleration of the universe. So we are left with the future event horizon but this choice encounters a causality problem [12]. It has been shown that with the inverse square root of the Ricci scalar, which is propertional to the curvature radius of spacetime, as the IR cutoff, vacuum energy can be a viable dark energy candidate, the Ricci dark energy [13] where
(1) 
where is the energy density of dark energy, is the Ricci scalar and is a constant parameter. The reduced Planck mass, , is set to be unity.
In order to study the dynamics of cosmological perturbations for Ricci dark energy, we follow the literature [14] to suppose that the above equation is also valid in the perturbed spacetime, which should be possible because the relation (1) can also be obtained from some arguments about quantum gravity effects [13]. In [14], the adiabaticity between dark energy and matter perturbations is assumed, i.e. , where , and is the density perturbations in dark energy, energy density of matter and density perturbations in matter respectively. For this case, the instability can occur in the future, and can be avoided if . In this work, we will check whether the adiabaticity between dark energy and matter is preserved in general, and study the dynamics of nonadiabatic perturbations between dark energy and matter as well as the instability in Ricci dark energy model. We suppose that the dynamics of Ricci dark energy can be presented in terms of the perfect fluid and the conservation of energymomentum.
In section 2, we derive some useful relations for the background dynamics of Ricci dark energy. The dynamics of density perturbations is studied in section 3, and the conclusions are given in section 4.
2 Background dynamics
Using the trace of Einstein equation, we can write the energy density and pressure in terms of Ricci scalar as
(2) 
Here, is the pressure of dark energy and , where and is the energy density of Cold Dark Matter (CDM) and baryon respectively. Substituting (1) into the above equation, we get
(3) 
so that the equation of state parameter of dark energy is given by
(4) 
where . The evolution equation for can be obtained from the energy conservation of matter and dark energy
(5) 
where is the Hubble parameter and is the cosmic scale factor. Using (4), the solution of the above differential equation is
(6) 
Here, and
(7) 
where is the present value of . We will use subscript to denote the present value of the variables through out this paper. Since , the evolution of can be presented in terms of as
(8) 
In order to obtain the present acceleration of the universe, must be negative, so that in the early epoch is nearly constant and given by
(9) 
Thus, it follows from (4) that dark energy evolves like a matter in the early epoch, i.e., .
In figure (1), the evolutions of for various values of are plotted. From the figure we see that when looking back in time reaches rather quickly, i.e. (9) is approximately valid during matter domination. Thus, during matter domination, the density parameter of matter and dark energy is given by
(10) 
Since should be in the range [14], (10) shows that is smaller than unity and is not small during matter domination. We will see in the next section that this feature greatly influences the ISW effect.
3 Evolution of perturbations
In order to study dynamics of density perturbations, we start by perturbing the trace of Einstein equation around FRW background as
(11) 
Since we suppose that (1) is also valid in the perturbed spacetime, the above equation can be written as
(12) 
Unlike the background case, the perturbations in CDM and baryon evolve differently, so that we have to analyze their effects on the perturbations in dark energy separately. We will express all perturbed quantities in Newtonian gauge in which the line element is
(13) 
According to the conservation of energy momentum tensor, the evolutions of density contrast and velocity perturbations of CDM obey [15, 16]
(14)  
(15) 
while the evolution equations for density contrast and velocity perturbation of baryon are given by [17, 18]
(16)  
(17) 
where a prime denotes derivative with respect to conformal time , a subscript denotes photon, is the sound speed of baryon, and is the Compton mean free path of a photon scattering off a baryon. Here, depends on the baryon sound speed and the coupling between baryon and photon. In the early epoch, photon and baryon are tightly coupled, so that [18]
(18)  
(19)  
where is the anisotropic perturbations in photon.
The evolution equations for the perturbations in dark energy can also be obtained from the conservation of energy momentum tensor. In the case where the anisotropic perturbations in dark energy is absent, one gets [15, 16]
(20)  
(21) 
where is the momentum density of dark energy.
Substituting (12) into the above equations yields
(22)  
(23) 
where . These equations can be solved if the evolution of is known. Thus, we use (14)  (17) to derive the evolution equation for as
(24)  
(25) 
Here, .
In our case, it is more interesting to describe the perturbations in dark energy in terms of the following variables
(26) 
It can be seen that the adiabaticity between dark energy and matter perturbations corresponds to . Using (22)  (25) The evolution equations for and can be written as
(27)  
From these equations, one can show that on superhorizon scales the evolution of is governed by
(29)  
and on subhorizon scales obeys
(30)  
Here, we keep the terms on the RHS of (29), because in the case of our interest can be very small so that its evolution on large scales can be affected by these terms.
To study how evolves, we start with the simplest case where the main contributions to the matter perturbations in (12) come from CDM perturbations, i.e., and , and the contribution from becomes negligible. In the radiation and matter eras, , so that the second and fourth terms on the RHS of (LABEL:du1) can be neglected. Hence, it follows from (27) and (LABEL:du1) that if the perturbations in dark energy and matter are initially adiabatic in the early epoch, i.e., , they will remain adiabatic at later time as long as both on superhorizon and subhorizon scales. Nevertheless, if the perturbations in dark energy and matter are not adiabatic in the early epoch, the nonadiabatic perturbations can grow both inside and outside the horizon. During the radiation and matter eras, the RHS of (29) can be neglected, so that the solutions of (29) are
(31) 
It can be seen that has a growing mode, and the constant adiabatic mode corresponds to . Similarly, during the radiation and matter eras, (30) becomes a homogeneous differential equation. To solve this equation, we use the variable , where and during radiation and matter domination respectively. In terms of this variable, (30) becomes
(32) 
where and during radiation and matter eras respectively. Since inside the horizon, the nonadiabatic perturbations grow inside the horizon as
(33) 
Again, the constant adiabatic mode corresponds to , and can grow inside the horizon. Since , is in the range during the radiation and matter eras. Thus, from small nonadiabatic perturbations between dark energy and matter in the early epoch, and also can grow by many orders of magnitude, and can become unphysically large at late time. This implies an unphysically large amplitude of and also together with the metric perturbation due to the back reaction which follows from the perturbed einstein equation [16],
(34)  
where a subscript denotes radiation including photon and neutrino, and a subscript denotes matter including CDM and also baryon in general case.
Roughly speaking, the existence of the unphysically large density perturbations, which implies an instability problem in dark energy, can be avoided if the nonadiabatic perturbations between dark energy and matter disappear or have a very small amplitude. However, near the end of matter domination, these nonadiabatic perturbations can be generated although they disappear in the early epoch. This is because becomes negative, so that (29) and (30) become inhomogeneous. Obviously, the constant adiabatic solution cannot be the particular solutions of these equations. On large scales, the terms on the RHS of (29) cannot influence the evolution of so much because their amplitude is small. Nevertheless, inside the horizon, to create the structure in the universe, must grow, consequently the inhomogeneous part of (30) leads to rapid growth of . As a result, , and the metric perturbation rapidly grow and become unphysically large. However, it is easy to see that the inhomogeneous part of (30) will always vanish if , and hence the generation of nonadiabatic mode after matter domination can be avoided for this value of . These qualitative considerations are in agreement with the numerical integration presented in figures (2) and (3) for the case where and respectively. In these figure, the initial conditions for and are chosen such that , and the perturbation variables are normalized such that the curvature perturbation defined below in (39) equals to 4. The numerical solution of the evolution equations for the perturbations and also the CMB power spectrum are obtained using CMBEASY [19]. From these figures one can see that the instability problem in Ricci dark energy can be avoided if .
It follows from the above consideration that if the contributions from baryon perturbations in (12) can be neglected, the instability problem in Ricci dark energy can be avoided if the perturbations in dark energy and matter are adiabatic initially, and . Let us now compute the initial conditions for the perturbations in dark energy and the other species needed for the numerical integration for density perturbations and CMB power spectrum. Since the energy density of Ricci dark energy depends on Ricci scalar, the amplitude of the density perturbation in this dark energy depends on the metric perturbations. The relation between the density and metric perturbations can be obtained from (1) as [14, 20]
(35) 
where the vanishing of the anisotropic perturbations in radiation is not assumed. Hence, on the superhorizon scales and in the radiation dominated epoch,
(36) 
To obtain the above equation, we suppose that is nearly constant on large scales because we use the adiabatic initial conditions for photon and neutrino, and the adiabatic perturbations between dark energy and matter is the constant mode. Using the adiabatic initial conditions derived in [16], the relation between the initial value of the density contrast of photon and metric perturbations can be written in Newtonian gauge as
(37) 
This implies that in the early epoch if the perturbations in matter and dark energy are adiabatic, i.e., , the perturbations in matter and photon cannot be adiabatic. This conclusion holds even though the anisotropic perturbations in neutrino are neglected, because, in this case, in the early epoch.
To study how the nonadiabaticity between radiation and matter perturbations influences OSW effect, we derive the evolution equation for the curvature perturbations on slices of uniform total energy density by ignoring the anisotropic perturbations in photon and neutrino , and treating photon and neutrino as a single perfect fluid. Furthermore, since in the early epoch dark energy evolves like matter and , we also treat dark energy and matter as a single fluid whose density contrast obeys (24). Here, a subscript denotes this combined fluid. The evolution of the density contrast of the total energy density , where and , is governed by the general evolution equation for the perfect fluid
(38) 
where is the total equation of state parameter, and is the total velocity perturbations. From the definition of the curvature perturbations on slices of uniform total energy density [16]
(39) 
and (38), one can show that on large scales
(40) 
where is the total adiabatic sound speed. It follows from the above equation that if the perturbation in radiation and fluid are adiabatic, will vanish. Consequently, is constant on large scales (cee e.g. [21]). In the radiation dominated epoch, the main contribution to comes from the curvature perturbation on the hypersurfaces of uniform radiation energy density, , where the subscript rd denotes the value in radiation dominated epoch. Similarly, during the matter domination, , where the subscript md denotes the value in the matter dominated epoch. Thus, the constancy of yields [22]. In the case of nonadiabatic perturbation, we have to know how the entropy perturbations in radiation and fluid evolve. According to the conservation of energy, obeys
(41) 
Using the above equation and (24), we get , where is the entropy perturbations between radiation and fluid. Hence, on large scales constant [22], so that until the matter domination. However, the detail analysis shows that the nonadiabatic perturbations between radiation and matter can grow in early epoch [23]. Since this growth of nonadiabatic perturbations between radiation and matter does not influence our rough estimation below, we assume that constant, and therefore
(42) 
which vanishes in the radiation dominated epoch and then becomes negative when , because in our notations is negative and is positive. This implies that increases during the transition from radiation domination to matter domination, so that and hence , which is in agreement with the numerical solution plotted in figure (4). Using the approximations and on large scales during last scattering epoch, one gets
(43) 
where is the temperature perturbations in Newtonian gauge [22] and the above expression is evaluated at last scattering. This equation shows that the nonadiabaticity between radiation and matter perturbations enhances the OSW effect compared with the adiabatic case in which .
The main contributions from the density perturbations in Ricci dark energy to the density perturbations in the universe come from the ISW effect which may be represented in terms of . On large scales and after last scattering, is governed by the perturbed Einstein equation
(44) 
Here, we have supposed that . Since and the density contrast have opposite signs in our cases, The contribution from can reduce because it cancels the contribution from the term . The contribution from increases if during matter era increases, and when the contribution from increases decreases. Usually OSW effect is partially canceled by ISW effect [22], a smaller leads to a smaller cancellation and therefore a large amplitude of CMB power spectrum at low multipoles. For the case of Ricci dark energy, is not small during the matter era, so that the contribution from can significantly reduce , leading to a large amplitude of the CMB spectrum at low multipoles. The reduction of due to the contribution from is shown in figure (4). In figure (5), the CMB power spectrum for the case where the instability problem in Ricci dark energy can be avoided, i.e., initially and , is plotted. It can be seen from the figure that the CMB power spectrum is greatly enhanced at small multipoles because the SW effect is modified due to the density perturbations in Ricci dark energy.
Finally, we include the baryon perturbations in in (12), so that the term on the RHS of (LABEL:du1) cannot be neglected. Using the same consideration as above, one can see that this term can generate a growing nonadiabatic mode inside the horizon, leading to an instability problem in Ricci dark energy model. Up to our trial, the generation of this nonadiabatic mode cannot be avoided by adjusting the initial conditions for the perturbations, special values of or can be canceled by the other terms, hence the instability problem is not easy to avoid in this case. The evolution of , and for this case is plotted in figure (6).
4 Conclusions
We study the evolution of density perturbations in Ricci dark energy model. We extend the analysis in [14] by allowing the nonadiabatic perturbations between matter and dark energy. In the case where the contributions from baryon perturbations are neglected, if the perturbations in matter and dark energy are adiabatic initially, they will remain adiabatic both on superhorizon and subhorizon scales until the equation of state parameter of dark energy becomes negative. When the equation of state parameter of dark energy becomes negative, the nonadiabatic perturbations between matter and dark energy perturbations will be generated even though they disappear in the early epoch. This nonadiabatic mode is a growing mode that can lead to an instability in Ricci dark energy, i.e., the density perturbations in dark energy, matter, and metric perturbations grow extremely large. The generation of this nonadiabatic mode can be avoided if and the perturbations in dark energy and matter are adiabatic initially. Nevertheless, the perturbations in matter and photon cannot be adiabatic in the early epoch, so that the OSW effect is modified. The main influence of the perturbations in dark energy on the density perturbations in the universe is the great modification of the ISW effect due to a large amount of dark energy during matter era. This modification of ISW effect greatly enhances the CMB spectrum at low multipoles. When the contributions from baryon perturbations are included, the coupling between baryon and photon can generate growing nonadiabatic mode that leads to the instability in Ricci dark energy. In this case, the generation of growing nonadiabatic mode cannot be avoided by simply adjusting the initial conditions for the density perturbation or spacial choices of .
Acknowledgments
The authors would like to thank A. Ungkitchanukit for comments on the manuscript. This work is supported by Thailand Research Fund (TRF) through grant RSA5480009.
References

[1]
A. G. Riess et al. [Supernova Search Team Collaboration],
2004 Astrophys. J. 607 665
[arXiv:astroph/0402512];
R. A. Knop et al., [Supernova Cosmology Project Collaboration], 2003 Astrophys. J. 598 102 [arXiv:astroph/0309368].  [2] C. L. Bennett et al., 2003 Astrophys. J. Suppl. 148 1 [arXiv:astroph/0302207].
 [3] M. Tegmark et al. [SDSS Collaboration], 2004 Phys. Rev. D69 103501 [arXiv:astroph/0310723].
 [4] J. Frieman, M. Turner and D. Huterer, 2008 Ann. Rev. Astron. Astrophys. 46 385 [arXiv:0803.0982 astroph].
 [5] C. W. Chen, JeAn Gu and P. Chen, 2009 Mod. Phys. Lett. A 24 1649 [arXiv:0903.2423 astroph.CO].
 [6] C. Wetterich, 1988 Nucl. Phys. B 302 668.
 [7] B. Ratra and J. Peebles, 1988 Phys. Rev. D 37 321.
 [8] R. R. Caldwell, R. Dave and P. J. Steinhardt, 1998 Phys. Rev. Lett. 80 1582 [arXiv:astroph/9708069].
 [9] A. G. Cohen, D. B. Kaplanand and A. E. Nelson, 1999 Phys. Rev. Lett. 82 4971 [arXiv:hepth/9803132].
 [10] W. Fischler and L. Susskind, 1999 JHEP. 9907 004 [arXiv:hepth/9806039].
 [11] M. Li, 2004 Phys. Lett. B 603 1 [arXiv:hepth/0403127].
 [12] R. G. Cai, 2007 Phys. Rev. D 76 123006 [arXiv:0707.4049v3 hepth].
 [13] C. Gao, X. Chen and Y. G. Shen, 2009 Phys. Rev. D 79 043511 [arXiv:0712.1394 astroph].
 [14] C. J. Feng and X. Z. Li, 2009 Phys. Lett. B 680 184187 [arXiv:0904.2972 hepth].
 [15] K. A. Malik and D. Wands, 2009 Phys. Rept. 475 151 [arXiv:0809.4944 astroph].
 [16] M. Doran and C. M. Müller, G. Schäfer and C. Wetterich, 2003 Phys. Rev. D 68 063505 [arXiv:astroph/0304212].
 [17] C. P. Ma and E. Bertschinger, 1995 Astrophys. J. 455 725 [arXiv:astroph/9506072].
 [18] M. Doran, 2005 JCAP 0506 011 [arXiv:astroph/0503277].
 [19] M. Doran, 2005 JCAP 0510 011 [arXiv:astroph/0302138].
 [20] A. Riotto, 2002 Trieste 2002, Astroparticle physics and cosmology* 317 [arXiv:hepph/0210162].
 [21] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, 2000 Phys. Rev. D 62 043527 [arXiv:astroph/0003278].
 [22] W. Hu and N. Sugiyama, 1995 Phys. Rev. D 51 2599 [arXiv:astroph/9411008].
 [23] I. A. Brown, A. J. Christopherson and K. A. Malik, 2011 [arXiv:1108.0639].