Dynamics of Bianchi I Universe in Extended Gravity with Scale Factors ()
1. Introduction
Einstein’s General Theory of Relativity (GR) is one of the revolutionary concepts in modern physics. This equation results out of the outstanding interplay
between matter and space-time. Based on Einstein equation
,
the presence of energy-momentum distribution has shown how gravity influenced space-time and how the curvature of space-time acts on gravity simultaneously. While (GR) is successful in several physical aspects, some astrophysical and cosmological issues still have not found their appropriate explanation by GR. Late time cosmic acceleration, is the most significant result that is indescribable in theory of GR. Of the proposed alternatives of the general theory of gravity, modified gravity theory is found to be one of the more appealing solutions. The recent cosmological observations of Type Ia supernovae (SN Ia) [1] [2] [3] , the cosmic ray Microwave Background (CMB) [4] [5] , Wilkinson Microwave Planck and the Sloan Digital Sky Survey [6] [7] strongly suggest the acceleration of the universe.
An unprecedented view of our dark universe has been shown by the recent discovery of the gravitational wave by LIGO (Laser Interferometer Gravitational- Wave Observatory). With the help of precise measurement of advanced LIGO, we will be able to test the differences between Einstein’s General Theory of Relativity and it’s extended theories [8] [9] .
Primary ideas of generalization of the GR have gone beyond the Einstein-Hilbert action. The generalization of interest here is including the term of
theory of gravity in the Einstein-Hilbert action. The presence of a late time cosmic acceleration of the universe can indeed be explained by
gravity [10] . Later on, Harko et al. [11] proposed some types of that
terms that can be included in Einstein-Hilbert action to describe the late time cosmology. After an essential suggestion of changing the action, there have been many other suggestions put forth in recent times that address critical issues in cosmology. A specific example is the
gravity [12] , which considered a relativistically co-variant model of interacting dark energy, based on the principle of the least action. Several cosmological models with different aspects have been constructed and analyzed in the last decade. Perlmutter et al. [13] has found the first evidence for the positive-cosmological-constant universe that shows expansion of universe based on observations of a set of 40 high-red-shift supernovae.
Bamba [14] has gone through different approaches of testing dark energy and alternative models of extended gravity by cosmography. An approach to work on this problem is to modify the regular part of the Einstein field equation where the standard Einstein-Hilbert action replaced by an arbitrary function of Ricci scalar R [15] . Yousaf [16] has discussed irregular energy density in the presence of imperfect fluid for
. Alves et al. [17] considered the extra polarization modes and found corresponding gravitational waves in
of the model. Zaregonbadi [18] solved the dark matter in galaxy scale in
gravity. Mishra et al. [19] obtained transition from early deceleration to late-time cosmic acceleration due to different deceleration parameter in the presence of hybrid scale factor.
Most studies on
gravity are with isotropic universes. However, Mishra et al. [20] [21] have studied the dynamical behavior of the anisotropic cosmological models in detail. A universe filled by perfect fluid in cylindrically symmetric space time in two classes of
gravity has been studied in [22] . In this respect, Halife [23] studied the anisotropic Bianchi type I universe in the presence of quark matter source in
with cosmological constant.
Naser [24] added new scale factor as
in Bianchi Type-V
Cosmology in
gravity. Several cosmological models were analysed on the perspective aspect of the anisotropic universe by considering recent observation of supernovae type I as a perfect standard candle and studied their spectrum [25] [26] . Harko et al. [11] considered three different classes of
gravity such as 1)
, 2)
and 3)
.
Simultaneously [27] investigated the role of cosmological constant to obtain the proper solution for relativistic stars by used of anisotropic fluid structure of the cosmic stellar filament . There is a provoking discussion in [28] , wherein the Lemaitre- Tolman- Bondi (LTB) space-time is applied as a powerful tool in describing the cosmic gravitational collapse. Hence it was suggested to look for the possibility of a viable polynomial model of
gravity for collapsing of an inhomogenous LTB space-time.
The arbitrariness in the choice of different functional forms of
to the need for constraints. Some constraints on
gravity from the energy conditions have been considered by Santos et al. [29] . In this direction, considerable research has been carried out in
gravity to find the emergence of relativistic compact stellar objects [30] and their stability in the background of some observations. The energy density has remained positive and approached to the maximum value at the core of stars. Yousaf et al. [31] have discussed the effect of
gravity on formation of spherically symmetric wormhole under distinct distribution of matter content as perfect , barotropic and anisotropic fluid only by using the validity of WEC (weak energy condition) and NEC(null energy condition). It is worth setting up the energy conditions of field equations in order to constraint the energy-momentum tensor selection in
gravity in an anisotropic universe, which is what this paper attempts to do. In the case all the energy conditions are satisfied then we can refer to
as the true sources of energy and momentum. Santos et. al. The study of dark
energy is possible through its equation of state parameter
. If the value of
this is equivalent to the cosmological constant
. For minimally
coupled scalar field such as quintessence (
) and phantom (
).
From results of supernovae data CMBR anisotopy, the values are
and
respectively.The notion of energy conditions have been used to explore viable models of gravity in [32] . The viable solutions are explained by assuming Starobinsky model
and specific shape function.
It is to be noted that perfect fluid is one of the matter of the universe which go through variation and contraction of the action. At the back drop of our mind, we would be interested to see that
would gradually be fitted into GR equation. The present paper is arranged as follows. In Section II, we present the field equations obtained from the substitution of the Bianchi type I space-time metric in the
gravity. Some cosmological parameters of this approach are also defined. In Section III, we have constructed the cosmological models with the power law cosmology and hybrid cosmology. The physical behaviors of the models along with the energy conditions are described in Section IV. Concluding remarks are given in Section V.
2. Basic Equations
In this section, we will derive the field equations of
gravity in the framework of Bianchi type I space-time and energy momentum tensor in the form of perfect fluid with the choice of the functional
, where
is an arbitrary constant.The action in the
theory of gravity can be written as .
(1)
with
be a function of R and T, g is the determinant of the metric and
is the matter Lagrangian density. The energy-momentum tensor can be
defined as
. It is to be mentioned that here the matter
Lagrangian is assumed to be dependent only on the metric tensor
and not on its derivatives. So, the stress-energy tensor can be expressed as :
(2)
Now, varying action (1) with respect to the metric component
, the field equations of
gravity can be obtained as
(3)
with
being the Ricci tensor, where
. Here,
,
and
, with
being the covariant derivative.
We consider the homogeneous and anisotropic Bianchi universe as
(4)
where the metric potentials
and
are function of the cosmic time. The energy-momentum tensor is
(5)
Here,
and p, respectively, denote the matter-energy density and proper pressure of the matter field and
is the four-velocity vector of the fluid in a co-moving coordinate system that satisfies
. Now, by assuming the matter Lagrangian as
, the field equations Equation (3) for the choice
can be reduced to
(6)
where,
. In the co-moving coordinate system, the field equations are
(7)
(8)
(9)
An over dot on the field variables represent the ordinary derivative with respect to the cosmic time. In order to understand the model physically, we will
consider the directional Hubble rates in the form
and
,
which subsequently yield the mean Hubble parameter as
, with (a) denotes the mean scale factor of the universe. So, the set of field Equations (7)-(9) can be expressed in term of the scale factor (a) as:
(10)
(11)
(12)
where we have defined
. By performing algebraic manipulation on Equations (10)-(12), the pressure (p) and energy density (
) of the model concerning the scale factor can be expressed respectively as
(13)
and
(14)
Subsequently, with respect to the scale factor, the equation of state (EoS) parameter
is obtained as
(15)
We also define an effective cosmological constant (ECC) for the model as
(16)
3. Cosmological Features of the Model
In this section, we will discuss the cosmological features of the model with the scale factor being in the form of the power law and hybrid law cosmologies.
3.1. Power Law Cosmology
In the previous section, we have expressed the physical parameters of the model with respect to the Hubble parameter. In order to study the dynamics of the universe, we need to know the parameters with respect to the cosmic time, hence
we consider here the power law function in the form
. This type of
assumption is quite helpful to understand the background cosmology, which can also be checked in the solutions for the standard cosmology Friedman equations. In order to employ some amount of anisotropic, we have assumed here
, where (k) is a constant. It can be seen that when
, the space time reduces to standard FRW model. With these assumptions, we obtain the
directional and mean Hubble parameters, respectively, as
,
and
. So, the pressure and energy density can be obtained as
(17)
(18)
The graphical behavior of the pressure and energy density have been represented in Figure 1 and Figure 2 respectively. The pressure remains entirely in the negative domain. It starts increasing at initial phase and after a specific time period, it remains constant; however it never touches the horizontal axis.The behavior of the variation of pressure with respect to cosmic time are
Figure 1. Plot of pressure (p) versus cosmic time for power law.
Figure 2. Plot of energy density (ρ) versus cosmic time for power law.
shown for
and various range of m: 1.42, 1.5, 1.6, 1.8, 2 and μ: 0.01, 1.15, 1.8, 0.07, 0.001 as anisotropy parameter and scaling constant where have shown with Red, Black, Orange, Green and Blue respectively. Pressure minimized on
and
.
Subsequently the EoS parameter and ECC can be obtained as
(19)
(20)
In Figure 3 and Figure 4, we have respectively graphically represented the EoS parameter and ECC. As it is provided in Figure 3, the EOS parameter is considered and represented as increasing function of time whereas the primary points of evolutionary process starts from 0.08 while it incrementally shifting to −0.6 while it remains on negative domain. It can be seen that the EoS parameter strongly depends on the anisotropic parameter k.
In Figure 4, ECC depend on the cosmic time (t). A function which decreased by time and approaches to its small positive value at present epoch. In the early universe, (μ) takes a very large value and after that decreasing as time increases. Therefore it explained that the energy density is predominance as a decreasing functional behavior in positive domain permanently.
3.2. Hybrid Scale Factor
The power-law and exponential law cosmologies are having constant deceleration parameter and only can be used to describe the epoch based evolution of the universe. These cosmologies do not exhibit the transition of the universe from deceleration to acceleration. So, in this section, in order to show the early deceleration and late time acceleration, we have considered a scale
factor in the form
known as hybrid scale factor [19] [21] [33] [34] .
Figure 3. Plot of EoS parameter (ω) versus k for power law.
Figure 4. Plot of ECC (Λ) versus cosmic time for power law.
This scale factor is having two components: one component behavior is like exponential universe and the other is the power law expansion. Moreover, in this scale factor the exponential component dominates at late phase while in the early phase of cosmic evolution the power law component dominates. It can be noted that when
, the scale factor reduces to a power law expansion whereas when for
, only the exponential law can be recovered. For this
scale factor, the Hubble rate is
,
,
. The pressure and energy density for the hybrid scale factor can be respectively obtained as:
(21)
(22)
In Figure 5 and Figure 6 we observe the fluid pressure as a function of time and variation of energy density. In the Figure 5 we assume representative value of the anisotropy i.e.
, for
In order to minimizing the pressure it is required to inequality of
Figure 5. Plot of pressure (p) versus cosmic time for hybrid scale factor.
Figure 6. Plot of energy density (ρ) versus cosmic time for hybrid scale factor.
the pressure raised up to its lowest negative value.
From Figure 6 we observe that the energy density always remain in the positive domain and it has decreased to its small value at late times. At the same time the energy density start decreasing from a high value and decreases with increase in the time period. On the other hand, the pressure raised up to its lowest negative value from some large negative value (≈0.023) and gradually vanishing to
at late times. We have obtained the EoS parameter ω and effective cosmological constant Λ can be obtained as
(23)
(24)
Figure 7 represents evolution of EoS with respect to cosmic time. In the present study for hybrid case ω remains in the negative domain favoring a quintessence phase. in the present epoch ω get the lowest value of −0.32 corresponding to
.
Figure 7. Plot of EoS parameter (ω) versus cosmic time for hybrid scale factor.
4. Physical Behavior of the Models
The scalar expansion and deceleration parameter of the model can be obtained
as
and
. We can infer that the scalar expansion
vanishes at late time and deceleration parameter remains constant at −1. It can be noted that when
, the model decelerates and accelerates for
.
The volume scale factor V and the ratio of anisotropy
can be respectively calculated as
and
. As usual for
, the model
reduces to an isotropic universe. The state finder diagnostic pair which give some insight to the geometrical part of the model. The pair for the power law
cosmology can be obtained as
,
.
For hybrid cosmology
and
. Figure 8 has
shown evolution of deceleration parameter in terms of cosmic time for transition from deceleration to acceleration phase. In this model we have
positive values of q in the range of
for decelerated universe and
facing to negative values of q for accelerated universe at the late time for
. Similarly for a hybrid cosmology, the pair can be obtained as
, and
Energy conditions as a sets of linear equations put some additional constraints on the model. However Equation (25) has demonstrated that the field equation has been physically approved. Energy conditions derived as:
Figure 8. Plot of Deceleration parameter (ω) versus cosmic time for hybrid scale factor.
(25)
Now, the energy conditions for power law cosmology can be obtained as:
(26)
(27)
(28)
In Figure 9, corresponding graphical behavior of energy conditions for power law cosmology were plotted for
and various values of m. In the present model all the energy conditions have been satisfied. The physical trend of all cases have remained almost the same for various choice of scaling constant μ.
Similarly, for a hybrid cosmology, the energy conditions for the model can be calculated
(29)
(30)
(31)
Figure 10 shows the graphical behavior of energy conditions for hybrid scale factor for
. It is implied by the Green(U), Blue(Q) and Red(G) energy frame which is expressed the dynamical trend for SEC, DEC, WEC. It is clearly shown that all energy conditions have been satisfied. Indeed it is assumed that typically SEC is the only known energy condition which we need to violate in order to find a bounce. An outstanding feature of this models is energy condition satisfaction for some values of the input parameters.
Another notable cases is evaluating the effectiveness of higher curvature on the dynamical variables of spherically compact stars which filled by an imperfect fluid as matter content has been investigated on [35] . One of the well-consistence form of
is
which is mostly coped with observational results. The selected non-linear complex model considered as
which has been compared with experiment result s of a neutron star named as 4U1820-30. Moreover the correspondence Raychaudhuri equations for shear and expansion scalar equations has established for shear and expansion scalar.
Figure 9. Plot of energy conditions versus cosmic time for power law cosmology.
Figure 10. Plot of energy conditions versus cosmic time for hybrid cosmology.
5. Conclusion
The work done in this paper is to investigate the cosmological model obtained in
gravity in Bianchi type I space-time with the functional
. The dynamical parameters are calculated and analyzed with a new mathematical formalism in power law cosmology and hybrid scale factor cosmology. We have shown that there exist various values for inputs parameters which the four energy conditions may be satisfied simultaneously. The state finder diagnostic pair and energy conditions are also studied. The results obtained in both the models are in accordance with the observational findings of accelerated expansion of the universe. Also the result of this work could be extended to a new approach of extended theories of gravity such as
that could lead to unveil the secret of dark energy in modern cosmology. Hence we addressed to some successful studies on this formalism as [36] which lightened the present of a higher derivatives model of
by using of energy conditions and evaluate their viability by making some constraints on coefficient of Ricci scalar and energy-momentum tensor. Also for finding more discussions on implication of energy conditions on formation of universe, it is good to present about the Gauss Bonnet invariant model of extended theory of gravity
in [37] which tries to find the constraint on cosmological parameters which essentially is needed for viability of the constructed model.
Acknowledgements
I would like to thank Professor.Bindu. A. Bambah for her expert advice and encouragement throughout this research paper. The author is thankful to the anonymous referees for constructive comments concerning our manuscript suggestions for the improvement of the paper.