Next Article in Journal
Outliers in Spectral Time Lag-Selected Gamma Ray Bursts
Next Article in Special Issue
Study of the Inflationary Spectrum in the Presence of Quantum Gravity Corrections
Previous Article in Journal
Exploring Anisotropic Lorentz Invariance Violation from the Spectral-Lag Transitions of Gamma-Ray Bursts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Phase Space Analysis of Interacting K-Essence Dark Energy Models in Loop Quantum Cosmology

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
*
Author to whom correspondence should be addressed.
Universe 2022, 8(10), 520; https://doi.org/10.3390/universe8100520
Submission received: 17 August 2022 / Revised: 23 September 2022 / Accepted: 28 September 2022 / Published: 6 October 2022
(This article belongs to the Special Issue Recent Advances in Quantum Cosmology)

Abstract

:
The present work deals with two kinds of k-essence dark energy models within the framework of loop quantum cosmology (LQC). The two kinds of k-essence models originates from two forms of Lagrangians, i.e., L 1 = F ( X ) V ( ϕ ) and L 2 = F ( X ) V ( ϕ ) , where F ( X ) and V ( ϕ ) stand for the kinetic term and potential of the scalar field ϕ , respectively. Two models are based on different phase variables settings, and the general form of autonomous dynamical system is deduced for each Lagrangian. Then, the dynamical stabilities of the critical points in each model are analysed in different forms of F ( X ) and V ( ϕ ) . Model I is a 3-dim system with four stable points, and Model II is a 4-dim system but reduced to a 3-dim system using the symmetry analysis, which has five stable points. Moreover, the corresponding cosmological quantities, such as Ω ϕ , w ϕ and q, are calculated at each critical point. To compare these with the case of the classical Einstein cosmology (EC), the dynamical evolutionary trajectories in the phase space and evolutionary curves of the cosmological quantities are drawn for both EC and LQC cases, which shows that the loop quantum gravity effects diminish in the late-time universe but are significant in the early time. Further, the effects of interaction Q = α H ρ m on the evolutions of the universe are discussed. With the loop quantum gravity effects, bouncing universe is achieved in both models for different initial values of ϕ 0 , ϕ ˙ 0 , H 0 , ρ 0 and coupling parameter α , which helps to avoid singularities. However, the interaction has little effect on bounce, although it is important to the stability of some critical points.

1. Introduction

As we know, the astronomical observations on the luminosity–redshift relation of distant Type Ia supernovas [1,2] in the late 90’s gave the world a striking surprise that the universe is experiencing an accelerated expansion. This mysterious power is the so-called dark energy (DE), which has the feature of negative pressure and its equation of state (EoS) w D E = p D E ρ D E < 1 / 3 , where p D E is its pressure and ρ D E is the density. It motivates the further observations and theoretical researches. Over the last twenty years, with the enhancement of observational accuracy [3,4,5], measurements indicate that the current energy density in the universe is composed of 68.3% dark energy, 26.8% dark matter (DM) and 4.9% baryons; meanwhile, the relative constraints on cosmological quantities were developed from different aspects [6,7,8,9]. Until now, the Λ CDM model is still a decent explanation to the accelerated expansion, which has a constant equation of state (EoS) w Λ = 1 as a candidate of DE in accordance with the observations. However, the Λ CDM model is not so perfect to unify all problems, such as reconciling the cosmological constant problem, the age problem and the coincidence problem [10,11,12,13,14]. Then, phenomenally a vast of scalar field models have been researched to replace the cosmological constant Λ , while generally in those models pressure and energy density are expressed by certain forms of kinetic and potential in terms of a scalar and its derivative to time, such as, quintessence, phantom, k-essence and Dirac-Born-Infeld models, etc. [15,16,17,18,19,20,21].
However, all the above models face the singularity problem, such as phantom field DE in classical Einstein cosmology, which leads to the Big Rip singularity for the fluids with w 1 [22]. Following the classical Einstein cosmology (EC), based on quantum geometry, loop quantum gravity (LQG) was developed as a background-independent, non-perturbative approach to quantum gravity, and has been applied as loop quantum cosmology (LQC), which adds a modified term ( 1 ρ ρ l c ) to the classical Friedmann equation [23], as well as some other forms, as in [24,25]. Within the framework of LQC, this provides us with a novel means of researching the big bang singularity and the black hole singularity [26,27]. In LQC, the future singularities in the classical phantom models are avoided, instead, by a recollapsing universe, where a bounce universe in an oscillatory regime occurs when the total energy density reaches a Plankton value [28,29,30,31,32,33,34]. In the case of constant EoS of DE [35,36], the bounce appears periodically forever, and the Hubble rate H, with its derivative H ˙ , are bounded; in the case of a certain form of scalar field models with a variable EoS of DE [29,30,37], the bounce has an infinite frequency of oscillations at a finite time in future, which leads to a diverge of H ˙ , and then the curvature singularity occurs for R = 6 H ˙ + 12 H 2 . The dynamical properties of the phase space are also studied [35,37,38]. In recent years, [32] investigated the inflation in the LQC-modified scalar-tensor theory and [39,40] focused on the generalised EoS in LQC with interactions between DM and DE.
Originating from k-inflation [41,42,43,44], the k-essence model could be viewed as a generalization of the quintessence model, which has a canonical Lagrangian L = X V ( ϕ ) . The k-essence model provides a variety of non-canonical Lagrangian terms, i.e., L = F ( X , ϕ ) [45,46,47], such as L = F ( X ) V ( ϕ ) , L = F ( X ) V ( ϕ ) , L = V ( ϕ ) f ( F ( X ) V ( ϕ ) ) and so on, which are used as the pressure of the dark energy [48,49,50,51,52]. The action is described by a single scalar field ϕ , and a canonical kinetic energy X 1 2 μ ϕ μ ϕ . Some modified kinetic terms have been discussed, such as F ( X ) = K X + L X 2 + and F ( X ) = 1 2 α 1 ( ( A X ) α 2 α α 0 A X ) [45,48,49,53,54,55,56,57]. The potential V ( ϕ ) takes many forms, such as V ( ϕ ) ϕ n , V ( ϕ ) e k ϕ , V ( ϕ ) s i n h 1 ( k ϕ ) and so on [58].
From the matter clustering properties, dark matter (DM) and dark energy (DE) are not the same substance; however there is research on the interactions between them, even some nonlinear interaction forms [59,60,61,62,63,64,65,66,67], which could provide a mechanism by which to generate acceleration and alleviate the coincidence problem. A theoretical explanation of the transition of the DE Eos from w > 1 to w < 1 has also been presented. With the interaction between DM and DE, ρ ϕ and ρ m do not separately satisfy independent conservation laws. Furthermore, a cosmological evolution system with DE, DM and unparticle with three kinds of interactions are researched in [68].
In this work, since the k-essence DE models are not widely discussed in the frame of LQC, on the basis of our previous work [69], we research the dynamical stabilities of the phase space of two kinds k-essence DE models in the framework of LQC, especially interactions between DM and DE. The first model is of L 1 = F ( X ) V ( ϕ ) f ( ϕ ) , where F ( X ) = X + X , V ( ϕ ) 1 / ϕ 2 , and f ( ϕ ) = 0 [67,69], together with a certain kind of interaction Q. Another is of L 2 = F ( X ) V ( ϕ ) , where F ( X ) = X η , V ( ϕ ) ϕ n [50,51], together with interaction Q. We investigate the possible cosmological behavior of these models in Friedmann–Robertson–Walker–Lemaître (FRWL) spacetime by performing a phase-space and stability analysis. The theory is based on [70,71], judging the stability of the critical points by their eigenvalues, whereas, in this model, for the convenience of calculation, the method using the determinant and trace of the Jacobian matrix of the autonomous differential equations is used [72]. 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, and the deceleration parameter q, to compare these with the case in Einstein cosmology without a loop quantum effect and the effect of interaction Q. The effect of interaction Q will be shown in the dynamical system for each model in LQC. In addition, for each model, a periodical bounce universe will be formed by setting suitable initial values in the dynamical system, in which both the Hubble rate and its derivative H ˙ are bounded. That is to say, both the Big Rip singularity and the curvature singularity are avoided. Throughout this paper, we work with a flat, homogeneous, and isotropic FRWL spacetime with a signature ( , + , + , + ) and in units c = 8 π G = 1 .
This paper is organized as follows: in the next section, we will discuss the phase space analysis of k-essence dark energy models with Lagrangian L 1 = F ( X ) V ( ϕ ) in the LQC framework with interaction Q = α H ρ m , where six critical points will be discussed. In the third section, we will choose another Lagrangian L 2 = F ( X ) V ( ϕ ) in the model and proceed with the discussion using another set of phase variables, where ten critical points are discussed. For both models, the critical points in LQC are in accordance with the EC, which could be viewed as a generalization of EC [51,69]. However, compared with EC, in LQC the evolutions of the cosmological quantities are quite different in the early time universe; additionally, the evolution of the Hubble rate is bounded and oscillates forever. Finally, we close with a few concluding remarks in the fifth section.

2. Model I: L 1 = F ( X ) V ( ϕ ) in Loop Quantum Cosmology

In a flat universe, the effectively modified Friedmann equation in the framework of LQC is given by [29,30]
3 H 2 = ρ ( 1 ρ ρ l c ) ,
where H is the Hubble parameter, ρ = ρ m + ρ ϕ is the total energy density and ρ ϕ , ρ m are the energy densities of dark energy and dark matter, respectively. The constant
ρ l c = 3 16 π 2 γ 3 G 2 = 4 3 γ 3 ,
is the critical loop quantum density, where γ is the dimensionless Barbero–Immirzi parameter. Usually, ρ l c = 1.5 [29,35], and the LQC goes back to EC when ρ ρ l c or ρ l c . Considering the interaction between DM and DE, the conservative equation of the total energy density ρ ˙ + 3 H ( p + ρ ) = 0 is represented as
ρ ˙ ϕ + 3 H ( ρ ϕ + p ϕ ) = Q ,
ρ ˙ m + 3 H ( ρ m + p m ) = Q .
From (1), one can obtain the modified Raychaudhuri equation as
H ˙ = 1 2 ( ρ m + ρ ϕ + p ϕ ) ( 1 2 ρ ρ l c ) .
In this section, against the background of LQC, we consider a k-essence dark energy model with Lagrangian
L 1 = p ϕ = F ( X ) V ( ϕ ) ,
as the pressure of the scalar field, while the energy density ρ ϕ = V ( ϕ ) [ 2 X F X F ] . Then, the density parameter Ω ϕ , the EoS parameters w ϕ , w t o t and the deceleration parameter q are, respectively, given by
Ω ϕ = V 3 H 2 [ 2 X F X F ] ,
w ϕ = F 2 X F X F ,
w t o t = F V 3 H 2 ( 1 z ) 1 ,
q = 1 H ˙ H 2 ,
where F X d F / d X and F X X d 2 F / d X 2 . And according to F ( X ) = X + X , it has 2 X F X F = X .
The equation of motion for the k-essence field is given by
ϕ ¨ ϕ ˙ [ F X + 2 X F X X ] V + ϕ ˙ [ 2 X F X F ] V ϕ + 6 X F X V H = Q ,
where V ϕ d V / d ϕ . Equations (3) and (5) are usually transformed into an autonomous dynamical system when performing the phase-space and stability analysis. By setting the phase variables:
x = ϕ ˙ , y = V 3 H , z = ρ ρ l c , s = V ϕ V 3 / 2 ,
where 0 z 1 by Equation (1), and using Equation (12), the modified Friedmann Equation (1) and Raychaudhuri Equation (5) become
1 = ( Ω m + y 2 [ 2 X F X F ] ) ( 1 z ) ,
H ˙ = 1 2 ( ρ m + V [ 2 X F X F ] ) ( 1 2 z ) .
Then, the autonomous dynamical system of the phase variables is:
x = 1 F X + 2 X F X X ( 3 x F X + 3 y s [ 2 X F X F ] + Q x V H ) ,
y = 3 2 x y 2 s + ( 3 2 y 1 z + 3 2 y 3 F ) ( 1 2 z ) ,
z = 3 z 3 z ( 1 z ) F y 2 ,
s = 3 Γ y x 3 3 2 x y s 2 ,
where the prime denotes for the derivative to N = l n a and the parameter Γ = V ϕ ϕ V 2 . In this work, we chose V = k ϕ 2 (k is a constant), then Γ = 6 k 1 and s = 2 k 1 2 simultaneously became constants. As a result, the 4-dim system became a 3-dim system made by { x , y , z } .
By choosing F = X + X and Q = α H ρ m , the autonomous system was:
x = 3 2 2 3 x 3 2 x 2 y s α x y ( 1 z ) + α 2 x ,
y = 3 2 x y 2 s + ( 3 2 y 1 z + 3 2 y 3 [ 2 2 x + 1 2 x 2 ] ) ( 1 2 z ) ,
z = 3 z 3 z ( 1 z ) y 2 [ 2 2 x + 1 2 x 2 ] ,
for x > 0 ; while, for x < 0 , the system was
x = 3 2 2 3 x 3 2 x 2 y s α x y ( 1 z ) + α 2 x ,
y = 3 2 x y 2 s + ( 3 2 y 1 z + 3 2 y 3 [ 2 2 x + 1 2 x 2 ] ) ( 1 2 z ) ,
z = 3 z 3 z ( 1 z ) y 2 [ 2 2 x + 1 2 x 2 ] ,
which has two parameters, i.e., the potential one s and the coupling one α . If z 0 , as no loop quantum gravity modification, the dynamical system goes back to the classical k-essence case, as in [69]. Meanwhile, the cosmological quantities (7)–(10) can be rewritten by the phase variables as:
Ω ϕ = 1 2 x 2 y 2 ,
w ϕ = 1 2 | x | 1 ,
w t o t = ( 2 2 | x | + 1 2 x 2 ) y 2 ( 1 z ) ,
q = 3 2 [ ( 1 z ) 1 + y 2 ( 2 2 | x | + 1 2 x 2 ) ] ( 1 2 z ) 1 .
The corresponding critical points { x c r i t , y c r i t , z c r i t } for x > 0 are
P 1 = { 3 6 s , 6 3 ( s 6 ) , 0 } ,
P 2 = { 3 6 + s , 6 3 ( s + 6 ) , 0 } ,
P 3 = { 2 ( α 3 ) 2 2 s 2 α + 9 + α 2 6 α , 6 6 2 s 2 α + α 2 6 α + 9 s ( α 3 ) , 0 } .
While, for x < 0 , the critical points are:
P 4 = { 3 6 + s , 6 3 ( s + 6 ) , 0 } ,
P 5 = { 3 s 6 , 6 3 ( 6 s ) , 0 } ,
P 6 = { 2 ( α 3 ) 2 2 s 2 α + 9 + α 2 6 α , 6 6 2 s 2 α + α 2 6 α + 9 s ( α 3 ) , 0 } .
The six critical points above are shown in Table 1, in which we also present the necessary conditions for their existences, as well as the corresponding cosmological quantities Ω ϕ , w ϕ and q in form of parameters s and α at each critical point. With these cosmological quantities, we could investigate the final state of the universe and discuss whether an acceleration phase exists. The existences are based on the physical meaning of x, y and z, and x > 0 for P 1 , P 2 and P 3 ; however, x < 0 for P 4 , P 5 and P 6 . They also require y > 0 and z > 0 according to the definition. Thus, P 1 and P 4 are excluded by the existence requirement. Meanwhile, 0 w ϕ 1 , and w ϕ < 1 / 3 are needed to accelerate expansion. From the viewpoint of LQC, all critical points have z c r i t = 0 , while x c r i t and y c r i t are as same as in EC [67,69], which implies that the loop quantum gravity effect vanishes and the EC dominates the whole universe in the late-time universe.
The stability of each point, except P 1 and P 4 , is analysed by the three eigenvalues λ 1 , λ 2 , λ 3 of the 3 × 3 Jacobian matrix as in [29,37,38] for the 3-dim system:
Point P 2 : 6 s , 3 6 2 s , 6 s 3 + α ,
Point P 3 : α 3 , 1 4 s ( α 3 ) 2 ( 2 α 3 s + 6 α s 3 + 9 α 2 s 27 s + Δ ) , 1 4 s ( α 3 ) 2 ( 2 α 3 s + 6 α s 3 + 9 α 2 s 27 s Δ ) ,
Point P 5 : 6 s , 3 + 6 2 s , + 6 s 3 + α ,
Point P 6 : α 3 , 1 4 s ( α 3 ) 2 ( 2 α 3 s + 6 α s 3 + 9 α 2 s 27 s + Δ ) , 1 4 s ( α 3 ) 2 ( 2 α 3 s + 6 α s 3 + 9 α 2 s 27 s Δ ) , where Δ = 36 α 2 s 6 2 α 7 + 36 α 5 s 2 108 α 3 s 4 + 42 α 6 459 α 4 s 2 + 648 α 2 s 4 378 α 5 + 2268 α 3 s 2 972 α s 4 + 1890 α 4 5346 α 2 s 2 5670 α 3 + 5832 α s 2 + 10,206 α 2 2187 s 2 10,206 α + 4374 .
Points P 2 and P 5 each have three real eigenvalues, which implies λ i < 0 for stability. However, points P 3 and P 6 , have a real eigenvalue λ 1 and a pair of quadratic eigenvalues λ 2 and λ 3 , which indicates the decomposition of an 1-dim subspace and a 2-dim subspace. The stability requires that the single real eigenvalue λ 1 < 0 ; meanwhile, the 2-dim subspace requires a determinant larger than zero and trace of less than zero, i.e., d e t = λ 2 · λ 3 > 0 and t r = λ 2 + λ 3 < 0 , without determining whether λ 2 and λ 3 are real numbers or complex ones. Additionally, constrained by the existence, stability, Ω ϕ and the accelerating expansion requirement w ϕ < 1 / 3 , the corresponding value ranges for parameters α and s are depicted in Figure 1a for P 2 and Figure 2a for P 3 . Meanwhile, with the observational constrains of α = 0 . 009 0.012 + 0.013 in [59], | α | < 0.01 in [60], and 1.06 < w ϕ < 0.996 in [5], together with the expressions of w ϕ for each stable point in Table 1, parameters α and s are constrained in certain values, i.e., 0.005 < s < 0.073 for P 2 and 0.498 ( α 3 ) 2 / α < s 2 < 0.53 ( α 3 ) 2 / α for P 3 . To make comparisons with different coupling parameter values α , we set | α | 0.02 in Model I. Since the properties for x > 0 are symmetrical to x < 0 , we only show the evolutional trajectories in the 3-dim phase space around P 2 ( s = 0.005 , α = 0.01 ) in Figure 1b and P 3 ( s = 16.5 , α = 0.015 ) in Figure 2b. With suitable parameter values, all trajectories around P 2 and P 3 drop down from z = 0.1 to z c r i t = 0 , and converge to the attractors P 2 = ( 0.709 , 1.996 , 0 ) and P 3 = ( 0.738 , 0.142 , 0 ) on x-y plane, respectively.
Then, we will make two kinds of comparisons: the first one is between LQC and EC from the perspective of trajectories in phase space, and the second one is the effects on the evolution of cosmological quantities for different values of s and α . First, compared with the former work in EC [69], there is no difference in the form of x c and y c for each critical point in LQC, and z c r i t = 0 for all critical points, which indicates that the final states of the universe are scarcely effected in the LQC. However, the early universe is effected in LQC model, which is shown by the projection of evolutional trajectories from the 3-dim phase space onto x-y phase plane for the past, i.e., l n a < 0 , a < 1 , as in Figure 3, where the trajectories of LQC further separate from the EC ones due to the effect of the LQC phase variable z as l n a . As a result, the evolution of each cosmological quantity of LQC varies from that of EC in the early universe, but reaches an identical and stable value in the future, as shown in Figure 4 and Figure 5. (It is interesting that the oscillation of w ϕ in Figure 5 is due to the long period spiral around P 3 , as in Figure 2b; however, w ϕ finally reach an identical and stable value in both LQC and EC for a longer period.) Second, Figure 4 shows that the coupling parameter α and the potential parameter s have little effect on the evolution of the cosmological quantities. From Figure 5, the effects of α and s on the cosmological quantities Ω ϕ and q are not obvious, but the EoS of DE w ϕ is sensitive to both the coupling parameter α and the potential parameter s. It is worth noting that, without the observational constrains on α and s, the evolution curves of each cosmological quantity obviously vary with larger values of α and s.
It is well known that the LQC could lead to a bouncing solution [35,36], which helps to avoid singularities. We drew the evolutions of H, ρ t o t and H ˙ in Figure 6 and Figure 7 for P 2 and P 3 , respectively, which are based on another autonomous dynamical system composed of the derivative of ϕ , ϕ ˙ , H and ρ with respect to t, as in [29,30]:
ϕ ˙ = ϕ ˙ ,
ϕ ¨ = ϕ ˙ 2 ϕ 1 + 3 2 H 6 H ϕ ˙ α H ρ ϕ 2 k ϕ ˙ + 1 2 α H ϕ ˙ ,
H ˙ = 1 2 [ ρ + ( 2 2 ϕ ˙ + ϕ ˙ 2 2 ) k ϕ 2 ] ( 1 2 ρ ρ l c ) ,
ρ ˙ = 3 H [ ρ + ( 2 2 ϕ ˙ + ϕ ˙ 2 2 ) k ϕ 2 ] .
We set the initial value of ϕ 0 , ϕ ˙ 0 , H 0 and ρ 0 under the conditions of the initial loop quantum gravity effect, as z 0 = ρ 0 ρ l c = 0.01 , i.e., ρ 0 = 0.015 . Then, H 0 = 0.07 from Equation (1). For P 2 , k = 16 as s = 0.5 were obtained by the definition; while, ϕ 0 = 20 and ϕ ˙ 0 = 0.75 were obtained for ( x 0 , y 0 ) = ( 0.75 , 1.6 ) around the stable point P 2 . For another stable point P 3 , k = 1 / 30 as s = 11 by the definition, while ϕ 0 = 9 and ϕ ˙ 0 = 0.9 for ( x 0 , y 0 ) = ( 0.9 , 0.17 ) around the stable point P 3 . In both cases, the bouncing is terminal, which indicates that the LQC effect acts not only in the early universe but also in future; then, the Big Rip singularity will not appear in future. In addition, the bounce appears periodically forever and the Hubble rate H with its derivative H ˙ are bounded; then, the curvature singularity is avoided for R = 6 H ˙ + 12 H 2 . From Figure 6b and Figure 7b, when ρ reaches the maximum ρ l c = 1.5 , the Hubble rate becomes zero for Equation (1), and H ˙ also reaches its maximum. On the left-hand side of the bouncing, H < 0 is the universe expanding backwards, while on the right-hand side, the universe keeps expanding for H > 0 . It can be concluded that, in future, although each of the cosmological quantities reach a stable value at the stable point, the total energy density ρ oscillates, which leads to the bounce in H and H ˙ as a recollapse universe.

3. Model II: L 2 = F ( X ) V ( ϕ ) in Loop Quantum Cosmology

In this section, we consider another kind of k-essence dark energy model with Lagrangian
L 2 = p ϕ = F ( X ) V ( ϕ ) ,
then the energy density and the corresponding cosmological quantities, such as the DE density parameter Ω ϕ , the decelerating parameter q, and the EoS parameters, are given by
ρ ϕ = V ( ϕ ) + [ 2 X F X F ] ,
Ω ϕ = ( 2 X F X F + V ) / ( 3 H 2 ) ,
q = 1 + 3 2 ( 2 X F X 3 H 2 + Ω m ) ( 1 2 ρ ρ l c ) ,
w t o t = F V 3 H 2 ( 1 ρ ρ l c ) ,
w ϕ = F V 2 X F X F + V .
According to the method in [50,51], in (40) we let ρ k = 2 X F X F , by setting the phase variables as
x = ρ k 3 H 2 , y = V 3 H 2 , z = ρ ρ l c , σ = 1 3 ρ k V ϕ ϕ ˙ V .
Then, the general autonomous dynamical system of the phase variables with three parameters γ k , Ξ , Γ and interaction Q is as follows,
x = 3 x 2 ( 1 1 z x 2 y 2 + γ k x 2 ) ( 1 2 z ) 3 2 ( γ k x σ y 2 ) Q 6 H 3 ,
y = 3 2 y [ ( 1 1 z x 2 y 2 ) + γ k x 2 ] ( 1 2 z ) 3 2 x y σ ,
z = 3 z + 3 z ( 1 z ) ( x 2 γ k x 2 + y 2 ) ,
σ = 3 σ 2 x ( 1 Γ ) + 3 2 σ γ k 3 2 σ 2 x 1 y 2 Q 6 H 3 3 σ 1 + 2 Ξ ( 1 σ x 1 y 2 γ k 1 ) ,
where γ k = 2 X F X 2 X F X F , Ξ = X F X X F X , Γ = V ϕ ϕ V V ϕ 2 . In [73], this provides a novel setting of phase variables based on a special F ( X ) , as in [57]. Whereas, in this work, by setting F ( X ) = A X η and V = B ϕ n , the three variables γ k = 2 η 2 η 1 , Ξ = η 1 and Γ = n 1 n become constants. The interaction is Q = α H ρ m , as in model I, then Q 6 H 3 = α 2 ( 1 1 z x 2 y 2 ) in (46). For z 0 , the above equations go back to the case in Einstein cosmology, as in [50,51]. By solving the equations above, we obtain eight critical points. For existence, x > 0 , y > 0 and 0 < z < 1 are required by definition. It is also constrained by 0 Ω ϕ 1 for physical meaning in the late-time universe, and w ϕ < 1 / 3 for accelerating expansion. The details are listed in Table 2. In this case, the cosmological quantities expressed by the phase variables are as follows:
Ω ϕ = x 2 + y 2 ,
q = 1 + 3 2 ( 1 1 z x 2 y 2 + γ k x 2 ) ( 1 2 z ) ,
w t o t = ( γ k x 2 x 2 y 2 ) ( 1 z ) ,
w ϕ = x 2 y 2 + γ k x 2 x 2 + y 2 .
The eight critical points for L 2 = A X η B ϕ n are:
P a = { 0 , 0 , 0 , 0 } , f o r α = 0 ,
P b = { 1 , 0 , 0 , 0 } ,
P c = { 0 , 1 , 0 , 0 } ,
P d = { 1 , 0 , 0 , 2 γ k 2 Γ 2 } ,
P e 1 = { α 3 ( 1 γ k ) , 0 , 0 , 9 γ k 3 36 γ k 2 + 45 γ k 18 + Δ e 12 α ( Γ γ k Γ γ k + 1 ) } ,
P e 2 = { α 3 ( 1 γ k ) , 0 , 0 , 9 γ k 3 36 γ k 2 + 45 γ k 18 Δ e 12 α ( Γ γ k Γ γ k + 1 ) } ,
P f 1 = { 3 α , 0 , α 2 + 9 γ k 9 9 ( γ k 1 ) , α γ k + 2 α + Δ f 12 ( Γ 1 ) } ,
P f 2 = { 3 α , 0 , α 2 + 9 γ k 9 9 ( γ k 1 ) , α γ k + 2 α Δ f 12 ( Γ 1 ) } ,
where Δ e = 8 α 2 ( γ k 1 ) ( Γ 1 ) ( α + 3 γ k 3 ) ( α + 3 3 γ k ) + 81 ( γ k 2 ) 2 ( γ k 1 ) 4 and Δ f = 72 γ k ( Γ 1 ) + α 2 ( γ k 2 ) 2 . P a only exists for the coupling parameter α = 0 . It is very interesting that, under the action of the interaction Q, the system has four extra critical points P e 1 , P e 2 , P f 1 and P f 2 , and they are only for α 0 as non-zero interaction between DE and DM. From the perspective of LQC, unlike other critical points, this has z c r i t 0 in P f 1 and P f 2 , which indicates that the loop quantum effect still exists in the final states, corresponding to these two points. However, from the requirement 0 < z < 1 , the existences of P f 1 and P f 2 contradict the physical meaning. That is to say, P f 1 and P f 2 are out of discussion, though P f 1 is mathematically stable, as discussed later.
The stability of the dynamical system around each critical point is analysed by the eigenvalues λ 1 , λ 2 , λ 3 and λ 4 of the 4 × 4 Jacobian matrix, as in [38]:
Point P a : 3 2 , 3 , 3 3 γ k 2 , 6 3 γ k 2 ;
Point P b : α + 3 γ k 3 , 3 γ k / 2 , 3 γ k , 3 γ k / 2 + 3 ;
Point P c : 3 , 0, 3 2 γ k 2 + 2 γ k 2 + 3 γ k 2 6 γ k + 4 2 γ k , 3 2 γ k 2 + 2 γ k 2 3 γ k 2 6 γ k + 4 2 γ k ;
Point P d : α + 3 γ k 3 , 3 ( 2 Γ γ k γ k 2 ) / ( 4 Γ 4 ) , 3 γ k , 3 γ k / 2 3 ;
Point P e 1 : α 2 9 γ k 2 + 18 γ k 9 6 ( γ k 1 ) , α 2 + 9 γ k 9 3 ( γ k 1 ) , Δ e 6 ( γ k 1 ) 2 , ( γ k 1 ) [ 4 α 2 ( Γ 1 ) + 9 ( γ k 1 ) ( 4 Γ + γ k 6 ) ] + Δ e 24 ( γ k 1 ) 2 ( Γ 1 ) ;
Point P e 2 : α 2 9 γ k 2 + 18 γ k 9 6 ( γ k 1 ) , α 2 + 9 γ k 9 3 ( γ k 1 ) , Δ e 6 ( γ k 1 ) 2 , ( γ k 1 ) [ 4 α 2 ( Γ 1 ) + 9 ( γ k 1 ) ( 4 Γ + γ k 6 ) ] Δ e 24 ( γ k 1 ) 2 ( Γ 1 ) ;
Point P f 1 : 3 4 α ( α γ k + Δ f ) , 3 4 α ( α γ k Δ f ) , 3 2 α Δ f , 3 ( α γ k 2 α Δ f ) 8 α ( Γ 1 ) ;
Point P f 2 : 3 4 α ( α γ k + Δ f ) , 3 4 α ( α γ k Δ f ) , 3 2 α Δ f , 3 ( α γ k 2 α + Δ f ) 8 α ( Γ 1 ) .
Among those critical points, only P d and P e 2 satisfy both the existence and stability properties ( x c r i t y c r i t z c r i t σ c r i t is supposed to be equivalently infinitesimal in the calculation of the eigenvalues of P a and P c ). We plotted the value region for parameters α , γ k and Γ in Figure 8a and Figure 9a, together with the dynamical evolution in the 3-dim phase space of { x , y , z } with each initial points and parameter values in Figure 8b and Figure 9b. In Model II, the stability condition and physical properties of some parameters move the value range of the parameter α out of the constrain | α | 0.02 for some stable points. Since P d is not effected by the coupling parameter α , the trajectories of α = 1 and α = 2 start from the same initial point and then converge to the same stable point P d = ( 1 , 0 , 0 , 30 7 ) , but the process differs according to the effect of interaction. Ω ϕ = 1 means the universe is dominated by DE, and ( x d = 1 , y d = 0 ) means the DE is composed by kinetics totally. As P e 2 is effected by the coupling parameter α , the trajectories with α = 0.01 and α = 0.3 start from the same initial point but converge to different stable points. The DE density parameter Ω ϕ = α 2 9 ( γ k 1 ) 2 is affected by the interaction between DE and DM, i.e., the universe is composed of both DE and DM.
By using the symmetry analysis in [50,51] for specific Lagrangian L = A X η B ϕ n , the 3-dim dynamical system is reduced to a 2-dim one in EC. In this work, the analogous analysis acts in the frame of LQC to reduce the 4-dim system to a 3-dim one. With d N = d l n a as an independent variable instesd of d t , the equations of H 2 , H ˙ , ρ ˙ m and ρ ˙ ϕ are rewritten as
H 2 = 1 3 [ ( 2 η 1 ) A X η + B ϕ n + ρ m ] ( 1 ( 2 η 1 ) A X η + B ϕ n + ρ m ρ l c ) ,
H d H d N = 1 2 [ 2 η A X η + ρ m ] ( 1 2 ( 2 η 1 ) A X η + B ϕ n + ρ m ρ l c ) ,
d ρ m d N = 3 ρ m + α ρ m ,
d d N [ ( 2 η 1 ) A X η + B ϕ n ] = 6 η A X η α ρ m .
If we consider ϕ , X, ρ m and ρ l c as the independent variables ( ρ ϕ is composed of ϕ and X), it has the transformation: ϕ ξ 2 η ϕ , X ξ 2 n X , ρ m ξ 2 n η ρ m , ρ l c ξ 2 n η ρ l c . Since X = 1 2 ( H d ϕ d N ) 2 , by the transformation above, the Hubble parameter should hold the transformation H ξ n 2 η H . However, to leave Equations (62) invariant requires H ξ n η H . The relation of n and η can be written as η = n 2 + n , then γ k = 2 η 2 η 1 = 2 n n 2 as a result; further, the symmetry allows for a reduction in the dynamical system from four variables into three, as follows:
σ = n B 1 n 2 3 ( A ( 2 η 1 ) ) 1 2 η ( x y ) 2 n = s ( x y ) 2 n ,
where s = n B 1 n 2 3 ( A ( 2 η 1 ) ) 1 2 η . As a result, the system (46) reduced to the following 3-dim form with three parameters of γ k , s and coupling parameter α as follows:
x = 3 x 2 ( 1 1 z x 2 y 2 + γ k x 2 ) ( 1 2 z ) 3 2 γ k x + 3 2 s x y ( y / x ) 2 γ k α 2 ( 1 1 z x 2 y 2 ) ,
y = 3 2 y [ ( 1 1 z x 2 y 2 ) + γ k x 2 ] ( 1 2 z ) 3 2 s x 2 ( y / x ) 2 γ k ,
z = 3 z + 3 z ( 1 z ) ( x 2 γ k x 2 + y 2 ) .
Then, there are two new critical points P g and P h :
P g = { σ γ k , 1 σ 2 γ k 2 , 0 } ,
P h = { 3 3 σ + α , α 2 + 3 α σ + 9 γ k 9 3 σ + α , 0 }
z c r i t = 0 is found in these two critical points, which indicates that the LQC effect disappears in the late-time universe. However, for both of the critical points containing σ with x / y inside, the expressions cannot be analytical. As a result, the stability analysis by the Jacobian matrix cannot proceed (the interaction Q = α H ρ m diminishes the coupling parameter α by substituting P g into (67)), instead we will take certain parameter values for later calculation. Using the expression (70), P g is not effected by the coupling parameter α , while the interaction effects P h , as in (71). These two critical points are in accordance with [51], in EC for the no interaction case, i.e., α = 0 . In LQC, by setting α = 0 , P h is lucky to have the analytical form P h 0 , which is expressed by parameters γ k and s, by substituting x h = 1 / σ and y h = γ k 1 / σ into (66):
P h 0 = { 1 / σ , γ k 1 / σ , 0 } = { s 1 ( γ k 1 ) γ k 2 2 γ k , s 1 ( γ k 1 ) γ k 1 γ k , 0 }
With the analytical form of P h 0 , the stability can be discussed using the eigenvalues of the 3 × 3 Jacobian matrix of the system, as shown below:
P h 0 : 3 , 3 4 s 1 ( γ k 1 ) 2 / γ k ( s ( γ k 1 ) 2 / γ k + Δ h 0 ) , 3 4 s 1 ( γ k 1 ) 2 / γ k ( s ( γ k 1 ) 2 / γ k Δ h 0 ) ,
in which Δ h 0 = ( γ k 1 ) 2 / γ k [ 8 γ k s 2 ( γ k 1 ) 2 / γ k 9 s 2 ( γ k 1 ) 2 / γ k 8 γ k ( γ k 1 ) 2 ] . We calculated the value ranges of parameters γ k and α under the condition of existence, physical meanings and stabilities for P h 0 , which are listed in Table 3 and depicted in Figure 10a. By setting s = 3 and γ k = 3 2 , from Figure 10b, which shows that P h 0 = ( 2 6 3 , 4 3 6 , 0 ) ( 0.374 , 0.265 , 0 ) has the same stability for both EC and LQC, which indicates that the late-time universe of LQC is in accordance with EC and the loop quantum effects vanish in the late-time universe.
After discussing the critical point P h 0 of the analytical form above, we research the general forms P g and P h , which are non-analytical, below. Here, we set η = 2 , n = 4 and s = 4 , i.e., F ( X ) V ( ϕ ) = A X 2 B ϕ 4 , as in [51], to compare EC and LQC. In this case, γ k = 4 3 and A B = 4 M p l 4 27 ; then, α becomes the only parameter. Numerically, point P g = ( x g , y g , z g ) with the eigenvalues λ i behaves as follows:
x g 0.994 , y g 0.109 , z g = 0 .
λ 1 3.95 , λ 2 0.95 + 0.99 α , λ 3 1.02 .
Compared with [51] in EC, x g and y g are the same as the ones in EC and, in this work, the stability condition requires α < 0.958 to make the critical point stable by λ 2 < 0; while, in [51] α = 0 of no interaction means that this point is a saddle point. Figure 11 shows the effect of coupling parameter α . Although α has no effects on the location of the critical point P g , it effects the stability, i.e., from Figure 11a–c P g becomes a stable point from a saddle point by changing the value of the coupling parameter α . Conversely, for P h , by setting α = 0 , this is the stable point P h 0 = ( 3 4 / 4 , 27 4 / 12 , 0 ) , in accordance with [51] in EC; then, by changing α from 0 to 1.6 , it becomes an unstable point. It can be concluded that coupling parameter α effects the stability of the critical points P g and P h . Further, from Table 2, the critical points P g and P h cannot be stable simultaneously under the stability condition.
Based on all four stable points P d , P e 2 , P g and P h in model II above, to delineate between the effect of LQC and classical EC, the evolution trajectories for both LQC and EC are projected onto the x-y phase plane, as shown in Figure 12a, Figure 13a, Figure 14a and Figure 15a. It can be seen that the dynamical evolutionary trajectories for both LQC and EC, starting from the same initial points around stable points, comply with each other in the late-time universe. However, the trajectories for the early universe as l n a < 0 bifurcate from each other. As a result, the values of cosmological quantities in LQC are varied from EC in the early universe, as the loop quantum gravity effect is obvious; however, for the late-time universe with the diminishing of loop quantum effect, the evolutions of the cosmological quantities show no difference from LQC to EC, as depicted in Figure 12b, Figure 13b, Figure 14b and Figure 15b.
Finally, we discuss the evolution of H ˙ , H and ρ , choosing certain parameter values and initial values in another dynamical system, formed as in Model I, as below:
ϕ ˙ = ϕ ˙ ,
ϕ ¨ = 6 H A η X η B n ϕ n 1 α H ( ρ A ( 2 η 1 ) X η B ϕ n ) ϕ ˙ 1 A η ( 2 η 1 ) X η 1 ,
H ˙ = 1 2 ( ρ + A X η B ϕ n ) ( 1 2 ρ ρ l c ) ,
ρ ˙ = 3 H ( ρ + A X η B ϕ n ) .
With the effect of LQC on cosmological evolution, we can predict that the future universe will experience an oscillating processes, which avoidd the Big Rip singularity and curvature singularity as H and H ˙ are bounded. Figure 16 is based on P e 2 for different initial values of H 0 and ρ 0 , from which the bouncing period is prolonged by larger values of H 0 and ρ 0 . Figure 17a is based on P h 0 for no interaction, i.e., α = 0 , while Figure 17b is based on P g for coupling parameter α = 1.2 , which shows that α has little effect on the bouncing. The initial values of ϕ 0 and ϕ ˙ 0 based on definition (45).

4. Conclusions

In summary, we have investigated the dynamical stability, cosmological evolution, and bouncing universe for the two kinds of k-essence DE models in depth, using the frame of LQC with the interaction Q = α H ρ m in FRWL spacetime. We not only numerically analyzed the dynamical stabilities for Model I with L 1 = F ( X ) V ( ϕ ) and Model II with L 2 = F ( X ) V ( ϕ ) , respectively, and obtained the four stable points among the six critical points for Model I, the five stable points among the eleven critical points for Model II, but also discussed the influences of coupling parameter and potential parameter on the evolutions of several cosmological quantities (such as the density parameter Ω ϕ , EoS of DE w ϕ and deceleration parameter q) around some stable points, especially the comparison between the cosmological evolutions in LQC with those in EC. The research results show that the EoS of DE w ϕ is more sensitive to both coupling parameter and potential parameter than the deceleration parameter q and the density parameter Ω ϕ . Moreover, we found that the evolutions of Ω ϕ , w ϕ and q in the last-time universe in LQC are in accordance with those in EC, although, in the early time, the differences from LQC to EC are great according to the quantum effects. It follows that the results given by us for the k-essence models in LQC could be viewed as a generalization of the results in EC. Finally, with the loop quantum gravity effects, we also achieved the bouncing universe in the two kinds of k-essence models for suitable initial values, which helps to avoid the future singularities and curvature singularity. Obviously, the bounce periods are evidently effected by the initial values, but the interaction does not effect the periods too much.

Author Contributions

Methodology, Y.W.; software, B.C., J.C., W.L. and Y.H.; supervision, Y.W.; writing—original draft, B.C.; writing—review and editing, Y.W. 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, 12175096, 12175095, 11865012, and LiaoNing Revitalization Talents Program grant number XLYC2007047.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare 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. Measurements of Ω and Λ from 42 High-Redshift Supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  3. Spergel, D.N. et al. [WMAP Collaboration]. 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. et al. [SDSS Collaboration]. 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. et al. [Planck Mission 2018 Collaboration]. Planck 2018 results VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [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. Relativ. 2001, 4, 1. [Google Scholar] [CrossRef] [PubMed]
  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. Dark Energy. Commun. Theor. Phys. 2011, 56, 525–604. [Google Scholar] [CrossRef] [Green Version]
  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. Tsujikawa, S. Quintessence: A Review. Class. Quant. Gravit. 2013, 30, 214003. [Google Scholar] [CrossRef] [Green Version]
  16. 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]
  17. Amendola, L.; Tsujikawa, S. Dark Energy: Theory and Observations; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  18. Goswami, G.K.; Pradhan, A.; Beesham, A. A Dark Energy Quintessence Model of the Universe. Mod. Phys. Lett. A 2020, 35, 2050002. [Google Scholar] [CrossRef] [Green Version]
  19. Ahmed, N.; Pradhan, A. Crossing the phantom divide line in universal extra dimensions. New Astron. 2020, 80, 101406. [Google Scholar] [CrossRef]
  20. Varshney, G.; Sharma, U.K.; Pradhan, A.; Kumar, N. Reconstruction of Tachyon, Dirac-Born-Infeld-essence and Phantom model for Tsallis holographic dark energy in f (R, T) gravity. Chin. J. Phys. 2021, 73, 56–73. [Google Scholar] [CrossRef]
  21. Dixit, A.; Maurya, D.C.; Pradhan, A. Phantom dark energy nature of bulk-viscosity universe in modified f (Q)-gravity. Int. J. Geom. Meth. Mod. Phys. 2022, 2250198. [Google Scholar] [CrossRef]
  22. Guo, Z.K.; Cai, R.G. Cosmological Evolution of Interacting Phantom Energy with Dark Matter. J. Cosmol. Astropart. Phys. 2005, 5, 2. [Google Scholar] [CrossRef]
  23. Bojowald, M. Loop Quantum Cosmology. Living Rev. Relativ. 2005, 8, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Li, B.F.; Singh, P.; Wang, A. Towards cosmological dynamics from loop quantum gravity. Phys. Rev. D 2018, 97, 084029. [Google Scholar] [CrossRef] [Green Version]
  25. Li, B.F.; Singh, P.; Wang, A. Qualitative dynamics and inflationary attractors in loop cosmology. Phys. Rev. D 2018, 98, 066016. [Google Scholar] [CrossRef] [Green Version]
  26. Ashtekar, A. Gravity, Geometry and Quantum. AIP Conf. Proc. 2006, 3, 861. [Google Scholar]
  27. Ashtekar, A.; Pawlowski, T.; Singh, P. Quantum Nature of the Big Bang. Phys. Rev. Lett. 2006, 96, 141301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Sami, M.; Singh, P.; Tsujikawa, S. Avoidence of future singularities in loop quantum cosmology. Phys. Rev. D 2007, 74, 043514. [Google Scholar] [CrossRef] [Green Version]
  29. Samart, D.; Gumjudpai, B. Phantom field dynamics in loop quantum cosmology. Phys. Rev. D 2007, 76, 043514. [Google Scholar] [CrossRef] [Green Version]
  30. Cailleteau, T.; Cardoso, A.; Vandersloot, K.; Wands, D. Singularity in loop quantum cosmology. Phys. Rev. Lett. 2008, 101, 251302. [Google Scholar] [CrossRef] [Green Version]
  31. Haro, J. Future singularity avoidence in semiclassical cosmology. Gen. Relativ. Gravit. 2012, 44, 2215–2238. [Google Scholar] [CrossRef]
  32. Kleidis, K.; Oikonomou, V.K. Loop quantum cosmology scalar field models. Int. J. Geom. Meth. Mod. Phys. 2018, 15, 1850071. [Google Scholar] [CrossRef] [Green Version]
  33. Novello, M.; Bergliaffa, S.E.P. Bouncing cosmologies. Phys. Rep. 2008, 463, 127–213. [Google Scholar] [CrossRef]
  34. Brandenberger, R.; Peter, P. Bouncing Cosmologies: Progress and Problems. Found. Phys. 2017, 47, 797–850. [Google Scholar] [CrossRef] [Green Version]
  35. Wu, P.X.; Zhang, S.N. Cosmological evolution of interacting phantom (quintessence) model in loop quantum cosmology. J. Cosmol. Astropart. Phys. 2008, 6, 7. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, S.B.; Wang, B.; Jing, J.L. Dynamics of interacting dark energy model in Einstein and loop quantum cosmology. Phys. Rev. D 2008, 78, 123503. [Google Scholar] [CrossRef]
  37. Fu, X.Y.; Yu, H.W.; Wu, P.X. Dynamics of interacting phantom in loop quantum cosmology. Phys. Rev. D 2008, 78, 063001. [Google Scholar] [CrossRef] [Green Version]
  38. Zonunmawia, H.; Khyllep, W.; Roy, N.; Dutta, J.; Tamanini, N. Extended phase space analysis of interacting DE models in loop quantum cosmology. Phys. Rev. D 2017, 96, 083527. [Google Scholar] [CrossRef] [Green Version]
  39. Oikonomou, V.K. The classical and loop quantum cosmology phase space of interacting dark energy and superfliud dark matter. Phys. Rev. D 2019, 99, 104042. [Google Scholar] [CrossRef] [Green Version]
  40. Oikonomou, V.K. Generalised logarithmic equation of state in classical and loop quantum cosmology dark energy-dark matter coupled system. Ann. Phys. 2019, 409, 167934. [Google Scholar] [CrossRef] [Green Version]
  41. Garriga, J.; Mukhanov, V.F. Perturbations in k-inflation. Phys. Lett. B 1999, 458, 219–225. [Google Scholar] [CrossRef] [Green Version]
  42. Srivastava, S.; Sharma, U.K.; Pradhan, A. New holographic dark energy in Bianchi-III universe with k-essence. New Astron. 2018, 68, 57–64. [Google Scholar] [CrossRef]
  43. Varshney, G.; Sharma, U.K.; Pradhan, A. Reconstructing the k-essence and the dilation field models of the THDE in f (R, T) gravity. Eur. Phys. J. Plus 2020, 135, 541. [Google Scholar] [CrossRef]
  44. Bhardwaj, V.K.; Pradhan, A.; Dixit, A. Compatibility between the scalar field models of tachyon, k-essence and quintessence in f (R, T) gravity. New Astron. 2021, 83, 101478. [Google Scholar] [CrossRef]
  45. Armendariz-Picon, C.; Damour, T.; Mukhanov, V. K-inflation. Phys. Lett. B 1999, 458, 209–218. [Google Scholar] [CrossRef] [Green Version]
  46. 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]
  47. Rendall, A.D. Dynamics of k-essence. Class. Quant. Gravit. 2006, 23, 1557–1570. [Google Scholar] [CrossRef] [Green Version]
  48. 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]
  49. 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]
  50. De-Santiago, J.; Cervantes-Cota, J.L.; Wands, D. Cosmological phase space analysis of the F (X)-V (ϕ) scalar field bouncing solutions. Phys. Rev. D 2013, 87, 023502. [Google Scholar] [CrossRef] [Green Version]
  51. De-Santiago, J.; Cervantes-Cota, J.L. Phase space analysis of the F (X)-V (ϕ) scalar field Lagrangian and scaling solutions in flat cosmology. J. Phys. Conf. Ser. 2014, 485, 012017. [Google Scholar] [CrossRef] [Green Version]
  52. Dutta, J.; Khyllep, W.; Tamanini, N. Cosmological dynamics of scalar fields with kinetic corrections: Beyond the exponential potential. Phys. Rev. D 2016, 93, 063004. [Google Scholar] [CrossRef] [Green Version]
  53. Scherrer, R.J. Purely kinetic k-essence as unified dark matter. Phys. Rev. Lett. 2004, 93, 011301. [Google Scholar] [CrossRef]
  54. Chiba, T.; Okabe, T.; Yamaguchi, M. Kinetically Driven Quintessence. Phys. Rev. D 2000, 62, 023511. [Google Scholar] [CrossRef] [Green Version]
  55. Chimento, L.P.; Feinstein, A. Power-law expansion in k-essence cosmology. Mod. Phys. Lett. A 2004, 19, 761–768. [Google Scholar] [CrossRef]
  56. Chimento, L.P. Extended tachyon field, Chaplygin gas and solvable k-essence cosmologies. Phys. Rev. D 2004, 69, 123517. [Google Scholar] [CrossRef] [Green Version]
  57. Chakraborty, A.; Ghosh, A.; Banerjee, N. Dynamical systems analysis of a k-essence model. Phys. Rev. D 2019, 99, 103513. [Google Scholar] [CrossRef] [Green Version]
  58. 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]
  59. Bolotin, Y.L.; Kostenko, A.; Lemets, O.A. Cosmological evolution with interaction between dark energy and dark matter. Int. J. Mod. Phys. D 2015, 24, 1530007. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, B. Dark Matter and Dark Energy Interactions: Theoretical Challenges, Cosmological Implications and Observational Signatures. Rep. Prog. Phys. 2016, 79, 096901. [Google Scholar] [CrossRef] [Green Version]
  61. Chimento, L.P. Linear and nonlinear interactions in the dark sector. Phys. Rev. D 2010, 81, 043525. [Google Scholar] [CrossRef] [Green Version]
  62. Goswami, G.K.; Pradhan, A.; Beesham, A. FLRW Accelerating Universe with Interactive Dark Energy. Pramana-J. Phys. 2019, 93, 89. [Google Scholar] [CrossRef] [Green Version]
  63. Pradhan, A.; Goswami, G.K.; Beesham, A.; Dixit, A. An FLRW interacting dark energy model of the Universe. New Astron. 2020, 78, 101368. [Google Scholar] [CrossRef]
  64. 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]
  65. 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]
  66. 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]
  67. 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]
  68. Chen, S.B.; Jing, J.L. Dark energy model interacting with dark matter and unparticle. Class. Quant. Gravit. 2009, 26, 155006. [Google Scholar] [CrossRef] [Green Version]
  69. Chen, B.H.; Wu, Y.B.; Xu, D.F.; Dong, W.; Nan, Z. Dynamical stability and geometrical dianostic of the power law k-essence dark energy model with interaction. Universe 2020, 6, 244. [Google Scholar] [CrossRef]
  70. 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]
  71. 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–777, 1–122. [Google Scholar]
  72. Leon, G.; Saridakis, E.N. Phase-space analysis of Horava-Lifshitz cosmology. J. Cosmol. Astropart. Phys. 2009, 11, 6. [Google Scholar] [CrossRef]
  73. Shi, J.L.; Wu, J.P. Dynamics of k-essence in loop quantum cosmology. Chin. Phys. C 2021, 4, 451–461. [Google Scholar] [CrossRef]
Figure 1. (a) The value ranges for parameters s and α to create and stablize the critical point P 2 when x > 0 ; (b) The 3-dim phase space for s = 0.005 , and α = 0.01 around the attractor P 2 = ( 0.709 , 1.996 , 0 ) when x > 0 .
Figure 1. (a) The value ranges for parameters s and α to create and stablize the critical point P 2 when x > 0 ; (b) The 3-dim phase space for s = 0.005 , and α = 0.01 around the attractor P 2 = ( 0.709 , 1.996 , 0 ) when x > 0 .
Universe 08 00520 g001
Figure 2. (a) The value ranges for parameters s and α to create and stabilize the critical point P 3 when x > 0 ; (b) The 3-dim phase space for s = 16.5 , and α = 0.015 around the attractor P 3 = ( 0.738 , 0.142 , 0 ) when x > 0 .
Figure 2. (a) The value ranges for parameters s and α to create and stabilize the critical point P 3 when x > 0 ; (b) The 3-dim phase space for s = 16.5 , and α = 0.015 around the attractor P 3 = ( 0.738 , 0.142 , 0 ) when x > 0 .
Universe 08 00520 g002
Figure 3. The projection of trajectories from 3-dim phase space onto x-y plane, and z-axis is vertical to the x-y plane. (a) The phase plane around P 2 = ( 0.709 , 1.996 , 0 ) . The green curves denote EC and the black curves denote LQC; (b) The phase plane around P 3 = ( 0.738 , 0.142 , 0 ) . The green curve denotes EC and the black curve denotes LQC.
Figure 3. The projection of trajectories from 3-dim phase space onto x-y plane, and z-axis is vertical to the x-y plane. (a) The phase plane around P 2 = ( 0.709 , 1.996 , 0 ) . The green curves denote EC and the black curves denote LQC; (b) The phase plane around P 3 = ( 0.738 , 0.142 , 0 ) . The green curve denotes EC and the black curve denotes LQC.
Universe 08 00520 g003
Figure 4. The comparison of the evolutions of Ω ϕ , w ϕ and q in EC (green) with the ones in LQC (black) around P 2 by the initial condition x 0 = 0.8 , y 0 = 1.5 , z 0 = 0.1 . (ac) correspond to different coupling parameters α = 0 , 0.005 , 0.02 and fixed potential parameter s = 0.004 ; alternately, (df) correspond to different potential parameters s = 0.005 , 0.003 , 0.001 and fixed coupling parameter α = 0.01 .
Figure 4. The comparison of the evolutions of Ω ϕ , w ϕ and q in EC (green) with the ones in LQC (black) around P 2 by the initial condition x 0 = 0.8 , y 0 = 1.5 , z 0 = 0.1 . (ac) correspond to different coupling parameters α = 0 , 0.005 , 0.02 and fixed potential parameter s = 0.004 ; alternately, (df) correspond to different potential parameters s = 0.005 , 0.003 , 0.001 and fixed coupling parameter α = 0.01 .
Universe 08 00520 g004
Figure 5. The comparison of the evolutions of Ω ϕ , w ϕ and q in EC (green) with those in LQC (black) around P 2 by the initial condition x 0 = 0.6 , y 0 = 0.15 , z 0 = 0.1 . (ac) correspond to different coupling parameters α = 0.012 , 0.015 , 0.018 and fixed potential parameter s = 16.5 ; alternately, (df) correspond to different potential parameters s = 13 , 15 , 17 and fixed coupling parameter α = 0.015 .
Figure 5. The comparison of the evolutions of Ω ϕ , w ϕ and q in EC (green) with those in LQC (black) around P 2 by the initial condition x 0 = 0.6 , y 0 = 0.15 , z 0 = 0.1 . (ac) correspond to different coupling parameters α = 0.012 , 0.015 , 0.018 and fixed potential parameter s = 16.5 ; alternately, (df) correspond to different potential parameters s = 13 , 15 , 17 and fixed coupling parameter α = 0.015 .
Universe 08 00520 g005
Figure 6. (a) The evolution of H, H ˙ and ρ in LQC around P 2 , by the initial condition of ϕ 0 = 20 , ϕ ˙ 0 = 0.75 , H 0 = 0.07 and ρ 0 = 0.015 ; (b) The local figure in (a) of t [ 410 , 420 ] where bouncing occurs.
Figure 6. (a) The evolution of H, H ˙ and ρ in LQC around P 2 , by the initial condition of ϕ 0 = 20 , ϕ ˙ 0 = 0.75 , H 0 = 0.07 and ρ 0 = 0.015 ; (b) The local figure in (a) of t [ 410 , 420 ] where bouncing occurs.
Universe 08 00520 g006
Figure 7. (a) The evolutions of H, H ˙ and ρ in LQC around P 3 , by the initial condition of ϕ 0 = 9 , ϕ ˙ 0 = 0.9 , H 0 = 0.07 and ρ 0 = 0.015 ; (b) The local figure in (a) of t [ 255 , 260 ] where bouncing occurs.
Figure 7. (a) The evolutions of H, H ˙ and ρ in LQC around P 3 , by the initial condition of ϕ 0 = 9 , ϕ ˙ 0 = 0.9 , H 0 = 0.07 and ρ 0 = 0.015 ; (b) The local figure in (a) of t [ 255 , 260 ] where bouncing occurs.
Universe 08 00520 g007
Figure 8. (a) The value ranges for parameters α , γ k and Γ to create and stabilize the critical point P d ; (b) The 3-dim phase space for fixed γ k = 8 7 and Γ = 11 10 , corresponding to two different values α = 1 (black) and α = 2 (green) around P d = ( 1 , 0 , 0 , 30 7 ) .
Figure 8. (a) The value ranges for parameters α , γ k and Γ to create and stabilize the critical point P d ; (b) The 3-dim phase space for fixed γ k = 8 7 and Γ = 11 10 , corresponding to two different values α = 1 (black) and α = 2 (green) around P d = ( 1 , 0 , 0 , 30 7 ) .
Universe 08 00520 g008
Figure 9. (a) The value ranges for parameters α , γ k and Γ to create and stabilize the critical point P e 2 ; (b) The 3-dim phase space around P e 2 with different coupling parameter α by fixing γ k = 6 5 and Γ = 6 5 . When α = 0.3 , the black trajectories converge to P e 2 = ( 1 2 , 0 , 0 , 4 ) , while for α = 0.01 , the green trajectories converge to P e 2 = ( 1 60 , 0 , 0 , 120 ) .
Figure 9. (a) The value ranges for parameters α , γ k and Γ to create and stabilize the critical point P e 2 ; (b) The 3-dim phase space around P e 2 with different coupling parameter α by fixing γ k = 6 5 and Γ = 6 5 . When α = 0.3 , the black trajectories converge to P e 2 = ( 1 2 , 0 , 0 , 4 ) , while for α = 0.01 , the green trajectories converge to P e 2 = ( 1 60 , 0 , 0 , 120 ) .
Universe 08 00520 g009
Figure 10. (a) Parameter region for P h 0 ; (b) The phase space around P h 0 = ( 2 6 3 , 4 3 6 , 0 ) ( 0.374 , 0.265 , 0 ) with s = 3 and γ k = 3 2 for EC(green) and LQC(black).
Figure 10. (a) Parameter region for P h 0 ; (b) The phase space around P h 0 = ( 2 6 3 , 4 3 6 , 0 ) ( 0.374 , 0.265 , 0 ) with s = 3 and γ k = 3 2 for EC(green) and LQC(black).
Universe 08 00520 g010
Figure 11. The stabilities corresponding to different values of coupling parameter α , which are shown by the projection from a 3-dim phase space onto to 2-dim x-y plane, and z-axis is vertical to the x-y plane. P g = ( 0.994 , 0.109 , 0 ) is displayed by a solid circle. (a) With α = 0 , trajectories converge to stable point P h = ( 0.329 , 0.190 , 0 ) , just as [51] in EC; (b) With α = 0.6 > 0.985 , trajectories still converge to stable point P h = ( 0.666 , 0.121 , 0 ) ; (c) With α = 1.2 < 0.985 , all trajectories converge to P g which is a stable point in this case, but P h = ( 1.231 , 0.113 , 0 ) becomes unstable.
Figure 11. The stabilities corresponding to different values of coupling parameter α , which are shown by the projection from a 3-dim phase space onto to 2-dim x-y plane, and z-axis is vertical to the x-y plane. P g = ( 0.994 , 0.109 , 0 ) is displayed by a solid circle. (a) With α = 0 , trajectories converge to stable point P h = ( 0.329 , 0.190 , 0 ) , just as [51] in EC; (b) With α = 0.6 > 0.985 , trajectories still converge to stable point P h = ( 0.666 , 0.121 , 0 ) ; (c) With α = 1.2 < 0.985 , all trajectories converge to P g which is a stable point in this case, but P h = ( 1.231 , 0.113 , 0 ) becomes unstable.
Universe 08 00520 g011
Figure 12. (a) The x-y phase plane projected from 4-dim system of {x,y,z, σ } for both LQC(black) and EC(green dash) around P d , corresponding to the parameters of γ k = 8 7 , Γ = 11 10 and α = 1 . The initial points are ( 0.8 , 0.1 , 0.1 , 4 ) for LQC, and ( 0.8 , 0.1 , 0 , 4 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases around P d .
Figure 12. (a) The x-y phase plane projected from 4-dim system of {x,y,z, σ } for both LQC(black) and EC(green dash) around P d , corresponding to the parameters of γ k = 8 7 , Γ = 11 10 and α = 1 . The initial points are ( 0.8 , 0.1 , 0.1 , 4 ) for LQC, and ( 0.8 , 0.1 , 0 , 4 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases around P d .
Universe 08 00520 g012
Figure 13. (a) The x-y phase plane projected from 4-dim system of {x, y, z, σ } for both LQC(black) and EC(green dash) around P e 2 , corresponding to the parameters of γ k = 6 5 , Γ = 6 5 and α = 0.01 . The initial points are ( 0.8 , 0.2 , 0.1 , 4 ) for LQC, and ( 0.8 , 0.2 , 0 , 4 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases around P e 2 .
Figure 13. (a) The x-y phase plane projected from 4-dim system of {x, y, z, σ } for both LQC(black) and EC(green dash) around P e 2 , corresponding to the parameters of γ k = 6 5 , Γ = 6 5 and α = 0.01 . The initial points are ( 0.8 , 0.2 , 0.1 , 4 ) for LQC, and ( 0.8 , 0.2 , 0 , 4 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases around P e 2 .
Universe 08 00520 g013
Figure 14. (a) The x-y phase plane projected from 3-dim system of {x,y,z} for both LQC(black) and EC(green) around P g , corresponding to the parameters of γ k = 4 3 , s = 4 and α = 1.2 . The initial points are ( 0.5 , 0.2 , 0.1 ) for LQC, and ( 0.5 , 0.2 , 0 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases for P g .
Figure 14. (a) The x-y phase plane projected from 3-dim system of {x,y,z} for both LQC(black) and EC(green) around P g , corresponding to the parameters of γ k = 4 3 , s = 4 and α = 1.2 . The initial points are ( 0.5 , 0.2 , 0.1 ) for LQC, and ( 0.5 , 0.2 , 0 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases for P g .
Universe 08 00520 g014
Figure 15. (a) The x-y phase plane projected from 3-dim system of {x,y,z} for both LQC(black) and EC(green) around P h , corresponding to the parameters of γ k = 4 3 , s = 4 and α = 0.6 . The initial points are ( 0.5 , 0.2 , 0.1 ) for LQC, and ( 0.5 , 0.2 , 0 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases for P h .
Figure 15. (a) The x-y phase plane projected from 3-dim system of {x,y,z} for both LQC(black) and EC(green) around P h , corresponding to the parameters of γ k = 4 3 , s = 4 and α = 0.6 . The initial points are ( 0.5 , 0.2 , 0.1 ) for LQC, and ( 0.5 , 0.2 , 0 ) for EC. The star stands for the final state of the evolution, while the solid circle stands for the initial point. (b) the evolution of the cosmological quantities of Ω ϕ , w ϕ and q for both LQC(black) and EC(green) cases for P h .
Universe 08 00520 g015
Figure 16. For P e 2 , the bouncing happens in LQC, with parameters η = 3 , n = 5 , A = 0.002 , B = 0.002 and coupling parameter α = 0.01 . (a) The initial values ϕ 0 = 0.6 , ϕ ˙ 0 = 1.94 , H 0 = 0.21 and ρ 0 = 0.15 . (b) The initial values ϕ 0 = 0.5 , ϕ ˙ 0 = 2.3 , H 0 = 0.3333 and ρ 0 = 0.5 .
Figure 16. For P e 2 , the bouncing happens in LQC, with parameters η = 3 , n = 5 , A = 0.002 , B = 0.002 and coupling parameter α = 0.01 . (a) The initial values ϕ 0 = 0.6 , ϕ ˙ 0 = 1.94 , H 0 = 0.21 and ρ 0 = 0.15 . (b) The initial values ϕ 0 = 0.5 , ϕ ˙ 0 = 2.3 , H 0 = 0.3333 and ρ 0 = 0.5 .
Universe 08 00520 g016
Figure 17. The bouncing happens in LQC, with parameters η = 2 , n = 4 , A = 0.5 , B = 0.3 , with the initial values ϕ 0 = 1.85 , ϕ ˙ 0 = 0.65 , H 0 = 0.212 and ρ 0 = 0.15 . (a) α = 0 for stable point P h 0 ; (b) α = 1.2 for stable point P g .
Figure 17. The bouncing happens in LQC, with parameters η = 2 , n = 4 , A = 0.5 , B = 0.3 , with the initial values ϕ 0 = 1.85 , ϕ ˙ 0 = 0.65 , H 0 = 0.212 and ρ 0 = 0.15 . (a) α = 0 for stable point P h 0 ; (b) α = 1.2 for stable point P g .
Universe 08 00520 g017
Table 1. The existence and stability conditions for six critical points; the cosmological quantities in the form of parameters s and α at each critical point.
Table 1. The existence and stability conditions for six critical points; the cosmological quantities in the form of parameters s and α at each critical point.
  Name ExistenceStability 0 Ω ϕ 1 w ϕ < 1 / 3 q < 0
P 1 none 0 < s < 6 , 6 s + α 3 < 0 1 6 s 3 3 6 s 2 1
P 2 6 s > 0 6 < s < 0 6 s + α 3 < 0 1 6 s + 3 3 6 s 2 1
P 3 ( α 3 ) 2 > 2 s 2 α s ( α 3 ) > 0 λ 1 = α 3 < 0 λ 2 λ 3 > 0 , λ 2 + λ 3 < 0 ( α 3 ) 2 6 s 2 2 α s 2 ( α 3 ) 2 1 α 2
P 4 none 0 < s < 6 , 6 s + α 3 < 0 1 6 s 3 3 6 s 2 1
P 5 6 s > 0 6 < s < 0 6 s + α 3 < 0 1 6 s + 3 3 6 s 2 1
P 6 ( α 3 ) 2 > 2 s 2 α s ( α 3 ) > 0 λ 1 = α 3 < 0 λ 2 λ 3 > 0 , λ 2 + λ 3 < 0 ( α 3 ) 2 6 s 2 2 α s 2 ( α 3 ) 2 1 α 2
Table 2. The existence and stability conditions for eight critical points, the cosmological quantities, and the range of the parameters γ k , Γ and α for acceleration.
Table 2. The existence and stability conditions for eight critical points, the cosmological quantities, and the range of the parameters γ k , Γ and α for acceleration.
   Name ExistenceStability 0 Ω ϕ 1 w ϕ < 1 / 3 q < 0
P a α = 0 unstable0none 3 2
P b alwaysunstable1 γ k 1 3 γ k 2 1
P c alwaysunstable1 1 1
P d always 0 < γ k < 2 , α + 3 γ k < 3 1 < Γ < 1 / γ k + 1 / 2 1 γ k 1 3 γ k 2 1
P e 1 α 1 γ k > 0 unstable α 2 9 ( γ k 1 ) 2 γ k 1 α 2 + 3 γ k 3 6 ( γ k 1 )
P e 2 α 1 γ k > 0 λ 1 < 0 , λ 2 < 0 , λ 3 < 0 , λ 4 < 0 α 2 9 ( γ k 1 ) 2 γ k 1 α 2 + 3 γ k 3 6 ( γ k 1 )
P f 1 nonestable 9 α 2 γ k 1 1
P f 2 noneunstable 9 α 2 γ k 1 1
Table 3. The existence and stability conditions for three critical points under the symmetric analysis, the cosmological quantities in form of the parameters γ k , Γ and α .
Table 3. The existence and stability conditions for three critical points under the symmetric analysis, the cosmological quantities in form of the parameters γ k , Γ and α .
NameExistenceStability Ω ϕ w tot w ϕ q
P g σ / γ k > 0 γ k 2 σ 2 > 0 , α σ + 3 σ 2 3 γ k < 0 1 σ 2 γ k γ k σ 2 γ k γ k 3 σ 2 2 γ k 2 γ k
P h γ k > 1 , s > 0 α σ + 3 σ 2 3 γ k > 0 , α γ k + 3 γ k σ 3 σ > 0 γ k σ 2 001/2
P h 0 α = 0 , γ k > 1 , s > 0 s 2 γ k ( γ k 1 ) ( γ k 2 ) γ k < 1 γ k σ 2 001/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.; Wu, Y.; Chi, J.; Liu, W.; Hu, Y. The Phase Space Analysis of Interacting K-Essence Dark Energy Models in Loop Quantum Cosmology. Universe 2022, 8, 520. https://doi.org/10.3390/universe8100520

AMA Style

Chen B, Wu Y, Chi J, Liu W, Hu Y. The Phase Space Analysis of Interacting K-Essence Dark Energy Models in Loop Quantum Cosmology. Universe. 2022; 8(10):520. https://doi.org/10.3390/universe8100520

Chicago/Turabian Style

Chen, Bohai, Yabo Wu, Jianan Chi, Wenzhong Liu, and Yiliang Hu. 2022. "The Phase Space Analysis of Interacting K-Essence Dark Energy Models in Loop Quantum Cosmology" Universe 8, no. 10: 520. https://doi.org/10.3390/universe8100520

APA Style

Chen, B., Wu, Y., Chi, J., Liu, W., & Hu, Y. (2022). The Phase Space Analysis of Interacting K-Essence Dark Energy Models in Loop Quantum Cosmology. Universe, 8(10), 520. https://doi.org/10.3390/universe8100520

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