Next Article in Journal
Supersymmetric HS Yang-Mills-like Models
Next Article in Special Issue
Taxonomy of Dark Energy Models
Previous Article in Journal
An Improved Framework for Quantum Gravity
Previous Article in Special Issue
Model Independent Diagnostics in Interacting Dark Energy Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Stability and Geometrical Diagnostic of the Power Law K-Essence Dark Energy Model with Interaction

1
Department of Physics, Liaoning Normal University, Dalian 116029, China
2
School of Liberal Arts and Sciences, North China Institute of Aerospace Engineering, Langfang 065000, China
3
Peng Huanwu Center for Fundamental Theory, Hefei 230026, China
*
Author to whom correspondence should be addressed.
Universe 2020, 6(12), 244; https://doi.org/10.3390/universe6120244
Submission received: 21 October 2020 / Revised: 11 December 2020 / Accepted: 13 December 2020 / Published: 18 December 2020
(This article belongs to the Special Issue Probing the Dark Universe with Theory and Observations)

Abstract

:
We investigate the cosmological evolution of the power law k-essence dark energy (DE) model with interaction in FRWL spacetime with the Lagrangian that contains a kinetic function F ( X ) = X + X . Concretely, the cosmological evolution in this model are discussed by the autonomous dynamical system and its critical points, together with the corresponding cosmological quantities, such as Ω ϕ , w ϕ , c s 2 , and q, are calculated at each critical point. The evolutionary trajectories are drawn in order to show the dynamical process on the phases plan around the critical points. The result that we obtained indicates that there are four dynamical attractors, and all of them correspond to an accelerating expansion of universe for certain potential parameter and coupling parameter. Besides that, the geometrical diagnostic by the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) of this scalar field model are numerically obtained by the phase components, as an extended null diagnostic for the cosmological constant. This diagnostic shows that both the potential parameter λ and interaction parameter α play important roles in the evolution of the statefinder hierarchy.

1. Introduction

There are two important stages in cosmology, the early inflation and the late time accelerated expansion. The inflation that we postulate is to explain some issues, such as the flat problem and the horizon problem, etc. While, the accelerating expansion is based on the observations on the luminosity-redshift relation of distant Ia supernovas [1,2], Cosmic Microwave Background [3], and Baryon Acoustic Oscillations [4], which indicate that the current energy density in the universe is composed by 68.3% dark energy (DE), 26.8% dark matter (DM), and 4.9% baryons [5] in order to drive the late time acceleration approximately. More details and cosmic constrain by observation are in [6,7,8,9]. Since its first observation in 1998, over the last twenty years, there have been many models to make explanation for the physical mechanism of this phenomenon. Among them, the simplest one is the Λ CDM model with a constant equation of state (EoS) w Λ = 1 , which provides the negative pressure for the expansion. Λ CDM model is in good accordance with the observation, but it has some crucial problems, such as the cosmological constant problem, the age problem [10,11,12,13,14], and the tensions on the parameters H 0 and σ 8 in the Λ CDM model in recent years [15,16,17,18,19,20]. Instead of the Λ CDM model, there is a class of phenomenological models with a scalar field to reconcile the problems above; for example, quintessence, phantom, quintom, tachyon, k-essence, and DBI models, etc. Among them, the quintessence model is one of the most popular one [21]. Its Lagrangian density L = X V ( ϕ ) , being the pressure, in the action is described by a single scalar field ϕ , and a canonical kinetic energy X 1 2 μ ϕ μ ϕ . To some extent, the quintessence model can describe a range of EoS necessary for an accelerated expansion, but the single-field quintessence with a canonical kinetic term limits the EoS w to be less than 1 . To this point, the quintom model with multiple fields [22] and some k-essence models with non-canonical kinetic terms are developed in order to realize the crossing of the phantom divide. Additionally, the combination of data sets Plank+R16+JLA support this point, such as in [15], for w = 1.079 0.057 + 0.072 .
It is known that the k-essence model is a better substitution for the quintessence model. Although it firstly originated from k-inflation coming from the string theory to explain the early inflation of the universe [23]; fortunately, it could also explain the late time accelerated expansion as DE [24], even in the form as a decent generalization of many scalar field dark energy models. Unlike the canonical kinetics that are defined by the quintessence model, the k-essence model provides a variety of non-canonical terms, i.e., L = F ( X , ϕ ) [25], which contains higher order terms of X. As a result, the k-essence model with non-canonical kinetic term in Largragian can also reconcile the H 0 tension problem by making the EoS smaller than 1 . The form of F ( X , ϕ ) was first discussed as a purely kinetic function F ( X ) , with a constant potential, V ( ϕ ) = c o n s t [26]. Further more, two common forms along with potential are F ( X , ϕ ) = F ( X ) V ( ϕ ) and F ( X , ϕ ) = F ( X ) V ( ϕ ) , which were widely investigated in [27,28]. Some of the modified kinetic terms were discussed, such as F ( X ) = K X + L X 2 [29], F ( X ) = 1 2 α 1 ( ( A X ) α 2 α α 0 A X ) [30], and F ( X ) = A X B X α [31] etc., which belongs to the class of power law k-essence dark energy model with the power law function F ( X ) = a i X b i . The approximation of the potential in scalar field dark energy models are discussed in [32], with both canonical and noncanonical kinetic terms.
In another aspect, from the matter clustering properties, dark matter (DM) and dark energy are not the same substance; however, there are researches regarding the interactions between them, even some nonlinear interaction forms [33,34,35,36,37], which can provide a mechanism for generating acceleration. By the recent observation, the interaction between DM and DE is too little to alleviate the coincidence problem, while, in our work, the k-essence model with the interaction between DM and DE can be a candidate, which helps to explain the H 0 tension and σ 8 tension between CMB and structure formation measurements [15,16,17].
The aim of this paper, which is based on the researches above, is to consider a model with F ( X , ϕ ) = F ( X ) V ( ϕ ) f ( ϕ ) , where F ( X ) = X + X , V ( ϕ ) 1 / ϕ 2 , and f ( ϕ ) = 0 [38], together with a certain kind of interaction Q. We investigate the possible cosmological behavior of this model in Friedmann–Robertson–Walker–Lemaître (FRWL) spacetime by performing a phase-space and stability analysis. The theory are based on [39,40] judging the stability of the critical points by the eigenvalues; whereas, in this model for the convenience of calculation, it prefers the method by the determinant and trace of the Jacobian matrix of the autonomous differential equations [41]. Some cosmological quantities will be calculated for each critical point, such as the dark energy density parameter Ω ϕ , the equation of state (EoS) parameter w ϕ of dark energy, the sound speed c s 2 , and the deceleration parameter q.
Finally, in order to distinguish this k-essence model from Λ CDM model, there are two main kinds of “null measure”: the Om diagnostic and the statefinder diagnostic. Om is constructed from the Hubble parameter H, and it provides a null test of the Λ CDM model [42,43]. In recent years, in order to distinguish those models from the best fitting model, the Λ CDM model, the statefinder hierarchy is used, which originate from statefinder diagnostic [44,45]. The statefinder pair { r , s } , is composed by the scale factor a ( t ) with its second and third derivatives; however, the statefinder hierarchy is based on the even higher derivatives [46]. In this paper, statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) are analyzed to the scalar dark energy model by the phase components (the auxiliary variables) in the autonomous equations, unlike the method that was mentioned before, which depends on the cosmological quantities. It shows that the hierarchies are varied from two parameters λ and α by the trajectories. This novel method is based on [47,48,49], which is used for the quintessence model and I ω C D M model, and it could be generated to a range of scalar DE models in the future.
This paper is organized, as follows: in the following section, we review k-essence dark energy models and its stability analysis. In the third section, we consider the dynamics of the k-essence scalar field with the interaction Q = α H ρ m (the coupling parameter α is a real arbitrary constant). In the fourth section, the statefinder diagnostic and statefinder hierarchy are analyzed in order to distinguish from Λ CDM model. Finally, we close with a few concluding remarks in the fifth section.

2. The Power Law K-Essence Dark Energy Model and Its Stability Analysis

We consider k-essence dark energy models with Lagrangian
L = p ϕ = G ( X ) U ( ϕ ) ,
where the kinetic term G ( X ) and potential term U ( ϕ ) are analytic functions of X and ϕ , respectively. Throughout this paper, we will work with a flat, homogeneous, and isotropic FRWL spacetime having signature ( , + , + , + ) and in units c = 8 π G = 1 . We are interested in the power law k-essence with a general form of kinetic term G ( X ) = K ( ϕ ) X + L ( ϕ ) X , which has been studied in [29,38,50]. Hence, the scalar field ϕ is redefined as the one in [29]. Consequently, Equation (1) is rewritten as p ϕ = F ( X ) V ( ϕ ), where the new kenetic term F ( X ) = X + X and new potential term V ( ϕ ) = ( K 2 / L ) U ( ϕ ) . Subsequently, the corresponding energy density ρ ϕ , the EoS parameter w ϕ and the effective sound speed c s 2 are, respectively, given by
ρ ϕ = V ( ϕ ) [ 2 X F X F ] = X V ,
w ϕ = F 2 X F X F = X X X ,
c s 2 = p / X ρ / X = F X F X + 2 X F X X = 1 1 2 X ,
where F X d F / d X and F X X d 2 F / d X 2 . The sound speed comes from the equation describing the evolution of linear adiabatic perturbations in a k-essence dominated universe [29,51] (a non-adiabatic perturbation of k-essence has been discussed in [52,53], here we only consider the case of adiabatic perturbation). From Equations (3) and (4), it has c s 2 = 1 + w ϕ 2 . Meanwhile, by considering the stability of solutions with respect to inhomogeneous perturbations as 1 c s 2 0 , it constrains the range 1 w ϕ 1 . It follows that the k-essence model in this paper does not permit phantom behaviour.
In the following discussion, we neglect baryonic matter ρ b and the radiation ρ r in the matter component. Subsequently, the Friedmann equations take the form
H 2 = 1 3 ( ρ m + ρ ϕ ) ,
H ˙ = 1 2 ( ρ m + ρ ϕ + p ϕ ) ,
where H = a ˙ / a is the Hubble parameter, ρ ϕ and ρ m are the DE and DM density, respectively. The equation of motion for the k-essence field is given by
( F X + 2 X F X X ) ϕ ¨ + 3 H F X ϕ ˙ + ( 2 X F X F ) V ϕ V = 0 ,
where V ϕ d V / d ϕ . Equations (5) and (6) are usually transformed into an autonomous dynamical system when performing the phase-space and stability analysis. Being derived from the Friedmann eqs., we obtain a ¨ a = 1 6 ( 3 p + ρ ) , which implies a continuous eq. ρ ˙ + 3 H ( p + ρ ) = 0 . Because, in this model, the density is composed by two parts, the dark energy density and the matter density, i.e., ρ = ρ m + ρ ϕ , with the interaction Q between DM and DE, ρ ϕ and ρ m do not separately satisfy conservation laws. Subsequently, the following two equations are conceived as:
ρ ˙ ϕ + 3 H ( ρ ϕ + p ϕ ) = Q ,
ρ ˙ m + 3 H ( ρ m + p m ) = Q .
Here, p m = 0 . For Q < 0 , there is a transfer of energy from dark energy to dark matter. The case of Q = 0 , as in no interaction, was discussed in the former paper [38]. While, in this paper, the interaction is chosen by Q = α H ρ m , which means that the transformation between dark energy and dark matter happen, to some degree, in that circumstance. By Equations (8) and (9), to keep the physical dimensionality, we have to set α as a dimensionless parameter [33].
By setting the phase components, the auxiliary variables are defined as
x = ϕ ˙ , y = V ( ϕ ) 3 H ,
in order to transform the cosmological Equations (5) and (6) into an autonomous dynamical system { x , y } , by considering (8) and (9), where the prime is the derivative with respect to N = l n a . Subsequently, after solving the eqs. { x c = 0 , y c = 0 } , the critical points X c = { x c , y c } are obtained. To discuss the stability of each critical point, we expand X = { x , y } around the critical points X c = { x c , y c } by setting { x , y } T = { x c , y c } T + U with the perturbational variables U (see, for example, Refs. [41,54,55,56,57,58]). Up to the first order we acquire U = M · U with the 2 × 2 matrix M determined by
M = x x x y y x y y .
The matrix M contains the coefficients of the perturbation equations, and thus its eigenvalues determine the stability of the critical points. In this 2-dim system, which has two eigenvalues of M, for hyperbolic critical points, all of the eigenvalues have real parts that are different from zero: sink for the negative real parts is stable, saddle for real parts of different sign is unstable, and source for positive real parts is unstable. However, for the convenience of calculation in this model, an alternative way to judge the stability of the critical points in a 2-dim system is given by the trace tr M < 0 and determinant det M > 0 .
For the more general linear interaction form Q = α H ρ m + β H ρ ϕ , the autonomous dynamical equations are derived, as follows:
x = 3 2 λ x 2 y 3 x + 3 2 2 + α x y 2 1 2 x α + x β 2 ,
y = 1 4 y 2 3 λ x y + 6 3 2 x y 2 + 3 x 2 y 2 ,
for x > 0 ; while, for x < 0 , the equations turn out to be
x = 3 2 λ x 2 y 3 x 3 2 2 + α x y 2 1 2 x α + x β 2 ,
y = 1 4 y 2 3 λ x y + 6 + 3 2 x y 2 + 3 x 2 y 2 ,
where λ V ϕ / V 3 2 = c o n s t , by assuming V ( ϕ ) = 4 λ 2 ϕ 2 in this paper. Because it is hard to obtain an analytic solution as a critical point from the equations above with three parameters λ , α , and β , it has to be simplified by two cases in the next step, i.e., Q 1 = α H ρ m and Q 2 = β H ρ ϕ , respectively. However, we only analyze Q 1 in this paper due to the unsatisfactory for Q 2 . Besides that, the values of α and β are constrained in certain ranges, which are imposed for the observations [18,19,20].

3. The Analysis of Stability for This Dark Energy Model with Interaction Q = α H ρ m

In this case of Q = α H ρ m , when x > 0 , Equations (5) and (6) are simplified, as follows:
x = 3 2 λ x 2 y 3 x + 3 2 2 + α x y 2 1 2 x α ,
y = 1 4 y 2 3 λ x y + 6 3 2 x y 2 + 3 x 2 y 2 ,
By (16) and (17), the corresponding critical points { x c , y c } are
P 1 = { x c , y c } = { 3 λ + 6 , 6 3 ( λ + 6 ) } ,
P 2 = { x c , y c } = { 3 λ + 6 , 6 3 ( λ + 6 ) } ,
P 3 = { x c , y c } = { 2 ( 9 + 6 α + α 2 ) 2 λ 2 α + 9 + α 2 + 6 α , 6 6 λ ( α + 3 ) ( 2 λ 2 α + 9 + α 2 + 6 α ) } .
While, when x < 0 , we have
x = 3 2 λ x 2 y 3 x 3 2 2 + α x y 2 1 2 x α ,
y = 1 4 y 2 3 λ x y + 6 + 3 2 x y 2 + 3 x 2 y 2 ,
Additionally, there are other three solutions as the critical points, as follows:
P 4 = { x c , y c } = { 3 λ + 6 , 6 3 ( λ + 6 ) } ,
P 5 = { x c , y c } = { 3 λ + 6 , 6 3 ( λ + 6 ) } ,
P 6 = { x c , y c } = { 2 ( 9 + 6 α + α 2 ) 2 λ 2 α + 9 + α 2 + 6 α , 6 6 λ ( α + 3 ) ( 2 λ 2 α + 9 + α 2 + 6 α ) } .
The corresponding density parameter, the EoS, the sound speed, and the deceleration parameter are reexpressed as, respectively,
Ω ϕ = 1 2 x 2 y 2 ,
w ϕ = 1 2 x 1 ,
c s 2 = 1 2 2 x 1 ,
q = 1 2 + 3 4 x 2 y 2 3 2 4 x y 2 .
Equations (16), (17), (21) and (22) form the self-autonomous dynamical systems, which are valid in the whole phase-space, not only at the critical points. The critical points { x c , y c } of the autonomous system are obtained by setting the left-hand sides of the equations to zero, namely by solving X = ( x , y ) T = 0 . Six critical points are obtained, as shown in Table 1, in which we also present the necessary conditions for their existences and stabilities, as well as the corresponding cosmological quantities, c s 2 , Ω ϕ , w ϕ , and q. With these cosmological quantities, we can investigate the final state of the universe and discuss whether there exists acceleration expansion or not. Physically, it requires 0 Ω ϕ 1 , so the auxiliary variables x and y are constrained as 0 1 2 x 2 y 2 1 . In order to comply with the accelerated expansion, it requires w ϕ < 1 3 . When considering the sound speed, it has to be 0 c s 2 1 . For the existence, it means y > 0 , and x < 0 or x > 0 for each case. The stability means det M > 0 and tr M < 0 , instead of the analysis by each eigenvalue.
For the power law k-essence dark energy presented in this paper, the specific expression of the M , det M , and tr M are as follows:
M = α x 2 y 2 1 2 α + 3 λ x y 3 2 α y 3 x + 3 2 λ x 2 1 4 y ( 2 3 λ y 3 2 y 2 + 6 x y 2 ) 3 λ x y + 3 2 9 2 4 x y 2 + 9 4 x 2 y 2 ,
det M = 9 4 λ 2 x 2 y 2 + 9 3 2 λ x y 15 6 8 λ x 2 y 3 + 3 3 2 λ x 3 y 3 9 2 + 27 2 4 x y 2 27 4 x 2 y 2 + 3 2 x y α λ + 9 2 8 x y 2 α 9 8 x 2 y 2 α 3 α 2 x 2 y 2 + 3 2 α 4 x ,
tr M = α x 2 y 2 α 2 3 2 9 2 4 x y 2 + 9 4 x 2 y 2 ,
for the case of x > 0 ; and,
M = α x 2 y 2 α 2 + 3 λ x y 3 2 α y 3 x + 3 2 λ x 2 1 4 y ( 2 3 λ y + 3 2 y 2 + 6 x y 2 ) 3 λ x y + 3 2 + 9 2 4 x y 2 + 9 4 x 2 y 2 ,
det M = 9 4 λ 2 x 2 y 2 + 9 3 2 λ x y + 15 6 8 λ x 2 y 3 + 3 3 2 λ x 3 y 3 9 2 27 2 4 x y 2 27 4 x 2 y 2 + 3 2 x y α λ 9 2 8 x y 2 α 9 8 x 2 y 2 α 3 α 2 x 2 y 2 3 2 α 4 x ,
tr M = α x 2 y 2 α 2 3 2 + 9 2 4 x y 2 + 9 4 x 2 y 2 ,
for the case of x < 0 .
According to the stability conditions of critical points by the determinant and trace, together with those cosmological quantities, we obtain the value range of λ and α in the parameter plane, which makes the critical points stable and causes accelerated expansion, as shown in Table 1. At first, P 2 and P 5 are excluded by using existence condition. Based on the range of parameters presented in Figure 1a, we plot the stable point P 1 and its evolutionary trajectory for λ = 0.5 and α = 0.01 in Figure 1b, as well as P 3 , P 4 , and P 6 (see Figure 2, Figure 3 and Figure 4) with some certain pairs of parameters λ and α , respectively, in order to have a visual understanding of the evolutionary behavior near critical points. Especially, P 3 and P 6 are spiral attractors, which have spiral evolutionary trajectories around them. Additionally, the evolutionary trajectories of the cosmological quantities are shown in Figure 5 and Figure 6. Below, we will analyze these stable points P 1 , P 3 , P 4 , and P 6 one-by-one.
For P 1 , it has Ω ϕ = 1 ; the universe will be dominated by k-essence dark energy. If λ = 0 , then the k-essence would behave like a cosmological constant. The deceleration parameter q = 1 + 6 2 λ . The final state of the universe depends on the potential parameter λ , i.e., the universe expansion would speed up if λ < 6 3 , it would expand at constant-speed if λ = 6 3 , and the universe expansion would slow down if 6 3 < λ .
For P 3 , Ω ϕ = ( 3 + α ) 2 6 λ 2 , the universe will be dominated by both k-essence dark energy and dark matter. When α and λ satisfy 2 α λ 2 ( 3 + α ) 2 = 1 , lying on the top edge of the grey region in Figure 2a, the k-essence will behave like a cosmological constant. The evolutionary trajectory in the phase space shown in Figure 2b will be spiral around P 3 , finally converging to the attractor P 3 . The deceleration parameter q = 1 2 + α 2 . The final state of the universe depends on the coupling parameter α , i.e., the expansion of the universe will speed up if 3 < α < 1 , will expand in a constant speed if α = 1 , and will slow down if 1 < α < 0 , respectively.
For x < 0 , the results are quite like the three critical points that are investigated above. Among the three critical points, P 4 and P 6 are attractors for some λ and α , as displayed in Figure 3 and Figure 4. Those two points are also physically meaningful; while, P 5 does not.
For P 4 , it has Ω ϕ = 1 , the universe will be dominated by k-essence dark energy. If λ = 0 , the k-essence will behave like a cosmological constant. The deceleration parameter q = 1 6 2 λ , which indicates that the final state of the universe depends on the potential, i.e., the universe expansion will speed up if λ > 6 3 , will expand with constant-speed if λ = 6 3 , and it will slow down if 6 3 > λ .
For P 6 , Ω ϕ = ( 3 + α ) 2 6 λ 2 , the universe will be dominated by both k-essence and dark matter. When α and λ satisfy 2 α λ 2 ( 3 + α ) 2 = 1 , the k-essence will behave like a cosmological constant. In the phase space, the evolutionary trajectory will be spiral around P 6 , and then finally converge to the attractor point. The deceleration parameter q = 1 2 + α 2 , which indicates that the final state of the universe depends on the dark matter: the universe will speed up if α < 1 , will expand with constant-speed if α = 1 , and it will slow down if α > 1 .
From Figure 5 and Figure 6, it is not difficult to see that the density parameter of DE Ω ϕ continuously varies, which indicates that there is the exchange of energy between dark energy and dark matter by the interaction Q, and, at present, the universe is composed by both DE and DM. Meanwhile, the deceleration parameter q evolves from positive to negative values, which shows that the universe experiences a decelerating expansion in the past, and then transforms to an accelerating expansion at present, and keeps on speeding up into the future. Especially, in Figure 5b and Figure 6b, the w ϕ > 0 happens in the early time, which means that, in early time, the dark energy performs in a relativistic matter, which provides positive pressure, and acts as the attraction force to enhance the structure formation. All four evolutional figures show the w ϕ < 1 / 3 in the late time universe, which indicates that this k-essence DE model could explain the accelerating expansion of the universe.

4. The Geometric Diagnostic of Statefinder Hierarchy

Because the Λ CDM model is the best fitting for observations until now, the statefinder pair is a way for distinguishing a certain model from the Λ CDM model, by showing the “distance” of trajectories in the s r plane from the spatially flat Λ CDM model scenario which is a fixed point { s , r } | Λ C D M = { 0 , 1 } . Beyond the Hubble parameter H = a ˙ a and the deceleration parameter q = a a ¨ a ˙ 2 , the third order derivative r = a a H 3 , together with a combination of r and q, which is s = r 1 3 ( q 1 2 ) , become the cosmological diagnostic pair { r , s } . In terms of Ω ϕ and w, the statefinder pair has the following form:
r = 1 + 9 2 Ω ϕ w ( 1 + w ) 3 2 Ω ϕ w ˙ H ,
s = 1 + w 1 3 w ˙ w H .
Further, the statefinder hierarchy is an extension of the statefinder pair, which comes from the view point of higher derivatives of the expansion factor d n a / d t n , in Taylor expanded: ( 1 + z ) 1 = a ( t ) a 0 = 1 + A n ( t 0 ) n ! [ H 0 ( t t 0 ) ] n , where A n = a ( t ) ( n ) a ( t ) H n , n N , and a ( t ) ( n ) = d n a ( t ) / d t n . By using A n , the series S n and S n ( 1 ) are constructed, as follows. In a spatially flat universe with pressureless matter and a cosmological constant, such as Λ C D M model, A n could be expressed by parameter q or Ω m , where Ω m = 2 3 ( 1 + q ) , as follows:
A 2 = 1 3 2 Ω m ,
A 3 = 1 ,
A 4 = 1 3 2 2 Ω m ,
A 5 = 1 + 3 Ω m + 3 3 2 Ω m 2 .
Then the statefinder hierarchy S n is defined as
S 2 = A 2 + 3 2 Ω m ,
S 3 = A 3 ,
S 4 = A 4 + 3 2 2 Ω m ,
S 5 = A 5 3 Ω m 3 3 2 Ω m 2 .
Finally, it derived the null diagnostic for the Λ C D M model, the S n ( 1 ) , as follows:
S 3 ( 1 ) = S 3 ,
S 4 ( 1 ) = A 4 + 3 ( 1 + q ) ,
S 5 ( 1 ) = A 5 2 ( 4 + 3 q ) ( 1 + q ) .
For the Λ C D M model, it always has S n | Λ C D M = 1 and S n ( 1 ) | Λ C D M = 1 . In this paper, we focus on how parameters λ and α effect the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) for the k-essence model with coupling Q = α H ρ m . We can obtain the following expressions with interaction Q:
S 3 ( 1 ) = 1 + 9 2 Ω ϕ w ϕ ( 1 + w ϕ ) 3 2 Ω ϕ w ϕ 3 w ϕ Q 2 H ρ ,
S 4 ( 1 ) = 1 9 4 w ϕ Ω ϕ 2 [ 3 w ϕ ( 1 + w ϕ ) w ϕ ] 3 4 Ω ϕ [ w ϕ ( 21 + 39 w ϕ + 18 w ϕ 2 ) ( 13 + 18 w ϕ ) w ϕ + 2 w ϕ ] + 3 w ϕ Q 2 H ρ ( 2 + 3 w ϕ ) 3 w ϕ Q 2 H ρ 3 w ϕ Q H ρ .
The former methods for analyzing the statefinder pair { r , s } or the statefinder hierarchy are mainly around purely kinetic k-essence dark energy models, by finding the relations among F X , w ϕ and a ( t ) . For the potential, in this paper, it is not as a constant as in PKK and the analytic relation cannot be derived directly; instead, the statefinder hierarchy should be represented by phase components { x , y } . Along with the numerical method presented in [49], after substituting (16) and (17) into (49) and (50) for x > 0 , and (21) and (22) into (49) and (50) for x < 0 , respectively, the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) could be expressed by variables x, y with parameters λ , α . We choose the initial points { x 0 , y 0 } around the attractors separately, and then adjust λ with a constant α , whereas adjust α with a constant λ , in order to see the effect on statefinder hierarchy for each by the figures.
For P 1 , from Figure 7, by changing the value of potential parameter λ with a fixed coupling parameter α = 0.01 , the difference of both S 3 ( 1 ) and S 4 ( 1 ) evolutionary curves are obvious. However, in Figure 8, there is nearly no difference in both S 3 ( 1 ) and S 4 ( 1 ) cases, with the same value of λ = 0.5 , but different values of α . It means that the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) are more sensible to the potential parameter λ than α around P 1 .
Oppositely, for P 3 , there is no difference by changing the potential parameter λ under the same value α = 0.02 for both S 3 ( 1 ) and S 4 ( 1 ) in Figure 9, while the statefinder hierarchy shows the sensibility to the coupling parameter α under the same value of λ = 11 in Figure 10. That is to say, curves with same value of α perform alike, while λ makes little effect around P 3 .
For the case of x < 0 , i.e., P 4 and P 6 , the evolutionary curves of the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) are analogous to the P 1 and P 3 , respectively. For all of cases above, the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) of the k-essence DE model in this paper can tell the difference from the Λ CDM model, which is a straight line in the figures.

5. Conclusions

In summary, we have deeply investigated the cosmological evolution, the dynamical stability, as well as the geometrical diagnostic of the power law k-essence dark energy model with the Lagrangian containing a kinetic function F ( X ) = X + X and interaction Q = α H ρ m in FRWL space time. Concretely, we have not only discussed the influences of the coupling parameter α and potential parameter λ on the evolution of several cosmological quantities (such as the density parameter Ω ϕ , EoS of dark energy w ϕ , the effective sound speed c s 2 , and deceleration parameter q), but also numerically analyzed the dynamical stability and showed that there are the four dynamical attractors in the phase space. In addition, the statefinder hierarchy S 3 ( 1 ) and S 4 ( 1 ) of this dark energy model have been numerically obtained, which shows that the potential parameter λ has more influence on the evolutions of S 3 ( 1 ) and S 4 ( 1 ) than one of the coupling parameter α around P 1 and P 4 ; while, for P 3 and P 6 , parameter α plays a more important role in S 3 ( 1 ) and S 4 ( 1 ) .

Author Contributions

B.-H.C.: search the references, plot the figures, write the manuscript; Y.-B.W.: corresponding author, review the manuscript, revise the manuscript; D.-F.X.: plot the figures, arrange the references, edit the manuscript; W.D.: plot the figures, explain the data, search the references; N.Z.: design the research plan, plot the figures, explain the data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant number 12075109, 11575075, 11705079 and 11865012.

Acknowledgments

Y. Wu has been supported by the National Natural Science Foundation of China under grants No.12075109 and No.11575075. W. Yang has been supported by the National Natural Science Foundation of China under grant No.11705079. J. Lu has been supported by the National Natural Science Foundation of China under grants No.11865012.

Conflicts of Interest

There is no conflict of interest.

References

  1. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef] [Green Version]
  2. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. [The Supernova Cosmology Project] Measurements of Ω and Λ from 42 High-Redshift Supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  3. Spergel, D.N.; Verde, L.; Peiris, H.V.; Komatsu, E.; Nolta, M.R.; Bennett, C.L.; Halpern, M.; Hinshaw, G.; Jarosik, N.; Kogut, A.; et al. First-year Wilkinson Microwave Anisotropy Probe (WMAP)* observations: Determination of cosmological parameters. Astrophys. J. Suppl. 2003, 148, 175. [Google Scholar] [CrossRef] [Green Version]
  4. Eisenstein, D.J.; Zehavi, I.; Hogg, D.W.; Scoccimarro, R.; Blanton, M.R.; Nichol, R.C.; Scranton, R.; Seo, H.J.; Tegmark, M.; Zheng, Z.; et al. Detection of the baryon acoustic peak in the large-scale correlation function of SDSS luminous red galaxies. Astrophys. J. 2005, 633, 560. [Google Scholar] [CrossRef]
  5. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. [Planck Collaboration] Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. J. 2020, 641, 1–67. [Google Scholar]
  6. Xu, L.X.; Lu, J.B.; Wang, Y.T. Revisiting generalized Chaplygin gas as a unified dark matter and dark energy model. Eur. Phys. J. C 2012, 72, 1883. [Google Scholar] [CrossRef]
  7. Xu, L.X.; Wang, Y.T.; Noh, H. Modified Chaplygin gas as a unified dark matter and dark energy model and cosmic constraints. Eur. Phys. J. C 2012, 72, 1931. [Google Scholar] [CrossRef]
  8. Yang, W.Q.; Li, H.; Wu, Y.B.; Lu, J.B. Cosmological implications of the dark matter equation of state. Int. J. Mod. Phys. D 2017, 26, 1750013. [Google Scholar] [CrossRef]
  9. Du, M.H.; Yang, W.Q.; Xu, L.X.; Pan, S.; Mota, D.F. Future constraints on dynamical dark-energy using gravitational-wave standard sirens. Phys. Rev. D 2019, 100, 043535. [Google Scholar] [CrossRef] [Green Version]
  10. Sahni, V.; Starobinsky, A.A. The case for a positive cosmological Λ-term. Int. J. Mod. Phys. D 2000, 9, 373. [Google Scholar] [CrossRef]
  11. Carroll, S.M. The cosmological constant. Living Rev. Rel. 2001, 4, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Carroll, S.M.; Hoffman, M.; Trodden, M. Can the dark energy equation-of-state parameter w be less than −1? Phys. Rev. D 2003, 68, 023509. [Google Scholar] [CrossRef] [Green Version]
  13. Li, M.; Li, X.D.; Wang, S.; Wang, Y. Commun. Dark Energy. Theor. Phys. 2011, 56, 525–604. [Google Scholar]
  14. Yang, R.J.; Zhang, S.N. The age problem in the ΛCDM model. Mon. Not. R. Astron. Soc. 2010, 407, 1835–1841. [Google Scholar] [CrossRef] [Green Version]
  15. Valentino, E.D.; Melchiorri, A.; Silk, J. Reconciling Planck with the local value of H0 in extended parameter space. Phys. Lett. B 2016, 761, 242. [Google Scholar] [CrossRef]
  16. Huang, Q.G.; Wang, K. How the dark energy can reconcile Planck with local determination of the Hubble constant. Eur. Phys. J. C 2016, 76, 506. [Google Scholar] [CrossRef] [Green Version]
  17. Vagnozzi, S. New physics in light of the H0 tension: An alternative view. Phys. Rev. D 2020, 102, 023518. [Google Scholar] [CrossRef]
  18. Kumar, S.; Nunes, R.C. Probing the interaction between dark matter and dark energy in the presence of massive neutrinos. Phys. Rev. D 2016, 94, 123511. [Google Scholar] [CrossRef] [Green Version]
  19. Pan, S.; Yang, W.Q.; Singha, C.; Saridakis, E.N. Observational constraints on sign-changeable interaction models and alleviation of the H0 tension. Phys. Rev. D 2019, 100, 083539. [Google Scholar] [CrossRef] [Green Version]
  20. Pan, S.; Yang, W.Q.; Valentino, E.D.; Saridakis, E.N.; Chakraborty, S. Interacting scenarios with dynamical dark energy: Observational constraints and alleviation of the H0 tension. Phys. Rev. D 2019, 100, 103520. [Google Scholar] [CrossRef] [Green Version]
  21. Tsujikawa, S. Quintessence: A Review. Class. Quant. Grav. 2013, 30, 214003. [Google Scholar] [CrossRef] [Green Version]
  22. Feng, B.; Wang, X.L.; Zhang, X.M. Dark Energy Constraints from the Cosmic Age and Supernova. Phys. Lett. B 2005, 607, 35–41. [Google Scholar] [CrossRef]
  23. Armendariz-Picon, C.; Damour, T.; Mukhanov, V. K-inflation. Phys. Lett. B 1999, 458, 209–218. [Google Scholar] [CrossRef] [Green Version]
  24. Armendariz-Picon, C.; Mukhanov, V.; Steinhardt, P.J. A Dynamical Solution to the Problem of a Small Cosmological Constant and Late-time Cosmic Acceleration. Phys. Rev. Lett. 2000, 85, 4438–4441. [Google Scholar] [CrossRef] [Green Version]
  25. Rendall, A.D. Dynamics of k-essence. Class. Quant. Grav. 2006, 23, 1557–1570. [Google Scholar] [CrossRef] [Green Version]
  26. Scherrer, R.J. Purely kinetic k-essence as unified dark matter. Phys. Rev. Lett. 2004, 93, 011301. [Google Scholar] [CrossRef] [Green Version]
  27. Bose, N.; Majumdar, A.S. A k-essence Model of Inflation, Dark Matter and Dark Energy. Phys. Rev. D 2009, 79, 103517. [Google Scholar] [CrossRef] [Green Version]
  28. Bose, N.; Majumdar, A.S. Unified Model of k-Inflation, Dark Matter and Dark Energy. Phys. Rev. D 2009, 80, 103508. [Google Scholar] [CrossRef] [Green Version]
  29. Chiba, T.; Okabe, T.; Yamaguchi, M. Kinetically Driven Quintessence. Phys. Rev. D 2000, 62, 023511. [Google Scholar] [CrossRef] [Green Version]
  30. Chimento, L.P.; Feinstein, A. Power-law expansion in k-essence cosmology. Mod. Phys. Lett. A 2004, 19, 761–768. [Google Scholar] [CrossRef]
  31. Chimento, L.P. Extended tachyon field, Chaplygin gas and solvable k-essence cosmologies. Phys. Rev. D 2004, 69, 123517. [Google Scholar] [CrossRef] [Green Version]
  32. Battye, R.A.; Pace, F. Approximation of the potential in scalar field dark energy models. Phys. Rev. D 2016, 94, 063513. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, B. Dark Matter and Dark Energy Interactions: Theoretical Challenges, Cosmological Implications and Observational Signatures. Rep. Prog. Phys. 2016, 79, 096901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Chimento, L.P. Linear and nonlinear interactions in the dark sector. Phys. Rev. D. 2010, 81, 043525. [Google Scholar] [CrossRef] [Green Version]
  35. Paliathanasis, A.; Pan, S.; Yang, W.Q. Dynamics of nonlinear interacting dark energy models. Int. J. Mod. Phys. D 2019, 28, 1950161. [Google Scholar] [CrossRef] [Green Version]
  36. Yang, W.Q.; Pan, S.; Barrow, J.D. Large-scale stability and astronomical constraints for coupled dark-energy models. Phys. Rev. D. 2018, 97, 043529. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, N.; Wu, Y.B.; Chi, J.N.; Yu, Z.; Xu, D.F. Diagnosing Tsallis holographic dark energy models with interactions. Mod. Phys. Lett. A 2020, 33, 2050044. [Google Scholar] [CrossRef] [Green Version]
  38. Yang, R.J.; Chen, B.H.; Li, J.; Qi, J.Z. The evolution of the power law k-essence cosmology. Astrophys. Space Sci. 2015, 356, 399–405. [Google Scholar] [CrossRef] [Green Version]
  39. Copeland, E.J.; Sami, M.; Tsujikawa, S. Dynamics of dark energy. Int. J. Mod. Phys. D. 2006, 15, 1753–1935. [Google Scholar] [CrossRef] [Green Version]
  40. Bahamonde, S.; Boehmer, C.G.; Carloni, S.; Copeland, E.J.; Fang, W.; Tamanini, N. Dynamical systems applied to cosmology: Dark energy and modified gravity. Phys. Rep. 2018, 775, 1–122. [Google Scholar] [CrossRef] [Green Version]
  41. Leon, G.; Saridakis, E.N. Phase-space analysis of Horava-Lifshitz cosmology. J. Cosmol. Astropart. Phys. 2009, 2009, 006. [Google Scholar] [CrossRef]
  42. Sahni, V.; Shafieloo, A.; Starobinsky, A.A. Two new diagnostics of dark energy. Phys. Rev. D 2008, 78, 103502. [Google Scholar] [CrossRef] [Green Version]
  43. Wu, Y.B.; Zhang, C.Y.; Lu, J.B.; Lu, J.W.; Zhang, X.; Qi, Q.M.; Fan, B. Analysis on modified Chaplygin gas as dark energy model. Mod. Phys. Lett. A 2015, 30, 1550005. [Google Scholar] [CrossRef]
  44. Sahni, V.; Saini, T.D.; Starobinsky, A.A.; Alam, U. Statefinder—A new geometrical diagnostic of dark energy. JETP Lett. 2003, 77, 201–206. [Google Scholar] [CrossRef]
  45. Alam, U.; Sahni, V.; Saini, T.D.; Starobinsky, A.A. Exploring the Expanding Universe and Dark Energy using the Statefinder Diagnostic. Mon. Not. R. Asron. Soc. 2003, 344, 1057. [Google Scholar] [CrossRef]
  46. Arabsalmani, M.; Sahni, V. The Statefinder hierarchy: An extended null diagnostic for concordance cosmology. Phy. Rev. D. 2011, 83, 043501. [Google Scholar] [CrossRef] [Green Version]
  47. Li, J.; Yang, R.J.; Chen, B.H. Discriminating dark energy models by using the statefinder hierarchy and the growth rate of matter perturbations. J. Cosmol. Astropart. Phys. 2014, 2014, 043. [Google Scholar] [CrossRef] [Green Version]
  48. Cui, J.L.; Yin, L.; Wang, L.F.; Li, Y.H.; Zhang, X. A closer look at interacting dark energy with statefinder hierarchy and growth rate of structure. J. Cosmol. Astropart. Phys. 2014, 2014, 024. [Google Scholar] [CrossRef]
  49. Liu, W.Z.; Liu, D.J. Statefinder Diagnostic for Quintessence with or without Thermal Interaction. Int. J. Mod. Phys. D 2009, 18, 43–52. [Google Scholar] [CrossRef] [Green Version]
  50. De-Santiago, J.; Cervantes-Cota, J.L. Generalizing a Unified Model of Dark Matter, Dark Energy, and Inflation with Non Canonical Kinetic Term. Phys. Rev. D 2011, 83, 063502. [Google Scholar] [CrossRef] [Green Version]
  51. Garriga, J.; Mukhanov, V.F. Perturbations in k-inflation. Phys. Lett. B 1999, 458, 219–225. [Google Scholar] [CrossRef] [Green Version]
  52. Unnikrishnan, S.; Sriramkumar, L. A note on perfect scalar fields. Phys. Rev. D 2010, 81, 103511. [Google Scholar] [CrossRef] [Green Version]
  53. Christopherson, A.J.; Malik, K.A. The non-adiabatic pressure in general scalar field system. Phys. Lett. B 2009, 675, 159–163. [Google Scholar] [CrossRef] [Green Version]
  54. Copeland, E.J.; Liddle, A.R.; Wands, D. Exponential potentials and cosmological scaling solutions. Phys. Rev. D 1998, 57, 4686–4690. [Google Scholar] [CrossRef] [Green Version]
  55. Yang, R.J.; Gao, X.T. Phase-space analysis of a class of k-essence cosmology. Class. Quantum Gravity 2011, 28, 065012. [Google Scholar] [CrossRef] [Green Version]
  56. Capozziello, S.; Nojiri, S.; Odintsov, S.D. Unified phantom cosmology: Inflation, dark energy and dark matter under the same standard. Phys. Lett. B 2006, 632, 597–604. [Google Scholar] [CrossRef] [Green Version]
  57. Carloni, S.; Elizalde, E.; Silva, P.J. An analysis of the phase space of Horava-Lifshitz cosmologies. Class. Quant. Grav. 2010, 27, 045004. [Google Scholar] [CrossRef]
  58. He, J.; Wu, Y.B.; Fu, M.H. Dynamical Attractor of Modified Chaplygin Gas. Chin. Phys. Lett. 2008, 25, 347. [Google Scholar]
Figure 1. (a) The value ranges for parameters λ and α to make the critical point P 1 exist and stable, which is also constrained by other cosmological quantities, when x > 0 . (b) The phase plane for λ = 0.5 and α = 0.01 around the attractor P 1 = ( 0.888 , 1.592 ) when x > 0 .
Figure 1. (a) The value ranges for parameters λ and α to make the critical point P 1 exist and stable, which is also constrained by other cosmological quantities, when x > 0 . (b) The phase plane for λ = 0.5 and α = 0.01 around the attractor P 1 = ( 0.888 , 1.592 ) when x > 0 .
Universe 06 00244 g001
Figure 2. (a) The value range for parameters λ and α to make the critical point P 3 exist and stable, also constrained by other cosmological quantities, when x > 0 . (b) The phase plane for λ = 11 and α = 0.02 around the attractor P 3 = ( 0.915 , 0.171 ) when x > 0 .
Figure 2. (a) The value range for parameters λ and α to make the critical point P 3 exist and stable, also constrained by other cosmological quantities, when x > 0 . (b) The phase plane for λ = 11 and α = 0.02 around the attractor P 3 = ( 0.915 , 0.171 ) when x > 0 .
Universe 06 00244 g002
Figure 3. (a) The value range for parameters λ and α to make the critical point P 4 exist and stable, also being constrained by other cosmological quantities, when x < 0 . (b) The phase plane for λ = 0.5 and α = 0.01 around the attractor P 4 = ( 0.888 , 1.592 ) when x < 0 .
Figure 3. (a) The value range for parameters λ and α to make the critical point P 4 exist and stable, also being constrained by other cosmological quantities, when x < 0 . (b) The phase plane for λ = 0.5 and α = 0.01 around the attractor P 4 = ( 0.888 , 1.592 ) when x < 0 .
Universe 06 00244 g003
Figure 4. (a) The value range for parameters λ and α to make the critical point P 6 exist and stable, also constrained by other cosmological quantities, when x < 0 . (b) The phase plane for λ = 11 and α = 0.02 around the attractor P 6 = ( 0.915 , 0.171 ) when x < 0 .
Figure 4. (a) The value range for parameters λ and α to make the critical point P 6 exist and stable, also constrained by other cosmological quantities, when x < 0 . (b) The phase plane for λ = 11 and α = 0.02 around the attractor P 6 = ( 0.915 , 0.171 ) when x < 0 .
Universe 06 00244 g004
Figure 5. (a) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.85 , 1.5 ) around P 1 corresponding to ( λ , α ) = ( 0.5 , 0.01 ) . (b) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.9 , 1.0 ) around P 3 corresponding to ( λ , α ) = ( 1.2 , 0.8 ) .
Figure 5. (a) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.85 , 1.5 ) around P 1 corresponding to ( λ , α ) = ( 0.5 , 0.01 ) . (b) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.9 , 1.0 ) around P 3 corresponding to ( λ , α ) = ( 1.2 , 0.8 ) .
Universe 06 00244 g005
Figure 6. (a) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.9 , 1.55 ) around P 4 corresponding to ( λ , α ) = ( 0.5 , 0.01 ) . (b) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.75 , 1.25 ) around P 6 corresponding to ( λ , α ) = ( 1 , 0.9 ) .
Figure 6. (a) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.9 , 1.55 ) around P 4 corresponding to ( λ , α ) = ( 0.5 , 0.01 ) . (b) The evolutions of Ω ϕ , w ϕ and the deceleration parameter q by the initial condition ( x 0 , y 0 ) = ( 0.75 , 1.25 ) around P 6 corresponding to ( λ , α ) = ( 1 , 0.9 ) .
Universe 06 00244 g006
Figure 7. (a) Graph of S 3 ( 1 ) for P 1 with different parameters λ and fixed α = 0.01 . (b) Graph of S 4 ( 1 ) for P 1 with different parameters λ and fixed α = 0.01 . The initial point is x 0 = 1 , y 0 = 1.2 .
Figure 7. (a) Graph of S 3 ( 1 ) for P 1 with different parameters λ and fixed α = 0.01 . (b) Graph of S 4 ( 1 ) for P 1 with different parameters λ and fixed α = 0.01 . The initial point is x 0 = 1 , y 0 = 1.2 .
Universe 06 00244 g007
Figure 8. (a) Graph of S 3 ( 1 ) for P 1 with different parameters α and fixed λ = 0.5 . (b) Graph of S 3 ( 1 ) for P 1 with different parameters α and fixed λ = 0.5 . The initial point is x 0 = 1 , y 0 = 1.2 .
Figure 8. (a) Graph of S 3 ( 1 ) for P 1 with different parameters α and fixed λ = 0.5 . (b) Graph of S 3 ( 1 ) for P 1 with different parameters α and fixed λ = 0.5 . The initial point is x 0 = 1 , y 0 = 1.2 .
Universe 06 00244 g008
Figure 9. (a) Graph of S 3 ( 1 ) for P 3 with different parameters λ and fixed α = 0.02 . (b) Graph of S 4 ( 1 ) for P 3 with different parameters λ and fixed α = 0.02 . The initial point is x 0 = 1 , y 0 = 0.3 .
Figure 9. (a) Graph of S 3 ( 1 ) for P 3 with different parameters λ and fixed α = 0.02 . (b) Graph of S 4 ( 1 ) for P 3 with different parameters λ and fixed α = 0.02 . The initial point is x 0 = 1 , y 0 = 0.3 .
Universe 06 00244 g009
Figure 10. (a) Graph of S 3 ( 1 ) for P 3 with different parameters α and fixed λ = 11 . (b) Graph of S 4 ( 1 ) for P 3 with different parameters α and fixed λ = 11 . The initial point is x 0 = 1 , y 0 = 0.3 .
Figure 10. (a) Graph of S 3 ( 1 ) for P 3 with different parameters α and fixed λ = 11 . (b) Graph of S 4 ( 1 ) for P 3 with different parameters α and fixed λ = 11 . The initial point is x 0 = 1 , y 0 = 0.3 .
Universe 06 00244 g010
Table 1. The existence and stability conditions for six critical points, and the cosmological quantities in form of the parameters λ and α in each critical point.
Table 1. The existence and stability conditions for six critical points, and the cosmological quantities in form of the parameters λ and α in each critical point.
NameExistenceStability 0 c s 2 1 0 Ω ϕ 1 ω ϕ < 1 3 q
P 1 λ < 6 6 λ 3 α < 0   and 3 6 λ 2 α 12 < 0 λ 6 1 6 λ 3 3 6 2 λ 1
P 2 none 6 λ + 3 + α < 0   and 3 6 λ + 2 α + 12 > 0 λ 6 1 6 λ + 3 3 1 6 2 λ
P 3 ( 3 + α ) 2 > 2 λ 2 α and   λ ( α + 3 ) > 0 λ ( ( 3 + α ) 2 6 λ 2 ) < 0   and 6 λ 2 α + 9 α 2 + 2 α 3 < 27 α λ 2 ( 3 + α ) 2 + 1 2 ( 3 + α ) 2 6 λ 2 2 α λ 2 ( 3 + α ) 2 1 2 + α 2
P 4 λ > 6 6 λ + 3 + α > 0   and 3 6 λ + 2 α + 12 > 0 λ 6 1 6 λ + 3 3 1 6 2 λ
P 5 none ( λ 3 6 α 6 ) ( 6 λ 6 ) < 0 and   3 6 λ 2 α 12 < 0 6 6 λ 1 6 λ 3 3 6 2 λ 1
P 6 ( 3 + α ) 2 > 2 λ 2 α and   λ ( α + 3 ) < 0 λ ( ( 3 + α ) 2 6 λ 2 ) > 0   and 6 λ 2 α + 9 α 2 + 2 α 3 < 27 α λ 2 ( 3 + α ) 2 + 1 2 ( 3 + α ) 2 6 λ 2 2 α λ 2 ( 3 + α ) 2 1 2 + α 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, B.-H.; Wu, Y.-B.; Xu, D.-F.; Dong, W.; Zhang, N. Dynamical Stability and Geometrical Diagnostic of the Power Law K-Essence Dark Energy Model with Interaction. Universe 2020, 6, 244. https://doi.org/10.3390/universe6120244

AMA Style

Chen B-H, Wu Y-B, Xu D-F, Dong W, Zhang N. Dynamical Stability and Geometrical Diagnostic of the Power Law K-Essence Dark Energy Model with Interaction. Universe. 2020; 6(12):244. https://doi.org/10.3390/universe6120244

Chicago/Turabian Style

Chen, Bo-Hai, Ya-Bo Wu, Dong-Fang Xu, Wei Dong, and Nan Zhang. 2020. "Dynamical Stability and Geometrical Diagnostic of the Power Law K-Essence Dark Energy Model with Interaction" Universe 6, no. 12: 244. https://doi.org/10.3390/universe6120244

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