Astrophysical observations have put unprecedentedly tight constraints on cosmological theories. The Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM model, mathematically simple and fits observational data sets well, is preferred for explaining the behavior of universe. But many basic features of the dark sectors are still unknown, which leaves room for various nonstandard cosmological hypotheses. As the pressure of the cosmological constant dark energy is unvarying, ignoring contributions from radiation and curvature terms at low redshift, the effective pressure keeps constant. In this paper, we propose two parametric models for a non-constant effective pressure in order to study the tiny deviation from Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM at low redshift. We recover our phenomenological models in the scenarios of quintessence and phantom fields, and we explore the behavior of the scalar field and potential. We constrain our model parameters with SNe Ia and BAO observations, and we detect subtle hints of ωde<-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}<-1$$\end{document} from the data-fitting results of both models, which indicates possibly a phantom dark energy scenario at present.
volume-issue-count12issue-article-count55issue-toc-levels0issue-typeRegularissue-online-date-year2015issue-online-date-month8issue-online-date-day27issue-pricelist-year2015issue-copyright-holderSIF and Springer-Verlag Berlin Heidelbergissue-copyright-year2015article-contains-esmNoarticle-numbering-styleContentOnlyarticle-toc-levels0article-registration-date-year2015article-registration-date-month6article-registration-date-day19article-grants-typeOpenChoicemetadata-grantOpenAccessabstract-grantOpenAccessbodypdf-grantOpenAccessbodyhtml-grantOpenAccessbibliography-grantOpenAccessesm-grantOpenAccessIntroduction
Since the discovery of the current acceleration of our universe expansion in 1998, maybe the greatest mystery in cosmology is the deceptive nature of the dark energy. Recent observational results [1] have put tight constraints on the properties of dark energy, but there is still no theoretical or observational indication pinning down its nature. On the one hand, although the simple cosmological constant Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document} can accommodate the accelerating expansion, it encounters two serious problems. The first one is the fine tuning problem: the measured energy of the vacuum is so much smaller than the estimated value ρvacobs≪ρvactheo\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho _{\mathrm{vac}}^{\mathrm{obs}}\ll \rho _{\mathrm{vac}}^{\mathrm{theo}}$$\end{document}, the famous 120-orders-of-magnitude discrepancy that makes the vacuum explanation suspect. On the other hand we may ask why there is dominance of the cosmological constant over the matter component at the present epoch. These two basic problems prompt us to propose some alternatives, which include an evolving scalar field called quintessence [2–8], a noncanonical scalar field (such as K-essence [9–11], phantom [7, 8, 12–18]), modified gravity [7, 8, 19–23], coupled dark energy [8, 24, 25] or decaying dark energy [26] models, and so on. On the other hand, the equation of state (EoS) parameter of the cosmological constant is precisely ωde=-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=-1$$\end{document}. Recent observations show that the EoS parameter of modeled dark energy is ωde=-1.006±0.045\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=-1.006\pm 0.045$$\end{document}, which slightly favors ωde<-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}<-1$$\end{document}. Anyhow, the small deviations from the cosmological constant Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document} allow one to consider models with ωde≠-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}\ne -1$$\end{document}. So one can make efforts to construct new models to explain the deviations which may be detectable at the precision of current and future observations.
Parameterization is an useful tool toward a more complete characterization of dark energy modeling and has been routinely employed to analyze data sets. Most parameterizations for dark energy models involve the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} for the dark energy behavior. Several well-known parameterizations for the EoS of dark energy have been proposed so far. We can write the parameterizations in polynomial form ωde(z)=∑n=0ωnxn(z)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}(z)=\sum \nolimits _{n=0}\omega _nx_n(z)$$\end{document} generally, where the expansions can be given in the following ways: (i) by redshift xn(z)=zn\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$x_n(z)=z^n$$\end{document}, (ii) by scale factor xn(z)=(1-aa0)n=(z1+z)n\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$x_n(z)=(1-\frac{a}{a_0})^n=(\frac{z}{1+z})^n$$\end{document}, (iii) logarithmic xn(z)=[ln(1+z)]n\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$x_n(z)=[\ln (1+z)]^n$$\end{document}. Parameterization (i) was proposed by Huterer and Turner [27] and Weller and Albrecht [28] with n≤1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$n\le 1$$\end{document}. Parameterization (ii) with n≤1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$n\le 1$$\end{document} was introduced by Chevalier, Polarski and Linder [29, 30], the famous Chevallier–Polarski–Linder (CPL) parameterization. ωde=ω0+ω1(1-a)=ω0+ω1z1+z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=\omega _0+\omega _1(1-a)=\omega _0+\omega _1\frac{z}{1+z}$$\end{document} behaves as ωde→ω0+ω1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}\rightarrow \omega _0+\omega _1$$\end{document} for z→∞\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z\rightarrow \infty $$\end{document} and ωde→ω0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}\rightarrow \omega _0$$\end{document} for z→0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z\rightarrow 0$$\end{document}. A more general form with ωde=ω0+ω1z(1+z)p\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=\omega _0+\omega _1\frac{z}{(1+z)^p}$$\end{document} was later proposed by Jassal et al. [32]. Parameterization (iii) with n≤1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$n\le 1$$\end{document} was introduced by Efstathiou [31]. In recent years, some new parameterizations have been proposed, such as using Padé parameterizations for the EoS of dark energy [33], namely ωde=ω0+ωa(1-a)1+ωb(1-a)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=\frac{\omega _0+\omega _a(1-a)}{1+\omega _b(1-a)}$$\end{document}, and ωde=ω0+ω1lna1+ω2lna\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=\frac{\omega _0+\omega _1\ln a}{1+\omega _2\ln a}$$\end{document}. It is worth mentioning that Sen proposed a parameterization for the pressure of the dark energy model [34, 35], PΛ=-P0+P1(1-a)+····\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_\Lambda =-P_0+P_1(1-a)+\cdot \!\cdot \!\cdot \cdot $$\end{document}, in order to study small deviations from the cosmological constant. Different from parameterizations which focused on the EoS of the dark energy mentioned above, in this paper we aim to make parameterizations for the relation between redshift and effective pressure of all energy components in the universe. In the following we propose two parametric models for the effective pressure in order to explore the late-stage evolution of the universe.
This paper is organized as follows: in Sect. 2, we propose two new parametric models for the effective pressure: P(z)=Pa+Pbz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P(z)=P_a+P_b z$$\end{document} and P(z)=Pc+Pd1+z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P(z)=P_c+\frac{P_d}{1+z}$$\end{document}. In Sect. 3, we relate our parametric models with the quintessence and phantom scalar fields, and the behavior of field and potential is then explored. In Sect. 4, we constrain our model parameters with SNe Ia and BAO observations. In Sect. 5, we end with discussions and conclusions.
Two parametric models
The Friedmann equations, the equation of energy conservation, and the equation of state constitute a closed system describing the background evolution of the universe. Substituting the EoS by a relation between the effective pressure P and the redshift z is also feasible, as the equation P=P(z)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P=P(z)$$\end{document} is not linearly dependent on the Friedmann equation and the equation of energy conservation. Also, the EoS can be recovered by inserting the P–z relation into the equation of energy conservation,ρ˙+3H(P+ρ)=0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \dot{\mathrm{\rho }}+3 H(P+\rho )=0, \end{aligned}$$\end{document}and integrating out the expression of ρ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho $$\end{document}. For example, the effective pressure for Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM at late stage is nearly constant, say P0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_0 $$\end{document}; accordingly, we can obtain from Eq. (1)ρ(a)=-P0+Ca-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho (a)=-P_0+C a^{-3} \end{aligned}$$\end{document}where C is an integration constant, and the two terms at the right side represent contributions from the cosmological constant and matter, respectively.
This is just an example of P parameterization; generally, we can have more complicated P–z relations. As the P–z relation is equivalent to the EoS, a parameterization on the effective pressure is equivalent to that of the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document}. Since ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} is the exponential of some component in EoS, the ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} parameterization requires a presupposition of the components in EoS; i.e., the physical mechanism of the possible deviation from Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM has to be dictated; we make parameterizations merely because we actually do not know the concrete mechanism behind the accelerative expansion. To illustrate, a deviation of Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM might come from the evolution of the EoS of the cosmological constant term, while an additional component might result in the same deviation. However, a parameterization of the effective pressure just circumvents this issue, and no knowledge of the concrete physical mechanism is required. We are able to directly study the deviation from the constant P–z relation without prejudice to a presupposition.
Model 1
In this subsection, we propose a model which readsP(z)=Pa+Pbz,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P(z)=P_a+P_b z, \end{aligned}$$\end{document}where Pa\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_a$$\end{document} and Pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_b$$\end{document} are free parameters.
For the scale factor a and the redshift z, we havea=a01+z=11+z,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} a=\frac{a_0}{1+z}=\frac{1}{1+z}, \end{aligned}$$\end{document}where a0=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a_0=1$$\end{document} corresponds to the value today. Substitute Eqs. (3) and (4) into Eq. (1); then the total energy density can be integrated asρ(a)=-(Pa-Pb)-32Pba-1+C1a-3,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho (a)=-(P_a-P_b)-\frac{3}{2}P_b a^{-1}+C_1 a^{-3}, \end{aligned}$$\end{document}where C1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_1$$\end{document} is an integration constant. If we set ρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho _0$$\end{document} to be the energy density today, the integration constant is then C1=ρ0+Pa+12Pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_1=\rho _0+P_a+\frac{1}{2}P_b$$\end{document}. In Eq. (5), we can interpret the inversely cubic term C1a-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_1 a^{-3}$$\end{document} as dust matter and the constant term -(Pa-Pb)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-(P_a-P_b)$$\end{document} as the cosmological constant in Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM. The term -32Pba-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-\frac{3}{2}P_b a^{-1}$$\end{document} does not appear in the Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM model, whose physical nature will be explored in the next section.
For convenience in data fitting, we introduce some dimensionless parameters. First, we define the dimensionless density and pressure asρ∗≡ρρ0=H2H02,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho ^*\equiv & {} \frac{\rho }{\rho _0}=\frac{H^2}{H_0^2} , \end{aligned}$$\end{document}P∗≡Pρ0.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P^*\equiv & {} \frac{P}{\rho _0} . \end{aligned}$$\end{document}The expressions of the total density Eq. (5) and the total pressure Eq. (3) can be rewritten asρ∗(a)=-(Pa∗-Pb∗)-32Pb∗a-1+C1∗a-3,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho ^*(a)= & {} -(P_a^*-P_b^*)-\frac{3}{2}P_b^*a^{-1}+C_1^*a^{-3} ,\end{aligned}$$\end{document}P∗(a)=(Pa∗-Pb∗)+Pb∗a-1,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P^*(a)= & {} (P_a^*-P_b^*)+P_b^*a^{-1} , \end{aligned}$$\end{document}where Pa∗≡Paρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_a^*\equiv \frac{P_a}{\rho _0}$$\end{document}, Pb∗≡Pbρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_b^*\equiv \frac{P_b}{\rho _0}$$\end{document}, and C1∗≡C1ρ0=1+Pa∗+12Pb∗\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C^*_1\equiv \frac{C_1}{\rho _0}=1+P_a^*+\frac{1}{2}P_b^*$$\end{document}.
We propose another parameterization asP(z)=Pc+Pd1+z,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P(z)=P_c+\frac{P_d}{1+z}, \end{aligned}$$\end{document}where Pc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_c$$\end{document} and Pd\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_d$$\end{document} are free parameters. Inserting Eqs. (4) and (17) into Eq. (1), we obtain the total energy density for model 2,ρ(a)=-Pc-34Pda+C2a-3,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho (a)=-P_c-\frac{3}{4}P_d a+C_2 a^{-3}, \end{aligned}$$\end{document}where C2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_2$$\end{document} is an integration constant. Setting the present energy density as ρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho _0$$\end{document}, then C2=ρ0+Pc+34Pd\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_2=\rho _0+P_c+\frac{3}{4}P_d$$\end{document}. Still, we can find the term C2a-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_2 a^{-3}$$\end{document} corresponding to dust matter, and the term -Pc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-P_c$$\end{document} corresponding to the cosmological constant. The difference between model 2 and model 1 rests in the rest term, -32Pba-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-\frac{3}{2}P_b a^{-1}$$\end{document} for model 2, whereas it is in -34Pda\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-\frac{3}{4}P_d a$$\end{document} for model 1. Their physical nature will be studied in the next section.
Like model 1, we need to introduce new model parameters in model 2. With Eqs. (6) and (7), we can obtain the expressions of the total density and the total pressure for model 2:ρ∗(a)=-Pc∗-34Pd∗a+C2∗a-3,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho ^*(a)= & {} -P_c^*-\frac{3}{4}P_d^*a+C_2^*a^{-3} , \end{aligned}$$\end{document}P∗(a)=Pc∗+Pd∗a,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P^*(a)= & {} P_c^*+P_d^*a , \end{aligned}$$\end{document}where Pc∗≡Pcρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_c^*\equiv \frac{P_c}{\rho _0}$$\end{document}, Pd∗≡Pdρ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_d^*\equiv \frac{P_d}{\rho _0}$$\end{document}, and C2∗≡C2ρ0=1+Pc∗+34Pd∗\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C^*_2\equiv \frac{C_2}{\rho _0}=1+P_c^*+\frac{3}{4}P_d^*$$\end{document}.
Deviations from the Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM in our models can be realized through different physical scenarios. Scalar fields are mainstream approaches to explain the acceleration of the universe’s expansion. In the scenarios of scalar fields, dark energy evolves with time. The scalar field dynamics has been studied in great detail (see Refs. [2–18]) and there are lots of issues involved such as (i) choosing the initial conditions for the scalar field; (ii) choosing the potential with solid theoretical motivation; (iii) the existence of the tracker field and so on. Generally, the evolution of a scalar field is studied over the cosmic history, and once the parameters of scalar field models are set they determine the entire cosmological evolution. So a more detailed analysis would involve studying the scalar field dynamics over cosmic history, and then comparing its evolution with that of a pressure parametrization model at low redshift. In this paper, we will merely compare the pressure and the energy density of a field with that of a model of pressure parametrization at low redshift and study the behavior of the field and the potential. The physical realization of the parameterizations through scalar fields means adjusting the behavior of the scalar fields to the dark energy term occurring in the parametric model. Specifically, we write the two equationsPeff=Pscalarfield,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&P_{\mathrm{eff}} = P_{\mathrm{scalar~field}},\end{aligned}$$\end{document}ρeff-ρm=ρscalarfield,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\rho _{\mathrm{eff}}-\rho _m = \rho _{\mathrm{scalar~field}}, \end{aligned}$$\end{document}as the mathematical definition of the realization.
In this section, we will take “quintessence” and “phantom” as two examples to realize our models.
Quintessence “Quintessence” denotes a canonical scalar field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} with a potential V1(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )$$\end{document} that does not interact with all the other components except standard gravity, whose EoS parameter ωde>-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}>-1$$\end{document}. Quintessence is described by the actionwhere κ2=8πG\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\kappa ^2=8\pi G$$\end{document}, R is the Ricci scalar, and SM\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$S_M$$\end{document} is the action of matter. The variation of the action Eq. (29) with respect to ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} givesϕ¨+3Hϕ˙+V1′(ϕ)=0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \ddot{\phi } + 3 H \dot{\phi } + V^\prime _1(\phi ) = 0, \end{aligned}$$\end{document}where V1(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )$$\end{document} is the potential of the quintessence field, the prime denotes the derivative with respect to ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. In a FLRW background, the energy density ρde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho _{\mathrm{de}}$$\end{document} and the pressure Pde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_{\mathrm{de}}$$\end{document} of the quintessence field areρde=12ϕ˙2+V1(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho _{\mathrm{de}}= & {} \frac{1}{2}\dot{\phi }^2+V_1(\phi ), \end{aligned}$$\end{document}Pde=12ϕ˙2-V1(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P_{\mathrm{de}}= & {} \frac{1}{2}\dot{\phi }^2-V_1(\phi ). \end{aligned}$$\end{document}Then the EoSωde=12ϕ˙2-V1(ϕ)12ϕ˙2+V1(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \omega _{\mathrm{de}}=\frac{\frac{1}{2}\dot{\phi }^2-V_1(\phi )}{\frac{1}{2}\dot{\phi }^2+V_1(\phi )}. \end{aligned}$$\end{document}Phantom The minimally coupled phantom model is also a possible realization, whose EoS parameter ωde<-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}<-1$$\end{document}. The action of the phantom field minimally coupled to gravity and matter sources iswhose variation with respect to ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} givesϕ¨+3Hϕ˙-V2′(ϕ)=0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \ddot{\phi } + 3 H \dot{\phi } - V^\prime _2(\phi ) = 0, \end{aligned}$$\end{document}where V2(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2(\phi )$$\end{document} is the potential of the phantom field, and the prime denotes the derivative with respect to ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. The energy density and pressure of the phantom are given by (assuming a flat FRW metric)ρde=-12ϕ˙2+V2(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \rho _{\mathrm{de}}= & {} -\frac{1}{2}\dot{\phi }^2+V_2(\phi ) ,\end{aligned}$$\end{document}Pde=-12ϕ˙2-V2(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} P_{\mathrm{de}}= & {} -\frac{1}{2}\dot{\phi }^2-V_2(\phi ) . \end{aligned}$$\end{document}The EoS of the phantom field is thenωde=--12ϕ˙2-V2(ϕ)-12ϕ˙2+V2(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \omega _{\mathrm{de}}=-\frac{-\frac{1}{2}\dot{\phi }^2-V_2(\phi )}{-\frac{1}{2}\dot{\phi }^2+V_2(\phi )}. \end{aligned}$$\end{document}So ωde<-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}<-1$$\end{document} for 12ϕ˙2<V2(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2<V_2(\phi )$$\end{document}.
Model 1
The EoS of the scalar fields for model 1 readsωde=Pscalarfieldρscalarfield=-1+13β(1+z)α+β(1+z).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \omega _{\mathrm{de}}=\frac{ P_{\mathrm{scalar~field}}}{\rho _{\mathrm{scalar~field}}}=-1+\frac{ \frac{1}{3}\beta (1+z)}{\alpha +\beta (1+z)}. \end{aligned}$$\end{document}Note that, in the above equation, there will be a singularity when z=-αβ-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z=-\frac{\alpha }{\beta }-1$$\end{document}. In this paper we only consider the universe at low redshift, so we do not need to worry about that situation. Besides, in Sect. 4 data fitting will support our argument.
In the quintessence scenario, assuming the cosmic components consist of matter and quintessence, comparing Eqs. (31) and (32) with Eqs. (3) and (5), we have-(Pa-Pb)-32Pba-1=12ϕ˙2+V1(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&-(P_a-P_b)-\frac{3}{2}P_b a^{-1} = \frac{1}{2}\dot{\phi }^2+V_1(\phi ) ,\end{aligned}$$\end{document}Pa-Pb+Pba-1=12ϕ˙2-V1(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&P_a-P_b+P_b a^{-1} = \frac{1}{2}\dot{\phi }^2-V_1(\phi ) . \end{aligned}$$\end{document}Simplify the above two equations, compare to Eqs. (6)–(11), replace model parameters (Pa\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_a$$\end{document}, Pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_b$$\end{document}) with the redefined parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}), and we obtain12ϕ˙2=16ρ0βa-1,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \frac{1}{2}\dot{\phi }^2= & {} \frac{1}{6}\rho _0\beta a^{-1} ,\end{aligned}$$\end{document}V1(ϕ)=ρ0α+56ρ0βa-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} V_1(\phi )= & {} \rho _0\alpha +\frac{5}{6}\rho _0\beta a^{-1} . \end{aligned}$$\end{document}From Eq. (43), it is easy to find that β>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta >0$$\end{document} in the scenario of quintessence. By Eqs. (43) and (44), one can construct the kinetic energy 12ϕ˙2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2$$\end{document} and the potential V1(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )$$\end{document} of the quintessence field with parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}) of model 1. In order to solve the above two equations, following [36], we choose the condition ϕa=1=MP\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi _{a=1}=M_P$$\end{document}, where MP\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M_P$$\end{document} is the reduced Planck mass. The Friedmann equation can then be rewritten asH2=13MP2ρ.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} H^2=\frac{1}{3M_P^2}\rho . \end{aligned}$$\end{document}Considering the dark energy domination at the present epoch in the universe, with the density parameter in the dark energy Ωde∼0.7\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{\mathrm{de}}\sim 0.7$$\end{document}, we define V0=ρ0=3MP2H02\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_0=\rho _0=3M_P^2H_0^2$$\end{document}. Simplifying Eqs. (43) and (44), we havedϕda=±MPβα+βa-1+(1-α-β)a-3a-32,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\frac{\mathrm{d}\phi }{\mathrm{d}a} = \pm M_P\sqrt{\frac{\beta }{\alpha +\beta a^{-1}+(1-\alpha -\beta ) a^{-3}}} a^{-\frac{3}{2}} , \end{aligned}$$\end{document}V1(ϕ)=V0α+56βa-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&V_1(\phi ) = V_0\left( \alpha +\frac{5}{6}\beta a^{-1}\right) . \end{aligned}$$\end{document}The symbol “±\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\pm $$\end{document}” in Eq. (46) corresponds to two solutions. Consider α=0.7,β=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7, \beta =0.05$$\end{document} for numerically solving the above two equations; the solutions are represented in Figs. 1 and 2, respectively. From Fig. 1, we find that ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} increases with a, the potential decreases with the increasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}, Eq. (47), implying that the potential will reach the minimum value V1(ϕ)=V0α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )=V_0\alpha $$\end{document} in the future. From Fig. 2, we can see that ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} decreases with a, and the potential decreases with decreasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}, Eq. (47) implies that the potential will reach the minimum V1(ϕ)=V0α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )=V_0\alpha $$\end{document} in the future. By Eqs. (14) and (15), we can obtain the expression of the density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} for model 1:Ωϕ=Ω1+Ω2=α+βa-1E2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Omega _\phi =\Omega _1+\Omega _2=\frac{\alpha +\beta a^{-1}}{E^2}. \end{aligned}$$\end{document}The evolution of the density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} in the scenario of quintessence is plotted in Fig. 3. From Fig. 3, we can see that until low redshift the energy density in the quintessence field becomes cosmologically dominant. Finally, the field comes to rest at the minimum of the potential V1(ϕ)=V0α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )=V_0\alpha $$\end{document}, and the universe eventually settles in a de Sitter phase [see Eq. (40)].
The solution of Eqs. (46) and (47) corresponding to a plus sign in Eq. (46). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values α=0.7,β=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =0.05$$\end{document}
The solution of Eqs. (46) and (47) corresponding to a minus sign in Eq. (46). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values α=0.7,β=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =0.05$$\end{document}
Evolution of the density parameters in the quintessence field (Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document}) and matter (Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document}) for model 1. Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} is indicated by solid line, and Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document} is indicated by a dashed line. We consider values α=0.7,β=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =0.05$$\end{document}
In the case of the phantom scenario, assuming the cosmological components to consist of matter and phantom, comparing Eqs. (37) and (38) with Eqs. (5) and (3), then we have-(Pa-Pb)-32Pba-1=-12ϕ˙2+V2(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&-(P_a-P_b)-\frac{3}{2}P_b a^{-1} = -\frac{1}{2}\dot{\phi }^2+V_2(\phi ) ,\end{aligned}$$\end{document}Pa-Pb+Pba-1=-12ϕ˙2-V2(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&P_a-P_b+P_b a^{-1} = -\frac{1}{2}\dot{\phi }^2-V_2(\phi ) . \end{aligned}$$\end{document}Replace the model parameters (Pa\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_a$$\end{document}, Pb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_b$$\end{document}) with the redefined parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}), we have12ϕ˙2=-16ρ0βa-1,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \frac{1}{2}\dot{\phi }^2= & {} -\frac{1}{6}\rho _0\beta a^{-1} ,\end{aligned}$$\end{document}V2(ϕ)=ρ0α+56ρ0βa-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} V_2(\phi )= & {} \rho _0\alpha +\frac{5}{6}\rho _0\beta a^{-1} . \end{aligned}$$\end{document}From Eq. (51), it is easy to find that in the scenario of phantom, β<0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta <0$$\end{document}. By Eqs. (51) and (52), one can construct the kinetic energy 12ϕ˙2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2$$\end{document} and potential V2(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2(\phi )$$\end{document} of the phantom field with the model parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}). Equations (51) and (52) can be rewritten asdϕda=±MP-βα+βa-1+(1-α-β)a-3a-32,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\frac{\mathrm{d}\phi }{\mathrm{d}a} = \pm M_P\sqrt{\frac{-\beta }{\alpha +\beta a^{-1}+(1-\alpha -\beta ) a^{-3}}} a^{-\frac{3}{2}} , \end{aligned}$$\end{document}V2(ϕ)=V0α+56βa-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&V_2(\phi ) = V_0\left( \alpha +\frac{5}{6}\beta a^{-1}\right) . \end{aligned}$$\end{document}Consider α=0.7,β=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7, \beta =-0.05$$\end{document} for numerically solving the above two equations, the two solutions are represented in Figs. 4 and 5, respectively. From Fig. 4, we can find ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} increases with a, and the potential increases with the increasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}, Eq. (54) implies that the potential will reach the maximum value V2(ϕ)=V0α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2(\phi )=V_0\alpha $$\end{document} in the future. From Fig. 5, ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} decreases with a, and the potential increases with decreasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}; in the future the potential will reach the maximum value V2(ϕ)=V0α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2(\phi )=V_0\alpha $$\end{document}. In Fig. 6, we plot the evolution of the density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} in the scenario of the phantom. Notice that the energy density in the phantom field becomes cosmologically dominant only in the recent past. In the future, the field comes to rest at the maximum of the potential and the universe eventually settles in a de Sitter phase.
The solution of Eqs. (53) and (54) corresponding to a plus sign in Eq. (53). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values α=0.7,β=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =-0.05$$\end{document}
The solution of Eqs. (53) and (54) corresponding to a minus sign in Eq. (53). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values α=0.7,β=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =-0.05$$\end{document}
Evolution of the density parameters in the phantom field (Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document}) and matter (Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document}) for model 1. Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} is indicated by a solid line, and Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document} is indicated by a dashed line. We consider the values α=0.7,β=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha =0.7,\beta =-0.05$$\end{document}
Model 2
We write down the EoS of the scalar fields for model 2:ωde=Pscalarfieldρscalarfield=-1-13δ(1+z)-1γ+δ(1+z)-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \omega _{\mathrm{de}}=\frac{ P_{\mathrm{scalar~field}}}{\rho _{\mathrm{scalar~field}}}=-1-\frac{\frac{1}{3}\delta (1+z)^{-1}}{\gamma +\delta (1+z)^{-1}}. \end{aligned}$$\end{document}It is obvious that only when parameters (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}) have opposite signs, there will be a singularity occurring when z=-1-δγ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z=-1-\frac{\delta }{\gamma }$$\end{document}. In Sect. 4, data-fitting results will show that such a singularity would not appear at low shift. In the quintessence scenario, compare Eqs. (31) and (32) with Eqs. (18) and (17), and we can obtain-Pc-34Pda=12ϕ˙2+V1(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&-P_c-\frac{3}{4}P_d a = \frac{1}{2}\dot{\phi }^2+V_1(\phi ) ,\end{aligned}$$\end{document}Pc+Pda=12ϕ˙2-V1(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&P_c+P_d a = \frac{1}{2}\dot{\phi }^2-V_1(\phi ) . \end{aligned}$$\end{document}Simplify the above two equations, referring to Eqs. (6), (7) and (19)–(22), replace the model parameters (Pc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_c$$\end{document}, Pd\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_d$$\end{document}) with the redefined parameters (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}), then12ϕ˙2=-16ρ0δa,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\frac{1}{2}\dot{\phi }^2 = -\frac{1}{6}\rho _0\delta a ,\end{aligned}$$\end{document}V1(ϕ)=ρ0γ+76ρ0δa.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&V_1(\phi ) = \rho _0\gamma +\frac{7}{6}\rho _0\delta a . \end{aligned}$$\end{document}From Eq. (58), it is easy to find that in the scenario of quintessence δ<0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta <0$$\end{document}. By Eqs. (58) and (59), the kinetic energy 12ϕ˙2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2$$\end{document} and potential V1(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1(\phi )$$\end{document} of the quintessence field are constructed with the parameters (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}) of model 2. Simplify these two equations, we havedϕda=±MP-δγ+δa+(1-γ-δ)a-3a-12,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\frac{\mathrm{d}\phi }{\mathrm{d}a}= \pm M_P\sqrt{\frac{-\delta }{\gamma +\delta a+(1-\gamma -\delta ) a^{-3}}} a^{-\frac{1}{2}} ,\end{aligned}$$\end{document}V1(ϕ)=V0γ+76δa,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&V_1(\phi ) = V_0\left( \gamma +\frac{7}{6}\delta a\right) , \end{aligned}$$\end{document}where V0=ρ0=3MP2H02\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_0=\rho _0=3M_P^2H_0^2$$\end{document}. Choose the parameters γ=0.7,δ=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =-0.05$$\end{document} for numerically solving the above two equations, the two solutions are represented in Figs. 7 and 8, respectively. From Fig. 7, we can find that ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} increases with a, and the potential decreases with increasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. From Fig. 8, we can find that ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} decreases with a, and the potential decreases with decreasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. Notice that since δ<0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta <0$$\end{document} in the scenario of quintessence, according to Eqs. (6), (21), and (45), the Friedmann equation is written as H2=13MP2ρ0[γ+δa+(1-γ-δ)a-3]\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$H^2=\frac{1}{3M_P^2}\rho _0[\gamma +\delta a+(1-\gamma -\delta )a^{-3}]$$\end{document}, which will not hold when the scale factor a is very large. Nevertheless at low redshift the relation is still feasible.
The solution of Eqs. (60) and (61) corresponding to a plus sign in Eq. (60). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values γ=0.7,δ=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =-0.05$$\end{document}
The solution of Eqs. (60) and (61) corresponding to a minus sign in Eq. (60). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_1$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values γ=0.7,δ=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =-0.05$$\end{document}
By Eqs. (23) and (24), we can obtain the expression of the density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} for model 2:Ωϕ=Ω1+Ω2=γ+δaE2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Omega _\phi =\Omega _1+\Omega _2=\frac{\gamma +\delta a}{E^2}, \end{aligned}$$\end{document}and the evolution curve of density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} in the scenario of quintessence has been plotted in Fig. 9, from which we see that the quintessence field begins to dominate at low redshift.
Evolution of the density parameters in the quintessence field (Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document}) and matter (Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document}) for model 2. Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} is indicated by a solid line, and Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document} is indicated by a dashed line. We consider values γ=0.7,δ=-0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =-0.05$$\end{document}
In order to realize model 2 in a phantom scenario, comparing Eqs. (37) and (38) with Eqs. (18) (17), we can obtain-Pc-34Pda=-12ϕ˙2+V2(ϕ),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&-P_c-\frac{3}{4}P_d a = -\frac{1}{2}\dot{\phi }^2+V_2(\phi ) ,\end{aligned}$$\end{document}Pc+Pda=-12ϕ˙2-V2(ϕ).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&P_c+P_d a = -\frac{1}{2}\dot{\phi }^2-V_2(\phi ) . \end{aligned}$$\end{document}Simplify and replace the model parameters (Pc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_c$$\end{document}, Pd\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P_d$$\end{document}) with the redefined parameters (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}), we have12ϕ˙2=16ρ0δa,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\frac{1}{2}\dot{\phi }^2 = \frac{1}{6}\rho _0\delta a , \end{aligned}$$\end{document}V2(ϕ)=ρ0γ+76ρ0δa.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&V_2(\phi ) = \rho _0\gamma +\frac{7}{6}\rho _0\delta a . \end{aligned}$$\end{document}From Eqs. (65) and (66), it is easy to find that in the scenario of phantom δ>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta >0$$\end{document}. By the above two equations, one can construct the kinetic energy 12ϕ˙2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2$$\end{document} and the potential V2(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2(\phi )$$\end{document} of the phantom field with parameters (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}) of model 2. Equations (65) and (66) can be rewritten asdϕda=±MPδγ+δa+(1-γ-δ)a-3a-12,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \frac{\mathrm{d}\phi }{\mathrm{d}a}=\pm M_P\sqrt{\frac{\delta }{\gamma +\delta a+(1-\gamma -\delta ) a^{-3}}} a^{-\frac{1}{2}}, \end{aligned}$$\end{document}V1(ϕ)=V0γ+76δa.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} V_1(\phi )=V_0\left( \gamma +\frac{7}{6}\delta a\right) . \end{aligned}$$\end{document}Choose the parameters γ=0.7,δ=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =0.05$$\end{document} for numerically solving the above two equations, the two solutions are represented in Figs. 10 and 11, respectively. From Fig. 10, we can find ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} increases with a, and the potential increases with the increasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. In Fig. 11, ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} decreases with a, and the potential increases with decreasing ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}. Notice that since δ>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta >0$$\end{document} in the scenario of phantom, the Friedmann equation can be written as H2=13MP2ρ0[γ+δa+(1-γ-δ)a-3]\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$H^2=\frac{1}{3M_P^2}\rho _0[\gamma +\delta a+(1-\gamma -\delta )a^{-3}]$$\end{document}, H→∞\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$H\rightarrow \infty $$\end{document} as a→∞\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a\rightarrow \infty $$\end{document}, which means there will be a “rip” in the future.
The solution of Eqs. (67) and (68) corresponding to a plus sign in Eq. (67). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values γ=0.7,δ=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =0.05$$\end{document}
The solution of Eqs. (67) and (68) corresponding to a minus sign in Eq. (67). The field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} as a function of a is depicted in the top panel, the potential V2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_2$$\end{document} as a function of ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document} is depicted in the bottom panel. The arrow indicates the direction of the evolution of the potential with respect to time. We consider values γ=0.7,δ=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =0.05$$\end{document}
In Fig. 12, we plot the evolution curve of density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} in the scenario of phantom. Note that the phantom becomes cosmologically dominant only in the recent past, finally the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} is less than -1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-1$$\end{document} [ωde=-43\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=-\frac{4}{3}$$\end{document}; see Eq. (55)] and the universe eventually settles in a “rip”.
Evolution of the density parameters in the phantom field (Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document}) and matter (Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document}) for model 2. Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} is indicated by a solid line, and Ωm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _m$$\end{document} is indicated by a dashed line. We consider values γ=0.7,δ=0.05\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma =0.7,\delta =0.05$$\end{document}
Astrophysical data constraints Type Ia supernovae
In this paper we use the Union2.1 SNe Ia data sets without systematic errors for data fitting, which compiles 580 SNe Ia covering the redshift range z=[0.015,1.4]\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z=[0.015,1.4]$$\end{document}. To perform the chi-square statistics, the theoretical distance modulus is defined asμth(zi)≡5log10DL(zi)+μ0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mu _{\mathrm{th}}(z_i)\equiv 5\log _{10}D_L(z_i)+\mu _0, \end{aligned}$$\end{document}where μ0≡42.39-5log10h\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu _0\equiv 42.39-5\log _{10}h$$\end{document} with h the Hubble parameter in units of 100km/s/Mpc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$100\,\mathrm{km/s/Mpc}$$\end{document},DL=(1+z)∫0zdz′E(z′;θ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} D_L= (1+z)\int ^z_0\frac{\mathrm{d}z'}{E(z';\theta )} \end{aligned}$$\end{document}is the Hubble-free luminosity distance in a spatially flat FRW universe, E(z;θ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$E(z;\theta )$$\end{document} is the dimensionless Hubble parameter, and θ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\theta $$\end{document} is for the model parameters.
The baryon acoustic oscillation (BAO) data sets are listed in Table 1. We use the parameter A to measure the BAO peak in the distribution of SDSS luminous red galaxies. In the following A is defined asA≡Ωm0E(zb)-131zb∫0zbdz′E(z′)23,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} A\equiv \sqrt{\Omega _{m0}}E(z_b)^{-\frac{1}{3}}\left[ \frac{1}{z_b}\int ^{z_b}_0\frac{\mathrm{d}z'}{E(z')}\right] ^{\frac{2}{3}}, \end{aligned}$$\end{document}where zb=0.35\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z_b=0.35$$\end{document}. The χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^2$$\end{document} for the BAO data isχBAO2=∑i=16[Aobs(zi)-Ath(zi;θ)]2σA2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \chi ^2_{\mathrm{BAO}}= \sum ^{6}_{i=1}\frac{[A_{\mathrm{obs}}(z_i)-A_{\mathrm{th}}(z_i;\theta )]^2}{\sigma ^2_A}. \end{aligned}$$\end{document}The total χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^2$$\end{document} is given byχ2=χ~SN2+χBAO2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \chi ^2 = \widetilde{\chi }^2_{\mathrm{SN}}+\chi ^2_{\mathrm{BAO}}. \end{aligned}$$\end{document}The fitting results and corresponding reduced χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^2$$\end{document} for model 1 and model 2 are listed in Table 2. The likelihoods of the parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}) and (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}) are shown in Figs. 13 and 14, respectively. Besides, the evolution of the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} with respect to the redshift z with 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} error propagation from data fitting (Table 2) are shown in Figs. 15 and 16, respectively.
1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\sigma $$\end{document} confidence ranges for parameter pair (α,β)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(\alpha ,\beta )$$\end{document} of model 1, constrained by SNe Ia and BAO data sets. The dotted straight line (β=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta =0$$\end{document}) corresponds to a Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM model. The blue dotted line and the red dotted line correspond to Ωm0=0.26\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{m0}=0.26$$\end{document} and Ωm0=0.32\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{m0}=0.32$$\end{document}, respectively
1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} and 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\sigma $$\end{document} confidence ranges for parameter pair (γ,δ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(\gamma ,\delta )$$\end{document} of model 2, constrained by SNe Ia and BAO data sets. The dotted straight line (δ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta =0$$\end{document}) corresponds to a Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM model. The blue dotted line and the red dotted line correspond to Ωm0=0.26\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{m0}=0.26$$\end{document} and Ωm0=0.32\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{m0}=0.32$$\end{document}, respectively
Evolution of the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} as a function of the redshift z with 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} error propagation, constrained by SNe Ia and BAO data sets for model 1. The solid line, the straight dotted line, and the light blue region represent the best-fit, ωde=-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=-1$$\end{document}(Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM), and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} region, respectively
Evolution of the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} as a function of the redshift z with 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} error propagation, constrained by SNe Ia and BAO data sets for model 2. The solid line, the straight dotted line, and the light blue region represent the best-fit, ωde=-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}=-1$$\end{document}(Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM), and 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} region, respectively
Six measurement points of the baryon acoustic oscillation data-sets
Parameters of model 1 and model 2 estimated by SNe Ia and BAO data sets with 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} errors
Since the observational confirmation on late-stage accelerative expansion of the universe many years ago, different models have been proposed to explain its source, among which parameterization is a widely used scheme to better characterize the dark energy and compare with observational results. In this paper, we studied two models parameterizing the effective pressure at low redshift, P(z)=Pa+Pbz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P(z)=P_a+P_b z$$\end{document} and P(z)=Pc+Pd1+z\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$P(z)=P_c+\frac{P_d}{1+z}$$\end{document}.
Deviations from the Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}CDM can be realized through different physical scenarios. Roughly speaking, there are two ways. One is to introduce some small but nonzero components besides the cosmological constant Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document}, such as imperfect fluid cosmology [38–41] and cosmic strings [35, 42]; whereas the other is to assume the cosmological constant Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document} exactly zero and the dark energy characterized by scalar fields evolving with time. In this paper, we pick the second way. We presented two parameterizations in the scenarios of quintessence and phantom fields, and accordingly expressed the kinetic energy term 12ϕ˙2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{1}{2}\dot{\phi }^2$$\end{document} and the potential term V(ϕ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V(\phi )$$\end{document} with the model parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}) and (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}), respectively. Then we reconstructed the density parameter Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _\phi $$\end{document} for quintessence and phantom evolving with redshift. In order to obtain a better physical understanding of the field and the potential, we numerically solved the field as a function of the scale factor a and the potential as a function of the field ϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\phi $$\end{document}.
We constrained the model parameters (α\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\alpha $$\end{document}, β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}) and (γ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma $$\end{document}, δ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta $$\end{document}) with the SNe Ia and BAO data sets. We reconstructed the evolution of the EoS parameter ωde\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}$$\end{document} in terms of the redshift z. For model 1, the value for the EoS parameter ωde0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}0} $$\end{document} is -1.027-0.0430.043\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-1.027^{0.043}_{ -0.043}$$\end{document} at present; for model 2, ωde0=-1.037-0.050+0.050\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}0} =-1.037^{+0.050}_{-0.050}$$\end{document}. These results show that model 1 and model 2 both slightly indicate that the EoS parameter of dark energy ωde<-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega _{\mathrm{de}}<-1$$\end{document}, which corresponds to a phantom dark energy scenario at present. Still, we cannot rule out a quintessence dark energy scenario or a Λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Lambda $$\end{document} dark energy scenario.
Different parameterizations possess their own advantages in addressing some particular problems, but their validity may not be ensured when applied to the explanation of the global evolution. For example, our two parameterizations of the effective pressure can estimate the deviation from the prediction of the standard model at a low redshift with a generality that does not depend on the concrete physical mechanism working in the background.
Acknowledgments
We are grateful for Jiaxin Wang’s instruction on model building and data analyzing. We also appreciate Prof. S. D. Odintsov’s recommendation of Refs. [7, 8, 18, 19] and Prof. V. K. Onemli’s recommendation of Refs. [13–15].
ReferencesP.A.R. Ade et al., [Plank Collaboration], arXiv:1502.01589v2 [astro-ph.CO]FujiiYPhys. Rev. D198226258068672110.1103/PhysRevD.26.25801982PhRvD..26.2580FFordLHPhys. Rev. D198735233910.1103/PhysRevD.35.23391987PhRvD..35.2339FWetterichCNucl. Phys. B198830266810.1016/0550-3213(88)90193-91988NuPhB.302..668WRatraBPeeblesPJEPhys. Rev. D198837340610.1103/PhysRevD.37.34061988PhRvD..37.3406RCaldwellRRDaveRSteinhardtPJPhys. Rev. Lett.199880158210.1103/PhysRevLett.80.15821998PhRvL..80.1582CNojiriSOdintsovSDInt. J. Geom. Methods Mod. Phys.20074115230381510.1142/S02198878070019281112.83047BambaKCapozzielloSNojiriSOdintsovSDAstrophys. Space Sci.201234215510.1007/s10509-012-1181-82012Ap&SS.342..155BChibaTOkabeTYamaguchiMPhys. Rev. D20006202351110.1103/PhysRevD.62.0235112000PhRvD..62b3511CArmendariz-PiconCMukhanovVFSteinhardtPJPhys. Rev. Lett.200085443810.1103/PhysRevLett.85.44382000PhRvL..85.4438AArmendariz-PiconCMukhanovVFSteinhardtPJPhys. Rev. D20016310351010.1103/PhysRevD.63.1035102001PhRvD..63j3510AR.R. Caldwell, Phys. Lett. B 545 (2002)OnemliVKWoodard ClassRPQuantum Gravity200219460710.1088/0264-9381/19/17/3111023.83025OnemliVKWoodardRPPhys. Rev. D20047010730110.1103/PhysRevD.70.1073012004PhRvD..70j7301OKahyaEOOnemliVKPhys Rev. D20077604351210.1103/PhysRevD.76.0435122007PhRvD..76d3512KSinghPSamiMDadhichNPhys. Rev. D20036802352210.1103/PhysRevD.68.0235222003PhRvD..68b3522SSamiMToporenskyAMod. Phys. Lett. A200419150910.1142/S02177323040139212004MPLA...19.1509SNojiriSOdintsovSDPhys. Lett. B2003562147198287110.1016/S0370-2693(03)00594-X2003PhLB..562..147N1027.83543NojiriSOdintsovSDPhys. Rev. D200368123512200396910.1103/PhysRevD.68.1235122003PhRvD..68l3512NMengX-HWangPClass. Quantum Gravity2003204949201953410.1088/0264-9381/20/22/0182003CQGra..20.4949M1054.83026X.-H. Meng, P. Wang, Quantum Gravity 21, 951 (2004)X.-H. Meng, P. Wang, Quantum Gravity 22, 23 (2005)HanDWangJ-XMengX-HEur. Phys. J. C2013222543WetterichCAstron. Astrophys.19953013211995A&A...301..321WAmendolaLPhys. Rev. D20006204351110.1103/PhysRevD.62.0435112000PhRvD..62d3511AP. Wang, X.-H. Meng, Class. Quantum Gravity 22, 283 (2005)HutererDTurnerMSPhys. Rev. D19996008130110.1103/PhysRevD.60.0813011999PhRvD..60h1301HWellerJAlbrechtAJPhys. Rev. D20026510351210.1103/PhysRevD.65.1035122002PhRvD..65j3512WChevallierMPolarskiDInt. J. Mod. Phys. D20011021310.1142/S02182718010008222001IJMPD..10..213CLinderEVPhys. Rev. Lett.20039009130110.1103/PhysRevLett.90.0913012003PhRvL..90i1301LEfstathiouGMon. Not. R. Astron. Soc.2000342810JassalHKBaglaJSPadmanabhanTMon. Not. R. Astron. Soc.2005356L1110.1111/j.1745-3933.2005.08577.x2005MNRAS.356L..11JWeiHYanX-PZhouY-NJCAP201401045319447410.1088/1475-7516/2014/01/0452014JCAP...01..045WSenAAPhys. Rev. D20087704350810.1103/PhysRevD.77.0435082008PhRvD..77d3508SS. Kumar, A. Nautiyal, A.A. Sen, Eur. Phys. J. C 73, 2562 (2013). arXiv:1207.4024v2CarrollSMHoffmanMPhys. Rev. D20036802350910.1103/PhysRevD.68.0235092003PhRvD..68b3509CBlakeCMon. Not. R. Astron. Soc.2011418170710.1111/j.1365-2966.2011.19592.x2011MNRAS.418.1707BMengX-HMaZ-YEur. Phys. J. C201272205310.1140/epjc/s10052-012-2053-72012EPJC...72.2053MMengX-HRenJMing-GuangHCommun. Theor. Phys.20074737938410.1088/0253-6102/47/2/0362007CoTPh..47..379MDouXMengX-HAdv. Astron.2011201182934010.1155/2011/8293402011AdAst2011E...1DWangJ-XMengX-HMod. Phys. Lett. A2014291450009316422210.1142/S02177323145000962014MPLA...2950009WR.J. Nemiroff, B. Patla, Am. J. Phys. 76, 265–276 (2008), arXiv:0703.739v2 (2007)