# Sequence of Routes to Chaos in a Lorenz-Type System

Authors: Fangyan Yang; Yongming Cao; Lijuan Chen; Qingdu Li
Journal: Discrete Dynamics in Nature and Society
 (2020)

---

## Abstract

This paper reports a new bifurcation pattern observed in a Lorenz-type system. The pattern is composed of a main bifurcation route to chaos (n
=
1) and a sequence of sub-bifurcation routes with n
=
3
,
4
,
5
,
…
,
14 isolated sub-branches to chaos. When n is odd, the n isolated sub-branches are from a period-n limit cycle, followed by twin period-n limit cycles via a pitchfork bifurcation, twin chaotic attractors via period-doubling bifurcations, and a symmetric chaotic attractor via boundary crisis. When n is even, the n isolated sub-branches are from twin period-n
/
2 limit cycles, which become twin chaotic attractors via period-doubling bifurcations. The paper also shows that the main route and the sub-routes can coexist peacefully by studying basins of attraction.


---

## Body

## 1. Introduction

Nonlinear phenomena in physics and engineering are often described by quadratic systems, such as the famous Lorenz system [1], the Nose-Hoover equations [2, 3], the brushless DC motor system [4], the pendulum system [5], etc. A fundamental requirement for such systems is to understand all their behaviors. However, it is a hard task, even for two dimensional (2D) cases (e.g., the Hilbert’s sixteenth problem). When the dimension is three or more, the existence of chaos makes it even harder to establish a general theory. Therefore, scientists from various backgrounds devoted themselves to exploiting nonlinear dynamics, especially chaos in Lorenz-type systems, in recent years [6–11].The most interesting feature of chaos is perhaps that it looks random but simple sequences and patterns are hidden behind. In the Sarkovskii theorem, the Sarkovskii ordering of the natural numbers implies the famous corollary by Li and York: if a continuous map defined on the real line has a period-3 point, then it must have periodic points of every other periods [12]. In the route to chaos, there is usually a sequence of period-doubling (PD) bifurcations as a parameter is adjusted. In Smale’s horseshoe [13] or a topological horseshoe [14–16], chaos can be explained by a sequence of operations: squishing, stretching, and folding. Near the quadratic homoclinic tangency, Gavrilov, Silnikov, and Newhouse proved the existence of long period sinks at an infinite sequence of parameter values [17, 18]. For better understanding of the complexity of chaos, whether there are more simple patterns or sequences is always a fundamental question to ask.In 1992, Petrov et al. found mixed-mode oscillations in a chemical system [19]. They found a pattern composed of a sequence of n isolated sub-branches in a one-parameter bifurcation diagram, where for each n
=
2
,
3
,
4
,
…, the n isolated sub-branches initiate from a period-n limit cycle and lead to chaos via period doubling bifurcations. Recently, similar pattern has also been found in hybrid dynamical system, i.e., the simplest walking model [20, 21]. More recent study suggests that this regular pattern may also exist in the compass-like walking model [22] and the simple walking model with upper body [23]. Although the period of the initial limit cycle increases with step size 1, the pattern is different from the period-adding phenomena reported in [24–26] since there is no parameter overlap between these branches and the period-doubling and period-adding happen in different parameter-changing directions. For this phenomenon, the maximal n found is 7, and the systems are all nonsymmetrical. Therefore, it is interesting to see the situation with a larger n for a symmetrical system.Therefore, this paper will study a Lorenz-type system and show a sequence sub-branch withn up to 14. When n is odd, the n isolated sub-branches are from a period-n limit cycle, otherwise the n sub-branches are from twin period-n
/
2 limit cycles. Although all these branches lead to chaos as a parameter increases, this pattern is completely different (see Figure 1) from the existing literature, according to the best of the authors’ knowledge.Figure 1
Sequence of routes withn sub-branches to chaos.
(a)
(b)The rest of paper is organized as follows: Section2 introduces the Lorenz-type system, studies its symmetry, equilibria and their stability, and briefly explores the nonlinear dynamics; Section 3 carries out detailed numerical research on the sequence of bifurcation branches; Section 4 draws the conclusion.


## 2. The Lorenz Type System

The Lorenz type system is given by following equations:(1)
x
˙
=
−
a
x
−
b
y
+
y
z
,
y
˙
=
−
c
x
−
x
z
,
z
˙
=
−
x
2
+
x
y
−
d
z
,wherex, y and z are the state variables, and a, b, c and d (when d
≠
0, this parameter can be eliminated with replacement x
→
d
x
̂, y
→
d
y
̂, z
→
d
z
̂, a
→
a
̂
d, b
→
b
̂
d, c
→
c
̂
d and t
→
t
̂
/
d) are positive parameters. Comparing with the famous Lorenz system, although there are some different terms, it has similar reflection symmetry, dissipation, equilibria, and even chaotic attractor. Obviously, system (1) is invariant under the same transformation,(2)
x
,
y
,
z
↔
−
x
,
−
y
,
z
.Therefore, ifx
t
,
y
t
,
z
t is a solution of (1), then −
x
t
,
−
y
t
,
z
t is also a solution. Since the Jacobian matrix of (1) at x
,
y
,
z is,(3)
J
=
−
a
z
−
b
y
−
c
−
z
0
−
x
y
−
2
x
x
−
d
.It is easy to see that the system is dissipative because there exists an exponential contraction rate,(4)
∇
V
=
∂
x
˙
∂
x
+
∂
y
˙
∂
y
+
∂
z
˙
∂
z
=
−
a
+
d
<
0
.SinceΔ
=
c
d
/
b
+
c
a
+
b
+
c
>
0, system (1) has three equilibria. The first one is the origin O, and the eigenvalues of the Jacobian matrix (3) at O are,(5)
λ
1
=
−
a
2
+
a
2
+
4
b
c
2
,
λ
2
=
−
a
2
−
a
2
+
4
b
c
2
and
λ
3
=
−
d
.Since all parameters are positive, we haveλ
1
>
0, λ
2
<
0 and λ
3
<
0, which means that O is unstable. The other two equilibria are,(6)
O
−
=
−
b
+
c
Δ
,
a
Δ
,
−
c
and
O
+
=
b
+
c
Δ
,
−
a
Δ
,
−
c
,which are obviously symmetric to each other under transformation (2). When the parameters take the following values,(7)
a
∈
3
,
30
,
b
=
3
,
c
=
63.5
,
and
d
=
5.They are both unstable. For example, whena
=
15, the eigenvalues at both O
− and O
+ are λ
1
,
2
≈
11.09
±
30.58
i, and λ
3
≈
−
39.18.To explore the nonlinear dynamics of this system, we compute Lyapunov exponents and a bifurcation diagram with the same initial condition (0, 0.01, 0.01) asa varies from 3 to 30. For Lyapunov exponents, we use the QR-based Jacobian method and the final time takes 10,000 to ensure a good accuracy. For the bifurcation diagram, we first take Π
=
{
x
,
y
,
z
|
x
=
0
and
x
˙
<
0
} as a Poincaré cross-section plane, and then define a Poincaré map h
:
Π
→
Π as follows: For each x
≜
y
,
z
∈
Π, h
x is taken to be the first return point in Π under the flow of system (1) with the initial condition x. To exclude transient behaviors of iterating h, we discard the first 1000 points, and then save the following 1000 points.The numerical results are shown in Figure2 with double axes where the left axis indicates the bifurcation diagram and the right axis indicates Lyapunov exponents. Since the last Lyapunov exponent is always a large negative number, it is not displayed. Figure 2 suggests that the Lyapunov exponents and the bifurcation diagram match very well and both show rich complex behaviors.Figure 2
The bifurcation diagram and the Lyapunov exponents by varyinga.When3
≤
a
<
9.612, system (1) has a limit cycle with period-1. A typical phase portrait at a
=
5 is illustrated as the first row in Figure 3 with three different view angles. As a increases from 9.621, there seems a period-doubling route to chaos. A typical chaotic attractor at a
=
11.2 is shown in the second row in Figure 3, and its three Lyapunov exponents are 2.90, 0.00, and −19.10, respectively. While a keeps increasing, the largest Lyapunov exponent also increases. When a
≈
11.36, the chaotic attractor suddenly becomes a large chaotic attractor. A typical large chaotic attractor at a
=
15 is shown in the last row in Figure 3, and its three Lyapunov exponents are 6.94, 0.00, and −26.94, respectively. When a increases from 15, the largest Lyapunov exponent may jump to zero occasionally, which indicates some periodic windows in the bifurcation diagram.Figure 3
Typical phase portraits ata
=
5
,
11.2
and
15.


## 3. Sequence of Routes to Chaos

If one carefully looks at the bifurcation diagram in Figure2, two unexpected issues will be found: (1) The period-doubling route seems broken near a
=
10.69; (2) There are many points located outside of the main bifurcation branch near a
=
10. The following detailed study will lead us to a sequence of sub-bifurcation routes to chaos, as shown in Figure 1(a).


### 3.1. The Main BranchB
1 to Chaos from the Period-1 Limit Cycle

For the first issue, we trace the bifurcation diagram by continuation. The results are shown in Figure4(e) with the left axis, where the bifurcation branches in magenta and blue indicate different attractors. When a
=
9
.5, there is only one limit cycle with period-1, as shown in Figure 4(a). At a
≈
9
.612, this cycle bifurcates into two different period-1 limit cycles. The new-born “twin” limit cycles are shown in Figure 4(b). Here “twin” means that they are the same under the transformation (2). As a increases, each limit cycle becomes a chaotic attractor via a sequence of period-doubling bifurcations. A typical phase portrait of the two chaotic attractors is shown in Figure 4(d). They are also symmetric to each other, and their Lyapunov exponents are 1.65, 0, and −17.65. They can peacefully coexist, and their basins of attraction are shown in Figure 4(f), which suggest that the boundaries are fractal. At last, the twin chaotic attractors merge into one chaotic attractor at a
≈
11.08 due to a boundary crisis. The merged attractor at a
=
11.2 has already been shown in the middle row of Figure 3.Figure 4
The detailed dynamics and bifurcations of the main branch.
(a)
(b)
(c)
(d)
(e)
(f)To examine how the original cycle bifurcates into the twin cycles, we compute the eigenvalues of the Jacobian of the Poincaré maph along the period-1 cycle. The eigenvalues are shown in Figure 4(e) with the right axis, where the blue dot line indicates the period-1 cycle after the bifurcation, and the green solid and dot lines indicate the absolute eigenvalues before and after the bifurcation, respectively. Since one eigenvalue crosses the unit circle from inside to outside at 1.0, it belongs to a pitchfork bifurcation, which is introduced by the reflection symmetry of the system.


### 3.2. Sub-BranchesB
3 to Chaos from a Period-3 Limit Cycle

For the second issue, we trace those points out of the main branch via continuation while adjustinga up and down, and find a sequence of sub-bifurcation branches consisted of new attractors, which can peacefully coexist with the attractors in the main branch.When9.6568
<
a
<
10.0903, although there are a pair of twin chaotic attractors in the main branch, we find three new sub-branches of bifurcation diagrams. One of them is shown in detail with Figure 5(e). It starts with a period-3 limit cycle, which is mirror symmetric under transformation (2). This period-3 cycle can peacefully coexist with the twin period-1 cycles in the main branch, as shown in Figure 5(a) for a
=
9
.7531, where the initial values are (0, 7.7231, −68.5635) and (0, ±12.5342, −78.6258), respectively. Like Figure 4(e), eigenvalues of the period-3 cycle are shown in Figure 5(e), which suggests that it is introduced by a fold bifurcation (also called a saddle-node bifurcation).Figure 5
The period-3 limit cycle, its bifurcations and basin of attraction.
(a)
(b)
(c)
(d)
(e)
(f)Whena increases, there is also a pitchfork bifurcation as shown in Figure 5(e). The period-3 cycle bifurcates into a pair of new twin period-3 cycles at a
≈
9
.8702. A phase portrait of the twin cycles is depicted in Figure 5(b) at a
≈
9
.9244. Note that the new twin period-3 cycles can also peacefully coexist with the twin period-1 cycles. Their basins of attraction of h are shown in Figure 5(f), where P
1
≈
12.5020
,
−
78.6634 and P
2
≈
6.0350
,
−
65.0803 are from the twin period-3 cycles, and P
3
≈
9.1301
,
−
71.6131 and P
4
≈
7.440
9
,
−
68.0668 are from the twin period-1 cycles. It is clear to see that all basins have fractional structures, and the boundaries are mixed together in many places, where an initial value could be trapped by any attractor.Like the main bifurcation branch, the twin period-3 cycles also become a pair of twin chaotic attractors via period-doubling bifurcations asa keeps increasing. A typical phase portrait of the twin chaotic attractors are shown in Figure 5(c) at a
=
10. When a
≈
10
.0145, the twin chaotic attractors suddenly merge into one chaotic attractor via a boundary crisis, as shown in Figure 5(d) at a
=
10.06.


### 3.3. Sub-BranchesB
4 to Chaos from a Pair of Twin Period-2 Limit Cycles

Whena takes a smaller range, i.e., 9.4370
<
a
<
9.5852, although there is a period-1 limit cycle in the main branch, we find four new sub-branches of bifurcations. Unlike the three sub-branches studied above, they are NOT from one limit cycle with period-4 but from a pair of twin limit cycles with period-2. A typical phase portrait of the twin cycles at a
=
9.4741 is shown in Figure 6(a). And a detailed branch of the bifurcation diagrams is shown in Figure 6(c) with double axes. It is easy to see that the twin limit cycles are also introduced by a fold bifurcation. The basins of attraction are shown in Figure 6(d), where P
1
≈
12
.
8113
,
−
84.3293 and P
2
=
9.0
434
,
−
75.4186 are from the twin cycles with period-2, and P
3
=
8.3636
,
−
69.7783 is from the main cycle with period-1. Obviously, the three limit cycles can coexist peacefully, and their basins have fractional structures too. When a keeps increasing from 9.4741, the twin period-2 cycles become a pair of twin chaotic attractors via period-doubling bifurcations. A phase portrait of the twin chaotic attractors at a
=
9
.5801 are shown in Figure 6(b), where the initial values are (0, ±12.9744, −84.9422). And they have similar basins shown in Figure 6(e).Figure 6
The period-2 limit cycle, its bifurcations and basin of attraction.
(a)
(b)
(c)
(d)
(e)


### 3.4. Sub-BranchesB
5 to Chaos from a Period-5 Limit Cycle

Whena takes another smaller range, i.e., 9.2401
<
a
<
9.2985, we find five sub-branches of bifurcations outside the main branch. Since five is an odd number like three, the bifurcation phenomenon is similar to the diagram from the limit cycles with period-3, as shown in Figure 7. First, there is a limit cycle (Figure 7(a)) with period-5, which is also introduced by a fold bifurcation (Figure 7(e)), then this cycle bifurcates into a pair of twin limit cycles (Figure 7(b)) with period-5 via a pitchfork bifurcation, after that they become a pair of twin chaotic attractors (Figure 7(c)) via period-doubling bifurcation, at last the twin attractors merge into one larger chaotic attractor (Figure 7(d)).Figure 7
Limit cycle with period-5 and its bifurcation.
(a)
(b)
(c)
(d)
(e)


### 3.5. Sub-BranchesB
6 to Chaos from a Pair of Twin Period-3 Limit Cycles

Whena keeps taking a smaller range, i.e., 9.0636
<
a
<
9.0895, we find six new sub-branches of bifurcations. Since six is an even number like four, the bifurcation phenomenon is similar to the four sub-branches from the twin limit cycles with period-2, as shown in Figure 8. First, there are a pair of twin limit cycles with period-3 (Figure 8(a)), which are also introduced by a fold bifurcation (Figure 8(c)), and then the twin cycles become a pair of twin chaotic attractors (Figure 8(b)) via period-doubling bifurcations.Figure 8
Limit cycle with period-6 and its bifurcations.
(a)
(b)
(c)


### 3.6. Sequence of Routes to Chaos with More Sub-Branches

According to the above observation, we have a sequence of sub-branchesB
n, n
=
3
,
4
,
5
,
6, which can peacefully coexist with the main branch B
1. When the integer n is odd, e.g., 1, 3, 5, the n sub-branches start with a limit cycle with period-n; otherwise, they start with a pair of twin limit cycles with period-n
/
2. This reminds us of the sequence of bifurcations found in the simplest passive walking model [19], where for each n
∈
{
3
,
4
,
5
,
6
,
7
}, there are n sub-branches of period-doubling bifurcations from the period-n gait, and they all can coexist peacefully with the main period-doubling branch from period-1.Although the bifurcation pattern here is more complex than [19], the sequence is regular enough to predict more sub-branches if they exist. Inspired by Feigenbaum constants, we suppose the following ratio tends to a limit, i.e.,(8)
lim
n
→
∞
r
n
−
r
n
−
1
r
n
−
1
−
r
n
−
2
=
η
,wherer
n represent the first bifurcation point in B
n. If n is odd, r
n is a pitchfork bifurcation point, otherwise a period doubling bifurcation. Then we predict an approximate value of the next bifurcation point by,(9)
r
n
+
1
≈
r
n
+
r
n
−
r
n
−
1
2
r
n
−
1
−
r
n
−
2
.By searching basins of attraction near this predicted point with the heterogeneous algorithm accelerated by GPU computing proposed in [18], many new sub-branches are found, as shown in Figure 1(a), which successfully extend the sequence up to 14. As listed in Table 1, for the n sub-branches B
n, x
n indicates an approximate point of the new limit cycle at a
=
a
n, which is induced by a fold bifurcation at a
=
s
n.Table 1
Sequence of bifurcation points.
n
x
n
a
n
s
n
r
n
r
n
−
r
n
−
1
r
n
−
1
−
r
n
−
2
r
n
−
s
n
r
n
−
1
−
s
n
−
1
3
(12.6421, −78.8718)
9.8005
9.655842
9.871091
—
—
4
(12.7825, −84.2335)
9.4646
9.436742
9.504714
—
3.1699
5
(12.5709, −87.4911)
9.2498
9.239972
9.266180
1.5359
2.5995
6
(8.4696, −76.3535)
9.0699
9.063520
9.075047
1.2481
2.2701
7
(8.2565, −76.6158)
8.9023
8.899685
8.905179
1.1246
2.0952
8
(11.7836, −92.3385)
8.7460
8.744490
8.747220
1.0754
2.0000
9
(7.9555, −77.0160)
8.5972
8.596478
8.597911
1.0577
1.9064
10
(11.3884, −93.8690)
8.4558
8.455405
8.456170
1.0533
1.8867
11
(2.1312, −50.9839)
8.3216
8.321398
8.321815
1.0551
1.8201
12
(1.8821, −49.7448)
8.1947
8.194608
8.194839
1.0581
1.8052
13
(1.6619, −48.5650)
8.0751
8.0750751
8.0752041
1.0614
1.7907
14
(13.6614 , −110.2840)
7.9627
7.9626959
7.9627687
1.0640
1.7917The detailed bifurcations are shown in Figure1(b), which all follows the known pattern: When n is odd, like the branch-3 shown in Figure 5, the n sub-branches first start with a period-n limit cycle introduced by a fold bifurcation, then bifurcate into a pair of twin limit cycles with the same period, after that the twin cycles become a pair of twin chaotic attractors via period doubling bifurcations, at last the twin chaotic attractors merge into one larger chaotic attractor. When n is even, like the branch-4 shown in Figure 6, the n sub-branches start with a pair of twin limit cycles with period-n
/
2, which are also introduced by a fold bifurcation, then they become a chaotic attractor via period-doubling bifurcations.Whenn get lager, the parameter range becomes smaller. For the n sub-branches in B
n, the parameter range of the initial limit cycles is r
n
−
s
n. The following range ratio(10)
r
n
−
s
n
r
n
−
1
−
s
n
−
1
,is listed in Table1. It seems to have a limit under a small error, which should be close to 1.79. This evidence suggests that the parameter range will get smaller very quickly as n grows. Due to unavoidable numerical errors and time-consuming computing, it becomes extremely hard to find new sub-branches with n
≥
15.


### 3.7. Connections with Periodic Windows of the Main Branch

The sequence of the sub-branches come from a sequence of the initial limit cycles with the following period,(11)
3
,
2
,
5
,
3
,
7
,
4
,
9
,
5
,
11
,
6
,
13
,
7
,
…This sequence is different from the Sarkovskii’s ordering. And it is also different from the period-adding phenomenon reported in the existing literature, according to the best of the authors’ knowledge. Therefore, it is significant to see where the limit cycles come from.Note that a period-n cycle of h is a fixed point of h
n. According to degree theory, as the parameter changes, the only possibilities for this fix point are: (a) coming/leaving through the boundary of the y
,
z domain, (b) being born (splitting)/disappearing (merging) in pairs. Therefore, we can track the limit cycles, with numerical continuation algorithms, back before the fold bifurcation.The results are shown in Figure9(a). It is interesting to see that all the cycles turn back to the region of chaos, which implies that the cycles might be born in the chaotic attractor. However, when we track the cycles with very small step size, it is surprising to see that they connect with some small periodic windows (see Figure 2) in the main bifurcation diagram. For example, after two-fold bifurcations, the period-3 cycle connects with three sub-branches in a period-3 windows. A detailed bifurcation diagram of one branch is shown in Figure 9(b), it looks like the reverse of the branch in Figure 4(e). Figure 9(c) also shows such a connection for n
=
5. When n is an even number, the situation is similar. For example, the twin period-n
/
2 cycles with n
=
4 or n
=
6 are also a periodic window after two-fold bifurcations. The window also contains a pair of twin limit cycles with period-n
/
2, which bifurcate into a pair of twin chaotic attractors via period-doubling bifurcations, as shown in Figure 9(e).Figure 9
Connections between the sub-branches and the periodic windows.
(a)
(b)
(c)
(d)
(e)The above evidence suggests that each branch in our sequence on the left could connect with a branch in a periodic window on the right. Such phenomenon is much different from the phenomenon reported in [19], where the cycles look passing through the main chaotic attractor bifurcated from the period-1 limit cycle.


## 3.1. The Main BranchB
1 to Chaos from the Period-1 Limit Cycle

For the first issue, we trace the bifurcation diagram by continuation. The results are shown in Figure4(e) with the left axis, where the bifurcation branches in magenta and blue indicate different attractors. When a
=
9
.5, there is only one limit cycle with period-1, as shown in Figure 4(a). At a
≈
9
.612, this cycle bifurcates into two different period-1 limit cycles. The new-born “twin” limit cycles are shown in Figure 4(b). Here “twin” means that they are the same under the transformation (2). As a increases, each limit cycle becomes a chaotic attractor via a sequence of period-doubling bifurcations. A typical phase portrait of the two chaotic attractors is shown in Figure 4(d). They are also symmetric to each other, and their Lyapunov exponents are 1.65, 0, and −17.65. They can peacefully coexist, and their basins of attraction are shown in Figure 4(f), which suggest that the boundaries are fractal. At last, the twin chaotic attractors merge into one chaotic attractor at a
≈
11.08 due to a boundary crisis. The merged attractor at a
=
11.2 has already been shown in the middle row of Figure 3.Figure 4
The detailed dynamics and bifurcations of the main branch.
(a)
(b)
(c)
(d)
(e)
(f)To examine how the original cycle bifurcates into the twin cycles, we compute the eigenvalues of the Jacobian of the Poincaré maph along the period-1 cycle. The eigenvalues are shown in Figure 4(e) with the right axis, where the blue dot line indicates the period-1 cycle after the bifurcation, and the green solid and dot lines indicate the absolute eigenvalues before and after the bifurcation, respectively. Since one eigenvalue crosses the unit circle from inside to outside at 1.0, it belongs to a pitchfork bifurcation, which is introduced by the reflection symmetry of the system.


## 3.2. Sub-BranchesB
3 to Chaos from a Period-3 Limit Cycle

For the second issue, we trace those points out of the main branch via continuation while adjustinga up and down, and find a sequence of sub-bifurcation branches consisted of new attractors, which can peacefully coexist with the attractors in the main branch.When9.6568
<
a
<
10.0903, although there are a pair of twin chaotic attractors in the main branch, we find three new sub-branches of bifurcation diagrams. One of them is shown in detail with Figure 5(e). It starts with a period-3 limit cycle, which is mirror symmetric under transformation (2). This period-3 cycle can peacefully coexist with the twin period-1 cycles in the main branch, as shown in Figure 5(a) for a
=
9
.7531, where the initial values are (0, 7.7231, −68.5635) and (0, ±12.5342, −78.6258), respectively. Like Figure 4(e), eigenvalues of the period-3 cycle are shown in Figure 5(e), which suggests that it is introduced by a fold bifurcation (also called a saddle-node bifurcation).Figure 5
The period-3 limit cycle, its bifurcations and basin of attraction.
(a)
(b)
(c)
(d)
(e)
(f)Whena increases, there is also a pitchfork bifurcation as shown in Figure 5(e). The period-3 cycle bifurcates into a pair of new twin period-3 cycles at a
≈
9
.8702. A phase portrait of the twin cycles is depicted in Figure 5(b) at a
≈
9
.9244. Note that the new twin period-3 cycles can also peacefully coexist with the twin period-1 cycles. Their basins of attraction of h are shown in Figure 5(f), where P
1
≈
12.5020
,
−
78.6634 and P
2
≈
6.0350
,
−
65.0803 are from the twin period-3 cycles, and P
3
≈
9.1301
,
−
71.6131 and P
4
≈
7.440
9
,
−
68.0668 are from the twin period-1 cycles. It is clear to see that all basins have fractional structures, and the boundaries are mixed together in many places, where an initial value could be trapped by any attractor.Like the main bifurcation branch, the twin period-3 cycles also become a pair of twin chaotic attractors via period-doubling bifurcations asa keeps increasing. A typical phase portrait of the twin chaotic attractors are shown in Figure 5(c) at a
=
10. When a
≈
10
.0145, the twin chaotic attractors suddenly merge into one chaotic attractor via a boundary crisis, as shown in Figure 5(d) at a
=
10.06.


## 3.3. Sub-BranchesB
4 to Chaos from a Pair of Twin Period-2 Limit Cycles

Whena takes a smaller range, i.e., 9.4370
<
a
<
9.5852, although there is a period-1 limit cycle in the main branch, we find four new sub-branches of bifurcations. Unlike the three sub-branches studied above, they are NOT from one limit cycle with period-4 but from a pair of twin limit cycles with period-2. A typical phase portrait of the twin cycles at a
=
9.4741 is shown in Figure 6(a). And a detailed branch of the bifurcation diagrams is shown in Figure 6(c) with double axes. It is easy to see that the twin limit cycles are also introduced by a fold bifurcation. The basins of attraction are shown in Figure 6(d), where P
1
≈
12
.
8113
,
−
84.3293 and P
2
=
9.0
434
,
−
75.4186 are from the twin cycles with period-2, and P
3
=
8.3636
,
−
69.7783 is from the main cycle with period-1. Obviously, the three limit cycles can coexist peacefully, and their basins have fractional structures too. When a keeps increasing from 9.4741, the twin period-2 cycles become a pair of twin chaotic attractors via period-doubling bifurcations. A phase portrait of the twin chaotic attractors at a
=
9
.5801 are shown in Figure 6(b), where the initial values are (0, ±12.9744, −84.9422). And they have similar basins shown in Figure 6(e).Figure 6
The period-2 limit cycle, its bifurcations and basin of attraction.
(a)
(b)
(c)
(d)
(e)


## 3.4. Sub-BranchesB
5 to Chaos from a Period-5 Limit Cycle

Whena takes another smaller range, i.e., 9.2401
<
a
<
9.2985, we find five sub-branches of bifurcations outside the main branch. Since five is an odd number like three, the bifurcation phenomenon is similar to the diagram from the limit cycles with period-3, as shown in Figure 7. First, there is a limit cycle (Figure 7(a)) with period-5, which is also introduced by a fold bifurcation (Figure 7(e)), then this cycle bifurcates into a pair of twin limit cycles (Figure 7(b)) with period-5 via a pitchfork bifurcation, after that they become a pair of twin chaotic attractors (Figure 7(c)) via period-doubling bifurcation, at last the twin attractors merge into one larger chaotic attractor (Figure 7(d)).Figure 7
Limit cycle with period-5 and its bifurcation.
(a)
(b)
(c)
(d)
(e)


## 3.5. Sub-BranchesB
6 to Chaos from a Pair of Twin Period-3 Limit Cycles

Whena keeps taking a smaller range, i.e., 9.0636
<
a
<
9.0895, we find six new sub-branches of bifurcations. Since six is an even number like four, the bifurcation phenomenon is similar to the four sub-branches from the twin limit cycles with period-2, as shown in Figure 8. First, there are a pair of twin limit cycles with period-3 (Figure 8(a)), which are also introduced by a fold bifurcation (Figure 8(c)), and then the twin cycles become a pair of twin chaotic attractors (Figure 8(b)) via period-doubling bifurcations.Figure 8
Limit cycle with period-6 and its bifurcations.
(a)
(b)
(c)


## 3.6. Sequence of Routes to Chaos with More Sub-Branches

According to the above observation, we have a sequence of sub-branchesB
n, n
=
3
,
4
,
5
,
6, which can peacefully coexist with the main branch B
1. When the integer n is odd, e.g., 1, 3, 5, the n sub-branches start with a limit cycle with period-n; otherwise, they start with a pair of twin limit cycles with period-n
/
2. This reminds us of the sequence of bifurcations found in the simplest passive walking model [19], where for each n
∈
{
3
,
4
,
5
,
6
,
7
}, there are n sub-branches of period-doubling bifurcations from the period-n gait, and they all can coexist peacefully with the main period-doubling branch from period-1.Although the bifurcation pattern here is more complex than [19], the sequence is regular enough to predict more sub-branches if they exist. Inspired by Feigenbaum constants, we suppose the following ratio tends to a limit, i.e.,(8)
lim
n
→
∞
r
n
−
r
n
−
1
r
n
−
1
−
r
n
−
2
=
η
,wherer
n represent the first bifurcation point in B
n. If n is odd, r
n is a pitchfork bifurcation point, otherwise a period doubling bifurcation. Then we predict an approximate value of the next bifurcation point by,(9)
r
n
+
1
≈
r
n
+
r
n
−
r
n
−
1
2
r
n
−
1
−
r
n
−
2
.By searching basins of attraction near this predicted point with the heterogeneous algorithm accelerated by GPU computing proposed in [18], many new sub-branches are found, as shown in Figure 1(a), which successfully extend the sequence up to 14. As listed in Table 1, for the n sub-branches B
n, x
n indicates an approximate point of the new limit cycle at a
=
a
n, which is induced by a fold bifurcation at a
=
s
n.Table 1
Sequence of bifurcation points.
n
x
n
a
n
s
n
r
n
r
n
−
r
n
−
1
r
n
−
1
−
r
n
−
2
r
n
−
s
n
r
n
−
1
−
s
n
−
1
3
(12.6421, −78.8718)
9.8005
9.655842
9.871091
—
—
4
(12.7825, −84.2335)
9.4646
9.436742
9.504714
—
3.1699
5
(12.5709, −87.4911)
9.2498
9.239972
9.266180
1.5359
2.5995
6
(8.4696, −76.3535)
9.0699
9.063520
9.075047
1.2481
2.2701
7
(8.2565, −76.6158)
8.9023
8.899685
8.905179
1.1246
2.0952
8
(11.7836, −92.3385)
8.7460
8.744490
8.747220
1.0754
2.0000
9
(7.9555, −77.0160)
8.5972
8.596478
8.597911
1.0577
1.9064
10
(11.3884, −93.8690)
8.4558
8.455405
8.456170
1.0533
1.8867
11
(2.1312, −50.9839)
8.3216
8.321398
8.321815
1.0551
1.8201
12
(1.8821, −49.7448)
8.1947
8.194608
8.194839
1.0581
1.8052
13
(1.6619, −48.5650)
8.0751
8.0750751
8.0752041
1.0614
1.7907
14
(13.6614 , −110.2840)
7.9627
7.9626959
7.9627687
1.0640
1.7917The detailed bifurcations are shown in Figure1(b), which all follows the known pattern: When n is odd, like the branch-3 shown in Figure 5, the n sub-branches first start with a period-n limit cycle introduced by a fold bifurcation, then bifurcate into a pair of twin limit cycles with the same period, after that the twin cycles become a pair of twin chaotic attractors via period doubling bifurcations, at last the twin chaotic attractors merge into one larger chaotic attractor. When n is even, like the branch-4 shown in Figure 6, the n sub-branches start with a pair of twin limit cycles with period-n
/
2, which are also introduced by a fold bifurcation, then they become a chaotic attractor via period-doubling bifurcations.Whenn get lager, the parameter range becomes smaller. For the n sub-branches in B
n, the parameter range of the initial limit cycles is r
n
−
s
n. The following range ratio(10)
r
n
−
s
n
r
n
−
1
−
s
n
−
1
,is listed in Table1. It seems to have a limit under a small error, which should be close to 1.79. This evidence suggests that the parameter range will get smaller very quickly as n grows. Due to unavoidable numerical errors and time-consuming computing, it becomes extremely hard to find new sub-branches with n
≥
15.


## 3.7. Connections with Periodic Windows of the Main Branch

The sequence of the sub-branches come from a sequence of the initial limit cycles with the following period,(11)
3
,
2
,
5
,
3
,
7
,
4
,
9
,
5
,
11
,
6
,
13
,
7
,
…This sequence is different from the Sarkovskii’s ordering. And it is also different from the period-adding phenomenon reported in the existing literature, according to the best of the authors’ knowledge. Therefore, it is significant to see where the limit cycles come from.Note that a period-n cycle of h is a fixed point of h
n. According to degree theory, as the parameter changes, the only possibilities for this fix point are: (a) coming/leaving through the boundary of the y
,
z domain, (b) being born (splitting)/disappearing (merging) in pairs. Therefore, we can track the limit cycles, with numerical continuation algorithms, back before the fold bifurcation.The results are shown in Figure9(a). It is interesting to see that all the cycles turn back to the region of chaos, which implies that the cycles might be born in the chaotic attractor. However, when we track the cycles with very small step size, it is surprising to see that they connect with some small periodic windows (see Figure 2) in the main bifurcation diagram. For example, after two-fold bifurcations, the period-3 cycle connects with three sub-branches in a period-3 windows. A detailed bifurcation diagram of one branch is shown in Figure 9(b), it looks like the reverse of the branch in Figure 4(e). Figure 9(c) also shows such a connection for n
=
5. When n is an even number, the situation is similar. For example, the twin period-n
/
2 cycles with n
=
4 or n
=
6 are also a periodic window after two-fold bifurcations. The window also contains a pair of twin limit cycles with period-n
/
2, which bifurcate into a pair of twin chaotic attractors via period-doubling bifurcations, as shown in Figure 9(e).Figure 9
Connections between the sub-branches and the periodic windows.
(a)
(b)
(c)
(d)
(e)The above evidence suggests that each branch in our sequence on the left could connect with a branch in a periodic window on the right. Such phenomenon is much different from the phenomenon reported in [19], where the cycles look passing through the main chaotic attractor bifurcated from the period-1 limit cycle.


## 4. Conclusions

We have studied a Lorenz-type system and revealed a sequence of route to chaos in the following new pattern. For an odd numbern
=
3
,
5
,
…
,
13, there are n isolated sub-bifurcation branches started with a period-n limit cycle and followed by twin period-n limit cycles via a pitchfork bifurcation, twin chaotic attractors via a sequence of period-doubling bifurcations and a symmetric chaotic attractor via boundary crisis. For even number n
=
4
,
6
,
…
,
14, there are n sub-branches started with a period-n
/
2 limit cycle, which becomes twin chaotic attractor via period-doubling bifurcations.In the above pattern, many chaotic attractors are born from limit cycles with high periods. This is much different from most existing results where there is just one chaotic attractor from a period-1 limit cycle. Since the sub-branches withn
≥
3 not only coexist with the main branch with n
=
1, but also connect with its periodic windows via two-fold bifurcations, the pattern revealed in this paper is also different from the period-adding phenomena. In comparison with the phenomenon in the walking models and the chemical system, the sequence here is much longer [21], which is more helpful for further study. In addition, since all equilibria are unstable for the given parameter in (3), we only observed self-excited attractors. Basin of attraction shows that there should be no hidden attractor like the other Lorenz-like system studied in [27, 28].We hope the studies carried out in the present paper may bring a new insight into the mechanism of chaos and coexisting attractors. For this purpose, it will also be interesting to find this pattern in real physical systems or simple 2D maps.



---

*Source: 3162170-2020-01-23.xml*