Advances in High Energy Physics

Volume 2016 (2016), Article ID 8543560, 8 pages

http://dx.doi.org/10.1155/2016/8543560

## Dynamics of an Anisotropic Universe in Theory

^{1}Department of Mathematics, Birla Institute of Technology and Science-Pilani, Hyderabad Campus, Hyderabad 500078, India^{2}Department of Physics, Indira Gandhi Institute of Technology, Sarang, Dhenkanal, Odisha 759146, India

Received 4 June 2016; Revised 15 August 2016; Accepted 25 August 2016

Academic Editor: Chao-Qiang Geng

Copyright © 2016 B. Mishra et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

Dynamics of an anisotropic universe is studied in gravity using a rescaled functional , where is the Ricci Scalar and is the trace of energy-momentum tensor. Three models have been constructed assuming a power law expansion of the universe. Physical features of the models are discussed. The model parameters are constrained from a dimensional analysis. It is found from the work that the anisotropic Bianchi type () model in the modified gravity generally favours a quintessence phase when the parameter is either or . We may not get viable models in conformity with the present day observation for .

#### 1. Introduction

In the past decade the idea of modifying gravity on cosmological scale has attracted a lot of attention. It gains momentum through the theoretical developments involving higher dimensional theories and in constructing renormalizable theories of gravity. Modified gravity represents an intriguing possibility for resolving the theoretical challenge posed by late time acceleration. However, it turns out to be extremely difficult to modify general relativity at low energy regime without violating observational constraints.

Observations from type Ia supernova [1–3], large scale structures [4–6], and cosmic microwave background [7, 8] confirm that the cosmic expansion is accelerating. The very cause behind this acceleration may be a mysterious energy source called dark energy which accounts for two-thirds of the total energy budget of the universe.

The role of modified theories to understand the mechanism behind the cosmic speed-up is more promising. The objective behind such theories is to explain the puzzling late time cosmic dynamics without the need of dark energy components in the field equations. Several modified theories of gravity have been proposed in recent times. Of these a few models such as gravity [9–11], gravity [12–16], and gravity [17–19] are studied widely. In theory, a more general function of the Ricci Scalar is used in place of whereas gravity is a generalized version of teleparallel gravity. Motivated by the success of cosmological constant as a simple and good candidate of dark energy, some matter field is also coupled with the function of in the geometry side of the action in some modified theories of gravity () theory).

Along the line of interest of incorporating some matter components in the action geometry, theory has been proposed by Harko et al. [20] which, of late, has been an interesting framework to investigate accelerating models. Moreover, the reconstruction of arbitrary FRW cosmologies is possible by an appropriate choice of the functional . Many authors have investigated the astrophysical and cosmological implications of the gravity [21–24]. Jamil et al. [25] have reconstructed some cosmological models for some specific forms of in this modified gravity. Shamir et al. [26] obtained exact solution of anisotropic Bianchi type I and type V cosmological models whereas Chaubey and Shukla [27] have obtained a new class of Bianchi cosmological models using special law of variation of parameter. Using a decoupled form of , that is, for Bianchi type universe, Ahmed and Pradhan [28] have studied the energy conditions of perfect fluid cosmological models and Yadav [29] obtained some string solutions.

Singh and Kumar [30] have investigated the effect of bulk viscosity in theory and suggested that inclusion of dissipative energy sources like bulk viscosity may be able to explain the early and late time accelerations of the universe. Mishra and Sahoo [31] have studied Bianchi type cosmological models assuming . In that work, Mishra and Sahoo have obtained exact solutions to the modified field equations by assuming a specific anisotropic relation among the metric potential. The models constructed in [31] are obtained to be accelerating for the model parameter . Samanta [32] has obtained exact solution of gravitational field equations in Kantowski-Sachs space time and Shamir [33] has constructed some cosmological models in Bianchi type V space time. In the framework of this modified gravity, recently, Mishra et al. [34] have presented the Einstein-Rosen nonstatic cosmological model with quadratic form of gravity. Sahoo et al. [35] have investigated the background cosmology of power law and exponential law of volumetric expansion in gravity in Kaluza-Klein model. Recently Yousaf et al. have investigated the irregularity energy density factor in theory responsible to disturb the stability of homogeneous universe [36] and also explored the evolutionary behaviour of compact objects [37].

In this work, we have constructed some cosmological models in theory for general Bianchi type () universe. The present work differs from the work of Mishra and Sahoo [31] in the sense that here we have used a rescaled functional instead of the one used in [31]. We have reconstructed some viable accelerating models in conformity with the recent observations from some physical basis. Our work is organised as follows. In Section 2, the basic formalism of the theory has been presented. The dynamics of anisotropic model has been presented in a general form. In Section 3, some specific models have been constructed to study the dynamics. We presented the conclusions of the works at the end in Section 4.

#### 2. Basic Formalism

The modified four-dimensional Einstein-Hilbert action in gravity theory with a specific choice of matter Lagrangian can be considered asHere, is an arbitrary function of and . Following the work of Harko et al. [20], we have considered the matter Lagrangian as , where is the pressure of the perfect cosmic fluid. The trace of the energy-momentum tensor is obtained from . is the rest energy density and is the four velocity vectors.

The field equations can be obtained from the action in (1) for the choice of the functional as [20] and are the partial differentiation of the respective functional with respect to their arguments. The functional can be chosen arbitrarily to get viable cosmological models. Here we chose the functional and to be linear in their arguments, and , so that . is a constant scaling factor.

Equation (2) now reduces toEquation (3) can now be recast as the usual Einstein Field Equation where can be identified with the cosmological constant that evolves with cosmic time. The role of the particular choice of the functional can be well understood in getting such a simplified approach to a modified theory which does not contain any dark sector component but can be instrumental in providing viable accelerating models.

For a spatially homogeneous and anisotropic Bianchi type () space time considered in the form(3) can be explicitly written aswhere the metric potentials , , and are functions of cosmic time . The constant exponent decides the behaviour of the model and can take integral values such as , and 1. An overdot on a field variable denotes differentiation with respect to time . We define the directional Hubble parameters along different directions as , , and . The mean Hubble parameter becomes . The field equations (5) can now be expressed aswhere .

The pressure and rest energy density can be obtained from (8) and (9) in general forms aswhere and .

From (11) and (12), we obtain the equation of state parameter and the effective cosmological constant as

Equations (11)–(13) provide the dynamical behaviour of the universe. However, the dynamics can only be assessed if the behaviour of these properties is known in terms of the directional Hubble rates for a given value of the exponent . In other words, the formalism as described above can help us to study a background cosmology for an assumed dynamics of the universe. In this context, we can consider the power law cosmology, where the cosmic expansion is governed through a volume scale factor of the form . is an arbitrary positive constant usually determined from the background cosmology. Power law cosmology has been widely studied in recent times because of its functional simplicity and ability to provide a first-hand information about the dynamics of the universe. For such an assumption, the radius scale factor can be . The deceleration parameter for power law expansion of the universe is a constant quantity: , which can be negative for and positive for . It is worth mentioning here that a positive describes a decelerating universe whereas a negative describes an accelerating universe. In order to keep a pace with the recent observational data favouring an accelerating universe, should be greater than 3.

Some other kinematical parameters of the universe are the scalar expansion , shear scalar , and the average anisotropy parameter defined, respectively, aswhere with . is a measure of deviation from isotropic expansion. One can get the isotropic behaviour of the model for . For the power law cosmology with , the scalar expansion becomes and consequently the mean Hubble rate becomes .

Geometrical analysis of dark energy models is usually performed through the statefinder pair and given by and . In the present model, we obtain these parameters as and . These parameters are constants of cosmic time and depend only on the exponent . It can be emphasized here that exact determination of these parameters from different observational bounds can constrain the exponent in narrow ranges.

#### 3. Dynamics of Anisotropic Universe

Once the cosmic expansion behaviour is known, it becomes simpler to study the background cosmology of the diagonal universe. However, the exponent in the metric can assume integral values such as , , and . Each value of corresponds to a different cosmological model with different dynamical behaviour. In view of this, in the following, we discuss the dynamical features of the three possible models in the framework of theory.

##### 3.1. Model-I

A substitution of in (10) yields , where the integration constant has been rescaled to unity. Assuming an anisotropic relationship , we can write the functionals and aswhere is an arbitrary positive constant. For a power law cosmology, we have , . Consequently the directional scale factors are and . Thus we can haveThe dynamical behaviour of the model is decided from the behaviour of the equation of state parameter and the effective cosmological constant . However, these two parameters depend on the time varying nature of functionals and which in turn depend on the parameters and . If , the terms within the square brackets in (16) dominate at late times of cosmic evolution whereas the terms containing dominate at early phase of cosmic evolution. Here we wish to adopt a dimensional analysis to get some idea into the general behaviour of these functionals. Since and are two dimensionless constants, it appears that the dimensionality of the time dependent factors for a given functional should remain the same. In other words, we can have so that and become

The equation of state parameter and the effective time varying cosmological constant are obtained from (13) as

It is obvious from the above expressions that the equation of state parameter is a constant quantity for a given value of scaling constant and the anisotropic parameter . However, the effective cosmological constant decreases quadratically with cosmic time. In order to get viable cosmological models in conformity with recent observations, the cosmological constant should be dynamically varying from large positive values at an initial epoch to vanishingly null values at late times of cosmic evolution. Similarly, the equation of state parameter should be negative with values less than at late times. This behaviour will enable us to constrain the parameter . In Figure 1, we have shown the variation of as a function of . Here we chose a negative value of the scaling constant ; that is, . Accordingly the model parameter is decided from the relation . The reason behind this particular choice is to look at the present problem from the backdrop of general relativity where the modified gravity field equation (3) appears as the Einstein Field Equation with a time varying effective cosmological constant . It can be inferred from the figure that the equation of state parameter increases almost linearly from a negative value for lower to zero at higher . One can note that the present model will collapse at and therefore we restrict the values of below 1. For , remains in the quintessence region. One can note that the results may be sensitive to the choice of the parameter . In Figure 2, we have shown the dynamical variation of the effective cosmological constant for some representative values of . As is required for an explanation to the late time cosmic acceleration, varies from large positive values in the beginning to vanishingly small values at late times. In Figure 2(a), the evolution of is shown as a function of cosmic time where as in Figure 2(b) its evolution is shown as function of redshift with the scale factor at present epoch being unity. As it appears from the figure, the evolutionary behaviour of is affected by the choice of the values of . At remote past, curves of with low values of remain below the curves with higher values of . However, at certain point of cosmic time corresponding to a redshift of , there occurs a reversal of the behaviour; that is, curves of with higher values of remain below the curves with lower values of . The choice of the model parameter does not affect the general time varying trend of the effective cosmological constant.