Cosmology with interaction in the dark sector
Unless some unknown symmetry in Nature prevents or suppresses a non-minimal coupling in the dark sector, the dark energy field may interact with the pressureless component of dark matter. In this paper, we investigate some cosmological consequences of a general model of interacting dark matter-dark energy characterized by a dimensionless parameter . We derive a coupled scalar field version for this general class of scenarios and carry out a joint statistical analysis involving SNe Ia data (Legacy and Constitution sets), measurements of baryon acoustic oscillation peak at (2dFGRS) and (SDSS), and measurements of the Hubble evolution . For the specific case of vacuum decay (), we find that, although physically forbidden, a transfer of energy from dark matter to dark energy is favored by the data.
As is well known, there is mounting observational evidence that our Universe is presently dominated by two exotic forms of matter or energy. Cold, nonbaryonic dark matter (CDM), which accounts for of the critical mass density and whose leading particle candidates are the axions and the supersymmetric charginos and neutralinos, was originally proposed to explain the general behavior of galactic rotation curves that differ significantly from the one predicted by Newtonian mechanics. Later on, it was also realized that the same concept is necessary for explaining the evolution of the observed structure in the Universe from density inhomogeneities of the size detected by a number of Cosmic Microwave Background (CMB) experiments. Dark energy or quintessence, which accounts for of the critical mass density and whose leading candidates are a cosmological constant and a relic scalar field , has been inferred from a combination of astronomical observations which includes distance measurements of type Ia supernovae (SNe Ia) indicating that the expansion of the Universe is speeding up, CMB anisotropy data suggesting , and clustering estimates providing (see, e.g., revde () for current reviews).
Another interesting aspect related to these two dark components is that, unless some unknown symmetry in Nature prevents or suppresses a non-minimal coupling between them (see carroll () for a discussion), the dark matter and dark energy fields may interact between themselves, giving rise to the so-called models of coupled quintessence (see cq (); cq1 (); jmaia1 (); saulo () and Refs. therein). From the observational side, no piece of evidence has been so far unambiguously presented against such an interaction, so that a weak coupling still below detection cannot be completely excluded. In the cosmological context, models of coupled quintessence are capable of explaining the current cosmic acceleration, as well as other recent observational results cq (); jmaia1 (); saulo (). From the theoretical viewpoint, however, critiques to these scenarios do exist and are mainly related to the fact that in order to establish a model and study its cosmological consequences, one needs first to specify a phenomenological coupling between the cosmic components.
In this paper, instead of adopting the traditional approach, we follow the qualitative arguments used in Refs. wm (); alc (); ernandes (); prd () and deduce the interacting law from a simple argument about the effect of the dark energy decay on the dark matter expansion rate. As a step towards fundamental physics motivation for this phenomenological dark matter/dark energy interaction, we study how this class of coupled cosmologies can be interpreted in terms of classical scalar field dynamics. We also discuss current observational constrains on the interacting term . To this end we use two different samples of SNe Ia data, namely, the SNLS collaboration sample of 115 data points Astier06 () and the most recent SNe Ia compilation, the so-called Constitution set of 397 SNe Ia covering a redshift range from to cs (). Given the complementarity of SNe Ia data with other cosmological observables, a joint analysis with measurements of baryon acoustic oscillation (BAO) peak at 2df () and bao (), and measurements of the Hubble evolution svj () is also performed.
Ii interacting dark matter/dark energy
If dark energy and dark matter interacts, the energy density of this latter component will dilute at a different rate compared to its standard evolution, , where is the cosmological scale factor. Thus, the deviation from the standard dilution may be characterized by a constant , such that
where we have set the present-day value of the cosmological scale factor . By assuming that the radiation and baryonic fluids are separately conserved, the energy conservation law for the two interacting components reads , where or, equivalently (in our notation, ),
By combining the above equations and considering that the dark energy component is described by an equation of state parameter , one can show that the energy density of the dark energy component evolves as prd ()
where the integration constant is the present-day fraction of the dark energy density. Clearly, in the absence of a coupling with the CDM component, i.e., , the conventional non-interacting quintessence scenario is fully recovered. For and , we may identify (the current value of the vacuum contribution), and the above expression reduces to the vacuum decaying scenario recently discussed in Refs. wm (); alc (); ernandes ().
Fig. (1a) shows the ratio as a function of the redshift parameter [Eqs. (1) and (3)] for some selected values of and and . It is shown that this fraction is a very sensitive function of the interacting parameter . For instance, at , the dark energy contribution relative to for is about 75% larger than the same contribution for . This amounts to saying that for relatively large and positive values of the energy density of the pressureless component is not negligible at high-. For the sake of comparison, the CDM case ( and ) is also shown (solid line).
Iii scalar field description
In order to find the scalar field counterpart for the flat interacting dark matter/energy model considered here (that shares the same dynamics and temperature evolution law), we use the procedure originally proposed in Ref. jmaia () (see also jmaia1 () for a general description with arbitrary curvature).
First, let us define the parameter , which is just another way of writing the field equations
where , and are the energy densities of baryons, cold dark matter and dark energy, respectively, whereas stands for the dark energy pressure. Thus, any cosmology dynamically equivalent to the coupled quintessence scenario considered here must have the same . Following standard lines, we replace the energy density and pressure in Eqs. (4) and (5) by the corresponding scalar field expressions, i.e., and , where and are, respectively, the energy density and pressure associated with the coupled scalar field whose potential is .
By defining a new parameter with , we can manipulate the above equations to obtain (see jmaia () for more details)
which link directly the field and its potential with the related quantities of the coupled quintesence case. From Eq. (6a), one can show that in terms of the scale factor the field is given by
with and . From now on, we consider to be a constant, which is equivalent to impose the condition that the scalar field version mimics exactly the interaction rate of its phenomenological coupled counterpart jmaia (). From Eqs. (6b) and (8), we also obtain
where and .
Figure (1b) shows the potential obtained from a numerical combination of Eqs. (7)-(9) by assuming and . Four different cases are shown, namely, quintessence () and phantom () EoS with and without a separately conserved baryonic component . As one may also check, the potential shown in Fig. (1) reduces to
when , where , and are constants ernandes (). In general, we can see that the effect of of baryons on the potential , as well as a changing in the EoS values, produces only a small shift relative to the case in which the baryonic content participates of the interacting process or to the case of vacuum decay, so that the general potential for these scenarios belongs to the same class of potentials given analytically by Eq. (10) ref ().
Iv Current observational constraints
The description of the interacting scenario discussed in Sections II and III clearly shows that it comprises a multitude of cosmological solutions for different combinations of , and . In this Section we investigate observational bounds on the parametric spaces and from statistical analyses involving three classes of cosmological observations. In our analysis we fix from WMAP results Sperg07 (), a value in good agreement with the constraints derived from primordial nucleosynthesis nucleo ().
We use current distance measurements to SNe Ia, namely, the SNLS collaboration sample of 115 points Astier06 () and the most recent SNe Ia compilation, the so-called Constitution set (CS) cs (). The SNLS sample includes 71 high- SNe Ia in the redshift range and 44 low- events compiled from the literature but analyzed from the same manner as the high- sample. The CS of 397 SNe Ia covers a redshift range from to , including 139 SNe Ia at , and constitutes the largest SNe Ia luminosity distance sample currently available.
and is given by Eqs. (1), (3) and (4). In our analysis, we also use this ratio along with 9 determinations of the Hubble evolution , as given in Ref. svj () (for more details on statistical analysis involving these data sets we refer the reader to Ref. sa ()).
The results of our statistical analyses are shown in Fig. (2). As usual, the total likelihood is written as , where the function takes into account the three data sets discussed above. By marginalizing over (), we can quantify how much the plane () can be constrained by the data. The contour levels for these analyses are displayed, respectively, in Panels (2a) and (2b). In both Panels, dashed lines correspond to 1 and 2 values of the joint analysis involving SNLS + BAO + H() measurements whereas the solid ones to CS + BAO + H() data. Although the bounds on from the latter analysis (CS) are tighter than those from the former one (SNLS), the important point worth emphasizing here is that both negative and positive values for the interacting parameter are allowed by these analyses. Physically, this amounts to saying that not only an energy flow from dark energy to dark matter () is observationally allowed but also a flow from dark matter to dark energy () [see Eq. (1)]. At 95.4% (C.L.), we have found and and and , respectively, for the SNLS + BAO + H() and CS + BAO + H() joint analyses.
In Fig. (2c) we show a similar analysis by fixing the dark energy EoS at , which is fully equivalent to the vacuum decay scenario discussed in Refs. wm (); alc (). As shown in Ref. alc (), for the vacuum decay process the interacting parameter is restricted by thermodynamical arguments to be positive ref2 (). By considering this physical constraint on we have found at 95.4% (C.L.) and [SNLS + BAO + H()] and and [CS + BAO + H()], which is equivalent to the standard CDM case. If we drop the thermodynamical constraint (), we find, respectively, and and and at 95.4% (C.L.).
V Final Remarks
We have discussed some cosmological consequences of an alternative mechanism of cosmic acceleration based on a general class of coupled quintessence scenarios whose interaction term is deduced from the effect of the dark energy on the CDM expansion rate. The resulting expressions for the model are parameterized by a small parameter and have many of the previous phenomenological approaches as a particular case.
Following the procedure proposed in Refs. jmaia (); jmaia1 (), we have also discussed a scalar field description for these cosmologies and found that this class of models is identified with a coupled quintessence field with the potential given by Eqs. (7)-(9) [see also Fig. (1)]. From the observational point of view, we have investigated the current bounds on the interacting parameter from the most recent data of SNe Ia (SNLS and Constitution samples), baryon acoustic oscillation peak at (2dFGRS) and (SDSS), and the Hubble evolution . We have shown that negative and positive values for the interacting parameter are observationally allowed, which means that both an energy flow from dark energy to dark matter as well as a flow from dark matter to dark energy are possible. It is worth emphasizing that in our analysis the interacting parameter has been considered constant whereas in a more realistic case it must vary with redshift. The theoretical and observational consequences of this more realistic interacting scenario, as well as a full comparison with the case discussed in the present analysis, will appear in a forthcoming communication new ().
Acknowledgements.FEMC is grateful to CAPES - Brazil for financial support. EMB and JSA thank CNPq - Brazil for the grants under which this work was carried out.
- (1) V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000); P. J. E. Peebles and B. Ratra Rev. Mod. Phys. 75, 559 (2003); T. Padmanabhan, Phys. Rept. 380, 235 (2003); E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D15, 1753 (2006); J. S. Alcaniz, Braz. J. Phys. 36, 1109 (2006); J. A. Frieman, AIP Conf. Proc. 1057, 87 (2008). arXiv:0904.1832 [astro-ph.CO].
- (2) S.M. Carroll, Phys. Rev. Lett. 81 3067 (1998).
- (3) M. zer and M. O. Taha, Phys. Lett. B 171, 363 (1986); O. Bertolami, Nuovo Cimento Soc. Ital. Fis., B93, 36 (1986); L. Amendola, Phys. Rev. D62 043511 (2000); W. Zimdahl, D. Pavón and L. Chimento, Phys. Lett. B 521, 133 (2001); L. Amendola and D.Tocchini-Valentini, Phys. Rev. D 64, 043509 (2001); I. L. Shapiro, J. Sola and H. Stefancic, JCAP 0501, 12 (2005); J. D. Barrow and T. Clifton, Phys. Rev. D 73, 103520 (2006); Q. Wu, Y. Gong, A. Wang and J. S. Alcaniz, Phys. Lett. B 659, 34 (2008); M. Quartin et al., JCAP 0805, 007 (2008).
- (4) I. Waga, Astrophys. J. 414, 436 (1993); J. M. Overduin and F. I. Cooperstock, Phys. Rev. D 58, 043506 (1998); S. Carneiro et al., Phys. Rev. D 74, 023532 (2006); L. Amendola, G. C. Campos and R. Rosenfeld, Phys. Rev. D 75, 083506 (2007); H. A. Borges, S. Carneiro and J. C. Fabris, Phys. Rev. D 78, 123522 (2008).
- (5) J. S. Alcaniz and J. M. F. Maia, Phys. Rev. D 67, 043502 (2003).
- (6) S. Carneiro et al., Phys. Rev. D 77, 083504 (2008).
- (7) P. Wang and X. Meng, Class. Quant. Grav. 22, 283 (2005).
- (8) J. S. Alcaniz and J. A. S. Lima, Phys. Rev. D 72, 063516 (2005).
- (9) F. E. M. Costa, J. S. Alcaniz and J. M. F. Maia, Phys. Rev. D 77, 083516 (2008)
- (10) J. F. Jesus et al. Phys. Rev. D 78, 063514 (2008).
- (11) P. Astier et al., Astron. Astrophys. 447, 31 (2006).
- (12) M. Hicken et al., arXiv:0901.4804 [astro-ph.CO].
- (13) W. J. Percival et al., MNRAS 381, 1053 (2007).
- (14) D. J. Eisenstein et al., Astrophys. J. 633, 560 (2005).
- (15) J. Simon, L. Verde and R. Jimenez, Phys. Rev. D 71, 123001 (2005).
- (16) J. M. F. Maia and J. A. S. Lima, Phys. Rev. D 65, 083513 (2002).
- (17) Double exponential potentials of the type (10) have been considered in the literature as viable examples of quintessence scenarios (see, e.g, dp ()). As discussed in Ref. dp1 (), a field potential given by the sum of two exponential terms is also motivated by dimensional reduction in M-theory with implications for the late-time accelerating behavior of the cosmic expansion.
- (18) D. N. Spergel et al., Astrophys. J. Supl. 170, 377 (2007).
- (19) G. Steigman, Ann. Rev. Nucl. Part. Sci, 57, 463 (2007).
- (20) As discussed in Ref. doran (), the BAO parameter can be considered as a good approximation for a large class of dark energy models. In our analysis, we do not use this quantity because it also uses the approximation that the sound travel distance is , which fails if dark energy has a significant contribution in the past. According to Fig. (1a), in our interacting scenarios it may happen for relatively large and positive values of (see saulo () for a discussion).
- (21) F. E. M. Costa and J. S. Alcaniz, in preparation (2009).
- (22) A. A. Sen and S. Sethi, Phys. Lett. B 532, 159 (2002).
- (23) L. Jarv, T. Mohaupt and F. Saueressig, JCAP 0408, 016 (2004).
- (24) M. Doran, S. Stern and E. Thommes, JCAP 0704, 015 (2007).
- (25) Note that the thermodynamical constraint seems not to apply to a general quintessence case (). See, e.g., Refs. t1 () and t2 () for a discussion.
- (26) S. H. Pereira and J. F. Jesus, Phys. Rev. D 79, 043517 (2009).
- (27) D. Pavon and B. Wang, Gen. Rel. Grav. 41, 1 (2009).
- (28) T. Padmanabhan and T.R. Choudhury, MNRAS 344, 823 (2003); S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 70, 043531 (2004); J.S. Alcaniz, Phys. Rev. D 69, 083521 (2004); T.R. Choudhury and T. Padmanabhan, Astron. Astrophys. 429, 807 (2005); Z. K. Guo, Z. H. Zhu, J.S. Alcaniz and Y. Z. Zhang, Astrophys. J. 646, 1 (2006); L. Samushia and B. Ratra, Astrophys. J. 650, L5 (2006); F. C. Carvalho et al., JCAP 0809, 008 (2008).