Global dynamics of Hoˇrava–Lifshitz cosmology with non-zero curvature and a wide range of potentials

The global dynamics of a cosmological model based on Hoˇrava–Lifshitz gravity in the presence of curvature is described by using the qualitative theory of differential equations.

With or without detailed-balance condition, few researchers have published several papers on Hořava-Lifshitz gravity through phase space analysis, see for instance [3][4][5][6][7]17]. Among these papers there are also some of them about the case of non-zero space curvature k. When the terms related to k are dominant, Friedmann solution will be generated. On the contrary, when the effective Hořava radiation dominates, the corresponding cosmology gives a solution of "radiation-dominated like" [6]. According to whether the cosmological constant is zero in flat (k = 0) and non-flat (k = 0) universe, Leon et al. [3][4][5]14] divided the Hořava-Lifshitz gravity into four cases: (1) k = 0, = 0; (2) k = 0, = 0; (3) k = 0, = 0; (4) k = 0, = 0, and analyzed these phase spaces. More precisely, they either studied the global dynamics of the planar case of Hořava-Lifshitz gravity with exponential potential by using the twodimensional Poincaré compactification, or discussed only local dynamics of this cosmological model with powerlaw potential but without investigate the dynamics close to infinity. The three-dimensional global dynamics of Hořavaa e-mail: gaofabao@sina.com (corresponding author) b e-mail: jllibre@mat.uab.cat Lifshitz cosmology with k = 0 and = 0 was studied in [7].
In this paper we will investigate the global dynamics of the Hořava-Lifshitz scalar field cosmology under the Friedmann-Lemaître-Robertson-Walker background spacetime in the presence of curvature and no cosmological constant term, i.e., k = 0 and = 0.

The cosmological equations
To describe the cosmological equations in this section we first give a brief review of the Hořava-Lifshitz gravitational theory as was proposed in [1]. In one common version of this theory [5,21], its field content can be given by a spatial scalar N , and a spatial vector N i . They are the 'lapse' and 'shift' variables usually found in general relativity. Then the full metric can be written as ds 2 = −N 2 dt 2 + g i j (dx i + N i dt)(dx j + N j dt), where g i j (i, j = 1, 2, 3) is the spatial metric. The scaling transformation of the coordinates is in the form of t → l 3 t, x i → lx i , under which both N and g i j are invari- Based on the detailed-balance condition [1], the full gravitational action of Hořava-Lifshitz is represented by and the extrinsic curvature The standard general covariant antisymmetric tensor is represented by i jm / √ g, and the indices are to raise and lower with the metric g i j . Besides λ is a dimensionless constant, μ, κ and w are also constants. For the potential V (φ), we consider the gravitational action term as follows and the metric N i = 0, is the dimensionless scale factor for the expanding universe, γ i j is a constant curvature metric of maximally symmetric. The curvature k = 1, 0, −1 corresponding to the closed, flat, and open space, respectively.
For simplicity we normalize κ 2 and N to the number one, then the cosmological equations can be described as where H =ȧ/a is the Hubble parameter. Note that the cosmological constant term is neglected in the present paper, i.e., = 0, so system (4) reduces to For simplicity the dimensionless transformation of the following variables [4,5] is given as Since the usual scalar field potential V (φ) admits multiple representations, the potential function f (s) also has various mathematical expressions (see the papers [4,17,22,23] and the references therein). In this paper we will investigate the monomial potential V (φ) = (μφ) 2n /2n with a positive constant μ and a natural number n. Thus the power-law potential is f (s) = −s 2 /(2n), and we obtain ds/dt = √ 6xs 2 /n. Furthermore it can be followed from the first three transformations in Eq. (6) and the first two equations in system (5) that x 2 + y 2 − z 2 = 1 and 3x 2 − 2z 2 = −H /H . Therefore the cosmological equations are reduced to the following autonomous dynamical system ds dt = √ 6 n xs 2 .
Our study will fully describes the global dynamics of system (7) in the region G = (x, z, s) ∈ R 3 : x 2 − z 2 ≤ 1 of physical interest adding its boundary at infinity. It should be noted from the derivation of system (7) that the space curvature k is implicit in the variables x and z. More details of system (7) can be found from Eqs. (113)-(115) of [4] or Eqs.

Phase portraits on the invariant planes and surface
In order to study the local phase portraits of the finite and infinite equilibrium points, and the global phase portraits of system (7) in the region G, which is the meaningful region for cosmology, see again [4] or [14]. We start discussing the phase portraits on its invariant planes and surface 3.1 The invariant plane z = 0 On this plane system (7) becomes The phase portraits of system (8) in the strip z = 0 and x 2 − z 2 ≤ 1, i.e., in z = 0, −1 ≤ x ≤ 1, it has been studied in [7], and the phase portraits is shown in Fig. 1, where the hyperbolic equilibrium point e 0 = (0, 0) is a saddle, both the semi-hyperbolic equilibrium points e 1 = (1, 0) and e 2 = (−1, 0) are saddle-nodes. In order to describe the dark energy or dark matter in the Hořava-Lifshitz universe more clearly, we introduce the density parameters of dark energy and dark matter (see [3,5] for more details): DE = −k 2 z 2 −u 2 , M = x 2 +y 2 , as well as the corresponding equation-of-state parameters: ω DE = (k 2 z 2 − 3u 2 )/(3k 2 z 2 + 3u 2 ), ω M = (x 2 − y 2 )/(x 2 + y 2 ).
x s Fig. 1 The phase portraits of the invariant plane z = 0 in the region −1 ≤ x ≤ 1 which has three equilibrium points e 0 = (0, 0), e 1 = (1, 0) and e 2 = (−1, 0). Here e 0 is a hyperbolic stable node with eigenvalues −3 and −2, both e 1 and e 2 are unstable hyperbolic nodes with eigenvalues 6 and 1. On the local chart U 1 (see Chapter 5 of [24] for more details on the Poincaré compactification) system (9) becomes Since this system vanishes at v = 0, all the points at infinity are equilibrium points. Taking the transformation with respect to time dτ = vdt yields This system has two hyperbolic points at infinity, e 3 = (− √ 6/2, 0) and e 4 = ( √ 6/2, 0), both of them are unstable hyperbolic saddle points with eigenvalues ± √ 6. On the local chart U 2 system (9) writes Rescaling the time of system (12) by letting dτ = vdt we obtain In view of (0, 0) is not an equilibrium point of system (13), we will not continue to study other infinite equilibrium points of system (13) because they have been studied in local chart U 1 . Therefore the global phase portraits of system (9) can be found in Fig. 2. In addition, similar to the previous section, the equilibrium points on this invariant plane and their characteristics are shown in Table 2.
First we prove that the surface x 2 − z 2 = 1 is invariant under the flow of system (7).
for some polynomial K , and this is the case with K = 2(3x 2 − 2z 2 − √ 6xs). On the surface x 2 − z 2 = 1 system (7) can be written as Then except for all the points on x = 0 which are equilibrium points, system (14) also admits two finite equilibrium points e 1 = (1, 0) and e 2 = (−1, 0). By using Theorem 2.19 of [24], we can find that both e 1 and e 2 are semi-hyperbolic saddle-nodes. On the local chart U 1 system (14) becomes It has two infinite equilibrium points e 5 = (0, 0) and e 6 = ( √ 6n/6, 0), where e 5 is a hyperbolic stable node with eigenvalues −1 of multiplicity two, and e 6 is a hyperbolic unstable saddle point with eigenvalues ±1.
On the local chart U 2 system (14) writes Let dτ = udt we obtain The origin e 7 = (0, 0) on the local chart U 2 is a hyperbolic stable node with eigenvalues − √ 6/n of multiplicity two. In short the global phase portraits of system (14) is shown in Fig. 3. Furthermore, similar to the previous Sects. 3.1 and  Table 3 Equilibrium points and their characteristics on the invariant plane 3.2, the equilibrium points and their characteristics on this invariant surface are collated in Table 3.

The finite equilibrium points
It is noted that system (7)

Phase portraits on the Poincaré sphere at infinity
In order to describe the dynamics of system (7) at infinity. We use the method of the three-dimensional Poincaré compactification (see [25] for more details) x = 1/z 3 , z = z 1 /z 3 , s = z 2 /z 3 , then the analytical vector field of system (7) on the In the different local charts of the Poincaré sphere, the infinity of R 3 corresponds to z 3 = 0, then system (18) has the equilibrium points listed in Table 4. Except that the equilibrium point u 31 denotes the origin of the local chart U 3 , the rest equilibrium points lie in the local chart U 1 . Moreover the straight line s = 0 of the local chart U 1 is filled with the equilibrium points u b0 for all b ∈ R.
For the case z 3 = 0 system (18) is reduced to After changing the of time dτ = z 2 dt system (19) becomes Then this system has three equilibrium points e i,1 , e i,2 and e i,3 with coordinates (−1, √ 6n/6), (1, √ 6n/6) and (0, √ 6n/(2(1 + n))), respectively. Here both e i,1 and e i,2 are unstable saddle points with eigenvalues √ 6/n and −2 √ 6. e i,3 is an unstable node with eigenvalues √ 6 and √ 6(1 + z s Fig. 4 The phase portrait of system (7) at infinity on local chart U 1 n)/n. The phase portrait on the Poincaré sphere at infinity on local chart U 1 is shown in Fig. 4. On the local chart U 2 in view of Poincaré compactification x = z 1 /z 3 , z = 1/z 3 , s = z 2 /z 3 , system (7) writes Since we want to study the infinity we take z 3 = 0, rescaling the time dτ = z 2 dt system (21) is equivalent to For any constant c, since (c, 0) is not the equilibrium point of the system (22), we will not continue to discuss other infinite equilibrium points of the system (22). Because this has been discussed in the case U 1 .
On the local chart U 3 the three-dimensional Poincaré compactification is x = z 1 /z 3 , z = z 2 /z 3 , s = 1/z 3 , then system (7) becomes Eur. Phys. J. C (2020) 80:137 For the case z 3 = 0 from system (23) we obtain It is noted that (0, 0) is a linearly zero equilibrium point, the topological index is zero from the Poincaré-Hopf theory (see Theorem 6.30 of [24] for more details). In order to study its local phase portrait we shall use vertical blow-ups (see for instance [26]) by letting w = z 2 /z 1 , then we obtain Eliminating the common factor z 1 of system (25) by changing the time dτ = z 1 dt it yields System (26) (25) and (26), so the direction of the trajectories in the local phase portraits of system (25) is opposite to that of Fig. 5a when z 1 < 0, see Fig. 5b for more details. In addition, all points on the w axis, i.e., z 1 = 0, are singularities of system (25). Thus the local phase portraits at (0, 0) of system (24) is shown in Fig. 5c. Then the phase portrait in the local chart U 3 is shown in Fig. 6.
In summary joining the previous information we obtain the global phase portraits at infinity in the Poincaré sphere in Fig. 7.

Phase portraits inside the Poincaré ball restricted to
It is noted that system (7) is invariant under the symmetries (x, z, s) → (−x, −z, −s) and (x, z, s) → (−x, z, −s), so it is invariant under the symmetry with respect to the origin and to the z-axis. Now we divide the Poincaré ball restricted to x 2 − z 2 ≤ 1 into the following four regions:  (26), (25) and (24), respectively x z Fig. 6 The phase portrait of system (24), i.e., the phase portrait at infinity in the local chart U 3 of system (7) Then due to the symmetries we only need to study the phase portraits of system (7) in the region R 1 .
Combining the phase portraits of the invariant surface x 2 − z 2 = 1 with the phase portraits of the planes z = 0, and s = 0, together with the phase portrait at infinity, we get the phase portraits on the boundary of the region R 1 as shown in Figs. 8, 9 and 10. Here we explain the definition of the three-dimensional cartesian coordinate system in this paper as follows: we regard the xz-plane as the horizontal plane in R 3 , where the direction of the z-axis is horizontal towards the right. If the z-axis axis rotates 90 • clockwise we get the xaxis. The s-axis is vertical upward, and then xzs constitutes a three-dimensional right-handed cartesian coordinate system.
In order to show the phase portraits more clearly, the boundary of the region R 1 is divided into six surfaces according to the orientation towards us, the back to us and the bottom area. See Fig. 11 for more details. It is noted that the North Pole u 31 of the Poincaré ball on the front boundary surfaces F 1 and F 2 is stable, and there is an elliptic sector and a stable

Dynamics in the interior of the region R 1
Note that the original system (7) admits the three finite equilibrium points p 0 , p 1 and p 2 in the three-dimensional cartesian coordinate system. The dynamical behavior of the system inside the region R 1 depends on the behavior of the flow in the following surfaces and planes These surfaces and planes divide the region R 1 into seven different subregions R 1i , i = (1, 2, . . . , 7), see Figs. 12, 13, 14, 15, 16, 17, 18 and 19 for more details. It should be noted that R 11 and R 13 represent parallel tunnels inside R 1 , respectively. It is easy to verify that h > 0 in the subregions R 12 , R 14 , R 15 and R 16 , and h < 0 in the subregions R 11 , R 13 and R 17 . Similarly, we can find that g > 0 in the subregions R 11 , R 12 , R 16 and R 17 , and g < 0 in the subregions R 13 , R 14 and R 15 . It should be noted that the dotted and solid lines in Figs. 13, 14, 15, 16, 17, 18 and 19 are consistent with those in Fig. 12.
As shown in the subregion R 11 (see Fig. 13) the front surface consists of three dashed lines and one solid line contained in the surface h = 0, and the surface on the back side (opposite to the above front surface) contained in the surface g = 0.
In Table 5 we describe the behavior ofẋ,ż andṡ in the seven subregions R 11 , . . . , R 17 . From this table we obtain that the variables x and z in the subregion R 11 decrease monotonically, and the variable s increases monotonically, so an orbit in the subregion R 11 either crosses the boundary h = 0 and enters into the subregion R 12 , or crosses the boundary g = 0 and goes to the subregion R 13 , or crosses the intersection curve of the subregions R 11 , R 12 and R 13 into the subregion R 14 with very low probability, i.e. an orbit in the subregion R 11 will not stay in the future in this region, but will exit through its boundaries into other subregions, and z s x Fig. 18 The subregion R 16 z s x Fig. 19 The subregion R 17 these orbits in backwards time come from the equilibrium point p 1 in the subregion R 11 . This process can be simply summarized as follows The front boundary of R 12 (see Fig. 14), i.e. the surface containing the equilibrium points u 12 , u 14 , u 31 and p 1 , is contained in the invariant surface x 2 − z 2 = 1 (x > 0). The surface on the back side of R 12 that contains the equilibrium points p 1 and u 31 consists of two parts: the upper part surface containing the infinite equilibrium point u 31 is the intersection of the subregions R 12 and R 14 , which is contained in the surface g = 0, and the lower part surface containing the finite equilibrium point p 1 , is the intersection of subregions R 11 and R 12 that is contained in the surface h = 0. The equilibrium points in the subregion R 12 are the finite equilibrium point p 1 , the infinite equilibrium points on the Poincaré sphere at s = 0 and the equilibrium point u 14 at infinity. However an orbit in the subregion R 12 does not start from these infinite equilibrium points at s = 0 or return to these points, but it comes from the finite equilibrium point p 1 or from the subregions R 11 and R 14 which have common boundaries with the subregion R 12 . This is due to the fact that the variable z is decreasing monotonically and the variables x and s are increasing monotonically inside the subregion R 12 according to Table 5. Moreover a trajectory in the subregion R 12 will tend to the equilibrium point u 14 at infinity. In summary, this dynamic behavior process can be represented as In the subregion R 13 (see Fig. 15) the front surface composed of three dashed lines and one solid line is contained in the surface g = 0, and its opposite surface is contained in the surface h = 0. From Table 5 an orbit in the subregion R 13 may only comes from the equilibrium points locate at the infinity on s = 0 in R 13 , or it comes from the subregions R 11 . Then the orbit passes through the intersection of subregions R 13 and R 14 (contained in the surface h = 0) and enters the subregion R 14 , or tends to the infinite equilibrium point u 31 , i.e., lies in the North Pole of the Poincaré sphere. So we obtain the following dynamics In subregion R 14 (see Fig. 16) the front triangular surface is contained in g = 0, and its opposite triangular surface is on the back side, which is included in the plane x = 0. According to Table 5 an orbit of the subregion R 14 may come from the infinite equilibrium points on x = 0, or come from the subregion R 13 through the surface h = 0, or from subregion R 15 traversing plane x = 0, then go through the surface g = 0 into the subregion R 12 , or directly goes to the infinite equilibrium point u 31 at the North Pole of the Poincaré sphere. This can be represented as It can be noted from Fig. 17 that the surface on the left of subregion R 15 (consisting of two dashed lines and one solid line) is contained in plane x = 0, and the opposite surface to it is included in the surface g = 0. According to Table 5 an orbit in the subregion R 15 can only come from the adjacent subregion R 16 , and then either enter into the subregion R 14 adjacent to it, or directly tend to the finite equilibrium point p 0 located at the center of the Poincaré ball. This dynamics can be denoted as We observe that the curved surfaces of the left and right sides of the subregion R 16 (see Fig. 18) are contained in the surfaces g = 0 and h = 0, respectively. It is known from Table 5 that an orbit in the subregion R 16 can only come from the subregion R 17 and crossing the surface g = 0, it enters into the subregion R 15 . Thus we get that The left and right surfaces of subregion R 17 (see Fig. 19) are included in the surface h = 0 and in the invariant surface x 2 − z 2 = 1 (x < 0), respectively. According to Table 5, it is known that all three variables in subregion R 17 are decreasing monotonically, so an orbit in this subregion must come from the infinite equilibrium point u 31 located at the North Pole of Poincaré sphere and eventually tends to the infinite equilibrium point U 13 in the invariant plane s = 0. Then we have that The dynamic behavior of the orbits inside the seven subregions of R 1 discussed above can be summarized as The above flow chart shows that the orbits of system (7) contained in the interior of the region R 1 have α-limit at the finite equilibrium point p 1 , and the some orbits on the boundary of the region R 1 have α-limit at the x = 0 in R 14 , or s = 0 in R 13 . Moreover the orbits have ω-limit either at the finite equilibrium point p 0 , or at the infinite equilibrium points u 14 and U 13 , where u 14 is located on the intersection curve of the Poincaré sphere and the invariant surface x 2 − z 2 = 1 (x > 0) at infinity, and U 13 is located on the intersection of the Poincaré sphere and the invariant surfaces s = 0, x 2 − z 2 = 1 (x < 0) at infinity (see Fig. 12). Furthermore the orbits of system (7) on the boundary of R 1 also have α-limit and ω-limit at the infinite equilibrium point u 31 , i.e. the North Pole of the Poincaré sphere.
For the sake of summary, similar to the previous Sects. 3.1-3.3, the equilibrium points on the region R 1 and their characteristics within this region are integrated in Table 6.
In this way the qualitative global dynamic behavior of system (7) is described.

Discussion and conclusions
From a cosmological point of view, the unstable finite equilibrium points p 1 and p 2 as shown in Figs. 8, 9 and 10 or 12 are consistent with the results of case 1 in the cosmological model in [3,5], corresponding to complete dark matter domination. It is known from the last flow chart in Sect. 6 that the physical importance of the finite equilibrium point p 0 is a late-time state of the universe if the initial condition is not on the invariant plane z = 0 or on the invariant surface x 2 − z 2 = 1. We also note that this late-time solution was already obtained in the strip − √ 6/2 < s < √ 6/2 for the cosmological model with exponential potential in [3][4][5], but our results hold for all values of s in the region G of the three-dimensional space. In addition, the main differences between the Ref. [5] and the present paper are: 1. The former specifically studied the case of exponential potential, our research is based on a wide range of potentials, and specifically considering the case of power-law potential. 2. The former is based on the planar cosmological model, and we are considering a three-dimensional one. 3. The former did not discussed the dynamics at infinity, but we studied the infinite equilibrium points, and found that the infinite equilibrium points u 14 and U 13 may also be the late-time state solutions of the universe.
See the last flow chart in Sect. 6 again, the equilibrium points u 14 and U 13 will be future attractors located at infinity. When the initial conditions are inside the region R 1 , they can be the late-time states of the universe. It is noted that only u 14 and U 13 are studied in the region R 1 because system (7) admits two aforementioned symmetries. Thus the results of u 15 and u 13 (see Fig. 7) are identical with those of u 14 and U 13 , respectively. It should also be noted that Fig. 5.5 in [3] or Fig. 3 in [5] can be obtained in Fig. 7 in the present paper when s = 0.6. Then it is easy to find that u 14 corresponds to P 6 in the Fig. 5.5 in [3] for the plane case when s → 0. However the infinite equilibrium point u 14 can be the latetime state of the universe from the above analysis, while P 6 corresponds to the universe dominated by dark matter due to the absence of and the effect of k will decrease with the increasing of the scale factor a(t).