Diverse Dynamic Behaviors and Firing Activities of the Modified Fractional-Order Hindmarsh–Rose Neuronal Model Induced by Fractional-Order

Yang, Xin; Zhang, GuangJun; Li, XueRen; Wang, Dong

doi:https://doi.org/10.1155/2021/8490695

Complexity

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest References Copyright Related Articles

Special Issue

Fractional-Order Chaotic Systems: Theory and Applications

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 8490695 | https://doi.org/10.1155/2021/8490695

Diverse Dynamic Behaviors and Firing Activities of the Modified Fractional-Order Hindmarsh–Rose Neuronal Model Induced by Fractional-Order

Xin Yang,¹GuangJun Zhang,²XueRen Li,¹and Dong Wang¹

Academic Editor: Guillermo Huerta Cuellar

Received28 Jun 2021

Revised09 Aug 2021

Accepted24 Oct 2021

Published12 Nov 2021

Abstract

It is important to investigate the firing activities of neurons, and previous experimental works have shown that fractional-order neuronal models depict the firing rate of neurons more verifiably. In this study, a modified fractional-order Hindmarsh–Rose neuronal model is proposed, and the dynamics and firing activities are investigated. Some novel phenomenon can be found. First, by analyzing numerically and theoretically, the Hopf bifurcation is found to occur when the external direct current stimulus is chosen appropriately. The effects of fractional-order on the bifurcation are also studied. Second, when injecting a direct current stimulus, compared with the integer-order model, the system has more varying dynamic behaviors and firing pattern transitions. Under different external current stimulus, periodic firing patterns and chaotic firing patterns occur when fractional-order changes, but the regions of chaotic firing patterns are different. In other words, the transition mode of periodic firing and chaotic firing induced by fractional-order is different under different external current stimulus. The two-dimensional colored diagram of firing patterns is also investigated. Finally, when injecting periodic current stimulus, regular/irregular bursting, multiple spiking, regular\irregular square wave bursting, and mixed firing mods are found by setting the appropriate fractional-order, amplitude, and frequency of the external current stimulus. Some firing patterns cannot be found in integer-order models. When the amplitude is chosen at appropriate values, the region of frequency when the system displays the mixed firing modes decreases with increasing fractional-order.

1. Introduction

The neuron was modeled and analyzed by a large amount of experimental data, and the results show that the firing behaviors of neurons are nonlinear processes. In 1952, Hodgkin and Huxley used equivalent circuits and large amounts of data from experiments to model and analyze the data. They constructed the Hodgkin–Huxley (HH) neuron model through theoretical derivation [1]. In 1961, FitzHugh simplified the variables in the HH model and constructed a low-dimensional model, the two-dimensional FitzHugh–Nagumo (FHN) model [2]. Morris and Lecar summarized the new neuron model (Morris–Lecar (ML) model) from the experimental data of the arctic goose muscle fiber, which is a further simplification of the HH model. In 1982, based on voltage clamp experiment data of snail nerve cells, Hindmarsh and Rose proposed the Hindmarsh–Rose (HR) model [3]. The modified HR model is called the polynomial model, which is the modified HR-type model [4]. This is also an important model for neurons, but research on this model’s firing activities is limited. The modified HR model has many dynamic characteristics and depicts bursting and spiking behaviors successfully. It is significant to investigate the dynamic behavior of the modified HR model. Firing activities are important phenomena in the neuronal system. A number of researchers have used different neuronal models to try to explain some of the phenomena observed in experiments. Studying the motor mechanism of the neuronal dynamical system helps us to understand relevant phenomena in the brain and contribute to the development of artificial intelligence.

Neuronal models have been studied numerically and theoretically [5, 6]. In [7], the HR neuronal model was used to develop a new type of a central pattern generator. Neuronal models describe the firing patterns in real neurons, for example [8]. For neuronal model firing, the effects of electromagnetic radiation [9, 10] and extracellular alternating-current feld [11] have been investigated. In [9], chaotic bursting and periodic firing were found in the HR model under an electric field, and Hopf bifurcation was verified in this system. In [10], the mode transition of electrical activities induced by electromagnetic radiation was studied. The two-dimensional discrete HR model was investigated in [12]. Diverse bifurcation points and dynamic behaviors were found. The neuronal system’s synchronization behaviors were investigated in [9, 13]. For collective behaviors in neuronal networks, [14] depicted the dynamics of the neurons and neuronal network, such as the mode transition in electrical activity, pattern formation, and selection in networks of neurons. In recent years, some researchers have investigated modified HR models by analyzing bifurcation points theoretically and firing pattern transitions [15–17]. Reference [15] analyzed the dynamic behaviors of a modified HR model with induced electromagnetic radiation and injected current stimulus, and diverse bifurcation points and firing activities were found. The firing pattern transitions with direct current stimulus and without current stimulus were investigated in [16, 17], respectively. At present, there are few studies on the modified HR model, but this model’s dynamic behaviors can be important to understand the experimental results.

The above studies studied the integer-order model only. In results of [18], after analyzing the dynamics of the firing rate with a range of stimulus dynamics, showed that the multiple time scale adaptation is consistent with fractional-order differentiation. Fractional computation is related to the intrinsic properties of dynamical systems [19]. Local memory, nonuniform diffusion effects, and many other aspects require fractional computation [19, 20]. Fractional-order differentiation is a fundamental and general computation that can contribute to efficient information processing, stimulus anticipation, and assessment of frequency-independent phase shifts of oscillatory neuronal firing. In this way, fractional-order differentiation has advantages over integer-order differentiation in the depiction of the firing characteristics of some types of neurons. The fractional-order model can give a more complete picture of neurons’ dynamic characteristics than the integer-order model. In the past, scholars have performed much research on fraction-order dynamical systems and applied them in many fields, such as financial systems [21], biomedicine [22, 23], and the spread of infectious diseases [24]. In terms of the fractional-order neuronal model, the single fractional-order HR neuronal model [25] was investigated, and the transitions of chaotic firing and periodical firing, spike firing, and bursting firing were observed. Other fractional-order models such as the fractional-order Lorenz [26, 27], piecewise-linear (PWL) hyperchaotic system [28], and Rossler system [27, 29] were also investigated, and diverse dynamic behaviors were found. The firing activities of neuronal models need to be further investigated, and modified fractional-order HR model may have more diverse firing patterns. However, there is no research studying the modified fractional-order HR model under external direct and periodic current stimulus. Like [21, 23], the Adomian decomposition methods (ADM) algorithm is used in the present study because of its great computational efficiency.

From this perspective, this study is organized as follows: first, the modified fractional-order HR model under an external current stimulus is proposed. The bifurcation analysis and influence of fractional-order on the bifurcation are studied numerically and theoretically. Then, the firing patterns, firing pattern transition, and bifurcation diagrams with a bifurcation parameter of fractional-order when injected with a direct current stimulus were investigated. Finally, when a periodic current stimulus is injected, diverse firing activities are studied.

2. Model Description and Bifurcation Analysis

2.1. Model Description

After Hindmarsh and Rose proposed the Hindmarsh–Rose model [3], a large number of studies have been constructed on the HR neuronal model. The HR neuronal model is described as follows:where is the membrane action potential, is a recovery variable, is a slow adaption current, and is the external stimulus current. are the parameters in the system, and they are usually set as . From previous studies [15–17], the integer-order modified HR model under a current stimulus is described as follows:where represents the membrane voltage of the model; and measure the gating dynamics of the ion channels and the dynamics of the cytosolic channels, respectively; are the parameters of the modified HR model, and in this study, they are set as . is a variable.

There are many definitions of fractional derivatives, and in practice, the most frequently used definitions are three: the Grunwald–Letnikov derivative, Riemann–Liouville derivative, and Caputo derivative. According to [25], under some conditions, the three definitions are equivalent and can be intertranslated. The Caputo derivative makes the Laplace transformation more concise, so it is simpler to solve the fractional-order derivative. Similar to most literature, the Caputo derivative is adopted in this study.

Definition 1. The Caputo derivative of function is defined aswhere , is the gamma function which is defined asEspecially, when ,The fractional-order modified HR neuronal model can be proposed aswhere is the differential operator defined by Caputo, is the fractional-order operator, and the other variables and parameters are the same as in the integer-order model. In this study, similar to the integer-order model, the parameters are set as follows: . and are the variables.

2.2. Bifurcation Analysis

The equilibrium point can be obtained by solving the equation system , which is

In this part, similar to the integer-order model in the above part, is variable, and the other parameters are set as . The equilibrium point can be described as follows after calculation:

From equation system (7) and the equilibrium point , the equation with one unknown quantity is as follows:

To solve the equation, the Cardan formula [5] is introduced. Now, the conversions are , , , , . Equation (9) can be simplified as follows:where , , and the Cardan discriminant is . After calculation, the Cardan discriminant , so there exists a single real root which is defined by

Now, we consider the existence of Hopf bifurcation. Suppose that the system has an equilibrium at which the following properties are satisfied:(1) has a simple pair of purely imaginary eigenvalues and other eigenvalues with negative real parts(2)

Then, the system has a Hopf bifurcation at the equilibrium . The eigenvalues can be analyzed by determining the Jacobian matrix associated with system (7), and the characteristic equation can be obtained as equation (12), where is an identity matrix of the same size as .in which

Using the same method as in equation (9), equation (9) can be simplified as follows:in which , . The three roots of equation (12) are

From the above analysis, when , there exists a simple pair of purely imaginary eigenvalues and a negative real root at equilibrium . Then, the results can be calculated as

The conclusion that there exists a Hopf bifurcation when can be drawn. For the fractional-order system, the stability can be studied using the following lemma:

Lemma 1. Consider the linear autonomous system , , where and . The system is asymptotically stable when and only when is satisfied for an arbitrary eigenvalue . The system is stable when and only when is satisfied for an arbitrary eigenvalue . The stable area of the q-order linear system is shown in Figure 1.

According to Lemma 1, the critical values of the external current stimulus for the fractional-order model to start firing at different fractional-orders can be calculated, and the results are shown in Figure 2.

2.3. Firing Behavior under Different Direct Current Stimulus and Fractional-Orders

In this section, the parameters are also set as . and are the variables. In [17], the firing activity of integer-order and the influence of parameters on the firing patterns were studied. In this section, the influence of fractional-order and direct current stimulus are investigated. In [17], when , the integer-order modified HR model shows a spiking firing pattern. Figure 3 shows the bifurcation diagram varying with the fractional-order when . Figure 3 shows that the firing patterns change with varying fractional-order.

Figure 4 shows the corresponding time series of and the phase diagram of . When , as shown in Figures 4(a) and 4(b), the neuronal model displays periodic-6 bursting. From Figures 4(c) and 4(d), the neuronal model displays periodic-5 bursting when q = 0.25. When q = 0.3 and q = 0.39, the neuronal model displays chaotic bursting (Figures 4(e), 4(f), 4(i), and 4(j)). The neuronal model displays periodic-4 bursting when q = 0.35 (Figures 4(g) and 4(h)), periodic-2 bursting when q = 0.415 (Figures 4(k) and 4(l)), and spiking when q = 0.55 (Figures 4(m) and 4(n)). From the above analysis, under the influence of the fractional-order, more complex dynamic behaviors occur. The fractional-order can induce a periodic firing pattern for the integer-order model transiting to the chaotic firing pattern.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

To prove that if fractional-order can induce a chaotic firing pattern for an integer-order model transitioning to a periodic firing pattern, the firing pattern varies with fractional-order when is set as an example. In [17], the integer-order neuronal model displays chaotic firing. As shown in Figure 5 (the bifurcation diagram when ), periodic bursting occurs when the fractional-order decreases. There exists a periodic window around .

When , the corresponding time series of and the phase diagram of are shown in Figure 6. When , as shown in Figures 6(a) and 6(b), the neuronal model displays chaotic firing. When q = 0.4, q = 0.3, q = 0.2, the neuronal model displays periodic-3, periodic-4, and periodic-5 bursting (Figures 6(c)–6(h)). Meanwhile, fractional-order can induce the chaotic firing pattern to transition to a periodic firing pattern.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

From Figure 7, the two-dimensional colored diagram of firing patterns, we can easily find the chaotic firing range and periodic firing range. We find that the chaotic range is larger when the direct current stimulus is smaller.

2.4. Firing Behavior under External Periodic Current Stimulus

In the above analysis, various firing activities are investigated when the modified fractional-order HR model is under direct external current. When the direct external current is replaced by a periodic external current stimulus, there are various firing activities. New firing patterns are found when the fractional-order changes. Spiking, bursting, and other firing patterns can occur at an appropriate fractional-order and by injecting a periodic current stimulus. The modified fractional-order HR model under periodic external current stimulus is as follows:where is the amplitude of the external current stimulus and is the angular frequency. In this section, similar to the above section, the parameters are set as , , , and the fractional-order , the amplitude , and angular frequency are considered as variables.

We set , and regular bursting and irregular bursting can be found. When , the system is in regular bursting, as shown in Figure 8(a). Irregular bursting can be found when (Figure 8(b)). After further investigating the firing activities, the critical value of fractional-order exists between regular bursting and irregular bursting. When the frequency is changed, we can also find irregular bursting and regular bursting. For example, when , the system is in regular bursting when (Figure 8(c)) and irregular bursting when (Figure 8(d)).

(a)

(b)

(c)

(d)

The amplitude of the external current stimulus is increased as . Multiple spiking consisting of different spiking frequencies, spiking, and square wave bursting can be found. Set , and the system is in regular bursting for (Figure 9(a)), multiple spiking when (Figure 9(b)), and square wave bursting when (Figure 9(c)). Setting , the system displays mixed firing modes that include bursting and spiking when (Figure 9(d)), irregular bursting when (Figure 9(e)), and another irregular bursting that is different from the above bursting activities when (Figure 9(f)). Furthermore, setting , the system displays regular spiking when (Figure 9(g)), mixed mode oscillations that consist of multiple numbers of spiking and subthreshold oscillations when (Figure 9(h)), and the system can also display regular bursting and mixed firing modes when (Figure 9(i)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

When , the firing activities are similar to the above firing patterns. However, a novel phenomenon can be found. First, setting , the system is in mixed firing modes when (Figure 10(a)), irregular square wave bursting when (Figure 10(b)), and regular square wave bursting when (Figure 10(c)). Then, setting , the system is in mixed firing modes when and (Figures 10(d) and 10(e)), and regular square wave bursting when (Figure 10(f)). When , we can find that the system displays mixed firing modes when , , and . From the above analysis, we can conclude that the region of frequency when the system displays the mixed firing modes decreases with increasing fractional-order.

(a)

(b)

(c)

(d)

(e)

(f)

3. Conclusion

In this study, the dynamic behaviors of the modified HR model were investigated, and diverse firing activities and some novel phenomena were found. First, according to the integer-order model and Cardan formula, the existence of Hopf bifurcation can be proven. The critical values of the direct external current stimulus for each fractional-order are also obtained. Then, when direct current stimulus is injected, the fractional-order model exhibits more complex dynamic behaviors compared with the integer-order model. In some direct current stimulus, the fractional-order can induce periodic firing for the integer-order model transitioning to chaotic firing, vice versa. In an appropriate direct current stimulus, from the bifurcation diagram with a bifurcation parameter of fractional-order, there is some transition of firing patterns such as chaotic bursting, periodic bursting (including different bursting modes), and periodic spiking. In addition, from the two-dimensional colored diagram, we can clearly find the chaotic firing range and periodic range. Finally, when periodic current stimulus is injected, at some appropriate fractional-order and parameters of periodic current stimulus, the fractional-order model can display regular/irregular bursting, regular/irregular spiking, multiple spiking, regular/irregular square wave bursting, and mixed firing modes. The region of frequency when the system displays the mixed firing modes decreases with the increasing fractional-order at some amplitude values of the external current stimulus.

Data Availability

The data used to support the findings of this study are included within the article, and other data used can be obtained from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

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: Publisher Site | Google Scholar
R. FitzHugh, “Impulses and physiological states in theoretical models of nerve membrane,” Biophysical Journal, vol. 1, no. 6, pp. 445–466, 1961.
View at: Publisher Site | Google Scholar
C. Morris and H. Lecar, “Voltage oscillations in the barnacle giant muscle fiber,” Biophysical Journal, vol. 35, no. 1, pp. 193–213, 1981.
View at: Publisher Site | Google Scholar
T. A. Krasimira, M. O. Hinke, R. Thorsten, and A. Sherman, “Full system bifurcation analysis of endocrine bursting models,” Journal of Theoretical Biology, vol. 264, no. 4, pp. 1133–1146, 2010.
View at: Google Scholar
C. S. Herrmann and A. Klaus, “Autapse turns neuron into oscillator,” International Journal of Bifurcation and Chaos, vol. 14, no. 2, pp. 623–633, 2004.
View at: Publisher Site | Google Scholar
A. Shilnikov and M. Kolomiets, “Methods of the qualitative theory for the Hindmarsh-Rose model: a case study-a tutorial,” International Journal of Bifurcation and Chaos, vol. 18, no. 8, pp. 2141–2168, 2008.
View at: Publisher Site | Google Scholar
D. G. Zhang, Q. Zhang, and X. Y. Zhu, “Exploring a type of central pattern generator based on Hindmarsh–Rose model: from theory to application,” International Journal of Neural Systems, vol. 25, no. 1, Article ID 1450028, 2015.
View at: Publisher Site | Google Scholar
H. Gu and B. Pan, “A four-dimensional 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
K. M. Wouapi, B. H. Fotsin, F. P. Louodop, K. F. Feudjio, Z. T. Njitacke, and T. H. Djeudjo, “Various firing activities and finite-time synchronization of an improved Hindmarsh-Rose neuron model under electric field effect,” Cognitive Neurodynamics, vol. 14, no. 3, pp. 375–397, 2020.
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
H. Yu, J. Wang, B. Deng, and X. Wei, “Firing patterns of map-based neuron under extracellular alternating-current feld,” Acta Biochimica et Biophysica Sinica, vol. 26, no. 10, pp. 907–918, 2010.
View at: Google Scholar
B. Li and Z. He, “Bifurcations and chaos in a two-dimensional discrete Hindmarsh-Rose model,” Nonlinear Dynamics, vol. 76, no. 1, pp. 697–715, 2014.
View at: Publisher Site | Google Scholar
H. Wang, Q. Wang, and Q. Lu, “Bursting oscillations, bifurcation and synchronization in neuronal systems,” Chaos, Solitons & Fractals, vol. 44, no. 8, pp. 667–675, 2011.
View at: Publisher Site | Google Scholar
B. Bao, Q. Yang, D. Zhu, Y. Zhang, Q. Xu, and M. Chen, “Initial-induced coexisting and synchronous firing activities in memristor synapse-coupled Morris-Lecar bi-neuron network,” Nonlinear Dynamics, vol. 99, no. 3, pp. 2339–2354, 2020.
View at: Publisher Site | Google Scholar
A. Mondal, R. K. Upadhyay, and J. Ma, “Bifurcation analysis and diverse firing activities of a modified excitable neuron model,” Cognitive Neurodynamics, vol. 13, no. 4, pp. 393–407, 2019.
View at: Publisher Site | Google Scholar
K. J. Wu, T. Q. Luo, H. W. Lu, and Y. Wang, “Bifurcation study of neuron firing activity of the modified Hindmarsh-Rose model,” Neural Computing & Applications, vol. 27, no. 3, pp. 739–747, 2016.
View at: Publisher Site | Google Scholar
K. J. Wu, W. Q. Li, and D. C. Wang, “Bifurcation of modified HR neural model under direct current,” Journal of Ambient Intelligence and Humanized Computing, vol. 11, no. 12, pp. 6211–6221, 2020.
View at: Publisher Site | Google Scholar
B. N. Lundstrom, M. H. Higgs, W. J. Spain, and A. L. Fairhall, “Fractional differentiation by neocortical pyramidal neurons,” Nature Neuroscience, vol. 1, no. 11, pp. 1335–1342, 2008.
View at: Publisher Site | Google Scholar
J. Ma, “Chaos theory and applications: the physical evidence, mechanism are important in chaotic systems,” Chaos Theory and Applications, vol. 4, no. 1, pp. 1–3, 2022.
View at: Google Scholar
J. A. T. Machado, V. Kiryakova, and F. Mainardi, “Recent history of fractional calculus,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, pp. 1140–1153, 2011.
View at: Publisher Site | Google Scholar
Z. Wang, X. Huang, and G. D. Shi, “Analysis of nonlinear dynamics and chaos in a fractional order financial system with time delay,” Computers & Mathematics With Applications, vol. 62, no. 3, pp. 1531–1539, 2011.
View at: Publisher Site | Google Scholar
Y. X. Fu, Y. M. Kang, and G. R. Chen, “Stochastic resonance based visual perception using spiking neural networks,” Frontiers in Computational Neuroscience, vol. 14, 2020.
View at: Google Scholar
T. J. Anastasio, “The fractional-order dynamics OF brain-stem vestibulooculomotor neurons,” Biological Cybernetics, vol. 72, no. 1, pp. 69–79, 1994.
View at: Publisher Site | Google Scholar
K. Rajagopal, N. Hasanzadeh, F. Parastesh, I. Hamarash, S. Jafari, and I. Hussain, “A fractional-order model for the novel coronavirus (COVID-19) outbreak,” Nonlinear Dynamics, vol. 101, no. 1, pp. 711–718, 2020.
View at: Publisher Site | Google Scholar
J. Dong, G. J. Zhang, Y. Xie, and Y. Hong, “Dynamic behavior analysis of fractional-order Hindmarsh-Rose neuronal model,” Cognitive Neurodynamics, vol. 8, no. 2, pp. 167–175, 2014.
View at: Google Scholar
S. B. He, K. H. Sun, X. Mei, B. Yan, and S. Xu, “Numerical analysis of a fractional-order chaotic system based on conformable fractional-order derivative,” European Physical Journal Plus, vol. 132, no. 1, p. 36, 2017.
View at: Publisher Site | Google Scholar
C. Letellier and L. A. Aguirre, “Dynamical analysis of fractional-order Rössler and modified Lorenz systems,” Physics Letters A, vol. 377, no. 28, pp. 1707–1719, 2013.
View at: Publisher Site | Google Scholar
L. M. Zhang, K. H. Sun, W. H. Liu, and S. B. He, “A novel color image encryption scheme using fractional-order hyperchaotic system and DNA sequence operations,” Chinese Physics B, vol. 26, no. 10, pp. 98–106, 2017.
View at: Publisher Site | Google Scholar
D. Cafagna and G. Grassi, “Hyperchaos in the fractional-order Rossler system with lowest-order,” International Journal of Bifurcation and Chaos, vol. 19, no. 1, pp. 339–347, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Xin Yang 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.

PDF Download Citation

Download other formats

Order printed copies

Views

211

Downloads

721

Citations