Memory Circuit Elements: Complexity, Complex Systems, and Applications
View this Special IssueResearch Article  Open Access
Bocheng Bao, Aihuang Hu, Han Bao, Quan Xu, Mo Chen, Huagan Wu, "ThreeDimensional Memristive Hindmarsh–Rose Neuron Model with Hidden Coexisting Asymmetric Behaviors", Complexity, vol. 2018, Article ID 3872573, 11 pages, 2018. https://doi.org/10.1155/2018/3872573
ThreeDimensional Memristive Hindmarsh–Rose Neuron Model with Hidden Coexisting Asymmetric Behaviors
Abstract
Since the electrical activities of neurons are closely related to complex electrophysiological environment in neuronal system, a novel threedimensional memristive Hindmarsh–Rose (HR) neuron model is presented in this paper to describe complex dynamics of neuronal activities with electromagnetic induction. The proposed memristive HR neuron model has no equilibrium point but can show hidden dynamical behaviors of coexisting asymmetric attractors, which has not been reported in the previous references for the HR neuron model. Mathematical model based numerical simulations for hidden coexisting asymmetric attractors are performed by bifurcation analyses, phase portraits, attraction basins, and dynamical maps, which just demonstrate the occurrence of complex dynamical behaviors of electrical activities in neuron with electromagnetic induction. Additionally, circuit breadboard based experimental results well confirm the numerical simulations.
1. Introduction
In the past three decades, numerous simplified neuron models had been fantastically extended from the classical Hodgkin–Huxley model [1] to reconstruct the main dynamical characteristics of neuronal electrical activities [2–8], among which the two and threedimensional Hindmarsh–Rose (HR) neuron models are effective and available for dynamical analysis in electrical activities of biological neurons [9, 10]. In the last few years, a wide variety of the HR neuron models, such as original threedimensional HR models [10–16], extended or nonlinear feedback coupled HR models [17–20], time delayed HR models [20–22], fractionalorder HR models [23, 24], and memristor based HR models under electromagnetic radiations [9, 25–27], have been proposed and further studied by bifurcation analysis methods for understanding the dynamics of electrical activities among neurons [8]. For this reason, bifurcation analysis theory plays an essential role in describing mode transitions between spiking and bursting in the neuronal electrical activities [9–27].
Inspired by the constructing approach of the threedimensional HR neuron model [3, 24], a novel threedimensional memristive HR neuron model is presented in this paper, which could be used to better describe complex dynamical characteristics of neuronal electrical activities with electromagnetic induction or further exhibit some undiscovered complex dynamical behaviors in neuronal electrical activities. Interestingly, our proposed memristive HR neuron model has no equilibrium point, which can be classified as a particular dynamical system with hidden oscillating patterns [28–31]. Furthermore, coexisting asymmetric attractors’ behavior can also be observed in such memristive HR neuron model as well, indicating the emergence of bistability dynamics, which has been found in some specified neuron models [32–35]. However, the phenomenon of hidden coexisting asymmetric attractors has not been previously reported for the HR neuron model.
This paper is organized as follows. In Section 2, based on the brief reviews on the HR neuron model, a threedimensional memristive HR neuron model is presented, upon which hidden coexisting asymmetric attractors are numerically revealed by phase portraits and time series and its bistability dynamics are confirmed by the attraction basins related to the initial values. In Section 3, hidden coexisting asymmetric attractors’ behaviors are demonstrated by bifurcation diagrams, Lyapunov exponents, and dynamical maps, from which numerous types of coexisting asymmetric attractors are easily observed. In addition, a physical implementation circuit is fabricated and breadboard experiments are carried out to confirm the hidden coexisting asymmetric attractors in Section 4. The conclusions are summarized in Section 5.
2. ThreeDimensional Memristive HR Neuron Model
2.1. Brief Reviews on the HR Neuron Model
Through simplifying the classical Hodgkin–Huxley model [1], the twodimensional Hindmarsh–Rose (HR) neuron model was proposed by Hindmarsh and Rose [2] in 1982, which is described by two firstorder ordinary differential equations asin which two variables and are the membrane potential and recovery variable (also called spiking variable), respectively, and a term is the externally applied current. The parameters , , , and are four positive constants, which are often assumed as , , , and , respectively [2, 11–14].
To permit numerous dynamical behaviors, for example, chaotic dynamics, for the membrane potential, an extra third equation was introduced by Hindmarsh and Rose [3] in 1984 to improve the twodimensional neuron model (1), which is expressed by three firstorder ordinary differential equations aswhere the variable is the bursting variable and the constant is the resting potential of the model. The newly added parameters and are two positive constants but is very small. Thus, a new variable , a slowly evolving current, is coupled into the first equation of the twodimensional model (1) to tune the externally applied current . If the threedimensional neuron model (2) is in its firing state, the value of increases [24].
2.2. Constructed Memristive HR Neuron Model
Motivated by the above constructing approach of the model (2), a threedimensional memristive HR neuron model with electromagnetic induction effect is proposed through introducing a fluxcontrolled ideal memristor into the first equation of the model (1), which can be mathematically modeled bywhere the new variable is the magnetic flux indicating the time integral of the membrane potential . The newly adding term represents the externally applied electromagnetic induction and is the strength of the electromagnetic induction.
It is important to stress that the introduced memristor in (3) is ideal and fluxcontrolled. According to the definitions of ideal memristor by statedependent Ohm’s laws between the terminal voltage and terminal current [36, 37], a fluxcontrolled ideal memristor is thereby given aswhere the memductance can be interpreted as the fluxdependent rate of change of charge. Therefore, the memductance utilized in (3) can be written byin which the coefficient is positive.
To exhibit three characteristic fingerprints of pinched hysteresis loop of the memristor modeled by (4) and (5) [38], a sinusoidal voltage source is connected at the input terminals of the memristor, where and are the amplitude and frequency, respectively. Let . When is maintained unchanged and is assigned as 0.1, 0.2, and 0.5, respectively, the plots are shown in Figure 1(a), while when is fixed and is determined as 3, 4, and 5, respectively, the plots are shown in Figure 1(b). It is seen from Figure 1 that the plots are the hysteresis loops pinched at the origin. The hysteresis loop is pinched regardless of the stimulus amplitudes but shrinks into a linear function at infinite frequency and its lobe area decreases with increasing the frequency. The numerical results in Figure 1 indicate that the memristor modeled by (4) can behave three fingerprints for distinguishing memristors [38].
(a)
(b)
In the next work, the threedimensional memristive HR neuron model given in (3) is considered. It should be remarked that the adjustable parameters of interest are and , and their regions are correspondent to the first quadrant of the parameter space ( and ). For any uncertain parameter , the existence of any equilibrium point is not allowed in the threedimensional memristive HR neuron model, neither stable nor unstable. Only if the applied current will the model show an equilibrium point, which is not in the considered parameter region. This case is often encountered in various kinds of nonlinear dynamical systems that are known to generate the specified hidden attractors [28–31].
2.3. Coexisting Asymmetric Attractors
When the original parameters are selected as , , , and , respectively [2], an example for model (3) with and is given as shown in Figure 2, where the orbits marked by the red and blue colors emerge from the initial values (0, 0, −2) and (0, 0, 2), respectively. In Figure 2(a), the phase portraits in the plane display the bistability phenomenon of hidden coexisting asymmetric attractors consisting of chaotic attractor and limit cycle in the memristive HR neuron model, whereas in Figure 2(b), the time series of the membrane potential demonstrate the coexistence of chaotic and periodic spikes in the memristive HR neuron model as well. Correspondingly, three Lyapunov exponents for the initial values (0, 0, −2) are 0.0782, 0, and −3.0684, respectively, while those for (0, 0, 2) are 0, −0.2717, and −2.8556, respectively. Remark that Wolf et al.’s method [39] with MATLAB ODE113 algorithm is here used to calculate three Lyapunov exponents.
(a)
(b)
For the coexisting asymmetric attractors shown in Figure 2(a), the corresponding attraction basins in the and planes of the initial values are drawn in Figures 3(a) and 3(b), where the attraction basins for chaotic attractors and periodic limit cycles are colored in the fuchsia and cyan regions, respectively. The results effectively indicate the emergence of bistability phenomenon in the memristive HR neuron model.
(a)
(b)
Particularly, the emerging coexisting asymmetric attractors do not associate with any equilibrium point, indicating that the memristive HR neuron model always operates in hidden oscillating patterns [28–31]. Additionally, it is interesting to note that, just like the selfexcited coexisting asymmetric attractors in hyperbolictype memristor based Hopfield neural network [32], such hidden coexisting asymmetric attractors in the memristive HR neuron model are induced by electromagnetic induction also, which illustrates the occurrence of complex dynamical behaviors of electrical activities in neuron with electromagnetic induction.
3. Hidden Coexisting Asymmetric Attractors’ Behavior
When the applied current and electromagnetic induction strength are considered as two bifurcation parameters, hidden coexisting asymmetric behaviors of the memristive HR neuron model are numerically studied by MATLAB ODE45 algorithm under two sets of the initial values (0, 0, −2) and (0, 0, 2).
3.1. Bifurcation Behaviors with Increasing
Figure 4 gives the bifurcation diagrams of and the first two Lyapunov exponents as = 1 and = 0.5~1.4, where in Figure 4(a) the orbits marked by the red and blue colors emerge from the initial values (0, 0, −2) and (0, 0, 2), respectively, and in Figure 4(b) the Lyapunov exponents marked by the red and fuchsia colors correspond to the initial values (0, 0, −2) and those marked by blue and green colors correspond to (0, 0, 2). It can be seen from Figure 4 that complex dynamics appear in the memristive HR neuron model, in which chaotic attractors with different topologies, limit cycles with different periodicities, perioddoubling bifurcation routes, tangent bifurcation routes, crisis scenarios, coexisting bifurcation modes, and so on can be found. Therefore, the electromagnetic induction by the introduced memristor induces numerous complex dynamics for the membrane potential , especially including hidden coexisting asymmetric behaviors.
(a)
(b)
The concernedly coexisting behaviors of asymmetric attractors mainly locate in two parameter regions and , in which some different types of hidden coexisting asymmetric attractors occur. When and and 0.81, respectively, the phase portraits in the plane for the other two types of hidden coexisting asymmetric attractors are depicted in Figure 5. In detail, Figure 5(a) exhibits the coexistence of hidden chaotic attractor and hidden limit cycle, and Figure 5(b) displays the coexistence of two hidden limit cycles with different periodicities.
(a)
(b)
3.2. Bifurcation Behaviors with Increasing
Figure 6 demonstrates the bifurcation diagrams of and the first two Lyapunov exponents as = 0~2.4 and = 0.9, where the initial values for different colored orbits in Figure 6(a) and different colored Lyapunov exponents in Figure 6(b) are consistent with those used in Figures 4(a) and 4(b). In the same way, it can be observed from Figure 6 that complex dynamics are coined in the memristive HR neuron model, reflecting the dynamical effect of the externally applied current in neurons.
(a)
(b)
The parameter region has the benefit for the coexisting behaviors of asymmetric attractors, in which several different types of hidden coexisting asymmetric attractors can be clearly found. When = 0.9 is fixed and is set to 1.15 and to 1.62, respectively, the phase portraits in the plane for the two types of hidden coexisting asymmetric attractors are plotted in Figure 7, where in Figure 7(a) the coexistence of hidden chaotic attractor and hidden period 1 limit cycle is displayed, and in Figure 7(b) the coexistence of hidden period 2 limit cycle and hidden chaotic attractor with large size is presented.
(a)
(b)
3.3. Coexisting Asymmetric Behaviors in the Parameter Space
For intuitively manifesting the coexisting behaviors of asymmetric attractors in the memristive HR neuron model, dynamical maps depicted by the largest Lyapunov exponent under two sets of the initial values are numerically plotted in the parameter space [40], as shown in Figures 8(a) and 8(b), where the luminous yellow, red, and black colored regions stand for the chaotic, periodic, and divergent behaviors, respectively. Figure 8 indicates how dynamical evolution in the electromagnetic induction strength and applied current affects the coexisting behaviors under different initial values being considered. When the two parameters and are evolved, some chaotic regions are embedded in the periodic regions; however, different chaotic regions appear on the parameter spaces of Figures 8(a) and 8(b), which are caused by the coexisting asymmetric attractors’ behaviors under different initial values. The numerical results in Figure 8 illustrate that the dynamical behaviors depicted by the largest Lyapunov exponent based dynamical maps are well agreed with those revealed by the bifurcation behaviors in Figures 4 and 6.
(a)
(b)
It should be remarkable for the dynamical maps in Figure 8 that the chaotic attractors in different locations of the chaotic region have different topologies and the limit cycles in different locations of the periodic region have different periodicities. Specifically, except for several types of coexisting asymmetric behaviors shown in Figures 2, 5, and 7, another type of coexisting asymmetric behaviors of chaotic attractor and divergent orbit can also be uncovered, which means that another form of bistability exists in the memristive HR neuron model.
4. Circuit Design and Breadboard Experiments
4.1. Physical Circuit Designs and Parameter Selections
The fluxcontrolled ideal memristor characterized by (4) and (5) and its constructing threedimensional memristive HR neuron model expressed by (3) can be physically realized by using an electronic circuit via analog multipliers and operational amplifiers connected with resistors and/or capacitors [41–43], as drawn in Figures 9(a) and 9(b), respectively. Of course, this threedimensional memristive HR neuron model also can be digitally implemented in fieldprogrammable gate arrays (FPGA) as well [44, 45].
(a)
(b)
The implementation circuit of the fluxcontrolled ideal memristor in Figure 9(a) contains an integrator with time constant , an inverter, a multiplier , and a resistor . For the input voltage and output current , the mathematical model for the memristor emulator can be easily given aswhere is the inner variable of the memristor emulator and is the gain of the multipliers , = , and .
The main circuit of the memristive HR neuron model in Figure 9(b) has two integrating channels for implementing the first and second equations of (3). According to Kirchhoff’s circuit laws and electrical properties of the circuit components, the circuit equations of Figure 9(b) are written aswhere and are two circuit variables, and are two applied voltages, and , , and are the gains of the multipliers , , and , respectively.
Considering that the dynamic amplitude of the recovery variable in the numerical simulations exceeds the linear operation ranges of operational amplifier and multiplier, the following linear transformationshould be utilized to reduce the dynamic voltage amplitude of in the circuit equations of the memristive HR neuron model. Thus, by comparing (8) with (3), there yieldsLet the time constant = 10 kΩ × 33 nF = 330 μs; that is, = 10 kΩ and = 33 nF and the multiplier gains = = 0.1 and = = 1. According to (9), the circuit parameters for the breadboard experiments of the memristive HR neuron model can be calculated, as listed in Table 1.

4.2. Results Captured from Breadboard Experiments
According to the circuit diagrams in Figure 9 and circuit parameters in Table 1, a hardware circuit using commercially discrete components can be welded on a breadboard. The operational amplifiers AD711JN and analog multipliers AD633JN supplied by ±15 V voltage modules are chosen. The DC voltages and are provided by Tektronix PWS 2326 DC Power Supply and the experimental results are measured by Tektronix TDS 3054C Digital Phosphor Oscilloscope. The photograph of the connectedly experimental prototype for the memristive HR neuron model is displayed in Figure 10.
For experimentally measuring the pinched hysteresis loops of the memristor emulator given in Figure 9(a), a sinusoidal voltage source generated by Tektronix AFG 3102C Function Generator is linked to the input terminals of the memristor emulator, where the physical frequency is calculated by . For = 1, the adjustable resistance = 10 kΩ. When the amplitudes and frequencies of the sinusoidal voltage source used during numerical simulations in Figure 1 are employed, the pinched hysteresis loops for the corresponding amplitudes and frequencies are captured, as shown in Figures 11(a) and 11(b), respectively, which experimentally validate the characteristic fingerprints of the memristor emulator. It should be addressed that for better observing the experimental results, all the output currents sensed by the current probe are magnified by ten times.
(a)
(b)
The circuit parameters listed in Table 1 are used and the different initial voltages of three capacitors are randomly sensed by repeatedly switching on and off the experimental power supply [46]. For the typical circuit parameters in Table 1, two adjustable circuit parameters of and correspond to the adjustable model parameters of = 1 and = 0.9. Corresponding to Figure 2, the phase portraits in the plane and time series of the variable that emerged from different initial voltages are experimentally obtained, as shown in Figure 12. The experimental results indicate that hidden coexisting asymmetric attractors also can be measured from the breadboard experiments of the memristive HR neuron model as well.
(a)
(b)
When the applied current = 1, that is, the resistance is fixed as 10 kΩ and the resistance is set to 1.36 kΩ and to 1.23 kΩ, respectively, the phase portraits in the plane are captured, as shown in Figures 13(a) and 13(b). Furthermore, when the electromagnetic induction strength = 0.9, that is, = 1.11 kΩ and is set to 8.70 kΩ and to 6.17 kΩ, respectively, the phase portraits in the plane are captured, as shown in Figures 13(c) and 13(d). Ignoring some tiny differences between numerical simulations and breadboard experiments due to the computational errors and parasitic circuit parameters, the experimental results are almost the same as the numerical simulations, which imply that the coexisting asymmetric attractors’ behaviors that emerged from the memristive HR neuron model can be validated experimentally.
(a)
(b)
(c)
(d)
5. Conclusions
This paper presents a novel threedimensional memristive HR neuron model to describe complex dynamics of neuronal activities with electromagnetic induction. The most prominent feature of this neuron model is that it does not contain any equilibrium point but can exhibit hidden coexisting behaviors of asymmetric attractors. Through executing bifurcation analyses, phase portraits, attraction basins, and dynamical maps, hidden coexisting asymmetric attractors are uncovered from the mathematical model and verified from the corresponding breadboard experiments. Thus, the proposed memristive HR neuron model can imitate the complex dynamical behaviors of electrical activities in neuron with electromagnetic induction. Further investigations will be performed in our future works.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the grants from the National Natural Science Foundations of China under Grant nos. 51777016, 51607013, 11602035, and 61601062 and the Natural Science Foundations of Jiangsu Province, China, under Grant no. BK20160282.
References
 A. L. Hodgkin and A. F. Huxley, “A quantitative description of membrane current and its application to conduction and excitation in nerve.,” The Journal of Physiology, vol. 117, no. 4, pp. 500–544, 1952. View at: Google Scholar
 J. L. Hindmarsh and R. M. Rose, “A model of the nerve impulse using two firstorder differential equations,” Nature, vol. 296, no. 5853, pp. 162–164, 1982. View at: Publisher Site  Google Scholar
 J. L. Hindmarsh and R. M. Rose, “A model of neuronal bursting using three coupled first order differential equations,” Proceedings of the Royal Society of London B: Biological Sciences, vol. 221, no. 1222, pp. 87–102, 1984. View at: Publisher Site  Google Scholar
 C. R. Laing and C. C. Chow, “A spiking neuron model for binocular rivalry,” Journal of Computational Neuroscience, vol. 12, no. 1, pp. 39–53, 2002. View at: Publisher Site  Google Scholar
 K. Tsumoto, H. Kitajima, T. Yoshinaga, K. Aihara, and H. Kawakami, “Bifurcations in MorrisLecar neuron model,” Neurocomputing, vol. 69, no. 46, pp. 293–316, 2006. View at: Publisher Site  Google Scholar
 H. Gu and B. Pan, “A fourdimensional neuronal model to describe the complex nonlinear dynamics observed in the firing patterns of a sciatic nerve chronic constriction injury model,” Nonlinear Dynamics, vol. 81, no. 4, pp. 2107–2126, 2015. View at: Publisher Site  Google Scholar
 F. Wu, C. Wang, and Y. Xu, “Model of electrical activity in cardiac tissue under electromagnetic induction,” Scientific Reports, vol. 6, article 28, 2016. View at: Publisher Site  Google Scholar
 J. Ma and J. Tang, “A review for dynamics in neuron and neuronal network,” Nonlinear Dynamics, vol. 89, no. 3, pp. 1569–1578, 2017. View at: Google Scholar  MathSciNet
 M. Ge, Y. Jia, Y. Xu, and L. Yang, “Mode transition in electrical activities of neuron driven by high and low frequency stimulus in the presence of electromagnetic induction and radiation,” Nonlinear Dynamics, vol. 91, no. 1, pp. 515–523, 2018. View at: Google Scholar  MathSciNet
 M. Storace, D. Linaro, and E. de Lange, “The HindmarshRose neuron model: bifurcation analysis and piecewiselinear approximations,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 18, no. 3, Article ID 033128, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 G. Innocenti, A. Morelli, R. Genesio, and A. Torcini, “Dynamical phases of the HindmarshRose neuronal model: studies of the transition from bursting to spiking chaos,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 17, no. 4, Article ID 043128, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 J. M. GonzálezMiranda, “Complex bifurcation structures in the HindmarshRose neuron model,” International Journal of Bifurcation and Chaos, vol. 17, no. 9, pp. 3071–3083, 2007. View at: Google Scholar
 J. M. GonzálezMiranda, “Observation of a continuous interior crisis in the HindmarshRose neuron model,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 13, no. 3, pp. 845–852, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 G. Innocenti and R. Genesio, “On the dynamics of chaotic spikingbursting transition in the HindmarshRose neuron,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 19, no. 2, Article ID 023124, 2009. View at: Google Scholar  MathSciNet
 T. Wei, B. Ou, J. Li et al., “Transcriptional profiling of rice early response to Magnaporthe oryzae identified OsWRKYs as important regulators in rice blast resistance,” PLoS ONE, vol. 8, no. 3, Article ID e59720, 2013. View at: Publisher Site  Google Scholar
 S. R. D. Djeundam, R. Yamapi, T. C. Kofane, and M. A. Azizalaoui, “Deterministic and stochastic bifurcations in the HindmarshRose neuronal model,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 23, no. 3, Article ID 033125, 2013. View at: Google Scholar  MathSciNet
 E. B. Ngouonkadi, H. B. Fotsin, P. Louodop Fotso, V. Kamdoum Tamba, and H. A. Cerdeira, “Bifurcations and multistability in the extended HindmarshRose neuronal oscillator,” Chaos, Solitons & Fractals, vol. 85, pp. 151–163, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 K. Wu, T. Luo, H. Lu, and Y. Wang, “Bifurcation study of neuron firing activity of the modified Hindmarsh–Rose model,” Neural Computing and Applications, vol. 27, no. 3, pp. 739–747, 2016. View at: Publisher Site  Google Scholar
 S. K. Thottil and R. P. Ignatius, “Nonlinear feedback coupling in Hindmarsh–Rose neurons,” Nonlinear Dynamics, vol. 87, no. 3, pp. 1879–1899, 2017. View at: Publisher Site  Google Scholar
 H. Wang, Y. Zheng, and Q. Lu, “Stability and bifurcation analysis in the coupled HR neurons with delayed synaptic connection,” Nonlinear Dynamics, vol. 88, no. 3, pp. 2091–2100, 2017. View at: Google Scholar  MathSciNet
 S. Lakshmanan, C. P. Lim, S. Nahavandi, M. Prakash, and P. Balasubramaniam, “Dynamical analysis of the HindmarshRose neuron with time delays,” IEEE Transactions on Neural Networks and Learning Systems, vol. 28, no. 8, pp. 1953–1958, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 H. Wang, Q. Wang, and Y. Zheng, “Bifurcation analysis for HindmarshRose neuronal model with timedelayed feedback control and application to chaos control,” Science China Technological Sciences, vol. 57, no. 5, pp. 872–878, 2014. View at: Publisher Site  Google Scholar
 D. Jun, Z. GuangJun, X. Yong, Y. Hong, and W. Jue, “Dynamic behavior analysis of fractionalorder HindmarshRose neuronal model,” Cognitive Neurodynamics, vol. 8, no. 2, pp. 167–175, 2014. View at: Publisher Site  Google Scholar
 E. Kaslik, “Analysis of two and threedimensional fractionalorder HindmarshRose type neuronal models,” Fractional Calculus and Applied Analysis, vol. 20, no. 3, pp. 623–645, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 M. Lv, C. Wang, G. Ren, J. Ma, and X. Song, “Model of electrical activity in a neuron under magnetic flow effect,” Nonlinear Dynamics, vol. 85, no. 3, pp. 1479–1490, 2016. View at: Publisher Site  Google Scholar
 M. Lv and J. Ma, “Multiple modes of electrical activities in a new neuron model under electromagnetic radiation,” Neurocomputing, vol. 205, pp. 375–381, 2016. View at: Publisher Site  Google Scholar
 L. Lu, Y. Jia, W. Liu, and L. Yang, “Mixed stimulusinduced mode selection in neural activity driven by high and low frequency current under electromagnetic radiation,” Complexity, vol. 2017, Article ID 7628537, pp. 1–11, 2017. View at: Google Scholar  MathSciNet
 S. Jafari, J. C. Sprott, and S. M. R. H. Golpayegani, “Elementary quadratic chaotic flows with no equilibria,” Physics Letters A, vol. 377, no. 9, pp. 699–702, 2013. View at: Publisher Site  Google Scholar
 V.T. Pham, C. Volos, S. Jafari, Z. Wei, and X. Wang, “Constructing a novel noequilibrium chaotic system,” International Journal of Bifurcation and Chaos, vol. 24, no. 5, Article ID 1450073, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 F. R. Tahir, S. Jafari, V.T. Pham, C. Volos, and X. Wang, “A novel noequilibrium chaotic system with multiwing butterfly attractors,” International Journal of Bifurcation and Chaos, vol. 25, no. 4, Article ID 1550056, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 A. P. Kuznetsov, S. P. Kuznetsov, E. Mosekilde, and N. V. Stankevich, “Coexisting hidden attractors in a radiophysical oscillator system,” Journal of Physics A: Mathematical and Theoretical, vol. 48, no. 12, Article ID 125101, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 B. Bao, H. Qian, Q. Xu, M. Chen, J. Wang, and Y. Yu, “Coexisting behaviors of asymmetric attractors in hyperbolictype memristor based hopfield neural network,” Frontiers in Computational Neuroscience, vol. 11, article 81, 2017. View at: Publisher Site  Google Scholar
 A. Shilnikov, R. L. Calabrese, and G. Cymbalyuk, “Mechanism of bistability: tonic spiking and bursting in a neuron model,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 71, no. 2, Article ID 056214, 2005. View at: Google Scholar  MathSciNet
 B. Bao, H. Qian, J. Wang et al., “Numerical analyses and experimental validations of coexisting multiple attractors in Hopfield neural network,” Nonlinear Dynamics, vol. 90, no. 4, pp. 2359–2369, 2017. View at: Publisher Site  Google Scholar
 T. Malashchenko, A. Shilnikov, and G. Cymbalyuk, “Six types of multistability in a Neuronal model based on slow calcium current,” PLoS ONE, vol. 6, no. 7, Article ID e21782, 2011. View at: Publisher Site  Google Scholar
 L. Chua, “Resistance switching memories are memristors,” Applied Physics A: Materials Science & Processing, vol. 102, no. 4, pp. 765–783, 2011. View at: Publisher Site  Google Scholar
 L. O. Chua, “The fourth element,” Proceedings of the IEEE, vol. 100, no. 6, pp. 1920–1927, 2012. View at: Publisher Site  Google Scholar
 S. P. Adhikari, M. P. Sah, H. Kim, and L. O. Chua, “Three fingerprints of memristor,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 60, no. 11, pp. 3008–3021, 2013. View at: Publisher Site  Google Scholar
 A. Wolf, J. B. Swift, and H. L. a. Swinney, “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena, vol. 16, no. 3, pp. 285–317, 1985. View at: Publisher Site  Google Scholar  MathSciNet
 P. C. Rech, “Periodadding and spiral organization of the periodicity in a Hopfield neural network,” International Journal of Machine Learning and Cybernetics, vol. 6, no. 1, 2013. View at: Publisher Site  Google Scholar
 H. Bao, N. Wang, B. Bao, M. Chen, P. Jin, and G. Wang, “Initial conditiondependent dynamics and transient period in memristorbased hypogenetic jerk system with four line equilibria,” Communications in Nonlinear Science and Numerical Simulation, vol. 57, pp. 264–275, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 B. C. Bao, H. Bao, N. Wang, M. Chen, and Q. Xu, “Hidden extreme multistability in memristive hyperchaotic system,” Chaos, Solitons & Fractals, vol. 94, pp. 102–111, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 K. Rajagopal, A. Akgul, S. Jafari, and B. Aricioglu, “A chaotic memcapacitor oscillator with two unstable equilibriums and its fractional form with engineering applications,” Nonlinear Dynamics, vol. 91, no. 2, pp. 957–974, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 K. Rajagopal, L. Guessas, A. Karthikeyan, A. Srinivasan, and G. Adam, “Fractional order memristor no equilibrium chaotic system with its adaptive sliding mode synchronization and genetically optimized fractional order PID synchronization,” Complexity, vol. 2017, Article ID 1892618, pp. 1–19, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 K. Rajagopal, A. Karthikeyan, and A. K. Srinivasan, “FPGA implementation of novel fractionalorder chaotic systems with two equilibriums and no equilibrium and its adaptive sliding mode synchronization,” Nonlinear Dynamics, vol. 87, no. 4, pp. 2281–2304, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 B. Bao, T. Jiang, G. Wang, P. Jin, H. Bao, and M. Chen, “Twomemristorbased Chua’s hyperchaotic circuit with plane equilibrium and its extreme multistability,” Nonlinear Dynamics, vol. 89, no. 2, pp. 1157–1171, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Bocheng Bao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.