Advances in High Energy Physics

Volume 2014 (2014), Article ID 878092, 13 pages

http://dx.doi.org/10.1155/2014/878092

Research Article

Interacting Quintessence Dark Energy Models in Lyra Manifold

M. Khurshudyan,¹ J. Sadeghi,² R. Myrzakulov,³ Antonio Pasqua,⁴ and H. Farahani⁵

¹Department of Theoretical Physics, Yerevan State University, 1 Alex Manookian, 0025 Yerevan, Armenia

²Department of Physics, Islamic Azad University, Ayatollah Amoli Branch, P.O. Box 678, Amol, Iran

³Department of Physics, Eurasian International Center for Theoretical Physics, Eurasian National University, Astana 010008, Kazakhstan

⁴Department of Physics, University of Trieste, Via Valerio, No. 2, 34127 Trieste, Italy

⁵Department of Physics, Mazandaran University, P.O. Box 47416-95447, Babolsar, Iran

Received 8 April 2014; Revised 8 August 2014; Accepted 13 August 2014; Published 1 September 2014

Academic Editor: Sally Seidel

Copyright © 2014 M. Khurshudyan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³.

Abstract.

We consider two-component dark energy models in Lyra manifold. The first component is assumed to be a quintessence field while the second component may be a viscous polytropic gas, a viscous Van der Waals gas, or a viscous modified Chaplygin gas. We also consider the possibility of interaction between components. By using the numerical analysis, we study some cosmological parameters of the models and compare them with observational data.

1. Introduction

Observations of high redshift type Supernovae Ia (SNeIa) [1–3] reveal the accelerated expansion of our universe, whose nature is not exactly clear until now. It is found that the density of matter is very much less than critical density [4]. Moreover, cosmic microwave background (CMB) radiation anisotropies observations indicate that the universe can be considered flat and the total energy density is very close to the critical value () [5]. Based on the experimental data, a component of the energy dubbed as dark energy is thought to be responsible for the physics of the accelerated expansion but it seems that it is not alone in the universe, so the mysterious matter component which is called dark matter should also exist. Dark energy can be described by a pressure sufficiently negative in order to drive the acceleration of the universe and by positive energy density. There are several different models proposed to explain the nature of dark energy. The cosmological constant is the simplest model which can be considered, but in presence of many research papers in these fields, the origin of dark energy and dark matter is still unknown, and the possible connection between them is also unknown as well as real role of the components in the history of the universe. This situation gives a lot of freedom to researchers and possibility of some simulations. The cosmological constant faced two main problems, that is, the absence of a fundamental mechanism which sets the cosmological constant zero or very small value (which is known as fine-tuning problem) and the problem known as cosmological coincidence problem, which asks why we are living in an epoch in which the densities of dark energy and matter are comparable. One of the interesting ways to solve the above mentioned problems is to consider interactions between components [6]. From observational point of view, no piece of evidence has been so far presented against such interactions. Indeed, possible interactions between the components of universe have been discussed in recent years. It is found that a suitable interaction can help to alleviate the coincidence problem. Different interacting models of dark energy have been investigated [7–14].

Alternative models of dark energy suggest a dynamical form of dark energy, which, at least in an effective level, can originate from a variable cosmological constant [15, 16] or from various fields, such as a canonical scalar field [17–21] (quintessence), a phantom field [20–22], or quintom [23–36]. Finally, interesting attempts to probe the nature of dark energy according to some basic quantum gravitational principles are the holographic dark energy paradigm [37–45] and agegraphic dark energy models [46, 47]. Among them a quintessence model is interesting in this paper as a component of dark energy. In that case, the dark energy may be dynamical approaching zero with time, or it may be slowly increasing. It is now dominating the universe because the reduction of mass and radiation energy density with the scale factor (which gives some information about the size of the universe) is greater than the decrease in dark energy density in the present epoch. In general, we would like the quintessence field to be decreasing with the scale factor and time at a smaller rate than the mass energy so that it will become dominant at redshifts less than one. The quintessence field has the property of being very weakly coupled to baryons but contributing a negative pressure to the equation of state. In the past it had a small contribution but with time it has decreased less quickly with the scale factor than the matter and radiation densities and is dominant now.

For the dark energy component we consider several models in this paper, including viscosity. Indeed, bulk viscosity is added to obtain more realistic models. However, viscous pressure can itself play the role of an agent that drives the present acceleration of the universe [48].

One of interesting dark energy models is the polytropic gas which was proposed to explain the accelerated expansion of the universe [49]. It was shown that the polytropic gas model in the presence of interaction can behave as phantom field [50]. It was pointed out that a polytropic scalar field can be reconstructed according to the evolutionary behaviors of the holographic and new agegraphic dark energy densities. The validity of the generalized second law of thermodynamics was also examined for the polytropic gas model in [51]. Another interesting model of dark energy may be Van der Waals gas which could be accounted as a fluid with unusual EoS or could be thought of as a fluid satisfying more general form of EoS; that is, [52]. There are also some important models to describe dark energy based on Chaplygin gas equation of state which were recently considered by several papers such as [53–56] and yield good agreement with observational data.

On the other hand, the Lyra geometry provides one of the possible alternatives in modification of the cosmological models. As we know the modification of the gravitational theory has long been famous, but the late-time cosmological acceleration caused more research in this field [57]. Now, we like to consider a universe filled with a two-component dark energy in Lyra manifold with possibility of interaction between components. The first component is assumed to be quintessence, while we have several choices for the second component such as viscous polytropic gas, viscous Van der Waals gas, or viscous Chaplygin gas. We suggest these as toy models to describe universe and compare our results with observational data to choose one of them as the best model.

This paper is organized as follows. In Section 2, we introduce our models. In Section 3, we recall the main properties of field equations. In Section 4, we give numerical results corresponding to constant . In Section 5, we give numerical results corresponding to varying . In Section 6, we obtain some observational constraints. Finally, in Section 7, we write the conclusions of this paper.

2. The Models

One of the well studied dark energy models is the quintessence model [17, 18], which is a scalar field model described by a field and a potential. It represents the simplest scalar field scenario without having theoretical problems such as the appearance of ghosts and Laplacian instabilities. The energy density and the pressure of the quintessence scalar field model are given, respectively, by

Canonical scalar field is not the unique solution. We can generalize it as follows [58]:

In the case of , (2) transform to the canonical scalar field model with rescaling of the field. Below, we would like to consider an interaction term between dark energy and dark matter described by

where and are positive constants with a typical value of 0.01–0.03. Nature of the interaction between dark energy and dark matter is not clear. If we believe that it has a quantum origin, then an absence of the final theory of quantum gravity leaves this question as an open problem. However if we believe that the link existing between components is due to the same origin of the dark energy and dark matter, then this approach does not give any exact solution, because the nature of the two components is not formulated and it is another open problem. Therefore, only phenomenological assumption is an appropriate approach. For the dark matter model we will consider once a viscous modified Chaplygin gas with the following equation of state (EoS):

where , , and are constants (with in general relativity).

For the second model we will use viscous polytropic fluid with EoS given by

where is the polytropic index and represents the viscous coefficient.

In the third model we would like to consider interaction between quintessence dark energy and a viscous Van der Waals gas of the general form:

where and are constants. Furthermore, we will consider two regimes: (1) is a numerical constant; (2) is a function of the cosmic time ; therefore it is a varying quantity. In particular, we choose the following form for the time-varying :

where is a positive constant and is the potential of the field which we consider as follows:

where is a constant parameter.

We investigate the behavior of the cosmological parameters like the Hubble parameter , deceleration parameter , and EoS parameters of the quintessence dark energy and an effective two-component fluid. Moreover, we perform stability analysis via the squared speed of the sound , which is defined as follows:

where and are, respectively, the pressure and the energy density of the effective fluid given by

where refers to one of the viscous fluids described above. We will finish our paper with the results obtained from observational constraints. Consideration of the statefinder analysis, different forms of interaction terms, and different and varying viscosity is possible and an interesting research. We hope to approach this question in future with forthcoming articles.

3. The Field Equations

The field equations governing our model are given by

Considering the content of the universe to be a perfect fluid, we have

where is the 4-velocity of the comoving observer, satisfying the relation . Let be a time-like vector field of displacement; then

where is a function of time alone and the factor is inserted in order to simplify the writing of all the following equations. By using FRW metric for a flat universe given by

field equations can be reduced to the following Friedmann equations:

where is the Hubble parameter, the dot stands for differentiation with respect to the cosmic time , , and represents the scale factor. The and parameters are the usual azimuthal and polar angles of spherical coordinates, with and . The coordinates () are called comoving coordinates.

The continuity equation is given by

The continuity equation given in (16) can be also rewritten in the compact form:

Comparing (16) and (17) we obtain the following link between and :

In order to introduce an interaction between dark energy and dark matter, we should mathematically split (17) into the two following equations:

The cosmological parameters of our interest are the EoS parameter of each fluid component , the EoS parameter of composed fluid,

and the deceleration parameter , which can be written as follows:

where index refers to the first components which is viscous modified Chaplygin gas or viscous polytropic fluid and index refers to the quintessence scalar field. A differential equation describing dynamics of the DE after some mathematics can be rewritten as

Taking into account the form of the varying from (7) for the Hubble parameter we will have

Hereafter, we will consider for mathematical simplicity.

4. Case of Constant

We found it reasonable to start our analysis from the models with constant . Without loss of generality, we would like to describe equations allowing us to find dynamics of the models. According to the assumption with constant (16) will be modified as follows:

and, using the expression , we will obtain that

The last equation can be integrated very easily and the result is the following:

where is the scale factor and is the integration constant. In our future calculations we will use as initial condition. For the Hubble parameter we will obtain

Concerning the form of the field equations, we need only to assume the form of and we will obtain the cosmological solutions. Concerning the mathematical hardness of the problem we will analyze models numerically and investigate graphical behavior of various important cosmological parameters with respect to the cosmic time. In the following subsections we consider three models with the particular form of with given forms of the EoS for the viscous dark matter fluids considered in the Introduction.

4.1. Model 1

The first toy model describes the dynamics of the universe within an effective fluid in case of the cosmological constant. The dynamics of the energy density of the viscous modified Chaplygin gas which will model dark matter in our universe and the differential equation describing the dynamics of it can be found to be

From (23) for the dynamics of the dark energy, we have

The best fit for the theoretical model of our consideration with the observational data we obtained for , , , , , , , and .

In Figures 1–3, we have chosen to plot some of the quantities derived for different values of the parameters involved. In particular, we have chosen , , , , and and we have chosen five different values of , that is, 0, 0.2, 0.3, 0.5, and 0.7.

(a)

(b)

Figure 1: Behavior of Hubble parameter and against for the constant . Model 1.

(a)

(b)

Figure 2: Behavior of EoS parameters and against for the constant . Model 1.

Figure 3: Squared sound speed against for the constant . Model 1.

Figure 1 shows that the Hubble parameter is decreasing with time to a constant at the late universe as expected and its value increased by . On the other hand, the value of the deceleration parameter decreased with . It is illustrated by the right plot of Figure 1. For we can see , while for we can see . Also acceleration to deceleration phase transition is seen in this model. At the late time we have being in agreement with observational data. As we know recent observations of type SNIa indicate that universe is accelerating with the deceleration parameter lying somewhere in the range .

Moreover, Figure 2 shows that the EoS parameter tends to −1 at the late time with , corresponding to a quintessence-like universe.

Unfortunately, analysis of squared sound speed (see Figure 3) shows that this model is not stable at the late time and will be considered only for the early universe.

4.2. Model 2

In the second model, after some mathematical calculations, we obtain the following differential equation to study dynamics of the model:

where is given as

The best fitted values of parameters are , , , , , , , and .

In Figures 4 and 5, we have chosen to plot some of the quantities derived for different values of the parameters involved. In particular, we have chosen , , , , and and we have chosen five different values of , that is, 0, 0.5, 1.0, 1.2, and 1.5.

(a)

(b)

Figure 4: Behavior of Hubble parameter and against for the constant . Model 2.

Figure 5: Squared sound speed against for the constant . Model 2.

Numerical results of the Hubble expansion parameter and the deceleration parameter show good behavior (see Figure 4), but the stability analysis illustrated in Figure 5 shows that this model also has instability in the late time and is only useful for the early universe. However, the effect of constant in this parameter is similar to the first model and the EoS parameter tends to −1 as before. In the case of , we can see that , which is not coincident with observational data. It tells that presence of may be necessary to obtain agreement with observations.

Figure 5 shows that instability of model may solve for the large value of the cosmological constant. Therefore model 2 is stable at the late time for the large value of the .

4.3. Model 3

In the third model we have the following expression for the pressure :

The constant assumption will lead us to the following expression for the Hubble parameter :

The differential equations describing the dynamics of the energy densities of both components are given by the following equations:

where is the EoS parameter of dark energy. The best fit for the theoretical model of our consideration with the observational data is obtained for , , , , , , , and .

The behaviors of the cosmological parameters are similar to the previous models and we can see late time instability of this model in Figure 6. This suggests considering varying to obtain more appropriate models.

Figure 6: Squared sound speed against for the constant . Model 3.

5. The Case of Varying

In this section we will consider three interacting fluid models and will investigate cosmological parameters like the Hubble parameter , deceleration parameter , and EoS parameters of the total fluid and dark energy . Based on numerical solutions, we will discuss graphical behaviors of the cosmological parameters. For the varying we take a phenomenological form which was considered by us recently [59]. The formula of is given as the following expression:

which is a function of the Hubble parameter, potential of the scalar field, and time derivative of the scalar field. For the potential we take a simple form ; therefore the form of can be written also in the following way as only a function of the filed :

Therefore, the dynamics of can be obtained from the following differential equation:

In forthcoming subsections within three different forms of we will investigate the dynamics of the universe. The question of the dynamics for the energy densities of the dark energy and dark matter is already discussed in a previous section; therefore we will not consider them here and we will start with the comments on the graphical behaviors of the cosmological parameters of the models. We will start with the model where

5.1. Model 4

Interacting viscous modified Chaplygin gas with the quintessence dark energy in the case of varying gives the following differential equation:

with the Hubble parameter obtained as

In Figure 7, we can see that the Hubble expansion parameter is a decreasing function of time , which yields a constant at the late time as expected. It is clear that and increase the value of the Hubble expansion parameter but decrease the value of the deceleration parameter. In order to obtain the deceleration parameter being in agreement with observational data, we should choose larger values of and . Also acceleration to deceleration phase transition happens in this model. We find an instability at the initial time but the model is completely stable at the late time.

(a)

(b)

(c)

(d)

Figure 7: Behavior of Hubble parameter and deceleration parameter against for varying . Model 4.

Figure 8 shows that the EoS parameters yield −1 at the late time, in agreement with observational data. Also effects of and are illustrated in the plots of Figure 8.

(a)

(b)

Figure 8: Behavior of the EoS parameters and against for the case corresponding to varying . Model 4.

In Figure 9, we can see behavior of the Hubble expansion parameter with the redshift which is also in agreement with observational data, since it is increasing function.

Figure 9: Hubble parameter against redshift for varying . Model 4.

5.2. Model 5

A polytropic fluid interacting with the quintessence dark energy yields the following differential equation:

with the Hubble parameter obtained as

In Figure 10, we can see behavior of the Hubble expansion parameter and the deceleration parameter with time and find that the value of decreases the value of the Hubble expansion parameter while it increases the value of the deceleration parameter. Also, acceleration to deceleration phase transition is seen in this model.

(a)

(b)

Figure 10: Behavior of Hubble parameter and deceleration parameter against for varying . Model 5.

However, squared sound speed plotted in Figure 11 shows that this model, similar to the constant , is completely instable at the late time and will be useful only for the early universe.

Figure 11: Squared sound speed against for varying . Model 5.

5.3. Model 6

Finally, in the third model for the case of varying where a viscous Van der Waals gas with EoS of the general form

interacts with the quintessence dark energy, the differential equations describing the dynamics of the energy densities of both components are given:

where is the EoS parameter of the dark energy. We obtain the behavior of important cosmological parameters as illustrated in Figure 12. Moreover, Figure 13 shows that this model is also instable at the late time.

(a)

(b)

Figure 12: Behavior of Hubble parameter and deceleration parameter against for varying . Model 6.

Figure 13: Squared sound speed against for varying . Model 6.

Therefore we can choose the first model (model 4) as an appropriate model for the late time cosmology. Apart from the instabilities discussed above for the second and third models we can construct models by using more observational data discussed in the next section.

6. Observational Constraints on Interacting Models with Varying

The SNIa test is based on the distance modulus which is related to the luminosity distance by

where is defined as

The quantities and denote the apparent and the absolute magnitudes, respectively. Baryonic acoustic oscillations have their origin in oscillations in the photon-baryon plasma at the moment of the decoupling at about . They can be characterized by the distance scale

The WiggleZ data [58] indicates the following information about and : , , and for the redshifts , , and , respectively. The key quantity of a statistical analysis is the parameter:

where is the theoretical evaluation of a given observable, depending on free parameters, is the corresponding observational value, and is the total number of observational data for the given test. There are many different data sets, obtained with different techniques. In some cases, these different samples may give very different results. The second point is the existence of two different calibration methods. One of them uses cosmological relations and takes into account with high , and the other one, using astrophysical methods, is suitable for small (MLCS2k2). Our observational analysis of the background dynamics uses the following three tests: the differential age of old objects based on the dependence as well as the data from and from . A fourth test could potentially be added: the position of the first peak of the anisotropy spectrum of the cosmic microwave background radiation (CMB). However, the CMB test implies integration of the background equations until which requires the introduction of the radiative component. But the inclusion of such radiative component considerably changes the structure of the equations and no analytic expression for is available. Hence, we will limit ourselves to the mentioned three tests for which a reliable estimation is possible.

In Tables 1 and 2 we fix parameters of three models by using mentioned observational data.

Table 1

4

5

6

Table 2

4

—
—

5
—
—
—

6

—
—
—

In Figure 14, we can see the behavior of in all models, which is approximately in agreement with observational data.

(a)

(b)

(c)

Figure 14: Observational data (SNeIa+BAO+CMB) for distance modulus versus our theoretical results for varying in models 4, 5, and 6.

7. Conclusion

In this paper we considered three different cosmological models for the universe based on Lyra geometry. First of all we introduced our models and then obtained field equations solved numerically to find behavior of cosmological parameters. In order to find effect of varying we also studied the case of constant and found that presence of is necessary to obtain results being in agreement with observational data. We assumed viscous modified Chaplygin gas (models 1 and 4), viscous polytropic gas (models 2 and 5), or viscous Van der Waals gas (models 3 and 6) to be a component which at the early universe plays the role of dark matter with , but at late times it tends to a cosmological constant. Moreover, we have a quintessence field, which will contribute to the dark energy sector including possibility of interaction between components. Easily we can check that and are of the same order. Also we considered case of varying and studied behavior of cosmological parameters numerically. We used observational data to fix parameters of the models and seen agreement with observational data by investigation of . By using stability analysis we concluded that the model is the best model considered in this paper to describe universe.

For the future work we can extend the present paper to include shear viscosity or varying bulk viscosity [53], and we can also consider cosmic Chaplygin gas versions [54] to obtain more general model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

A. G. Riess, A. V. Filippenko, P. Challis, et al., “Observational evidence from supernovae for an accelerating universe and a cosmological constant,” The Astronomical Journal, vol. 116, no. 3, p. 1009, 1998.

S. Perlmutter, G. Aldering1, G. Goldhaber, et al., “Measurements of and from 42 high-redshift supernovae,” The Astrophysical Journal, vol. 517, p. 565, 1999.

R. Amanullah, C. Lidman, D. Rubin, et al., “Spectra and hubble space telescope light curves of six type Ia supernovae at and the union2 compilation,” The Astrophysical Journal, vol. 716, p. 712, 2010.

A. C. Pope, T. Matsubara, A. S. Szalay, et al., “Cosmological parameters from eigenmode analysis of sloan digital sky survey galaxy redshifts,” The Astrophysical Journal, vol. 607, no. 2, pp. 655–660, 2004.

D. N. Spergel, L. Verde, H. V. Peiris, et al., “First year Wilkinson Microwave Anisotropy Probe (WMAP) observations: determination of cosmological parameters,” The Astrophysical Journal Supplement, vol. 148, p. 175, 2003.

M. Jamil, K. Yesmakhanova, D. Momeni, and R. Myrzakulov, “Phase space analysis of interacting dark energy in f(T) cosmology,” Central European Journal of Physics, vol. 10, no. 5, pp. 1065–1071, 2012.

W. Zimdahl, “Models of interacting dark energy,” AIP Conference Proceedings, vol. 147, p. 51, 2012.

M. Khurshudyan, B. Pourhassan, and E. Chubaryan, “A Universe with a generalized ghost dark energy and Van der Waals fluid interacting with a fluid,” http://arxiv.org/abs/1302.1819

Y. L. Bolotin, A. Kostenko, O. A. Lemets, and D. A. Yerokhin, “Cosmological evolution with interaction between dark energy and dark matter,” http://arxiv.org/abs/1310.0085

M. Khurshudyan, B. Pourhassan, and E. O. Kahya, “Interacting two-component fluid models with varying EoS parameter,” International Journal of Geometric Methods in Modern Physics, vol. 11, no. 6, Article ID 1450061, 14 pages, 2014.

P. P. Avelino and H. M. R. da Silva, “Effective dark energy equation of state in interacting dark energy models,” Physics Letters B, vol. 714, no. 1, pp. 6–10, 2012.

J. Sadeghi, M. Khurshudyan, A. Movsisyan, and H. Farahani, “Interacting ghost dark energy models in the higher dimensional cosmology,” http://arxiv.org/abs/1401.6649

M. Khurshudyan, J. Sadeghi, A. Pasqua, S. Chattopadhyay, R. Myrzakulov, and H. Farahani, “Interacting Ricci dark energy models with an effective -term in Lyra manifold,” http://arxiv.org/abs/1402.5678

M. Khurshudyan, J. Sadeghi, M. Hakobyan, H. Farahani, and R. Myrzakulov, “Interaction between modified Chaplygin gas and ghost dark energy in the presence of extra dimensions,” The European Physical Journal Plus, vol. 129, article 119, 2014.

J. Sola and H. Stefancic, “Effective equation of state for dark energy: Mimicking quintessence and phantom energy through a variable Λ,” Physics Letters B, vol. 624, no. 3-4, pp. 147–157, 2005.

I. L. Shapiro and J. Solà, “On the possible running of the cosmological ‘constant’,” Physics Letters B, vol. 682, no. 1, pp. 105–113, 2009.

B. Ratra and P. J. E. Peebles, “Cosmological consequences of a rolling homogeneous scalar field,” Physical Review D, vol. 37, no. 12, pp. 3406–3427, 1988.

M. Khurshudyan, E. Chubaryan, and B. Pourhassan, “Interacting quintessence models of dark energy,” International Journal of Theoretical Physics, vol. 53, no. 7, pp. 2370–2378, 2014.

Z. K. Guo, N. Ohta, and Y. Z. Zhang, “Parametrizations of the dark energy density and scalar potentials,” Modern Physics Letters A, vol. 22, no. 12, p. 883, 2007.

S. Dutta, E. N. Saridakis, and R. J. Scherrer, “Dark energy from a quintessence (phantom) field rolling near a potential minimum (maximum),” Physical Review D, vol. 79, Article ID 103005, 2009.

E. N. Saridakis and S. V. Sushkov, “Quintessence and phantom cosmology with nonminimal derivative coupling,” Physical Review D, vol. 81, Article ID 083510, 2010.

R. R. Caldwell, M. Kamionkowski, and N. N. Weinberg, “Phantom energy: dark energy with causes a cosmic doomsday,” Physical Review Letters, vol. 91, Article ID 071301, 2003.

B. Feng, X. L. Wang, and X. M. Zhang, “Dark energy constraints from the cosmic age and supernova,” Physics Letters B, vol. 607, no. 1-2, pp. 35–41, 2005.

E. Elizalde, S. Nojiri, and S. D. Odintsov, “Late-time cosmology in a (phantom) scalar-tensor theory: dark energy and the cosmic speed-up,” Physical Review D, vol. 70, Article ID 043539, 2004.

Z. K. Guo, “Cosmological evolution of a quintom model of dark energy,” Physics Letters B, vol. 608, pp. 177–182, 2005.

M.-Z. Li, B. Feng, and X.-M. Zhang, “A single scalar field model of dark energy with equation of state crossing −1,” Journal of Cosmology and Astroparticle Physics, vol. 2005, p. 002, 2005.

B. Feng, M. Li, Y.-S. Piao, and X. Zhang, “Oscillating quintom and the recurrent universe,” Physics Letters B, vol. 634, no. 2-3, pp. 101–105, 2006.

S. Capozziello, S. Nojiri, and S. D. Odintsov, “Unified phantom cosmology: inflation, dark energy and dark matter under the same standard,” Physics Letters, Section B, vol. 632, no. 5-6, pp. 597–604, 2006.

W. Zhao and Y. Zhang, “Quintom models with an equation of state crossing −1,” Physical Review D, vol. 73, Article ID 123509, 2006.

Y. F. Cai, T. Qiu, Y. S. Piao, M. Li, and X. Zhang, “Bouncing universe with quintom matter,” Journal of High Energy Physics, vol. 710, article 71, 7 pages, 2007.

E. N. Saridakis and J. M. Weller, “A Quintom scenario with mixed kinetic terms,” Physical Review D, vol. 81, no. 12, Article ID 123523, 11 pages, 2010.

Y. F. Cai, T. Qiu, R. Brandenberger, Y. S. Piao, and X. Zhang, “On perturbations of quintom bounce,” Joint Center for Artificial Photosynthesis, vol. 2008, no. 03, article 013, 2008.

M. R. Setare and E. N. Saridakis, “Coupled oscillators as models of quintom dark energy,” Physics Letters B, vol. 668, no. 3, pp. 177–181, 2008.

M. R. Setare and E. N. Saridakis, “Quintom cosmology with general potentials,” International Journal of Modern Physics D, vol. 18, no. 4, pp. 549–557, 2009.

Y. F. Cai, E. N. Saridakis, M. R. Setare, and J. Xia, “Quintom cosmology: theoretical implications and observations,” Physics Reports, vol. 493, no. 1, pp. 1–60, 2010.

T. Qiu, “Theoretical aspects of Quintom models,” Modern Physics Letters A, vol. 25, no. 11-12, pp. 909–921, 2010.

Q. G. Huang and M. Li, “The holographic dark energy in a non-flat universe,” Journal of Cosmology and Astroparticle Physics, vol. 8, 2004.

M. Ito, “Holographic-dark-energy model with non-minimal coupling,” Europhysics Letters, vol. 71, no. 5, pp. 712–716, 2005.

X. Zhang and F. Q. Wu, “Constraints on holographic dark energy from type Ia supernova observations,” Physical Review D, vol. 72, Article ID 043524, 2005.

S. Nojiri and S. D. Odintsov, “Unifying phantom inflation with late-time acceleration: scalar phantom—non-phantom transition model and generalized holographic dark energy,” General Relativity and Gravitation, vol. 38, pp. 1285–1304, 2006.

E. Elizalde, S. Nojiri, S. D. Odintsov, and P. Wang, “Dark energy: vacuum fluctuations, the effective phantom phase, and holography,” Physical Review D, vol. 71, Article ID 103504, 2005.

H. Li, Z. K. Guo, and Y. Z. Zhang, “A tracker solution for a holographic dark energy model,” International Journal of Modern Physics D, vol. 15, p. 869, 2006.

E. N. Saridakis, “Restoring holographic dark energy in brane cosmology,” Physics Letters B, vol. 660, no. 3, pp. 138–143, 2008.

E. N. Saridakis, “Holographic dark energy in braneworld models with moving branes and the crossing,” Journal of Cosmology and Astroparticle Physics, vol. 2008, no. 4, article 020, 2008.

E. N. Saridakis, “Holographic dark energy in braneworld models with a Gauss-Bonnet term in the bulk. Interacting behavior and the crossing,” Physics Letters B, vol. 661, no. 5, pp. 335–341, 2008.

R.-G. Cai, “A dark energy model characterized by the age of the universe,” Physics Letters B, vol. 657, no. 4-5, pp. 228–231, 2007.

H. Wei and R. G. Cai, “Interacting agegraphic dark energy,” The European Physical Journal C, vol. 59, pp. 99–105, 2009.

O. F. Piattella, J. C. Fabris, and W. Zimdahl, “Bulk viscous cosmology with causal transport theory,” Journal of Cosmology and Astroparticle Physics, vol. 2011, no. 5, article 029, 2011.

K. Karami, Z. Safari, and S. Asadzadeh, “Cosmological constraints on polytropic gas model,” http://arxiv.org/abs/1209.6374

K. Karami, S. Ghaffari, and J. Fehri, “Interacting polytropic gas model of phantom dark energy in non-flat universe,” European Physical Journal C, vol. 64, no. 1, pp. 85–88, 2009.

K. Karami and S. Ghaffari, “The generalized second law of thermodynamics for the interacting polytropic dark energy in non-flat FRW universe enclosed by the apparent horizon,” Physics Letters B, vol. 688, p. 125, 2010.

G. M. Kremer, “Cosmological models described by a mixture of van der Waals fluid and dark energy,” Physical Review D, vol. 68, no. 12, Article ID 123507, 2003.

H. Saadat and B. Pourhassan, “Effect of varying bulk viscosity on generalized chaplygin gas,” International Journal of Theoretical Physics, vol. 53, no. 4, pp. 1168–1173, 2014.

B. Pourhassan, “Viscous modified cosmic chaplygin gas cosmology,” International Journal of Modern Physics D, vol. 22, no. 9, Article ID 1350061, 2013.

H. Saadat and B. Pourhassan, “FRW bulk viscous cosmology with modified Chaplygin gas in flat space,” Astrophysics and Space Science, vol. 343, no. 2, pp. 783–786, 2013.

H. Saadat and B. Pourhassan, “FRW bulk viscous cosmology with modified cosmic Chaplygin gas,” Astrophysics and Space Science, vol. 344, no. 1, pp. 237–241, 2013.

V. K. Shchigolev, “Cosmology with an effective Λ-term in Lyra Manifold,” Chinese Physics Letters, vol. 30, Article ID 119801, 2013.

M. Jamil, S. Ali, D. Momeni, and R. Myrzakulov, “Bianchi type I cosmology in generalized Saez-Ballester theory via Noether gauge symmetry,” European Physical Journal C, vol. 72, no. 4, article 1998, pp. 1–6, 2012.

M. Khurshudyan, “Toy models of Universe with an effective varying -term,” http://arxiv.org/abs/1403.0109