# The Effect of de-Sitter Like Background on Increasing the Zero Point Budget of Dark Energy

Authors: Haidar Sheikhahmadi; Ali Aghamohammadi; Khaled Saaidi
Journal: Advances in High Energy Physics
 (2016)

---

## Abstract

During this work, using subtraction renormalization mechanism, zero point quantum fluctuations for bosonic scalar fields in a de-Sitter like background are investigated. By virtue of the observed value for spectral index,n
s
(
k
), for massive scalar field the best value for the first slow roll parameter, ϵ, is achieved. In addition, the energy density of vacuum quantum fluctuations for massless scalar field is obtained. The effects of these fluctuations on other components of the universe are studied. By solving the conservation equation, for some different examples, the energy density for different components of the universe is obtained. In the case which all components of the universe are in an interaction, the different dissipation functions, Q
~
i, are considered. The time evolution of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) shows that Q
~
=
3
γ
H
(
t
)
ρ
m has the best agreement in comparison to observational data including CMB, BAO, and SNeIa data set.


---

## Body

## 1. Introduction

For a de-Sitter like background, to investigate the effects of the quantum fluctuations on the energy budget of the universe both massive and massless scalar fields are considered. Because of the appearance of the first slow roll parameter, the relation between power spectral index and slow roll parameters leads to a way to connect theoretical results with observations. In fact one can estimate the best value of the first slow roll parameter based on observations, for instance, Planck 2013 [1]. From this perspective, for zero point quantum fluctuations both energy density and pressure are calculated. Interestingly it is observed that the contribution of zero point energy is increased against the results for the normal de-Sitter ones [2, 3]. The importance of this note is that whereas the zero point contribution of dark energy potentially is detectable, therefore the possibility of dark energy detection is increased. Also it should be stressed, due to time dependency of Hubble parameter which appeared in energy density of zero point quantum fluctuations, energy can transfer between different components of the universe. By considering this concept one can propose different manners to investigate interaction between zero point fluctuations and other sectors. In a general case where all components are in an interaction, three different dissipation functions Q
~ are considered. For this special case time evolution of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) shows that Q
~
=
3
γ
H
(
t
)
ρ
m is in best agreement in comparison to observations including CMB, BAO, and SNeIa data set. Furthermore, energy density of matter, ρ
m, has some deviations in comparison to the result of ordinary de-Sitter one. In fact in ρ
m, beside a
-
3, some extra terms appeared, where these new terms could be interpreted as new source for cold dark matter risen from interaction between matter and quantum fluctuations. In fact as it was discussed in [3] zero point quantum fluctuations could be considered as subdark energy, and these extra terms in question can be proposed as subdark matter. The scheme of this paper is as follows.In Section2 the general framework of this work including the mathematical calculations is discussed. In Section 3 the cosmological role for zero point energy density is investigated, and the results of this work are compared with previous works. In Section 4 to estimate the amount of subdark energy and also subdark matter, the bounds which are risen from time evolution of dark energy are considered. And at last, we have conclusion.


## 2. Massive Scalar Field and Slow Roll Parameters

To study the effect of zero point quantum fluctuations let us consider a real minimally coupled bosonic scalar fieldΦ in a semiclassical general relativity mechanism [4]. In such scenario the geometry is not quantized but the energy momentum tensor related to the scalar field is calculated by means of the vacuum expectation value concept. To begin we consider the action (1)
S
=
S
m
+
∫
d
4
x
-
g
κ
2
R
-
2
Λ
+
κ
2
L
Φ
,where S
m is the matter action, g is the determinant of the metric, R is the Ricci scalar, Λ is the Einstein’s cosmological constant, and L
Φ is defined as (2)
L
Φ
=
-
1
2
g
μ
ν
∇
μ
Φ
∇
ν
Φ
-
V
Φ
,where V
(
Φ
) is the potential of the model. Variation of the action (1) with respect to metric yields (3)
R
μ
ν
-
1
2
g
μ
ν
R
+
g
μ
ν
Λ
=
κ
2
T
μ
ν
M
+
0
T
μ
ν
Φ
0
r
e
n
,where T
μ
ν
M
=
(
-
2
/
-
g
)
×
(
δ
S
m
/
δ
g
μ
ν
) and (4)
T
μ
ν
Φ
=
∂
μ
Φ
∂
ν
Φ
-
1
2
g
μ
ν
∂
α
Φ
∂
α
Φ
+
2
V
Φ
.Using the above equation, the 00 and i
i components of energy-momentum tensor read (5)
T
0
Φ
0
=
-
1
2
Φ
˙
2
+
Φ
,
i
2
a
2
t
+
2
V
Φ
,
1
3
T
i
Φ
i
=
1
2
Φ
˙
2
-
1
3
Φ
,
i
2
a
2
t
-
2
V
Φ
.It is also obvious that variation (1) with respect to Φ yields (6)
□
Φ
-
V
,
Φ
Φ
=
0
,where □ is d’Alembert operator. As it is mentioned above in semiclassical approach one can quantize the scalar field and therefore the vacuum expectation value of energy density and pressure could be obtained. For this purpose, the quantized scalar field is defined as (7)
Φ
t
,
x
=
∫
d
3
k
2
π
3
2
k
a
k
ϕ
t
e
-
j
k
·
x
+
a
k
†
ϕ
∗
t
e
j
k
·
x
,where ϕ
(
t
) is a function which should be determined, a
k (a
k
†) is annihilation (creation) operators and -
1
=
j [5, 6]. Now introducing (7) to (6), it expresses (8)
ϕ
¨
+
3
H
t
ϕ
˙
+
k
2
a
2
ϕ
+
V
,
ϕ
=
0
,where H
=
a
˙
/
a, and a
(
t
) is the scale factor of the universe, overdot denotes derivation with respect to the cosmic time, and the k
2
ϕ
/
a
2 term appeared due to scalar field’s spatial dependency. For more convenience one can use conformal time η
=
∫
d
t
/
a
(
t
), and therefore (8) could be rewritten as (9)
ϕ
′
′
+
2
H
ϕ
′
+
k
2
ϕ
+
a
2
V
′
ϕ
ϕ
′
=
0
,in which overprime indicates derivation with respect to η, and H
=
a
′
/
a
=
a
˙ is the conformal, comoving, Hubble parameter. Now, by defining χ
≔
ϕ
/
a, using a power law potential, V
(
ϕ
)
=
m
2
ϕ
2
/
2, relation (9) reads (10)
χ
′
′
+
k
2
+
m
2
a
2
-
a
′
′
a
χ
=
0
.It should be stressed for de-Sitter like background in question that the scale factor is taken as a
(
t
)
=
t
1
/
ϵ, where ϵ
=
-
H
˙
/
H
2
≪
1 is the first slow roll parameter [7, 8]. Using definitions of scale factor (i.e., a
(
η
)
≃
-
(
1
+
ϵ
)
/
H
η) and second slow roll parameter, τ
=
m
2
/
3
H
2, one can attain (11)
-
m
2
a
2
+
a
′
′
a
≃
1
η
2
2
+
3
ϵ
-
τ
.By introducing ν
~
≔
(
ϵ
-
τ
)
+
3
/
2 and keeping only up to first orders of the first and second slow roll parameters (10) could be rewritten as (12)
χ
′
′
+
k
2
-
ν
~
2
-
1
/
4
η
2
χ
=
0
,where ν
~
2
≃
3
(
ϵ
-
τ
)
+
9
/
4. Solving the above Bessel like equation, the magnitude of ϕ can be achieved as (13)
ϕ
≃
H
1
-
ϵ
k
k
a
H
1
-
ϵ
3
/
2
-
ν
~
.To estimate the best value of ϵ, one can use the power spectrum concept for instance in light of Planck 2013 [1]. It is obvious that, by considering the Fourier transformation for an arbitrary function as (14)
Φ
~
k
t
,
r
=
∫
d
3
k
2
π
3
2
k
e
j
k
·
r
ϕ
~
k
t
,where k is comoving momentum and r is spatial vector, the power spectrum can be defined as (15)
ϕ
~
k
t
ϕ
~
k
′
∗
t
=
2
π
2
k
2
ρ
S
2
π
3
δ
3
k
-
k
′
,where ρ
S is power spectrum in question and ϕ
~
k
(
t
)
ϕ
~
k
′
∗
(
t
) indicates the mean square value of ϕ
~
k
(
t
). By combining (7), (14), and (15) the power spectrum can be achieved as (16)
ρ
S
=
k
2
2
π
2
ϕ
t
2
.To achieve this result the relation between annihilation and creation operators reads [
a
k
,
a
k
′
†
]
=
(
2
π
)
3
δ
3
(
k
-
k
′
). At last if one considers the relation (13), the power spectrum could be rewritten as (17)
ρ
S
=
1
-
2
ϵ
H
2
k
a
H
1
-
ϵ
3
-
2
ν
~
.Also another important quantity is notable that is the spectral index which plays a crucial role in inflationary investigations and is defined as (18)
n
S
-
1
=
d
ln
⁡
ρ
S
d
ln
⁡
k
.Substituting (17) into the above equation, one has (19)
n
S
-
1
=
1
-
2
ϵ
3
-
2
ν
~
,and at last after some manipulations the spectral index is obtained as n
S
-
1
=
-
2
ϵ
+
2
τ. By means of observed value of spectral index, n
S
<
0.9675, risen from Planck 2013 [1], the best value of the first slow roll parameter is ϵ
≃
0.02.


### 2.1. Typical Example: Massless Scalar Field

In this case we want to consider massless scalar field. Therefore if, in action (1), one assumes V
(
Φ
)
=
0, solving (9), for positive modes, is attained as (20)
ϕ
η
=
1
a
η
1
-
i
ξ
k
η
e
-
i
k
η
,where ξ
=
(
2
+
3
ϵ
)
/
2. In this case, as massive ones, the first slow roll parameter is considered only up to the first order. For this typical case we want to estimate the zero point quantum fluctuation contribution in the energy budget of the universe. To begin we have to calculate the vacuum expectation value of the scalar field energy-momentum tensor. By means of (5) and (20), one has (21)
ρ
v
a
c
=
0
T
00
ϕ
0
=
1
2
∫
0
a
t
Λ
c
d
3
k
2
π
3
2
k
ϕ
˙
2
+
k
2
a
2
t
ϕ
2
,
P
v
a
c
=
1
3
Σ
i
0
T
i
ϕ
i
0
=
1
2
∫
0
a
t
Λ
c
d
3
k
2
π
3
2
k
ϕ
˙
2
-
k
2
3
a
2
t
ϕ
2
,where ρ
v
a
c and P
v
a
c, indicate the energy density and pressure of the vacuum quantum fluctuations, respectively. According to quantum field theory the cutoff Λ
c should be considered greater than physical momenta k
/
a
(
t
). By considering (20) and definition of scale factor for a de-Sitter like background, one has (22)
ϕ
˙
2
+
k
2
a
2
t
ϕ
2
=
1
a
4
2
k
2
+
a
2
H
2
1
+
9
ϵ
+
O
ϵ
2
,
ϕ
˙
2
-
k
2
3
a
2
t
ϕ
2
=
2
3
a
4
k
2
-
a
2
H
2
2
1
-
3
ϵ
+
O
ϵ
2
.By virtue of (22), solving (21) leads to (23)
ρ
v
a
c
=
Λ
c
4
16
π
2
+
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
,
P
v
a
c
=
Λ
c
4
48
π
2
-
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
.The first terms in (23) are the contribution of the energy density and pressure for Minkowskian space time; and because of the cutoff dependency the latter terms are well known bare quantities. To get rid of quartic divergencies the subtraction mechanism is a good suggestion, which is close to Casimir approach [9]. The base of the Casimir effect is on the subtraction mechanism. In this effect the energy that gives rise to observable effects is the difference between vacuum energy risen from two parallel conducting plates, for example, for electromagnetic field and vacuum energy computed in Minkowskian space-time in an infinite volume. Therefore when two infinite values for energy of space time and plates are subtracted one can obtain a finite quantity. Therefore by means of Arnowitt-Deser-Misner (ADM) approach [2] and subtraction mechanism one concludes that the vacuum energy to an asymptotically flat space-time with metric g
μ
ν can be achieved as E
=
H
G
R
(
g
μ
ν
)
-
H
G
R
(
η
μ
ν
), where H
G
R refers to the Hamiltonian which is calculated in general relativity. This equation indicates that flat space-time does not gravitate and the contribution of the energy which is obtained in Minkowski space-time can be subtracted from related quantity in curved background [2, 9–11]. Thence, because flat space-time does not gravitate one is able to subtract the contribution of quartic terms in Minkowski space from the same terms in FriedmannLemaitreRobertsonWalker (FLRW) space-time. Also it should be stressed that the results which were obtained for de-Sitter like scenario have some notable differences with normal de-Sitter ones. The first is, for de-Sitter like investigations, the Hubble parameter is not a constant and this causes the appearance of the first slow roll parameter in the model which could be considered to investigate the accuracy of this model. As a second note, the coefficients which appeared in energy density and pressure cause the increase of the zero point quantum fluctuations contribution in the dark energy. Now let us investigate the bare quantities which are (24)
ρ
b
a
r
e
=
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
,
P
b
a
r
e
=
-
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
.Following [3] one can introduce counter terms for energy density and pressure, respectively, as follows:(25)
ρ
c
o
u
n
t
=
-
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
+
ρ
Z
,
(26)
P
c
o
u
n
t
=
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
+
P
Z
,where (27)
ρ
Z
=
1
+
9
ϵ
H
2
t
M
2
16
π
2
,
(28)
P
Z
=
-
1
-
3
ϵ
H
2
t
M
2
48
π
2
,where subscript Z refers to the zero point and M is in order of Planck mass. The interesting result which could be realized from (27) and (28) is that they appear to contain a mix of ultraviolet (UV) cutoff, M as Planck mass, and infrared (IR) cutoff, and H
(
t
) as the inverse of horizon size. Based on quantum field theory concepts, we know the widely separated energy scales are decoupled; but one can find that the origin of this mix of UV-IR cutoffs arises from classical subtraction procedure, that is, E
=
H
G
R
(
g
μ
ν
)
-
H
G
R
(
η
μ
ν
). In fact the origin of this mix goes back to (21), where it involves T
μ
ν
ϕ and 0
T
μ
ν
ϕ
0 as local and global quantities, respectively. For more discussion we refer the reader to [2, 11, 12]. Using (27) and (28), the equation of state (EoS) parameter for the vacuum fluctuations could be expressed as (29)
ω
Z
=
P
Z
ρ
Z
=
-
1
3
+
4
ϵ
.This relation indicates that the EoS of zero point quantum fluctuations is dependent on the first slow roll parameter. It should be noted also that whereas this approach is similar to Casimir mechanism both positive and negative signs for energy density are acceptable. To consider this fact one can consider coefficient σ
=
±
1, for (27), and redefines it as (30)
ρ
Z
=
σ
1
+
9
ϵ
H
2
t
M
2
16
π
2
.The positive sign causes an attractive force and negative ones are related to the repulsive case. Now for more discussions about time dependency of energy density of vacuum fluctuations, one is able to redefine it based on critical energy density of the universe. Hence considering definition of critical energy density (ρ
c
r
i
=
3
H
2
(
t
)
/
8
π
G) and by means of definition of Planck mass M
P
l
=
1
/
G, (G is the Newtonian constant), the energy density of zero point fluctuations could be rewritten as (31)
ρ
Z
=
Ω
Z
ρ
c
r
i
+
O
ϵ
2
=
β
1
+
9
ϵ
ϵ
2
a
-
2
ϵ
t
,where Ω
Z
=
β
(
1
+
9
ϵ
)
/
M
P
l
2 and β
=
σ
M
2
/
16
π
2. Because of the time dependency of ρ
Z
(
t
), the conservation equation for ρ
Z is not satisfied. Hence using ρ
˙
Z
=
-
2
ϵ
H
(
t
)
ρ
Z and (31) one has (32)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
Q
,where Q is dissipation function and it could be obtained as (33)
Q
=
2
H
t
ρ
Z
1
+
5
ϵ
.Therefore the energy density of quantum fluctuations is capable of exchanging energy with other components of the universe. To investigate the transformation of energy we consider some different cases as follows.


## 2.1. Typical Example: Massless Scalar Field

In this case we want to consider massless scalar field. Therefore if, in action (1), one assumes V
(
Φ
)
=
0, solving (9), for positive modes, is attained as (20)
ϕ
η
=
1
a
η
1
-
i
ξ
k
η
e
-
i
k
η
,where ξ
=
(
2
+
3
ϵ
)
/
2. In this case, as massive ones, the first slow roll parameter is considered only up to the first order. For this typical case we want to estimate the zero point quantum fluctuation contribution in the energy budget of the universe. To begin we have to calculate the vacuum expectation value of the scalar field energy-momentum tensor. By means of (5) and (20), one has (21)
ρ
v
a
c
=
0
T
00
ϕ
0
=
1
2
∫
0
a
t
Λ
c
d
3
k
2
π
3
2
k
ϕ
˙
2
+
k
2
a
2
t
ϕ
2
,
P
v
a
c
=
1
3
Σ
i
0
T
i
ϕ
i
0
=
1
2
∫
0
a
t
Λ
c
d
3
k
2
π
3
2
k
ϕ
˙
2
-
k
2
3
a
2
t
ϕ
2
,where ρ
v
a
c and P
v
a
c, indicate the energy density and pressure of the vacuum quantum fluctuations, respectively. According to quantum field theory the cutoff Λ
c should be considered greater than physical momenta k
/
a
(
t
). By considering (20) and definition of scale factor for a de-Sitter like background, one has (22)
ϕ
˙
2
+
k
2
a
2
t
ϕ
2
=
1
a
4
2
k
2
+
a
2
H
2
1
+
9
ϵ
+
O
ϵ
2
,
ϕ
˙
2
-
k
2
3
a
2
t
ϕ
2
=
2
3
a
4
k
2
-
a
2
H
2
2
1
-
3
ϵ
+
O
ϵ
2
.By virtue of (22), solving (21) leads to (23)
ρ
v
a
c
=
Λ
c
4
16
π
2
+
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
,
P
v
a
c
=
Λ
c
4
48
π
2
-
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
.The first terms in (23) are the contribution of the energy density and pressure for Minkowskian space time; and because of the cutoff dependency the latter terms are well known bare quantities. To get rid of quartic divergencies the subtraction mechanism is a good suggestion, which is close to Casimir approach [9]. The base of the Casimir effect is on the subtraction mechanism. In this effect the energy that gives rise to observable effects is the difference between vacuum energy risen from two parallel conducting plates, for example, for electromagnetic field and vacuum energy computed in Minkowskian space-time in an infinite volume. Therefore when two infinite values for energy of space time and plates are subtracted one can obtain a finite quantity. Therefore by means of Arnowitt-Deser-Misner (ADM) approach [2] and subtraction mechanism one concludes that the vacuum energy to an asymptotically flat space-time with metric g
μ
ν can be achieved as E
=
H
G
R
(
g
μ
ν
)
-
H
G
R
(
η
μ
ν
), where H
G
R refers to the Hamiltonian which is calculated in general relativity. This equation indicates that flat space-time does not gravitate and the contribution of the energy which is obtained in Minkowski space-time can be subtracted from related quantity in curved background [2, 9–11]. Thence, because flat space-time does not gravitate one is able to subtract the contribution of quartic terms in Minkowski space from the same terms in FriedmannLemaitreRobertsonWalker (FLRW) space-time. Also it should be stressed that the results which were obtained for de-Sitter like scenario have some notable differences with normal de-Sitter ones. The first is, for de-Sitter like investigations, the Hubble parameter is not a constant and this causes the appearance of the first slow roll parameter in the model which could be considered to investigate the accuracy of this model. As a second note, the coefficients which appeared in energy density and pressure cause the increase of the zero point quantum fluctuations contribution in the dark energy. Now let us investigate the bare quantities which are (24)
ρ
b
a
r
e
=
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
,
P
b
a
r
e
=
-
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
.Following [3] one can introduce counter terms for energy density and pressure, respectively, as follows:(25)
ρ
c
o
u
n
t
=
-
1
+
9
ϵ
H
2
t
Λ
c
2
16
π
2
+
ρ
Z
,
(26)
P
c
o
u
n
t
=
1
-
3
ϵ
H
2
t
Λ
c
2
48
π
2
+
P
Z
,where (27)
ρ
Z
=
1
+
9
ϵ
H
2
t
M
2
16
π
2
,
(28)
P
Z
=
-
1
-
3
ϵ
H
2
t
M
2
48
π
2
,where subscript Z refers to the zero point and M is in order of Planck mass. The interesting result which could be realized from (27) and (28) is that they appear to contain a mix of ultraviolet (UV) cutoff, M as Planck mass, and infrared (IR) cutoff, and H
(
t
) as the inverse of horizon size. Based on quantum field theory concepts, we know the widely separated energy scales are decoupled; but one can find that the origin of this mix of UV-IR cutoffs arises from classical subtraction procedure, that is, E
=
H
G
R
(
g
μ
ν
)
-
H
G
R
(
η
μ
ν
). In fact the origin of this mix goes back to (21), where it involves T
μ
ν
ϕ and 0
T
μ
ν
ϕ
0 as local and global quantities, respectively. For more discussion we refer the reader to [2, 11, 12]. Using (27) and (28), the equation of state (EoS) parameter for the vacuum fluctuations could be expressed as (29)
ω
Z
=
P
Z
ρ
Z
=
-
1
3
+
4
ϵ
.This relation indicates that the EoS of zero point quantum fluctuations is dependent on the first slow roll parameter. It should be noted also that whereas this approach is similar to Casimir mechanism both positive and negative signs for energy density are acceptable. To consider this fact one can consider coefficient σ
=
±
1, for (27), and redefines it as (30)
ρ
Z
=
σ
1
+
9
ϵ
H
2
t
M
2
16
π
2
.The positive sign causes an attractive force and negative ones are related to the repulsive case. Now for more discussions about time dependency of energy density of vacuum fluctuations, one is able to redefine it based on critical energy density of the universe. Hence considering definition of critical energy density (ρ
c
r
i
=
3
H
2
(
t
)
/
8
π
G) and by means of definition of Planck mass M
P
l
=
1
/
G, (G is the Newtonian constant), the energy density of zero point fluctuations could be rewritten as (31)
ρ
Z
=
Ω
Z
ρ
c
r
i
+
O
ϵ
2
=
β
1
+
9
ϵ
ϵ
2
a
-
2
ϵ
t
,where Ω
Z
=
β
(
1
+
9
ϵ
)
/
M
P
l
2 and β
=
σ
M
2
/
16
π
2. Because of the time dependency of ρ
Z
(
t
), the conservation equation for ρ
Z is not satisfied. Hence using ρ
˙
Z
=
-
2
ϵ
H
(
t
)
ρ
Z and (31) one has (32)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
Q
,where Q is dissipation function and it could be obtained as (33)
Q
=
2
H
t
ρ
Z
1
+
5
ϵ
.Therefore the energy density of quantum fluctuations is capable of exchanging energy with other components of the universe. To investigate the transformation of energy we consider some different cases as follows.


## 3. Transformation of Energy between Different Components of the Universe

### 3.1. Transformation of Energy between Zero Point Fluctuations and Matter

In this case one has(34)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
2
H
t
ρ
Z
1
+
5
ϵ
,
ρ
˙
m
+
3
H
t
ρ
m
=
-
2
H
t
ρ
Z
1
+
5
ϵ
,and therefore the combination of the two sections of the above equation yields (35)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
+
ρ
˙
m
+
3
H
t
ρ
m
=
0
,which indicates that the conservation equation in general is satisfied. Therefore by means of (31) and (35), ρ
m could be achieved as (36)
ρ
m
=
Ψ
a
-
2
ϵ
t
+
Ψ
~
a
-
3
t
,where Ψ
=
-
2
β
(
1
+
44
ϵ
/
3
)
/
3
ϵ
2 and Ψ
~ is integration constant. From (36) it is realized that in our model the matter density equation is modified, where the first term indicates the matters which are risen from interaction of quantum fluctuations with matter and the latter indicates the remaining contribution of matter, namely, ordinary cold dark matter.


### 3.2. Transformation of Energy between Zero Point Energy and the Remnant Components of Dark Energy

Assume there is an internal interaction between different components of dark energy, namely,ρ
Λ and ρ
Z, where ρ
Λ indicates energy density of cosmological constant. In this case, one can suppose that ω
Λ
=
-
1 and therefore the conservation equation reads (37)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
+
ρ
˙
Λ
=
0
.By virtue of ρ
˙
Z
=
-
2
β
(
1
+
3
ϵ
)
/
(
ϵ
2
a
3
ϵ
(
t
)
) and definition of scale factor in de-Sitter like background, one has (38)
ρ
Λ
=
β
ϵ
1
+
14
ϵ
H
2
t
+
C
0
,where C
0 is integration constant. In addition it is obvious that because ρ
˙
Z is proportional to ρ
˙
Λ the Big Bang Nucleosynthesis (BBN) constraint which has been discussed in [2] could be considered to estimate the upper bound on Ω
Z
=
β
(
1
+
9
ϵ
)
/
M
P
l
2. Also it should be stressed that, by comparing ρ
Z with the one in normal de-Sitter model, it is clear that the coefficient (
1
+
9
ϵ
) causes the increase of the magnitude of zero point energy density.


### 3.3. Transformation of Energy between All Components of the Universe

In this stage, one can consider a general case where all components of the universe are in an interaction. Therefore the conservation equations could be written as(39)
ρ
˙
m
+
3
H
t
ρ
m
=
Q
~
i
,
(40)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
Q
~
i
,where Q
~
i are dissipation functions and are defined as follows: (i)
Q
~
1
=
3
κ
H
(
t
)
ρ
Λ,
(ii)
Q
~
2
=
3
γ
H
(
t
)
ρ
m,
(iii)
Q
~
3
=
3
θ
H
(
t
)
ρ
Z.A also κ, γ, and θ indicate the strength of the interaction between different components of the universe [13].


#### 3.3.1. Solving Conservation Equation forQ
~
1

By virtue of (40) and dissipation function Q
~
1, the conservation equation for dark energy components of the universe is as follows:(41)
ρ
˙
Λ
+
ρ
˙
Z
+
3
κ
H
t
ρ
Λ
+
3
H
t
ρ
Z
1
+
ω
Z
=
0
.Based on (31) and (29), (41) could be rewritten as (42)
ρ
˙
Λ
+
3
κ
H
t
ρ
Λ
+
2
β
1
+
14
ϵ
H
3
t
=
0
,and hence solving this differential equation yields (43)
ρ
Λ
=
-
D
˘
3
κ
-
2
ϵ
a
-
2
ϵ
t
+
c
~
a
-
3
κ
,where D
˘
=
2
(
1
+
14
ϵ
)
β
/
ϵ
2 and c
~ is integration constant. By substituting (43) into (39) one can attain ρ
m as follows (44)
ρ
m
=
-
D
˘
3
-
2
ϵ
1
-
1
/
κ
a
-
2
ϵ
t
+
κ
c
~
1
-
κ
a
-
3
κ
+
B
˘
a
-
3
t
.In the above equation B
˘ is integration constant. From this relation one can conclude that, in matter equation, only ordinary cold dark matter does not appear; rather, an extra term appears which is risen from interaction of matter and quantum fluctuations, namely, subdark matter. In the following, we will come back to this issue.


#### 3.3.2. Solving Conservation Equation forQ
~
2

By rewriting (39) and (40) for Q
~
2, one arrives to(45)
ρ
˙
m
+
3
H
t
ρ
m
=
3
γ
H
t
ρ
m
,
(46)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
γ
H
t
ρ
m
,and thus (45) could be considered as (47)
ρ
˙
m
+
3
H
t
1
-
γ
ρ
m
=
0
,and after solving, one can attain (48)
ρ
m
=
ρ
m
0
a
-
3
1
-
γ
t
.Substituting (48) into (46) yields (49)
ρ
Λ
=
β
1
+
14
ϵ
ϵ
H
2
t
+
ρ
m
0
γ
3
1
-
γ
a
-
3
1
-
γ
+
M
,where M is the integration constant.


#### 3.3.3. Solving Conservation Equation forQ
~
3

In this stage, one can suppose that the interaction between different components of the universe is determined by virtue ofQ
~
3. Thence (39) and (40) are rearranged as (50)
ρ
˙
m
+
3
H
t
ρ
m
=
3
θ
H
t
ρ
Z
,
(51)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
θ
H
t
ρ
Z
.Using definition of scale factor in de-Sitter like background, one is allowed to rewrite (50) as (52)
ρ
˙
m
+
3
ϵ
ρ
m
a
-
ϵ
-
3
θ
ϵ
β
1
+
3
ϵ
ϵ
2
a
-
3
ϵ
=
0
,and solving this equation for ρ
m yields (53)
ρ
m
=
A
˘
a
-
2
ϵ
+
K
˘
a
-
3
t
,where A
˘
≃
3
θ
β
(
1
+
8
ϵ
)
/
ϵ
2 and K
˘ is the integration constant. Hence, to attain ρ
Λ, from (51) one has (54)
ρ
˙
Λ
+
β
ϵ
3
2
+
3
θ
+
27
ϵ
1
+
θ
a
-
3
ϵ
t
=
0
,and solution of this equation for ρ
Λ reads (55)
ρ
Λ
=
3
β
ϵ
2
+
3
θ
+
27
ϵ
1
+
θ
H
2
t
.


## 3.1. Transformation of Energy between Zero Point Fluctuations and Matter

In this case one has(34)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
2
H
t
ρ
Z
1
+
5
ϵ
,
ρ
˙
m
+
3
H
t
ρ
m
=
-
2
H
t
ρ
Z
1
+
5
ϵ
,and therefore the combination of the two sections of the above equation yields (35)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
+
ρ
˙
m
+
3
H
t
ρ
m
=
0
,which indicates that the conservation equation in general is satisfied. Therefore by means of (31) and (35), ρ
m could be achieved as (36)
ρ
m
=
Ψ
a
-
2
ϵ
t
+
Ψ
~
a
-
3
t
,where Ψ
=
-
2
β
(
1
+
44
ϵ
/
3
)
/
3
ϵ
2 and Ψ
~ is integration constant. From (36) it is realized that in our model the matter density equation is modified, where the first term indicates the matters which are risen from interaction of quantum fluctuations with matter and the latter indicates the remaining contribution of matter, namely, ordinary cold dark matter.


## 3.2. Transformation of Energy between Zero Point Energy and the Remnant Components of Dark Energy

Assume there is an internal interaction between different components of dark energy, namely,ρ
Λ and ρ
Z, where ρ
Λ indicates energy density of cosmological constant. In this case, one can suppose that ω
Λ
=
-
1 and therefore the conservation equation reads (37)
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
+
ρ
˙
Λ
=
0
.By virtue of ρ
˙
Z
=
-
2
β
(
1
+
3
ϵ
)
/
(
ϵ
2
a
3
ϵ
(
t
)
) and definition of scale factor in de-Sitter like background, one has (38)
ρ
Λ
=
β
ϵ
1
+
14
ϵ
H
2
t
+
C
0
,where C
0 is integration constant. In addition it is obvious that because ρ
˙
Z is proportional to ρ
˙
Λ the Big Bang Nucleosynthesis (BBN) constraint which has been discussed in [2] could be considered to estimate the upper bound on Ω
Z
=
β
(
1
+
9
ϵ
)
/
M
P
l
2. Also it should be stressed that, by comparing ρ
Z with the one in normal de-Sitter model, it is clear that the coefficient (
1
+
9
ϵ
) causes the increase of the magnitude of zero point energy density.


## 3.3. Transformation of Energy between All Components of the Universe

In this stage, one can consider a general case where all components of the universe are in an interaction. Therefore the conservation equations could be written as(39)
ρ
˙
m
+
3
H
t
ρ
m
=
Q
~
i
,
(40)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
Q
~
i
,where Q
~
i are dissipation functions and are defined as follows: (i)
Q
~
1
=
3
κ
H
(
t
)
ρ
Λ,
(ii)
Q
~
2
=
3
γ
H
(
t
)
ρ
m,
(iii)
Q
~
3
=
3
θ
H
(
t
)
ρ
Z.A also κ, γ, and θ indicate the strength of the interaction between different components of the universe [13].


### 3.3.1. Solving Conservation Equation forQ
~
1

By virtue of (40) and dissipation function Q
~
1, the conservation equation for dark energy components of the universe is as follows:(41)
ρ
˙
Λ
+
ρ
˙
Z
+
3
κ
H
t
ρ
Λ
+
3
H
t
ρ
Z
1
+
ω
Z
=
0
.Based on (31) and (29), (41) could be rewritten as (42)
ρ
˙
Λ
+
3
κ
H
t
ρ
Λ
+
2
β
1
+
14
ϵ
H
3
t
=
0
,and hence solving this differential equation yields (43)
ρ
Λ
=
-
D
˘
3
κ
-
2
ϵ
a
-
2
ϵ
t
+
c
~
a
-
3
κ
,where D
˘
=
2
(
1
+
14
ϵ
)
β
/
ϵ
2 and c
~ is integration constant. By substituting (43) into (39) one can attain ρ
m as follows (44)
ρ
m
=
-
D
˘
3
-
2
ϵ
1
-
1
/
κ
a
-
2
ϵ
t
+
κ
c
~
1
-
κ
a
-
3
κ
+
B
˘
a
-
3
t
.In the above equation B
˘ is integration constant. From this relation one can conclude that, in matter equation, only ordinary cold dark matter does not appear; rather, an extra term appears which is risen from interaction of matter and quantum fluctuations, namely, subdark matter. In the following, we will come back to this issue.


### 3.3.2. Solving Conservation Equation forQ
~
2

By rewriting (39) and (40) for Q
~
2, one arrives to(45)
ρ
˙
m
+
3
H
t
ρ
m
=
3
γ
H
t
ρ
m
,
(46)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
γ
H
t
ρ
m
,and thus (45) could be considered as (47)
ρ
˙
m
+
3
H
t
1
-
γ
ρ
m
=
0
,and after solving, one can attain (48)
ρ
m
=
ρ
m
0
a
-
3
1
-
γ
t
.Substituting (48) into (46) yields (49)
ρ
Λ
=
β
1
+
14
ϵ
ϵ
H
2
t
+
ρ
m
0
γ
3
1
-
γ
a
-
3
1
-
γ
+
M
,where M is the integration constant.


### 3.3.3. Solving Conservation Equation forQ
~
3

In this stage, one can suppose that the interaction between different components of the universe is determined by virtue ofQ
~
3. Thence (39) and (40) are rearranged as (50)
ρ
˙
m
+
3
H
t
ρ
m
=
3
θ
H
t
ρ
Z
,
(51)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
θ
H
t
ρ
Z
.Using definition of scale factor in de-Sitter like background, one is allowed to rewrite (50) as (52)
ρ
˙
m
+
3
ϵ
ρ
m
a
-
ϵ
-
3
θ
ϵ
β
1
+
3
ϵ
ϵ
2
a
-
3
ϵ
=
0
,and solving this equation for ρ
m yields (53)
ρ
m
=
A
˘
a
-
2
ϵ
+
K
˘
a
-
3
t
,where A
˘
≃
3
θ
β
(
1
+
8
ϵ
)
/
ϵ
2 and K
˘ is the integration constant. Hence, to attain ρ
Λ, from (51) one has (54)
ρ
˙
Λ
+
β
ϵ
3
2
+
3
θ
+
27
ϵ
1
+
θ
a
-
3
ϵ
t
=
0
,and solution of this equation for ρ
Λ reads (55)
ρ
Λ
=
3
β
ϵ
2
+
3
θ
+
27
ϵ
1
+
θ
H
2
t
.


## 3.3.1. Solving Conservation Equation forQ
~
1

By virtue of (40) and dissipation function Q
~
1, the conservation equation for dark energy components of the universe is as follows:(41)
ρ
˙
Λ
+
ρ
˙
Z
+
3
κ
H
t
ρ
Λ
+
3
H
t
ρ
Z
1
+
ω
Z
=
0
.Based on (31) and (29), (41) could be rewritten as (42)
ρ
˙
Λ
+
3
κ
H
t
ρ
Λ
+
2
β
1
+
14
ϵ
H
3
t
=
0
,and hence solving this differential equation yields (43)
ρ
Λ
=
-
D
˘
3
κ
-
2
ϵ
a
-
2
ϵ
t
+
c
~
a
-
3
κ
,where D
˘
=
2
(
1
+
14
ϵ
)
β
/
ϵ
2 and c
~ is integration constant. By substituting (43) into (39) one can attain ρ
m as follows (44)
ρ
m
=
-
D
˘
3
-
2
ϵ
1
-
1
/
κ
a
-
2
ϵ
t
+
κ
c
~
1
-
κ
a
-
3
κ
+
B
˘
a
-
3
t
.In the above equation B
˘ is integration constant. From this relation one can conclude that, in matter equation, only ordinary cold dark matter does not appear; rather, an extra term appears which is risen from interaction of matter and quantum fluctuations, namely, subdark matter. In the following, we will come back to this issue.


## 3.3.2. Solving Conservation Equation forQ
~
2

By rewriting (39) and (40) for Q
~
2, one arrives to(45)
ρ
˙
m
+
3
H
t
ρ
m
=
3
γ
H
t
ρ
m
,
(46)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
γ
H
t
ρ
m
,and thus (45) could be considered as (47)
ρ
˙
m
+
3
H
t
1
-
γ
ρ
m
=
0
,and after solving, one can attain (48)
ρ
m
=
ρ
m
0
a
-
3
1
-
γ
t
.Substituting (48) into (46) yields (49)
ρ
Λ
=
β
1
+
14
ϵ
ϵ
H
2
t
+
ρ
m
0
γ
3
1
-
γ
a
-
3
1
-
γ
+
M
,where M is the integration constant.


## 3.3.3. Solving Conservation Equation forQ
~
3

In this stage, one can suppose that the interaction between different components of the universe is determined by virtue ofQ
~
3. Thence (39) and (40) are rearranged as (50)
ρ
˙
m
+
3
H
t
ρ
m
=
3
θ
H
t
ρ
Z
,
(51)
ρ
˙
Λ
+
ρ
˙
Z
+
3
H
t
ρ
Z
1
+
ω
Z
=
-
3
θ
H
t
ρ
Z
.Using definition of scale factor in de-Sitter like background, one is allowed to rewrite (50) as (52)
ρ
˙
m
+
3
ϵ
ρ
m
a
-
ϵ
-
3
θ
ϵ
β
1
+
3
ϵ
ϵ
2
a
-
3
ϵ
=
0
,and solving this equation for ρ
m yields (53)
ρ
m
=
A
˘
a
-
2
ϵ
+
K
˘
a
-
3
t
,where A
˘
≃
3
θ
β
(
1
+
8
ϵ
)
/
ϵ
2 and K
˘ is the integration constant. Hence, to attain ρ
Λ, from (51) one has (54)
ρ
˙
Λ
+
β
ϵ
3
2
+
3
θ
+
27
ϵ
1
+
θ
a
-
3
ϵ
t
=
0
,and solution of this equation for ρ
Λ reads (55)
ρ
Λ
=
3
β
ϵ
2
+
3
θ
+
27
ϵ
1
+
θ
H
2
t
.


## 4. Bounds Which Are Risen from Time Evolution of Dark Energy

In this section we want to compare the results of this work with results which are risen from standardΛCDM model. For this end, one can start from the Friedmann equation. Therefore, the ratio of dark energy density and critical energy density in the standard model as a function of red shift parameter, z, is obtained as (56)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Λ
1
+
z
3
1
+
ω
Λ
Ω
M
1
+
z
3
+
Ω
Λ
.It should be emphasized in the above equation that one can get Ω
Λ
=
0.73, Ω
M
=
0.23, and ω
Λ
=
-
0.98, which are obtained from a combination of CMB, BAO, and SNeIa data sets [2, 14–16]. Whereas, in this work, the components of dark energy are quantum fluctuations and cosmological constant, thence the Friedmann equation is obtained as (57)
3
H
2
t
8
π
G
=
ρ
m
t
+
ρ
D
E
=
ρ
m
t
+
ρ
Λ
t
+
Ω
Z
ρ
c
r
i
t
,where subscript DE indicates dark energy. By substituting (31) into (57), the Friedmann equation could be rewritten as (58)
H
2
t
=
H
0
2
t
1
-
Ω
z
Ω
m
t
+
Ω
Λ
t
,where Ω
i
(
t
)
=
ρ
i
(
t
)
/
ρ
c
r
i
(
0
), and i refers to m, Λ, and Z, respectively; in addition ρ
c
r
i
(
0
) denotes the critical energy density in present epoch. It should be noted that the energy densities of curvature and radiation are neglected. From (31) and definition of dimensionless energy density parameters, one gets (59)
ρ
D
E
=
Ω
Z
ρ
c
r
i
t
+
ρ
c
r
i
0
Ω
Λ
.By virtue of definition of ρ
c
r
i
(
t
), (58) could be rearranged as (60)
ρ
c
r
i
t
=
ρ
c
r
i
0
1
-
Ω
Z
Ω
m
t
+
Ω
Λ
.By dividing ρ
D
E and ρ
c
r
i
(
t
), one arrives to (61)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Z
+
1
-
Ω
Z
Ω
Λ
Ω
m
t
+
Ω
Λ
.Whereas different equations forρ
m (based on different conditions for interaction) are attained, Ω
m
(
t
) in (60) gets different forms. Thus by considering (36), (44), (53), and (61) it could be rewritten, respectively, as follows.(a)
From (36), consider (62)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Z
+
1
-
Ω
Z
Ω
D
E
-
Ω
Z
Ω
m
1
+
z
3
+
Ω
~
m
1
+
z
2
ϵ
+
Ω
D
E
-
Ω
Z
,
whereΩ
~
m
=
Ψ
/
ρ
c
r
i
(
0
) and Ω
m
=
Ψ
~
/
ρ
c
r
i
(
0
). In Figure 1, the evolution of (62) versus z parameter shows a deviation in comparison to the ordinary ρ
D
E
(
t
)
/
ρ
c
r
i
(
t
), (56) that illuminates the effects of Ω
Z and ϵ in the evolution of this function.
(b)
From (44), consider(63)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Z
+
1
-
Ω
Z
Ω
D
E
-
Ω
Z
Ω
m
1
+
z
3
+
Ω
-
m
1
+
z
2
ϵ
+
Ω
m
∗
1
+
z
3
κ
+
Ω
D
E
-
Ω
Z
,
whereΩ
-
m
=
-
D
˘
/
ρ
c
r
i
(
0
)
(
3
-
2
ϵ
(
1
-
1
/
κ
)
) and Ω
m
∗
=
κ
c
~
/
ρ
c
r
i
(
0
)
(
1
-
κ
). Considering the above equation, it is found out that the Q
~
1 case does not lead to physical result, because the equation for ρ
m can not satisfy observations as well.
(c)
From (48), consider (64)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Z
+
1
-
Ω
Z
Ω
D
E
-
Ω
Z
Ω
m
1
+
z
3
1
-
γ
+
Ω
D
E
-
Ω
Z
,
whereΩ
m
(
1
+
z
)
3
(
1
-
γ
) indicates dark matter component of the universe. Figure 2 shows a deviation with respect to ordinary ρ
D
E
(
t
)
/
ρ
c
r
i
(
t
) (56) where in comparison to Figure 1 it is realized that Q
2 function has better results.
(d)
From (53), consider (65)
ρ
D
E
t
ρ
c
r
i
t
=
Ω
Z
+
1
-
Ω
Z
Ω
D
E
-
Ω
Z
Ω
m
1
+
z
3
+
Ω
m
♯
1
+
z
2
ϵ
+
Ω
D
E
-
Ω
Z
,
whereΩ
m
=
K
˘
/
ρ
c
r
i
(
0
) and Ω
m
♯
=
A
˘
/
ρ
c
r
i
(
0
). The energy density parameter in this case behaves as the case related to (38). This manner is similar to case (a) that is (62), and therefore it does not need to plot its evolution. It should be noted that the behaviour of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) based on different quantities for Ω
Z is plotted in Figure 3.As it is obvious, Figure 3 illuminates the evolution of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) versus z (62) for three different values Ω
Z
=
0.16,0.05, and -
0.16. The results are shown with green dashed, red dotted, and blue solid lines, respectively. It is clear that ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) decreases with decreasing Ω
Z.Figure 1
ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) versus z have been shown with dashed-blue line (62) and solid-red line, respectively (56). The auxiliary parameters are Ω
Z
=
0.16, Ω
m
=
0.17, Ω
~
m
=
0.1, Ω
D
E
=
0.73, ϵ
=
0.02, and ω
Λ
=
-
0.98.Figure 2
ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) versus z have been shown with dotted line (64) and solid line (56). The auxiliary parameters are Ω
Z
=
0.16, Ω
m
=
0.23, Ω
D
E
=
0.73, γ
=
0.03, and ω
Λ
=
-
0.98.Figure 3
ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) versus z (62) for three different values Ω
Z
=
0.16, Ω
Z
=
0.05, and Ω
Z
=
-
0.16. They have been shown with dashed-green, dotted-red, and blue (solid) lines, respectively.


## 5. Conclusion and Discussion

Vacuum quantum fluctuations in a de-Sitter like background for both massive and massless bosonic scalar fields have been investigated. In light of Planck database2013 we have estimated the best value for the first slow roll parameter and using such quantity the different components of the Universe’s energy budget have been calculated. It should be stressed that the scalar fields have been quantized in a FLRW framework and it has been shown that the contribution of vacuum fluctuations has increased in such background in comparison with normal de-Sitter case. It should be emphasized that the subtraction approach has been used to eliminate the infinities which appeared in the calculations. Incidentally, using the physical energy density of zero point quantum fluctuation, it has been realized that this component of the universe has to have an interaction with other components of the universe. In addition, when the energy density of matter is achieved, it has been found that beside ordinary dark matter there exist matter components which were created due to interaction with zero point quantum fluctuations. Also whereas zero point energy density was time dependent, the transformation of energy between different ingredients of the universe has been investigated. It was considerable that, for the state in which all components of the universe exchange energy between themselves, time evolution of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) has shown that Q
2
=
3
γ
H
ρ
m is in best agreement in comparison with observational database and also the interaction term Q
~
1
=
3
κ
H
(
t
)
ρ
Λ had not any physical results. At last for more details, the bounds which have risen from time evolution of dark energy density in comparison to standard Λ cosmology have been investigated. To compare the results of this work with observational data, we have regarded the time evolution of ρ
D
E
(
z
)
/
ρ
c
r
i
(
z
) which was concluded from a combination of CMB, BAO, and SNeIa data sets. From Figures 2 and 3, the evolution of (62) and (64) versus z in comparison to observational results has been illustrated.



---

*Source: 2594189-2016-06-08.xml*