# Unified model of Baryonic matter and dark components

###### Abstract

We investigate an interacting two-fluid cosmological model and introduce a scalar field representation by means of a linear combination of the individual energy densities. Applying the integrability condition to the scalar field equation we show that this ”exotic quintessence” is driven by an exponential potential and the two-fluid mixture can be considered as a model of three components. These components are associated with baryonic matter, dark matter and dark energy respectively. We use the determination of the redshift dependence of the Hubble parameter to constrain the current density parameters of this model. With the best fit density parameters we obtain the transition redshift between non accelerated and accelerated regimes and the time elapsed since the initial singularity . We study the perturbation evolution of this model and find that the energy density perturbation decreases with the cosmological time.

## I Introduction

Astrophysical data suggests that the Universe is accelerating Perlmutter:1998np (); Spergel:2003cb (). This acceleration may be explained by very different models, among them, the simplest one is the CDM Padmanabhan:2002ji (); Sahni:1999gb (). It assumes a cosmological constant arising from the energy density of the zero point fluctuations of the quantum vacuum and cold dark matter in form of pressureless dust. While it fits rather well all the observational constraints, the small positive value of the energy density of the vacuum remains as an explanatory challenge for physics today. See also Refs. Triay:2005jr (); Mersini-Houghton:2006ue () for objections to this interpretation. The next step is to propose a dark energy component that may vary with time and that is generally modeled by a scalar field. Most of these models assume that dark matter and scalar field components evolve independently. Again, this is not the solution because in the analysis of these models through SNIa or WMAP data, best fit models for one set data alone is usually ruled out by the other set at a large confidence limit Jassal:2005qc (). As these conclusions are valid only for standard models, where dark energy and dark matter are decoupled, many papers have been devoted to interacting models Amendola:2006ku (). Certain models conceive the interaction as a time-variable dark mass, evolving with an inverse power law potential or an exponential potential, Farrar:2003uw (), Amendola:2002kd (), Neupane:2007mu (), Axenides:2004kb ().

In this paper we clarify this point and establish an exactly solvable model with a smooth transition from a matter dominated phase to a period of accelerated expansion. We introduce an interacting two-fluid cosmological model and investigate the effects of imposing the integrability condition of the whole equation of conservation. It forces the dynamic of the model to be governed by a modified Friedmann equation with three components. One of them, associated with an exponential potential, drives an exotic scalar field, (exotic quintessence). Here, we show that the problem of an accelerating universe can be realized in a comparatively simple manner within the framework of general relativity. Finally, the perturbation evolution of the model is investigated.

Our paper is organized as follows. In section II we consider the general interacting two-fluid cosmological model and introduce the exotic quintessence. There, we obtain the evolution equation of the exotic field, find their implicit solutions, build the modified Friedmann equation and show the asymptotic behavior of the scale factor by using stability analysis. In section III we introduce baryonic matter, dark matter and dark energy components and show that dark components satisfy separately an effective equation of conservation with variable equation of state. In section IV we find confidence regions for the parameters of the model by using the Hubble function H(z) data, the age of universe and the redshift of the transition from non-accelerated to accelerated regime for the best fit model. In section V we present the equations governing the perturbations of the model. In section VI we express the conclusions.

## Ii Exotic quintessence

Cosmological models are based in the Einstein equations of the gravitational field where the source includes different kinds of matter known, for instance protons, neutrons, photons, neutrinos, etc., as well as non-relativistic non-baryonic cold dark matter and dark energy. We embark in a less ambitious project by considering a model consisting of two perfect fluids with an energy momentum-tensor . Here , where and are the energy density and the equilibrium pressure of fluid and is the four-velocity. Assuming that the two fluids interact between them in a spatially flat homogeneous and isotropic Friedmann-Robertson-Walker (FRW) cosmological model, the Einstein’s equations reduce to two algebraically independent equations:

(1) | |||

(2) |

where a(t) is the FRW scale factor and is the Hubble expansion rate. We introduce an equation of state for each fluid component and for simplicity we assume that are constants and . This simplified model leads to a reduction of the number of fundamental parameters required to describe observations. It can be considered as an advantage from the computational point of view. We choose units such as the gravitational constant is set to and .

The whole equation of conservation (2) shows the interaction between both fluid components allowing the mutual exchange of energy and momentum, meaning that, there will be no local energy-momentum conservation for these fluids separately. Then, we assume an overall perfect fluid description with an effective equation of state, , where and . So that, from Eqs. (1)-(2) we get

(3) |

To avoid an eventual super acceleration of the universe, which could lead to a “big rip” singularity, we choose and crossing (). Hence, and there is no big rip singularity. These models can be investigated by introducing a scalar field representation of the interacting two-fluid mixture

(4) |

with . The dynamic equation for the scalar field is obtained from the equation of conservation (2)

(5) |

It can be integrated by setting the interaction between the two fluids by

(6) |

where is a new constant parameter of the model. Integrating Eq. (6), we find that the energy density of the second fluid can be associated with an exponential potential

(7) |

where is a positive integration constant and , are the present values of Hubble constant and scalar field.

From Eqs. (4), (5) and (7) we obtain the total energy density and pressure of the fluid mixture and the dynamical equation for the scalar field

(8) | |||

(9) | |||

(10) |

with . These equations are different than the conventional ones describing quintessence, in contrast, define an exotic scalar field. When the interacting two-fluids system is related to the scalar field in the form and , with equations of state and , meaning that (stiff matter) and (vacuum energy), the exotic scalar field reduces to quintessence. Then, due to the interactions between the two-fluid components the energy-momentum tensor conservation of the system, as a whole, is equivalent to the Klein-Gordon equation. For any other interacting two-fluid mixture, the cosmological model contains an exotic quintessence field driven by an exponential potential.

Using the integrability condition (6) in the field equation (5), its first integral is given by

(11) |

and

(12) |

where is an arbitrary integration constant, is the redshift parameter and is the present scale factor. This model is finally closed when the Eq. (11) is inserted into the energy density (8) 2c () . Hence, the Friedmann equation (1) reads

(13) |

As a consequence of the linear term in the expansion rate , this equation can be seen as amodified Friedmann equation. Its solution gives the scale factor and the model we propose, containing exotic quintessence, could be formally solved.

In order to obtain the asymptotic behavior of the scale factor it will be useful to find the constant solutions of the dynamical equation for the overall equation of state

(14) |

and investigate their asymptotic stability. This equation has two stationary solutions and (Eq. (3) excludes the solution ). Assuming that , we find that is an unstable solution while becomes asymptotically stable. Essentially, the evolution of the geometry is dictated by , meaning that the universe begins to evolve from an unstable phase as it were dominated by the first fluid at early times, , and ends in a stable expanding phase dominated by the exponential potential, . The latter becomes an expanding accelerated phase when the slope of the potential satisfy the inequality .

## Iii Dynamic of baryonic and dark components

The integrability condition (6) can be considered as an effective equation of conservation for the second fluid. This allows us to identify with the energy density of the dark energy component and with its effective equation of state. Expressing the latter in term of the exotic field we get

(15) |

and the relation

(16) |

linking the overall and dark energy equations of state. Also, for convenience we write with

(17) |

The linear term in the Eq. (13) is adequate to describe non-relativistic non-baryonic cold dark matter components whose energy-momentum tensor is approximately dust-like. Finally, the baryonic matter is introduced by setting in the third term of the Eq. (13). Making these identifications the modified Friedmann equation (13) becomes

(18) |

where

(19) | |||

(20) | |||

(21) |

are the present dark matter, dark energy and baryonic matter density parameters respectively. As these are constrained according to , we conclude that

(22) |

During the accelerated epoch of the universe , the SEC is violated and leads to . Hence, one finds that the conditions and are consistent with the requirement of having positive density parameters.

The exotic quintessence is essentially based on the integrability condition (6) under which the conservation equation (2) and the exotic field equation (5) can be integrated. Actually, from the above condition along with Eq. (11), we get for an expanding universe. Then, is a Liapunov function and the solution , of Eq. (6), is asymptotically stable. Also, from Eqs. (4), (7) and (11) the energy density of the first fluid has a vanishing limit in the remote future. Hence, this general model is viable and it does not contradict basic cosmological conjectures.

From Eqs. (19) and (21) we can express the parameter and the integration constant in terms of the present density parameters ,

(23) | |||

(24) |

besides,

(25) |

Finally the original problem of the interacting two-fluid mixture governed by system equations (1)-(2) is equivalent to an effective model with a “three-fluid” mixture. So that, the effective dynamical equations of our model read

(26) | |||

(27) | |||

(28) | |||

(29) |

where

(30) | |||

(31) | |||

(32) |

are the effective energy densities of dark and baryonic components and

(33) |

is the effective equation of state of the dark energy. Also, we find the following relation

(34) |

between the effective equations of state of dark components. So that, the knowledge of them determines the remaining ones including . To obtain the above effective dynamical equations of the model we have taken into account that . Hence, we have identified the equation of conservation (6) with (28) to express, after using (8), the effective for dark energy in terms of the observed density parameters, the present Hubble expansion rate and the redshift parameter.

In our model the stationary solutions at early times and at late times along with Eqs. (16) and (34) lead to the stationary solutions , and , . As the observed density parameter satisfies the condition , then by using Eq. (23). In this case, we find that is an unstable solution at early times and becomes asymptotically stable at late times. Here, the evolution of the geometry represents a universe that begins to evolve as it were matter dominated at early times and ends in an accelerated phase dominated by the dark energy component. As the accelerated epoch begins at or , the corresponding redshift is given by the expression

(35) |

On the other hand, the Eq. (33) allow us to find the elapsed time from the creation

(36) |

where we have used that the scale factor behaves as , at early times.

## Iv Observational Constraints

In Ref. Lazkoz:2007zk () it was used the recently published Hubble function data [SV&J(2005)]Simon:2004tf (), extracted from differential ages of passively evolving galaxies. This is interesting for, among other reasons, the function is not integrated over, in contrast to standard candle luminosity distances or standard ruler angular diameter distances Since the Hubble parameter depends on the differential age of the Universe as a function of z in the form , it can be directly measured through a determination of dz/dt. In the procedure of calculating the differential ages, Simon et al. have employed the new released Gemini Deep Deep Survey Abraham:2004ra () and archival data Treu:2001hq (), Nolan:2003bt () to determine the 9 numerical values of in the range , and their errors [see Table 1]. These data will be inserted in our Eqs. (17) and (18) to derive restrictions on the range of possible values for the density parameters (19)-(21).

0.09 | 69 | |
---|---|---|

0.17 | 83 | |

0.27 | 70 | |

0.40 | 87 | |

0.88 | 117 | |

1.30 | 168 | |

1.43 | 177 | |

1.53 | 140 | |

1.75 | 202 |

We adopt the prior for . It is exactly the mean value of the results from the Hubble Space Telescope key project Freedman:2000cf () and consistent with the one from WMAP 3-year result Spergel:2006hy ().

The parameters of the model can be determined by minimizing the function

(37) |

where

(38) |

is the predicted value for the Hubble parameter, obtained from Eq. (18), and is calculated from Eqs. (12) and (17).

is the observed value of at the redshift , is the corresponding uncertainty, and the summation is over the 9 observational data points at redshift Wei:2006ut (). Also, we adopt the prior Astier:2005qq ().

We find a local minimum of for and () and present the observational data in Fig. 1 with error bars and the theoretical line corresponding to the best fit parameter. The Fig. 2 shows that the Universe begins to accelerate on . A similar result can be obtained from CDM flat cosmology when the density parameters are and Copeland:2006wr (). In Fig. 3 we plot confidence regions in the - plane. The true values of those parameters are inside the inner ellipse or between both ellipses with or percent of probability respectively. In Fig. 4 we plot the equations of state for dark matter , dark energy and overall fluid . According to the previous stability analysis predictions we see that the asymptotic behavior of the overall state parameter varies from (near cold matter behavior) in the far past to in the far future (dark energy behavior). In Fig. 5 we show the energy densities corresponding to the baryonic and dark components in terms of the redshift .

We calculate the age of the universe (36) with the best fit parameters and find that . In Fig. 6 we plot the time elapsed (in units of ), since the initial singularity to present days, for our model and the flat CDM model, as a function of the matter density. We also show the border coming from the bound of the oldest stellar ages. The age of the universe in this coupled scenario tends to be much higher when compared with the ACDM case Franca:2003zg ()

## V Linear perturbations

Cosmological models with two interacting fluids have been investigated with the purpose of describe the evolution of dark components. There the energy momentum tensor of the interacting components is not separately conserved. Usually these cosmological model are presented with interacting matter species have a non constant equation of state parameter koi () or with DE having a constant equation of state parameter coupled to DM ladw (),diego (). However, in the model we are investigating the interaction between the two fluids is setting by Eq. (6). For this choice the field equation can be integrated, generalizing the case of quintessence driven by the exponential potential. Then, it will be interesting to investigate the evolution of the density perturbation.

In the synchronous gauge the line element is given by:

(39) |

where the commoving coordinate are related to the proper time and position by , and is the metric perturbation. The scalar mode of is described by the two fields and in the Fourier space,

(40) |

with . The Einstein equations to linear order in k-space, expressed in terms of and , are given by the following four equations Ma ():

(41) | |||

(42) | |||

(43) | |||

(44) |

Here, the quantities and are defined as , and denotes the traceless component of the tensor . In addition, is the divergence of the fluid velocity and means .

Let us consider a fluid moving with a small coordinate velocity , then, can be treated as a perturbation of the same order as energy density, pressure and metric perturbations. Hence, to linear order in the perturbations, the energy-momentum tensor, with vanishing anisotropic shear perturbation , is given by

(45) | |||

(46) | |||

(47) |

For a fluid with equation of state , the perturbed part of energy-momentum conservation equations in the k-space leads to the equations

(48) | |||

(49) |

where and . Besides, using equations (41), (43), (45) and (47) we arrive at

(50) |

We have showed that our interacting two-fluid model can be associated with an overall perfect fluid description based in an effective equation of state . Hence, we investigate the asymptotic regimes at early and late times assuming nearly constant equations of state and respectively.

At early time, when the overall fluid has , the effective fluid perturbations evolve similar to those of ordinary dust with , and from Eqs. (48-50) we obtain

(51) |

and , where and are arbitrary integration constants. In this dust dominated era the perturbation grows as showing an initial unstable phase and compatible with the observation that the primordial universe would have tiny perturbations which seed the formation of structures in the later universe.

At late times, we are interested to find the evolution of the linear scalar perturbations for any mode . To this end we write the second order differential equation for the density perturbation and the first order differential equation for the divergence of the fluid velocity , evaluating them on the asymptotically stable equation of state . In this case, from Eqs. (48-50) we get:

(52) | |||

(53) |

Taking into account that in the late time regime the scale factor behaves as we can calculate the conformal time , and

(54) | |||

(55) | |||

(56) |

From Eqs. (52) and (53) the perturbation evolution becomes mode dependent with the term, and for low energy modes their solutions can be obtained assuming a power law dependence of the perturbations with the scale factor, and . In this case the approximate solutions for are given by

(57) | |||

(58) |

where , and are integration constants while is a function of and . This shows that the coupling to in Eq. (52) can be neglected for all scales we are interested. Finally, expressing the Eq. (52) in term of conformal time we get

(59) |

The general solution of the latter equation in terms of the Bessel functions is

(60) |

with

(61) |

At late times, it can be approximated by the two first terms of the Eq. (58) showing that the energy density perturbation decreases for large cosmological times for modes satisfying the condition . For high energy modes, , the perturbation

(62) |

decreases but slowly that the low energy modes. This results can be understood writing the Eq. (52) as the equation of motion of a dissipative mechanical system by using the analogy with the classical potential problem

(63) |

where

(64) | |||

(65) | |||

(66) |

The potential has an extreme at , it is maximum for or a minimum for . On the other hand, assuming that the perturbation depends on the scalar factor in the form , we find that . Hence, for any mode the perturbation begins to grow at early times for , while at late times for , the function inside the square bracket in Eq. (63) is a Liapunov function and the perturbation decreases asymptotically reaching in the limit .

## Vi Conclusions

We have shown an interacting two-fluid cosmological model that allows us to reproduce the accelerated behavior of our universe and its probable age . The model gives rise to an exotic scalar field dubbed exotic quintessence which reduces to quintessence when one fluid is associated with stiff matter and the other with vacuum energy. Setting the interaction between the two fluids by Eq. (6), the field equation is integrated, generalizing the case of quintessence driven by the exponential potential and, the equation governing the scale factor (13) looks like a modified Friedmann equation.

We have obtained the evolution equation for the overall equation of state of the model and showed the asymptotic behavior of the scale factor i.e., the universe begins from an unstable phase dominated by the first fluid, , and ends in a stable expanding phase dominated by the exponential potential, . The latter becomes accelerated when the exponential potential slope satisfy . Setting , the scale factor interpolates between pressureless matter and dark energy phases.

Using the Hubble function data from Table 1 we minimize the function (37), obtaining the best fit densities parameters and with a reduced . These results are consistent with those found in the literature, see for instance Ref. Amendola:2006dg () for null coupling () or with the result obtained in Ref. Feldman:2003nu (), () through mean relative peculiar velocity measurements for pairs of galaxies. With our best densities parameters and the priors for and , we obtain the theoretical function, plotted together with the SVJ(2005) experimental data in Fig. 1. In Fig. 2 we plot the acceleration of the model as a function of the redshift and find the transition from the non accelerated phase to the accelerated one around Copeland:2006wr (). Our value agrees with the result obtained by a nearly model independent characterization of dark energy properties as a function of redshift, (,Daly:2006ax ()). The problem of why an accelerated expansion should occur now in the long history of the universe seems to be naturally dressed in our model. Considering the age of the universe, we take into account that the age of the oldest stellar objects have been constrained for instance, by using a distance-independent method Jimenez:1996 (), ( for Globular clusters in the Milky Way) and the white dwarfs cooling sequence method Hansen:2002ij () ( for the globular cluster M4). Then, the age of universe needs to satisfy the lower bound . This condition is fulfilled by our model with , as it can be seen in Fig. 5.

The energy density perturbation of the model grows in the first stage of the universe showing that initial instabilities in the primordial universe could leads to the formation of structure in the later universe. At late times we have found a Liapunov function which indicates that the perturbation decreases asymptotically reaching in the limit .

## Acknowledgments

This work was partially supported by the University of Buenos Aires and Consejo Nacional de Investigaciones Científicas y Técnicas under Projects X224 and 5169.

## References

- (1) S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) [arXiv:astro-ph/9812133].
- (2) D. N. Spergel et al. [WMAP Collaboration], “First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Astrophys. J. Suppl. 148, 175 (2003) [arXiv:astro-ph/0302209].
- (3) T. Padmanabhan, Phys. Rept. 380, 235 (2003) [arXiv:hep-th/0212290].
- (4) V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000) [arXiv:astro-ph/9904398].
- (5) R. Triay, Int. J. Mod. Phys. D 14 (2005) 1667 [arXiv:gr-qc/0510088].
- (6) L. Mersini-Houghton, arXiv:gr-qc/0609006.
- (7) H. K. Jassal, J. S. Bagla and T. Padmanabhan, Phys. Rev. D 72, 103503 (2005) [arXiv:astro-ph/0506748].
- (8) L. Amendola, M. Gasperini and F. Piazza, Phys. Rev. D 74, 127302 (2006) [arXiv:astro-ph/0610574]. L. Amendola, G. Camargo Campos and R. Rosenfeld, “Consequences of dark matter - dark energy interaction on cosmological Phys. Rev. D 75, 083506 (2007) [arXiv:astro-ph/0610806]. L. P. Chimento, A. S. Jakubi, D. Pavon and W. Zimdahl, Phys. Rev. D 67, 083513 (2003) [arXiv:astro-ph/0303145]. L. P. Chimento, Phys. Lett. B 633, 9 (2006) [arXiv:gr-qc/0503049]. L. P. Chimento and D. Pavon, Phys. Rev. D 73, 063511 (2006) [arXiv:gr-qc/0505096]. E. Majerotto, D. Sapone and L. Amendola, arXiv:astro-ph/0410543.
- (9) G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 (2004) [arXiv:astro-ph/0307316].
- (10) L. Amendola, Mon. Not. Roy. Astron. Soc. 342, 221 (2003) [arXiv:astro-ph/0209494].
- (11) I. P. Neupane, [arXiv:0706.2654 [hep-th]].
- (12) M. Axenides and K. Dimopoulos, JCAP 0407 (2004) 010 [arXiv:hep-ph/0401238].
- (13) L.P Chimento and R. Lazkoz Phys. Lett. B639, 591-595 (2006).
- (14) L.P. Chimento Class. and Quantum Grav. 15, 965-974 (1998).
- (15) R. Lazkoz and E. Majerotto, arXiv:0704.2606 [astro-ph].
- (16) J. Simon, L. Verde and R. Jimenez, Phys. Rev. D 71, 123001 (2005) [arXiv:astro-ph/0412269].
- (17) R. G. Abraham et al., “The Gemini Deep Deep Survey: I. Introduction to the Survey, Catalogs and Astron. J. 127, 2455 (2004) [arXiv:astro-ph/0402436].
- (18) T. Treu, M. Stiavelli, P. Moller, S. Casertano and G. Bertin, “The properties of field elliptical galaxies at intermediate redshift. II: Mon. Not. Roy. Astron. Soc. 326, 221 (2001) [arXiv:astro-ph/0104177].
- (19) P. L. Nolan, W. F. Tompkins, I. A. Grenier and P. F. Michelson, Astrophys. J. 597, 615 (2003) [arXiv:astro-ph/0307188].
- (20) W. L. Freedman et al., Astrophys. J. 553 (2001) 47 [arXiv:astro-ph/0012376].
- (21) D. N. Spergel et al. [WMAP Collaboration], arXiv:astro-ph/0603449.; L. Page et al. [WMAP Collaboration], arXiv:astro-ph/0603450.; G. Hinshaw et al. [WMAP Collaboration], arXiv:astro-ph/0603451.; N. Jarosik et al. [WMAP Collaboration], arXiv:astro-ph/0603452.
- (22) H. Wei and S. N. Zhang, Phys. Lett. B 644 (2007) 7 [arXiv:astro-ph/0609597].
- (23) P. Astier et al. [The SNLS Collaboration], Astron. Astrophys. 447 (2006) 31 [arXiv:astro-ph/0510447].
- (24) E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15 (2006) 1753 [arXiv:hep-th/0603057].
- (25) U. Franca and R. Rosenfeld, Phys. Rev. D 69, 063517 (2004) [arXiv:astro-ph/0308149].
- (26) T. Koivisto, Phys. Rev.D 72, 043516 (2005).
- (27) W. Zimdahl, D. Pavon, and L. P. Chimento, Phys. Lett. B521, 133 (2001); L. P. Chimento, A. S. Jakubi, D. Pavon, and W. Zimdahl, Phys. Rev. D 67, 083513 (2003).
- (28) G. Olivares, F. Atrio-Barandela, and D. Pavon, Phys. Rev.D 74, 043521 (2006).
- (29) C. P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995).
- (30) L. Amendola, G. Camargo Campos and R. Rosenfeld, Phys. Rev. D 75 (2007) 083506 [arXiv:astro-ph/0610806].
- (31) H. A. Feldman et al., Astrophys. J. 596 (2003) L131 [arXiv:astro-ph/0305078].
- (32) R. A. Daly and S. G. Djorgovski, arXiv:astro-ph/0609791.
- (33) R. Jimenez, P. Thejll, U. Jorgensen, J. MacDonald and B. Pagel, MNRAS, 282, 926 (1996).
- (34) B. M. S. Hansen et al., Astrophys. J. 574 (2002) L155 [arXiv:astro-ph/0205087]. H. B. Richer et al., Astrophys. J. 574 (2002) L151 [arXiv:astro-ph/0205086].