In a large scale structure, the universe is highly homogeneous and isotropic, shown in . These properties are also supported by the work of Benet et al.  and Spergel et al.  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  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 (where, , 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 has been introduced by Ananda and Bruni  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  used inhomogeneous Hubble parameter term to discuss the variation of dark energy universe relative to different equation of state. In 2006, Capozzelio  pointed out that the observational constraint impact along with quadratic EoS on dark energy. Nojiri et al.  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.  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. , Firoze et al.  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  found a new exact solution of Einstein-Maxwell equation in terms of quadratic EoS. Unification of vacuum, dark energy and radiation  was also brought under active consideration through a cosmological model based on quadratic EoS.
Later Chavanis  considered early inflation, intermediate decelerating expansion, late accelerating expansion as quadratic EoS and developed a model. In 2013, Sharma and Ratanpal’s work  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 , Malaver  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.
2. Metric and Field Equations
The Bianchi type-I line element is given by
where A, B and C are scale factors and are functions of time t only.
The Einstein field equations, in natural limits ( and ), are
where is the Ricci tensor, R is the Ricci scalar and is the energy-momentum tensor.
The energy-momentum tensor for the perfect fluid is given by
where is the energy density, p is the pressure and is the four velocity vector satisfying .
Here, we have assumed an equation of state (EoS) in the general form for the matter distribution.
We have considered it in the quadratic form as follows
where , and are the constant and strictly 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:
where overhead dot ( ) denotes differentiation with respect to time t.
The vanishing divergence of Einstein tensor which leads to , yields the following energy conservation equation
We have defined the spatial volume V and average scale factor “a” for Bianchi type-I space-time as
The mean Hubble parameter H for Bianchi type-I universe is defined as
where , , are the directional Hubble parameters in the directions of x, y and z axes respectively.
3. Solution of the Field Equations
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 of the shear tensor in the model that leads to
Subtracting (7) from (8), we have
From equations (12) and (13), we get
which can be reduced to the following form
Integrating Equation (15), we obtain
Equation (16) can be written as follows
which is a first order linear differential equation in and t. Hence the general solution of Equation (17) is found as
where and 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,
where n is real number. Then (18) yields
Choosing without loss of generality and hence we obtain
From Equation (12), (19) and (21), we get
The metric of the model can be written in the form
Using the directional scale factors for the model (23) in (10), we compute the spatial volume V as follows
For the spatial volume to be positive, we must have . It is seen that the spatial volume is zero at .
Choosing and in Equation (4) as studied in    and from Equations (5) and (6) by using Equations (19), (21) and (22), we obtain the energy density for the model (23) as
From (4) and (25), we obtain the pressure as
Again, choosing and in Equation (4) as in  and from Equations (5) and (6) by using Equations (19), (21) and (22), we obtain the energy density for the model (23) as
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 and deceleration parameter q which are defined and found to be
4. Results and Discussions
It has been observed from Equation (25) and Figure 1 that the energy density is a decreasing function of time and while diverges at the initial epoch of the cosmological evolution for and . Also, from Equation (27) and Figure 1, it is clearly seen that the energy density is a decreasing function of time and while diverges at the initial epoch of the cosmological evolution for and . 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    . In an agreement with the energy dominated universe, the model, with and , shows a negative pressure corresponding to Figure 2 and Equation (26).
Figure 1. Density vs Cosmic Time using Equations (25) and (27) ( ).
Figure 2. Pressure vs Cosmic Time using Equations (26) and (28) ( ).
Again, Figure 2, corresponding to Equation (28) for and , depicts the evolution of pressure p and shows that p is positive and when 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 . 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 when . The expansion scalar and shear scalar start with infinitely large value at and tend to zero as 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 , but for , it tends to infinity. Also, since except at , 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   as well as CMBR   are in favour of accelerating models ( ), but both do not entirely rule out the decelerating ones which are also realistic with these observations supported by the work of Vishwakarma . 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.
 Bennett, C.L., et al. (2003) First-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results. Astrophysical Journal Supplement Series, 148, 1-27.
 Spergel, D.N., et al. (2003) First-Year Wilkinson Microwave Anisotropy Probe (WMAP)* Observations: Determination of Cosmological Parameters. The Astrophysical Journal Supplement Series, 148, 175-194.
 Ananda, K.N. and Bruni, M. (2006) Cosmo-Dynamics and Dark Energy with a Quadratic EoS: Anisotropic Models, Large-Scale Perturbations and Cosmological Singularities. Physical Review D, 74, Article ID: 023523.
 Nojiri, S. and Odintsov, S.D. (2005) Inhomogeneous Equation of State of the Universe: Phantom Era, Future Singularity, and Crossing the Phantom Barrier. Physical Review D, 72, Article ID: 023003.
 Capozziello, S., Cardone, V.F., Elizalde, E., Nojiri, S. and Odintsov, S.D. (2006) Observational Constraints on Dark Energy with Generalized Equations of State. Physical Review D, 73, Article ID: 043512.
 Rahaman, F., Jamil, M. and Chakraborty, K. (2010) Revisiting the Classical Electron Model in General Relativity. Astrophysics and Space Science, 331, 191-197.
 Chavanis, P.H. (2013) A Cosmological Model Based on a Quadratic Equation of State Unifying Vacuum Energy, Radiation, and Dark Energy. Journal of Gravity, 2013, Article ID: 682451.
 Chavanis, P.H. (2015) A Cosmological Model Describing the Early Inflation, the Intermediate Decelerating Expansion, and the Late Accelerating Expansion of the Universe by a Quadratic Equation of State. Universe, 1, 357.
 Sharma, R. and Ratanpal, B.S. (2013) Relativistic Stellar Model Admitting a Quadratic Equation of State. International Journal of Modern Physics D, 22, Article ID: 1350074.
 Ahmed, N., et al. (2019) Transition from Decelerating to Accelerating Universe with Quadratic Equation of State in f(R,T) Gravity. NRIAG Journal of Astronomy and Geophysics, 8, 198-203.
 Reddy, D.R.K., Adhav, K.S. and Purandare, M.A. (2015) Bianchi Type-I Cosmological Model with Quadratic Equation of State. Astrophysics and Space Science, 357, 20.
 Riess, A.G., et al. (1998) Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. The Astronomical Journal, 116, 1009-1038.
 Vishwakarma, R.G. (2000) A Study of Angular Size-Redshift Relation for Models in Which Λ Decays as the Energy Density. Classical and Quantum Gravity, 17, 3833.