USTC-ICTS-07-03
A String-Inspired Quintom Model Of Dark Energy
Yi-Fu Caia, Mingzhe Lib,c, Jian-Xin Lud, Yun-Song Piaoe, Taotao
Qiua,
Xinmin Zhanga
aInstitute of High Energy Physics, Chinese Academy of Sciences, P.O. Box 918-4, Beijing 100049, P. R. China
bInstitut fЭr Theoretische Physik, Philosophenweg 16, D-69120 Heidelberg, Germany
cFakultДt fЭr Physik, UniversitДt Bielefeld, D-33615 Bielefeld, Germany
dInterdisciplinary Center for Theoretical Study, University of Science and Technology of China, Hefei, Anhui 230026, China
eCollege of Physical Sciences, Graduate School of Chinese Academy of Sciences,
YuQuan Road 19A, Beijing 100049, China
We propose in this paper a quintom model of dark energy with a single scalar field φ given by the
lagrangian
= -V (φ)
. In the limit of β′ →0 our model reduces to the
effective low energy lagrangian of tachyon considered in the literature. We study the cosmological
evolution of this model, and show explicitly the behaviors of the equation of state crossing the
cosmological constant boundary.
I. INTRODUCTION
The current data from type Ia supernovae, cosmic microwave background (CMB) radiation, and other cosmological
observations[1–4] have provided strong evidences for a spatially flat and accelerated expanding universe at
the present time. Within the framework of the standard cosmology, this acceleration can be understood by
introducing a mysterious component, dubbed dark energy (DE). The simplest candidate for DE is a minor
positive cosmological constant, but it suffers from the fine-tuning and coincidence problems. As a possible
solution to these problems various dynamical models of DE have been proposed, such as quintessence. In the
recent years with the accumulated astronomical observational data it becomes possible to probe the current and
even early behavior of DE. Although the current fits to the data in combination of the 3-year WMAP[5], the
recently released 182 SNIa Gold sample[6] and also other cosmological observational data show the consistence of
the cosmological constant, it is worth noting that the dynamical DE models are not excluded and a class of
dynamical models with equation of state across -1, dubbed quintom, is mildly favored (for recent references see e.g.
[7, 8]).
Theoretically it is a big challenge to the model building of the quintom dark energy. The Ref.[9] is the first
paper pointing out this challenge and showing explicitly the difficulty of realizing w crossing over -1 in the
quintessence and phantom like models. In general with a single fluid or a single scalar field with a lagrangian of
form
=
(φ,∂μφ∂μφ) it has been proved [10] (see also [11]) that the dark energy perturbation would be
divergent as the equation of state (EOS) w approaches to -1. This “no-go” theorem forbids the dynamical models
widely studied in the literature with a single scalar field such as quintessence, phantom and k-essence to make
the EOS cross over the cosmological constant boundary. The quintom scenario of dark energy is designed to
understand the nature of dark energy with w across -1. The quintom models of dark energy differ from the
quintessence, phantom and k-essence and so on in the determination of the cosmological evolution and the fate of the
universe.
To realize a viable quintom scenario of dark energy it needs to introduce extra degree of freedom to the conventional theory
with a single fluid or a single scalar field. The first model of quintom scenario of dark energy is given in Ref.[9] with two scalar
fields. This model has been studied in detail later on. In the recent years there have been active studies on models of quintom
like dark energy such as models with high derivative term[12], vector field[13], or even extended theory of gravity[14] and so
on, see e.g. [15].
In this paper, we propose a new type of quintom model inspired by the string theory. We will demonstrate in this paper our
model can realize the equation of state crossing -1 naturally. This paper is organized as follows. In section 2 we present our
model and study its properties especially on the conditions required for the model parameters when w crosses over -1. By
solving numerically the model we will study the evolution of the equation of state. The section 3 is the summary of our
paper.
II. OUR MODEL
Our model is given by the following action We have adopted signature (+,-,-,-) in this paper.
This model generalizes the usual“Born-Infeld-type” action for the effective description of tachyon dynamics by adding a term
φ□φ to the usual ∇μφ∇μφ in the square root. It is known that the 4D effective action of tachyon dynamics to the lowest order
in ∇μφ∇μφ around the top of the tachyon potential can be obtained by the stringy computations for either a
D3 brane in bosonic theory [16, 17] or a non-BPS D3 brane in supersymmetric theory [18]. To this order,
we have no need to include the operator φ□φ in the action since it is equivalent to the usual ∇μφ∇μφ term.
However, we cannot exclude its existence in an action such as the “Born-Infeld-type” one when incorporating an
infinite number of higher derivative terms since now the two terms are in general different dynamically. For
example, we cannot simply replace the φ□φ term in the action (1) by the ∇μφ∇μφ. Further this term in the
above generalized action has new cosmlogical consequence as will be shown in this paper. Including an infinite
number of higher derivative terms that has a significant cosmological consequence has also been discussed in the
context of p-adic string recently in [19]. The two parameters α′ and β′ in (1) can also be made arbitrary when
the background flux is turned on [20]. Without further reasoning, we will take the action (1) as our starting
point.
As will be demonstrated, the β′ term in (1) is crucial to realize the w across -1. There we have defined α′ = α∕M4 and
β′ = β∕M4 with α and β being the dimensionless parameters, respectively and M an energy scale used to make the “kinetic
energy terms” dimensionless. V (φ) is the potential of scalar field φ (e.g., a tachyon) with dimension of [mass]4 with an
expected tachyon potential behavior in general, i.e., bounded and reaching its minimum asymptotically, unless specifically
stated. Note that, □ =
∂μ
gμν
∂ν, therefore, in (1) the terms ∇μφ∇μφ and φ□φ both involve two fields and two
derivatives.
The model with operator φ□φ for the realization of w crossing -1 has been proposed in [12]. However in general
for a model with lagrangian as a sum of operators with a polynomial function of the scalar field φ and its
derivatives, the operator φ□φ can be rewritten as a total derivative term which makes no contribution after
integration and a term which renormalizes the canonical kinetic term as discussed above. So if one considers a
renormalizable lagrangian, the operator φ□φ will not be included. Ref.[12] considered a dimension-6 operator as (□φ)2.
However in the present model, the operator φ□φ appears at the same order as the operator ∇μφ∇μφ does
in the “Born-Infeld-type” action. As discussed above, this model appears more natural than the one used in
[12].
With (1) we obtain the equation of motion of the scalar field φ, where f =
and V φ = dV∕dφ. Following the convention in Ref.[12], the energy-momentum tensor
Tμν is given by the standard definition: δgμνS ≡-∫
d4x
Tμνδg
μν, where ψ ≡
= -
.
For a flat Friedman-Robertson-Walker (FRW) universe and a homogenous scalar field φ, the equation of motion in (2) can
be solved equivalently by the following two equations where we have made use of the ψ as defined before and the first equation above is just the defining equation for ψ in terms of
φ and its derivatives. H = ȧ∕a is the Hubble parameter.
One can see from equations above the β term plays a role in determining the evolution of the scalar field φ. We can read the
energy density from (3) as and similarly the pressure
With the Friedman equations H2 =
ρ and Ḣ = -4πG(ρ + p), we now study the cosmological evolution of equation of
state for the present model. Given w = p∕ρ, we have ρ + p = (1 + w)ρ. To explore the possibility of the w across -1, we need
to check if
(ρ + p)
0 can be held when w →-1.
Using (6) and (7) as well as making use of the defining equation for ψ, we have and from which, Eq. (8) implies that we have either (i)
= 0 or (ii)-![αV(φ)
f](article5c19x.png)
= β
[
] when w →-1.
Let us assume
= 0 first when w →-1. With this and Eq. (9), we have
(ρ + p) = β
[
] = -
β2φ2
□φ.
Therefore the conditions for having the w across over -1 are (i-1)φ
0; (i-2)
0; (i-3)
□φ
0 in addition to the
= 0. Since the
characteristic behavior of V (φ) is to have a finite maximum value at a finite φ and to reach its (exponentially) vanishing minimum
asymptotically
where one expects
= 0 and
□φ = 0 but a non-zero
, so realizing the above crossing over conditions must happen before
reaching the potential minimum asymptotically. This implies once the crossing over conditions are met, the field φ must
continue to run away as it should be since we have
0.
Let us turn to the second case (ii) when w →-1. We now have -
= β
(
) = β
(
)+βφ
(
) which is equivalent to saying
Y ≡ (α + β)![V
f-](article5c49x.png)
+ βφ
(
) = 0. For now, we have ![φ˙](article5c53x.png)
0. Then the above can be further expressed as
[(α + β)lnφ + β ln(
)] = 0
if φ
0.
Then from (9), the further condition for crossing over is
(ρ + p)
0 which implies Ẏ
0 or
[(α + β)lnφ + β ln(
)]
0. In
summary, the following conditions are required for the crossing over: ii-1)Y = 0 and ii-2) Ẏ
0 in addition to the ![˙φ](article5c66x.png)
0
and φ
0. Given the definition of the above Y and the characteristic behavior of V (φ) discussed previously,
one expects that both Y and Ẏ vanish asymptotically while
can reach a non-vanishing value. This implies
that the crossing over if occurring at all must occur before the V (φ) reaches its minimum asymptotically as
anticipated.
Before we demonstrate numerically that the crossing over of the present model can indeed be realized using specific
examples, let us remark one salient feature of the present model that the β term in the present model (1) is the key for
realizing the crossing over. One can check from (8) that if β = 0 the ω →-1 is possible only for
= 0. Then from (9), we will
have
(ρ + p) = 0, the impossibility of crossing over.
The analysis above shows various possibilities of our model in realizing the EOS w crossing -1. Now we consider some
specific examples for numerical calculations of the evolution of the EOS. In Figures 1 we take V (φ) = V 0e-λφ2
and plot the
behavior of EOS. In the numerical calculations we have normalized the values of the scalar field and V 0, respectively, by the
energy scale M. In Figure 1 our model predicts the EOS crossing -1 during the evolution and a big-rip singularity for the
fate of the universe. Numerically we have checked that
= 0 when w crosses over the cosmological constant
boundary.
In Figure 2 we take a different potential for numerical calculations. One can see that the EOS crosses over -1 during the
evolution. When we take β positive Figure 3 shows the EOS starts with w < -1, crosses over -1 into the region of w > -1,
then transits again to w < -1.
The potentials used for the numerical calculations in Figures 1-3 are well motivated by the string theory. However as
a phenomenological study of our model as a quintom dark energy we in Figure 4 plot the evolution of the
EOS w with a potential which is a sum of eλφ and e-λφ. This type of potential does not have the general
behavior of the tachyon potential, however has been used for the study of phenomenological models of dark
energy. One can see from this Figure that the EOS evolves from the region where w > -1 to w < -1, and
stays there for a period of time then comes back to w > -1. At late time our model gives rise to the de-sitter
phase.
III. CONCLUSION AND DISCUSSION
The current cosmological observations indicate the possibility that the acceleration of the universe is driven by
dark energy with EOS across -1, which if confirmed further in the future will challenge the theoretical model
building of the dark energy. In this paper we have proposed a string-inspired model of dark energy through
modifying the usual effective “Born-Infeld-type” description of tachyon dynamics. As shown in the present
work, this modification by including a β term in the action (1) is the key for the EOS crossing -1 during the
evolution.
Compared to other models with w across -1 in the literature so far the present one has a motivation inspired from
string theory consideration and is also economical in the sense that it involves a single scalar field such as a
tachyon.
Acknowledgments
We thank Bo Feng, Gongbo Zhao, Hong Li and Junqing Xia for useful discussions. This work is supported in part by National
Natural Science Foundation of China under Grant Nos. 90303004, 10533010 and 19925523. The author M.L.
would like to acknowledge the support by Alexander von Humboldt Foundation. JXL acknowledges support by
grants from the Chinese Academy of Sciences and grants from the NSF of China with Grant Nos: 10588503 and
10535060.
_______________________________________________
[1] S. Perlmutter et al., Astrophys. J. 483, 565 (1997); A. G. Riess et al., Astrophys. J. 116, 1009 (1998).
[2] D. N. Spergel et al., Astrophys. J. Suppl. 148, 175 (2003).
[3] A. G. Riess et al., Astrophys. J. 607, 665 (2004).
[4] U. Seljak et al., Phys. Rev. D71, 103515 (2005).
[5] D. N. Spergel et al., astro-ph/0603449.
[6] A. G. Riess et al., astro-ph/0611572.
[7] H. Li, M. Su, Z. Fan, Z. Dai, and X. Zhang, astro-ph/0612060.
[8] G.-B. Zhao et al, astro-ph/0612728; H. Wei, N. Tang, and S. N. Zhang, astro-ph/0612746; J. Zhang, X. Zhang, and
H. Liu, astro-ph/0612642; U. Alam, V. Sahni, and A. A. Starobinsky, astro-ph/0612381; Y.-G. Gong, and A. Wang,
astro-ph/0612196; V. Barger, Y. Gao, D. Marfatia, astro-ph/0611775.
[9] B. Feng, X. Wang and X. Zhang, Phys. Lett. B607, 35 (2005).
[10] G.-B. Zhao, J.-Q. Xia, M. Li, B. Feng, and X. Zhang, Phys. Rev. D72, 123515 (2005).
[11] R. R. Caldwell, M. Doran, Phys. Rev. D72 043527 (2005); A. Vikman, Phys. Rev. D71, 023515 (2005); W. Hu, Phys.
Rev. D71, 047301 (2005);
[12] M. Li, B. Feng, and X. Zhang, JCAP 0512, 002 (2005); X.-F. Zhang, and T.-T. Qiu, Phys. Lett. B642, 187 (2006).
[13] H. Wei, and R.-G. Cai, Phys. Rev. D73, 083002 (2006).
[14] R.-G. Cai, H.-S. Zhang, and A. Wang, Commun. Theor. Phys. 44, 948 (2005); P. S. Apostolopoulos, and N. Tetradis,
Phys. Rev. D74, 064021 (2006); H.-S. Zhang, and Zong-Hong Zhu, astro-ph/0611834.
[15] B. Feng, M. Li, Y. Piao, and X. Zhang, Phys. Lett. B634, 101 (2006); X.-F. Zhang, H. Li, Y.-S. Piao, and X. M. Zhang,
Mod. Phys. Lett. A21, 231 (2006); Z. Guo, Y. Piao, X. Zhang, and Y.-Z. Zhang, Phys. Lett. B608, 177 (2005); H. Li, B.
Feng, J.-Q. Xia, and X. Zhang, Phys. Rev. D73, 103503 (2006); I. Y. Aref’eva, A. S. Koshelev, and S. Yu. Vernov, Phys.
Rev. D72, 064017 (2005); Z.-K. Guo, Y.-S. Piao, X. Zhang, and Y.-Z. Zhang, astro-ph/0608165; Y.-F. Cai, H. Li, Y.-S.
Piao, and X. Zhang, gr-qc/0609039; W. Zhao, and Y. Zhang, Phys. Rev. D73, 123509 (2006).
[16] A. A. Gerasimov and S. L. Shatashvili, JHEP 0010, 034 (2000).
[17] D. Kutasov, M. Marino and G. W. Moore, JHEP 0010, 045 (2000).
[18] D. Kutasov, M. Marino and G. W. Moore, arXiv:hep-th/0010108.
[19] N. Barnaby, T. Biswas and J. M. Cline, arXiv:hep-th/0612230.
[20] P. Mukhopadhyay and A. Sen, JHEP 0211, 047 (2002).
[21] A. Sen, JHEP 0204, 048 (2002); A. Sen, JHEP 0207, 065 (2002); A. Sen, Mod. Phys. Lett. A17, 1797 (2002); A. Sen,
Int. J. Mod. Phys. A18, 4869 (2003).
[22] G. Gibbons, Phys. Lett. B537, 1 (2002); S. Mukohyama, Phys. Rev. D66, 024009 (2002); D. Choudhury, D. Ghoshal, D.
P. Jatkar, and S. Panda, Phys. Lett. B544, 231 (2002); T. Padmanabhan, Phys. Rev. D66, 021301 (2002); J. Hao, and
X. Li, Phys. Rev. D66, 087301 (2002); J. S. Bagla, H. K. Jassal, and T. Padmanabhan, Phys. Rev. D67, 063504 (2003);
E. J. Copeland, M. R. Garousi, M. Sami, and S. Tsujikawa, Phys. Rev. D71, 043003 (2005); K.-F. Zhang, W. Fang, and
H.-Q. Lu, Int. J. Theor. Phys. 45, 1296, 2006.