32 Page 6 of 26 Eur. Phys. J. C (2020) 80 :32
p
DE
=−
3
2
1
8π G
0
˙
G
G
σ. (16)
This effective anisotropic pressure (in the presence of shear)
induced by the Jordan field in the Jordan frame is absent in
the Einstein frame (see, e.g., Ref. [132]).
3 Exact solution compatible with standard CDM
model
We followthe method, given in [100], for obtaining exactcos-
mological solutions of a scalar-tensor gravity theory compat-
ible with the CDM model, albeit we consider LRS Bianchi
type I spacetime rather than spatially flat Robertson–Walker
spacetime.
6,7
To do so, we first re-express the set of differen-
tial equations given by Eqs. (4)–(7)using
dz
dt
=−H (1 + z),
where H =
˙
S/S is the average Hubble parameter and z is
the cosmic redshift we define in terms of S as z =−1 +
1
S
.
Accordingly, we re-express the energy density equation (4)
as
3H
2
1 −
2ω
3
(1 + z)
2
ϕ
2
ϕ
2
− 2(1 + z)
ϕ
ϕ
−
σ
2
2
− 2ωM
2
=
4
ϕ
2
ρ, (17)
and obtain the pressure equation as
−
dH
2
dz
1 + z − (1 + z)
2
ϕ
ϕ
+3H
2
1 + (1 + z)
2
2
3
ϕ
ϕ
+
2
3
(1 + ω)
ϕ
2
ϕ
2
−
2
3
(1 + z)
ϕ
ϕ
+
σ
2
2
− 2ωM
2
=−
4
ϕ
2
p, (18)
6
We note that in the GR limit (say, ϕ = const.) the anisotropic model
under consideration here, (4)–(7), reduces mathematically to the stan-
dard CDM+stiff matter model in the RW framework [133]. However,
the role of the stiff matter with positive energy density in [133]isplayed
by the shear scalar here, so that one can straightforwardly utilize the
rich class of solutions given in [133] for the GR limit in our model.
This presents a good example that one could find various solutions of
the model under consideration here but yet we focus on the solution
obtained by extending the method given in [100] to LRS Bianchi type
I metric.
7
We followed the method given in [100] for obtaining an exact cosmo-
logical solution of BD gravity theory for the expansion rate H (z),with
its CDM model counterpart up to a large redshift, viz., for pressure-
less matter in a spatially flat, homogeneous and isotropic universe. A
more general exact solution of the same setup was given in [134](much
earlier than [100]), was nevertheless given for the scale factor in cosmic
time and very complicated for extracting an exact H (z) required for
observational analyses.
using Eqs. (5) and (6), the shear propagation equation as
σ
σ
−
3
1 + z
+ 2
ϕ
ϕ
= 0, (19)
by subtracting Eq. (5) from Eq. (6), and finally re-express the
scalar field equation (7)as
−
dH
2
dz
ω(1 + z)
2
ϕ
ϕ
+
3
2
(1 + z)
+H
2
−2ω
ϕ
ϕ
(1 + z)
2
+ 4ω(1 + z)
ϕ
ϕ
+ 6
+
σ
2
2
− 2ωM
2
= 0, (20)
where
denotes derivative with respect to redshift (d/dz).
We note that Eqs. (18) and (20) have exactly the same
mathematical form, that is, a first order linear differential
equation in H
2
, such as
A
i
(z)
dH
2
dz
+ B
i
(z)H
2
+ C
i
(z) + D
i
= 0, i = 1, 2, (21)
provided that we set p = p
m
= 0, viz., the universe is filled
with only pressureless matter. We next note that constants
D
1
= D
2
=−2ωM
2
, and that the shear scalar σ
2
is the term
that differs in our system of equations from the one given in
[100]. Fortunately, it contributes to (18) and (20)inthesame
way, namely, as C
1
(z) = C
2
(z) =
σ
2
2
. In accordance with
these points, we assume A
1
(z) = A
2
(z) in (18) and (20) and
then solve for the rate of change of Jordan field in z as
ϕ
ϕ
=−
1
2(1 +ω)
1
1 + z
, (22)
which in turn renders the coefficients of H
2
identical, i.e,
B
1
(z) = B
2
(z), as well. Thus, integrating (22), it turns out
that the solution of the Jordan field is
ϕ
2
= ϕ
2
0
(1 + z)
−
1
1+ω
, (23)
where ϕ
0
= ϕ(z = 0) is the present time value of the Jordan
field, and which in turn gives
G = G
0
(1 + z)
1
1+ω
. (24)
The integration of the shear propagation equation given in
Eq. (19)gives
σ
2
= σ
2
0
(1 + z)
6
ϕ
ϕ
0
−4
, (25)
123