On a Family of Hamilton–Poisson Jerk Systems

: In this paper, we construct a family of Hamilton–Poisson jerk systems. We show that such a system has infinitely many Hamilton–Poisson realizations. In addition, we discuss the stability and we prove the existence of periodic orbits around nonlinearly stable equilibrium points. Particularly, we deduce conditions for the existence of homoclinic and heteroclinic orbits. We apply the obtained results to a family of anharmonic oscillators


Introduction
Jerk is the rate of change of acceleration, the third derivative of position with respect to time [1].A jerk equation ... x = j(x, ẋ, ẍ) and the corresponding jerk system, which is a three-dimensional system given by can model processes characterized by changes in acceleration.Despite their simple form, jerk systems provide examples of chaotic behavior (see, e.g., [2][3][4]).Bifurcations in the dynamics of jerk systems are also analyzed (see, e.g., [5][6][7][8]).
In this paper, we study how Hamilton-Poisson jerk systems can be constructed.Roughly speaking, a three-dimensional system is a Hamilton-Poisson system if it has two independent constants of motion (for details on Hamilton-Poisson mechanics, see e.g., [9]).Using such functions, we obtain a family of Hamilton-Poisson jerk systems, given by ... x + f ′ (x) ẋ = 0, which are in fact jerk versions of the system with one degree of freedom ẍ + f (x) = 0 [10].For instance, the equations of the harmonic oscillator, the mathematical pendulum, the Duffing oscillator, and other anharmonic oscillators are of this form.
Oscillatory systems, characterized by repetitive patterns or cycles, are found in various biological phenomena such as circadian rhythms (see, e.g., [11]), neuronal activity (see, e.g., [12]), and even in cellular processes like metabolic oscillations (see, e.g., [13]).Population dynamics in predator-prey relationships often display cyclic behavior, where changes in predator and prey populations exhibit periodic patterns (see, e.g., [14]).Modeling changes in population sizes or ecological systems often involves sudden shifts or rapid changes in growth rates, which can be compared analogously to jerk-like behavior in dynamic systems.Moreover, in neural systems, sudden changes in firing rates or neuronal activities might indirectly relate to rapid changes in behavior akin to jerk-like dynamics.
The paper is organized as follows: in Section 2, we recall some notions regarding Hamilton-Poisson systems and then we give some conditions for which system (1) is of this type.Using the integrable deformation method (see [15] and references therein), we construct a family of Hamilton-Poisson jerk systems.Also, we give Hamilton-Poisson realizations of such a system.In Section 3, we analyze some dynamical properties of the obtained system, namely, the stability of the equilibrium points, the existence of the periodic orbits around some nonlinearly stable equilibria, and the existence of homoclinic or heteroclinic orbits.In Section 4, we apply these results to a family of anharmonic oscillators [16].

A Family of Hamilton-Poisson Jerk Systems
In this section, we construct a family of jerk systems that have Hamilton-Poisson realizations.
Recall that the three-dimensional dynamical system is a Hamilton-Poisson system on R 3 if there are the smooth functions ν, H, C such that on R 3 (see, e.g., [17,18]).The function ν is called the rescaling function.In addition, the Hamiltonian function H and the Casimir function C are constants of motion of the above system.In fact, a Hamilton-Poisson system on R 3 is a triple (R 3 , Π, H), where Π is a Poisson structure, and in this case it is given by the matrix where we have denoted C x = ∂C ∂x .Such a system writes ( ẋ, ẏ, ż) t = Π • ∇H.Details on Hamiltonian mechanics can be found, for example, in [19].
In the following, we consider ν = 1.One of our goals is to obtain jerk systems that can be written in the form ( ẋ, ẏ, ż) T = ∇H × ∇C, that is, to determine functions H, C such that and which are constants of motion of system (1), that is, We note that achieving this goal appears to be complicated for a general jerk function j.However, the next result holds.Theorem 1.If j z ̸ = 0, then jerk system (1) cannot have a Hamilton-Poisson formulation ( ẋ, ẏ, ż) T = ∇H × ∇C.
Instead of starting with a function j and checking for the existence of the functions H and C, we can construct Hamilton-Poisson jerk systems using integrable deformation method [15].
Consider the jerk equation ... x = 0 and the corresponding jerk system It is easy to see that the functions are constants of motion for system (3).Moreover, system (3) writes ( ẋ, ẏ, ż) T = ∇H × ∇C; thus, it has the Hamilton-Poisson realization (R 3 , Π, H) with the Hamiltonian H and the Poisson structure given by the matrix Now, we alter the above Hamiltonian and Casimir functions, that is, we consider the functions where α and β are smooth, and g 1 , g 2 are real parameters.Then, an integrable deformation of system (3) is given by ( ẋ, ẏ, ż) that is System ( 7) is jerk only if and (8) turns into thus, β y = 0 and, consequently, α y = 0. Therefore, Consequently, if the functions α = α(x) and β = β(x) satisfy the relation then we constructed the following family of Hamilton-Poisson jerk systems The corresponding jerk equation is given by .
The jerk versions of the most known oscillators (the harmonic oscillator ẍ + x = 0, the mathematical pendulum ẍ + sin x = 0, and the Duffing oscillator ẍ + x 3 − x = 0) belong to the above family of jerk equations.
In the following, we give Hamilton-Poisson realizations of system (10).For this purpose, using (2), the functions give the matrices and respectively.
Proof.Using ( 12)-( 14), it is easy to see that system (10) writes The sum of the matrices Π 1 and Π 2 is a Poisson structure.Therefore, Π 1 and Π 2 are compatible Poisson structures, and ( 10) is a bi-Hamiltonian system.
As a consequence, we obtain the next result.
the Poisson structure is defined by , and a Casimir of the Poisson structure is Since ∇C(x, y, z) = (g 2 β ′ (x), 0, 1) ̸ = (0, 0, 0), for all (x, y, z) ∈ R 3 , every level set of the Casimir function C is a regular surface.We denote such a level set by The regular symplectic leaves associated with the Poisson structure Π 1,0 are given by the connected components corresponding to pre-images of regular values of the Casimir function C. Therefore, O c is the regular symplectic leaf of the Poisson structure Π 1,0 corresponding to the regular value c ∈ R of C. In addition, the dynamics of the Hamilton-Poisson system R 3 , Π 1,0 , H are foliated by these symplectic leaves.Moreover, the restriction of system (10) to a regular leaf O c is the following completely integrable Hamiltonian system (O c , ω = dp ∧ dq, H| O c ), where the Hamiltonian H| O c = H(p, q) is given by The reduced equations are where W ′ (q) = c − g 2 β(q).Thus, on each level set O c the dynamics of system (10) are given by system (16) or Equation (17), representing a nonlinear oscillator with the kinetic energy T = 1  2 p 2 and the potential energy V(q) = −W(q) = g 1 α(q) + g 2 qβ(q) − cq (for details about the system ẍ = f (x), see, e.g., [10,20]).

Some Dynamical Properties
In this section, we study the stability of system (10) and we prove the existence of some periodic orbits.Also, we obtain sufficient conditions for the existence of heteroclinic and homoclinic orbits.
(iii) If β ′ (M) = 0 and there is a neighborhood V ⊂ R of M such that F(x) < 0, for all x ∈ V \ {M}, then the equilibrium point e M is nonlineary stable.
(iv) If β ′ (M) = 0 and there is a neighborhood V = (a, b) ⊂ R of M such that F(x) > 0, for all x ∈ (a, M) or x ∈ (M, b), then the equilibrium point e M is unstable.
Proof.The Jacobian matrix of system ( 10) at e M is with the characteristic polynomial and eigenvalues (21), it results that one of the eigenvalues is a positive number.
Therfore, e M is an unstable equilibrium point.
If g 2 = 0, system (10) becomes (3), and it has the solution Thus, e M is an unstable equilibrium point.
From the Arnold stability test, it results that the equilibrium point e M is nonlineary stable for g 2 β ′ (M) > 0. (iii) Let U ⊂ R 3 be a neighborhood of (M, 0, 0) such that {x|(x, 0, 0) ∈ U} = V.We consider the function and we prove that it is a Lyapunov function.By the condition L(x, y, z) = 0, we obtain Then, using the hypothesis, we deduce that x = M and y = 0. Therefore, L(x, y, z) = 0 on U if and only if x = M, y = z = 0, that is, L given by ( 22) is a positive definite function on U.Moreover, by (10) we obtain L = ∇L • ( ẋ, ẏ, ż) = 0; thus, L is a Lyapunov function.Therefore, the equilibrium point (M, 0, 0) is nonlineary stable.
Remark 1.If g 2 β is an increasing function such that β ′ (M) = 0, then the function F fulfills the hypothesis given in Theorem 4 (iii); thus, the equilibrium point (M, 0, 0) is nonlinearly stable.
The next result shows the existence of a family of periodic orbits around some nonlinearly stable equilibrium points.
In the following, we study the existence of homoclinic and heteroclinic orbits of system (10).
Let us consider an arbitrary unstable equilibrium point (M, 0, 0), M ∈ R of system (10), which is a saddle, that is, g 2 β ′ (M) < 0. A homoclinic or heteroclinic orbit is given by the intersection of the level sets C(x, y, z) = C(M, 0, 0) and H(x, y, z) = H(M, 0, 0), provided it exists.In this case, we can reduce system (10) to where the smooth function F is given by (18).We have Considering only the level set C(x, y, z) = C(M, 0, 0), system (10) reduces to for which an equilibrium point is (q * , 0), if F ′ (q * ) = 0.The above system writes ẍ = F ′ (x), and it is given by the Hamiltonian H = 1 2 y 2 + V(x), where V(x) = −F(x) is the potential energy.Therefore, for a given function F, "a look at the graph of the potential energy is enough for a qualitative analysis of such an equation" [10].In addition, "if there are two saddle points with the same energy level, corresponding to two maxima of V(x), with no higher maximum between them, then they must be connected by heteroclinic orbits" [20].
The motion of the particle is confined to the region F(x) ≥ 0, and the points with the property F(x) = 0 determine the bounds for the motion.Because heteroclinic and homoclinic orbits are bounded, we assume there is b > M such that F(b) = 0 and F(x) > 0 for all x ∈ (M, b).Since F(M) = 0, we obtain that F has at least a local maximum N ∈ (M, b); hence, (N, 0, 0) is an equilibrium point.Moreover, F is concave in a neighborhood of N; thus, (N, 0, 0) is a nonlinearly stable equilibrium point (via Theorem 4 (ii); note that F ′′ (x) is the same for all equilibrium points).Then, we obtain the next result.Theorem 6.Let α, β be smooth functions such that g 1 α ′ (x) = −g 2 xβ ′ (x), where g 1 , g 2 ∈ R and (M, 0, 0), M ∈ R represent an arbitrary unstable equilibrium point of system (10) such that g 2 β ′ (M) < 0. We consider the function F defined in (18), that is, F( Assume there is b > M such that F(b) = 0, F(x) > 0 for all x ∈ (M, b), and the function F does not have local minima on (M, b).(i) If F ′ (b) = 0 and g 2 β ′ (b) < 0, then a heteroclinic orbit HE(t) = (x(t), y(t), z(t)) given by (25) exists, which connects the unstable equilibrium points (M, 0, 0) and (b, 0, 0).(ii) If F ′ (b) ̸ = 0, then a homoclinic orbit H(t) = (x(t), y(t), z(t)) given by (25) exists, which connects the unstable equilibrium point (M, 0, 0) with itself.
Remark 2. The above theorem also holds for b < M.

The Anharmonic Oscillator
In this section, we apply the obtained results to the jerk version of the anharmonic oscillator given by the equation ẍ + δx n = 0, where δ ̸ = 0 and n > 1 integer.
Therefore, system (28) belongs to the considered family of Hamilton-Poisson jerk systems (10) if The constants of motion are given by The stability of the equilibrium points follows by Theorem 4.
(i) If δM n−1 < 0, then the equilibrium point e M is unstable.
(ii) If δM n−1 > 0, then the equilibrium point e M is nonlinearly stable.
Around some nonlinearly stable equilibrium points, there is a family of periodic orbits of the considered system.More precisely, by Theorem 5 we deduce the next result.Proposition 2. Let e M = (M, 0, 0) be a nonlineary stable equilibrium point of system (28) in the case δM n−1 > 0.Then, for each sufficiently small ϵ ∈ R * + , any integral surface As we have seen in Theorem 6, some homoclinic or heteroclinic orbits can exist in the considered dynamics.Proposition 3. Let e M = (M, 0, 0) be an unstable equilibrium point of system (28).If n is even and δM < 0, then a homoclinic orbit H exists that connects the unstable equilibrium point (M, 0, 0) with itself.Moreover, the heteroclinic orbits cannot exist in this case.
Proof.Let nδM n−1 < 0. We consider the function F defined in (18), namely, Using F ′′ (x) and F ′ (x), we deduce the following: (a) Let δ < 0 and M > 0.Then, there is an unique > 0 for all x ∈ (b, M), and F(x) < 0 otherwise.Using Theorem 6, a homoclinic orbit H exists that connects the unstable equilibrium point (M, 0, 0) with itself.Moreover, the heteroclinic orbits cannot exist.(b) If δ > 0 and M < 0, then we obtain the same result on (M, b), which finishes the proof.Remark 3. If n is odd, then (M, 0, 0), M ̸ = 0 is an unstable equilibrium point for δ < 0. In this case, the above-mentioned function F has the property F(M) = 0 and F(x) > 0 otherwise.Thus, the motion of system (28) is unbounded.
As a particular case, we consider n = 2, that is, Let δ > 0. Thus, the equilibrium point e M = (M, 0, 0), M > 0 is nonlinearly stable, and a family of periodic orbits of the above system surrounds it (white curves in Figure 1).Choosing initial conditions farther and farther from e M , these periodic orbits approach the unstable equilibrium point e −M , that is, they tend towards the homoclinic orbit that connects the unstable equilibrium point e −M with itself (the pink curve in Figure 1).After that, the unbounded curves appear in the dynamics of system (31) (yellow curves in Figure 1).periodic orbits (white) around the stable equilibrium point (M, 0, 0), M > 0 (red), a homoclinic orbit (pink) that connects the unstable equilibrium point (−M, 0, 0) (green) with itself, and unbounded curves (yellow).

Theorem 3 .
Let a, b, c, d ∈ R such that ad − bc = 1.Then, system (10) admits infinitely many Hamilton-Poisson realizations R 3 , Π a,b , H c,d , where the Hamiltonian H c,d is given by

contains at least one periodic orbit γ e M ϵ
of system (28) whose period is close to