In a large scale structure, the universe is highly homogeneous and isotropic, shown in [1]. These properties are also supported by the work of Benet et al. [2] and Spergel et al. [3] via the cosmic microwave background radiation. Aside homogeneity property, isotropy throws a challenge that anisotropic models cannot be isotropized appropriately because of their time evolvement in future. However this is to note that inflation is a suitable method to solve this problem.

Ananda and Bruni’s [4] works in relativistic dynamics study of Robertson Walker models using a non-linear quadratic equation of state (EoS) are excellent to enhance to open new door of research. The authors considered an equation of state as of the form p = p 0 + α ρ + β ρ 2 (where, p 0 , α and β are parameters), in particular the quadratic term in the above equation and revealed that anisotropy at the singularity existing in brane world scenario can be reproduced in the context of general relativity.

Depending on the appearance or look of the different forms of quadratic EoS, there are extensive uses of it in many contemporary research areas. For instances, EoS of the form p = α ρ + ( ρ 2 / ρ c ) has been introduced by Ananda and Bruni [5] to study the effect of quadratic EoS in homogeneous and inhomogeneous cosmological models to isotropize the universe at early age when initial singularity is continued.

In a firm study, Nojiri and Odintsov [6] used inhomogeneous Hubble parameter term to discuss the variation of dark energy universe relative to different equation of state. In 2006, Capozzelio [7] pointed out that the observational constraint impact along with quadratic EoS on dark energy. Nojiri et al. [6] works reveal that quadratic EoS may be nice tool to describe dark energy or unified dark matter in a more generous way. In 2009, Rahaman et al. [8] showed how electron can be modeled precisely instead of spherically symmetric charged perfect fluid distribution of matter characterized by quadratic EoS. Followed by the works of Rahaman et al. [8], Firoze et al. [9] introduced “Charged anisotropic matter” into quadratic EoS where the approach results in a new class of static spherically symmetric models of relativistic star. Einstein and Maxwell equations were brought under the application of quadratic EoS. For charged anisotropic matter distribution, Maharaj and Takisa [10] found a new exact solution of Einstein-Maxwell equation in terms of quadratic EoS. Unification of vacuum, dark energy and radiation [11] was also brought under active consideration through a cosmological model based on quadratic EoS.

Later Chavanis [12] considered early inflation, intermediate decelerating expansion, late accelerating expansion as quadratic EoS and developed a model. In 2013, Sharma and Ratanpal’s work [13] results in a favorable solution focusing the interior of a static symmetric spherically asymmetric compact anisotropic star and showed that the model admits EoS in quadratic form. After Maharaj and Tasika [10], Malaver [14] solved Maxwell in Maxwell system of equations using some elementary function throughout the behavior of compact anisotropic relativistic objects in quadratic equation of state.

It reveals from the above discussion that quadratic EoS plays an important and interesting role in many areas including dark energy and dynamics of different models in general relativity. In the study of anisotropic problem, in particular, Bianchi type-I homogeneous cosmological model containing perfect fluid is an interesting idea that lead us to continue research using quadratic EoS to generate an expanding and decelerating universe.

We have organized this paper as follows: In Section 2, we have applied Bianchi type-I metric to obtain field equations. Explicit solution of the equations is presented in Section 3 followed by relevant interpretations of the results in Section 4. We summarize this paper by giving a conclusion in Section 5.

The Bianchi type-I line element is given by

d s 2 = d t 2 − A 2 d x 2 − B 2 d y 2 − C 2 d z 2 (1)

where A, B and C are scale factors and are functions of time t only.

The Einstein field equations, in natural limits ( 8 π G = 1 and c = 1 ), are

R i j − 1 2 R g i j = − T i j (2)

where R i j is the Ricci tensor, R is the Ricci scalar and T i j is the energy-momentum tensor.

The energy-momentum tensor T i j for the perfect fluid is given by

T i j = ( ρ + p ) u i u j − p g i j (3)

where ρ is the energy density, p is the pressure and u i is the four velocity vector satisfying g i j u i u j = 1 .

Here, we have assumed an equation of state (EoS) in the general form p = p ( ρ ) for the matter distribution.

We have considered it in the quadratic form as follows

p = α ρ 2 + β ρ + γ (4)

where α , β and γ are the constant and strictly α ≠ 0 which preserves the quadratic nature of the equation of state.

In co-moving coordinate system, the Einstein field Equations (2) for the metric (1) with the help of Equation (3) reduce to following set of equations:

A ˙ B ˙ A B + B ˙ C ˙ B C + A ˙ C ˙ A C = ρ (5)

B ¨ B + C ¨ C + B ˙ C ˙ B C = − p (6)

A ¨ A + C ¨ C + A ˙ C ˙ A C = − p (7)

A ¨ A + B ¨ B + A ˙ B ˙ B C = − p (8)

where overhead dot ( ⋅ ) denotes differentiation with respect to time t.

The vanishing divergence of Einstein tensor which leads to T ; j i j = 0 , yields the following energy conservation equation

ρ ˙ + ( A ˙ A + B ˙ B + C ˙ C ) ( ρ + p ) = 0 (9)

We have defined the spatial volume V and average scale factor “a” for Bianchi type-I space-time as

V = a 3 = A B C (10)

The mean Hubble parameter H for Bianchi type-I universe is defined as

H = a ˙ a = 1 3 ( H x + H y + H z ) (11)

where H x = A ˙ A , H y = B ˙ B , H z = C ˙ C are the directional Hubble parameters in the directions of x, y and z axes respectively.

The field Equation (5)-(8) are four independent equations in five unknowns A, B, C, ρ and p. For the exact explicit solutions to these Einstein’s field equations, we have assumed the expansion scalar ( θ ) as proportional to the component σ 1 1 of the shear tensor σ i j in the model that leads to

A = ( B C ) m (12)

Subtracting (7) from (8), we have

B ¨ B − C ¨ C + A ˙ A ( B ˙ B − C ˙ C ) = 0 (13)

From equations (12) and (13), we get

B ¨ B − C ¨ C + m ( B ˙ B − C ˙ C ) ( B ˙ B + C ˙ C ) = 0 (14)

which can be reduced to the following form

d d t ( B ˙ B − C ˙ C ) B ˙ B − C ˙ C + ( m + 1 ) ( B ˙ B + C ˙ C ) = 0 (15)

Integrating Equation (15), we obtain

( B ˙ B − C ˙ C ) ( B C ) m + 1 = l 1 (16)

Equation (16) can be written as follows

d d t ( B m + 1 ) − ( m + 1 ) C ˙ C B m + 1 = ( m + 1 ) l 1 C m + 1 (17)

which is a first order linear differential equation in B m + 1 and t. Hence the general solution of Equation (17) is found as

B m + 1 = ( m + 1 ) C m + 1 ∫ l 1 C 2 ( m + 1 ) d t + l 2 C m + 1 (18)

where l 1 and l 2 are integrating constants.

It is to be noted that in this case the solution of Einstein’s field equations reduces to the integration of (18) if the explicit form of the scale factor C is known.

We obtain particular solution of (18) for a simple choice of the function C. For this, we choose,

C = t n (19)

where n is real number. Then (18) yields

B m + 1 = ( m + 1 ) l 1 1 − 2 ( m + 1 ) n t 1 − ( m + 1 ) n + l 2 t ( m + 1 ) n (20)

Choosing l 2 = 0 without loss of generality and hence we obtain

B = { ( m + 1 ) l 1 1 − 2 ( m + 1 ) n } 1 m + 1 t 1 m + 1 − n (21)

From Equation (12), (19) and (21), we get

A = { ( m + 1 ) l 1 1 − 2 ( m + 1 ) n } m m + 1 t m m + 1 (22)

The metric of the model can be written in the form

d s 2 = d t 2 − { ( m + 1 ) l 1 1 − 2 ( m + 1 ) n } 2 m m + 1 t 2 m m + 1 d x 2 − { ( m + 1 ) l 1 1 − 2 ( m + 1 ) n } 2 m + 1 t 2 m + 1 − 2 n d y 2 − t 2 n d z 2 (23)

Using the directional scale factors for the model (23) in (10), we compute the spatial volume V as follows

V = ( m + 1 ) l 1 1 − 2 ( m + 1 ) n t (24)

For the spatial volume to be positive, we must have ( m + 1 ) l 1 1 − 2 ( m + 1 ) n t > 0 . It is seen that the spatial volume is zero at t = 0 .

Choosing β = 1 and γ = 0 in Equation (4) as studied in [15] [16] [17] and from Equations (5) and (6) by using Equations (19), (21) and (22), we obtain the energy density for the model (23) as

ρ = 1 α { − n + 1 m + 1 − ( m − n ) 2 } 1 t (25)

From (4) and (25), we obtain the pressure as

p = { − n + 1 m + 1 − ( m − n ) 2 } 1 t 2 + 1 α { − n + 1 m + 1 − ( m − n ) 2 } 1 t (26)

Again, choosing β = − 1 and γ = 0 in Equation (4) as in [18] and from Equations (5) and (6) by using Equations (19), (21) and (22), we obtain the energy density for the model (23) as

ρ = 1 α { m − 1 ( m + 1 ) 2 + m n − m 2 + 2 m + 1 m + 1 − m ( m + 1 ) − n ( n − 1 ) } 1 t (27)

p = { m − 1 ( m + 1 ) 2 + m n − m 2 + 2 m + 1 m + 1 − m ( m + 1 ) − n ( n − 1 ) } 1 t 2     − 1 α { m − 1 ( m + 1 ) 2 + m n − m 2 + 2 m + 1 m + 1 − m ( m + 1 ) − n ( n − 1 ) } 1 t (28)

Some other physically important parameters for the model (23) such as the mean Hubble parameter H, the scalar expansion θ , the average anisotropy parameter Δ , the shear scalar σ 2 and deceleration parameter q which are defined and found to be

H = 1 3 ( A ˙ A + B ˙ B + C ˙ C ) = 1 3   t (29)

θ = 3 H = 1 t (30)

Δ = 1 3 ∑ i = 1 3 ( Δ H i H ) 2 = 2 ( m 2 − m + 1 ) + 6 n ( m + 1 ) ( m n + n − 1 ) ( m + 1 ) 2 (31)

σ 2 = 3 2 Δ H 2 = ( m 2 − m + 1 ) + 3 n ( m + 1 ) ( m n + n − 1 ) 3 ( m + 1 ) 2 t 2 (32)

q = d d t ( 1 H ) − 1 = 2 (33)

It has been observed from Equation (25) and Figure 1 that the energy density ρ is a decreasing function of time and ρ > 0 while diverges at the initial epoch of the cosmological evolution for α = − 1 and β = 1 . Also, from Equation (27) and Figure 1, it is clearly seen that the energy density ρ is a decreasing function of time and ρ > 0 while diverges at the initial epoch of the cosmological evolution for α = 1 and β = − 1 . Moreover, we have investigated through using Matlab tools that it is obvious to require a fine tuning between α and β to obtain physically viable cosmological model [15] [16] [17] [18]. In an agreement with the energy dominated universe, the model, with α = − 1 and β = 1 , shows a negative pressure corresponding to Figure 2 and Equation (26).

Again, Figure 2, corresponding to Equation (28) for α = 1 and β = − 1 , depicts the evolution of pressure p and shows that p is positive and p → 0 when t → ∞ as in matter dominated universe.

From Equation (24), it has been observed that initially, the spatial volume V vanishes and linearly increases as time t increases and becomes infinitely large for t → ∞ . This shows that the universe starts expanding with initial zero volume and keeps expanding with increase of cosmic time t. From Equation (29), it has been observed that Hubble parameter is a decreasing function of time and H → 0 when t → ∞ . The expansion scalar and shear scalar start with infinitely large value at t = 0 and tend to zero as t → ∞ corresponding to the Equations (30) and (32) which confirms the expansion of the universe. The mean anisotropy parameter, as in Equation (31), is independent of time and is uniform throughout the whole expansion of the universe when m ≠ − 1 , but for m = − 1 , it tends to infinity. Also, since ( σ 2 / θ 2 ) = constant except at m ≠ − 1 , the model does not approach to isotropy.

The deceleration parameter (q) is an important parameter to determine whether the model inflates or not. The negative value of q confirms inflation whereas the positive value of q corresponds to standard deceleration. Our model shows decelerating behavior as the deceleration parameter is a positive constant. The remarkable thing is that though the recent observations of SNeIa [19] [20] as well as CMBR [2] [3] are in favour of accelerating models ( q < 0 ), but both do not entirely rule out the decelerating ones which are also realistic with these observations supported by the work of Vishwakarma [21]. Moreover, the study of the decelerating model of the universe is important to make a bridge between early inflation and late time acceleration.

We investigated spatially homogeneous and anisotropic Bianchi type-I cosmological models with quadratic equation of state (EoS) in the context of general relativity. The important physical and kinematical parameters in the discussion of cosmological model have been determined and studied. The limitation in the choices of the parameters of quadratic equation of state has been explored and discussed. We also showed, the universe will not be shear free and will remain anisotropic throughout the cosmic evolution with standard deceleration.

The authors declare no conflicts of interest regarding the publication of this paper.

Alam, M.M., Ullah, M.A. and Chowdhury, SMEK. (2020) Anisotropic Standard Decelerating Cosmological Model with Quadratic Equation of State. Journal of Applied Mathematics and Physics, 8, 1990-1998. https://doi.org/10.4236/jamp.2020.89148