Abstract

We use one of the simplest forms of the K-essence theory and we apply it to the classical anisotropic Bianchi type I cosmological model, with a barotropic perfect fluid () modeling the usual matter content and with cosmological constant . Classical exact solutions for any and are found in closed form, whereas solutions for are found for particular values in the barotropic parameter. We present the possible isotropization of the cosmological model Bianchi I using the ratio between the anisotropic parameters and the volume of the universe. We also include a qualitative analysis of the analog of the Friedmann equation.

1. Introduction

In recent times, some attempts to unify the description of dark matter, dark energy, and inflation, by means of a scalar field with nonstandard kinetic term, have been conducted [16]. The K-essence theory is based on the idea of a dynamical attractor solution which causes it to act as a cosmological constant only at the onset of matter domination. Consequently, K-essence overtakes the matter density and induces cosmic acceleration at about the present epoch. Usually K-essence models are restricted to the Lagrangian density of the form [4, 710] where the canonical kinetic energy is given by . K-essence was originally proposed as a model for inflation and then as a model for dark energy, along with explorations of unifying dark energy and dark matter [7, 11, 12]. Other motivations to consider the action above originate from string theory [13, 14]. For more details for K-essence applied to dark energy you can see [15] and references therein.

In this framework, gravitational and matter variables have been reduced to a finite number of degrees of freedom. For homogenous cosmological models the metric depends only on time and gives a model with a finite dimensional configuration space, called minisuperspace. In this work, we use this formulation to obtain classical solutions to the anisotropic Bianchi type I cosmological model with a perfect fluid. This class of models was considered initially in this formalism by Chimento and Forte [16]. The first step is to write the theory for the Bianchi type I model in the usual manner; that is, we calculate the corresponding energy-momentum tensor to the scalar field and give the equivalent Lagrangian density. Next, by means of a Legendre transformation, we proceed to obtain the canonical Lagrangian , from which the classical Hamiltonian can be found.

One of the simplest K-essence models, without self-interaction, has the following Lagrangian density: where is the scalar curvature and is an arbitrary function of the scalar field.

From the Lagrangian (2) we can build the complete action: where is the matter Lagrangian, is the cosmological constant Lagrangian, and is the determinant of the metric tensor. The field equations for this theory arewhere we work in units with and, as usual, the semicolon means a covariant derivative and a subscripted denotes differentiation with respect to .

The same set of (4a), (4b) is obtained if we consider the scalar field as part of the matter content; that is to say, with the corresponding energy-momentum tensor:

Considering the energy-momentum tensor of a barotropic perfect fluid, where is the four-velocity, which satisfies the relation , is the energy density, and is the pressure of the fluid. For simplicity we consider a comoving perfect fluid. The pressure, the energy density, and the four-velocity, corresponding to the energy-momentum tensor of the field , become thus, the barotropic parameter is and we notice that the case of a constant barotropic index (with the exception ) can be obtained by the function We have the following states in the evolution of our universe in this formalism: The mathematical analysis for the last two cases is very complicated in both regimes, classical and quantum one. For quantum radiation case, the resulting Wheeler-DeWitt equation appears as fractionary differential equation, and the results will be reported elsewhere. In [6], the authors present the analysis to radiation era using dynamical systems obtaining bouncing solutions. It is worthy to mention that (4b) was especialized in [17] to flat FRW cosmological model and solved using an extended tachyon field from which an extended Chaplygin gas can arise, which interpolates between a power law phase and a de Sitter phase. The author found the general solutions for linear expansion to the functional and particular and showed from the kinematic point of view that this model and the quintessence scalar field one driven by an exponential potential are the same, but they are dynamically nonequivalent because the K-field and the scalar potential are linked by the Einstein equation, similar to our next equations (19). However our results contain other functions ; then we argued that explanation to the corresponding observational data and the kinematic properties, since each formalism will have big differences.

1.1. Anisotropic Cosmological Bianchi Class A Models, = Constant

Considering the cosmological anisotropic Bianchi class A models with metric (23), (4b) in terms of becomes (here and all where appear the means, , with the usual cosmic time) and its corresponding solution becomes with a constant.

Note that (12) gives us the possible solution , as a function of the scale factor, and therefore the behavior of all physical properties of the K-essence (like , ) is completely determined by the function and does not depend on the evolution of the other types of energy density. The only dependence of the K-essence component on other components enters through in Bianchi class A cosmological models.

In the following we present the analysis when is a constant and generic function of the field assuming a Bianchi type I metric, which is the anisotropic generalization of flat FRW cosmological model, and we present the solution in quadrature form.

1.1.1. Quintessence Radiation Like Case: and = Constant

Using the equation for the energy kinetic we have the form then the field has the solution

1.1.2. Quintessence Stiff Fluid Like Case: and Constant

The field equations for this particular case areand the energy-momentum tensor (5) has the following form: In this new line of reasoning, action (3) can be rewritten as a geometrical part and matter content (usual matter plus a term that corresponds to the exotic scalar field component of the K-essence theory).

The equation of motion for the field (16b) has the following property, using the metric of the Bianchi type I model (however, this is satisfied by all cosmological Bianchi class A models): which can be integrated at once with the following result: where is an integration constant and has the same sign as . Considering the particular form of or with and being constants, the classical solutions for the field in quadrature are In the particular gauge , (20) is simplified to (remember that )

The complete solution to the quintessence scalar field depends strongly on the behavior of the scale factor of the cosmological model under consideration and particular form to the function ; Bianchi type model appears mimetic in the function, as in (15) and (19), except when one chooses a particular gauge which includes this function, as in (21).

In our particular case, it is evident that the contribution of the scalar field is equivalent to a stiff fluid with a barotropic equation of state . This is an instance of the results of the analysis of the energy-momentum tensor of a scalar field (17) by Madsen [18] for general relativity with scalar matter and by Pimentel [19] for the general scalar tensor theory. In both works a free scalar field is equivalent to a stiff matter fluid. In this way, we write action (3) in the usual form and, consequently, the classical equivalence between the two theories. We can infer that this correspondence is also satisfied in the quantum regime, so we can use this structure for the quantization program, where the ADM formalism is well known for different classes of matter [20].

This work is arranged as follows. In Section 2 we construct the Lagrangian and Hamiltonian densities for the anisotropic Bianchi type I cosmological model. In Section 3 we present some ideas as the anisotropic cosmological model can obtain its isotropization via the mean volume function and next we obtain the classical exact solution for all values in the gamma parameter. Finally, Section 3 is devoted to some final remarks.

2. Hamiltonian for the Bianchi Type I Cosmological Model

Let us recall here the canonical formulation in the ADM formalism of the diagonal Bianchi class A models. The metric has the form where is a scalar, is the lapse function, is a diagonal matrix, , are one-forms that characterize each cosmological Bianchi type model and obey , and are the structure constants of the corresponding invariance group. For the Bianchi type I model we have The total Lagrangian density then for this metric becomes using the standard definition of the momenta, , where , we obtain and, introducing them into the Lagrangian density, we obtain the canonical Lagrangian as , where we have used energy-momentum conservation, the law for a perfect fluid; ; we assumed an equation of state , so the corresponding Hamiltonian density is with .

2.1. Classical Equations

The corresponding Einstein field equations (16a) and (16b) for the anisotropic cosmological model Bianchi type I are the following (remember that the prime is the derivative over the time ): The solution of this last equation was putted in (19), The combination between the second and third equations gives us the solution for the anisotropic function ; also the sum of third and fourth equations, putting the solution, gives us the form of the function, where are integration constants. So, (29) are rewritten as

2.2. Isotropization

The current observations of the cosmic background radiation set a very stringent limit to the anisotropy of the universe [21]; therefore it is important to consider the anisotropy of the solutions. Recalling the Friedmann equation (constraint equation), we can see that isotropization is achieved when the terms with go to zero or are negligible with respect to the other terms in the differential equation. We find in the literature the criteria for isotropization, among others, , , that are consistent with our above remark. In the present case the comparison with the density should include the contribution of the scalar field. We define an anisotropic density that is proportional to the shear scalar: and will compare it with , , and . From the Hamilton Jacobi analysis we know that and the ratios are

Here we see that for expanding a universe the anisotropic density is dominated by the fluid density (with the exception of the stiff fluid) or by the term and then at late times the isotropization is obtained if the expansion goes to infinity. Hence it is necessary to determine when we have an ever expanding universe. Equation (34) in the new variable is

That is equivalent to the equation of motion in the coordinate of a particle under the potential with energy , where We can now have a qualitative idea of the different solutions from energy diagrams. We assume that is nonnegative since it is proportional to the energy density of the fluid. On the other hand, and are real and can take positive, null, or negative values.

From Figure 1 that is qualitatively correct for , we see that for negative all the expanding solutions will recollapse eventually, regardless of the sign of . For positive an expanding solution will expand forever. We also note that, when , that is, when the ghost contribution is dominant, there are contracting solutions that reach a minimum and then expand forever; these solutions do not have a big bang singularity.

2.2.1. Radiation Case,

Reproducing the set of equation for this case, we have the fact that the Friedmann-like equation is and, making the same analysis that is in the previous case, we have Also the relation between the anisotropic function and (32) is satisfied. In the last equation we have used (14).

Also, we can follow the same structure for the matter case and, following the Hamilton Jacobi analysis, we know that and the corresponding ratios are

Hence it is necessary to determine when we have an ever expanding universe. Equation (41) in the new variable

That is equivalent to the equation of motion in the coordinate of a particle under the potential with energy , where For this case, the qualitative analysis is the same when .

In the following we obtain exact solutions in order to give the volume function to each case to .

2.3. Exact Classical Solutions

In order to find the solutions for the remaining minisuperspace variables we employ the Einstein-Hamilton-Jacobi equation, which arises by making the identification in the Hamiltonian constraint , which results in In order to solve the above equation, we assume a solution of the form which results in the following set of ordinary differential equations: Here the are separation constants satisfying the relation , are real, should have the same signs as , and for consistency with (19) we have . Recalling the expressions for the momenta we can obtain solutions for equations (47)–(48) in quadrature; in particular

We already know the solution for (48). As can be seen from (50), in order to obtain solutions for one needs to find a solution for , which can be obtained from (49).

Equation (49) does not have a general solution; however, it is possible to find solutions for particular values of the barotropic parameter with .

(1)   and . Equation (49) can be written as When we consider the time transformations , and the change of variable , this equation has the solution where and . With this, the time transformation becomes And the closed form is [22] where is a hypergeometric function. We also have The anisotropy functions and the scalar field are given by As a concrete example we consider the particular value and then and becomes If isotropization is possible, this volume function would not make it quick, and the integral becomes So, the classical solutions for the anisotropic function and field are

(2)   and . In this case (49) is with that we assume positive. By integrating we obtain If isotropization is possible, this volume function would not make it quick. For the anisotropic functions we have and the scalar field is given by

(3)   and . Equation (49) has the form where : Solving for The inverse volume function is If isotropization is possible, the corresponding volume function would make it quick and its integral becomes And, as the anisotropic function is dependent on this integral, so Also, the field uses this integral (see (20)), and the corresponding solutions become with the condition .

Also we can consider the special case, when or : We have the following expression: The isotropization is possible for this case, because the corresponding volume function would make it quick and its integral becomes So, the anisotropic function and the function become for this case

(4)   and . For this case, (49) is with and .

The function becomes and we have Also, the isotropization is possible, because the corresponding volume function would make it quick. For the anisotropic function we have and for the function

(5)   and . Equation (49) becomes where . In this case also we have two possible solutions depending on the value of the cosmological constant.

(i) . The solution becomes so the function is For this case the isotropization is possible and the corresponding volume function would make it quick. For the anisotropic function we have and for the field

(ii) . The corresponding solutions are The function is As the volume has an oscillatory behavior, the isotropization does not yield for this case, and for completeness we calculate The anisotropic function , And the field

3. Final Remarks

In this work we present the study of the classical cosmological anisotropic Bianchi type I in the K-essence formalism. In previous work made by Chimento and Forte [16], they present the possible isotropization of this model. Our goal in this work is that we obtain the corresponding classical solutions for a barotropic perfect fluid and cosmological term which mimics the scalar field in (1). In the case of and we obtain the solutions in closed form. With these solutions we can validate our qualitative analysis on isotropization of the cosmological model, implying that this occurs when the volume is large in the corresponding time evolution. So, only one solution does not present the large volume, when in stiff matter era in the ordinary matter content. We include a qualitative analysis to Friedmann equation when it is written as an equation for the volume that is equivalent to the equation of motion of a particle under a potential and we conclude the same about the isotropization of this anisotropic model, considering the stiff matter and radiation cases. In [17] the author shows from the kinematic point of view the model when is linear (as in our case) and the particular value to the function , using our notation, and the quintessence scalar field one driven by an exponential potential are the same, but they are dynamically nonequivalent because the K-field and the scalar potential are linked by the Einstein equation. However our results contain other functions ; then we argued that explanation to the corresponding observational data and the kinematic properties, since each formalism will have big differences. All of these comments occur because the complete solution to the quintessence scalar field depends strongly on the behavior of the scale factor of the cosmological model under consideration and particular form to the function ; Bianchi type model appears mimetic in the function, as in (15) and (19), except when one chooses a particular gauge which includes this function, as in (21). In the quantum analysis for this model, considering the scalar field, the solutions are similar to those found in the Bianchi type IX cosmological model [23]; you can see (32). For quantum radiation case, the resulting Wheeler-DeWitt equation appears as fractionary differential equation, and the results will be reported elsewhere.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was partially supported by CONACYT Grants 167335 and 179881 and PROMEP Grants UGTO-CA-3 and UAM-I-43. This work is part of the collaboration within the Instituto Avanzado de Cosmología and Red PROMEP: Gravitation and Mathematical Physics under project Quantum Aspects of Gravity in Cosmological Models, Phenomenology, and Geometry of Space-Time. Many calculations were done by Symbolic Program REDUCE 3.8.