Coexisting Oscillation and Extreme Multistability for a Memcapacitor-Based Circuit

The coexisting oscillations are observedwith amemcapacitor-based circuit that consists of two linear inductors, two linear resistors, and an active nonlinear charge-controlled memcapacitor. We analyze the dynamics of this circuit and find that it owns an infinite number of equilibrium points and coexisting attractors, which means extreme multistability arises. Furthermore, we also show the stability of the infinite many equilibria and analyze the coexistence of fix point, limit cycle, and chaotic attractor in detail. Finally, an experimental result of the proposed oscillator via an analog electronic circuit is given.


Introduction
Memristor is known as the fourth basic circuit element that was firstly postulated by Chua in 1971 [1] and has attracted worldwide immense attention from both theory and applications even since a solid-state device called the memristor was fabricated by the HP Lab in 2008 [2].Then, Di Ventra et al. extend the notion of memristor to memcapacitor and meminductor [3], whose properties depend on the state and history of them in circuit.These nanoscale elements, that store information without need for an internal power supply, can be used for nonvolatile memories, nonlinear circuits, neural networks, and so on.
In the field of nonlinear dynamics, memristors can be as fundamental elements for the designs of new nonlinear circuits substituting other nonlinear devices, whereupon many memristor-based chaotic oscillators were presented by using HP memristor and the other memristor models of piecewise linear, quadric, and cubic - (flux-charge) functions.Itoh and Chua firstly derived several nonlinear oscillators based on Chua's oscillators by replacing Chua's diodes with piecewise linear models of memristor [4].The authors in [5] designed a periodically forced memristive Chua's circuit using the flux-controlled memductance model () = − + |()|.HP memristor model is a first model of actual memristive device, which was implemented by a memristor emulator in [6].The paper [7] introduced a complete mathematical model for the HP memristor which takes into consideration all boundary situations providing the interrelation between memristance charge and flux.Based on HP memristor and other memristor models, some chaotic oscillators were presented [8][9][10].Recently, a new memristorbased chaotic system was designed based on HP memristor in [11] and was realized using FPGA (Field Programmable Gate Array) technology.A novel digital-analog hybrid chaotic system with generalized memristor was constructed for the production of random number [12].However, the memcapacitor and meminductor have received a little attention since the actual solid-state devices of them have not been successfully achieved.Because of their potential application values, memcapacitor and meminductor have still attracted more and more attentions.
The first meminductor emulator was designed in [13] whose inductance can be varied by an external current source without employing any memristive system.In [14], a flux-controlled memristive emulator using light-dependent resistor (LDR) was proposed and the mutator for transferring memristor into a flux-controlled meminductor is described.

Mathematical Problems in Engineering
Active memcapacitor A mathematical model of a flux-controlled meminductor and its equivalent circuit model for exploring the properties of the meminductor in a nonlinear circuit were presented in [15].It would be worthwhile to prospectively study effective memcapacitor models and its applications.So several memcapacitor models, including piecewise linear, quadric and cubic function models, memristor-based memcapacitor models, and memcapacitor emulators, were proposed in [16][17][18][19][20][21], and a mathematical memcapacitor model and a corresponding circuit model are established in [22].Some special phenomena such as hidden attractors, coexistence attractors, and extreme multistability were found in memcapacitorbased chaotic oscillators [23,24] and memristor-based chaotic oscillators [25][26][27].In fact, multistability and coexisting attractors have caught the attention of researcher in general chaotic systems [28][29][30].
This paper introduces a new chaotic oscillator based on a charge-controlled memcapacitor model and its dynamical behaviors are analyzed.The most important properties of the memcapacitor-based circuit are that it possesses an infinite number of equilibrium points and coexisting attractors and displays the coexisting attractors and stability of the infinite many equilibria, called extreme multistability.The rest of this paper is organized as follows.In Section 2, the chaotic oscillator is described and the typical chaotic attractors are given.In Section 3, the dissipativity and equilibrium stability of the system are studied.In Section 4, coexisting attractors with different initial (0) and (0) are described.In Section 5, an analog circuit is designed to realize the memcapacitorbased oscillator.Finally, some conclusions of this paper are given in Section 6.

Memcapacitor-Based Chaotic Oscillator Circuit
The chaotic oscillator circuit based on a charge-controlled memcapacitor is shown in Figure 1, which contains two linear inductors,  1 and  2 , two linear resistors,  1 and  2 , and an active memcapacitor.The active memcapacitor is an active 2terminal circuit consisting of a linear active resistor − 0 and a nonlinear memcapacitor in series, in which the active resistor − 0 provides energy for the circuit.Hence, the series circuit, that is, the active 2-terminal circuit as a whole, exhibits an active feature, so called the active memcapacitor.
The circuit dynamics can be described by where   and V  are the charge through the memcapacitor and the voltage across the memcapacitor, respectively.The relationship between the charge   and voltage V  of the memcapacitor is defined as [19] V where  +  2  is the inverse memcapacitance and   is the time integral of   ; namely,  = ∫ ( Typical chaotic attractors obtained by simulating (3) are shown in Figure 2 in terms of the parameters listed in Table 1, with the initial condition (0.01, 0.01, 0.01, 0.01).
The corresponding Lyapunov exponents are calculated as LE 1 = 0.3787, LE 2 = 0.0092, LE 3 = −0.0094,and LE 4 = −0.9786.The Poincare mapping on  = 0 and the timedomain waveforms of system (3) are described in Figure 3.  Mathematical Problems in Engineering Obviously, time-domain waveforms are nonperiodic continuous chaotic sequences, and the curve of Poincare map is a continuous curve, and the time-domain waveform is not periodic, which indicate that these attractors are chaotic attractors.

Dynamical Characteristic Analysis
3.1.Dissipativity and Stability of Equilibrium Set.Dissipative systems will shrink to a phase space of a limited area, called the attractor.The dissipativity of the system is described as When the parameters are set as  = 2.6 and  = −3.4,the exponential constraint rate satisfies It means that the volume of phase space will be contracted to zero and all trajectories of the system are confined to a subset of zero volume.
The equilibrium points of system (3) can be calculated as where  1 ,  2 , and  3 are arbitrary constants.Obviously, the system has two equilibrium sets, and every set has an infinite number of equilibria.
If the parameters are  = ,  =  or  ̸ = ,  ̸ = , the Jacobian matrix of the circuit at equilibrium sets  1 and  2 is given by where  2 and  3 have the same effect on the system, so they can be replaced by .The characteristic equation is written as According to the Routh-Hurwitz criterion, the necessary and sufficient condition for stability of the system is that all the principal minors are greater than zero; that is, where whereas it is unstable and has the possibility of generating chaos under the proper parameters.If we fix  = 5.8,  = 1.5,  = 2.6,  = 0.1,  = −3.4, = 0.2,  = 6.8, and  = 2.8, the stability of equilibrium sets  1 and  2 depends on the constant . Figure 4 shows the change rule of Δ 1 , Δ 2 , Δ 3 , and Δ 4 with the constant .Here, Δ 1 is slightly greater than zero and Δ 4 is equal to zero, so the black curve overlaps so much with the green curve in Figure 4.Note that the stability of an infinite number of equilibria is called extreme multistability.
The following analysis will show that the stability of the all equilibrium points depends on the parameter , resulting in chaos when  ∈ [−1.5, 1.5].

Coexisting Oscillations.
In this section, we show that the dynamics of the system depend on not only its parameters but also its initial conditions, which leads to a special characteristic, called the coexisting oscillation in this paper.In the coexisting oscillation, many coexisting attractors appear with different initial conditions, such as coexistence of limit circle or fixed point and chaotic attractor, coexistence of chaotic and hyperchaotic attractors, and coexistence of two limits or two fixed points.Multistability is a common phenomenon in many nonlinear systems, which corresponds to the coexistence of more than one stable attractor for the same set of system parameters [31].If infinitely many attractors coexist for the same set of system parameters, the phenomenon is called extreme multistability [32], showing that various coexisting attractors have essentially extreme multistability.

Coexisting Attractors with Initial
Variable (0).If we fix the parameters as Table 1 and vary initial value (0) in the initial conditions ( 0 , 0.01, 0.01, 0.01), the bifurcation diagram and corresponding Lyapunov exponent spectrum with respect to (0) can be obtained as shown in Figure 5. System (3), with the increase of initial state (0) in Figure 5, undergoes a complex bifurcation from chaos to period, or from period to chaos, implying the coexistence of periodic and chaotic attractors.
When (0) is varied, many coexisting attractors appear, as shown in Figure 6, in which every subdiagram represents a pair of coexisting attractors; the blue attractor and the red attractor start from different initial values of (0).The blue attractors start from (0) ∈ (−0.85, 0) and the red attractors start from (0) ∈ (0, 0.84).From Figure 6, we can see that there exist different coexistence attractors, including coexistence of a pair of chaotic attractors, coexistence of a pair of quasiperiodic limit cycles, and coexistence of a pair of limit cycles.Various coexisting attractors and their corresponding initial values are listed in Table 2.

Coexisting Attractors with Initial
Variable (0).Similar to the abovementioned, the bifurcation diagrams and Lyapunov exponent spectrums with respect to the initial value (0) are described in Figure 7, in which we also fix parameters as Table 1 and vary parameter (0) in the initial condition: (0.01, 0.01, 0.01,  0 ).Figures 7(a Coexisting attractors with respect to initial value (0) are shown in Figure 8, and the corresponding regimes and initial conditions for all attractors are listed in Table 3.Note that Figure 8(a) shows the coexistence of the three attractors, that is, fixed point, two different limited circles.
Note that the phenomenon of coexisting oscillation is different with the initial value sensitivity of usual chaotic systems.For the usual chaotic system, trajectories starting from different initial values evolve along with the time and eventually converge to the same chaotic attractor.For the proposed oscillator, however, trajectories starting from different initial values can converge eventually to different attractors, including chaotic attractors, periodic attractors, or fixed points.Hence, it seems that trajectories in the coexisting attractors are more sensitive to initial values not only for chaotic orbits but also for periodic orbits.

Dimension Reduction Analysis
In order to solve the 3D equations (1), we must add an equation of   which is an internal variable of the memcapacitor and obtain the 4D equations (3).However, this makes it difficult to solve the 4D equations.In general, if there is one memory element increased into a circuit, there will be two state equations added to the system equations, in which one is the equation corresponding to a circuit state variable and the other corresponds to the internal variable for the memory element.
In order to reduce the complexity of solving (3), we integrate directly (1) on both sides and then the 3D equations can be obtained: where  0 is the initial charge of the memcapacitor and  1 ,  2 ,   , and   are the time integrals of  1 ,  2 , V  , and   , respectively.For convenience, let  Figure 11(a) shows a complex bifurcation process, which has typical period 3 window, and also possesses novel period 1 and period 2 windows.Furthermore, there exist perioddoubling and inverse period-doubling bifurcations in both sides of the period 2 window.Figures 12(a)-12(d) give the typical attractors in the - plane with  0 = 0.01,  0 = 0.05,  0 = 0.11, and  0 = 0.29.

Analog Circuit Experiment
An experimental circuit is designed to realize the memcapacitor-based chaotic oscillator.The designed circuit for realizing (12) is shown in Figure 13, and its corresponding state equations are written as Note that the circuit shown in Figure 13 and used in the experiment is designed based on (12), which is a set of dimension reduction (three-dimensional) equations which is obtained from five-dimensional equation (3).Compared with other memcapacitor-based chaotic oscillators, the advantages of this experiment are that the experimental circuit is simple and the attractors or time-domain waveforms with respect to charge ( 1 and  2 ) through the inductor ( 1 and  2 ) can be observed directly from the experimental oscilloscope.

Conclusion
In this paper, a new chaotic oscillator based on chargecontrolled memcapacitor is designed.The dynamical characteristics with respect to the variations of initial conditions are investigated, and then the complex dynamic characteristics, including equilibrium set or an infinite number of equilibrium points, extreme multistability of the equilibrium set, and coexisting attractors, are found.Finally, chaotic attractors are captured experimentally by the analog circuit of implementing this oscillator.The most prominent feature of this system is the coexisting oscillation and coexisting attractors.It is concluded that one coexisting oscillation system can be as multi-random signal seeds to design different chaotic pseudorandom number generators, which can be used in secret communications, information encryptions, and various pseudorandom number generation.