Exploring Limit Cycle Bifurcations in the Presence of a Generalized Heteroclinic Loop

: This work revisits the number of limit cycles (LCs) in a piecewise smooth system of Hamiltonian with a heteroclinic loop generalization, subjected to perturbed functions through polynomials of degree m . By analyzing the asymptotic expansion (AE) of the Melnikov function with ﬁrst-order M ( h ) near the generalized heteroclinic loop (HL), we utilize the expansions of the corresponding generators. This approach allows us to establish both lower and upper bounds for the quantity of limit cycles in the perturbed system. Our analysis involves a combination of expansion techniques, derivations, and divisions to derive these ﬁndings


Introductory Notes
The theory of limit cycle (LC) bifurcation in piecewise smooth integrable differential systems finds widespread application in numerous practical problems, spanning domains like electronic engineering, neural networks, automatic control, and biological mathematics, among others [1,2].Additionally, this theory bears relevance to Hilbert's 16th problem, adding further significance to its exploration.As a result, the investigation of LC bifurcation within differential systems (piecewise smooth) has emerged as a prominent and dynamic area of research in recent years [3][4][5].
The prevalence of natural laws and diverse influencing factors has led to the extensive application of LC bifurcation theory in piecewise smooth integrable differential systems across various practical fields.Engineers often rely on this theory to address complex challenges in electronic engineering, while researchers in neural networks and automatic control also find it instrumental in their work [6,7].
Beyond its wide-ranging applicability, the theory's connection to Hilbert's 16th problem and Weakened Hilbert's 16th problem underscores its fundamental importance in mathematics.Consequently, researchers have increasingly directed their attention toward the study of LC bifurcation within differential systems, recognizing the richness and relevance of this subject.The growing body of work in this area is evident in recent literature, with various studies exploring different aspects and applications [8,9].As a result, the pursuit of understanding and harnessing the dynamics of LCs in such systems has become a vibrant and active research pursuit, attracting researchers from diverse backgrounds and interests.
Considered the following perturbed piecewise smooth Hamiltonian systems [10,11]: where where a δ,ζ , b δ,ζ are perturbation parameters.As long as ε = 0, the associated Hamiltonian functions to (1) can be written as and The unperturbed version of (1) | ε=0 has a family of orbits with periodicity as comes next: where h ∈ (0, 1  2 ).The dynamics of the system are characterized by the presence of pe- riodic orbits situated between two key points.Firstly, there is the origin, which serves as an elementary center of the parabolic-parabolic type, see [12,13] for more.Secondly, we observe a generalized heteroclinic loop (HL) encircling the origin, with two saddle points located at (±1, 0).This configuration of periodic orbits is visually represented in Figure 1.Noting that the unperturbed piecewise smooth Hamiltonian [14,15] systems find application in various areas of physics and engineering where dynamical systems exhibit discontinuous behaviors or interactions.These systems can describe physical scenarios where the underlying dynamics switch between different regimes due to certain thresholds or conditions.Examples include mechanical systems with impacts or friction, electrical circuits with switches, and biological systems with abrupt transitions.The unperturbed nature of the system refers to its behavior without any external disturbances, allowing for the study of the inherent dynamics and the effects of the discontinuities or switching on the system's behavior [16,17].The system's behavior is intriguing as it exhibits periodic orbits positioned in a delicate balance between the elementary center at the origin and the encircling generalized HL [18].At the origin, the dynamics manifest as a parabolic-parabolic type, contributing to the system's unique characteristics.Meanwhile, the presence of saddle points at (±1, 0) within the generalized HL introduces an additional layer of complexity to the system's behavior.The interplay between these fundamental elements gives rise to the observed periodic orbits, which showcase captivating and intricate trajectories in the phase space.To gain further insights into the system's dynamics and the significance of these periodic orbits, refer to Figure 1, providing a visual representation of this intriguing phenomenon.Figure 1 encapsulates the captivating phenomenon occurring within the system under study.The figure highlights the existence of periodic orbits that reside in the region between the origin and a generalized HL surrounding the origin.While the origin itself serves as an elementary center of parabolic-parabolic type, the loop features saddle points at (±1, 0), adding a distinct dynamic dimension to the system.As a result, the periodic orbits gracefully navigate this delicate balance between the attractive forces of the origin and the more complex interactions near the saddle points.Understanding the underlying principles governing these orbits is vital for comprehending the system's overall behavior and offers valuable insights into its practical applications and implications.By the main results in [19,20], we know that the first order Melnikov function (FOMF) associated with the periodic orbits of (1) as comes next ( A fascinating connection exists between the LCs and zeros of the FOMF in the context of system (1).Specifically, the LCs of the system are associated with isolated zeros of the first non-vanishing Melnikov functions (4), see fore more [21,22].
In a related work by Liang and Han [23], they delved into the study of Hopf bifurcation and generalized heteroclinic bifurcation, managing to derive an upper bound for the count of LCs in system (1) up to the first order in ε.To facilitate their calculations, the authors made simplifications by selecting specific forms for f + (s, z) and f − (s, z) in each bifurcation case.For the Hopf bifurcation, they used , and d ± δ,ζ are the perturbation coefficients.On the other hand, for the generalized heteroclinic bifurcation, they employed and 2016) expanded upon the Melnikov method to address a specific category of perturbed hybrid piecewise smooth systems featuring two switching manifolds in a planar configuration.Their study delved into the continuity of periodic orbits, particularly when the unperturbed system itself possesses such orbits [24].Other new results about periodic orbits of piecewise smooth systems can be traced in [25][26][27].
In this paper, we take a fresh approach to revisit the problem of LC bifurcation in (1) by considering the more general case of perturbed polynomials f ± (s, z) and g ± (s, z).Unlike the previous study, we do not impose specific constraints on the form of these polynomials, allowing for a more comprehensive exploration of the system's dynamics.
This paper unfolds with a systematic investigation into the intricate dynamics of the FOMF M(h) in Section 2. This section delves into the algebraic properties and underlying structure of M(h), shedding light on its essential characteristics and revealing valuable insights into its behavior.Building upon this foundation, Section 3 embarks on an exploration of the asymptotic expansion (AE) of M(h) near pertinent regions of interest.This critical analysis provides a deeper understanding of the function's behavior and unveils valuable information about its asymptotic properties.More results are obtained because of more general perturbations (general perturbed polynomial differential certain system).Finally, the paper culminates the conclusion in Section 4, summarizing the key findings and implications of the study.Together, these sections present a comprehensive investigation into the intriguing relationship between LCs and the FOMF, contributing to the broader understanding of dynamical systems and bifurcation phenomena.

The Algebraic Structure of the FOMF M(h)
Our major result consists of the theorem presented below.
Theorem 1.The count of LCs arising from system (1), originating from the period annulus near the origin and encompassing all valid combinations of f ± (s, z) and g ± (s, z) adhering to ( 2 Toward such a goal, we build several different theoretical findings.Now we investigate the algebraic structure of FOMF M(h) for (1).When h ∈ (0, 1  2 ), we denote Noting that with respect to the x-axis (s-axis), the orbits Γ ± h are symmetric.Therefore, I δ,2ζ+1 (h) ≡ 0. The following lemma is now stated and proved.Lemma 1.For h ∈ (0, 1  2 ) and when m ∈ N, we have ) are constants and independent of each other.Here, [•] stands for the integer part.
Proof.Using the Green's Formula, one has the following and By one obtains ( 7) and ( 8), one obtains Hence, M(h) can be given by It is easy to verify that ξ δ,ζ might be chosen as free constants.To continue the proof, we establish the following relations: where ) can be taken as free parameters.Indeed, upon differentiating both sides of H + (s, z) = h with respect to z, one attains z − s ∂s ∂z Upon multiplying (11) by the one-form s δ z ζ−1 dz and subsequently integrating, we arrive at the following relation: Similarly, multiplying H + (s, z) = h both sides by s δ−2 z ζ dz and integrating over Γ + h , we obtain another relation Elementary manipulations reduce Equations ( 12) and (13) to and We will now establish the claim using induction on m.To preserve the generality of the argument, we can make the following assumption and move on.We will focus on proving the claim for even values of m (similarly, the claim can be proven for odd values of m).To begin, let us perform a direct computation using the two equalities mentioned above: Hence, for m = 2, 3 one has Indeed, the claim has been verified for m = 2, 3. Let us now make the assumption that the claim holds for all values of δ 2) and (δ, ζ) = (m, 0) into ( 16), we obtain: . . .
Thus, by leveraging the hypothesis of induction along with (17), one can deduce: where σδ , θδ , σ δ and θ δ are constants and (2.6) is proved in view of I 0,0 (h) = −2 √ 2h.Next, we establish the independence between ).To begin, we utilize the induction hypothesis, which ensures that σδ ) are already established as independent.In other words, the determinant of the following Jacobian matrix governs their independence and is not zero ,l 0 .
Here, the summation for ξ δ,ζ is less or equal than m − 1.This follows directly from the calculations, and we can express the Jacobian matrix as follows: ,l 0 , ξ 0,m

m+1
, where B is a column vector and 0 is a row vector.We now obtain ) are independent of each other.This completes the proof.
Before ending this section, it is recalled that the work [28] introduces an extension of the Melnikov function method and the averaging method to study the presence of LCs in piecewise smooth near-integrable systems, even in higher dimensional spaces.Their work demonstrates the equivalence of these methods, which were previously known for planar analytic systems, and establishes their applicability for piecewise systems.That lies in deriving the formula for the second-order Melnikov function in planar piecewise near-Hamiltonian systems.

The Asymptotic Expansion of M(h)
In this section, our focus shifts towards providing the essential asymptotic expansions of the FOMF M(h) in close proximity to both the generalized HL Γ and the origin.These expansions offer valuable insights into the behavior of M(h) within these critical regions and lay the groundwork for a deeper understanding of the system's dynamics.
To begin, let us explore the AE of M(h) near the generalized HL Γ.By employing appropriate mathematical techniques, we derive expressions that shed light on the intricate relationship between M(h) and Γ.The obtained expansions reveal how the Melnikov function behaves as h approaches Γ and provide valuable information about the system's response to perturbations in this particular region.Understanding the behavior of M(h) near Γ is crucial for comprehending the system's stability and the emergence of LCs in this vicinity.
Next, we direct our attention to the AE of M(h) near the origin.Employing rigorous mathematical analysis, we unveil the behavior of M(h) as it approaches the origin.The obtained expansions offer valuable insights into the system's dynamics in the vicinity of this critical point and provide a deeper understanding of the underlying mechanisms governing LC bifurcations.By unraveling the intricate relationship between M(h) and the origin, we gain essential knowledge about the stability and behavior of the system for perturbations close to the origin.
Following the presentation of the asymptotic expansions, we introduce a key lemma that contributes to our understanding of the system's behavior in the region where 0 < 1 2 − h 1.The lemma provides an explicit AE of I 1,0 (h), a vital component of the Melnikov function, as h approaches the neighborhood of 1  2 .This expansion, expressed in terms of real constants ρ 1 and ρ 1k , with k ≥ 0, illuminates the intricate characteristics of I 1,0 (h) and its relationship with h in this specific region.Such detailed insights into the behavior of I 1,0 (h) contribute significantly to the overall understanding of the Melnikov function and its role in LC bifurcations near 1  2 .
Proof.Some direct computations show that Let s = √ 2l + 1 and u = h − 1 2 , then one can rewrite (20) into Notice that −1 < l −2 u < 0 for 0 < h < 1 2 and 0 < s < 1 − √ 1 − 2h, and the following AE where (21), one can find that where It is easy to check that the series in ρ 11 is convergent.This ends the proof.

Proof. A direct calculation gives
√ h where .
where c δ and d δ are independent of each other, moreover Proof.Substituting ( 19) and ( 24) into (12) gives (26).A direct computation shows that the determinant of the following Jacobian matrix , is 1, which means that the coefficients c δ and d δ can be arbitrarily selected.With this, we have concluded the proof.
At this moment, we are in a position to establish the major finding of this work using the obtained theoretical findings in what follows.

Proof of Theorem 1
By Proposition 1, we know that c ) are independent of each other.We can select them in such a way that they read the conditions below: ] + 1 = m simple and distinct zeros in (0, 1  2 ).Consequently, system (1) can exhibit m LCs for 0 < h < 1 2 .Now, our objective is to determine the upper bound for the number of LCs of system (1).To achieve this, we begin with some preliminary considerations.

Concluding Summaries
In this work, we have investigated the bifurcations and stability of LCs in a general perturbed polynomial differential system (1).Our main focus was on understanding the impact of perturbation parameters on the existence and number of LCs in the system.
In Section 2, we studied the algebraic structure of the FOMF M(h) for the system (1).By appropriately selecting the perturbation coefficients, we showed that the number of LCs can be controlled to be m, which represents an exciting result in the study of perturbed polynomial differential systems.We provided a clear theoretical foundation for this result and proved it rigorously.
In Section 3, we analyzed the AE of the FOMF M(h) in two critical regions: near the generalized HL Γ and near the origin.These expansions provided valuable insights into the behavior of M(h) in these regions and its relation to the underlying dynamics of the system.Additionally, we introduced a key Lemma 2 that presented an explicit AE of I 1,0 (h), which plays a significant role in the Melnikov function near the value 1  2 .The obtained results contributed significantly to our understanding of the system's stability and the emergence of LCs in these critical regions.
In conclusion, our work sheds light on the rich dynamical behavior of perturbed polynomial differential systems and provides a deeper understanding of their stability and bifurcation properties.The findings presented in this work contribute to the existing body of knowledge in the field of dynamical systems theory and can be applied to a wide range of real-world applications, including engineering, physics, and biology.
Future research in this area could explore the behavior of LCs and stability in higherdimensional perturbed polynomial systems, as well as investigate the effects of other types of perturbations on the system's dynamics.Additionally, numerical simulations could be conducted to validate and further explore the theoretical findings presented in this paper.