[ProQuest: [...] denotes non US-ASCII text; see PDF]
Bogdan Danila 1 and Tiberiu Harko 2,3 and Man Kwong Mak 4,5 and Praiboon Pantaragphong 6 and Sorin V. Sabau 7
Academic Editor:Douglas Singleton
1, Astronomical Observatory, 19 Ciresilor Street, Cluj-Napoca, Romania
2, Department of Physics, Babes-Bolyai University, Kogalniceanu Street 400084, Cluj-Napoca, Romania
3, Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
4, School of Science and Technology, The Open University of Hong Kong, Homantin, Kowloon, Hong Kong
5, Departamento de Física, Facultad de Ciencias Naturales, Universidad de Atacama, Copayapu 485, Copiapó, Chile
6, Mathematics Department, King Mongkut's Institute of Technology, Ladkrabang, Bangkok 10520, Thailand
7, School of Science, Department of Mathematics, Tokai University, Sapporo 005-8600, Japan
Received 23 July 2016; Accepted 18 September 2016; 26 October 2016
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 SCOAP3 .
1. Introduction
A large number of cosmological observations, obtained initially from distant Type Ia Supernovae, have convincingly proven that the Universe has undergone a late-time accelerated expansion [1-4]. In order to explain these observations, a deep change in our paradigmatic understanding of the cosmological dynamics is necessary, and many ideas have been put forward to address them. The "standard" explanation of the late-time acceleration is based on the assumption of the existence of a mysterious component, called dark energy, which is responsible for the observed characteristics of the late-time evolution of the Universe.
On the other hand, the combination of the results of the observations of high redshift supernovae and of the WMAP and of the recently released Planck data indicates that the location of the first acoustic peak in the power spectrum of the Cosmic Microwave Background Radiation is consistent with the prediction of the inflationary model for the density parameter Ω , according to which Ω = 1 . The cosmological observations also provide strong evidence for the behavior of the parameter w = p / ρ of the equation of state of the cosmological fluid, where p is the pressure and ρ is the density, as lying in the range - 1 <= w = p / ρ < - 1 / 3 [5].
In order to explain the observed cosmological dynamics, it is assumed usually that the Universe is dominated by two main components: cold (pressureless) dark matter (CDM) and dark energy (DE) with negative pressure, respectively. CDM contributes Ω m ~ 0.3 [6], and its introduction is mainly motivated by the necessity of theoretically explaining the galactic rotation curves and the large scale structure formation. On the other hand, DE represents the major component of the Universe, providing Ω D E ~ 0.7 . Dark energy is the major factor determining the recent acceleration of the Universe, as observed from the study of the distant Type Ia Supernovae [5]. Explaining the nature and properties of dark energy has become one of the most active fields of research in cosmology and theoretical physics, with a huge number of proposed DE models (for reviews, see, e.g., [7-11]).
One interesting possibility for explaining DE is cosmological models containing a mixture of cold dark matter and quintessence, representing a slowly varying, spatially inhomogeneous component [12]. From a theoretical as well as a particle physics point of view, the idea of quintessence can be implemented by assuming that it is the energy associated with a scalar field Q , having a self-interaction potential V ( Q ) . When the potential energy density V ( Q ) of the quintessence field is greater than the kinetic one, it follows that the pressure p = Q 2 / 2 - V ( Q ) associated with the quintessence Q -field is negative. The properties of the quintessential cosmological models have been actively considered in the physical literature (for a recent review, see [13]). As opposed to the cosmological constant of standard general relativity, the equation of state of the quintessence field changes dynamically with time [14]. Alternative models, in which the late-time acceleration can be driven by the kinetic energy of the scalar field, called k -essence models, have also been proposed [15-19].
Scalar fields [varphi] that are minimally coupled to gravity via a negative kinetic energy can also explain the recent acceleration of the Universe. Interestingly enough, they allow values of the parameter of the equation of state with w < - 1 . These types of scalar fields, known as phantom fields, have been proposed in [20]. For phantom scalar fields, the energy density and pressure are given by ρ [varphi] = - [varphi] 2 / 2 + V ( [varphi] ) and p [varphi] = - [varphi] 2 / 2 - V ( [varphi] ) , respectively. The interesting properties of phantom cosmological models for dark energy have been investigated in detail in [21-28]. Recent cosmological observations show that at some moment during the cosmological evolution the value of the parameter w may have crossed the standard value w = - 1 , corresponding to the general relativistic cosmological constant Λ . This cosmological situation is called the phantom divide line crossing [28]. In the case of scalar field models with cusped potentials, the crossing of the phantom divide line was investigated in [27]. Another alternative way of explaining the phantom divide line crossing is to model dark energy by a scalar field, which is nonminimally coupled to gravity [27].
Scalar fields are also assumed to play a fundamental role in the evolution of the very early Universe, playing a major role in the inflationary scenario [29, 30]. Originally, the idea of inflation was proposed to provide solutions to the singularity, flat space, horizon, and homogeneity problems, to the absence of magnetic monopoles, and to the problem of large numbers of particles [31, 32]. However, presently it is believed that the most important feature of inflation is the generation of both initial density perturbations and the background cosmological gravitational waves. These important cosmological parameters can be determined in many different ways, like, for example, through the study of the anisotropies of the microwave background radiation, the analysis of the local (peculiar) velocity galactic flows, of the clustering of galaxies, and the determination of the abundance of gravitationally bound structures of different types, respectively [32].
In many inflationary models, the dynamical evolution of the early Universe is driven by a single scalar field, called the inflaton, with the inflaton rolling in some underlying self-interaction potential [29-32]. One common approximation in the study of the inflationary evolution is the slow-roll approximation, which can be successfully used in two separate contexts. The first situation is in the study of the classical inflationary dynamics of expansion in the lowest order approximation. Hence, this implies that the contribution of the kinetic energy of the inflaton field to the expansion rate is ignored. The second situation is represented by the calculation of the perturbation spectra. The standard expressions deduced for these spectra are valid in the lowest order in the slow-roll approximation [33]. Finding exact inflationary solutions of the gravitational field equations for different types of scalar field potentials is also of great importance for the understanding of the dynamics of the early Universe. Such exact solutions have been found for a large number of inflationary potentials. Moreover, the potentials allowing a graceful exit from inflation have been classified [34].
Hence, the theoretical investigation of the scalar field models is an essential task in cosmology. Among the various methods used to study the properties of scalar fields, the methods based on the applications of the mathematical formalism of the qualitative study of dynamical systems are of considerable importance.
The usefulness of dynamical systems formulation of physical models is mainly determined by their powerful predictive power. This predictive power is essentially determined by the stability of their solutions. In a realistic physical system, due to the limited precision of the measurements, some uncertainties in the initial conditions always exist. Therefore, a physically meaningful mathematical model must also offer detailed and useful information on the evolution of the deviations of the possible trajectories of the dynamical system from a given reference trajectory. Hence, an important requirement in mathematical modelling is the understanding of the local stability of the physical and cosmological processes. This information on the system behavior is as important as the understanding of the late-time deviations. The global stability of the solutions of the dynamical systems described by systems of nonlinear ordinary differential equations is analyzed in the framework of the well-known mathematical theory of Lyapunov stability. In this mathematical approach, the fundamental quantities that measure exponential deviations from a given trajectory are the so-called Lyapunov exponents [35, 36]. It is usually very difficult to determine the Lyapunov exponents analytically for a given dynamical system, and thus one must resort to numerical methods. On the other hand, the important problem of the local stability of the solutions of dynamical systems, described by ordinary differential equations, is less understood.
Cosmological models have been intensively investigated by using methods from dynamical systems and Lyapunov stability theory [37-53]. In particular, phase space analysis proved to be a very useful method for the understanding of the cosmological evolution. When studying the evolution of cosmological models, the dynamical equations can be represented by an autonomous dynamical system, described by a set of coupled (usually strongly nonlinear) differential equations for the physical parameters. This representation allows the study of the Lyapunov stability of the model, without explicitly solving the field equations for the basic variables. Furthermore, the importance of the Lyapunov analysis is related to the fact that stationary points of the dynamical system correspond to exact or approximate analytic solutions of the field equations. Thus, the dynamical system formulation provides a useful tool for obtaining exact or approximate solutions of the field equations in cosmologically interesting situations.
Even though the mathematical methods of the Lyapunov stability analysis are well established, the study of the stability of the dynamical systems from different points of view is extremely important. The comparison of the results of the alternative approach with the corresponding Lyapunov exponents analysis can provide a deeper understanding of the stability properties of the system. An alternative and very powerful method for the study of the systems of the ordinary differential equations is represented by the so-called Kosambi-Cartan-Chern (KCC) theory, which was initiated in the pioneering works of Kosambi [54], Cartan and Kosambi [55], and Chern [56], respectively. The KCC theory was inspired and influenced by the geometry of the Finsler spaces (for a recent review of the KCC theory, see [57]). From a mathematical point of view, the KCC theory is a differential geometric theory of the variational equations for the deviations of the whole trajectory with respect to the nearby ones [58]. In the KCC geometrical description of the systems of ordinary differential equations, one associates a nonlinear connection and a Berwald-type connection with the system of equations. With the use of these geometric quantities, five geometrical invariants are obtained. The most important invariant is the second invariant, also called the curvature deviation tensor, which gives the Jacobi stability of the system [57-60]. The KCC theory has been applied for the study of different physical, biochemical, or technical systems (see [59-77]).
An alternative geometrization method for dynamical systems, with applications in classical mechanics and general relativity, was proposed in [78, 79] and further investigated in [80-84]. The Henon-Heiles system and Bianchi type IX cosmological models were also investigated within this framework. In particular, in [80] a theoretical approach based on the geometrical description of dynamical systems and of their chaotic properties was developed. For the base manifold, a Finsler space was introduced, whose properties allow the description of a wide class of physical systems, including those with potentials depending on time and velocities, for which the Riemannian approach is unsuitable.
It is the purpose of the present paper to consider a systematic investigation of the Jacobi stability properties of the flat homogeneous and isotropic general relativistic cosmological models. By starting from the standard Friedmann equations, we perform, as a first step in our analysis, a "second geometrization" of these equations, by associating with them a nonlinear connection and a Berwald connection, respectively. This procedure allows obtaining the so-called KCC invariants of the Friedmann equations. The second invariant, called the curvature deviation tensor, gives the Jacobi stability properties of the cosmological model. The KCC theory can be naturally applied to systems of second-order ordinary differential equations. The Friedmann equations can be formulated as second-order differential equations, similarly to the Klein-Gordon equation describing the scalar field. Therefore, the KCC theory can be applied to matter and scalar field dominated cosmological models. We obtain the general condition for the Jacobi stability of scalar fields, which is described by two inequalities involving the second and the first derivatives of the scalar field potential, the energy density of the field, and the time derivative of the field itself. The geodesic deviation equations describing the time variation of the deviation vector are also obtained. As an application of the developed formalism, we investigate the stability properties of the scalar field cosmological models with exponential and Higgs type potentials, respectively. It turns out that the exponential potential scalar field is Jacobi unstable during its entire evolution, while the time evolution of the scalar field cosmological models with Higgs potential shows complicated dynamics with alternating stable and unstable Jacobi phases. The Jacobi stability properties of the Higgs type models are determined by the numerical value of the ratio of the self-coupling constant and the square of the mass of the Higgs particle.
The Lyapunov stability properties of the scalar field cosmological models are usually investigated by reformulating the evolution equation as a set of three first-order ordinary differential equations. In order to apply the KCC theory to such systems, they must be lifted to the tangent bundle. From mathematical point of view, this requires taking the time derivative of the first-order equations, so that their "second geometrization" can be easily performed. We consider in detail the Jacobi stability properties of the phantom quintessence and tachyon scalar field cosmological models. We study in detail the Jacobi stability condition of these models, and we find that they are Jacobi unstable during the entire expansionary cosmological evolution.
The present paper is organized as follows. We review the basic ideas and the mathematical formalism of the KCC theory in Section 2. The Jacobi stability analysis of the homogenous isotropic flat cosmological models by using the second-order formulation of the dynamics is performed in Section 3. We consider both cases of the matter dominated and scalar field dominated cosmological models. As an application of the developed formalism, we investigate in detail the Jacobi stability of the scalar fields with exponential potential and Higgs potential, respectively. The Jacobi stability of the first-order dynamical system formulation of scalar field cosmological models is considered in Section 4, in which the KCC geometrization of the phantom quintessence and tachyonic scalar field models is analyzed in detail. We discuss and conclude our results in Section 5. The KCC geometric quantities giving the geometric description of the phantom quintessence and tachyon scalar field cosmologies are presented in Appendix A and Appendix B, respectively.
2. Kosambi-Cartan-Chern (KCC) Theory and Jacobi Stability
In the present section, we briefly summarize the basic concepts and results of the KCC theory (for a detailed presentation, see [57, 58]).
2.1. Dynamical Systems as Paths of a Semispray
In the following, we denote by M a real, smooth n -dimensional manifold and by T M its tangent bundle. On an open connected subset Ω of the Euclidian ( 2 n + 1 ) dimensional space R n × R n × R 1 , we introduce a 2 n + 1 -dimensional coordinates system ( x i , y i , t ) , i = 1,2 , ... , n , where ( x i ) = ( x 1 , x 2 , ... , x n ) , ( y i ) = ( y 1 , y 2 , ... , y n ) , and t is the time t . The coordinates y i are defined as [figure omitted; refer to PDF]
We assume that the time t is an absolute invariant, and therefore the only admissible change of coordinates will be [figure omitted; refer to PDF]
Definition 1 (see [85]).
A deterministic dynamical system is a formal set of rules that describe the evolution of points in some set S with respect to an external, discrete, or continuous time parameter t running in another set T .
In a rigorous mathematical formulation, a dynamical system is a map [85]: [figure omitted; refer to PDF] which satisfies the condition [varphi] ( t , · ) [composite function] [varphi] ( s , · ) = [varphi] ( t + s , · ) for all times t , s ∈ T . In order to model realistic dynamical systems or physical processes additional structures must be added to the above definition.
In many cases, the equations of motion of a dynamical system can be derived from Lagrangian L via the Euler-Lagrange equations: [figure omitted; refer to PDF] where F i , i = 1,2 , ... , n , is the external force. For regular Lagrangian L , the Euler-Lagrange equations given by (4) are equivalent to a system of second-order differential equations [86]: [figure omitted; refer to PDF] where each function G i ( x j , y j , t ) is C ∞ in a neighborhood of some initial conditions ( ( x ) 0 , ( y ) 0 , t 0 ) in Ω .
A vector field S on T M of the form [figure omitted; refer to PDF] is called a semispray [87-89]. The functions G i ( x i , y i ) are the local coefficients of the semispray, and they are defined on domains of local charts. In the particular case in which the coefficients G i = G i ( x i , y i ) are homogeneous of degree two in y i , the vector field S is called a spray .
Definition 2.
A path of the semispray S is defined as a curve c : t [arrow right] x i ( t ) on M , with the property that its lift c [variant prime] : t [arrow right] ( x i ( t ) , x i ( t ) ) to T M is an integral curve of S , that is, a curve satisfying the equation [figure omitted; refer to PDF]
Conversely, for any system of ordinary differential equations of form (7), which is globally defined, the functions G i define a semispray on T M [88, 89].
More generally, one can start from an arbitrary system of second-order differential equations of form (5) , where no a priori given Lagrangian function is assumed, and study the behavior of its trajectories by analogy with the trajectories of the Euler-Lagrange equations.
2.2. The KCC Geometrization of Dynamical Systems
To associate a geometrical structure with the dynamical system defined by (5), we introduce first a nonlinear connection N on M , with coefficients N j i , defined as [86] [figure omitted; refer to PDF]
The nonlinear connection can be understood in terms of a dynamical covariant derivative ∇ N [85]: for two vector fields v , w defined over M , we introduce the covariant derivative ∇ N as [figure omitted; refer to PDF]
For N j i ( x , y ) = Γ j l i ( x ) y l , (9) reduces to the definition of the covariant derivative for the special case of a linear connection.
For nonsingular coordinate transformations given by (2), the KCC-covariant differential of an arbitrary vector field θ i ( x ) on the open subset Ω ⊆ R n × R n × R 1 is defined as [57, 59] [figure omitted; refer to PDF]
For θ i = y i , we obtain [figure omitted; refer to PDF] where the contravariant vector field [...] i on Ω is called the first KCC invariant.
2.2.1. The Curvature Deviation Tensor
Let us now vary the trajectories x i ( t ) of system (5) into nearby ones according to [figure omitted; refer to PDF] where ( η ) is a small parameter and ξ i ( t ) are the components of some contravariant vector field defined along the path x i ( t ) . Substituting (12) into (5) and taking the limit η [arrow right] 0 , we obtain the variational equations [58-61] [figure omitted; refer to PDF]
By using the KCC-covariant differential, we can write (13) in the covariant form: [figure omitted; refer to PDF] where we have denoted [figure omitted; refer to PDF] and we have introduced the Berwald connection G j l i , defined as [57-61, 86] [figure omitted; refer to PDF]
The tensor P j i is called the second KCC invariant, or the deviation curvature tensor , and (14) is called the Jacobi equation , respectively. In either Riemann or Finsler geometry, when system (5) describes the geodesic equations in the given geometry, (14) represents the Jacobi field equation.
The trace P of the curvature deviation tensor is obtained as [figure omitted; refer to PDF]
The third, fourth, and fifth invariants of system (5) are defined as [58] [figure omitted; refer to PDF]
The third invariant P j k i can be interpreted as a torsion tensor. The fourth and fifth invariants P j k l i and D j k l i are called the Riemann-Christoffel curvature tensor and the Douglas tensor, respectively [57, 58]. In Berwald spaces, these tensors always exist. Generally, they can be used to describe the geometrical properties of systems of second-order differential equations.
2.3. The Definition of Jacobi Stability
In many physical applications, the behavior of the trajectories of the dynamical system (5) in a vicinity of a point x i ( t 0 ) is of major interest. In the following, for simplicity, we take t 0 = 0 . The trajectories x i = x i ( t ) can be considered as curves in the Euclidean space ( R n , ( · , · ) ) , where ( · , · ) defines the canonical inner product of R n . As for the deviation vector ξ , we assume that it satisfies the natural initial conditions ξ ( 0 ) = O and ξ ( 0 ) = W ≠ O , where O ∈ R n is the null vector [57-60].
We describe the focusing tendency of the trajectories around t 0 = 0 as follows: if ( ξ ( t ) ) < t 2 , t [approximate] 0 + , the trajectories are bunching together. On the other hand, if ( ξ ( t ) ) > t 2 , t [approximate] 0 + , the trajectories are dispersing [57-60]. Alternatively, the focusing tendency of the trajectories can be described in terms of the deviation curvature tensor: the trajectories of the system of second-order differential equations (5) are bunching together for t [approximate] 0 + if and only if the real parts of the eigenvalues of the deviation curvature tensor P j i ( 0 ) are strictly negative. On the other hand, the trajectories are dispersing if and only if the real parts of the eigenvalues of the tensor P j i ( 0 ) are strictly positive [57-60].
Based on the above intuitive considerations, we can define the Jacobi stability for a second-order system of differential equations as follows [57-60].
Definition 3.
If the system of second-order differential equations (5) satisfies the initial conditions [figure omitted; refer to PDF] with respect to the norm ( · ) induced by a positive definite inner product, then we call the trajectories of (5) Jacobi stable if and only if the real parts of the eigenvalues of the deviation tensor P j i are strictly negative everywhere. Otherwise, the trajectories are called Jacobi unstable.
The focusing/dispersing behavior of the trajectories of a system of second-order ordinary differential equations near the origin is represented in Figure 1.
Figure 1: Behavior of the trajectories near zero.
[figure omitted; refer to PDF]
2.4. The Deviation Curvature Tensor for Two- and Three-Dimensional Dynamical Systems
In the important two-dimensional case, the curvature deviation tensor can be written in a matrix form as [figure omitted; refer to PDF] with the eigenvalues given by [figure omitted; refer to PDF]
The eigenvalues of the curvature deviation tensor can be obtained as solutions of the quadratic equation: [figure omitted; refer to PDF]
A very powerful algebraic method to obtain the signs of the eigenvalues of the curvature deviation tensor is represented by the Routh-Hurwitz criteria [90]. According to these criteria, all of the roots of the polynomial P ( λ ) are negatives or have negative real parts if the determinants of all Hurwitz matrices det [...] H j , j = 1,2 , ... , n , are strictly positive. For n = 2 , corresponding to the case of (22), the Routh-Hurwitz criteria take the simple form [figure omitted; refer to PDF]
The curvature properties along a given geodesic are described by the eigenvalues of the deviation curvature tensor λ ± , which are invariant functions on the tangent space. Moreover, once they are known, we can introduce two quantities that can characterize the way the geodesic explores the base manifold. They are defined as the (half) of the Ricci curvature scalar along the flow, κ , and the anisotropy θ , given by [figure omitted; refer to PDF] respectively.
In the case of a three-dimensional dynamical system with n = 3 , the characteristic equation of the matrix of the curvature deviation tensor becomes [figure omitted; refer to PDF] Therefore, the conditions of the Jacobi stability for a three-dimensional system of second-order ordinary differential equations can be formulated as follows: [figure omitted; refer to PDF]
3. Jacobi Stability Analysis of Isotropic Matter Dominated and Scalar Field Cosmologies
In the present section, we use the KCC approach for the study of the dynamical properties of the matter dominated and scalar field cosmologies. We explicitly obtain the nonlinear and Berwald connections and the deviation curvature tensors. The eigenvalues of the deviation curvature tensor are also obtained, and we study their properties in the equilibrium points of the matter dominated model. The study of the sign of the eigenvalues allows us to obtain the Jacobi stability properties of the fixed points of the matter field dominated cosmological models. Next, we proceed to a detailed analysis of the scalar field cosmologies in the framework of the KCC theory. A full "second geometric" description is introduced, and the time evolution of the relevant physical and geometrical parameters is obtained. In particular, the nature of the Jacobi stability is analyzed in detail. In the present study, we restrict our analysis to the case of the flat Friedmann-Robertson-Walker metric, given by [figure omitted; refer to PDF] where a is the scale factor.
3.1. Jacobi Stability of Matter Dominated Cosmological Models
For a Universe filled with pressureless dust and radiation only, the cosmological expansion is described mathematically by the Friedmann equations, which take the well-known form [figure omitted; refer to PDF] where a dot denotes the derivative with respect to the time t . In (28), H = a / a denotes the Hubble function, ρ m is the baryonic matter energy density, ρ r is the energy density of the radiation, and Λ is the cosmological constant.
3.1.1. Friedmann Equations as an Autonomous Dynamical System
In order to reformulate the cosmological evolutions equations as a dynamical system, we need to introduce first the density parameters ( Ω m , Ω r , Ω Λ ) of the matter, radiation and cosmological constant, defined as [figure omitted; refer to PDF] The density parameters satisfy the normalization relation [figure omitted; refer to PDF] As the basic variables ( x , y ) in the phase space we adopt the quantities x ≡ Ω r and y ≡ Ω Λ , respectively [52]. Then, the density parameter of the matter is given by Ω m = 1 - x - y , and the range of the variables ( x , y , Ω m ) is 0 <= x <= 1 , 0 <= y <= 1 , and 0 <= Ω m <= 1 . To describe the cosmological dynamics, we define the physically significant phase space as Φ = ( ( x , y ) : x + y <= 1 , 0 <= x <= 1 , 0 <= y <= 1 ) . Next, we take the time derivatives of x and y with respect to the time t , and, after introducing the new time variable τ = ln [...] a ( t ) , we can formulate the Friedmann equations as an autonomous dynamical system given by [52] [figure omitted; refer to PDF] The critical points of system (31) and (32) in the phase space region defined by Ψ are obtained by solving the algebraic equations x ( 1 - x + 3 y ) = 0 and ( 3 + x - 3 y ) y = 0 , respectively, and are given by [figure omitted; refer to PDF]
The Lyapunov stability properties of the system follow from the study of the Jacobian matrix: [figure omitted; refer to PDF]
The critical points of system (31) and (32) have a clear cosmological interpretation. Thus, the critical point P d S = ( 0,1 ) , with Ω Λ = 1 , has a ( t ) [proportional, variant] e x p ( Λ / 3 t ) and is associated with an accelerated de Sitter type expansion, being a future attractor [53]. The critical point P r = ( 1,0 ) , with Ω r = 1 and a ( t ) [proportional, variant] t , corresponds to the radiation-dominated era in the cosmological evolution of the Universe and is a source point or a past attractor. Finally, the critical point P m ( 0,0 ) , with Ω m = 1 and a ( t ) [proportional, variant] t 2 / 3 , corresponds to the decelerating matter dominated phase of the cosmological expansion. It turns out that P m is a saddle critical point [53].
3.1.2. The KCC Geometrization and the Jacobi Stability of the Friedmann Equations
In order to apply the KCC theory to the cosmological dynamical system given by (31) and (32), we first relabel the variables as x ≡ x 1 and y ≡ x 2 . We also denote y 1 = d x 1 / d τ and y 2 = d x 2 / d τ , respectively. Hence, we obtain [figure omitted; refer to PDF]
Next, we take the derivative with respect to τ of (35) and (36), respectively, thus obtaining the following lift on the tangent bundle of the cosmological dynamical system: [figure omitted; refer to PDF]
By comparison with (5), we obtain immediately [figure omitted; refer to PDF]
Hence, the components of the nonlinear connection ( N ) = N i j associated with the matter dominated cosmological dynamical system in the presence of a cosmological constant are obtained as [figure omitted; refer to PDF]
For the components of the deviation tensor P = P i j , we obtain [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the Hessian of f and H g is the Hessian of g . Explicitly, the curvature deviation tensor for the Friedmann cosmological dynamical system can be obtained as [figure omitted; refer to PDF]
Evaluating P at the critical points gives [figure omitted; refer to PDF] where A = ( J ) ( y 1 = 0 , y 2 = 0 , x 1 = x c r 1 , x 2 = x c r 2 ) .
3.1.3. Jacobi Stability of the Critical Points of the Matter Dominated Cosmological Models
In order to obtain the Jacobi stability of the critical points of the cosmological model described by (31) and (32), we need to compute the numerical values of the curvature deviation tensor at the critical points. At the critical point P d S = ( x 1 = 0 , x 2 = 1 ) , the curvature deviation tensor takes the form [figure omitted; refer to PDF] and has the eigenvalues λ 1 = 4 and λ 2 = 9 . Hence, it follows that the critical point P d S of the standard Λ CDM cosmological model is Jacobi unstable . For the critical point P r = ( x 1 = 1 , x 2 = 0 ) , we find [figure omitted; refer to PDF] with the corresponding eigenvalues of the curvature deviation tensor obtained as λ 1 = 16 and λ 2 = 1 / 4 . Hence, the critical point P r is also Jacobi unstable . Finally, for the last critical point P m = ( 0,0 ) , we obtain [figure omitted; refer to PDF] with the eigenvalues λ 1 = 9 and λ 2 = 1 / 4 . Hence, we obtain the result that all the critical points of the dynamical system equivalent to the cosmological Friedmann equations are Jacobi unstable .
3.2. The Cosmological Evolution Equations in the Presence of Minimally Coupled Scalar Fields
Let us consider a rather general class of scalar field models, minimally coupled to the gravitational field, for which the Lagrangian density in the Einstein frame reads [figure omitted; refer to PDF] where R is the curvature scalar, [varphi] is the scalar field, V ( [varphi] ) is the self-interaction potential, and κ = 8 π G / c 4 is the gravitational coupling constant, respectively. In the following, we use natural units with c = 8 π G = [Planck's over 2pi] = 1 , and we adopt as the signature for the metric ( + 1 , - 1 , - 1 , - 1 ) , as is common in particle physics.
For a flat FRW scalar field dominated Universe, the evolution of a cosmological model is determined by the system of the field equations [figure omitted; refer to PDF] and the evolution equation for the scalar field [figure omitted; refer to PDF] where the overdot denotes the derivative with respect to the time-coordinate t and the prime denotes the derivative with respect to the scalar field [varphi] , respectively. By substituting a / a from (49) into (50) and (51), we can reformulate the dynamics of the scalar field cosmological models in terms of two second-order nonlinear ordinary differential equations, given by [figure omitted; refer to PDF] respectively.
3.2.1. The Nonlinear and Berwald Connections and the KCC Invariants of the Scalar Field Cosmological Models
In the following, we introduce a new notation for the dependent variables a and [varphi] and for their time derivatives, respectively, as [figure omitted; refer to PDF]
The cosmological dynamics of scalar field dominated Universes can be formulated as a second-order differential system, given by two second-order differential equations of the form [figure omitted; refer to PDF]
From (52) and (53), it follows immediately that [figure omitted; refer to PDF] respectively. Therefore, we first obtain the components of the nonlinear connection as [figure omitted; refer to PDF]
For the nonzero components of the Berwald connection, defined as [figure omitted; refer to PDF] we obtain [figure omitted; refer to PDF]
The components of the first KCC invariant of the minimally coupled scalar field model are obtained as [figure omitted; refer to PDF] respectively.
The components of the curvature deviation tensor for minimally coupled scalar field cosmological models are given by [figure omitted; refer to PDF] For the trace of the curvature deviation tensor, we obtain [figure omitted; refer to PDF] while for χ = P 1 1 P 2 2 - P 2 1 P 1 2 = P 1 1 P 2 2 we have [figure omitted; refer to PDF]
Therefore, if the conditions P 1 1 + P 2 2 < 0 and P 1 1 P 2 2 - P 2 1 P 1 2 > 0 are simultaneously satisfied, the scalar field cosmological model is Jacobi stable . These conditions allow us to formulate the following.
Jacobi Stability Condition of Isotropic and Homogeneous Scalar Field Cosmological Models . If the parameters ( V ( [varphi] ) , ρ [varphi] , p [varphi] , [varphi] ) of a homogeneous scalar field in an isotropic flat FRW geometry simultaneously satisfy the conditions [figure omitted; refer to PDF] the corresponding cosmological model is Jacobi stable, but Jacobi unstable otherwise.
For the variational differential equations determining the deviation vector ξ , we obtain [figure omitted; refer to PDF] respectively.
In the case of the isotropic cosmological scalar field models, the third, fourth, and fifth KCC invariants, as defined by (18), are identically equal to zero.
3.2.2. Applications: The Case of the Scalar Field with Exponential Self-Interaction Potential
As an application of the KCC geometrization of the scalar field cosmological models, we will consider the case of a scalar field with an exponential self-interaction potential of the form [figure omitted; refer to PDF] where V 0 and λ are constants. The cosmological equations describing the time evolution of this scalar field model are [figure omitted; refer to PDF] respectively. By introducing a new time variable τ = V 0 t , it turns out that the system of (70) and (71) can be written as [figure omitted; refer to PDF]
The system of (72)-(73) must be integrated with the initial conditions a ( 0 ) = a 0 , a ( 0 ) = a 0 H 0 , [varphi] ( 0 ) = [varphi] 0 , and [varphi] ( 0 ) = [varphi] 0 , respectively.
The variations of the scale factor and of the scalar field with exponential potential are presented, for different values of the parameter λ , in Figure 2. In order to numerically solve the field equations, we have used the initial conditions a ( 0 ) = 0 , [varphi] ( 0 ) = 0.001 , a ( 0 ) = 1 , and [varphi] ( 0 ) = - 0.001 , respectively.
Figure 2: Variation of the scale factor of the Universe (a) and of the scalar field (b) as a function of the dimensionless time τ for the scalar field with exponential potential, for different values of the parameter λ : λ = - 1 (solid curve), λ = - 2 (dotted curve), λ = - 3 (dashed curve), and λ = - 2 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
As one can see from Figure 2, the scale factor is a monotonically increasing function of time, while the scalar field also increases during the cosmological evolution. The time variation of the scalar field potential is presented in Figure 3.
Figure 3: Variation of the scalar field potential V = exp [...] ( - λ [varphi] ) as a function of the dimensionless time τ for different values of the parameter λ : λ = - 1 (solid curve), λ = - 2 (dotted curve), λ = - 3 (dashed curve), and λ = - 2 (dashed dotted curve).
[figure omitted; refer to PDF]
The scalar field potential is a monotonically decreasing function of time, which tends, in the large time limit, to zero. The time variations of the KCC invariants 2 κ = P 1 1 + P 2 2 and χ = P 1 1 P 2 2 - P 2 1 P 2 1 are represented in Figure 4.
Figure 4: Time variations of the KCC invariants 2 κ = P 1 1 + P 2 2 (a) and χ = P 1 1 P 2 2 - P 2 1 P 2 1 (b) for a Universe filled with a scalar field with exponential potential, for different values of the parameter λ : λ = - 1 (solid curve), λ = - 2 (dotted curve), λ = - 3 (dashed curve), and λ = - 2 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
In order for the considered model of the Universe to be stable, the invariants must simultaneously satisfy the conditions 2 κ < 0 and χ > 0 , respectively. As one can see from the figures, from the chosen set of parameters, these conditions are not satisfied for any interval of time during the cosmological evolution. Therefore, it follows that during its entire evolution an exponential potential scalar field Universe is in a Jacobi unstable state. Finally, in Figure 5 we present the time variation of the components of the deviation vector ξ i , i = 1,2 .
Figure 5: Time variations of the deviation vector components ξ 1 (a) and ξ 2 (b) for a Universe filled with a scalar field with exponential potential, for different values of the parameter λ : λ = - 1 (solid curve), λ = - 2 (dotted curve), λ = - 3 (dashed curve), and λ = - 2 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
The component ξ 1 of the deviation vector increases exponentially in time, indicating that the trajectories do diverge exponentially near the origin. Thus, this result also confirms the presence of the Jacobi instability for the exponential potential scalar field cosmology.
3.2.3. Scalar Fields with Higgs Potential
As a next case of the investigation of the Jacobi stability of a cosmological scalar field model, we consider that the Universe is filled with a Higgs-like field, with self-interaction potential given by [figure omitted; refer to PDF] where V 0 is a constant and M 2 < 0 is related to the mass of the Higgs boson by the relation m H = V [variant prime] [variant prime] 2 , where V 2 = - M 2 / λ gives the minimum of the potential. The Higgs self-coupling constant λ [approximate] 1 / 8 [91], a value inferred based on the determination of m H from accelerator experiments. By introducing a new dimensionless time variable τ = M t , the basic equations determining the cosmological evolution are given by [figure omitted; refer to PDF] where v 0 = V 0 / M 2 and η = λ / 4 M 2 > 0 , respectively. The time variations of the scale factor of the Higgs field filled Universe and of the scalar field are represented in Figure 6. To numerically integrate the gravitational field equations, we have fixed the value of the constant v 0 to v 0 = 0.005 , and we have adopted as initial conditions the values of the initial conditions a ( 0 ) = a 0 , a ( 0 ) = a 0 H 0 , [varphi] ( 0 ) = [varphi] 0 , and [varphi] ( 0 ) = [varphi] 0 , respectively.
Figure 6: Variation of the scale factor of the Universe (a) and of the scalar field (b) as a function of the dimensionless time τ for the scalar field with Higgs type self-interaction potential, for different values of the parameter η : η = 10 (solid curve), η = 20 (dotted curve), η = 30 (dashed curve), and η = 40 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
The Universe filled with a Higgs type scalar field is expanding, with the scale factor monotonically increasing in time. In the initial phases of the expansion, the scale factor can be approximated by a linearly increasing function of time. The Higgs scalar field [varphi] keeps a constant value in the initial stages of the evolution, followed by a rapid increase associated with an oscillatory behavior with a decreasing amplitude, associated with the energy dissipation. The variation of the Higgs potential is represented in Figure 7.
Figure 7: Variation of the Higgs type potential of the scalar field as a function of the dimensionless time τ for different values of the parameter η : η = 10 (solid curve), η = 20 (dotted curve), η = 30 (dashed curve), and η = 40 (dashed dotted curve).
[figure omitted; refer to PDF]
After an initial phase in which the potential is constant, rapid oscillations with decreasing amplitude do follow, and the potential decays in time. The variation of the KCC invariants 2 κ and χ is represented in Figure 8.
Figure 8: Time variations of the KCC invariants 2 κ = P 1 1 + P 2 2 (a) and χ = P 1 1 P 2 2 - P 2 1 P 2 1 (b) for a Universe filled with a Higgs type scalar field, for different values of the parameter η : η = 10 (solid curve), η = 20 (dotted curve), η = 30 (dashed curve), and η = 40 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
Jacobi stability of the cosmological model with Higgs type scalar field requires the conditions 2 κ < 0 and χ > 0 to be simultaneously satisfied. The Jacobi stable/unstable regions formed during the cosmological expansion of scalar field cosmologies with Higgs self-interacting potential are presented in Figure 9.
Figure 9: The Jacobi stability regions for the scalar field cosmological model with Higgs potential for different values of η : η = 10 (a), η = 20 (b), η = 30 (c), and η = 40 (d). The blue area shows the time intervals with χ > 0 , and the black area shows the time intervals with κ < 0 . Where they overlap, both conditions χ > 0 and κ < 0 are simultaneously satisfied, and the evolution is Jacobi stable.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF]
In Figure 10, we present the time variation of the components of the deviation vector ξ i , i = 1,2 .
Figure 10: Time variations of the deviation vector components ξ 1 (a) and ξ 2 (b) for a Universe filled with a Higgs type scalar field, for different values of the parameter η : η = 10 (solid curve), η = 20 (dotted curve), η = 30 (dashed curve), and η = 40 (dashed dotted curve).
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
4. Jacobi Stability Analysis of the Phantom Quintessence and Tachyonic Cosmological Models
In the present section, we will consider the Jacobi stability analysis of the standard dynamical system formulation of scalar field cosmological models. In this formulation, the Friedmann equations are rewritten in the equivalent form of a three-dimensional first-order dynamical system. Geometrically, we can describe the solutions of a first-order dynamical system as a flow [straight phi] t : D ⊂ R n [arrow right] R n , or, more generally, [straight phi] t : D ⊂ M [arrow right] M , where M is a smooth n -dimensional manifold. The canonical lift of [straight phi] t to the tangent space T M is given by [straight phi] ^ t : T M [arrow right] T M , [straight phi] ^ t ( u ) = ( [straight phi] t ( u ) , [straight phi] t ( u ) ) , where [straight phi] t ( u ) = ∂ [straight phi] ( t , u ) / ∂ t . Therefore, in order to apply the KCC theory to first-order dynamical systems, we must first lift the equations to the tangent bundle. Mathematically, this is equivalent to simply taking the time derivative of the dynamical system.
4.1. First-Order Dynamical System Formulation of Quintessence and Phantom Quintessence Scalar Field Cosmological Models
In the following, we assume that the energy density and pressure of the scalar field can be generally represented as [figure omitted; refer to PDF] where ζ = ± 1 . The case ζ = + 1 corresponds to the quintessence fields, while ζ = - 1 describes the phantom scalar fields. We assume that the Universe is filled with ordinary matter with energy density ρ m and pressure p m and scalar fields. All models are characterized by their dimensionless density parameters Ω m and Ω [varphi] , defined as [figure omitted; refer to PDF] and satisfying the constraint Ω m + Ω [varphi] = 1 . The cosmological equations describing this Universe model take the form [figure omitted; refer to PDF] respectively. In the following, we assume that the matter obeys a linear barotropic equation of state of the form p m = ( γ m - 1 ) ρ m , where γ m is constant and 1 <= γ m <= 2 . From the second Friedmann equation, we immediately obtain the relation [figure omitted; refer to PDF] As basic variables in the first-order dynamical description of cosmological dynamics, we introduce the quantities x and y , defined as [figure omitted; refer to PDF]
In these variables, the density parameter of the scalar field is given by Ω [varphi] = x 2 + y 2 , while the density parameter of the matter can be written as Ω m = 1 - ( x 2 + y 2 ) . Moreover, the energy density and pressure of the scalar field are p [varphi] = 3 H 2 ( ζ x 2 + y 2 ) and p [varphi] = 3 H 2 ( ζ x 2 - y 2 ) , respectively. Instead of the ordinary time variable t , we introduce the new independent variable τ , defined as τ = ln [...] a ( t ) , giving d / d t = H ( d / d τ ) . Then, by taking the time derivative of x and y , we obtain after some simple transformations [53] [figure omitted; refer to PDF] The above system is not a closed system since one more equation involving V [variant prime] / V is still lacking. Therefore, we introduce a third variable z , defined as z = V [variant prime] / V . By taking its derivative with respect to the time τ , we obtain [figure omitted; refer to PDF] where Γ = V V [variant prime] [variant prime] / V [variant prime] 2 . In order to lift this dynamical system to the tangent bundle, we relabel the coordinates as [figure omitted; refer to PDF] Then, by taking the derivatives with respect to τ of the first-order cosmological dynamical system, we obtain the equivalent second-order system: [figure omitted; refer to PDF] The KCC geometric quantities (nonlinear connection, deviation curvature tensor, and geodesic deviation equations) describing the Jacobi stability properties of the phantom quintessence cosmological model are presented in Appendix A.
4.2. The Phantom Scalar Field with Power Law Potential
We assume that the potential of the phantom quintessence scalar field is given by a power law function, so that V ( [varphi] ) = V 0 [varphi] α , where α is a constant. Then, it follows immediately that z = V [variant prime] / V = α / [varphi] , Γ = ( α - 1 ) / α = c o n s t a n t , and Γ [variant prime] ( z ) = 0 . An important cosmological indicator is represented by the parameter of the total equation of state of the matter w , which for dust is defined as [figure omitted; refer to PDF]
The variations of the density parameter of the phantom quintessence field Ω [varphi] and of the parameter of the total equation of state are represented as a function of the dust ( γ m = 1 ) cosmological matter density parameter in Figure 11. The initial conditions used to integrate the cosmological dynamical system are x ( 0 ) = 0.28 , y ( 0 ) = 0.47 , and z ( 0 ) = 0.1 , respectively.
Figure 11: Variation of the density parameter of the phantom quintessence field as a function of the density parameter of the dust cosmological matter (a) and of the parameter w of the total equation of state (b) for different values of α : α = 0.8 (solid curve), α = 1 (dotted curve), α = 1.2 (dashed curve), and α = 1.4 (dashed dotted curve), respectively.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
In this model, the Universe starts its evolution in the large time limit from a matter dominated phase, with Ω m = 1 . During the cosmological expansion the role of the scalar field becomes dominant, and the Universe reaches the present day in a state of accelerated expansion, with Ω [varphi] = 0.75 and Ω m = 0.25 , respectively. On the other hand, the de Sitter phase with w = - 1 is reached only for vanishing ordinary matter density, when the Universe is fully dominated by the phantom quintessence field. It is also important to note that the cosmological evolution is basically independent of the numerical values of the exponent α in the scalar field potential.
The conditions of the Jacobi stability of the phantom quintessence cosmological model in its three-dimensional dynamic system representation are given by the four conditions that must be satisfied by the quantities ( Σ , Φ , Ψ , Ω ) , which are functions of the components of the deviation curvature tensor and which are presented in (26), for different values of α . The time variations of ( Σ , Φ , Ψ , Ω ) are represented in Figure 12.
Figure 12: Time variation of the quantities ( Σ , Φ , Ψ , Ω ) , giving the Jacobi stability of the phantom quintessence field cosmological model with power law potential for different values of α : α = 0.8 (solid curve), α = 1 (dotted curve), α = 1.2 (dashed curve), and α = 1.4 (dashed dotted curve), respectively.
[figure omitted; refer to PDF]
As one can see from Figure 12, the Jacobi stability condition Ω > 0 is not satisfied during the entire time evolution of this cosmological model. Therefore, phantom quintessence scalar field cosmological models with power law potential are Jacobi unstable for all time intervals.
4.3. Jacobi Stability of Tachyon Field Cosmological Models
Tachyon scalar fields have been proposed as possible candidates to explain both inflation and the late acceleration of the Universe [92, 93]. It is also possible for the tachyonic field to trigger the inflationary expansion and at a later time generate a nonrelativistic fluid, which could account for the existence of dark matter. A tachyonic scalar field is described by the Lagrangian L = - V ( [varphi] ) 1 - [varphi] 2 , and it can be described in terms of effective energy density and pressure, given by [figure omitted; refer to PDF] The Friedmann equations are given by [figure omitted; refer to PDF] while the Klein-Gordon type equation satisfied by the scalar field is [figure omitted; refer to PDF]
The parameter of the equation of state of the tachyon field is given by [figure omitted; refer to PDF] By adopting again for the matter the linear barotropic equation of state p m = ( γ m - 1 ) ρ m , 1 <= γ m <= 2 , from the Friedmann equation we obtain first [figure omitted; refer to PDF] To formulate the cosmological evolution equation as a dynamical system, we introduce a set of variables ( τ , x , y , z ) defined as [figure omitted; refer to PDF] In these variables, the density parameters of the scalar field and of the matter are given by [figure omitted; refer to PDF]
Then, by taking the derivative of x and y with respect to τ , we obtain [figure omitted; refer to PDF] where Γ ( z ) = ( 2 / 3 ) V V [variant prime] [variant prime] / V [variant prime] 2 . By lifting the cosmological field equations to the tangent bundle, we obtain [figure omitted; refer to PDF]
The KCC geometric quantities (nonlinear connection, deviation curvature tensor, and geodesic deviation equations) describing the Jacobi stability properties of the tachyon scalar field cosmological model are presented in Appendix B.
4.3.1. Power Law Potential Tachyonic Scalar Field Cosmological Models
In the following, we assume again that the potential of the tachyonic scalar field is power law type, with V ( [varphi] ) = V 0 [varphi] α , where V 0 and α are constants. This choice immediately gives Γ = ( 2 / 3 ) ( α - 1 ) / α . The parametric dependence of the density parameter Ω [varphi] of the tachyonic scalar field with power law potential on the matter energy density Ω m is represented in Figure 13(a). The parametric variation of the matter density parameter Ω m as a function of the parameter of the total equation of state of matter w is shown in Figure 13(b).
Figure 13: Variation of the density parameter of the tachyonic scalar field with power law potential as a function of the density parameter of the dust cosmological matter ( γ m = 1 ) (a) and of the parameter w of the total equation of state (b) for different values of α : α = 0.8 (solid curve), α = 1 (dotted curve), α = 1.2 (dashed curve), and α = 1.4 (dashed dotted curve), respectively.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
In order to numerically integrate the dynamical system corresponding to the tachyon scalar field model, we have used the initial conditions x ( 0 ) = 0.28 , y ( 0 ) = 0.45 , and z ( 0 ) = 0.1 , respectively, and we have assumed that the ordinary matter in the Universe is in the form of zero pressure dust, with γ m = 1 . There is a linear relation between Ω [varphi] and Ω m . The Universe starts its cosmological evolution in a matter dominated phase, with Ω m = 1 and Ω [varphi] = 0 , with decelerating expansion. The energy density of the tachyon field increases in time, and the Universe enters in a de Sitter phase, with the matter density becoming negligibly small. The parameter of the total equation of state satisfies the condition w < 0 during the entire cosmological evolution. It is also interesting to note that the time evolution of this model is basically independent of the numerical values of the exponent α in the power law potential. The time variations of ( Σ , Φ , Ψ , Ω ) describing the Jacobi stability of the tachyon scalar field cosmological model are represented in Figure 14.
Figure 14: Time variation of the quantities ( Σ , Φ , Ψ , Ω ) , giving the Jacobi stability of the tachyon scalar field cosmological model with power law potential for different values of α : α = 0.8 (solid curve), α = 1 (dotted curve), α = 1.2 (dashed curve), and α = 1.4 (dashed dotted curve), respectively.
[figure omitted; refer to PDF]
As one can see from Figure 14, the Jacobi stability condition Ω > 0 is not satisfied during the entire cosmological evolution period of this scalar field model. Therefore, it follows that tachyon scalar field cosmological models with power law potential are Jacobi unstable for all time intervals. This result is independent of the numerical values of the exponent α in the tachyon scalar field potential.
5. Discussions and Final Remarks
In the present paper, we have investigated the Jacobi stability properties of the scalar field cosmological models by using the KCC theory, which represents a powerful mathematical method for the analysis of dynamical systems. Scalar field cosmological models represent a nontrivial testing object for studying nonlinear effects in the framework of general relativity. From a mathematical point of view, the Jacobi (in)stability represents a natural generalization of the (in)stability of the geodesic flow on a differentiable manifold, endowed with a Riemannian or Finslerian type metric to a nonmetric setting . The KCC theory can be applied to scalar field cosmological models that can be formulated mathematically as sets of second-order ordinary nonlinear differential equations. Then, the geometric invariants associated with this system (nonlinear and Berwald connections) and the deviation curvature tensor, as well as its eigenvalues, can be explicitly obtained. The time evolution of the components of the deviation vector can also be obtained by explicitly solving the geodesic deviation equations.
The Jacobi stability, and its theoretical foundation, the KCC theory, offers an alternative approach to the "classical" Lyapunov approach, by investigating the deviations of the entire trajectory of the cosmological evolution equations with respect to the nearby ones under the effects of a small perturbation. In the framework of general relativity, we may call the applications of the KCC theory to the study of the gravitational fields as a "second geometrization," in which already geometric quantities are supplemented by additional geometric structures. Hence, general relativistic cosmological models can be described in geometric terms originating from their dynamical system structure, with these new geometric structures fully determined by the underlying Riemannian geometry and the physical properties of the scalar fields (their self-interaction potential). The stability properties of the perturbations of a given trajectory describing the cosmological evolution are determined by the properties of the curvature deviation tensor, a geometric quantity constructed from the connections (nonlinear and Berwald) associated with the dynamical system describing the cosmological evolution. It is important to note that the KCC theory can be directly applied to systems of second-order differential equations, which can be interpreted geometrically as the paths (or geodesics) associated with a semispray. In investigating the Jacobi stability of cosmological models, we have followed two approaches. Since the cosmological evolution equations (the Friedmann equations) are second-order differential equations, the KCC theory can be naturally and directly applied to study the stability of the cosmic evolution. As a first step, one obtains the two nonzero components of the nonlinear connection, with the N 2 1 component depending on the product of the scale factor and of the time derivative of the field, while the N 2 2 component depends on the energy density of the scalar field, as well as of the scalar field potential. After obtaining the components of the deviation curvature tensor, we have formulated the general condition of the stability of the scalar field cosmological models, which is determined by two inequalities involving the second and the first derivative of the scalar field potential, the energy density of the field, and the time variation of the scalar field itself.
As an application of the developed formalism, we have considered two scalar field models, both being relevant for the study of both the early and the late stages of the cosmological evolution. The first case we did consider is the scalar field with exponential potential. We have studied in detail the KCC geometric properties of this model. It turns out that the Jacobi stability condition which can be expressed in terms of the components of the deviation curvature tensor is not satisfied during the cosmological evolution and that the Universe described by the exponential potential scalar field is in a Jacobi unstable state. This result is independent of the numerical values of the parameter λ , describing the properties of the potential, and can also be inferred from the behavior of the components of the deviation vector ξ i , with ξ 1 diverging exponentially in time. As a second case, we have considered the case of the Higgs type potential. For this potential, the KCC geometric quantities show a complex behavior. After a period in which the scalar field and the potential are almost constant, the field starts to oscillate, with the amplitude of the oscillations decreasing in time. This behavior of the Higgs field is also reflected in the behavior of the components of the deviation curvature tensor, which are also some oscillating functions. The Jacobi stability of this cosmological model strongly depends on the numerical value of the parameter η = λ / 4 M 2 . For small values of η , the Universe evolves between successive Jacobi stable and unstable states. With the increase of the numerical value of η , the time intervals in which the Universe is Jacobi stable decrease quickly, and for large values of η the Universe is in a Jacobi unstable state during its entire cosmological evolution.
As a second approach for the study of the Jacobi stability of scalar field cosmologies, we have considered the first-order dynamical system formulation of scalar field evolution equations. In this approach, by introducing a new set of variables, expressed in terms of the square root of the potential, the time derivative of the scalar field, and the Hubble function, respectively, the Friedmann equations in the presence of scalar fields can be reformulated as a first-order dynamical system, consisting of three highly nonlinear ordinary differential equations. In order to apply the KCC theory, this dynamical system must be lifted to the tangent bundle and formulated as a second-order differential system. We have analyzed, by using this approach, two specific scalar field models, the phantom quintessence and the tachyon scalar field with power law potentials. It turns out that for this choice of the potential both scalar field models are Jacobi unstable. The power law potential gives a very simple form for the function Γ ( z ) , which takes constant values during the cosmological evolution. This situation is similar to the case of the exponential potential and leads to a significant simplification of the mathematical formalism.
We have started our study of the applications of the KCC theory to cosmological problems with the investigation of the standard matter dominated cosmological models in their dynamical system formulation. We have studied the Jacobi stability of the critical points of different models, and we have shown that they are Jacobi unstable. A full comparison between the Jacobi and Lyapunov properties of the critical points for second-order systems was given in [59, 60], and hence we will not discuss in detail this relation. However, this study of the critical points of matter dominated cosmological models also shows the fundamental differences between the Lyapunov stability and KCC theories: while Lyapunov stability is mostly restricted to the study of critical points, the KCC theory has the potential to investigate the deviation of the full trajectory during the entire period of the cosmological evolution. Therefore, we can consider a Lyapunov stability analysis of steady states (called the linear analysis) and a "Lyapunov type" stability analysis of the whole trajectory (the KCC or Jacobi stability analysis), and these two methods are complementary but distinct to each other.
The KCC theory also introduces the first set of KCC invariants [...] i , i = 1 , ... , n , giving the contravariant KCC derivative of the vector field y i . The first KCC invariant can be interpreted as an external force. We did not study in detail the time evolution of the first KCC invariant, since its properties are not directly related to the stability issues that were our main points of interest.
In the present paper, we have performed a stability analysis of the scalar field cosmological models, in which we have considered a description of the deviations of the whole trajectories of the differential system describing the cosmological dynamics, and we have provided some basic theoretical and computational tools for this study. Further investigations of the Jacobi stability properties of cosmological models may provide some methods for discriminating between different evolutionary scenarios, as well as for better understanding of some other fundamental processes, like, for example, structure formation that played an essential role in the evolution of our Universe.
[1] A. G. Riess, A. V. Filippenko, P. Challis, "Observational evidence from supernovae for an accelerating universe and a cosmological constant," The Astronomical Journal , vol. 116, no. 3, 1998.
[2] S. Perlmutter, G. Aldering, G. Goldhaber, "Measurements of Ω and Λ from 42 high-redshift supernovae," The Astrophysical Journal , vol. 517, no. 2, pp. 565, 1999.
[3] R. A. Knop, G. Aldering, R. Amanullah, "New constraints on ΩM , ΩΛ , and w from an independent set of 11 high-redshift supernovae observed with the hubble space telescope," The Astrophysical Journal , vol. 598, no. 1, pp. 102-137, 2003.
[4] R. Amanullah, C. Lidman, D. Rubin, "Spectra and Hubble Space Telescope light curves of six type Ia supernovae at 0.511 < z < 1.12 and the Union2 compilation," The Astrophysical Journal , vol. 716, no. 1, pp. 712, 2010.
[5] D. H. Weinberg, M. J. Mortonson, D. J. Eisenstein, C. Hirata, A. G. Riess, E. Rozo, "Observational probes of cosmic acceleration," Physics Reports , vol. 530, no. 2, pp. 87-255, 2013.
[6] P. A. R. Ade, N. Aghanim, M. Arnaud, "Planck 2015 results. XVIII. Background geometry & topology," Astronomy & Astrophysics , vol. 594, article A18, 2015.
[7] P. J. Peebles, B. Ratra, "The cosmological constant and dark energy," Reviews of Modern Physics , vol. 75, no. 2, pp. 559-606, 2003.
[8] T. Padmanabhan, "Cosmological constant--the weight of the vacuum," Physics Reports , vol. 380, no. 5-6, pp. 235-320, 2003.
[9] S. Nojiri, S. D. Odintsov, "Unified cosmic history in modified gravity: from F(R) theory to Lorentz non-invariant models," Physics Reports , vol. 505, no. 2-4, pp. 59-144, 2011.
[10] M. Li, X.-D. Li, S. Wang, Y. Wang, "Dark energy: a brief review," Frontiers of Physics , vol. 8, no. 6, pp. 828-846, 2013.
[11] M. J. Mortonson, D. H. Weinberg, M. White, "Dark Energy: a short review," https://arxiv.org/abs/1401.0046
[12] R. R. Caldwell, R. Dave, P. J. Steinhardt, "Cosmological imprint of an energy component with general equation of state," Physical Review Letters , vol. 80, no. 8, pp. 1582-1585, 1998.
[13] S. Tsujikawa, "Quintessence: a review," Classical and Quantum Gravity , vol. 30, no. 21, 2013.
[14] L. P. Chimento, A. S. Jakubi, D. Pavón, "Enlarged quintessence cosmology," Physical Review D , vol. 62, no. 6, 2000.
[15] T. Chiba, T. Okabe, M. Yamaguchi, "Kinetically driven quintessence," Physical Review D , vol. 62, no. 2, 2000.
[16] C. Armendariz-Picon, V. Mukhanov, P. J. Steinhardt, "Dynamical solution to the problem of a small cosmological constant and late-time cosmic acceleration," Physical Review Letters , vol. 85, no. 21, pp. 4438-4441, 2000.
[17] C. Armendariz-Picon, V. F. Mukhanov, P. J. Steinhardt, "Essentials of k-essence," Physical Review D , vol. 63, no. 10, 2001.
[18] N. Arkani-Hamed, H.-C. Cheng, M. A. Luty, S. Mukohyama, "Ghost condensation and a consistent infrared modification of gravity," Journal of High Energy Physics , vol. 2004, article 074, 2004.
[19] F. Piazza, S. Tsujikawa, "Dilatonic ghost condensate as dark energy," Journal of Cosmology and Astroparticle Physics , vol. 2004, no. 7, article 004, 2004.
[20] R. R. Caldwell, "A phantom menace? Cosmological consequences of a dark energy component with super-negative equation of state," Physics Letters B , vol. 545, no. 1-2, pp. 23-29, 2002.
[21] S. M. Carroll, M. Hoffman, M. Trodden, "Can the dark energy equation-of-state parameter w be less than -1?," Physical Review D , vol. 68, no. 2, 2003.
[22] P. Singh, M. Sami, N. Dadhich, "Cosmological dynamics of a phantom field," Physical Review D , vol. 68, no. 2, 2003.
[23] M. Sami, A. Toporensky, "Phantom field and the fate of the universe," Modern Physics Letters A , vol. 19, no. 20, pp. 1509-1517, 2004.
[24] J. M. Cline, S. Jeon, G. D. Moore, "The phantom menaced: constraints on low-energy effective ghosts," Physical Review D , vol. 70, no. 4, 2004.
[25] E. Elizalde, S. Nojiri, S. D. Odintsov, "Late-time cosmology in a (phantom) scalar-tensor theory: dark energy and the cosmic speed-up," Physical Review D , vol. 70, no. 4, 2004.
[26] E. Elizalde, S. Nojiri, S. D. Odintsov, D. Sáez-Gómez, V. Faraoni, "Reconstructing the universe history, from inflation to acceleration, with phantom and canonical scalar fields," Physical Review D , vol. 77, no. 10, 2008.
[27] A. Y. Kamenshchik, "Quantum cosmology and late-time singularities," Classical and Quantum Gravity , vol. 30, no. 17, 2013.
[28] U. Alam, V. Sahni, T. D. Saini, A. A. Starobinsky, "Is there supernova evidence for dark energy metamorphosis?," Monthly Notices of the Royal Astronomical Society , vol. 354, no. 1, pp. 275-291, 2004.
[29] B. A. Bassett, S. Tsujikawa, D. Wands, "Inflation dynamics and reheating," Reviews of Modern Physics , vol. 78, no. 2, pp. 537-589, 2006.
[30] A. Maleknejad, M. M. Sheikh-Jabbari, J. Soda, "Gauge fields and inflation," Physics Reports , vol. 528, no. 4, pp. 161-261, 2013.
[31] A. D. Linde Particle Physics and Inflationary Cosmology , Harwood Academic, Reading, UK, 1990.
[32] A. R. Liddle, "Inflation and the cosmic microwave background," Physics Reports , vol. 307, no. 1-4, pp. 53-60, 1998.
[33] E. J. Copeland, E. W. Kolb, A. R. Liddle, J. E. Lidsey, "Reconstructing the inflaton potential: perturbative reconstruction to second order," Physical Review D , vol. 49, no. 4, pp. 1840-1844, 1994.
[34] F. E. Schunk, E. W. Mielke, "Method of generating exact inflationary solutions," Physical Review D , vol. 50, no. 8, pp. 4794-4806, 1994.
[35] A. M. Mancho, D. Small, S. Wiggins, "A tutorial on dynamical systems concepts applied to Lagrangian transport in oceanic flows defined as finite time data sets: theoretical and computational issues," Physics Reports , vol. 437, no. 3-4, pp. 55-124, 2006.
[36] G. Boffetta, M. Cencini, M. Falcioni, A. Vulpiani, "Predictability: a way to characterize complexity," Physics Reports , vol. 356, no. 6, pp. 367-474, 2002.
[37] J. Wainwright, G. F. R. Ellis Dynamical Systems in Cosmology , Cambridge University Press, Cambridge, UK, 1997.
[38] A. A. Coley Dynamical Systems and Cosmology , Kluwer Academic Publishers, Dordrecht, The Netherlands, 2003.
[39] J. Miritzis, S. Cotsakis, E. Papantonopoulos, "Introduction to cosmological dynamical systems," Cosmological Crossroads: An Advanced Course in Mathematical, Physical and String Cosmology , vol. 592, of Lecture Notes in Physics, pp. 111-143, Springer, Berlin, Germany, 2002.
[40] H. van Elst, C. Uggla, J. Wainwright, "Dynamical systems approach to G2 cosmology," Classical and Quantum Gravity , vol. 19, no. 1, pp. 51-82, 2002.
[41] J. Carot, M. M. Collinge, "Scalar field cosmologies: a dynamical systems study," Classical and Quantum Gravity , vol. 20, no. 4, pp. 707-728, 2003.
[42] M. Szydlowski, A. Krawiec, "Dynamical system approach to cosmological models with a varying speed of light," Physical Review D , vol. 68, no. 6, 2003.
[43] L. Lara, M. Castagnino, "Brane cosmology as a dynamical system," Physical Review D , vol. 70, no. 4, 2004.
[44] L. Arturo Urena-Lopez, "Scalar phantom energy as a cosmological dynamical system," Journal of Cosmology and Astroparticle Physics , vol. 2005, no. 9, pp. 13, 2005.
[45] N. Goheer, J. A. Leach, P. K. S. Dunsby, "Dynamical systems analysis of anisotropic cosmologies in R n -gravity," Classical and Quantum Gravity , vol. 24, no. 22, pp. 5689-5708, 2007.
[46] A. Bonanno, S. Carloni, "Dynamical system analysis of cosmologies with running cosmological constant from quantum einstein gravity," New Journal of Physics , vol. 14, 2012.
[47] D. J. E. Marsh, E. R. M. Tarrant, E. J. Copeland, P. G. Ferreira, "Cosmology of axions and moduli: a dynamical systems approach," Physical Review D , vol. 86, no. 2, 2012.
[48] G. Leon, C. R. Fadragas Cosmological Dynamical Systems: And Their Applications , LAP Lambert Academic Publishing, 2012.
[49] S. Carloni, T. Koivisto, F. S. N. Lobo, "Dynamical system analysis of hybrid metric-Palatini cosmologies," Physical Review D , vol. 92, 2015.
[50] N. Mahata, S. Chakraborty, "A dynamical system analysis of three fluid cosmological model," http://arxiv.org/abs/1512.07017
[51] A. Stachowski, M. Szydlowski, "Dynamical system approach to running Λ cosmological models," http://arxiv.org/abs/1601.05668
[52] C. G. Boehmer, N. Chan, "Dynamical systems in cosmology," https://arxiv.org/abs/1409.5585
[53] R. García-Salcedo, T. Gonzalez, F. A. Horta-Rangel, I. Quiros, D. Sanchez-Guzmán, "Introduction to the application of dynamical systems theory in the study of the dynamics of cosmological models of dark energy," European Journal of Physics , vol. 36, no. 2, 2015.
[54] D. D. Kosambi, "Parallelism and path-spaces," Mathematische Zeitschrift , vol. 37, no. 1, pp. 608-618, 1933.
[55] E. Cartan, D. D. Kosambi, "Observations sur le mémoire précédent," Mathematische Zeitschrift , vol. 37, no. 1, pp. 619-622, 1933.
[56] S. S. Chern, "Sur la géométrie d'un système d'équations différentielles du second order," Bulletin des Sciences Mathematiques , vol. 63, pp. 206-212, 1939.
[57] C. G. Böhmer, T. Harko, S. V. Sabau, "Jacobi stability analysis of dynamical systems--applications in gravitation and cosmology," Advances in Theoretical and Mathematical Physics , vol. 16, no. 4, pp. 1145-1196, 2012.
[58] P. L. Antonelli Handbook of Finsler Geometry , vol. 1, Kluwer Academic, Dordrecht, The Netherlands, 2003.
[59] S. V. Sabau, "Some remarks on Jacobi stability," Nonlinear Analysis: Theory, Methods & Applications , vol. 63, no. 5-7, pp. e143-e153, 2005.
[60] V. S. Sabau, "Systems biology and deviation curvature tensor," Nonlinear Analysis. Real World Applications , vol. 6, no. 3, pp. 563-587, 2005.
[61] P. L. Antonelli, "On y -Berwald connections and Hutchinson's ecology of social interactions," Tensor, New Series , vol. 52, no. 1, pp. 27-36, 1993.
[62] T. Yajima, H. Nagahama, "K{CC}-theory and geometry of the Rikitake system," Journal of Physics A , vol. 40, no. 11, pp. 2755-2772, 2007.
[63] T. Yajima, H. Nagahama, "KCC-theory and geometry of the Rikitake system," Journal of Physics A: Mathematical and Theoretical , vol. 40, no. 11, pp. 2755-2772, 2007.
[64] T. Harko, V. S. Sabau, "Jacobi stability of the vacuum in the static spherically symmetric brane world models," Physical Review D , vol. 77, no. 10, 2008.
[65] C. G. Böhmer, T. Harko, "Nonlinear stability analysis of the emden-fowler equation," Journal of Nonlinear Mathematical Physics , vol. 17, no. 4, pp. 503-516, 2010.
[66] T. Yajima, H. Nagahama Acta Mathematica Academiae Paedagogicae Nyíregyháziensis , vol. 24, 2008.
[67] T. Yajima, H. Nagahama, "Tangent bundle viewpoint of the Lorenz system and its chaotic behavior," Physics Letters A , vol. 374, no. 11-12, pp. 1315-1319, 2010.
[68] H. Abolghasem, "Liapunov stability versus Jacobi stability," Journal of Dynamical Systems and Geometric Theories , vol. 10, no. 1, pp. 13-32, 2012.
[69] H. Abolghasem, "Jacobi stability of circular orbits in a central force," Journal of Dynamical Systems and Geometric Theories , vol. 10, no. 2, pp. 197-214, 2012.
[70] H. Abolghasem, "Stability of circular orbits in Schwarzschild spacetime," International Journal of Differential Equations and Applications , vol. 12, no. 3, pp. 131-147, 2013.
[71] H. Abolghasem, "Jacobi stability of hamiltonian systems," International Journal of Pure and Applied Mathematics , vol. 87, no. 1, pp. 181-194, 2013.
[72] T. Harko, C. Y. Ho, C. S. Leung, S. Yip, "Jacobi stability analysis of the Lorenz system," International Journal of Geometric Methods in Modern Physics , vol. 12, no. 7, 2015.
[73] T. Harko, P. Pantaragphong, S. Sabau, "A new perspective on the Kosambi-Cartan-Chern theory, and its applications," https://arxiv.org/abs/1509.00168
[74] M. J. Lake, T. Harko, "Dynamical behavior and Jacobi stability analysis of wound strings," The European Physical Journal C , vol. 76, article 311, 2016.
[75] M. K. Gupta, C. K. Yadav, "Jacobi stability analysis of Rikitake system," International Journal of Geometric Methods in Modern Physics , vol. 13, no. 7, pp. 1650098-1650818, 2016.
[76] T. Yajima, K. Yamasaki, "Jacobi stability for dynamical systems of two-dimensional second-order differential equations and application to overhead crane system," International Journal of Geometric Methods in Modern Physics , vol. 13, no. 4, 2016.
[77] T. Harko, P. Pantaragphong, S. V. Sabau, "Kosambi-Cartan-Chern (KCC) theory for higher-order dynamical systems," International Journal of Geometric Methods in Modern Physics , vol. 13, no. 2, 2016.
[78] M. Pettini, "Geometrical hints for a nonperturbative approach to Hamiltonian dynamics," Physical Review E , vol. 47, no. 2, pp. 828-850, 1993.
[79] H. E. Kandrup, "Geometric interpretation of chaos in two-dimensional Hamiltonian systems," Physical Review E , vol. 56, no. 3, pp. 2722-2732, 1997.
[80] M. Di Bari, D. Boccaletti, P. Cipriani, G. Pucacco, "Dynamical behavior of Lagrangian systems on Finsler manifolds," Physical Review E , vol. 55, no. 6, pp. 6448-6458, 1997.
[81] P. Cipriani, M. Di Bari, "Finsler geometric local indicator of chaos for single orbits in the hénon-heiles hamiltonian," Physical Review Letters , vol. 81, no. 25, pp. 5532-5535, 1998.
[82] M. Di Bari, P. Cipriani, "Geometry and chaos on Riemann and Finsler manifolds," Planetary and Space Science , vol. 46, no. 11-12, pp. 1543-1555, 1998.
[83] L. Casetti, M. Pettini, E. G. D. Cohen, "Geometric approach to Hamiltonian dynamics and statistical mechanics," Physics Reports , vol. 337, no. 3, pp. 237-341, 2000.
[84] G. Ciraolo, M. Pettini, "Geometry of chaos in models of stellar dynamics," Celestial Mechanics and Dynamical Astronomy , vol. 83, no. 1-4, pp. 171-190, 2002.
[85] R. Punzi, M. N. R. Wohlfarth, "Geometry and stability of dynamical systems," Physical Review E , vol. 79, no. 4, 2009.
[86] R. Miron, D. Hrimiuc, H. Shimada, S. V. Sabau The geometry of Hamilton and Lagrange spaces , vol. 118, of Fundamental Theories of Physics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2001.
[87] P. L. Antonelli, I. Bucataru, "New results about the geometric invariants in KCC-theory," Analele Stiintifice ale Universitatii 'Al. I. Cuza' Iasi, Mathematics, New Series , vol. 47, pp. 405-420, 2001.
[88] P. L. Antonelli, R. S. Ingarden, M. Matsumoto The Theory of Sprays and Finsler Spaces with Applications in Physics and Biology , vol. 58, of Fundamental Theories of Physics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1993.
[89] I. Bucataru, O. Constantinescu, M. F. Dahl, "A geometric setting for systems of ordinary differential equations," International Journal of Geometric Methods in Modern Physics , vol. 8, no. 6, pp. 1291-1327, 2011.
[90] Q. I. Rahman, G. Schmeisser Analytic Theory of Polynomials. London Mathematical Society Monographs , of New Series 26, Oxford University Press, Oxford, UK, 2002.
[91] G. Aad, B. Abbott, J. Abdallah, "Combined measurement of the Higgs boson mass in pp collisions at s = 7 and 8 TeV with the ATLAS and CMS experiments," Physical Review Letters , vol. 114, 2015.
[92] M. Sami, P. Chingangbam, T. Qureshi, "Aspects of tachyonic inflation with an exponential potential," Physical Review D , vol. 66, no. 4, 2002.
[93] G. Calcagni, A. R. Liddle, "Tachyon dark energy models: dynamics and constraints," Physical Review D , vol. 74, no. 4, 2006.
Appendix
A. The Nonlinear Connection, the Curvature Deviation Tensor, and the Geodesic Deviation Equations for the Phantom Quintessence Scalar Field Cosmological Model
The coefficients of the nonlinear connection associated with the dynamical system describing the phantom quintessence cosmological model are given by [figure omitted; refer to PDF]
The first KCC invariants [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , are given by [figure omitted; refer to PDF]
All Berwald connection components are zero here. The components of the curvature deviation tensor are [figure omitted; refer to PDF]
The geodesic deviation equations are obtained in the form [figure omitted; refer to PDF]
B. The Nonlinear Connection, the Curvature Deviation Tensor, and the Geodesic Deviation Equations for the Tachyon Scalar Field Cosmological Model
The coefficients of the nonlinear connection associated with the dynamical system describing the tachyon scalar field cosmological model are given by [figure omitted; refer to PDF]
The first KCC invariants [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , are [figure omitted; refer to PDF] All Berwald connection components are zero here. The components of the curvature deviation tensor are [figure omitted; refer to PDF]
The geodesic deviation equations are obtained in the form [figure omitted; refer to PDF]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2016 Bogdan Danila 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 SCOAP3 .
Abstract
We study the stability of the cosmological scalar field models by using the Jacobi stability analysis, or the Kosambi-Cartan-Chern (KCC) theory. In this approach, we describe the time evolution of the scalar field cosmologies in geometric terms, by performing a "second geometrization" and considering them as paths of a semispray. By introducing a nonlinear connection and a Berwald-type connection associated with the Friedmann and Klein-Gordon equations, five geometrical invariants can be constructed, with the second invariant giving the Jacobi stability of the cosmological model. We obtain all the relevant geometric quantities, and we formulate the condition for Jacobi stability in scalar field cosmologies. We consider the Jacobi stability properties of the scalar fields with exponential and Higgs type potential. The Universe dominated by a scalar field exponential potential is in Jacobi unstable state, while the cosmological evolution in the presence of Higgs fields has alternating stable and unstable phases. We also investigate the stability of the phantom quintessence and tachyonic scalar field models, by lifting the first-order system to the tangent bundle. It turns out that in the presence of a power law potential both of these models are Jacobi unstable during the entire cosmological evolution.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer