Next Article in Journal
A Novel Three-Step Numerical Solver for Physical Models under Fractal Behavior
Next Article in Special Issue
Klein-Gordon Theory in Noncommutative Phase Space
Previous Article in Journal
Solvent Effects and Metal Ion Recognition in Several Azulenyl-Vinyl-Oxazolones
Previous Article in Special Issue
Quantum Geometry of Spacetime and Quantum Equilibrium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Jacobi and Lyapunov Stability Analysis of Circular Geodesics around a Spherically Symmetric Dilaton Black Hole

1
Faculty of Mathematics and Computer Sciences, Babes-Bolyai University, 400084 Cluj-Napoca, Romania
2
Department of Mathematics and Computer Sciences, Babes-Bolyai University, 400084 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Department of Physics, Babes-Bolyai University, 400084 Cluj-Napoca, Romania.
Symmetry 2023, 15(2), 329; https://doi.org/10.3390/sym15020329
Submission received: 9 December 2022 / Revised: 4 January 2023 / Accepted: 17 January 2023 / Published: 24 January 2023

Abstract

:
We analyze the stability of the geodesic curves in the geometry of the Gibbons–Maeda–Garfinkle–Horowitz–Strominger black hole, describing the space time of a charged black hole in the low energy limit of the string theory. The stability analysis is performed by using both the linear (Lyapunov) stability method, as well as the notion of Jacobi stability, based on the Kosambi–Cartan–Chern theory. Brief reviews of the two stability methods are also presented. After obtaining the geodesic equations in spherical symmetry, we reformulate them as a two-dimensional dynamic system. The Jacobi stability analysis of the geodesic equations is performed by considering the important geometric invariants that can be used for the description of this system (the nonlinear and the Berwald connections), as well as the deviation curvature tensor, respectively. The characteristic values of the deviation curvature tensor are specifically calculated, as given by the second derivative of effective potential of the geodesic motion. The Lyapunov stability analysis leads to the same results. Hence, we can conclude that, in the particular case of the geodesic motion on circular orbits in the Gibbons–Maeda–Garfinkle–Horowitz–Strominger, the Lyapunov and the Jacobi stability analysis gives equivalent results.

1. Introduction

The analysis of the global stability of the solutions of the systems of strongly nonlinear, first or second order differential equations, describing the temporal evolution of complicated dynamical systems, is usually performed in the framework of a rigorous mathematical approach, with the help of the well known, and mathematically extensively investigated Lyapunov stability theory. In this commonly accepted approach to the problem of the stability of the solutions of differential equations, the basic mathematical quantities to be studied are the Lyapunov exponents. The Lyapunov exponents determine the exponential departures of the trajectories, obtained as solutions of the differential equations describing a given dynamical system, as compared to a standard trajectory, taken as reference [1,2]. However, when one attempts to obtain the Lyapunov exponents, one must understand that generally it is a very complicated task to find them in an exact analytical form. For that reason, in order to evaluate the basic Lyapunov exponents, and their properties, one should apply mostly complicated numerical methods. Currently, a large number of powerful numerical methods have been established for their calculation. Such numerical methods are extensively used for the mathematical, as well as physical description of the temporal evolution near the critical points of the differential equations describing dynamical systems [3,4].
The mathematical approach based on the Lyapunov linear stability analysis is well developed from a mathematical point of view. Moreover, it provides an intuitive and clear comprehension of some of the stability properties of the systems of differential equations. Nevertheless, to obtain a more profound understanding of the time evolution and general properties of the natural and mathematical structures, alternative approaches for the investigation of the stability problems must also be proposed, examined, and fully investigated. Hence, after a new method for the study of the stability is developed, one could contrast exhaustively the results, predictions, and the possible shortcomings of the newly introduced method with the consequences and conclusions gained by adopting the linear Lyapunov stability analysis of the considered system of differential equations.
One of the alternative approaches to the problem of stability that could provide valuable information for the investigation of the stability of the systems of ordinary second order differential equations is given by what one could designate as the geometro-dynamical approach. Historically, one of the first illustrations of such a different mathematical investigation of the stability problem of the systems of differential equations is the Kosambi–Cartan–Chern (KCC) theory. In its early formulation, the KCC theory was developed in the major investigations of Kosambi [5], Cartan [6] and Chern [7], respectively, which led to a first rigorous mathematical investigation of the geometric properties of dynamical systems. From a strictly mathematical perspective, the KCC theory is strongly influenced by the geometry of the Finsler spaces, which largely represents, and significantly contributed, to its theoretical justification. The KCC theory is built up by using the basic conjecture according to which there is a mathematical analogy between autonomous or non-autonomous dynamical system, formulated in terms of second order differential equations, and the equations of the geodesics in a Finsler space. The Finsler type geodesic equations can always be associated with a given dynamical system, described by second order, usually strongly nonlinear systems of differential equations (for an in depth presentation of the KCC theory, and of its applications, see [8]).
The KCC theory represents essentially a geometric way of thinking about the variational differential equations describing the divergence/convergence of a bunch of trajectories of a dynamical system, or of a system of differential equations, as compared to the neighbouring ones [9]. The KCC theory proposes a geometrical type characterization of the systems of second order differential equations, considered as geodesic curves. Within the framework of this description, one can construct for each system of differential equations two connections, which are essentially geometric quantities. The first connection is the nonlinear connection N j i , while a Berwald type connection G j k i is also considered for the description of the dynamical system. With the use of these two connections, which can be generally defined, five important geometrical invariants can be constructed rigorously. Of these five geometrical invariants, the most relevant, from a mathematical and physical point of view, is the second invariant P j i , which is named the deviation curvature tensor. From the general perspective of the scientific, engineering and mathematical applications, its importance is given by its main property of determining the so-called Jacobi stability of the considered system of second order nonlinear differential equation [8,9,10,11,12]. Various engineering, physical, chemical, biochemical, or medical systems have been thoroughly investigated with the help of the KCC stability theory [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27].
The KCC theory has also found applications in the study of the gravitational phenomena. Thus, in [14], the static, spherically symmetric structure equations of the static vacuum in the brane world models were analyzed, from the point of view of their stability, by applying both the Jacobi stability analysis, and the linear (Lyapunov) stability analysis. It was shown that the trajectories that are unstable on the static brane with spherical symmetry behave chaotically, which implies that, after a bunch of particles travel a restricted range of a radial distance, it would not be possible to differentiate the trajectories that were extremely close to each other at an initial time, and at an initial point. Thus, the KCC theory together with the Jacobi stability analysis represent a very powerful method for giving some important constraints on the physical properties of the vacuum on the four-dimensional brane.
The stability of the radial solutions of the Lane–Emden semilinear elliptic equation Δ u + u n = 0 , with the initial conditions u ( 0 ) = 1 and u ( 0 ) = 0 , respectively, were studied on the positive real line in [15], by using the Lyapunov standard linear stability analysis, the Jacobi stability approach, and the Lyapunov function method, respectively. By using the KCC theory, one can obtain for the stability of a polytropic star the criterion E i / E g < n ρ ( r ) / ρ ¯ , where E i is the internal energy, E g is the gravitational energy, ρ ( r ) is the density of the stellar matter, and ρ ¯ ( r ) is the mean mass density, respectively. An investigation of the stability properties of the inflationary cosmological models in the presence of scalar fields, by using the Jacobi stability analysis, was performed by considering a “second geometrization” of the models, and interpreting them as paths of a semispray, in [23]. The KCC stability properties of the cosmological models in the presence of scalar fields with exponential and Higgs type potentials were considered in detail. The KCC theory was used to investigate the Jacobi type stability/instability of the string equations in [25]. Moreover, by using this approach, precise bounds on the geometrical and physical parameters that guarantee dynamical stability of the windings were determined. It was found that, for the same initial conditions, and in higher dimensions, the topology and the curvature of the internal space have significant influences on the microscopic behavior of the string. On the other hand, it turns out, surprisingly, that the macroscopic behavior of the string is not sensitive to the details of the physical motion in the compact space. The stability of the circular restricted three-body problem, which considers the motion of a particle with a very small mass due to the gravitational attraction of two massive stellar type objects moving on circular orbits about their common center of mass, was analyzed in [27], by using the KCC) theory. It was found that, from the geometric perspective of the KCC theory, the five Lagrangian equilibrium points of the restricted three body problem are all unstable.
Black holes are important observational astrophysical objects, as well as fundamental testing grounds for the theories of gravitation. In particular, the motion of particles around the central black hole can give essential information on the physical processes taking place in the cosmic environments. There are many known black hole solutions, obtained in the framework of the different gravitational theories. Exact black hole type solutions of string theory play an important role for the confrontation of the predictions of the theory with observations. Static, spherically symmetric charged black hole solutions in the low energy limit of string theory were found in [28,29], respectively. These solutions are characterized by three independent parameters: their gravitational mass, electric charge, and the asymptotic value of the scalar dilaton field, respectively, whose presence has important physical consequences. A particular class of solutions, the extremely charged “black holes”, represent, from a geometric point of view, geodesically complete spacetimes, without event horizons and singularities.
It is the goal of the present paper to study the stability properties of the geodesic trajectories in the charged dilatonic solution of the low energy limit of string theory, as obtained in [28,29], respectively. After obtaining the geodesic equations of motion, and the expression of the effective potential, the stability of the trajectories is analyzed by using both the linear Lyapunov and the KCC theory based Jacobi stability approaches. It turns out that, in the case of the string theory inspired dilatonic black hole solution, the predictions of both stability methods coincide.
The present paper is organized as follows: We briefly review the Lyapunov and the Jacobi stability approaches in Section 2. The charged black hole solution of the dilatonic low energy string theory is written down in Section 3, where the equations of the geodesic lines are also obtained. The stability of the trajectories of the particles moving around the black hole is investigated, by using both the Lyapunov and the Jacobi approaches, in Section 4. We discuss and conclude our results in Section 5.

2. Brief Review of Lyapunov Stability and of the KCC Theory

In the present section, we quickly introduce, in a concise but rigorous way, the fundamental ideas, the basic concepts, and the results of the linear Lyapunov stability theory, and of the KCC theory, respectively. For in depth discussions of the mathematical aspects of the linear (Lyapunov) stability theory, and of its applications in astrophysics and cosmology, see [30,31,32,33].
Furthermore, we present the notations of the geometrical and physical quantities used in the present investigation, and we introduce the basic definitions of the relevant geometric objects met in KCC theory (for a detailed presentation, we refer the reader to [8,9]).

2.1. Linear Stability of the Systems of Ordinary Differential Equations

We start by concisely presenting, mostly using the approach introduced in [34], some key consequences of the Lyapunov stability analysis of arbitrary dynamical systems, expressed by general systems of first order Ordinary Differential Equations (ODEs). In the beginning of our presentation, we would like to mention that the stability of a system of first ordinary differential equations is determined, in general, by the roots of its characteristic polynomial. To prove this property, we consider the system of autonomous first order ordinary differential equations [34],
d x 1 d t = f 1 x 1 , x 2 , , x n , d x 2 d t = f 2 x 1 , x 2 , , x n , , d x n 1 d t = f n 1 x 1 , x 2 , , x n , d x n d t = f n x 1 , x 2 , , x n ,
where f 1 , f 2 , , f n are, by definition, n smooth functions, possessing derivatives of all orders in their domain of definition. Now, let us linearize the system (1) about one of its steady states (fixed point, or equilibrium point) x 0 1 , x 0 2 , , x 0 n , by associating with the arbitrary system (1) the linear system
d x 1 d t d x 2 d t d x n d t = A x 1 x 2 x n ,
where the Jacobian matrix A of the system (1), defined as A = J f 1 , f 2 , , f n x 0 1 , x 0 2 , , x 0 n , is evaluated at the steady state,
A = f 1 x 1 f 1 x 2 f 1 x n f n x 1 f n x 2 f n x n x 0 1 , x 0 2 , , x 0 n .
The solutions of Equation (2) can be obtained as [34]
x 1 x 2 x n = C 1 e λ 1 t C 2 e λ 2 t C n e λ n t ,
where, by C 1 , C 2 , , C n , we have denoted the n components of an arbitrary vector quantity having constant components, while λ 1 , λ 2 , , λ n denote the proper values (eigenvalues) of the matrix A, which are determined as the algebraic roots of the characteristic polynomial p ( λ ) , which are defined according to the relation,
det A λ I n = 0 ,
where, by I n , we have denoted the identity matrix, defined in the standard way.
Definition 1
([34]). Let us assume that a solution of the system of the ordinary differential Equation (1) is known. The solution is called stable if and only if all the roots λ 1 , λ 2 , , λ n , of the characteristic polynomial p ( λ ) are located on the left-hand side of the complex plane, that is, for all roots λ, the condition Re λ < 0 is satisfied.
Assuming that the condition Re λ < 0 is satisfied, then it follows that, as t , x i ( t ) = e λ i t , i = 1 , 2 , , n tend exponentially to zero for all i. Hence, the point x 1 , x 2 , , x n = 0 , 0 , , 0 is stable with respect to small (linear) perturbations of the system of differential equations.
In the case of n-dimensional system of ordinary differential equations, the characteristic polynomial is given by
p ( λ ) = λ n + a 1 λ n 1 + + a n = 0 ,
where the coefficients a i are all real numbers, a i R n , i = 1 , 2 , , n . Moreover, without any loss of generality, it is possible to consistently consider a n 0 , since otherwise we would obtain λ = 0 , and thus we would have a characteristic polynomial of order n 1 , having the coefficient of the zeroth order non-vanishing.
The important and fundamental necessary and sufficient conditions for the polynomial p λ to have all algebraic solutions satisfying the condition Re λ < 0 may be generally formulated as [34]
a n > 0 , D 1 = a 1 > 0 , D 2 = a 1 a 3 1 a 2 > 0 ,
D 3 = a 1 a 3 a 5 1 a 2 a 4 0 a 1 a 3 > 0 , , D k = a 1 a 3 1 a 2 a 4 0 a 1 a 3 0 1 a 2 0 0 a k > 0 ,
for all k = 1 , 2 , , n .
An alternative and important method for the study of the linear, Lyapunov stability question is to investigate the algebraic relations involving the non-zero solutions of the characteristic polynomial p λ . Thus, we obtain
s : = i = 1 n λ i = a 1 ,
μ 1 : = i , j , i j n λ i λ j = a 2 ,
p : = λ 1 λ n = 1 n a n .
From the values of these coefficients, and by taking into account their algebraic properties, we can determine essential and significant information on the stability properties of the system of ordinary differential Equation (1). in the following, for the sake of completeness of our discussion, we also present the
Remark 1
(Descartes’ Rule of Signs) [34]). Let us consider the characteristic polynomial (5) of the system of differential Equation (1), with the coefficient a n satisfying the condition a n > 0 . Let us denote by m the number of changes in sign in the sequence of the coefficients a n , a n 1 , , a 0 , ignoring any coefficients in the sequence that are zero. Then, the polynomial p ( λ ) has at most m roots that are real and positive. In addition, m, m 2 , m 4 , …, real positive roots of the polynomial do exist [34].
By taking ω : = λ , Descartes’ Rule of Signs gives essential clues about the potential existence of real negative roots of the characteristic polynomial, a piece of information that is crucial for the study of the stability of the systems of strongly nonlinear differential equations.
If the proper values of the Jacobian A associated with a system of differential equations, evaluated at the equilibrium point x 0 : = x 0 1 , x 0 2 , , x 0 n , are known, with Equation (3), we obtain the behavior of solution near x 0 . For example, if we consider a two-dimensional autonomous differential system, we obtain the following classification of the equilibrium points: if the eigenvalues of A have negative real parts, then in the phase plane, all solutions are converging towards the steady state (equilibrium point) x 0 . The point x 0 is named a hyperbolic sink (stable point). Moreover, if the real parts of the proper values (eigenvalues) of A are greater then zero, then all integral curves diverge from the equilibrium point, and x 0 is named a hyperbolic source (unstable point). If one eigenvalue is positive, and the other one is negative, the fixed point x 0 is a saddle point (unstable). If the proper values (eigenvalues) of A are complex conjugate pairs, and Re λ 0 , then the equilibrium point is a spiral point (stable if Re λ < 0 , and unstable otherwise). If the proper values (eigenvalues) of A are purely imaginary values ( Re λ = 0 ), the fixed point is named a center.
If the eigenvalues of the linearized system (2) evaluated in x 0 have nonzero real parts, then the equilibrium point is said to be hyperbolic. Otherwise, it is called nonhyperbolic. The relation between the linear stability of the system (1) and its linearization (2) at equilibrium points is given by the following Theorem.
Theorem 1
((Hartman–Grobman) [20]). Let us consider a system of ordinary differential equations x ˙ = f ( x ) , x R n , with the vector field f being C 1 . Let us assume that x ¯ is a hyperbolic fixed point of the considered system of differential equations. Then, a neighborhood of the point x ¯ , on which the flow is topologically equivalent to the flow of the linearization of the system of differential equations at x ¯ , does always exist.
For the sake of completeness, we mention that the linear stability of the system (1) near a nonhyperbolic point could be investigated with the aid of the Lyapunov function introduced through the following Theorem and definition.
Theorem 2
((Lyapunov stability theorem) [20]). Let us consider that a vector field x ˙ = f ( x ) , x R n is given. Let x ¯ denote an equilibrium point of the vector field x ˙ . Moreover, let Σ : U R be a C 1 function, defined on some neighborhood U of x ¯ , and having the properties
(i) Σ x ¯ = 0 and Σ ( x ) > 0 , if x x ¯ ;
(ii) Σ ˙ ( x ) 0 in U x ¯ .
Then, if the conditions (i) and (ii) are satisfied, the point x ¯ is stable. Furthermore, if
(iii) Σ ˙ ( x ) < 0 in U x ¯ , the point x ¯ is asymptotically stable.
The function Σ ( x ) from the above Theorem is called the Lyapunov function. It has the property that, near the equilibrium point x 0 , the integral curves f ( x ) are tangent to the surface levels of Σ , or they cross the surface level oriented towards their interior. In other words, Σ ( x 0 ) · f ( x 0 ) 0 .

2.2. KCC Stability Theory

In the following, we introduce the basic ideas, and geometric concepts, of the KCC theory. Our presentation mostly follows the similar expositions of the theory in [8,25], respectively.

2.2.1. Geometrization of Arbitrary Dynamical Systems

We introduce first a set of dynamical variables x i , i = 1 , 2 , , n , assumed to be defined on a real, smooth n-dimensional manifold M . We denote in the following by T M the tangent bundle of M . Usually, M is considered as R n , M = R n , and therefore T M = T R n = R n [24].
Let us consider now a particular subset Ω of the ( 2 n + 1 ) dimensional Euclidian space defined as R n × R n × R 1 . On Ω , we assume the existence of a 2 n + 1 dimensional coordinate system denoted x i , y i , t , i = 1 , 2 , , n , where we have also introduced the notations x i = x 1 , x 2 , , x n , and y i = y 1 , y 2 , , y n , respectively. By t, we denote the ordinary time coordinate. We define the coordinates y i as
y i = d x 1 d t , d x 2 d t , , d x n d t .
A fundamental supposition in the KCC theory refers to the coordinate t, which can be interpreted physically as the time variable, and which is assumed to be an absolute invariant, which does not change in the coordinate transformations. Thus, by taking into account this assumption, on the base manifold M , the only allowed transformations of the coordinates are of the general form [24],
t ˜ = t , x ˜ i = x ˜ i x 1 , x 2 , , x n , i = 1 , 2 , , n .
In various mathematical applications of scientific importance, and interest, the equations of motion describing the evolution of natural or engineering systems are obtainable from a Lagrangian function L, which describes the state of the system, and is an application L : T M R . The dynamical evolution of the system can be obtained with the help of the Euler–Lagrange equations, given by [24]
d d t L y i L x i = F i , i = 1 , 2 , , n .
For the specific case of mechanical systems, the quantities F i , i = 1 , 2 , , n give the components of the external force F , which cannot be derived from a potential. By assuming that the Lagrangian L is regular, by using some simple calculations, we obtain the fundamental results that the Euler–Lagrange equations Equation (12) are equivalent mathematically to a complicated system of second-order ordinary, strongly nonlinear differential equations, given by [8,24]
d 2 x i d t 2 + 2 G i x j , y j , t = 0 , i = 1 , 2 , , n .
As for the functions G i , we assume that each function G i x j , y j , t is C in a neighborhood of some initial conditions x 0 , y 0 , t 0 , defined in Ω .
The key conjecture of the KCC theory is the following. Let us assume that an arbitrary system of strongly nonlinear second-order ordinary differential equations of the form (13) is defined in a general form. Even if the Lagrangian function for the system is not known a priori, we still can investigate the evolution and the behavior of its trajectories by using techniques suggested by the differential geometry of the Finsler spaces. This analysis can be carried out due to the existence of a close similarity between the paths of the Euler–Lagrange system, and the geodesics in a Finsler geometry.

2.2.2. The Nonlinear and Berwald Connections, and the KCC Invariants Associated with a Dynamical System

To investigate from a geometrical perspective, the mathematical properties of the dynamical system described by the system of differential equations Equation (13), we first introduce the nonlinear connection N, defined on the base manifold M , and having the coefficients N j i given by [35]
N j i = G i y j .
From a general geometric perspective, the nonlinear connection N j i can be characterized with the help of a dynamical covariant derivative N by using the following procedure. Let us assume that two vector fields v and w, respectively, are given, with both vector fields defined over a manifold M . Next, we define the covariant derivative N of w as [36]
v N w = v l x l w i + N j i ( x , y ) w j x i .
For N i j ( x , y ) = Γ i l j ( x ) v l , from Equation (15), we directly reobtain the definition of the standard covariant derivative for the particular case of the Levi–Civita linear connection, as introduced usually in the Riemannian geometry [36].
Now, we consider the open subset Ω R n × R n × R 1 on which the system of differential Equation (13) is defined, together with the coordinate transformations, defined by Equation (11), and assumed to be non-singular. On this mathematical structure, we can define the KCC-covariant derivative of an arbitrary vector field ξ i ( x ) by means of the definition [9,10,11,12],
D ξ i d t = d ξ i d t + N j i ξ j .
By taking ξ i = y i , we obtain
D y i d t = N j i y j 2 G i = ϵ i .
The contravariant vector field ϵ i is defined on the subset Ω of the Euclidian space, as one can see immediately from the above equation. It represents the first KCC invariant. From a general physical perspective, and within the mathematical formalism of the classical Newtonian mechanics, ϵ i , the first KCC invariant, could be understood as an external force, not derivable from a potential, and acting on the dynamical system.
As a next step in our discussion of the KCC theory, we consider the infinitesimal variations of the trajectories x i ( t ) of the dynamical system (13) into neighbouring ones, with the variations defined according to the prescriptions,
x ˜ i t = x i ( t ) + η ξ i ( t ) , y ˜ i t = y i ( t ) + η d ξ i ( t ) d t ,
where, by η < < 1 , we have denoted a small infinitesimal quantity, while ξ i ( t ) denotes the components of an arbitrary contravariant vector field ξ . The vector field ξ is defined along the trajectory x i ( t ) of the system of the differential equations under consideration. After the substitution of Equation (18) into Equation (13), and by considering the limit η 0 , we arrive to the deviation, or Jacobi, equations, representing the central mathematical result of the KCC theory, and which are given by [9,10,11,12]
d 2 ξ i d t 2 + 2 N j i d ξ j d t + 2 G i x j ξ j = 0 .
With the help of the KCC-covariant derivative, as introduced in Equation (16), we can reformulate Equation (19) in a fully covariant form as
D 2 ξ i d t 2 = P j i ξ j ,
where we have denoted
P j i = 2 G i x j 2 G l G j l i + y l N j i x l + N l i N j l + N j i t .
In Equation (21), we have defined the important tensor G j l i , given by [8,9,10,11,12,35]
G j l i N j i y l ,
and which in the KCC theory is named the Berwald connection.
The tensor P j i is called the second KCC-invariant. It is the fundamental quantity in the KCC theory, and in the Jacobi stability investigations. Alternatively, P j i is also called the deviation curvature tensor, by indicating its essential geometric nature. Furthermore, we will call in the following Equation (20) the Jacobi equation. This equation exists in both Riemann or Finsler geometries. If one assumes that the system of Equation (13) corresponds to the geodesic motion of a physical system, then Equation (20) gives the so-called Jacobi field equation, which can always be introduced in the considered geometry.
The trace P = P i i of the curvature deviation tensor, constructed from P j i , is a scalar invariant. It can be calculated from the relation
P = P i i = 2 G i x i 2 G l G i l i + y l N i i x l + N l i N i l + N i i t .
Other important invariants can also be constructed in the KCC theory. The most commonly used invariants are the third, fourth and fifth invariants associated with the given dynamical system, whose evolution and behavior is represented by the second order system of nonlinear Equation (13). These invariants are introduced according to the definitions [9]
P j k i 1 3 P j i y k P k i y j , P j k l i P j k i y l , D j k l i G j k i y l .
From a geometrical perspective, P j k i , the third KCC invariant, can be described as a torsion tensor. P j k l i , the fourth KCC invariant, represents the equivalent of the Riemann–Christoffel curvature tensor, while D j k l i , the fifth KCC invariant, is called the Douglas tensor [8,9]. Let us point out now that, in a Berwald geometry, these three tensors can always be defined. It is also important to mention that, in the KCC theory, the five invariants defined above are the fundamental mathematical quantities that describe the geometrical properties, and interpretation, of a dynamical system whose time evolution and behavior are represented by an arbitrary system of second-order strongly nonlinear differential equations.

2.2.3. Jacobi Stability of Dynamical Systems

In a large number of scientific investigations, the analysis of the stability of biological, chemical, biochemical, engineering, medical or physical systems, as well as the study of the trajectories of a system of differential equations, as given, for example, by Equation (13), in the vicinity of a given point x i t 0 , is of fundamental significance for the understanding of its properties. Moreover, the study of the stability can give important information on the temporal evolution of a general dynamical system.
For simplicity, in the following, we adopt as the origin of the time variable t the value t 0 = 0 . Moreover, we define . , . as representing the canonical inner product of R n . We also introduce the null vector O defined in R n , O R n . Next, we assume that the trajectories x i = x i ( t ) of the system of differential Equation (13) represent smooth curves in the Euclidean space R n , endowed with the canonical inner product . , . . Moreover, to completely characterize the deviation vector ξ , we assume, as a general property, that it satisfies the set of two essential initial conditions ξ 0 = O and ξ ˙ 0 = W O , respectively [8,9,10,11].
To describe the dispersing/focusing trend of the trajectories of a dynamical system around t 0 = 0 , we introduce the following intuitive and simple mathematical picture. Let us assume first that the deviation vector ξ satisfies the condition ξ t < t 2 , t 0 + [8,24]. If this is the case, then it turns out that all the trajectories of the dynamical system are focusing together, and converge towards the origin. Let us assume now that the deviation vector ξ has the property that the condition ξ t > t 2 , t 0 + is satisfied. In this case, it follows that all the solutions of the system of differential Equation (13) have at infinity a dispersing behavior [8,9,10,11].
The dispersing/focusing behavior of the solutions of a given system of second order ordinary differential equations can be also characterized, from the geometrical perspective introduced by the KCC theory, by considering the algebraic properties of the deviation curvature tensor P i j . This can be achieved in the following way. For t 0 + , the trajectories of the system of strongly nonlinear equations Equation (13) are bunching/focusing together if and only if the real parts of the characteristic values (eigenvalues) of the deviation tensor P j i 0 are strictly negative. Otherwise, the trajectories of the dynamical system have a dispersing behavior if and only if the real parts of the characteristic values (eigenvalues) of P j i 0 are strictly positive [8,9,10,11].
By taking into account the qualitative discussion presented above, we present now the rigorous mathematical definition of the notion of Jacobi stability for an arbitrary dynamical system, described by a system of ordinary second order differential equations, which is given by the following [8,9,10,11]:
Definition 2.
Let us consider that the general system of differential equations Equation (13), describing the time evolution of a dynamical system, satisfies the initial conditions:
| | x i t 0 x ˜ i t 0 | | = 0 , | | x ˙ i t 0 x ˜ ˙ i t 0 | | 0 ,
defined with respect to the norm . induced by a positive definite inner product.
If these conditions are satisfied, the trajectories of the dynamical system, given by Equation (13), are designated as Jacobi stable if and only if the real parts of the characteristic values (eigenvalues) of the curvature deviation tensor P j i are everywhere strictly negative.
On the other hand, if the real parts of the characteristic values (eigenvalues) of the curvature deviation tensor P j i are strictly positive everywhere, the trajectories of the dynamical system are designated as unstable in the Jacobi sense.
This definition allows us to straightforwardly investigate the stability of the systems of second order differential equations, and of the associated dynamical systems, as an alternative to the standard Lyapunov linear stability method.

2.3. The Correlation between Lyapunov and Jacobi Stability for a Two-Dimensional Autonomous Differential System

Let us recall that the Lyapunov stability is determined by the nature and sign of the eigenvalues of the Jacobian matrix evaluated at an equilibrium point (fixed point). On the other hand, the Jacobi stability is given by the sign of the real part of the eigenvalues of the curvature deviation tensor P j i , calculated at the same point.
The Jacobi stability of the two-dimensional systems of first order differential equations was explored in [10,11], where the authors considered a system of two arbitrary differential equations written in the general form,
d u d t = f ( u , v ) , d v d t = g ( u , v ) .
Moreover, it was assumed that the point ( 0 , 0 ) is a fixed point of the system (25), that is, f ( 0 , 0 ) = g ( 0 , 0 ) = 0 . If the equilibrium point is ( u 0 , v 0 ) ( 0 , 0 ) , with the change of the variable u ˜ = u u 0 , and v ˜ = v v 0 , respectively, the equilibrium point ( u 0 , v 0 ) is moved to the origin ( 0 , 0 ) .
In the approach pioneered in [10,11], after relabelling the variables by denoting v as x, and g ( u , v ) as y, and by also assuming that the condition g u | ( 0 , 0 ) 0 is satisfied by the function g ( u , v ) , it turns out that it is possible to eliminate from the system of the two equations (25) the variable u. By taking into account that the point ( u , v ) = ( 0 , 0 ) is a fixed point, from the Theorem of the Implicit Functions, it follows that, in the neighbourhood of the point ( x , y ) = ( 0 , 0 ) , the algebraic equation g ( u , x ) y = 0 has a unique solution u = u ( x , y ) . By taking into account that x ¨ = g ˙ = g u f + g v y , where the subscripts denote the partial derivatives with respect to u and v, respectively, one obtains finally an autonomous one-dimensional second order equation, equivalent mathematically to the system (25), and which is obtained in the general form as
x ¨ 1 + g 1 ( x , y ) = 0 ,
where
g 1 ( x , y ) = g u ( u ( x , y ) , x ) f ( u ( x , y ) , x ) g v ( u ( x , y ) , x ) y .
Hence, the Jacobi stability properties of a system of two arbitrary first order differential equations can be studied in detail via the equivalent Equation (26) by using the methods of the KCC theory [10,11]. Thereupon, the KCC stability properties of a first order system of differential equations of the form (25) can also be determined easily, thus allowing an in depth comparison between the Jacobi and Lyapunov stability properties of the two-dimensional dynamical systems, which can be performed in a straightforward manner [8].
Let us recall some results from [8] which are applied in this paper. The Jacobian matrix of (25) is
J ( u , v ) = f u f v g u g v .
The characteristic equation is given by
λ 2 ( tr A ) λ + det A = 0 ,
where tr A = f u + g v , and det A = f u g v g u f v are the trace and the determinant of the Jacobian matrix A = J | ( 0 , 0 ) , respectively.
The signs of the discriminant Δ = ( f u g v ) 2 + 4 f v g u , and the trace and the determinant of A, give the Lyapunov (linear) stability of the fixed point ( 0 , 0 ) .
The system (25) is equivalent from a mathematical point of view to the second order differential Equation (26). Performing the Jacobi stability analysis for this last equation in [8], assuming that g u | ( 0 , 0 ) 0 , the authors obtained the result that the curvature deviation tensor P 1 1 at the fixed point ( 0 , 0 ) is given by
4 P 1 1 | ( 0 , 0 ) = 4 g , 1 1 | ( 0 , 0 ) + ( g ; 1 1 ) 2 | ( 0 , 0 ) = Δ ,
where
Δ : = ( t r A ) 2 4 det A
is the discriminant of the characteristic Equation (29) and A = J | ( 0 , 0 ) . Therefore, the authors of [8] concluded their analysis by formulating the following
Theorem 3
([8]). Let us consider the system of two first order ordinary differential equations, given by Equation (25), with the fixed point P ( 0 , 0 ) , such that g u | ( 0 , 0 ) 0 . Then, the trajectory v = v ( t ) is Jacobi stable if and only if Δ < 0 .
Boehmer et al. [8] emphasized also that, if f v | ( 0 , 0 ) 0 , eliminating the variable v and relabeling u as x, we obtain a similar result that the trajectory u = u ( t ) is Jacobi stable if and only if Δ < 0 .
In other words, they demonstrated that: if one considers the system of ordinary differential Equation (25) with the fixed point P ( 0 , 0 ) located in the origin of the coordinate system, and satisfying the condition g u | ( 0 , 0 ) 0 , then the Jacobian matrix J evaluated at the point P has complex proper values (eigenvalues) if and only if P satisfies the condition of being a stable point in the Jacobi sense.
Let us recall that the condition of stability in the Lyapunov sense for a solution of the dynamical system (25) is that Re λ < 0 for all roots of the characteristic (eigenvalue) Equation (29). The condition for the Jacobi stability thus requires that the discriminant of the same characteristic (eigenvalue) equation to take negative values. Consequently, Lyapunov stability is in general not equivalent with Jacobi stability, and it is worth finding out the equilibrium points where a system is stable in both Jacobi and Lyapunov sense.
In the next section, we will investigate the relation between Jacobi and Lyapunov stability of the circular orbits around a GMGHS black hole.

3. Black Hole Solutions in Dilaton Gravity

Theoretical models based on the field equations obtained from the string-like action [28,29]
S = 1 2 d 4 x g R 2 γ ϕ 2 + e 2 α ϕ F 2 ,
where α = constant > 0 , and γ = constant > 0 , have been extensively investigated in the physical literature. In Equation (31), ϕ denotes the dilaton field, while F represents the Maxwell two-form field, having the components defined as F μ ν = μ A ν ν A μ . The action (31) can be obtained from the string frame low energy effective action, given by
S ^ = d 4 x g ^ e 2 α ϕ R ^ 2 γ ^ ϕ 2 + F 2 ,
by using the conformal transformation g ^ μ ν = e 2 α ϕ g μ ν . As particular cases, the models constructed from the action (32) include the Einstein–Maxwell theory, corresponding to α = γ = 0 , and low energy string theory, which is obtained by taking α = 2 , and γ = 2 , respectively.
The gravitational field equations derived from the action (32) by varying the metric are given by [29]
μ e 2 ϕ F μ ν = 0 , 2 ϕ + 1 2 e 2 ϕ F 2 = 0 ,
and
R μ ν = 2 μ ϕ ν ϕ + 2 e 2 ϕ F μ ρ F ν ρ 1 2 g μ ν e 2 ϕ F 2 ,
respectively. By assuming for the static, spherically symmetric metric an ansatz of the form,
d s 2 = λ 2 d t 2 + d r 2 λ 2 + R 2 d Ω ,
where λ and R are functions of r only, Gibbons and Maeda [28], and, independently, three years later, Garfinkle, Horowitz and Strominger [29] found that the unique static charged black hole solution corresponding to the action (31) is given by
d s 2 = 1 2 M r d t 2 + 1 1 2 M r d r 2 + r r Q 2 M ( d θ 2 + sin 2 θ d φ 2 )
where the electric field strength, and the dilaton field is given by
F r t = Q r 2 , e 2 α ϕ = 1 Q 2 M r .
This form of the metric is different from the usual standard form of a spherically symmetric metric, like, for example, the Schwarzschild metric. The above black hole solution (36) is known as a GMGHS black hole. For Q 2 < 2 M 2 , the black hole has an event horizon, and if Q 2 = 2 M 2 , the solution corresponds, from a physical point of view, to a naked singularity. This later case is called the extremal GMGHS black hole solution.
An external observer cannot see the region inside the event horizons. The main physical characteristics of a black hole can be acquired by investigating the behavior of matter and light outside the event horizons. A massive test particle moves along time-like geodesics, and photons move on null geodesics. The geodesics of a GMGHS black hole were studied extensively, for both the extremal and the non-extremal cases [37,38,39,40,41,42]. In what follows, we will consider the Jacobi stability of circular time-like geodesics around a GMGHS black hole. The existence and the Lyapunov stability of the circular orbits around GMGHS black holes was already investigated in [41].

Geodesic Equations in a Static-Charged Black Hole Geometry

We derive now the geodesics equations in the GMGHS space-time by using the mathematical formalism based on the Euler–Lagrange Equation (12). The Lagrangian for the metric (36) is:
2 L = 1 2 M r t ˙ 2 + r ˙ 2 1 2 M r + r r Q 2 M θ ˙ 2 + sin 2 θ φ ˙ 2
where by a dot we have denoted the differentiation with respect to τ —an affine parameter defined along the geodesic line. The parameter τ is chosen so that the Lagrangian L satisfies the condition 2 L = 1 on a time-like geodesics, 2 L = 0 on a null geodesics, and 2 L = 1 on a space-like geodesics, respectively. Moreover, all the functions F i , i = 1 , 2 , n , from the right-hand side of Equation (12) are zero.
We recall that the coordinates t and φ are cyclic. Thus, we obtain that
1 2 M r t ˙ = constant = E ,
is the energy integral, where E is the total energy of the particle, and
2 sin 2 θ · r 1 2 M r φ ˙ = constant = L ,
is the integral of the angular momentum.
The Euler–Lagrange equation for θ is
d d τ r r Q 2 M θ ˙ = r r Q 2 M sin θ cos θ · φ ˙ 2 .
We note that, if θ = π / 2 , then θ ˙ 0 , θ ¨ 0 , and therefore θ = π / 2 is located on the geodesic curve. Thus, the motion of a massive particle is planar. On the other hand, the angular momentum integral (40) takes the form,
r 1 2 M r φ ˙ = L ,
where L represents, from a physical point of view, the angular momentum of the particle, oriented in the direction of an axis perpendicular to the plane in which the motion of the particle takes place.
Being complicated, the Euler–Lagrange equation for r is replaced with the constancy of the Lagrangian,
1 2 d r d τ 2 + 1 2 1 2 M r L 2 r r Q 2 M ϵ = E 2 ,
where the constant ϵ takes the numerical values ϵ = 1 for time-like geodesics, ϵ = 0 for null geodesics, and ϵ = + 1 for space-like geodesics, respectively. The second term on the left-hand side of Equation (43) is the effective potential. Hence, Equation (43) can be written as
1 2 d r d τ 2 + V ( r ) = E 2
A massive test particle moving freely around a GMGHS black hole describes a time-like geodesic line. For them ϵ = 1 , and the effective potential, becomes
V ( r ) = 1 2 1 2 M r L 2 r r Q 2 M + 1 .
By introducing the new variable η , defined as r = 2 M η , the effective potential (45) becomes
V ( η ) = 1 2 1 1 η l 2 η 2 1 q 2 η + 1 ,
where we have denoted l 2 = L 2 / 4 M 2 , and q 2 = Q 2 / 2 M 2 , respectively. The variation of the potential V ( η ) is presented, for different values of l and q, in Figure 1.
The first derivative of the potential can be obtained immediately as
V ( η ) = 1 2 η 2 1 + l 2 η 2 ( 1 2 η ) q 2 2 ( 1 3 2 η ) η 1 q 2 η 2 .
The variation of the derivative V ( η ) of the potential as a function of η is represented in Figure 2.
For the second derivative of the potential, we obtain
V ( η ) = 2 η 3 1 l 2 η 4 3 ( 1 2 η ) η 2 + ( 1 3 η ) q 4 3 ( 1 8 3 η ) η q 2 1 q 2 η 3 .
The variation of V ( η ) is presented, for different values of l and q, in Figure 3.
The dependence of the real solution η 0 of the equation V ( η 0 ) = 0 on the parameters l 2 and q 2 of the GMGHS black hole solution is represented in Figure 4.

4. Stability of the Circular Orbits of the Free Test Particles in a GMGHS Spacetime

Next, we will consider the stability in Lyapunov and Jacobi sense of the circular orbits on which massive test particles move freely around a GMGHS black hole.

4.1. Lyapunov Stability Analysis

Equation (43) is the starting point in the dynamical systems approach for the analysis of the stability of the geodesic curves in the GMGHS geometry. Differentiating Equation (43) with respect to τ and dividing the result with r ˙ , we obtain the following second order differential equation
r ¨ = d V d r .
Equation (49) corresponds to following system of first order differential equations:
d r d τ = p d p d τ = d V d r .
The Jacobian matrix of the system (50) is
J = 0 1 V ( r ) 0
where means differentiation with respect to r.
The characteristic equation is
λ 2 + V ( r ) = 0 ,
and the proper values (eigenvalues) of the Jacobian matrix associated with the system (50) are given by
λ = ± V ( r ) ,
and so a simple fixed point ( r 0 , 0 ) of (50) is a saddle point if ( V ( r 0 ) < 0 ) and a center if ( V ( r 0 ) > 0 ) .
In [41], the existence and stability in the sense of Liapunov of circular time-like geodesics around a GMGHS black hole were explored. The study showed that, for certain values of the parameters, outside the black hole, there are two circular geodesics. In other words, outside the black hole, we can find two values r 0 > 2 M so that V ( r 0 ) = 0 , one for a minimum of the effective potential ( V ( r 0 ) < 0 ) and the other for a maximum of the potential ( V ( r 0 ) > 0 ).
If V ( r 0 ) < 0 , the eigenvalues (53) of the linearization of the system (50) are real and have opposite signs; therefore, the fixed point ( r 0 , 0 ) of the system is a saddle point, which is Liapunov unstable. A saddle point is a hyperbolic point, and so, based on the Hartman–Grobman theorem, the fixed point ( r 0 , 0 ) of the system (50) is Liapunov unstable.
If V ( r 0 ) > 0 , the values of λ from (53) are purely complex conjugate, and the study of the linear stability of system at the fixed point ( r 0 , 0 ) begins with the search of a Liapunov function for the system (50). We note that the function
V ( r , p ) = p 2 2 + V ( r )
has the property that V ( r , p ) · f ( r , p ) = 0 , ( r , p ) , where f ( r , p ) = p , d V / d r , meaning that the function (54) could be chosen as a Liapunov function. Furthermore, we have to check if V fulfills the condition of the Liapunov stability theorem. In other words, we have to verify if the fixed point is a local minimum of V . Therefore, we compute the Hessian matrix of V
H V = V ( r ) 0 0 1 .
We note that the matrix (55) is positive definite when V ( r ) > 0 , meaning that the fixed point ( r 0 , 0 ) of the linearized system corresponding to (50) is a center. In addition, thus, we finally conclude that the circular orbits around a GMGHS black hole, r = r 0 = constant , are Lyapunov stable when V ( r 0 ) > 0 , and Lyapunov unstable when V ( r 0 ) < 0 .

4.2. Jacobi Stability Analysis

In this section, we will first perform the study of the Jacobi stability analysis of Equation (49), giving the geodesic trajectories of a massive particle in the GMGHS geometry. Then, we will consider the stability of the circular orbits in both Lyapunov and KCC approaches.

4.2.1. Stability of the GMGHS Orbits in the General Case

The affine parameter τ in the geodesic Equation (49) is an absolute invariant, and hence all the results of the KCC theory can be applied to this case. By introducing the dimensionless radial coordinate η , the geodesic equation of motion in the GMGHS geometry takes the form
d 2 η τ d τ 2 + 1 4 M 2 d V η d η = 0 ,
or, equivalently,
d 2 η ( τ ) d τ 2 + 1 4 M 2 { 1 η 2 ( τ ) 1 + l 2 η 2 ( τ ) 1 q 2 η ( τ ) 1 1 η ( τ ) 2 l 2 η 3 ( τ ) 1 q 2 η ( τ ) + l 2 q 2 η 4 ( τ ) 1 q 2 η ( τ ) 2 } = 0 .
By denoting η τ = x 1 , and η τ = y 1 , Equation (56) takes the form
d 2 x 1 d τ 2 + 2 G 1 x 1 = 0 ,
where
G 1 ( x 1 ) = 1 8 M 2 d V x 1 d x 1 ,
or
G 1 x 1 = 1 8 M 2 { 1 x 1 2 1 + l 2 x 1 2 1 q 2 x 1 1 1 x 1 2 l 2 x 1 2 1 q 2 x 1 + l 2 q 2 x 1 2 1 q 2 x 1 2 } .
The nonlinear connection associated with Equation (58) is obtained as
N 1 1 = G 1 x 1 y 1 0 .
For the Berwald connection, we have
G 11 1 = N 1 1 y 1 0 .
For the curvature deviation tensor, we obtain now the simple expression
P 1 1 = 2 G 1 x 1 x 1 = 1 4 M 2 V ( x 1 ) ) .
Hence, the geodesic trajectories of a massive test particle in the spherically symmetric GMGHS black hole are Jacobi stable if the condition V x 1 x 1 = x 0 1 < 0 . On the other hand, we obtain a geometric interpretation of the second derivative of the potential, giving the deviation curvature tensor of the geodesic trajectories. On the other hand, for the first KCC invariant of the system, we obtain
ϵ 1 = 2 G 1 ( x 1 ) = 1 4 M 2 d V x 1 . d x 1
Hence, the first derivative of the potential, representing physically the force acting on the particle, has a geometric interpretation as the first KCC invariant. Moreover, the geodesic deviation equations take the form
d 2 ξ 1 d τ 2 + 1 4 M 2 V ( x 1 ) ξ 1 = 0 .
The variation of the deviation curvature tensor of the GMGHS black hole is represented as a function of the solution parameters l and q, for a fixed value of η , in Figure 5. The contour plot of the deviation tensor is also represented.

4.2.2. Stability of the Circular Orbits

In this section, we will perform the Jacobi stability of the system (50) using the results obtained by Boehmer et al. in [8]. For the system (50), we consider the Jacobian matrix of the system and evaluate its trace and determinant
tr J = 0 det J = V ( r ) .
Thus, the discriminant of the characteristic equation becomes
Δ = tr ( J ) 2 4 det ( J ) = 4 V ( r ) .
Based on the Theorem which makes the link between the discriminant of the characteristic equation and the Jacobi stability, we can conclude that the circular orbit of a free particle moving in an GMGHS spacetime, r = r 0 , is Jacobi stable when V ( r 0 ) > 0 and Jacobi unstable for V ( r 0 ) < 0 .
We note that we have found the same condition for Jacobi stability as for the for Lyapunov stability, meaning that, for the circular orbits on which massive test particle move around a GMGHS black hole, the two types of stability coincide.
We consider now some particular cases of stability, corresponding to some specific values of the parameters l and q of the GMGHS black hole. By taking l = 8 , and q = 0.95 , it turns out that the equation V ( η ) = 0 has two real solutions satisfying the condition η > 1 , given by η 1 = 1.1996 and η 2 = 127.8522 . Evaluating the second derivative of the potential in these points gives V ( η ) η = 1.19996 = 185.5041 < 0 , and V ( η ) η = 127.8522 = + 2.39 × 10 7 > 0 , respectively. Hence, we can conclude that the circular trajectory located at r = 2 M η 1 is both Lyapunov and Jacobi unstable (left panel in Figure 6), while the circular trajectory located at r = 2 M η 2 is both Lyapunov and Jacobi stable (right panel in Figure 6).
Let us now consider the Lyapunov stability corresponding to the value l = 3 / 2 of the parameter l of a GMGHS black hole. By taking q { 0.55 , 0.6512 , 0.75 } , the equation V ( η ) = 0 has no real solution, two real and equal solutions, and two different real solutions, satisfying the condition η > 1 , meaning that, outside the black hole, there are no circular orbits, one circular orbit and two circular orbits, respectively. For q = 0.6512 , η 1 = η 2 = 2.5243 , by evaluating the second derivative of the potential, we obtain V ( η ) η = 2.5243 = 0 , the point η = 2.5243 is an inflection point for the potential, which leads to a cusp in the phase diagram (see the middle panel in Figure 7). For q = 0.75 , η 1 = 1.8333 and η 2 = 3.3836 and by evaluating the second derivative of the potential in these points, we obtain V ( η ) η = 1.8333 = 0.111 < 0 and V ( η ) η = 3.3836 = + 0.0075 > 0 , respectively. Therefore, we can conclude that the circular trajectory located at r = 2 M η 1 is Lyapunov unstable, and the circular trajectory located at r = 2 M η 2 is Lyapunov stable (see the right panel in Figure 7).

5. Conclusions

In the present paper, we have first revisited, and carefully investigated, two methods of stability analysis: Lyapunov stability and the Jacobi stability approaches, respectively. The Lyapunov stability analysis is conducted by the linearization of the system of differential equations describing the dynamical system at the fixed points. On the other hand, and Jacobi stability involves the perturbations of a whole set of trajectories. Intuitively, the Jacobi stability indicates how the trajectories bunch together, or disperse, when approaching the fixed point.
As an application of the two stability methods, we have comparatively investigated the behavior of the trajectories of the solutions for a specific charged black hole solution that originates in the low energy limit of string theory, called the GMGHS solution. The study of the properties of the geodesic curves in black hole geometries is an important field of investigation [37,38,39,40,41,42], which could lead not only to a better understanding of the theoretical properties of these objects, but could also open new perspectives on their observational detection. Moreover, black holes represent a fertile testing ground of modified gravity theories.
An analysis of the stability of the orbits in the Schwarzschild geometry was performed in [20], by using both Lyapunov and Jacobi stability approaches. As a result of this study, it was shown that stable circular orbits do exist at a radius r = R + , where R + = L 2 + L 4 12 L 2 M 2 / 2 M , while unstable circular orbits exist at r = R , where R = L 2 L 4 12 L 2 M 2 / 2 M . A similar analysis of the motion of the particles in Newtonian mechanics in the presence of a central force field f ( r ) was carried out in [19]. In nonrelativistic mechanics, circular orbits in a central field do exist if the equation
V ( r ) = f ( r ) L 2 M r 3 = 0 ,
has real roots. Furthermore, if r = r 0 is a root of V ( r ) = 0 , the circular orbit is stable if the condition
V ( r ) r = r 0 = f ( r ) r = r 0 + 3 L 2 M r 0 4 > 0 ,
is satisfied.
In a realistic astrophysical environment, massive general relativistic objects, like, for example, black holes or neutron stars, are often enclosed by an accretion disk. Accretion disks around compact objects can be the basis of physical models that could convincingly provide explanations for many astrophysical phenomena, like, for example, active galactic nuclei and X-ray binaries. The disks can be described theoretically by assuming that they are composed of massive test particles (baryons) that evolve in the gravitational field of the central massive and compact astrophysical object. The disks cool down through the electromagnetic radiation emission from their surface, and this form of energy emission represents an efficient physical mechanism for avoiding the extreme heating of the disk surface [43]. The disk has an inner edge, which is located at the marginally stable orbits of the gravitational, potential created by the central massive compact object. Hence, in higher orbits, the motion of the gas in the disk is Keplerian.
Therefore, the problem of the determination of the position of circular orbits, and of their stability, is fundamental from an astrophysical point of view. The electromagnetic emissivity properties of the accretion disks provide distinct observational signatures for different classes of astrophysical objects, including black holes, neutron, quark or other types of exotic stars. The parameters l and q of the GMGHS black hole can also be constrained from the physical properties of the accretion disks. The condition of the stability of the particle trajectories in the disk also imposes strong constraints on the parameters of central object. In this respect, the results obtained via the applications of the concept of Jacobi and Lyapunov stability may prove to be essential for the understanding of the nature of the black holes, or other types of compact objects.
For example, in [44], it was shown that the equation governing the vertical perturbations of the trajectories of the test particles in the equatorial orbits around massive general relativistic objects is given by
d 2 δ z d s 2 + ν d δ z d s + ω 2 δ z = ξ z g ; z ,
where δ z is the perturbation of the coordinate z, ν is a constant, ξ z g ; z is the external force, and
ω 2 = Γ t t , z z + 2 Γ t ϕ , z z Ω c + Γ ϕ ϕ , z z Ω c 2 u t 2 ,
respectively. In the above equation, Γ μ ν λ denote the Christoffel symbols of a Riemannian metric, ω is the azimuthal angular velocity, while u t is the temporal component of the four-velocity of the particles in the disk. To obtain Equation (70), it was assumed that the particles in the disk move along the geodesic lines. The presence of a viscous dissipation and of an external (stochastic force) was also taken into account. Since the vertical perturbations of the disk are described by a second order differential equations, the study of the stability of the disk around the GMGHS black holes can be analyzed by using the theoretical concepts discussed in the present work.
To conclude, in the present study, we have carried out an independent investigation of the stability of the geodesic trajectories of massive, baryonic test particles moving in GMGHS geometry, by using both the Lyapunov and the Jacobi methods for the stability analysis. We have obtained the basic mathematical result that the condition for Jacobi stability of circular orbits in a GMGHS spacetime is the same as the condition for Lyapunov stability, meaning that, in this case, these two types of stability are equivalent. This result is also a consequence of the two-dimensional nature of the system of differential equations, corresponding to the geodesic motion in the GMGHS geometry in spherical static symmetry. For higher dimensional dynamical systems, and in the presence of a complex behavior, the predictions of the Jacobi and Lyapunov stability theories may be different, thus allowing for a better explanation of the physical and mathematical properties of these systems on both qualitative and quantitative levels.

Author Contributions

Conceptualization, C.B.; methodology, P.B.; formal analysis, T.H.; writing—review and editing, T.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the three anonymous reviewers for comments and suggestions that helped us to improve our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boffetta, G.; Cencini, M.; Falcioni, M.; Vulpiani, A. Predictability: A way to characterize complexity. Phys. Rep. 2002, 356, 367–474. [Google Scholar] [CrossRef] [Green Version]
  2. Mancho, A.M.; Small, D.; Wiggins, S. A tutorial on dynamical systems concepts applied to Lagrangian transport in oceanic flows defined as finite time data sets: Theoretical and computational issues. Phys. Rep. 2006, 437, 55–124. [Google Scholar] [CrossRef]
  3. Motter, A.E.; Gruiz, M.; Károlyi, G.; Tél, T. Doubly Transient Chaos: Generic Form of Chaos in Autonomous Dissipative Systems. Phys. Rev. Lett. 2013, 111, 194101. [Google Scholar] [CrossRef] [Green Version]
  4. Donetti, L.; Hurtado, P.I.; Munoz, M.A. Entangled Networks, Synchronization, and Optimal Network Topology. Phys. Rev. Lett. 2005, 95, 188701. [Google Scholar] [CrossRef] [Green Version]
  5. Kosambi, D.D. Parallelism and path-spaces. Math. Z. 1933, 608, 608–618. [Google Scholar]
  6. Cartan, E. Observations sur le méemoir précédent. Math. Z. 1933, 37, 619–622. [Google Scholar] [CrossRef]
  7. Chern, S.S. Sur la géométrie d’un systéme d’equations differentialles du second ordre. Bull. Des Sci. Math. 1939, 63, 206–212. [Google Scholar]
  8. Boehmer, C.G.; Harko, T.; Sabau, S.V. Jacobi stability analysis of dynamical systems—Applications in gravitation and cosmology. Adv. Theor. Math. Phys. 2012, 16, 1145–1196. [Google Scholar] [CrossRef]
  9. Antonelli, P.L. (Ed.) Handbook of Finsler Geometry; Kluwer Academic: Dordrecht, HI, USA, 2003; Volume 1. [Google Scholar]
  10. Sabau, S.V. Systems biology and deviation curvature tensor. Nonlinear Anal. Real World Appl. 2005, 6, 563–587. [Google Scholar]
  11. Sabau, S.V. Some remarks on Jacobi stability. Nonlinear Anal. 2005, 63, 143–153. [Google Scholar] [CrossRef]
  12. Antonelli, P.L.; Bucataru, I. New results about the geometric invariants in KCC theory. An. St. Univ. ”Al.I. Cuza” Iasi. Mat. (N.S.) 2001, 47, 405–420. [Google Scholar]
  13. Yajima, T.; Nagahama, H. KCC-theory and geometry of the Rikitake system. J. Phys. A Math. Theor. 2007, 40, 2755–2772. [Google Scholar]
  14. Harko, T.; Sabau, V.S. Jacobi stability of the vacuum in the static spherically symmetric brane world models. Phys. Rev. D 2008, 77, 104009. [Google Scholar] [CrossRef] [Green Version]
  15. Boehmer, C.G.; Harko, T. Nonlinear Stability Analysis of the Emden-Fowler Equation. J. Nonlinear Math. 2010, 17, 503–516. [Google Scholar]
  16. Yajima, T.; Nagahama, H. Nonlinear dynamical systems and KCC-theory. Acta Math. Acad. Nyíregyházi. 2008, 24, 179–189. [Google Scholar]
  17. Yajima, T.; Nagahama, H. Tangent bundle viewpoint of the Lorenz system and its chaotic behavior. Phys. Lett. A 2010, 374, 1315–1319. [Google Scholar] [CrossRef]
  18. Abolghasem, H. Liapunov stability versus Jacobi stability. J. Dyn. Syst. Geom. 2012, 10, 13–32. [Google Scholar]
  19. Abolghasem, H. Jacobi Stability of Circular Orbits in a Central Force. J. Dyn. Syst. Geom. 2012, 10, 197–214. [Google Scholar] [CrossRef]
  20. Abolghasem, H. Stability of circular orbits in Schwarzschild spacetime. Int. J. Differ. Equ. Appl. 2013, 12, 131–147. [Google Scholar]
  21. Abolghasem, H. Jacobi stability of Hamiltonian systems. Int. J. Pure Appl. Math. 2013, 87, 181–194. [Google Scholar]
  22. Harko, T.; Ho, C.Y.; Leung, C.S.; Yip, S. Jacobi stability analysis of the Lorenz system. Int. J. Geom. Methods Mod. Phys. 2015, 12, 1550081. [Google Scholar] [CrossRef] [Green Version]
  23. Harko, T.; Pantaragphong, P.; Sabau, S. A new perspective on the Kosambi–Cartan–Chern theory, and its applications. arXiv 2015, arXiv:1509.00168. [Google Scholar]
  24. Harko, T.; Pantaragphong, P.; Sabau, S. Kosambi–Cartan–Chern (KCC) theory for higher-order dynamical systems. Int. J. Geom. Methods Mod. Phys. 2016, 13, 1650014. [Google Scholar]
  25. Danila, B.; Harko, T.; Mak, M.K.; Pantaragphong, P.; Sabau, S. Jacobi stability analysis of scalar field models with minimal coupling to gravity in a cosmological background. Adv. High Energy Phys. 2016, 2016, 7521464. [Google Scholar]
  26. Lake, M.J.; Harko, T. Dynamical behavior and Jacobi stability analysis of wound strings. Eur. Phys. J. C 2016, 76, 311. [Google Scholar]
  27. Blaga, C.; Blaga, P.A.; Harko, T. Jacobi stability analysis of the classical restricted three body problem. Rom. Astron. J. 2021, 31, 101–112. [Google Scholar]
  28. Gibbons, G.W.; Maeda, K. Black holes and membranes in higher-dimensional theories with dilaton fields. Nucl. Phys. B 1988, 298, 741–775. [Google Scholar]
  29. Garfinkle, T.; Horowitz, G.A.; Strominger, A. Charged black holes in string theory. Phys. Rev. D 1991, 43, 3140–3143. [Google Scholar]
  30. Wainwright, J.; Ellis, G.F.R. Dynamical Systems in Cosmology; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  31. Boehmer, C.G.; Chan, N.; Lazkoz, R. Dynamics of dark energy models and centre manifolds. Phys. Lett. B 2012, 714, 11–17. [Google Scholar]
  32. Boehmer, C.G.; Chan, N. Dynamical systems in cosmology. arXiv 2014, arXiv:1409.5585. [Google Scholar]
  33. Garcia-Salcedo, R.; Gonzalez, T.; Horta-Rangel, F.A.; Quiros, I.; Sanchez-Guzmán, D. Introduction to the application of dynamical systems theory in the study of the dynamics of cosmological models of dark energy. Eur. J. Phys. 2015, 36, 025008. [Google Scholar]
  34. Murray, J.D. Mathematical Biology; Springer: New York, NY, USA, 1993. [Google Scholar]
  35. Miron, R.; Hrimiuc, D.; Shimada, H.; Sabau, V.S. The Geometry of Hamilton and Lagrange Spaces; Kluwer Academic Publishers: Dordrecht, The Netherlands; Boston, MA, USA, 2001. [Google Scholar]
  36. Punzi, R.; Wohlfarth, M.N.R. Geometry and stability of dynamical systems. Phys. Rev. E 2009, 79, 046606. [Google Scholar] [CrossRef] [Green Version]
  37. Blaga, C.; Blaga, P.A. On the geodesics for a spherically symmetric dilaton black hole. Serb. Astron. J. 1998, 158, 55–59. [Google Scholar]
  38. Fernando, S. Null geodesics of charged black holes in string theory. Phys. Rev. D 2012, 85, 024033. [Google Scholar] [CrossRef] [Green Version]
  39. Pradhan, P.P. Horizon straddling ISCOs in spherically symmetric string black holes. Int. J. Mod. Phys. D 2015, 24, 1550086. [Google Scholar] [CrossRef] [Green Version]
  40. Olivares, M.; Villanueva, J.R. Massive neutral particles on heterotic string theory. Eur. Phys. J. C 2013, 73, 2659. [Google Scholar] [CrossRef] [Green Version]
  41. Blaga, C. Circular time-like geodesics around a charged spherically symmetric dilaton black hole. Automat. Comp. Appl. Math. 2013, 22, 41–49. [Google Scholar]
  42. Blaga, C. Timelike geodesics around a charged spherically symmetric dilaton black hole. Serb. Astron. J. 2015, 190, 41–48. [Google Scholar]
  43. Shahidi, S.; Harko, T.; Kovács, Z. Distinguishing Brans-Dicke-Kerr type naked singularities and black holes with their thin disk electromagnetic radiation properties. Eur. Phys. J. C 2020, 80, 162. [Google Scholar]
  44. Harko, T.; Mocanu, G.R. Stochastic oscillations of general relativistic discs. Mon. Not. R. Astron. Soc. 2012, 421, 3102–3110. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Variation as a function of η of the effective potential V ( η ) of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Figure 1. Variation as a function of η of the effective potential V ( η ) of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Symmetry 15 00329 g001
Figure 2. Variation as a function of η of the derivative V ( η ) of the effective potential of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Figure 2. Variation as a function of η of the derivative V ( η ) of the effective potential of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Symmetry 15 00329 g002
Figure 3. Variation as a function of η of the second derivative V ( η ) of the effective potential of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Figure 3. Variation as a function of η of the second derivative V ( η ) of the effective potential of the GMGHS black hole for q = 0.95 , and l = 16 (solid curve), l = 14 (dotted curve), l = 12 (short dashed curve), l = 10 (dashed curve), and l = 8 (long dashed curve), respectively (left panel), and for l = 10 , and q = 0.95 (solid curve), q = 0.85 (dotted curve), q = 0.75 (short dashed curve), q = 0.65 (dashed curve), and q = 0.55 (long dashed curve, respectively (right panel).
Symmetry 15 00329 g003
Figure 4. Variation as a function of l 2 and q 2 of the solution η 0 of the algebraic equation V ( η 0 ) = 0 .
Figure 4. Variation as a function of l 2 and q 2 of the solution η 0 of the algebraic equation V ( η 0 ) = 0 .
Symmetry 15 00329 g004
Figure 5. Variation of the curvature deviation tensor P 1 1 η , l , q as a function of l and q for η = 1.75 (left panel), and the contour plot of the curvature deviation tensor for η = 1.25 (right panel).
Figure 5. Variation of the curvature deviation tensor P 1 1 η , l , q as a function of l and q for η = 1.75 (left panel), and the contour plot of the curvature deviation tensor for η = 1.25 (right panel).
Symmetry 15 00329 g005
Figure 6. The phase plane for the specific values of the parameters l = 8 and q = 0.95 . The behavior of the solution near point η 1 = 1.1996 , represented in the left panel, shows that it is a saddle point. The behavior of the solution near the second point η 2 = 127.8522 , from the right panel, shows that the point is a center.
Figure 6. The phase plane for the specific values of the parameters l = 8 and q = 0.95 . The behavior of the solution near point η 1 = 1.1996 , represented in the left panel, shows that it is a saddle point. The behavior of the solution near the second point η 2 = 127.8522 , from the right panel, shows that the point is a center.
Symmetry 15 00329 g006
Figure 7. The phase plane for the specific values of the parameters l = 3 / 2 and q { 0.55 , 0.6512 , 0.75 } . The left panel corresponds to q = 0.55 , in which case the first derivative of the potential has no solution for η > 1 , thus there are no critical points in the phase diagram. The middle panel is for q = 0.6512 , in which case the first derivative of the potential has a double solution for η > 1 , meaning that, in the phase plane, there is a cusp in the point ( 2.5243 , 0 ) . The right panel corresponds to q = 0.75 , in which case there are two circular orbits around the black hole, corresponding to η 1 = 1.8333 and η 2 = 3.3836 , respectively, with the first generating a saddle point, while the second a center.
Figure 7. The phase plane for the specific values of the parameters l = 3 / 2 and q { 0.55 , 0.6512 , 0.75 } . The left panel corresponds to q = 0.55 , in which case the first derivative of the potential has no solution for η > 1 , thus there are no critical points in the phase diagram. The middle panel is for q = 0.6512 , in which case the first derivative of the potential has a double solution for η > 1 , meaning that, in the phase plane, there is a cusp in the point ( 2.5243 , 0 ) . The right panel corresponds to q = 0.75 , in which case there are two circular orbits around the black hole, corresponding to η 1 = 1.8333 and η 2 = 3.3836 , respectively, with the first generating a saddle point, while the second a center.
Symmetry 15 00329 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Blaga, C.; Blaga, P.; Harko, T. Jacobi and Lyapunov Stability Analysis of Circular Geodesics around a Spherically Symmetric Dilaton Black Hole. Symmetry 2023, 15, 329. https://doi.org/10.3390/sym15020329

AMA Style

Blaga C, Blaga P, Harko T. Jacobi and Lyapunov Stability Analysis of Circular Geodesics around a Spherically Symmetric Dilaton Black Hole. Symmetry. 2023; 15(2):329. https://doi.org/10.3390/sym15020329

Chicago/Turabian Style

Blaga, Cristina, Paul Blaga, and Tiberiu Harko. 2023. "Jacobi and Lyapunov Stability Analysis of Circular Geodesics around a Spherically Symmetric Dilaton Black Hole" Symmetry 15, no. 2: 329. https://doi.org/10.3390/sym15020329

APA Style

Blaga, C., Blaga, P., & Harko, T. (2023). Jacobi and Lyapunov Stability Analysis of Circular Geodesics around a Spherically Symmetric Dilaton Black Hole. Symmetry, 15(2), 329. https://doi.org/10.3390/sym15020329

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop