New Chaotic Systems with Two Closed Curve Equilibrium Passing the Same Point: Chaotic Behavior, Bifurcations, and Synchronization

: In this work, we introduce a chaotic system with inﬁnitely many equilibrium points laying on two closed curves passing the same point. The proposed system belongs to a class of systems with hidden attractors. The dynamical properties of the new system were investigated by means of phase portraits, equilibrium points, Poincaré section, bifurcation diagram, Kaplan–Yorke dimension, and Maximal Lyapunov exponents. The anti-synchronization of systems was obtained using the active control. This study broadens the current knowledge of systems with inﬁnite equilibria.

According to a new classification of chaotic dynamics [16], there are two kinds of attractors: self-excited attractors and hidden attractors. Recall that an attractor is referred to as being self-excited if its basin of attraction intersects any arbitrarily small open neighborhood of an equilibrium, otherwise it is called a hidden attractor. The basin of attraction for a hidden attractor is not connected with any unstable fixed point. For example, hidden attractors are observed in the systems without fixed points, with no unstable fixed points, or with one stable fixed point. A system with infinitely many equilibrium points can be classified as one system with hidden attractors, for the reason that we do not know which part of the equilibria may be used to localize the hidden attractors, which should be treated in detail (see, e.g., [17]). Recent important investigations and developments in the study of chaotic dynamical systems with practical problems and challenges have been asked to satisfy at least one of the following criteria as Sprott mentioned in [18]: (S1) The system should credibly model some important unsolved problem in nature and shed light on that problem; (S2) the system should exhibit some behavior previously unobserved; (S3) the system should be simpler than all other known examples exhibiting the observed behavior. An important ongoing research topic is dedicated to discovering and developing new and novel chaotic systems with different shapes of closed curve equilibrium.
The main goal of this work is to present a new system with infinitely many equilibrium points arranged on two closed curves passing the same point, which extends the general knowledge about such systems. Our new chaotic system (see Section 2 below for details) is meaningful for satisfying two of the three conditions, (S1), (S2), and (S3), as well as there being a certain novelty value in this work. In Section 2, some dynamical properties of the proposed system, which have been studied using a bifurcation diagram, phase portrait, Poincaré section, maximal Lyapunov exponents, and Kaplan-Yorke dimension, are presented. The ability of anti-synchronization of the new system is also discussed in Section 3.

A New Family with Two Closed Curve Equilibrium
In this work, motivated by the known dynamic systems mentioned above, we study the following general model given byu where u, v and w are three state variables, f (u, v, w) and g(u, v) are two nonlinear functions. The equilibrium points in model (1) can be obtained by calculating It is obtained that the equilibrium points locate on a curve described by g(u, v) = 0 in the plane w = 0. In fact, by selecting appropriate functions f and g, some known systems, both chaotic and with different closed curve equilibrium, can be constructed.

(Example C)
Very recently, Zhu and Du [13] discovered and studied a new family of systems with different equilibrium (as shown in Figure 3) described bẏ where k ∈N. In fact, the chaotic system (5) can be obtained by putting f (u, v, w) = αv + βv 2 + uw and g(u, v) = |u| k + |v| k − 1 into model (1). In [13], Zhu and Du analyzed the dynamical properties of their proposed systems using the methods of equilibrium points, eigenvalues, phase portraits, maximal Lyapunov exponents, and Kaplan-Yorke dimension; see [13] for more details. The results established in [11,13,19] are very important for indicating the existence of chaotic systems with different shapes of equilibrium points (see Table 1). Note that the first two equations of the systems proposed in [13,19] are the same as in [11]. The difference is the third equation. When we choose a different third equation, we can get different systems to display new features, such as different shapes of equilibrium point and other dynamic properties.
To the best of our knowledge, there is no paper devoted to the study of chaotic dynamical systems with eye-shaped curve equilibrium. Therefore, this study is an important ongoing research topic. In this paper, motivated and inspired by this, two functions, f (u, v, w) and g(u, v), are chosen in the following forms where α and β are two positive parameters. Substituting (6) into system (1), our new system is described asu It is verified that system (7) has infinitely many equilibrium points (u * , v * ). These equilibrium points are located on the curve in the coordinate plane described by It means that the new system (7) has eye-shaped curve equilibrium as shown in Figure 4. Note that the eye-shaped curve is different from some other shapes reported, such as line, square, circle, or cloud-shaped [11,19], and is symmetric about the u-axis, v-axis, and origin. Furthermore, system (7) has hidden attractors [17]. Above all, investigating system (7) will strengthen our understanding of hidden attractors. For α = 5, β = 30, and initial conditions (0.06, 0.01, 0.01), the new system (7) has chaotic attractors (see Figures 5 and 6). For the simulation, we used the Wolf et al. method to calculate the Lyapunov exponents [20], the time of computation was 1000, and we obtained the Lyapunov exponents (0.0424, 0, −0.2484). The method of Wolf et al. is rooted conceptually in a previously developed technique that could only be applied to analytically defined model systems to monitor the long-term growth rate of small volume elements in an attractor. In addition, the corresponding Kaplan-Yorke dimension of system (7) is 2.1707. Poincaré return maps are often used to transform complicated behavior of a dynamic system in phase space to discrete maps in a lower dimensional space to reveal the complicated behaviors. Poincaré return maps corresponding with phase portraits in Figure 6 are presented in Figure 7; there are some dense points in the Poincaré section, and it can be determined that the motion is a chaotic state. These results reveal that the system is chaotic.   Gradually changing the value of the parameter β or α, the bifurcation plot of the system can be discovered in Figure 8. Figures 9 and 10   The new system with eye-shaped equilibrium has periodic behavior in the range 36 ≤ β ≤ 48. For instance, the system can display period-1 behavior for α = 5, β = 45, period-2 behavior for α = 5, β = 38, and period-4 behavior for α = 5, β = 36 (see Figure 11a-c, respectively).

Anti-Synchronization of New Systems
Synchronization of chaos is a phenomenon that may occur when two, or more, dissipative chaotic systems are coupled. Since the pioneering work of Pecora and Carroll related to synchronization in chaotic systems [21], some methods of chaotic synchronization have been presented related to complete, generalized, lag, and imperfect phase synchronization [22]. Many papers on applications of chaos synchronization for cryptographic [23], kinetics [24], physiology [25], neural networks [26], and economics [27] have appeared.
In the following, we consider the anti-synchronization of the systems with eye-shaped equilibrium related to the driver-response system. The driver system with eye-shaped closed curve equilibrium is as follows:u where u, v, and w are are three state variables, and the value of α = 5, β = 30.
The response system is described aṡ where the control is In order to reveal the difference between the driver system (7) and the response system (9), the state errors can be defined as and we obtainė 1 =u +u 1 , e 2 =v +v 1 , e 3 =ẇ +ẇ 1 .
Combining (7), (9), (10), and (11), we get the state errors systeṁ e 1 = e 3 + h 1 , We choose the control proposed by where k i > 0 (i = 1, 2, 3) are the positive gain constants used to control the rate of anti-synchronization. By substituting (12) into (11), we get the state errors systeṁ Obviously, the eigenvalues (−k 1 , −k 2 , −k 3 ) of the Jacobian matrix of the state errors system are negative. Then, the complete anti-synchronization between the driver system (7) and the response system (9) is proved.
In numerical simulation, we assume the initial values of the driver system (7) and the response system (9) The positive gain constants here are selected as k 1 = k 2 = k 3 = 3. It is obvious in Figure 12 that there exists anti-synchronization of the respective states of the new systems with two closed curve equilibrium (7) and (9). The time history of the synchronization errors e 1 , e 2 , e 3 is shown in Figure 13 which plots the anti-synchronization of the driver-response system.

Conclusions
In this work, we propose and study the following new system: u = w, v = −w(αv + βv 2 + uw), w = u 2 − |u| + |v| + v 2 , with eyed-shaped equilibrium points which are located on two closed curves passing the same point. In Section 2, some dynamical properties of the proposed system are presented, which were investigated using bifurcation diagram, phase portrait, maximal Lyapunov exponents, and Kaplan-Yorke dimension. Furthermore, the anti-synchronization of systems is obtained by using active control in Section 3. This study will broaden the current knowledge of chaotic systems with infinitely many equilibria.
Author Contributions: Both authors contributed equally to this work. Both authors read and approved the final manuscript.