Next Article in Journal
On the Existence, Uniqueness and a Numerical Approach to the Solution of Fractional Cauchy–Euler Equation
Previous Article in Journal
I-Convergence Sequence Paranormed Spaces of Order (α, β)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System

1
School of Information, Zhejiang Guangsha Vocational and Technical University of Construction, Dongyang 322100, China
2
Department of Big Data Science, School of Science, Zhejiang University of Science and Technology, Hangzhou 310023, China
3
College of Sustainability and Tourism, Ritsumeikan Asia Pacific University, Beppu 874-8577, Oita, Japan
4
School of Electronic and Information Engineering, Taizhou University, Taizhou 318000, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 625; https://doi.org/10.3390/axioms13090625
Submission received: 27 July 2024 / Revised: 8 September 2024 / Accepted: 9 September 2024 / Published: 12 September 2024
(This article belongs to the Section Mathematical Analysis)

Abstract

:
Aiming to explore the subtle connection between the number of nonlinear terms in Lorenz-like systems and hidden attractors, this paper introduces a new simple sub-quadratic four-thirds-degree Lorenz-like system, where x ˙ = a ( y x ) , y ˙ = c x x 3 z , z ˙ = b z + x 3 y , and uncovers the following property of these systems: decreasing the powers of the nonlinear terms in a quadratic Lorenz-like system where x ˙ = a ( y x ) , y ˙ = c x x z , z ˙ = b z + x y , may narrow, or even eliminate the range of the parameter c for hidden attractors, but enlarge it for self-excited attractors. By combining numerical simulation, stability and bifurcation theory, most of the important dynamics of the Lorenz system family are revealed, including self-excited Lorenz-like attractors, Hopf bifurcation and generic pitchfork bifurcation at the origin, singularly degenerate heteroclinic cycles, degenerate pitchfork bifurcation at non-isolated equilibria, invariant algebraic surface, heteroclinic orbits and so on. The obtained results may verify the generalization of the second part of the celebrated Hilbert’s sixteenth problem to some degree, showing that the number and mutual disposition of attractors and repellers may depend on the degree of chaotic multidimensional dynamical systems.

1. Introduction

As the fourteenth mathematical problem of the twenty-first century collected by Smale [1], revealing the nature of the Lorenz attractor has continued to be a hot topic of ongoing research in nonlinear science since its introduction [2,3,4,5,6,7,8]. As part of this ongoing effort, when studying the chaos of three-dimensional quadratic autonomous differential systems using the contraction map and boundary problem, Shilnikov et al. introduced the following classification: chaos of the Shilnikov homoclinic orbit, or heteroclinic orbit, or homoclinic and heteroclinic orbits hybrid orbit type, etc. [7]. Combining numerical technique and theoretical analysis, Kokubu et al. gave some explanation of Lorenz-like attractors from the viewpoint of the collapse of singularly degenerate heteroclinic cycles [9,10]. Llibre and Zhang applied homogeneous-weight polynomials and the method of characteristic curves to solve the linear partial differential equations in order to study invariant algebraic surface of the Lorenz system, the collapse of which may generate strange attractors [11]. Liao et al. argued that the existence of a global attractive compact set and having at least one positive Lyapunov exponent are the two sufficient conditions of a continuous system exhibiting chaos, and verified their findings using the chaotic Lorenz family [12,13]. Based on the algebraic structure and topological characterization, Letellier et al. introduced the concepts of Lorenz-like systems and attractors [8]. Chen separated the vector fields of Lorenz chaotic family into the linear and nonlinear parts [14], while others divided them into the conservative, dissipative and the external force field part [15,16,17].
However, Kuznetsov et al. turned their attention to the relationship between the degree of the considered model and Lorenz-like attractors and generalized the second part of the celebrated Hilbert’s sixteenth problem: the degree may control the number and mutual disposition of attractors and repellers [18,19]. Zhang and Chen reasserted this conjecture and coined two coexisting two-scroll Lorenz attractors from a cubic Lorenz system [20]. Motivated by that, Wang et al. guessed that decreasing the powers of x of the second and third equations of the quadratic Lorenz-like system [21] may widen the range of the parameter c for hidden attractors, and verified this via two sub-quadratic Lorenz-like analogues with degrees of four-thirds and six-fifths [22,23].
Now, one can not help but wonder what happens when decreasing the powers of x of the cross products x z and x y of the Lorenz-like system [21], especially for the self-excited and hidden attractors. To the best of our knowledge, little attention seems to have been paid to this problem. Furthermore, this newly reported Lorenz-like system also satisfies the second criterion of Sprott [24], i.e., the main contribution of this study, validating the generalization of the second part of the Hilbert’s sixteenth problem to some degree: the decrease in the powers of nonlinear terms may narrow or even eliminate the scope of some certain parameters for hidden Lorenz-like attractors, but enlarge it for self-excited attractors. This compelled us to carry out the research detailed here.

2. New Four-Thirds-Degree Lorenz-like System and Its Main Dynamics

By replacing the nonlinear term c x 3 in the sub-quadratic Lorenz-like system [22] with the linear one c x , we formulate the analogue as follows:
x ˙ = a ( y x ) , y ˙ = c x x 3 z , a 0 , ( b , c ) R 2 , z ˙ = b z + x 3 y .
In order to distinguish system (1) from the systems in [21,22,23], we must first present its basic dynamics in the following propositions. We have done this indentation.
Proposition 1.
(i) 
If b = 0 (resp. b 0 and b c 0 ), then E z = { ( 0 , 0 , z ) | z R } is the non-isolated equilibria (resp. a single equilibrium point) of system (1).
(ii) 
If b c > 0 , then E ± = ( ± ( b c ) 3 2 , ± ( b c ) 3 2 , b c 2 ) is a pair of nontrivial equilibrium points in system (1) beside E 0 .
Remark 1.
As in the systems in [21,22,23], a generic (resp. degenerate) pitchfork bifurcation at E 0 (resp. E z ) occurs in system (1) when b 0 (resp. c 0 ) and c (resp. b) passes through the zero value and b c > 0 .
Proposition 2.
For a 0 and ( b , c ) R 2 (resp. b = 0 and z 0 ), Table 1 (resp. Table 2) lists the local dynamics of E 0 (resp. E z ).
Remark 2.
As for the system in [22], using a linear analysis, one can easily obtain the characteristic equations of E 0 and E z : ( λ + b ) ( λ 2 + a λ a c ) = 0 and λ ( λ 2 + a λ a ( c z 3 x 2 3 ) ) , where x 0 , and from which Proposition 2 follows.
In the next proposition, let us discuss the local dynamics of E ± .
Proposition 3.
Make S = { ( a , b , c ) | a 0 , b c > 0 } , S 1 = { ( a , b , c ) S : a + b > 0 , a b + b c 2 a c 3 > 0 , 2 a b c 3 > 0 } , S 2 = S S 1 and
S 1 1 = { ( a , b , c ) S 1 : a b ( a + b ) c [ a ( 2 a + b ) 3 b 2 ] < 0 } , S 1 2 = { ( a , b , c ) S 1 : c = 3 a b ( a + b ) ( a b ) ( 2 a + 3 b ) } , S 1 3 = { ( a , b , c ) S 1 : a b ( a + b ) c [ a ( 2 a + b ) 3 b 2 ] > 0 } .
Then, E ± is unstable (resp. asymptotically stable) when ( a , b , c ) S 1 1 (resp. S 1 3 ). However, when ( a , b , c ) S 1 2 , system (1) undergoes Hopf bifurcation at E ± .
As stated in [21] (Proposition 2.4, p. 2567) (resp. Proposition 3), the non-trivial equilibria E ± of the following quadratic Lorenz-like system is
x ˙ = a ( y x ) , y ˙ = c x x z , a 0 , ( b , c ) R 2 , z ˙ = b z + x y ,
(resp. system (1)) is asymptotically stable when 0 < c < a ( a + b ) a b (resp. 0 < c < 3 a b ( a + b ) ( a b ) ( 2 a + 3 b ) ). Due to 3 a b ( a + b ) ( a b ) ( 2 a + 3 b ) < a ( a + b ) a b , system (1) may experience chaotic behaviors coexisting with the unstable origin and stable E ± in a narrower range of the parameter c in contrast to the quadratic one (2).
Likewise, for ( a , b ) = ( 4 , 1 ) , the E ± of system (1) (resp. (2)) is asymptotically stable when 0 < c < 20 11 (resp. 0 < c < 20 3 ), and Figure 1 shows the periodic behavior rather than chaotic attractors displayed in system (2) [22] (Figure 3, p. 363).
For ( a , b ) = ( 3 , 1.5 ) and c [ 0.1 , 599.1 ] , the quadratic Lorenz-like system (2) mainly experiences periodic behaviors, whereas system (1) mainly experiences chaotic ones, as shown in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6.
Therefore, compared with another two sub-quadratic Lorenz-like analogues [22] (Figures 1–2, p. 362), [23] (Property, Figures 2–4, p. 2450071-5-7) and Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, one may obtain the convincing argument:
Property. A decrease in the powers of nonlinear terms of the quadratic Lorenz-like system (2) may narrow or even eliminate the range of the parameter c for hidden attractors, but enlarge it for self-excited attractors.
Meanwhile, unlike most of other Lorenz-like systems [9,10,22,23,25,26,27,28], the collapse of the singularly degenerate heteroclinic cycles in the system in (1) makes it hard to create strange attractors, as shown in the following numerical result.
Numerical Result . 2 . 1 According to the dynamics of E z in Table 2, for a > 0 , z 1 < 0 and t , the one-dimensional unstable manifold W u ( E z 1 ) ( E z 1 = ( 0 , 0 , z 1 ) ) tending towards the stable E z 2 = ( 0 , 0 , z 2 ) with z 2 > 0 creates singularly degenerate heteroclinic cycles, as shown in Figure 3a. Moreover, a tiny perturbation in b > 0 may change singularly degenerate heteroclinic cycles to limit cycles, as depicted in Figure 3b.
Remark 3.
When ( a , b ) = ( 3 , 1.5 ) , E ± is asymptotically stable (resp. unstable) when 0 < c < 3.8571 (resp. c > 3.8571 ). However, when ( a , c , b ) = ( 3 , 3.8571 , 1.5 ) , the system in (1) undergoes Hopf bifurcation at E ± .
Finally, similarly to [21,22,23,26,28,29,30,31,32,33,34,35], we will discuss the heteroclinic orbits of the system in (1) and present it in the following proposition.
Proposition 4.
If c > 0 and 3 b 4 a > 0 , then (a) the ω-limit of any one trajectory of system (1) is an equilibrium point; (b) system (1) has a pair of heteroclinic orbits to E 0 and E ± but no homoclinic orbits.
Remark 4.
Generally speaking, the equilibria of heteroclinic orbits are all saddles or saddle-foci, or unstable nodes [7]. As heteroclinic orbits [21,22,23,26,28,29,30,31,32,33,34,35], those discussed in Proposition 4 are heteroclinic wiggles [7] (Figure 14.3.2, p. 439), which connect the stable E ± and unstable E 0 .
The rest of the paper is arranged as follows. Section 3 studies the stability and Hopf bifurcation of E ± followed by the proof of Proposition 3. Section 4 discusses the heteroclinic orbits and the proof of Proposition 4 is outlined. A conclusion is drawn and the subject of future work is discussed in Section 5, particularly related to the the relationship between the degree and the Lorenz-like attractors.

3. Hopf Bifurcation

Using the theory of Hopf bifurcation, we then sketched a proof for Proposition 3.
Proof of Proposition 3.
First of all, the characteristic equation of E ± is calculated:
λ 3 + ( a + b ) λ 2 + [ a b + b c 2 a c 3 ] λ + 2 a b c 3 = 0 .
Next, based on the Routh–Hurwitz criterion and Equation (3), we derived the stability of E ± ; however, we have omitted its proof here.
( a , b , c ) S 1 2 , λ 1 , 2 = ± ω i = ± a b [ 1 + ( 3 b 2 a ) ( a + b ) ( a b ) ( 2 a + 3 b ) ] i and λ 3 = ( a + b ) are a pair of conjugate purely imaginary roots and one negative real root for Equation (3), respectively. Moreover, one has
d Re ( λ 1 ) d c | c = c * = a b 3 b 2 + 2 a 2 6 [ ω 2 + ( a + b ) 2 ] = ( 3 b + 2 a ) ( a b ) 6 [ ω 2 + ( a + b ) 2 ] 0 ,
from which the transversal condition is verified. Therefore, Hopf bifurcation happens at E ± . □
Next, we applied the project method [36,37] to compute the Lyapunov coefficients, aiming to determine the nondegeneracy (or stability) of the Hopf bifurcation at E ± .
Firstly, based on the time and coordinate transformation
( x , y , z , t ) ( x 3 , y , z , 1 3 x 2 3 t ) ,
system (1) can be transformed to the equivalent one:
x ˙ = a ( y x 3 ) , y ˙ = 3 x 2 ( c x 3 x z ) , z ˙ = 3 x 2 ( b z + x y ) .
The E ± of system (1) corresponds to E 1 , 2 = ( ± b c , ± ( b c ) 3 , b c 2 ) in system (4). One can verify the transversality of Hopf bifurcation at E 1 , 2 .
In fact, the characteristic equation at E 1 , 2 is
λ 3 + 3 b c ( a + b ) λ 2 + 3 ( b c ) 2 [ ( 3 b 2 a ) c + 3 a b ] λ + 18 a ( b c ) 4 = 0 ,
with λ 1 , 2 = ± ω i = ± 3 ( b c * ) 2 [ ( 3 b 2 a ) c * + 3 a b ] i and λ 3 = 3 b c * ( a + b ) < 0 when ( a , b , c ) S 1 2 . We then obtained the following derivative
d Re ( λ 1 ) d c | c = c * = Δ 2 [ ω 2 + 9 ( b c * ) 2 ( a + b ) 2 ] 0 ,
where Δ = 3 b ( a + b ) ω 2 + 72 a b 4 c * 3 3 a b ( a + b ) [ 9 b 2 c * 2 ( 3 b 2 a ) + 18 a b 3 c * ] , which thus verifies the transversality of Hopf bifurcation of E 1 , 2 .
Then, the following transformation
( x , y , z ) ( x + b c * , y + ( b c * ) 3 , z + b ( c * ) 2 ) ,
converts system (4) into the resulting one
x ˙ y ˙ z ˙ = 3 a b c * a 0 6 b 2 c * 3 0 3 ( b c * ) 3 3 ( b c * ) 3 3 ( b c * ) 2 3 b ( b c * ) 3 x y z + 3 a b c * x 2 7 c * ( b c * ) 3 x 2 3 b c * x z 3 ( b c * ) 5 x 2 + 4 ( b c * ) 3 x y 3 b 2 c * x z + a x 3 9 b c * ( 3 c * x 3 x 2 z ) 9 ( b c * b x 2 z + ( b c * ) 2 x 3 + 2 b c * x 2 y ) + 0 3 x 3 ( 5 c * b c * x z ) 3 x 3 ( 4 b c * y b z + ( b c * ) 3 x ) + 0 3 c * x 5 3 x 4 y .
To compute the first Lyapunov coefficient l 1 , one has to distill the following multi-linear symmetric functions from system (4)
B ( x , y ) = 6 a b c * x 1 y 1 42 c * ( b c * ) 3 x 1 y 1 9 b c * ( x 1 y 3 + x 3 y 1 ) 18 ( b c * ) 5 x 1 y 1 + 12 ( b c * ) 3 ( x 1 y 2 + x 2 y 1 ) 9 b 2 c * ( x 1 y 3 + x 3 y 1 ) ,
C ( x , y , z ) = 6 a x 1 y 1 z 1 162 b c * c * x 1 y 1 z 1 18 b c * ( x 3 y 1 z 1 + x 1 y 3 z 1 + x 1 y 1 z 3 ) 18 b b c * ( x 3 y 1 z 1 + x 1 y 3 z 1 + x 1 y 1 z 3 ) + 54 ( b c * ) 2 x 1 y 1 z 1 + 36 b c * ( x 2 y 1 z 1 + x 1 y 2 z 1 + x 1 y 1 z 2 )
When l 1 = 0 , one needs to compute the other multi-linear symmetric functions to get the second Lyapunov exponent or the third or even higher order ones.
Due to complex algebraic structure of system (6) itself, it is difficult to obtain the explicit form of l 1 . However, one can easily calculate it for a concrete problem, e.g., ( a , c , b ) = ( 4 , 20 11 , 1 ) . Here, E 1 , 2 = ( ± 20 11 , ± ( 20 11 ) 3 , 400 121 ) , whose eigenvalues are λ 1 , 2 = ± ω i = ± 5.3713 i and λ 3 = 27.2727 , and the transversality condition holds: d Re ( λ 1 ) d c | c = c * = 20 11 0.5251 > 0 . Moreover, the l 1 of E 1 , 2 is discussed in the following proposition.
Proposition 5.
For ( a , c , b ) = ( 4 , 20 11 , 1 ) , system (6) undergoes a Hopf bifurcation at E 1 , 2 , for which the first Lyapunov coefficient is l 1 91.2608 < 0 , and thus E 1 , 2 are both weakly unstable foci. Because of d Re ( λ 1 ) d c 0.5251 > 0 , the Hopf bifurcation at E 1 , 2 is supercritical. In a word, for c > c * = 20 11 when it is close to c * = 20 11 , there is at least a pair of stable close orbits around the unstable E 1 , 2 .
Proof. 
Based on the method proposed in [36,37], one can easily obtain the following expressions
p = 0.0318 + 0.0349 i 0.0173 + 0.0233 i 0.00019 0.02346 i ,   q = 4 21.8181 + 5.3712 i 23.5357 15.9337 i ,   h 11 = 108.4504 462.1017 435.16077 ,
h 20 = 55.35369 + 20.7412 i 375.6721 + 261.7939 i 750.3901 236.67872 i ,   G 21 = 182.5216 1713.4399 i and l 1 = 1 2 G 21 = 91.2608 . Since d Re ( λ 1 ) d c 0.5251 > 0 , the Hopf bifurcation at E 1 , 2 is supercritical. Namely, when c = 1.8182 > c * , there exists a pair of stable close orbits around the unstable E 1 , 2 = ( ± 1.3484 , ± 2.4517 , 3.3059 ) for system (4), i.e., E ± = ( ± 2.4517 , ± 2.4517 , 3.3059 ) of system (1), as illustrated in Figure 7. This finishes the proof. □
In the next section, we will study the existence of heteroclinic orbits in system (1). For the sake of our argument, we list the following symbols:
(1) ψ ( t ; ψ 0 ) = ( x ( t ; x 0 ) , y ( t ; y 0 ) , z ( t ; z 0 ) ) : a solution for system (1) with the initial condition ψ 0 = ( x 0 , y 0 , z 0 ) .
(2) γ = { ψ ( t ; ψ 0 ) | ψ ( t ; ψ 0 ) = ( x + ( t ; x 0 ) , y + ( t ; y 0 ) , z + ( t ; z 0 ) ) W u ( E 0 ) , t R } (resp. γ + = { ψ + ( t ; ψ 0 ) | ψ + ( t ; ψ 0 ) = ( x + ( t ; x 0 ) , y + ( t ; y 0 ) , z + ( t ; z 0 ) ) W + u ( E 0 ) , t R } ): the negative (resp. positive) branch of W u ( E 0 ) with x + < 0 (resp. x + > 0 ) when t .

4. Existence of Heteroclinic Orbit

In this section, as in [21,22,23,26,28,29,30,31,32,33,34,35], with a suitable choice of Lyapunov functions, the proof of Proposition 4 can be divided into two stages: (1) 3 b 4 a > 0 , (2) 3 b 4 a = 0 .

4.1. 3 b 4 a > 0

This subsection introduces the first Lyapunov function
V 1 ( ψ ( t ; ψ 0 ) ) = 1 2 [ b ( b 4 a 3 ) ( y x ) 2 + ( b z + x 4 3 ) 2 + 3 b 4 a 6 a ( b c x 2 3 + x 4 3 ) 2 + 3 b 4 a 12 a ( b 2 c 2 + x 4 3 ) 2 ]
and the following assertions:
Lemma 1.
If c > 0 and 3 b 4 a > 0 , then we have the following results:
1. 
Assume t 1 , 2 , t 1 < t 2 and V 1 ( ψ ( t 1 ; ψ 0 ) ) = V 1 ( ψ ( t 2 ; ψ 0 ) ) . ψ 0 is one of the stationary points.
2. 
If lim t ψ ( t ; ψ 0 ) = E 0 and x ( t 3 ; x 0 ) < 0 , t 3 R , then we arrive at V 1 ( E 0 ) > V 1 ( ψ ( t ; ψ 0 ) ) and x ( t ; x 0 ) < 0 , t R . Namely, ψ 0 γ .
Proof. 
(1) By taking the derivative of V 1 with respect to ψ ( t ; ψ 0 ) , we arrive at
d V 1 ( ψ ( t ; ψ 0 ) ) d t | ( 1 ) = a b ( b 4 a 3 ) ( y x ) 2 b ( b z + x 4 3 ) 2 ,
and derive
y ( t ; y 0 ) x ( t ; x 0 ) , b z ( t ; z 0 ) x 4 3 ( t ; x 0 ) ,
under the condition of (1), t ( t 1 , t 2 ) .
Based on system (1) and Equation (8), the identities x ˙ ( t ; x 0 ) y ˙ ( t ; y 0 ) z ˙ ( t ; z 0 ) 0 hold, t ( t 1 , t 2 ) . Namely, ψ 0 is a fixed point.
(2) Now, t R , we prove the fact V 1 ( E 0 ) > V 1 ( ψ ( t ; ψ 0 ) ) . If not, t R , V 1 ( E 0 ) V 1 ( ψ ( t ; ψ 0 ) ) . In fact, the first assertion suggests that ψ 0 is just an equilibrium point, contradicting the assumption that lim t ψ ( t ; ψ 0 ) = E 0 and x ( t 3 ; x 0 ) < 0 . In a word, V 1 ( E 0 ) > V 1 ( ψ ( t ; ψ 0 ) ) , t R .
Next, let us prove x ( t ; x 0 ) < 0 , t R . Otherwise, x ( t 4 ; x 0 ) 0 , t 4 R . Since x ( t 3 ; x 0 ) < 0 , t 3 R , one arrives at x ( t 5 ; x 0 ) = 0 , t 5 R . Due to V 1 ( E 0 ) > V 1 ( ψ ( t ; ψ 0 ) ) , t R , one gets ψ ( t 5 ; ψ 0 ) ( x , y , z ) | V 1 ( x , y , z ) < V 1 ( E 0 ) ( x , y , z ) | x = 0 . Further, the following statement is derived, { ( x , y , z ) | V 1 ( x , y , z ) < V 1 ( E 0 ) } { ( x , y , z ) | x = 0 } = { ( 0 , y , z ) | 1 2 [ b ( b 4 a 3 ) y 2 + b 2 z 2 + ( 3 b 4 a ) b 4 c 4 12 a ] < ( 3 b 4 a ) b 4 c 4 24 a } = , which is impossible. As a result, one arrives at x ( t ; x 0 ) < 0 , for all t R . This ends the proof. □
Lemma 2.
When c > 0 and 3 b > 4 a > 0 , all of the solutions for system (1) tend towards an equilibrium point such as t . Namely, system (1) has no closed orbits.
Proof. 
On the basis of Equation (7), one derives lim t + V 1 ( ψ ( t ; ψ 0 ) ) = Φ ( ψ 0 ) , 0 V 1 ( ψ ( t ; ψ 0 ) ) V 1 ( ψ ( 0 ; ψ 0 ) ) = V 1 ( ψ 0 ) , t 0 , and obtains that x ( t ; x 0 ) , y ( t ; y 0 ) and z ( t ; z 0 ) are all bounded, t [ 0 , + ) , i.e., ψ ( t ; ψ 0 ) is bounded.
Denoting the ω -limit set of ψ ( t ; ψ 0 ) by Ω ( ψ 0 ) . For w Ω ( ψ 0 ) , { t n } , we have
lim n + t n = + , lim n + ψ ( t n , ψ 0 ) = w .
Next, for all t R , ψ ( t ; w ) = lim n + ψ ( t ; ψ ( t n ; ψ 0 ) ) = lim n + ψ ( t + t n ; ψ 0 ) suggests V 1 ( ψ ( t ; w ) ) = V 1 [ lim n + ψ ( t ; ψ ( t n ; ψ 0 ) ) ] = lim n + V 1 ( ψ ( t + t n ; ψ 0 ) ) = Φ ( ψ 0 ) . As a result, w { E + , E , E 0 } . Because of connectedness of Ω ( ψ 0 ) , one only deduces Ω ( ψ 0 ) = { E + } , Ω ( ψ 0 ) = { E } or Ω ( ψ 0 ) = { E 0 } , yielding that ψ ( t ; ψ 0 ) approaches an equilibrium point such as t + . The proof is completed. □
Lastly, one considers the existence of heteroclinic orbits with the help of above two lemmas.
Theorem 1.
When c > 0 and 3 b > 4 a > 0 , the following two statements hold.
1. 
System (1) has no homoclinic orbits.
2. 
A pair of heteroclinic orbits at E ± and E 0 exists in system (1).
Proof. 
For c > 0 and 3 b > 4 a > 0 , one firstly shows that homoclinic orbits at E ± or E 0 do not exist in system (1). If not, a homoclinic orbit at E 0 , E + or E can be denoted by ψ ( t ) , i.e., lim t ± ψ ( t ) = e ± , where e = e + { E , E + , E 0 } .
From Equation (7), one obtains
V 1 ( e ) V 1 ( ψ ( t ) ) V 1 ( e + ) .
In both cases, we obtain the fact that V 1 ( e ) = V 1 ( e + ) and further arrive at V 1 ( ψ ( t ) ) V 1 ( e + ) . From the first statement in Lemma 1, ψ ( t ) is only an equilibrium point. Namely, homoclinic orbits at E 0 , E + or E are non-existent.
Then, let us prove that γ is a heteroclinic orbit at E 0 and E , i.e., lim t + p ( t ) = E . From the concept of γ and Lemma 1, we only arrive at x + ( t ) < 0 , for all t R , yielding lim t + ψ ( t ) E 0 . As a result, lim t + ψ ( t ) = E is true.
Finally, let us prove the uniqueness of γ .
Let ψ 1 ( t ) = ( x 1 ( t ) , y 1 ( t ) , z 1 ( t ) ) be a solution for system (1) such that lim t ± ψ 1 ( t ) = e 1 ± , and { e 1 , e 1 + } = { E 0 , E } . Like Equation (9), one gets V 1 ( e 1 ) V 1 ( ψ 1 ( t ) ) V 1 ( e 1 + ) , t R , from Equation (8). Due to V 1 ( E 0 ) > V 1 ( E ) , one derives e 1 = E 0 and e 1 + = E , i.e.,
lim t ψ 1 ( t ) = E 0 , lim t + ψ 1 ( t ) = E ,
which leads to ψ 1 ( t ) γ from the second assertion of Lemma 1. Since system (1) is axis-symmetrical with respect to the z-axis, a single heteroclinic orbit γ + , i.e., the ones at E 0 and E + , also exists in system (1). This completes the proof. □

4.2. 3 b 4 a = 0

Firstly, we introduce another Lyapunov function
V 2 ( ψ ( t ; ψ 0 ) ) = 1 2 [ ( y x ) 2 + 3 8 a 2 ( 4 a c 3 x 2 3 + x 4 3 ) 2 + 3 16 a 3 ( 16 a 2 c 2 9 + x 4 3 ) 2 ]
from which the following lemma is deduced:
Lemma 3.
When c > 0 and 3 b = 4 a > 0 , one arrives at the following four assertions.
(i) If lim t ψ ( t ; ψ 0 ) is bounded, then Q ( ψ ( t ; ψ 0 ) ) = z ( t ; z 0 ) 3 4 a x 4 3 ( t ; x 0 ) = 0 .
(ii) If 4 a z ( t ; z 0 ) = 3 x 4 3 ( t ; x 0 ) , then d V 2 ( ψ ( t ; ψ 0 ) ) d t | ( 1 ) = a ( y x ) 2 0 .
(iii) If 4 a z ( t ; z 0 ) = 3 x 4 3 ( t ; x 0 ) and V 2 ( ψ ( t 1 ; ψ 0 ) ) = V 2 ( ψ ( t 2 ; ψ 0 ) ) , t 1 , 2 , t 1 < t 2 , then ψ 0 is an equilibrium point.
(iv) If lim t ψ ( t ; ψ 0 ) = E 0 and x ( t 3 ; x 0 ) < 0 , t 3 R , then V 2 ( E 0 ) > V 2 ( ψ ( t ; ψ 0 ) ) and x ( t ; x 0 ) < 0 , t R . Namely, ψ 0 γ .
Proof. 
(i) When c > 0 and 3 b = 4 a > 0 , the derivative of Q ( ψ ( t ; ψ 0 ) ) = z ( t ; z 0 ) 3 4 a x 4 3 ( t ; x 0 ) is calculated as follows: d Q ( ψ ( t ; ψ 0 ) ) d t | ( 1 ) = 4 a 3 Q ( ψ ( t ; ψ 0 ) ) , i.e.,
Q ( ψ ( t ; ψ 0 ) ) = Q ( ψ ( τ ; ψ 0 ) ) e 4 a 3 ( t τ ) , τ , t R .
Since lim τ ψ ( τ ; ψ 0 ) is bounded, Equation (10) yields Q ( ψ ( t ; ψ 0 ) ) 0 , i.e., z ( t ; z 0 ) 3 4 a x 4 3 ( t ; x 0 ) .
(ii) The fact that 4 a z ( t ; z 0 ) 3 x 4 3 ( t ; x 0 ) and system (1) result in Conclusion (ii) of Lemma 3.
(iii) Based on assumed conditions and the above statement, we obtain d V 2 ( ψ ( t , ψ 0 ) ) d t | ( 1 ) = 0 , for all t ( t 1 , t 2 ) , i.e.
y ( t ; y 0 ) x ( t ; x 0 ) .
In virtue of x ˙ , Equation (11) and 4 a z ( t ; z 0 ) 3 x 4 3 ( t ; x 0 ) , one deduces
x ˙ ( t ; x 0 ) y ˙ ( t ; y 0 ) z ˙ ( t ; z 0 ) 0 , t ( t 1 , t 2 ) .
Hence, ψ 0 is only an equilibrium point.
(iv) First of all, one shows that V 2 ( E 0 ) > V 2 ( ψ ( t ; ψ 0 ) ) , t R . If not, V 2 ( E 0 ) V 2 ( ψ ( t 0 ; ψ 0 ) ) , t 0 R . In addition, the first, second and third assertions yield that ψ 0 is an equilibrium point, which contradicts lim t ψ ( t ; ψ 0 ) = 0 and x ( t 3 ; x 0 ) < 0 . Namely, we find that V 2 ( E 0 ) > V 2 ( ψ ( t ; ψ 0 ) ) holds for t R .
Then, we prove x ( t ; x 0 ) < 0 , t R . If not, t 4 R , such that x ( t 4 ; x 0 ) 0 . Due to x ( t 3 ; x 0 ) < 0 , t 5 R , such that x ( t 5 ; x 0 ) = 0 . Since V 2 ( E 0 ) > V 2 ( ψ ( t ; ψ 0 ) ) , t R , we obtain ψ ( t 5 ; ψ 0 ) { ( x , y , z ) | V 2 ( x , y , z ) < V 2 ( E 0 ) } { ( x , y , z ) | x = 0 } . On the other hand, { ( x , y , z ) | V 2 ( x , y , z ) < V 2 ( E 0 ) } { ( x , y , z ) | x = 0 } = { ( 0 , y , z ) | 1 2 [ y 2 + 16 a c 4 27 ] < 8 a c 4 27 } = . As such, a contradiction happens. Therefore, x ( t ; x 0 ) < 0 holds, t R . □
Lemma 4.
Set c > 0 and 3 b = 4 a > 0 . If ψ ( t , ψ 0 ) is bounded when t , then lim t ψ ( t , ψ 0 ) E 0 , or E ± . In a word, there are no closed orbits in system (1).
Proof. 
From the first and second assertion of Lemma 3, lim t V 2 ( ψ ( t ; ψ 0 ) ) = Ψ ( ψ 0 ) exists. Assume h α ( ψ 0 ) , i.e., { t n } , we have lim t n ψ ( t n ; ψ 0 ) = h , n + . For all t R ,
ψ ( t ; h ) = lim n + ψ ( t ; ψ ( t n ; ψ 0 ) ) = lim n + ψ ( t + t n ; ψ 0 )
leads to
ψ ( t ; h ) is bounded , t R , V 2 ( ψ ( t ; h ) ) = lim n + V 2 ( ψ ( t + t n ; ψ 0 ) ) = Ψ ( ψ 0 ) .
On the basis of Lemma 3, one obtains h { E , E 0 , E + } . Therefore,
α ( ψ 0 ) { E , E 0 , E + } .
Because α ( ψ 0 ) is connected, one derives α ( q 0 ) = { E } , or α ( q 0 ) = { E 0 } , or α ( q 0 ) = { E + } , suggesting that lim n + ψ ( t ; ψ 0 ) approaches an equilibrium point. The proof is completed. □
Theorem 2.
When 3 b = 4 a > 0 and c > 0 , we derive the statements as follows.
(i) There are no homoclinic orbits in system (1).
(ii) A pair of heteroclinic orbits E ± and E 0 exist in system (1).
Proof. 
(i) When 3 b = 4 a > 0 and c > 0 , we are able to prove the non-existence of homoclinic orbits connecting E + , E and E 0 in system (1). Otherwise, we assume that ψ ( t ) = ( x ( t ) , y ( t ) , z ( t ) ) is a homoclinic orbit to E + , E or E 0 , i.e.,
lim t ± ψ ( t ) = e ± , e = e + { E , E 0 , E + } .
Based on Lemma 3 and V 2 ( e ) = V 2 ( e + ) , we find that ψ ( t ) is only a stationary point.
As such, homoclinic orbits to E ± or E 0 are non-existent in system (1).
(ii) Then, we prove the uniqueness of γ , i.e., the heteroclinic orbit at E 0 and E .
Suppose ψ 1 ( t ) = ( x 1 ( t ) , y 1 ( t ) , z 1 ( t ) ) is a solution of system (1) such that
lim t ± ψ 1 ( t ) = e 1 ± , { e 1 , e 1 + } = { E 0 , E } .
For all t R , the first and second assertions of Lemma 3 yield
V 2 ( e 1 ) V 2 ( ψ 1 ( t ) ) V 2 ( e 1 + ) .
Due to V 2 ( E 0 ) > V 2 ( E ) , one derives e 1 = E 0 and e 1 + = E , i.e.,
lim t + ψ 1 ( t ) = E and lim t ψ 1 ( t ) = E 0 ,
which leads to ψ 1 ( t ) γ based on Lemma 3.
Finally, one shows that γ is just the heteroclinic orbit to E and E 0 ; that is, lim t + ψ ( t ) = E . According to Lemma 3, we deduce the following:
z + ( t ; z 0 ) 3 4 a x + 4 3 ( t ; x 0 ) , d V 2 ( ψ ( t ) ) d t | ( 1 ) = a ( y + ( t ; y 0 ) x + ( t ; x 0 ) ) 2 , V 2 ( ψ ( t ) ) < V 2 ( E 0 ) , x + ( t ) < 0 , t R .
Based on the second part of Equation (14), one arrives at lim t V 2 ( ψ ( t ) ) = v 2 . Again, Equation (14) suggests that x + ( t ) , y + ( t ) and z + ( t ) are all bounded, and also shows the boundedness of ψ ( t ) , for all t [ 0 , + ) . We denote the ω -limit set of ψ ( t ) as Ω ; that is, for all h Ω , { t n } , we obtain lim n + t n = + and lim n + ψ ( t n ) = h . In a word, for all t R ,
ψ ( t ; h ) = lim n + ψ ( t ; ψ ( t n ) ) = lim n + ψ ( t + t n ) ,
ψ ( t ; h ) is bounded , for all t R , V 2 ( ψ ( t ; h ) ) = lim n + V 2 ( ψ ( t + t n ) ) = v 2 ,
and Lemma 3 all lead to h { E , E 0 } . Therefore, Ω { E , E 0 } . Because of the connectedness of Ω , the relation Ω = E or Ω = E 0 holds. Based on Lemma 3 and Equation (14), one derives Ω E 0 . In a word, Ω = E , i.e., lim n + ψ + ( t ) = E . Therefore, there exists a single heteroclinic orbit to E and E 0 . Due to the symmetry, there exists another unique heteroclinic orbit to E + and E 0 in system (1), as shown in Figure 8. The proof is finished. □

5. Conclusions

Combining a theoretical analysis and numerical simulation, this paper investigates a newly reported 3D sub-quadratic four-thirds-degree Lorenz-like system, and reveals most of inherent dynamics of the Lorenz system family, i.e., self-excited attractors, Hopf bifurcation, generic and degenerate pitchfork, heteroclinic orbits, singularly degenerate heteroclinic cycle, an invariant algebraic surface, etc.
In contrast to the existing quadratic and sub-quadratic Lorenz-like analogues, we may find a new property: decreasing the powers of nonlinear terms may narrow and even eliminate the range of some certain parameters for hidden attractors, but enlarge it for self-excited attractors. This may verify the generalization of the second part of the celebrated Hilbert’s sixteenth problem to some degree: the number and mutual disposition of attractors and repellers may depend on the degree of polynomials of chaotic multidimensional dynamical systems. However, previous studies mainly emphasized that the aforementioned dynamics may explain the forming mechanism of strange attractors.
In future, we will expect other researchers to test this property through more Lorenz-like analogues, and clarify the relationship between other complex dynamics and the degrees, shedding light on the nature of chaos and providing reference for chaos-based applications.

Author Contributions

Conceptualization, G.K.; Methodology, J.P.; Software, G.K., J.P. and H.W.; Validation, H.W.; Investigation, G.K. and J.P.; Writing—original draft, G.K. and J.P.; Writing—review & editing, H.W.; Visualization, G.K., J.P., F.H. and H.W.; Supervision, J.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part Natural Science Foundation of Zhejiang Guangsha Vocational and Technical University of construction under Grant 2022KYQD-KGY, in part National Natural Science Foundation of China under Grant 12001489, in part Zhejiang Public Welfare Technology Application Research Project of China Grant LGN21F020003, in part Natural Science Foundation of Taizhou University under Grant T20210906033.

Data Availability Statement

There are no data because the results obtained in this paper can be reproduced based on the information given in this paper.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Smale, S. Mathematical problems for the next century. Math. Intell. 1998, 20, 7–15. [Google Scholar] [CrossRef]
  2. Lorenz, E.N. Deterministic nonperiodic flow. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  3. Afraimovich, V.S.; Bykov, V.V.; Shilnikov, L.P. The origin and structure of Lorenz attractor. Sov. Phys. Dokl. 1977, 22, 253–255. [Google Scholar]
  4. Tucker, W. The Lorenz attractor exists. Comptes Rendus l’Académie Sci. Ser. I Math. 1999, 328, 1197–1202. [Google Scholar] [CrossRef]
  5. Viana, M. What’s new on Lorenz strange attractors? Math. Intell. 2000, 22, 6–19. [Google Scholar] [CrossRef]
  6. Stewart, I. Mathematics: The Lorenz attractor exists. Nature 2000, 406, 948–949. [Google Scholar] [CrossRef] [PubMed]
  7. Shilnikov, L.P.; Shilnikov, A.L.; Turaev, D.V.; Chua, L.O. Methods of Qualitative Theory in Nonlinear Dynamics Part I, II; World Scientific: Singapore, 2001. [Google Scholar]
  8. Letellier, C.; Mendes, E.M.A.M.; Malasoma, J. Lorenz-like systems and Lorenz-like attractors: Definition, examples, and equivalences. Phys. Rev. E 2023, 108, 044209. [Google Scholar] [CrossRef]
  9. Kokubu, H.; Roussarie, R. Existence of a singularly degenerate heteroclinic cycle in the Lorenz system and its dynamical consequences: Part I. J. Dyn. Differ. Equ. 2004, 16, 513–557. [Google Scholar] [CrossRef]
  10. Messias, M. Dynamics at infinity and the existence of singularly degenerate heteroclinic cycles in the Lorenz system. J. Phys. A Math. Theor. 2009, 42, 115101. [Google Scholar] [CrossRef]
  11. Llibre, J.; Zhang, X. Invariant algebraic surfaces of the Lorenz system. J. Math. Phys. 2002, 43, 1622–1645. [Google Scholar] [CrossRef]
  12. Liao, X.; Yu, P.; Xie, S.; Fu, Y. Study on the global property of the smooth Chua’s system. Int. J. Bifurc. Chaos 2006, 16, 2815–2841. [Google Scholar] [CrossRef]
  13. Liao, X. New Research on Some Mathematical Problems of Lorenz Chaotic Family; Huazhong University of Science & Technology Press: Wuhan, China, 2017. [Google Scholar]
  14. Chen, G. Generalized Lorenz systems family. arXiv 2020, arXiv:2006.04066. [Google Scholar]
  15. Pasini, A.; Pelino, V. A unified view of Kolmogorov and Lorenz systems. Phys. Lett. A 2000, 275, 435–446. [Google Scholar] [CrossRef]
  16. Pelino, V.; Maimone, F.; Pasini, A. Energy cycle for the lorenz attractor. Chaos Solitons Fractals 2014, 64, 67–77. [Google Scholar] [CrossRef]
  17. Liang, X.; Qi, Q. Mechanical analysis of Chen chaotic system. Chaos Solitons Fractals 2017, 98, 173–177. [Google Scholar] [CrossRef]
  18. Leonov, G.A.; Kuznetsov, N.V. On differences and similarities in the analysis of Lorenz, Chen, and Lu systems. Appl. Math. Comput. 2015, 256, 334–343. [Google Scholar] [CrossRef]
  19. Kuznetsov, N.V.; Mokaev, T.N.; Kuznetsova, O.A.; Kudryashova, E.V. The Lorenz system: Hidden boundary of practical stability and the Lyapunov dimension. Nonlinear Dyn. 2020, 102, 713–732. [Google Scholar] [CrossRef]
  20. Zhang, X.; Chen, G. Constructing an autonomous system with infinitely many chaotic attractors. Chaos 2017, 27, 071101. [Google Scholar] [CrossRef]
  21. Liu, Y.; Yang, Q. Dynamics of a new Lorenz-like chaotic system. Nonl. Anal. RWA 2010, 11, 2563–2572. [Google Scholar] [CrossRef]
  22. Wang, H.; Ke, G.; Pan, J.; Hu, F.; Fan, H. Multitudinous potential hidden Lorenz-like attractors coined. Eur. Phys. J. Spec. Top. 2022, 231, 359–368. [Google Scholar] [CrossRef]
  23. Wang, H.; Pan, J.; Ke, G. Revealing more hidden attractors from a new sub-quadratic Lorenz-like system of degree 6 5 . Int. J. Bifurc. Chaos 2024, 34, 2450071. [Google Scholar] [CrossRef]
  24. Sprott, J.C. A proposed standard for the publication of new chaotic systems. Int. J. Bifurc. Chaos 2011, 21, 2391–2394. [Google Scholar] [CrossRef]
  25. Llibre, J.; Messias, M.; Silva, P.R. On the global dynamics of the Rabinovich system. J. Phys. A Math. Theor. 2008, 41, 275210. [Google Scholar] [CrossRef]
  26. Wang, H.; Ke, G.; Pan, J.; Su, Q. Conjoined Lorenz-like attractors coined. Miskolc Math. Notes 2023. [Google Scholar]
  27. Ke, G. Creation of three-scroll hidden conservative Lorenz-like chaotic flows. Adv. Theory Simulations 2024. [Google Scholar] [CrossRef]
  28. Wang, H.; Ke, G.; Hu, F.; Pan, J.; Dong, G.; Chen, G. Pseudo and true singularly degenerate heteroclinic cycles of a new 3D cubic Lorenz-like system. Results Phys. 2024, 56, 107243. [Google Scholar] [CrossRef]
  29. Wang, H.; Pan, J.; Ke, G. Multitudinous potential homoclinic and heteroclinic orbits seized. Electron. Res. Arch. 2024, 32, 1003–1016. [Google Scholar] [CrossRef]
  30. Wang, H.; Pan, J.; Ke, G.; Hu, F. A pair of centro-symmetric heteroclinic orbits coined. Adv. Cont. Discr. Mod. 2024, 2024, 14. [Google Scholar] [CrossRef]
  31. Chen, Y.; Yang, Q. Dynamics of a hyperchaotic Lorenz-type system. Nonlinear Dyn. 2014, 77, 569–581. [Google Scholar] [CrossRef]
  32. Li, T.; Chen, G.; Chen, G. On homoclinic and heteroclinic orbits of the Chen’s system. Int. J. Bifurc. Chaos 2006, 16, 3035–3041. [Google Scholar] [CrossRef]
  33. Tigan, G.; Constantinescu, D. Heteroclinic orbits in the T and the Lü system. Chaos Solitons Fractals 2009, 42, 20–23. [Google Scholar] [CrossRef]
  34. Liu, Y.; Pang, W. Dynamics of the general Lorenz family. Nonlinear Dyn. 2012, 67, 1595–1611. [Google Scholar] [CrossRef]
  35. Tigan, G.; Llibre, J. Heteroclinic, homoclinic and closed orbits in the Chen system. Int. J. Bifurc. Chaos 2016, 26, 1650072. [Google Scholar] [CrossRef]
  36. Kuzenetsov, Y.A. Elements of Applied Bifurcation Theory, 3rd ed.; Springer: New York, NY, USA, 2004; Volume 112. [Google Scholar]
  37. Sotomayor, J.; Mello, L.F.; Braga, D.C. Lyapunov coefficients for degenerate Hopf bifurcations. arXiv 2007, arXiv:0709.3949. [Google Scholar]
Figure 1. For ( a , b ) = ( 4 , 1 ) , c [ 0 , 5 ] and ( x 0 1 , y 0 1 , z 0 1 ) = ( 0.13 , 1.3 , 1.6 ) × 10 7 ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (1). In contrast to system (2) [22] (Figure 3, p. 363), these figures suggest that the solutions for the system in (1) display stable equilibria and period orbits, rather than the self-excited and hidden attractors shown in the system in (2).
Figure 1. For ( a , b ) = ( 4 , 1 ) , c [ 0 , 5 ] and ( x 0 1 , y 0 1 , z 0 1 ) = ( 0.13 , 1.3 , 1.6 ) × 10 7 ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (1). In contrast to system (2) [22] (Figure 3, p. 363), these figures suggest that the solutions for the system in (1) display stable equilibria and period orbits, rather than the self-excited and hidden attractors shown in the system in (2).
Axioms 13 00625 g001
Figure 2. For ( a , b ) = ( 3 , 1.5 ) , c [ 0.1 , 599.1 ] and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (2). The subfigures (ac) are consistent with the subfigure (d), showing that system (2) mainly experiences periodic behaviors.
Figure 2. For ( a , b ) = ( 3 , 1.5 ) , c [ 0.1 , 599.1 ] and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (2). The subfigures (ac) are consistent with the subfigure (d), showing that system (2) mainly experiences periodic behaviors.
Axioms 13 00625 g002
Figure 3. For ( a , b ) = ( 3 , 1.5 ) , c [ 0.1 , 599.1 ] and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (1). In contrast with Figure 2, the four subfigures show that system (1) mainly behaves in a similar way to self-excited attractors, verifying the introduced property, i.e., a decrease in powers of nonlinear terms of the quadratic Lorenz-like system (2) may narrow or even eliminate the range of the parameter c for hidden attractors, but enlarge it for self-excited attractors.
Figure 3. For ( a , b ) = ( 3 , 1.5 ) , c [ 0.1 , 599.1 ] and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) ; (ac) bifurcation diagrams; (d) Lyapunov exponents versus c of system (1). In contrast with Figure 2, the four subfigures show that system (1) mainly behaves in a similar way to self-excited attractors, verifying the introduced property, i.e., a decrease in powers of nonlinear terms of the quadratic Lorenz-like system (2) may narrow or even eliminate the range of the parameter c for hidden attractors, but enlarge it for self-excited attractors.
Axioms 13 00625 g003
Figure 4. Phase portraits of system (1) for ( a , c , b ) = ( 3 , 100 , 1.5 ) and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) illustrating the existence of two-scroll self-excited attractor suggested in Figure 3.
Figure 4. Phase portraits of system (1) for ( a , c , b ) = ( 3 , 100 , 1.5 ) and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) illustrating the existence of two-scroll self-excited attractor suggested in Figure 3.
Axioms 13 00625 g004aAxioms 13 00625 g004b
Figure 5. Poincaré cross-sections of system (1) for ( a , c , b ) = ( 3 , 100 , 1.5 ) and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) showing the geometrical structure of the Lorenz-like attractor depicted in Figure 4.
Figure 5. Poincaré cross-sections of system (1) for ( a , c , b ) = ( 3 , 100 , 1.5 ) and ( x 0 2 , y 0 2 , z 0 2 ) = ( 1.314 , 2.236 , 4.669 ) showing the geometrical structure of the Lorenz-like attractor depicted in Figure 4.
Axioms 13 00625 g005
Figure 6. Phase portraits of system (1) for ( a , c ) = ( 1 , 36 ) , (a) b = 0 , (b) b = 0.07 , and ( x 0 1 , 3 , y 0 1 , 3 , z 0 3 ) = ( ± 1.3 × 10 8 , ± 1.3 × 10 7 , 1 ) . Both figures imply that collapsing singularly degenerate heteroclinic cycles in system (1) create limited cycles rather than strange attractors.
Figure 6. Phase portraits of system (1) for ( a , c ) = ( 1 , 36 ) , (a) b = 0 , (b) b = 0.07 , and ( x 0 1 , 3 , y 0 1 , 3 , z 0 3 ) = ( ± 1.3 × 10 8 , ± 1.3 × 10 7 , 1 ) . Both figures imply that collapsing singularly degenerate heteroclinic cycles in system (1) create limited cycles rather than strange attractors.
Axioms 13 00625 g006
Figure 7. Phase portraits of system (1) for ( a , c , b ) = ( 4 , 1.8182 , 1 ) and (a) ( x 0 1 , 3 , y 0 1 , 3 , z 0 4 ) = ( ± 0.13 , ± 1.3 , 1.6 ) × 10 7 , ( x 0 3 , 4 , y 0 3 , 4 , z 0 4 ) = ( ± 2.4 , ± 2.41 , 3.3 ) and ( x 0 5 , 6 , y 0 5 , 6 , z 0 5 ) = ( ± 2.39 , ± 2.4 , 3.32 ) , (b) ( x 0 3 , y 0 3 , z 0 4 ) = ( 2.4 , 2.41 , 3.3 ) , ( x 0 5 , y 0 5 , z 0 5 ) = ( 2.39 , 2.4 , 3.32 ) , (c) ( x 0 4 , y 0 4 , z 0 4 ) = ( 2.4 , 2.41 , 3.3 ) and ( x 0 6 , y 0 6 , z 0 5 ) = ( 2.39 , 2.4 , 3.32 ) , showing at least five limit cycles for system (1) when Hopf bifurcation occurs at E ± , i.e., two around E + , two around E and one around E ± .
Figure 7. Phase portraits of system (1) for ( a , c , b ) = ( 4 , 1.8182 , 1 ) and (a) ( x 0 1 , 3 , y 0 1 , 3 , z 0 4 ) = ( ± 0.13 , ± 1.3 , 1.6 ) × 10 7 , ( x 0 3 , 4 , y 0 3 , 4 , z 0 4 ) = ( ± 2.4 , ± 2.41 , 3.3 ) and ( x 0 5 , 6 , y 0 5 , 6 , z 0 5 ) = ( ± 2.39 , ± 2.4 , 3.32 ) , (b) ( x 0 3 , y 0 3 , z 0 4 ) = ( 2.4 , 2.41 , 3.3 ) , ( x 0 5 , y 0 5 , z 0 5 ) = ( 2.39 , 2.4 , 3.32 ) , (c) ( x 0 4 , y 0 4 , z 0 4 ) = ( 2.4 , 2.41 , 3.3 ) and ( x 0 6 , y 0 6 , z 0 5 ) = ( 2.39 , 2.4 , 3.32 ) , showing at least five limit cycles for system (1) when Hopf bifurcation occurs at E ± , i.e., two around E + , two around E and one around E ± .
Axioms 13 00625 g007
Figure 8. For ( a , c ) = ( 6 , 100 ) , b = 8 , 9 and ( x 0 1 , 3 , y 0 1 , 3 , z 0 3 ) = ( ± 0.13 , ± 1.3 , 1.6 ) × 10 7 , phase portraits of system (1), verifying the existence of a pair of heteroclinic orbits to unstable E 0 and stable E ± when c > 0 and 3 b 4 a > 0 .
Figure 8. For ( a , c ) = ( 6 , 100 ) , b = 8 , 9 and ( x 0 1 , 3 , y 0 1 , 3 , z 0 3 ) = ( ± 0.13 , ± 1.3 , 1.6 ) × 10 7 , phase portraits of system (1), verifying the existence of a pair of heteroclinic orbits to unstable E 0 and stable E ± when c > 0 and 3 b 4 a > 0 .
Axioms 13 00625 g008
Table 1. The dynamical behaviors of E 0 .
Table 1. The dynamical behaviors of E 0 .
bacProperty of E 0
<0<0<0A 1D W l o c s and a 2D W l o c u
>0A 3D W l o c u
>0<0A 2D W l o c s and a 1D W l o c u
>0A 1D W l o c s and a 2D W l o c u
>0<0<0A 2D W l o c s and a 1D W l o c u
>0A 1D W l o c s and a 2D W l o c u
>0<0A 3D W l o c s
>0A 2D W l o c s and a 1D W l o c u
Table 2. The dynamical behaviors of E z .
Table 2. The dynamical behaviors of E z .
zaProperty of E z
>0<0A 1D W l o c s , a 1D W l o c c and a 1D W l o c u
>0A 2D W l o c s and a 1D W l o c c
<0<0A 2D W l o c u and a 1D W l o c c
>0A 1D W l o c s , a 1D W l o c c and a 1D W l o c u
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ke, G.; Pan, J.; Hu, F.; Wang, H. Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System. Axioms 2024, 13, 625. https://doi.org/10.3390/axioms13090625

AMA Style

Ke G, Pan J, Hu F, Wang H. Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System. Axioms. 2024; 13(9):625. https://doi.org/10.3390/axioms13090625

Chicago/Turabian Style

Ke, Guiyao, Jun Pan, Feiyu Hu, and Haijun Wang. 2024. "Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System" Axioms 13, no. 9: 625. https://doi.org/10.3390/axioms13090625

APA Style

Ke, G., Pan, J., Hu, F., & Wang, H. (2024). Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System. Axioms, 13(9), 625. https://doi.org/10.3390/axioms13090625

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