We investigate the Shilnikov sense homoclinicity in a 3D system and consider the dynamical behaviors in vicinity of the principal homoclinic orbit emerging from a third order simplified system. It depends on the application of the simplest normal form theory and further evolution of the Hopf-zero singularity unfolding. For the Shilnikov sense homoclinic orbit, the complex form analytic expression is accomplished by using the power series of the manifolds surrounding the saddle-focus equilibrium. Then, the second order Poincaré map in a generally analytical style helps to portrait the double pulse dynamics existing in the tubular neighborhood of the principal homoclinic orbit.

1. Introduction

In recent years, there has been a great deal interest in understanding the dynamics and bifurcations in problems governed by 3D autonomous differential systems [13]. These problems include sequences of periodic doubling, existence of strange attractors, and bifurcations of the multipulse orbits. The common adapted tool for analyzing multipulse behaviors of the 3D nonlinear dynamical system is the well-known Poincaré map [4]. It starts by a standard analysis of the local and global cross-sections to the principal homoclinic orbit of the equilibrium, where the stable and unstable manifolds can intersect along the tubular neighborhood of the principal homoclinic loop, giving rise to Poincaré homoclinic structure. Among the literature, Glendinning and Sparrow [5] studied the local behavior of systems' near homoclinic orbits to the stationary equilibrium of saddle-focus type. They described the process of periodic orbit approached homoclinicity and discussed the resulting global patterns of bifurcations, such as the emergency of subsidiary homoclinic manifolds. More explicit work is offered by Gonchenko et al. [6]. They reported the research of multipulse homoclinic loops in two-parameter families of vector fields surrounding the principal homoclinic tangency and then proved the high structural instability in the neighborhood of a Shilnikov type saddle-focus.

The prosperities of the Poincaré map reflect the topological and stability of the actual flow and are basically adapted in this paper to concern the dynamics of double-pulse orbit at the tubular field. Besides that we complexly quantify the principal orbit before the analysis of the return map with the expectation to generalize the expression and obtain higher computational precise. Actually, the so-called principal homoclinicity has been frequently mentioned [79], but the orbit itself and its usage for the return map have not been broadly detected. Li and Zhu constructed the dimensional series expression of the stable and unstable manifolds [10], which formed the homoclinic orbit near the hyperbolic singular point of the Lorenz system.

With a further extension of the methodology, the Shilnikov sense principal homoclinic loop will be analytically obtained from the composition of series expression of the invariant manifolds surrounding the saddle-focus equilibrium. These quantified manifolds intersect the cross-sections locally and globally produce more explicit coordinates up to a desired order and also characterize the Poincaré homoclinic map of the system.

When talking about the complex dynamics, we acknowledge the relevance of the splitting homoclinic connection with the unfolding of the Hopf-zero singularity [11]. Meanwhile, the local dynamics presented by a normal form approach (NF) [1215], a standard result of Birkhoff NF theory, is also applied in the initial simplification.

So, in this paper, we try to construct a second order Poincaré map based on the Hopf-zero unfolding in a common 3D system. In Section 2, it is supposed that the system exhibit the Hopf-zero singularity at the fixed parameter value of the linear part of the vector field. In Section 3, the explicit expression of the principal homoclinic loop H is given by using the invariant manifold method. In Section 4, analytical series expressions motivate our interest in chasing a more general and higher order Poincaré map from its quantified stable and unstable manifolds. The new map can be easily degenerated to the original result presented by Wiggins [16]. And it gives further facility of studying the complex dynamics, such as subsidiary or multipulse homoclinicity surrounding the principal homoclinic orbit.

2. The Shilnikov Unfolding Based on the Hopf-Zero Normal Form

Consider the following 3D system: where is regarded as the orbit splitting parameter and are , and vanish along with their derivates at . As mentioned in [11], the unfolding of the Hopf-zero local bifurcation contains the global Shilnikov sense homoclinicity emerging from the saddle-focus equilibrium. Given that (1) enjoys Hopf-zero singularity near the origin, the normal form (NF) can be written as in which , an over-bar denotes complex conjugate; for example, denotes the complex conjugate of ; and are the NF coefficients, given in the complex and real forms, respectively. The final expressions have different forms depending on whether the total order n of the NF is an odd or even number; that is,(i) , if n is odd;(ii) , if n is even.

Here we only use the third order NF in (2) because the system is terminated in the third order [11]:

The dynamics of bifurcation near the structurally unstable homoclinic orbits in the Shilnikov sense requires the description on a two-parameter group . It takes the function as the distance between the point of the intersection of principal homoclinic H with some surface of the cross-section and the line of intersection of stable manifold with the same surface of section. On the other hand, , defined in its complex form. Further discussion is restricted to the condition .

Firstly, it is needed to quantify the analytical structure of the principal homoclinic orbit in the complex form at . And then constructing a second order Poincaré map to prohibit a clear insight of the complex dynamics surrounding the principal orbit is carried out.

3. Analytical Expression of the Principal Homoclinic Orbit

The analytical expression of the principal homoclinic orbit is the combination of the 1D manifold associated with the eigenvalue and the 2D manifold associated with . So, the 1D manifold for (3) is where , and are the coefficients to be determined. is in the complex form and denotes its complex conjugate. The is a real number. Differentiate (4) with respect to time t and substitute it into (3) and then balance the same order terms of on both sides of the resulting equations. That yields

The 2D manifolds subject to the complex enginvalues are where , and are the coefficients to be determined, . can be expressed in the following complex form: . Differentiate (6) with respect to time and substitute it into (3) and then balance the same order terms of on both sides of the resulting equation. It gives the coefficients of the 2D manifold:

These foregoing 1D and 2D invariant manifolds finally form the Shilnikov type principal homoclinic orbit H in its analytical form. It benefits the construction of the second order Poincaré map shown in Section 4.

4. The Second Order Poincaré Map

The Poincaré map is one of the most commonly used approach to study the dynamics surrounding the Shilnikov homoclinic loop including double or multipulse (circuit) orbits. In a general process, the first step to construct a return map is built from the composition of the first hit maps between the surfaces of section [4, 17]. The first hit maps are locally defined since at the surfaces of the section intersect the homoclinic cycle transversally. So, with the obtained analytical manifolds, the local linear stable and unstable manifolds are given by

The orbit intersects the complex plane at a countable set of points at . So, we select as the starting point of (8) in its coordinates and construct a two-dimensional cross-section as shown in Figure 2. Following the orbit, there emerges another point belonging to the second auxiliary cross-section . Let be a rectangle lying in the plane and a rectangle parallel to the plane at . So, the Poincaré map along the orbit should be constructed by dividing the flow close to the stationary point of its unstable manifold into two parts: a small neighborhood of the stationary point where the flow is essentially linear and a global flow which takes trajectories close to the unstable manifold away from and then back to the linear region, as shown in Figure 1.

Starting from the point on to reach on , it produces the first local map . According to (8), this local map can be written as

From the equation of , the time of flight for points starting on to reach can be calculated from

Substituting the flight time into map , it is found that the expression for the first hit map can be expressed as where . Since the time of flight is bounded in this progress, which generates another global map according to (6), where include the remaining components of , and these equations portrait the global behavior outside the equilibrium. So, the combination of the local and global map finally forms the second order Poincaré map:

That is a more general expression of the return map because the former results as given by Gonchenko et al. and Zhou et al. [6, 7], can be easily obtained by letting , and ignoring the second order components, while the new map makes itself in a broader usage because of the general expression and better computational precision.

5. Conclusions

In this paper, the second order Poincaré map in a common 3D nonlinear system is constructed. Shilnikov sense subsidiary homoclinic dynamics is also researched. Firstly, we find the strategy to analytically construct the homoclinic orbit surrounding the saddle-focus equilibrium point in terms of the complex SNF and method of invariant manifolds. By using these quantified manifolds, we can form the expression of the first and the second order Poincaré map, making their appearance in a more general style with the capacity to enhance computational precise. The further possible application of the subsidiary homoclinicity is shown in Figure 2 which is a diagram for subsidiary homoclinic orbit; the return map can portrait the dynamics unfolding from the principal orbit.

However, the analytical work lags behind the numerical results in great extent. So, to employ the higher order Poincaré map and then to unfold the successive bifurcation behaviors more simply and exactly, in a narrow parameter interval, will be the topics for further discussion.


The SNF coefficients of the Hopf-zero bifurcation up to the third order are as follows:

The second order coefficients of the two-dimensional manifold are as follows

Conflict of Interests

The authors declare that there is no conflict of interests.


This work was supported by the National Natural Science Foundation of China (Grant no. 11272229) and Science Foundation of Tianjin education committee (Grant no. 20120902).