Next Article in Journal
Identification of Inhomogeneities: The Selected Coordinate Descent Method Applied in the Drilling Area
Next Article in Special Issue
On the Double-Zero Bifurcation of Jerk Systems
Previous Article in Journal
The Mean Value Theorem in the Context of Generalized Approach to Differentiability
Previous Article in Special Issue
Exploring Limit Cycle Bifurcations in the Presence of a Generalized Heteroclinic Loop
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Hopf and Fold Bifurcations of Jerk Systems

by
Cristian Lăzureanu
1,2,*,† and
Jinyoung Cho
1,2,†
1
Department of Mathematics, Politehnica University Timişoara, P-ta Victoriei 2, 300006 Timişoara, Romania
2
Department of Mathematics, Faculty of Mathematics and Computer Science, West University of Timişoara, Parvan Blv. 4, 300223 Timişoara, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(20), 4295; https://doi.org/10.3390/math11204295
Submission received: 11 September 2023 / Revised: 12 October 2023 / Accepted: 13 October 2023 / Published: 15 October 2023
(This article belongs to the Special Issue Application of Mathematical Method and Models in Dynamic System)

Abstract

:
In this paper we consider a jerk system x ˙ = y , y ˙ = z , z ˙ = j ( x , y , z , α ) , where j is an arbitrary smooth function and α is a real parameter. Using the derivatives of j at an equilibrium point, we discuss the stability of that point, and we point out some local codim-1 bifurcations. Moreover, we deduce jerk approximate normal forms for the most common fold bifurcations.
MSC:
70K20; 70K45; 70K50

1. Introduction

A jerk equation is a third-order ordinary differential equation of the form
x = j ( x , x ˙ , x ¨ ) ,
where the rate of change of acceleration of the motion, x ( t ) , is called jerk [1]. Its equivalent form is the jerk system
x ˙ = y , y ˙ = z , z ˙ = j ( x , y , z ) .
Jerk systems have triggered wide interest. On one hand, they have a simple form but display complex behavior, and on the other hand, jerk systems model, among other things, some oscillators, in particular chaotic oscillators (see, e.g., [2,3,4,5]). Jerk systems display different dynamic behaviors, which depend on the total number of terms and parameters of j, particularly on the number of nonlinearities (see, e.g., [6] and references therein). The variation in the parameters leads to some bifurcations and an asymptotically stable or periodic orbit can becomes chaotic. Bifurcations of particular jerk systems were analyzed (see, e.g., [7,8,9,10,11]). Moreover, new chaotic jerk systems were recently reported (see, e.g., [12,13,14]).
In this paper, we analyze the jerk system
x ˙ = y y ˙ = z z ˙ = j ( x , y , z , α ) ,
where j is smooth, and α is a real parameter under the assumption that there is α = α 0 such that E ( x 0 , 0 , 0 ) is an equilibrium point of (1), that is,
j ( x 0 , 0 , 0 , α 0 ) = 0 .
Of course, E can be an equilibrium point for every α R . For instance, if j is odd, i.e., j ( x , y , z , α ) = j ( x , y , z , α ) , for every x , y , z , α , then j ( 0 , 0 , 0 , α ) = 0 , for every α .
We denote
j x ( x 0 , 0 , 0 , α ) = a ( α ) , j y ( x 0 , 0 , 0 , α ) = b ( α ) , j z ( x 0 , 0 , 0 , α ) = c ( α ) ,
and
a ( α 0 ) = a , b ( α 0 ) = b , c ( α 0 ) = c .
The Jacobian matrix of system (1) at the equilibrium E is
J ( x 0 , 0 , 0 , α ) = 0 1 0 0 0 1 a ( α ) b ( α ) c ( α )
and the characteristic equation is given by
λ 3 c ( α ) λ 2 b ( α ) λ a ( α ) = 0 .
For α = α 0 , it becomes
λ 3 c λ 2 b λ a = 0 ,
which allows us to discuss the stability of E and some codim-1 local bifurcations, namely Hopf and fold bifurcations.
A fold bifurcation is a local bifurcation associated with the appearance of a zero eigenvalue. In the one-dimensional case, the generic fold bifurcation, also called the saddle-node bifurcation, has the normal form given by the equation
η ˙ = β ± η 2 ,
where β is a real parameter (see, e.g., [15]). The other common fold bifurcations are the transcritical bifurcation and the pitchfork bifurcation, and their normal forms are also given by one-dimensional equations. In higher dimensions, the same normal form for the saddle-node bifurcation is obtained by reduction on the center manifold. One of the goals of this work is to obtain (approximate) normal forms for these fold bifurcations that continue being jerk systems, that is, the simplest jerk systems that display such bifurcations. To obtain such normal forms, we first study the stability of an arbitrary jerk system, and we point out the framework that ensures the existence of the above-mentioned bifurcations.
The paper is organized as follows: In Section 2, we discuss the local stability of an equilibrium point E of system (1), taking into account the sign of a , b , c and a + b c . In Section 3, using again the sign of a , b , c and the behavior of the function h ( α ) = a ( α ) + b ( α ) c ( α ) , we prove that system (1) experiences a Hopf bifurcation. In Section 4, we establish conditions of existence for each above-mentioned fold bifurcation. In Section 5, we deduce approximate normal forms for the saddle-node bifurcation, the transcritical bifurcation, and the pitchfork bifurcation. These normal forms are also jerk systems. We study the stability of the proposed approximate normal form using the reduction on the center manifold. Also, we give bifurcation diagrams and some orbits of these normal forms. In Section 6, we present some conclusions.

2. Stability of a Jerk System

The equilibria of system (1) are of the form E ( x * , 0 , 0 ) , where x * fulfills the condition j ( x * , 0 , 0 , α 0 ) = 0 .
Theorem 1.
Consider jerk system (1) with property (2) and notations (3) and (4).
(a) 
For a 0 , the equilibrium point E ( x 0 , 0 , 0 )
(i) 
 is asymptotically stable if a < 0 , b < 0 , c < 0 , and b c + a > 0 ;
(ii) 
  is unstable if a > 0 , or b 0 , or c 0 , or b c + a < 0 .
(b) 
For a = 0 , the equilibrium point E ( x 0 , 0 , 0 ) is unstable if b > 0 or c > 0 .
Proof. 
Taking into account (6), the conclusions follow from the Routh–Hurwitz theorem (see, e.g., [16]). □
Remark 1.
For a 0 , the stability of the equilibrium E is discussed in all cases except a < 0 , b < 0 , c < 0 , b c + a = 0 . In this case, there is a pair of pure imaginary eigenvalues. Therefore, reduction on the center manifold will be necessary. In particular, if the first Lyapunov coefficient l 1 ( 0 ) is computed, then we can conclude that E is weakly asymptotically stable if l 1 ( 0 ) < 0 and unstable if l 1 ( 0 ) > 0 .
If a = 0 , the follow cases remain: b < 0 , c < 0 (a zero eigenvalue), b = 0 , c < 0 (a pair of zero eigenvalues), b < 0 , c = 0 (a zero eigenvalue and a pair of pure imaginary eigenvalues), and b = c = 0 (all eigenvalues are zero). In the first two cases, the stability of E results from the reduced form of the center manifold.
Some of the above-mentioned cases are related to some local codim-1 bifurcations that can occur in the dynamics of the considered system. Such bifurcations will be discussed in the next sections.

3. Hopf Bifurcation

In this section, we give conditions for which system (1) experiences a Hopf bifurcation.
Theorem 2.
Consider jerk system (1) with property (2) and notations (3) and (4), and let h ( α ) = a ( α ) + b ( α ) c ( α ) . Assume that E ( x 0 , 0 , 0 ) is an equilibrium point of system (1) for every α in the neighborhood of α 0 .
If  a < 0 , b < 0 ,   c < 0 ,   h ( α 0 ) = 0 ,  and  h ( α 0 ) 0 ,  then system (1displays a Hopf bifurcation at the point  E ( x 0 , 0 , 0 )  when α passes through the critical value  α 0 .
Proof. 
For h ( α 0 ) = 0 , that is, a = b c , the roots of (6) are λ 1 = c < 0 and λ 2 , 3 = ± i b .
It remains to check the transversality condition from Hopf’s bifurcation theorem (see, e.g., [17]). Using (5) and the implicit function theorem, we have
d λ d α = c ( α ) λ 2 + b ( α ) λ + a ( α ) 3 λ 2 2 c ( α ) λ b ( α )
and
R e d λ d α ( α 0 ) = a ( α 0 ) + b ( α 0 ) c + b c ( α 0 ) 2 ( b c 2 ) 0 ,
which finishes the proof. □
Remark 2.
In some particular cases, we can calculate the first Lyapunov coefficient l 1 ( 0 ) , which allows to see if the Hopf bifurcation is non-degenerate ( l 1 ( 0 ) 0 ) or degenerate ( l 1 ( 0 ) = 0 ) .
Remark 3.
Since the Hopf bifurcation involves a change in the stability, the conditions a < 0 ,   b < 0 , c < 0 are necessary conditions. In this case, the change in the stability is given by h ( α 0 ) = 0 and h ( α 0 ) 0 (the transversality condition), where h ( α ) = a ( α ) + b ( α ) c ( α ) (3).
Taking into account the above results, we notice that one of the simplest nonlinear jerk system which experiences the Hopf bifurcation at O ( 0 , 0 , 0 ) when the parameter α passes through the critical value α = 0 is given by
x ˙ = y y ˙ = z z ˙ = ( α b c ) x + b y + c z + y z ,
with b < 0 , c < 0 . In fact, it is easy to see that any jerk system (1) with
j ( x , y , z , α ) = ( α b c ) x + b y + c z + ψ ( x , y , z )
has the same property, where ψ is smooth, ψ ( 0 , 0 , 0 ) = 0 , d ψ ( 0 , 0 , 0 ) = 0 , b < 0 , c < 0 .

4. Fold Bifurcations

In this section, we study the bifurcation of the considered system associated with the appearance of a zero eigenvalue, namely the fold bifurcation or zero bifurcation. The non-degenerate fold bifurcation is called the saddle-node bifurcation. The most known degenerate fold bifurcations are the transcritical bifurcation and the pitchfork bifurcation.
Let α R be the bifurcation parameter such that at α = α 0 , system (1) has the non-hyperbolic equilibrium E ( x 0 , 0 , 0 ) with only one zero eigenvalue. Using (6), it follows that a = 0 and b 0 .
In the following, we point out fold bifurcations of system (1) using Sotomayor’s theorem [18] (also see [17,19]). First, we compute the expressions that are used in Sotomayor’s theorem.
Proposition 1.
Let  F ( x , y , z , α ) = ( y , z , j ( x , y , z , α ) )  given by jerk system (1with property (2and notations (3and (4). If a = 0 , then the following assertions hold:
1. 
D F x 0 , 0 , 0 , α 0 has a simple eigenvalue λ 1 = 0 with right eigenvector v = ( 1 , 0 , 0 ) t and left eigenvector w = b , c , 1 t ;
2. 
F α x 0 , 0 , 0 , α 0 = 0 , 0 , j α ( x 0 , 0 , 0 , α 0 ) t and w t F α x 0 , 0 , 0 , α 0 = j α ( x 0 , 0 , 0 , α 0 ) ;
3. 
D F α x 0 , 0 , 0 , α 0 ( v ) = 0 , 0 , 2 j x α ( x 0 , 0 , 0 , α 0 ) t and
w t D F α x 0 , 0 , 0 , α 0 ( v ) = 2 j x α ( x 0 , 0 , 0 , α 0 ) ;
4. 
D 2 F x 0 , 0 , 0 , α 0 ( v , v ) = 0 , 0 , 2 j x 2 ( x 0 , 0 , 0 , α 0 ) t and
w t D 2 F x 0 , 0 , 0 , α 0 ( v , v ) = 2 j x 2 ( x 0 , 0 , 0 , α 0 ) ;
5. 
D 3 F x 0 , 0 , 0 , α 0 ( v , v , v ) = 0 , 0 , 3 j x 3 ( x 0 , 0 , 0 , α 0 ) t and
w t D 3 F x 0 , 0 , 0 , α 0 ( v , v , v ) = 3 j x 3 ( x 0 , 0 , 0 , α 0 ) .
Now, using Proposition 1 and Sotomayor’s theorem, we obtain the next results.
Theorem 3.
Let  α R  be the bifurcation parameter of system (1with the critical value α 0 R If
SN1. 
j ( x 0 , 0 , 0 , α 0 ) = 0 , j x ( x 0 , 0 , 0 , α 0 ) = a = 0 , j y ( x 0 , 0 , 0 , α 0 ) = b 0 ,
j z ( x 0 , 0 , 0 , α 0 ) = c 0 ,
SN2. 
j α ( x 0 , 0 , 0 , α 0 ) 0 ,
SN3. 
2 j x 2 ( x 0 , 0 , 0 , α 0 ) 0 ,
then system (1displays a saddle-node bifurcation at the point  E ( x 0 , 0 , 0 ) .
More precisely, there are two hyperbolic equilibrium points of the system near E for α < α 0 (or α > α 0 ), which collide for α = α 0 into the non-hyperbolic equilibrium E and then disappear. In addition, the two equilibrium points have stable manifolds of dimensions k and k + 1 respectively, where k is the number of eigenvalues of D F ( x 0 , 0 , 0 , α 0 ) with a negative real part.
Remark 4.
If we consider b < 0 , c < 0 in the above theorem, then the Jacobian matrix D F ( x 0 , 0 , 0 , α 0 ) has k = 2 eigenvalues with negative real part. Therefore, the two equilibrium points that collide are a saddle and a stable node respectively (as in the one-dimensional case of the saddle-node bifurcation). If b < 0 , c > 0 , then k = 0 and the points are an unstable node and a saddle. Finally, if b > 0 , then k = 1 and both equilibrium points are saddles.
Theorem 4.
Let  α R  be the bifurcation parameter of system (1with the critical value α 0 R If
T1. 
j ( x 0 , 0 , 0 , α 0 ) = 0 , j x ( x 0 , 0 , 0 , α 0 ) = a = 0 , j y ( x 0 , 0 , 0 , α 0 ) = b 0 , j z ( x 0 , 0 , 0 , α 0 ) = c 0 ,
T2. 
j α ( x 0 , 0 , 0 , α 0 ) = 0 , 2 j x α ( x 0 , 0 , 0 , α 0 ) 0 ,
T3. 
2 j x 2 ( x 0 , 0 , 0 , α 0 ) 0 ,
then system (1displays a transcritical bifurcation at the point  E ( x 0 , 0 , 0 ) .
Theorem 5.
Let  α R  be the bifurcation parameter of system (1with the critical value  α 0 R If
P1. 
j ( x 0 , 0 , 0 , α 0 ) = 0 , j x ( x 0 , 0 , 0 , α 0 ) = a = 0 , j y ( x 0 , 0 , 0 , α 0 ) = b 0 , j z ( x 0 , 0 , 0 , α 0 ) = c 0 ,
P2. 
j α ( x 0 , 0 , 0 , α 0 ) = 0 , 2 j x α ( x 0 , 0 , 0 , α 0 ) 0 ,
P3. 
2 j x 2 ( x 0 , 0 , 0 , α 0 ) = 0 , 3 j x 3 ( x 0 , 0 , 0 , α 0 ) 0 ,
then system (1displays a pitchfork bifurcation at the point  E ( x 0 , 0 , 0 ) .
Remark 5.
The condition S N 1 ( T 1 , P 1 ) implies that E ( x 0 , 0 , 0 ) is an equilibrium point of the considered jerk system with a simple zero eigenvalue and no other eigenvalues on the imaginary axis. The genericity conditions of the saddle-node bifurcation are the transversality condition S N 2 and the nondegeneracy condition S N 3 . Transcritical and pitchfork bifurcations are obtained violating some of these genericity conditions.

5. Approximate Normal Forms for Fold Bifurcations

In this section, we derive approximate normal forms for each of the above fold bifurcations. We mention that these normal forms are also jerk systems.
In the following, we consider x 0 = 0 and α 0 = 0 . Then, for the fold bifurcation we have
j ( 0 , 0 , 0 , 0 ) = 0 ,   a = j x ( 0 , 0 , 0 , 0 ) = 0 ,   b = j y ( 0 , 0 , 0 , 0 ) 0 ,   c = j z ( 0 , 0 , 0 , 0 ) 0 .

5.1. Approximate Normal Forms for the Saddle-Node Bifurcation

We follow the steps used for the one-dimensional case (see, e.g., [15]). Assume that conditions S N 2 and S N 3 from Theorem 3 are fulfilled.
Consider a Taylor expansion of the function j with respect to ( x , y , z ) at ( 0 , 0 , 0 ) ,
j ( x , y , z ) = j 0 ( α ) + j 1 ( α ) x + 1 2 j 2 ( α ) x 2 + b ( α ) y + c ( α ) z + B ( α ) x y + C ( α ) x z + O ,
where O is the remainder of the expansion, and
j 0 ( α ) = j ( 0 , 0 , 0 , α ) , j 1 ( α ) = j x ( 0 , 0 , 0 , α ) , j 2 ( α ) = 2 j x 2 ( 0 , 0 , 0 , α ) , b ( α ) = j y ( 0 , 0 , 0 , α ) , c ( α ) = j z ( 0 , 0 , 0 , α ) , B ( α ) = 2 j x y ( 0 , 0 , 0 , α ) , C ( α ) = 2 j x z ( 0 , 0 , 0 , α ) .
We do not explicitly use the term y z in the above Taylor expansion because such a term and its derivatives vanish at ( 0 , 0 , 0 , α ) .
Step 1. Using the translation x = ε δ , y = y , z = z , where δ = δ ( α ) , system (1) is written as
ε ˙ = y y ˙ = z z ˙ = j 0 ( α ) j 1 ( α ) δ + 1 2 j 2 ( α ) δ 2 + O ( δ 3 ) + j 1 ( α ) j 2 ( α ) δ + O ( δ 2 ) ε + 1 2 j 2 ( α ) + O ( δ ) ε 2 + ( b ( α ) B ( α ) δ ) y + ( c ( α ) C ( α ) δ ) z + B ( α ) ε y + C ( α ) ε z + O .
We find δ such that the coefficient of ε vanishes. Let
F ( α , δ ) = j 1 ( α ) j 2 ( α ) δ + ψ ( α , δ ) δ 2 .
Then,
F ( 0 , 0 ) = j 1 ( 0 ) = j x ( 0 , 0 , 0 , 0 ) = 0 ( S N 1 ) ,
and
F δ ( 0 , 0 ) = j 2 ( 0 ) = 2 j x 2 ( 0 , 0 , 0 , 0 ) 0 ( S N 3 ) .
Therefore, by using the implicit function theorem, there is a function δ = δ ( α ) with δ ( 0 ) = 0 and F ( α , δ ( α ) ) = 0 in the neighborhood of α = 0 . Moreover, since
F α ( 0 , 0 ) = j 1 α ( 0 ) = 2 j x α ( 0 , 0 , 0 , 0 ) : = a ( 0 ) ,
we obtain
δ = a ( 0 ) j 2 ( 0 ) α + O ( α 2 ) .
Using δ and the following Taylor expansions
j 0 ( α ) = j α ( 0 , 0 , 0 , 0 ) α + O ( α 2 ) = j 0 ( 0 ) α + O ( α 2 ) , j 1 ( α ) = 2 j x α ( 0 , 0 , 0 , 0 ) α + O ( α 2 ) , j 2 ( α ) = j 2 ( 0 ) + O ( α 2 ) , b ( α ) = b + O ( α ) , c ( α ) = c + O ( α ) , B ( α ) = B ( 0 ) + O ( α ) , C ( α ) = C ( 0 ) + O ( α ) ,
system (7) becomes
ε ˙ = y y ˙ = z z ˙ = j 0 ( 0 ) α + O ( α 2 ) + 1 2 j 2 ( 0 ) + O ( α ) ε 2 + ( b + O ( α ) ) y + ( c + O ( α ) ) z + ( B ( 0 ) + O ( α ) ) ε y + ( C ( 0 ) + O ( α ) ) ε z + O .
Step 2. We denote
j 0 ( 0 ) α + ψ ( α , δ ) α 2 = μ .
Let
G ( μ , α ) = j 0 ( 0 ) α + ψ ( α , δ ) α 2 μ
We have G ( 0 , 0 ) = 0 and
G α ( 0 , 0 ) = j 0 ( 0 ) = j α ( 0 , 0 , 0 , 0 ) 0 ( S N 2 ) .
Then, there is a function α = α ( μ ) such that α ( 0 ) = 0 and F ( μ , α ( μ ) ) = 0 in the neighborhood of μ = 0 . Consequently, system (8) becomes
ε ˙ = y y ˙ = z z ˙ = μ + 1 2 j 2 ( 0 ) + O ( μ ) ε 2 + ( b + O ( μ ) ) y + ( c + O ( μ ) ) z + ( B ( 0 ) + O ( μ ) ) ε y + ( C ( 0 ) + O ( μ ) ) ε z + O .
Step 3. We denote
A ( μ ) = 1 2 j 2 ( 0 ) + O ( μ ) .
It follows that
A ( 0 ) = 1 2 j 2 ( 0 ) = 1 2 2 j x 2 ( 0 , 0 , 0 , 0 ) 0 ( S N 3 ) ,
thus, A ( μ ) 0 in the neighborhood V of μ = 0 . Moreover,
s i g n ( A ( μ ) ) : = s = sign ( 2 j x 2 ( 0 , 0 , 0 , 0 ) )
on V .
Using the change in the parameter and variables given by
β = s A ( μ ) μ , X = s A ( μ ) ε , Y = s A ( μ ) y , Z = s A ( μ ) z ,
with μ = 2 s j 2 ( 0 ) β + O ( β 2 ) , system (9) becomes
X ˙ = Y Y ˙ = Z Z ˙ = β + s X 2 + ( b + O ( β ) ) Y + ( c + O ( β ) ) Z + ( 2 s B ( 0 ) j 2 ( 0 ) + O ( β ) ) X Y + ( 2 s C ( 0 ) j 2 ( 0 ) + O ( β ) ) X Z + O ,
where the remainder O was properly changed, but it does not contain terms in the form of those already written in the third equation of (10).
In conclusion, we have proven the following result.
Theorem 6.
Consider jerk system (1with smooth j such that  O ( 0 , 0 , 0 )  is an equilibrium at  α = 0  with a zero eigenvalue, that is,  j ( 0 , 0 , 0 , 0 ) = 0 j x ( 0 , 0 , 0 , 0 ) = 0 .  Suppose that
SN1. 
b = j y ( 0 , 0 , 0 , 0 ) 0 , c = j z ( 0 , 0 , 0 , 0 ) 0 ,
SN2. 
j α ( 0 , 0 , 0 , 0 ) 0 ,
SN3. 
2 j x 2 ( 0 , 0 , 0 , 0 ) 0 .
Then, there are invertible coordinate and parameter changes transforming system  (1)  into system (10).
Let s = 1 . An equilibrium point of system (10) has the form E ( x 0 , 0 , 0 ) , where x 0 is a solution of the equation β + X 2 + O ( X 3 ) = 0 . Since the fold bifurcation takes place at β = 0 in O, we are concerned with equilibria close to O, thus x 0 is near 0. Recall that for any sufficiently small β < 0 , the above equation has two solutions near the origin, which are close to the solutions of the equation β + x 2 = 0 for the same parameter value (see, e.g., the proof of Lemma 3.1 [15]). Therefore, in order to obtain a normal form for the fold bifurcation in O ( 0 , 0 , 0 ) at β = 0 , we can truncate system (10). Thus, we obtain the following truncated normal form for the generic fold bifurcation of a jerk system:
X ˙ = Y Y ˙ = Z Z ˙ = β + s X 2 + b Y + c Z + 2 s B ( 0 ) j 2 ( 0 ) X Y + 2 s C ( 0 ) j 2 ( 0 ) X Z .
Moreover, we can consider the following approximate normal form for the generic fold bifurcation of a jerk system:
x ˙ = y y ˙ = z z ˙ = α + s x 2 + b y + c z ,
where α R is the parameter, s = ± 1 , and b , c R are fixed such that b , c 0 .
We note that the last two systems have the same equilibria and the corresponding characteristic polynomials are similar.
We will study system (12) near the origin and α in the neighborhood of 0 .
Let s = 1 . Then, system (12) has two equilibria E ± ( ± α , 0 , 0 ) for α < 0 ; one equilibrium point O ( 0 , 0 , 0 ) for α = 0 ; and no equilibria for α > 0 .
Proposition 2.
Let O ( 0 , 0 , 0 ) and E ± ( ± α , 0 , 0 ) be the equilibria of system (12), where α 0 such that | α | is sufficiently small.
1. 
Assume b < 0 and c < 0 .
(a) 
If α < 0 , then E is asymptotically stable and E + is unstable.
(b) 
If α = 0 , then O = E is unstable.
2. 
If b > 0 or c > 0 , then O and E ± are unstable for any α , provided they exist.
Proof. 
According to Theorem 1, it only remains to study the case α = 0 with b < 0 , c < 0 . Because the equilibrium point has a zero eigenvalue, we will reduce system (12) on the center manifold (see, e.g., [17,20]). Using the eigenvectors of the Jacobian matrix at O, we deduce the following transformation
u = c + c 2 + 4 b 2 b c 2 + 4 b c + c 2 + 4 b 2 y + z v = c c 2 + 4 b 2 b c 2 + 4 b c c 2 + 4 b 2 y + z w = x + c b y 1 b z .
Then, system (12) becomes
u ˙ = c c 2 + 4 b 2 u + c + c 2 + 4 b 2 b c 2 + 4 b ( u + v + w ) 2 v ˙ = c + c 2 + 4 b 2 v c c 2 + 4 b 2 b c 2 + 4 b ( u + v + w ) 2 w ˙ = 1 b ( u + v + w ) 2 .
Consider the center manifold
W loc O = ( u , v , w ) | u = F ( w ) , v = G ( w ) , | w | 1 , F ( 0 ) = F ( 0 ) = G ( 0 ) = G ( 0 ) = 0 ,
where F ( w ) = A w 2 + O ( w 3 ) , G ( w ) = B w 2 + O ( w 3 ) . Using (13), we obtain
A = ( c + c 2 + 4 b ) 2 4 b 2 c 2 + 4 b , B = ( c c 2 + 4 b ) 2 4 b 2 c 2 + 4 b ,
and the reduced equation on W l o c O , given by
w ˙ = 1 b w 2 + O ( w 2 ) .
Thus, O is unstable, and the proof is finished. □
For s = 1 , we obtain a similar result.
Remark 6.
We notice that the results presented in Remark 4 are valid for the approximate normal form (12).
In Proposition 2, we discuss the stability of system (12) for α close to zero. The equilibrium points can be stable only if b , c < 0 . In this case, when α < 0 moves away from zero we deduce that E + remains unstable and E changes its stability for α = b 2 c 2 4 (see Theorem 1). Moreover, a Hopf bifurcation occurs in the dynamics of system (12) (Theorem 2), and a stable limit cycle is born (see Figure 1c–e). The saddle-node and Hopf bifurcations are also pointed out in the bifurcation diagram obtained by numerical simulations using MatCont/MatLab (Figure 1a).

5.2. Approximate Normal Forms for the Transcritical Bifurcation

The transcritical bifurcation is a non-generic fold bifurcation, and the method used in the above paragraph is useless.
Assume that hypotheses T 1 T 3 of Theorem 4 are satisfied. Therefore, we have the following Taylor expansion of the function j with respect to ( x , y , z , α ) at ( 0 , 0 , 0 , 0 )
j ( x , y , z , α ) = b y + c z + a 11 x 2 + a 12 x y + a 13 x z + a 14 α x + O ,
where b , c , a 11 , a 14 0 . Let s = s i g n ( a 11 ) . Using the transformation
X = s · a 11 x , Y = s · a 11 y , Z = s · a 11 z , β = a 14 α ,
system (1) becomes
X ˙ = Y Y ˙ = Z Z ˙ = β X + s X 2 + b Y + c Z + s a 12 a 11 X Y + s a 13 a 11 X Z + O .
An equilibrium point of system (14) has the form E ( x 0 , 0 , 0 ) , where x 0 is a solution of the equation β X + s X 2 + O ( X 3 ) = 0 . A solution is x 1 = 0 , and the others satisfy F ( β , X ) = 0 near ( 0 , 0 ) , where F ( β , X ) = β + s X + X 2 φ ( X , β ) . From the implicit function theorem, there is a unique function x = x ( β ) with x ( 0 ) = 0 and F ( β , x ( β ) ) = 0 in the neighborhood of β = 0 . Moreover, x ( 0 ) = s , thus, x ( β ) = s β + O ( β 2 ) near 0, that is, for any sufficiently small β , the equation F ( β , X ) = 0 has a unique solution X = X ( β ) near the origin, which is close to the solution of the equation β + s x = 0 for the same parameter value. In conclusion, near the origin, the equilibria of system (14) are close to the equilibria of its truncated normal forms
X ˙ = Y Y ˙ = Z Z ˙ = β X + s X 2 + b Y + c Z + s a 12 a 11 X Y + s a 13 a 11 X Z
and
x ˙ = y y ˙ = z z ˙ = α x + s x 2 + b y + c z ,
where α R is the parameter, s = ± 1 , and b , c R are fixed such that b , c 0 .
The characteristic polynomials at the equilibria of the last two systems are similar. Therefore, we consider that system (16) is an approximate normal form for the transcritical bifurcation of a jerk system.
Let s = 1 . Then, system (16) has two equilibria for every α 0 , namely, O ( 0 , 0 , 0 ) and E ( α , 0 , 0 ) , which collide for α = 0 .
Proposition 3.
Let O and E be the equilibria of system (16) and α R with sufficiently small | α | .
1. 
Assume b < 0 and c < 0 .
(a) 
If α < 0 , then O is asymptotically stable and E is unstable.
(b) 
If α = 0 , then O = E is unstable.
(c) 
If α > 0 , then E is asymptotically stable and O is unstable.
2. 
If b > 0 or c > 0 , then O and E are unstable for any α .
Proof. 
We immediately obtain the results 1 ( a ) , 1 ( c ) , and 2 via Theorem 1. For α = 0 ( b < 0 , c < 0 ), we proceed as in the proof of Proposition 2. In this case, the reduced equation on the center manifold is given by
w ˙ = 1 b w 2 + O ( w 3 ) .
Thus, O is unstable, and the proof is finished. □
For s = 1 , we obtain a similar result.
Remark 7.
We notice that the case b < 0 , c < 0 explains the transcritical bifurcation, also named the exchange of stability: two equilibria collide when the parameter passes through the critical value and then split, exchanging their stability. In the other cases, even if both equilibria are unstable, there is an exchange of the dimensions of the stable manifolds at the bifurcation value.
The bifurcation diagram of system (16) is presented in Figure 2. Again, we notice a Hopf bifurcation, which follows from Theorem 2. In this case, the orbits look similar to those depicted in Figure 1.

5.3. Approximate Normal Forms for the Pitchfork Bifurcation

In this case, we work under the hypothesis P 1 P 3 of Theorem 5. It results in the following Taylor expansion of the function j with respect to ( x , y , z , α ) at ( 0 , 0 , 0 , 0 )
j ( x , y , z , α ) = b y + c z + a 11 x 3 + a 12 x y + a 13 x z + a 14 α x + O ,
where b , c , a 11 , a 14 0 . Using the transformation
X = p x , Y = p y , Z = p z , β = a 14 α ,
where p = s a 11 , s = s i g n ( a 11 ) , system (1) writes
X ˙ = Y Y ˙ = Z Z ˙ = β X + s X 3 + b Y + c Z + a 12 p X Y + a 13 p X Z + O .
An equilibrium point of system (17) has the form E ( x 0 , 0 , 0 ) , where x 0 is a solution of the equation β X + s X 3 + O ( X 4 ) = 0 . A solution is x 1 = 0 , and the others satisfy β + s X 2 + O ( X 3 ) = 0 near ( 0 , 0 ) , which is the same equation as in the case of the saddle-node bifurcation. Thus, we can truncate the above expansion. Consequently, we obtain the following approximate normal forms for the pitchfork bifurcation of a jerk system:
X ˙ = Y Y ˙ = Z Z ˙ = β X + s X 3 + b Y + c Z + a 12 p X Y + a 13 p X Z
and
x ˙ = y y ˙ = z z ˙ = α x + s x 3 + b y + c z ,
where α R is the parameter, s = ± 1 , and b , c R are fixed such that b , c 0 .
Proposition 4.
Let s = 1 . Then, O ( 0 , 0 , 0 ) is an equilibrium point of system (19for every α R , and other equilibria are E ± ( ± α , 0 , 0 ) if α > 0 .
Let α R with sufficiently small | α | .
1. 
Assume b < 0 and c < 0 .
(a) 
If α < 0 , then O is asymptotically stable.
(b) 
If α = 0 , then O = E + = E is stable.
(c) 
If α > 0 , then E ± are asymptotically stable and O is unstable.
2. 
If b > 0 or c > 0 , then O and E ± are unstable for any α, provided they exist.
Proof. 
For α = 0 ( b < 0 , c < 0 ), we proceed as in the proof of Proposition 2. In this case, the reduced equation of the center manifold is given by
w ˙ = 1 b w 3 + O ( w 4 ) ,
thus O is stable. The other conclusions result from Theorem 1. □
Similarly, we obtain:
Proposition 5.
Let s = 1 . Then, O ( 0 , 0 , 0 ) is an equilibrium point of system (19for every α R , and other equilibria are E ± ( ± α , 0 , 0 ) if α < 0 .
Let α R with sufficiently small | α | .
1. 
Assume b < 0 and c < 0 .
(a) 
If α < 0 , then O is asymptotically stable and E ± are unstable.
(b) 
If α = 0 , then O = E + = E is unstable.
(c) 
If α > 0 , then O is unstable.
2. 
If b > 0 or c > 0 , then O and E ± are unstable for any α, provided they exist.
Proof. 
For α = 0 , b < 0 , c < 0 the reduced equation on the center manifold is given by
w ˙ = 1 b w 3 + O ( w 4 ) ,
thus O is unstable. □
Remark 8.
The pitchfork bifurcation shows the following local phenomenon in the dynamics of a system: at the point of bifurcation, the stability of an equilibrium changes, and a pair of equilibria appears (see Figure 3) or disappears (see Figure 4). This situation is pointed out in the case b < 0 , c < 0 . More precisely, for s = 1 we have obtained the supercritical pitchfork bifurcation, that is, a stable equilibrium point becomes unstable and a pair of stable equilibria is born from it at the branch point BP (Figure 3). For s = 1 , a stable equilibrium point collides with a pair of unstable equilibria and becomes unstable, which indicates a subcritical pitchfork bifurcation (Figure 4).
The pitchfork bifurcation is specific to systems with symmetries, more precisely to those systems that are invariant under the transformation ( x , y , z ) ( x , y , z ) . We notice that system (19) is such a system. Moreover, such systems can have symmetrical orbits and symmetrical chaotic attractors (see, e.g., [2,10]), as shown in Figure 5.

6. Conclusions

In this paper, we have studied the local stability and codim-1 bifurcations of an arbitrary one-parameter jerk system. We have established conditions for which Hopf or fold bifurcations occur in the dynamics of a jerk system. It is known that for an n dimensional dynamical system, a normal form of a local bifurcation is obtained by reduction on the center manifold. The normal forms of the most well-known fold bifurcations, namely the saddle-node bifurcation, the transcritical bifurcation, and the pitchfork bifurcation, are given by one-dimensional equations. Our goal was to derive three-dimensional (approximate) normal forms, which are themselves jerk systems. We have obtained such a normal form for each above-mentioned fold bifurcation. Particularly, because the saddle-node bifurcation is a generic fold bifurcation, we have derived such an approximate normal form using invertible coordinate and parameter changes.

Author Contributions

Conceptualization, C.L.; methodology, C.L.; formal analysis, C.L. and J.C.; writing—original draft preparation, C.L. and J.C.; writing—review and editing, C.L.; visualization, C.L. and J.C.; supervision, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the referees very much for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schot, S.H. Jerk: The time rate of change of acceleration. Am. J. Phys. 1978, 46, 1090–1094. [Google Scholar] [CrossRef]
  2. Coullet, P.H.; Tresser, C.; Arneodo, A. Transition to stochasticity for a class of forced oscillators. Phys. Lett. A 1979, 72, 268–270. [Google Scholar] [CrossRef]
  3. Gottlieb, H.P.W. Question 38. What is the simplest jerk function that gives chaos? Am. J. Phys. 1996, 64, 525. [Google Scholar] [CrossRef]
  4. Sprott, J.C. Some simple chaotic jerk functions. Am. J. Phys. 1997, 65, 537–543. [Google Scholar] [CrossRef]
  5. Wei, Z.; Sprott, J.C.; Chen, H. Elementary quadratic chaotic flows with a single non-hyperbolic equilibrium. Phys. Lett. A 2015, 379, 2184–2187. [Google Scholar] [CrossRef]
  6. Singh, P.P.; Roy, B.K. Pliers shaped coexisting bifurcation behaviors in a simple jerk chaotic system in comparison with 21 reported systems. IFAC PapersOnLine 2022, 55, 920–926. [Google Scholar] [CrossRef]
  7. Tang, D.; Zhang, S.; Ren, J. Dynamics of a general jerky equation. J. Vib. Control 2019, 25, 922–932. [Google Scholar] [CrossRef]
  8. Sang, B.; Huang, B. Zero-Hopf Bifurcations of 3D Quadratic Jerk System. Mathematics 2020, 8, 1454. [Google Scholar] [CrossRef]
  9. Braun, F.; Mereu, A.C. Zero-Hopf bifurcation in a 3D jerk system. Nonlinear Anal. Real World Appl. 2021, 59, 103245. [Google Scholar] [CrossRef]
  10. Lăzureanu, C. Dynamical Properties, Deformations, and Chaos in a Class of Inversion Invariant Jerk Equations. Symmetry 2022, 14, 1318. [Google Scholar] [CrossRef]
  11. Lyu, W.; Li, S.; Chen, Z.; Bi, Q. Bursting Dynamics in a Singular Vector Field with Codimension Three Triple Zero Bifurcation. Mathematics 2023, 11, 2486. [Google Scholar] [CrossRef]
  12. Zhang, S.; Zeng, Y. A simple Jerk-like system without equilibrium: Asymmetric coexisting hidden attractors, bursting oscillation and double full Feigenbaum remerging trees. Chaos Solitons Fractals 2019, 120, 25–40. [Google Scholar] [CrossRef]
  13. Li, F.; Zeng, J. Multi-Scroll Attractor and Multi-Stable Dynamics of a Three-Dimensional Jerk System. Energies 2023, 16, 2494. [Google Scholar] [CrossRef]
  14. Vaidyanathan, S.; Kammogne, A.S.T.; Tlelo-Cuautle, E.; Talonang, C.N.; Abd-El-Atty, B.; Abd El-Latif, A.A.; Kengne, E.M.; Mawamba, V.F.; Sambas, A.; Darwin, P.; et al. A Novel 3-D Jerk System, Its Bifurcation Analysis, Electronic Circuit Design and a Cryptographic Application. Electronics 2023, 12, 2818. [Google Scholar] [CrossRef]
  15. Kuznetsov, Y.A. Elements of Applied Bifurcation Theory; Springer: New York, NY, USA, 1998. [Google Scholar]
  16. Pontryagin, L.S. Ordinary Differential Equations; Addison-Wesley Publishing Company Inc.: Reading, MA, USA, 1962. [Google Scholar]
  17. Perko, L. Differential Equations and Dynamical Systems, 3rd ed.; Texts in Applied Mathematics 7; Springer: New York, NY, USA, 2001. [Google Scholar]
  18. Sotomayor, J. Generic Bifurcations of Dynamical Systems. In Dynamical Systems; Peixoto, M.M., Ed.; Academic Press: New York, NY, USA, 1973; pp. 561–582. [Google Scholar]
  19. Guckenheimer, J.; Holmes, P. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Applied Mathematical Sciences; Springer: Berlin/Heidelberg, Germany, 1983. [Google Scholar]
  20. Osipenko, G. Center Manifolds. In Encyclopedia of Complexity and Systems Science; Meyers, R., Ed.; Springer: New York, NY, USA, 2009. [Google Scholar] [CrossRef]
Figure 1. The normalform (12), s = 1 , b = 1 , c = 1 : (a) the bifurcation diagram (H stands for Hopf bifurcation and LP for the limit point, that is, the fold bifurcation); (b) an asymptotically stable orbit ( a = 0.2 ); (c) a stable limit cycle ( a = 0.25 ); (d,e) unstable orbits attracted by the limit cycle ( a = 0.28 and a = 0.25 , respectively); (f) an unbounded unstable orbit ( a = 0 ). The initial conditions are close to E in (d) and close to O in the other cases.
Figure 1. The normalform (12), s = 1 , b = 1 , c = 1 : (a) the bifurcation diagram (H stands for Hopf bifurcation and LP for the limit point, that is, the fold bifurcation); (b) an asymptotically stable orbit ( a = 0.2 ); (c) a stable limit cycle ( a = 0.25 ); (d,e) unstable orbits attracted by the limit cycle ( a = 0.28 and a = 0.25 , respectively); (f) an unbounded unstable orbit ( a = 0 ). The initial conditions are close to E in (d) and close to O in the other cases.
Mathematics 11 04295 g001aMathematics 11 04295 g001b
Figure 2. The bifurcation diagram of the normal form (16), s = 1 , b = 1 , c = −1 (BP stands for the branch point, that is, it highlights the transcritical bifurcation in this case).
Figure 2. The bifurcation diagram of the normal form (16), s = 1 , b = 1 , c = −1 (BP stands for the branch point, that is, it highlights the transcritical bifurcation in this case).
Mathematics 11 04295 g002
Figure 3. The bifurcation diagram of the normal form (19), s = 1 , b = 1 , c = −1.
Figure 3. The bifurcation diagram of the normal form (19), s = 1 , b = 1 , c = −1.
Mathematics 11 04295 g003
Figure 4. The bifurcation diagram of the normal form (19), s = 1 , b = 1 , c = −1.
Figure 4. The bifurcation diagram of the normal form (19), s = 1 , b = 1 , c = −1.
Mathematics 11 04295 g004
Figure 5. Chaotic attractors in the dynamics of system (19) for s = 1 , b = 1 , c = 1 : (a) α = 1.37 ; (b) α = 1.47 .
Figure 5. Chaotic attractors in the dynamics of system (19) for s = 1 , b = 1 , c = 1 : (a) α = 1.37 ; (b) α = 1.47 .
Mathematics 11 04295 g005
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

Lăzureanu, C.; Cho, J. On Hopf and Fold Bifurcations of Jerk Systems. Mathematics 2023, 11, 4295. https://doi.org/10.3390/math11204295

AMA Style

Lăzureanu C, Cho J. On Hopf and Fold Bifurcations of Jerk Systems. Mathematics. 2023; 11(20):4295. https://doi.org/10.3390/math11204295

Chicago/Turabian Style

Lăzureanu, Cristian, and Jinyoung Cho. 2023. "On Hopf and Fold Bifurcations of Jerk Systems" Mathematics 11, no. 20: 4295. https://doi.org/10.3390/math11204295

APA Style

Lăzureanu, C., & Cho, J. (2023). On Hopf and Fold Bifurcations of Jerk Systems. Mathematics, 11(20), 4295. https://doi.org/10.3390/math11204295

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