It is shown that many neurological diseases are caused by the changes of firing patterns induced by bifurcations. Therefore, the bifurcation control may provide a potential therapeutic method of these neurodegenerative diseases. In this paper, we investigate the Hopf bifurcation control of the Morris-Lecar (ML) model with Homoclinic (Hc) bifurcation type by introducing a dynamic state-feedback control. The results indicate that the linear term can change the ML model from Hc bifurcation type to SNIC bifurcation type without changing the firing patterns. The cooperation of linear and cubic term can transform the ML model from the Hc bifurcation type to the Hopf bifurcation type, resulting in the transformation of firing patterns from type I to type II. Besides, we utilize the Poincare Birkhoff (PB) normal form method to derive the analytical expression of the bifurcation stability index for the controlled ML model with Hc bifurcation type, and the results show that the cubic term can regulate the criticality of the Hopf bifurcation. Numerical simulation results are consistent with the theoretical analysis.
Citation: Qinghua Zhu, Meng Li, Fang Han. Hopf bifurcation control of the ML neuron model with Hc bifurcation type[J]. Electronic Research Archive, 2022, 30(2): 615-632. doi: 10.3934/era.2022032
Related Papers:
[1]
Xueyong Zhou, Xiangyun Shi .
Stability analysis and backward bifurcation on an SEIQR epidemic model with nonlinear innate immunity. Electronic Research Archive, 2022, 30(9): 3481-3508.
doi: 10.3934/era.2022178
[2]
Fengrong Zhang, Ruining Chen .
Spatiotemporal patterns of a delayed diffusive prey-predator model with prey-taxis. Electronic Research Archive, 2024, 32(7): 4723-4740.
doi: 10.3934/era.2024215
[3]
Zhili Zhang, Aying Wan, Hongyan Lin .
Spatiotemporal patterns and multiple bifurcations of a reaction- diffusion model for hair follicle spacing. Electronic Research Archive, 2023, 31(4): 1922-1947.
doi: 10.3934/era.2023099
[4]
Li Li, Zhiguo Zhao .
Inhibitory autapse with time delay induces mixed-mode oscillations related to unstable dynamical behaviors near subcritical Hopf bifurcation. Electronic Research Archive, 2022, 30(5): 1898-1917.
doi: 10.3934/era.2022096
[5]
Mengting Sui, Yanfei Du .
Bifurcations, stability switches and chaos in a diffusive predator-prey model with fear response delay. Electronic Research Archive, 2023, 31(9): 5124-5150.
doi: 10.3934/era.2023262
[6]
Rina Su .
Dynamic analysis for a class of hydrological model with time delay under fire disturbance. Electronic Research Archive, 2022, 30(9): 3290-3319.
doi: 10.3934/era.2022167
[7]
Qixiang Wen, Shenquan Liu, Bo Lu .
Firing patterns and bifurcation analysis of neurons under electromagnetic induction. Electronic Research Archive, 2021, 29(5): 3205-3226.
doi: 10.3934/era.2021034
[8]
Ruizhi Yang, Dan Jin .
Dynamics in a predator-prey model with memory effect in predator and fear effect in prey. Electronic Research Archive, 2022, 30(4): 1322-1339.
doi: 10.3934/era.2022069
[9]
Weiyu Li, Hongyan Wang .
Dynamics of a three-molecule autocatalytic Schnakenberg model with cross-diffusion: Turing patterns of spatially homogeneous Hopf bifurcating periodic solutions. Electronic Research Archive, 2023, 31(7): 4139-4154.
doi: 10.3934/era.2023211
[10]
Miao Peng, Rui Lin, Zhengdi Zhang, Lei Huang .
The dynamics of a delayed predator-prey model with square root functional response and stage structure. Electronic Research Archive, 2024, 32(5): 3275-3298.
doi: 10.3934/era.2024150
Abstract
It is shown that many neurological diseases are caused by the changes of firing patterns induced by bifurcations. Therefore, the bifurcation control may provide a potential therapeutic method of these neurodegenerative diseases. In this paper, we investigate the Hopf bifurcation control of the Morris-Lecar (ML) model with Homoclinic (Hc) bifurcation type by introducing a dynamic state-feedback control. The results indicate that the linear term can change the ML model from Hc bifurcation type to SNIC bifurcation type without changing the firing patterns. The cooperation of linear and cubic term can transform the ML model from the Hc bifurcation type to the Hopf bifurcation type, resulting in the transformation of firing patterns from type I to type II. Besides, we utilize the Poincare Birkhoff (PB) normal form method to derive the analytical expression of the bifurcation stability index for the controlled ML model with Hc bifurcation type, and the results show that the cubic term can regulate the criticality of the Hopf bifurcation. Numerical simulation results are consistent with the theoretical analysis.
1.
Introduction
Individual neurons are the basic units of the biological nervous system, they process information by working together in neuronal circuits with specific synaptic connectivity. In response to an applied current, the signal can be propagated from one neuron to the next neuron by generating electrical impulses or spikes[1]. Hodgkin and Huxley (HH) firstly proposed the mathematical description of the neural model based on a large number of biological experiments[2], which built a new field of neuronal electrophysiology and laid the foundation for theoretical neuroscience. However, the HH model is a complex nonlinear system, and is difficult to be analyzed theoretically. Subsequently, researchers have successively proposed many simplified neuron models such as the Fitzhugh-Nagumo (FN) model [3], Morris-Lecar (ML) model[4], Hindmarsh-Rose (HR) model[5], and the modified forms of them[6,7,8], to explore the physiological mechanisms of neurons.
Previous studies have shown that many neurological dynamic diseases such as Epilepsy, Parkinson's disease, Alzheimer's disease and Schizophrenia may be induced by bifurcations caused by the changes of the regulatory parameters in the neural system, which leads to the disorders in the firing patterns of nervous systems[9,10,11,12]. In addition, the changes of firing patterns in neuron systems may cause neurons to produce a kind of plasticity, which finally affects learning, memory and other neural activities that rely on neuron-plasticity[13]. Therefore, a method is expected to be proposed to control the bifurcation mechanism of the neuron model, which can change the firing patterns from abnormal to normal under certain conditions.
Bifurcation control has attracted increasing attention due to its potential applications in many fields, such as engineering, biomedicine, and meteorology systems etc[14]. The main idea of bifurcation control is to change the inherent bifurcation characteristics of the original system by designing the controller without changing the parameters, thereby obtaining some desirable bifurcation behaviors. This includes delaying or advancing the onset of the bifurcation points, modifying the shape or type of bifurcations, and so on. Many bifurcation control methods have been provided in previous studies[15,16,17,18,19]. Xie et al. [10] studied the dynamic behaviors of Hopf bifurcation under the control of the washout filter in HH model. The results show that the linear term of the controller can control the position of the Hopf bifurcation point, and the cubic term can change the criticality of the Hopf bifurcation point. The washout filter dynamic feedback controller is introduced into the two-dimensional HR model to change the firing pattern from type I to type II by delaying the location of its inherent Hopf bifurcation[12]. Huang et al. [19] studied the Hopf bifurcation control of the ML model with type I under the control of washout filter-aided dynamic feedback controller. However, the firing pattern cannot be changed from type I to type II by using the washout filter-aided dynamic feedback controller in ML model. Therefore, a new controller should be introduced to address this issue. In addition, the analytical expression for determining the stability index of Hopf bifurcation criticality has not been obtained. Moreover, none of these studies involved the bifurcation control of the ML model with Hc bifurcation type.
Motivated by above analysis, in the present paper, we introduce a new controller proposed by Nguyen et al.[20] into the ML model to study the influence of the controller on the ML model with Hc bifurcation type. In addition, we give the analytical expression of stability index of the system to determine the criticality of Hopf bifurcation.
The structure of the paper is organized as follows. In Sect. 2, we reviewed the bifurcation dynamics and the conditions for the emergence of Hopf bifurcations of the ML model. In Sect. 3, the dynamic state-feedback control law was applied to the ML neuron model with Hc bifurcation type. Also, we derived the analytical expression of the bifurcation stability index for the controlled ML model with Hc bifurcation type. Conclusions are given in Sect. 4.
2.
Reviews of the ML model
2.1. ML model
The ML neuron model is a simplified form of HH model and its mathematical description can be expressed by the following nonlinear differential equations:
where V is the membrane potential, and N is the recovery variable. I is the current stimulation and is treated as the bifurcation parameter. C is the membrane capacitance. Vi(i=K,Ca,L) represent the reversal potential, and gi(i=K,Ca,L) represent the maximum conductance. Vi(i=1,2,3,4) are constant potentials, and ϕ represents the temperature factor. Here, we named system (1) as the original system.
2.2. The bifurcation dynamics of the ML model
There are three types of bifurcations in ML model by choosing different parameter sets, including Hopf bifurcation, saddle node bifurcation on invariant cycle (SNIC) and homoclinic (Hc) bifurcation (see Table 1). In this paper, we mainly focus on the Hopf bifurcation control of the ML model with Hc bifurcation type. The bifurcation diagram of the ML model with Hc bifurcation type is shown in Figure 1. From Figure 1 we can see that as the bifurcation parameter I increases, the ML model with Hc bifurcation type is divided into six types of dynamic regions. In region I (when I<−9.949), the system has only one stable equilibrium point in the lower branch; In region II (when −9.949<I<34.94), the system contains three equilibrium points, the stable equilibrium point in the upper branch, the unstable equilibrium point in the middle branch and the stable equilibrium point in the lower branch; In region III (when 34.94<I<36.32), the system contains a stable limit cycle and three equilibrium points; In region IV (when 36.32<I<39.96), the system contains a stable limit cycle, an unstable limit cycle and three equilibrium points; In region V(when 39.96<I<40.56), the system contains a stable limit cycle, an unstable limit cycle and a stable equilibrium point; In region VI(when I>40.56), the system contains a stable equilibrium point. The corresponding phase diagrams of six typical dynamical behaviors are shown in Figure 2. Compared it with [19], the ML model with Hc bifurcation type has more dynamical behaviors than that with SNIC bifurcation type. In particular, Hc bifurcation occurs in the system when I=34.94, and the corresponding phase diagram is shown in Figure 3.
Table 1.
Parameter setting in the simulation and experiment.
Figure 1.
(a) Bifurcation diagram of the ML model with Hc bifurcation type. (b) The enlarged representation of (a). Here, the red and black lines represent stable equilibrium point and unstable equilibrium point, respectively. The green and blue circles represent stable limit cycle and unstable limit cycle, respectively. SNi(i=1,2) represent saddle node bifurcation points, H represents Hopf bifurcation point, and Hc represents homoclinic bifurcation point. I−VI represent six types of dynamical behavior regions, respectively. SN1: I=−9.949μA/cm2, Hc: I=34.94μA/cm2, H: I=36.32μA/cm2, SN2: I=39.96μA/cm2.
Figure 2.
Six types of phase diagrams of the ML model with Hc bifurcation type. Here, the blue and red curves represent N-nullcline and V-nullcline, respectively. The black solid line and the black dashed line represent stable limit cycle (SLC) and unstable limit cycle (ULC), respectively. The red circle and the black circle represent the stable equilibrium point (SEP) and the unstable equilibrium point (UEP), respectively. (a) I=−30μA/cm2, (b) I=30μA/cm2, (c) I=35.5μA/cm2, (d) I=37.5μA/cm2, (e) I=40μA/cm2, (f) I=45μA/cm2.
Some researchers have proposed that the firing patterns can be determined by the types of bifurcation[21,22]. Specifically, type I firing pattern occurs when the system undergoes SNIC bifurcation or Hc bifurcation, while type II firing pattern occurs when the system undergoes Hopf bifurcation. For type I, the system can fire at any low frequency; for type II, the system fires from a positive minimum frequency. Figure 4 shows the firing patterns of the ML model under the three types of bifurcations. From Figure 4 we can see that the ML model with Hopf bifurcation type shows the type II firing pattern (see Figure 4b), and the corresponding time evolutions are shown in Figure 4a; the ML model with SNIC bifurcation or Hc bifurcation type both show the type I firing pattern (see Figure 4d and 4f), and the corresponding time evolutions are shown in Figure 4c and 4e, respectively.
Figure 4.
(a) (c) (e) Time evolution diagrams of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively. (b) (d) (f) Frequency-I curves of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively.
2.4. Conditions for the emergence of Hopf bifurcation
Consider the following n-dimensional nonlinear system:
dxdt=f(x,μ),
(3)
where x∈Rn represents the state variables, μ∈R represents the bifurcation parameter. Assume that the system has a fixed point x0 when the bifurcation parameter μ=μ0 and the corresponding eigenvalue λμ0=αμ0+iβμ0. The Hopf bifurcation occurs at μ=μ0 when the system satisfies the following conditions: βμ0≠0, αμ0=0, and ddμ(αμ)∣μ=μ0>0. However, the eigenvalues of the system cannot be obtained easily when the system is a high-dimensional system. Therefore, the following Routh-Hurwitz stability criterion can be used to determine whether the Hopf bifurcation occurs[23]. Assume that the corresponding characteristic polynomial of system (3) is as follows:
Based on the Routh-Hurwitz stability criterion, the Hopf bifurcation occurs at μ=μ0 when the system satisfies the following conditions:
(I) Eigenvalue crossing condition:
pn(μ0)>0,Di(μ0)>0,i=1,⋯,n−2,Dn−1(μ0)=0.
(7)
(II) Transversality condition:
dDn−1(μ)dμ∣μ=μ0≠0.
(8)
3.
Hopf bifurcation control of the ML model
3.1. The dynamic state-feedback controller
In this section, we provide a dynamic state-feedback controller to control the Hopf bifurcation of the original system. Consider the following nonlinear system:
u=u(x,y),˙y=g(x,y).
(9)
where y∈Rm(1≤m≤n) represents the state variables of the controller, u(x,y) represents the feedback control and g(x,y) is a smooth function. In this paper, we utilize the following feedback control law proposed by Nguyen and Hong[20].
ui(xi,yi)=k1ixi+k3i(xi−x0i)3−liyi,˙yi=ui(xi,yi).
(10)
where x0i(i=1,2,...,m) are the equilibrium points of the system at the Hopf bifurcation point, k1i and k3i represent the linear term control gain and the cubic term control gain, respectively. li are constant parameters. Then we obtain the following controlled system:
According to Eq (10), we have ui=0 when ˙y=ui=0. So, if x0 is the equilibrium point of the original system, (x0,y0) is the equilibrium point of the controlled system, where y0=(y01,y02,⋯,y0m) and then y0i=(k1ix0i+k3i(x0i−x01i)3)/li(i=1,2,⋯,m). That is, the equilibrium points of the original system are the equilibrium points of the controlled model. The reason for using state feedback control to conduct the Hopf bifurcation control is that the equilibrium structure of the original system and the controlled system remains unchanged during the control process. In the present paper, we add the controller to the membrane potential V, and then we obtain the following controlled system:
where V0 is the equilibrium membrane potential of the original system at the Hopf bifurcation point.
3.2. Hopf bifurcation control of the ML model with Hc bifurcation type
The Hopf bifurcation control of neuron models has been widely studied, including HH model, HR model, ML model with Hopf bifurcation type or SNIC bifurcation type[10,11,12,19,20]. However, the Hopf bifurcation control of the ML model with Hc bifurcation type has not been studied. In addition, the controllers used in previous researches only studied the system with single Hopf bifurcation point. Therefore, in the present paper, we utilize the dynamic state-feedback control proposed by Nguyen et al. to make the Hopf bifurcation point change from one to two. The parameter values of the model are ϕ=0.23,gCa=4,V3=12,V4=17.4,VCa=120,VK=−84,VL=−60,gK=8,gL=2,V1=−1.2,V2=18,C=20,l=0.1, and I is treated as the bifurcation parameter.
3.2.1. The transition from Hc bifurcation type to SNIC bifurcation type
The Hopf bifurcation point is located at I=36.32 in the original ML model with Hc bifurcation type (Figure 1). To change the model from the Hc bifurcation type to the SNIC bifurcation type, the Hopf bifurcation point needs to be moved forward. Here, we change the Hopf bifurcation point from I=36.32 to I0=70. The equilibrium point of the original system is (V0=6.7697,N0=0.35407). Therefore, the equilibrium point of the controlled system is (V0=6.7697,N0=0.35407,y0=67.697k1), and then we obtain the Jacobian matrix of the controlled system at the equilibrium point:
Solving Eq (16), we have k1=0.0893. Then, we have:
dD2(I)dI∣I=I0=−0.2032×10−4≠0.
(17)
which satisfy the transversal condition (II).
The bifurcation diagram of the controlled system when k1=0.0893,k3=0 is shown in Figure 5a. As expected, the Hopf bifurcation point of the ML model with Hc bifurcation type moves forward from I=36.32 to I=70. Figure 5b shows the two-parameter bifurcation diagram of bifurcation parameter I with respect to k1, which further verifies that Hopf bifurcation occurs at I=70 when k1=0.0893. Now, the model shows a change from the Hc bifurcation type to the SNIC bifurcation type (see Figure 5a), but does not change the firing pattern of the system (see Figure 5c and 5d).
Figure 5.
(a) Bifurcation diagram of the controlled system when k1=0.0893,k3=0,I0=70. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled model with different bifurcation parameter I when k1=0.0893,k3=0. (d) Frequency-I curves of the controlled system when k1=0.0893,k3=0.
The results of this section show that the ML model can be transformed from the Hc bifurcation type to the SNIC bifurcation type by introducing a dynamic state-feedback control. In addition, there is only one Hopf bifurcation point in the system when k1≠0,k3=0 (Figure 5a), that is, k1 can change the location of the Hopf bifurcation point of the ML neuron model with Hc bifurcation type, resulting in the change of bifurcation type from Hc type to SNIC type.
3.2.2. The transition from Hc bifurcation type to Hopf bifurcation type
Previous study has investigated the transition of ML model from SNIC bifurcation type to Hopf bifurcation type[18]. Here, we will study the transition of the ML model from Hc bifurcation type to Hopf bifurcation. First, the inherent Hopf bifurcation point in the Hc bifurcation type needs to be moved forward, and then a new Hopf bifurcation point is constructed in front of the moved inherent Hopf bifurcation point. Here, we expect to change the inherent Hopf bifurcation point from I=36.32 to I1=200, and then a new Hopf bifurcation point is constructed at I2=60.
When I1=200, the equilibrium point of the controlled system is (V0=12.94,N0=0.52697,y0=129.4k1), and then we obtain the Jacobian matrix of the controlled system at the equilibrium point:
Solving Eq (20), we have k1=0.34966383. Then, substituting k1=0.34966383 into the controlled system when I2=60, we obtain the equilibrium point of the controlled system (V0=6.1368,N0=0.33762,y0=21.46+3148.761k3). The Jacobian matrix of the controlled system at the equilibrium point is as follows:
Solving Eq (23), we have k3=−0.00205656. Figure 6 shows the dynamic behaviors of the controlled system when k1=0.34966383,k3=−0.00205656. From Figure 6a we can see that the Hopf bifurcation point of the ML neuron model with Hc bifurcation type is moved forward from I=36.32 to I1=200, and a new Hopf bifurcation point is constructed at I2=60. Figure 6b shows the two-parameter bifurcation diagram of bifurcation parameter I with respect to k1, which further verifies that the model undergoes the Hopf bifurcation at I1=200,I2=60 when k1=0.0893. Figure 6c shows the time evolutions of the controlled system with different bifurcation parameter I. From Figure 6c we can see that the system fires from a positive minimum frequency, which is further verified in Figure 6d. Now, the model shows a change from Hc bifurcation type to Hopf bifurcation type (see Figure 6a), and the firing pattern changes from type I to type II (see Figure 6d).
Figure 6.
(a) Bifurcation diagram of the controlled system with k1=0.34966383 and k3=−0.00205656. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled system with different bifurcation parameter I when k1=0.34966383 and k3=−0.00205656. (d) Frequency-I curves of the controlled system when k1=0.34966383 and k3=−0.00205656.
The results in this section show that the ML model can be transformed from the Hc bifurcation type to the Hopf bifurcation type by introducing dynamic state-feedback control. In addition, the model performs type II firing pattern when k1≠0,k3≠0, that is, the combination of k1 and k3 can change the firing patterns of the ML neuron model with Hc bifurcation type.
3.2.3. Theoretical analysis of the criticality of Hopf bifurcation
As mentioned above, the criticality of the Hopf bifurcation can be changed by a dynamic state-feedback controller. In this section, we derive the stability coefficient β2 of the Hopf bifurcation by utilizing the Poincare Birkhoff (PB) normal form method. The Hopf bifurcation is subcritical when β2>0, resulting in an unstable limit cycle; the Hopf bifurcation is supercritical when β2<0, resulting in a stable limit cycle.
Consider the following controlled system:
˙x=f(x,μ),
(24)
where x=(V,N,y)T, f=(f1,f2,f3)T. The equilibrium point of the system at μ0: x0(μ0)=(V0(μ0),N0(μ0),y0(μ0))T. Rewrite the controlled system as: ˙x=f(x,μ)=Jx+g(x), where J is the Jocabian matrix of the system at the equilibrium point, g(x) is the nonlinear term:
Then, we conduct the linear transformation of the controlled system: x=Py+x0(μ0), where P=[Re(v1),Im(v1),v3], Re(v1) and Im(v1) are the real and imaginary parts of the eigenvector v1 corresponding to eigenvalues λ1, and v3 is the eigenvector corresponding to the eigenvalue λ3. Then, we have:
where λ3 is the real characteristic root of the Jacobian matrix, ω0 is the imaginary part of the pure virtual characteristic root of the Jacobian matrix, and ω0>0. Solving Eqs (39) and (40) we obtain:
Because Ay is the first-order function with respect to y and the second or more derivative of y is required in the following calculations, the derivative of y is omitted. We only calculate the following equation: F=P−1f(Py+x0(μ0)).
Next, we calculate the stability coefficient β2 of the controlled system when k1=0.3497, k3=−0.0021, μ0=I1=200. The equilibrium point of the controlled system is (V0=12.94,N0=0.52697,y0=45.2465), and the corresponding Jacobian matrix is:
J=(0.2122−38.7760−0.10.0066−0.230100.34970−0.1).
(46)
Then, we obtain the eigenvalues and eigenvectors: λ1=0.4938i, λ2=−0.4938i, λ3=−0.1179.
Finally, substituting the parameters into the stability coefficient formula we have β2=−9.9169×10−4<0, that is, the controlled system undergoes the supercritical Hopf bifurcation at (V0,N0)=(12.94,0.52697) when I1=200. The membrane potential curve and phase trajectory of the controlled system at the bifurcation point are shown in Figure 7a and 7b. From Figure 7b we can see that the phase trajectory of the controlled system converges to a stable limit cycle when I1=200. However, the phase trajectory of the original system converges to a stable equilibrium point when I1=200 (see Figure 7c and 7d). The similar results can be obtained when I2=60, so we omit them here.
Figure 7.
Time evolutions and phase diagrams of the controlled system (a-b) and original system (c-d) when I1 = 200μA/cm2.
In this paper, we investigate the Hopf bifurcation control of the ML model with Hc bifurcation type by utilizing the method of dynamic state-feedback control. The results show that (1) the location of the Hopf bifurcation point can be changed by introducing the state-feedback control without changing the model parameters; (2) the control gain k1 transforms the ML model from the Hc bifurcation type to the SNIC bifurcation type, but does not change the firing pattern of the model; (3) the cooperation of control gain k1 and k3 transforms the ML model from the Hc bifurcation type to the Hopf bifurcation type, and the firing pattern shows a change from Type I to Type II. In addition, we find that the criticality of the Hopf bifurcation of the system can be changed by introducing the control gain k3, that is, the criticality of the Hopf bifurcation is changed from the subcritical Hopf bifurcation to the supercritical Hopf bifurcation. Many neurodegenerative diseases are caused by the change of firing patterns. Therefore, the results of this paper may have potential applications for the treatment of these diseases.
Compared to previous studies, the conclusions of this paper further enrich the research results. Some scholars have studied the Hopf bifurcation control based on the ML model with SNIC bifurcation type by introducing a washout filter-aided dynamic feedback controller[19]. In their research, they only studied how the location and criticality of the Hopf bifurcation changed under different control conditions. However, they cannot perform the change of the Hopf bifurcation point from one to two. Therefore, it is impossible to study how to control the firing pattern of the system changes from type I to type II. Here, we utilize a new dynamic state feedback controller to control the Hopf bifurcation of the ML model with Hc bifurcation type. This controller not only reproduces the previous results, but also realizes the control of two Hopf bifurcation points. We further study the effects of the controller on firing patterns, providing potential application value for the treatment of some dynamic diseases. In addition, the stability index to determine the criticality of the Hopf bifurcation is given in our paper. The mutual verification of theoretical analysis and numerical simulation improves the integrity of this research.
Acknowledgments
This work was supported by the national Natural Science Foundation of China under Grant 11972115 and the Fundamental Research Funds for the Central Universities. The authors would like to thank the reviewers for their comments.
Conflict of interest
The authors declare there is no conflicts of interest.
References
[1]
G. B. Ermentrout, D. H. Terman, Mathematical foundations of neuroscience, Springer-Verlag, New York, 2010.
[2]
A. L. Hodgkin, A. F. Huxley, A quantitative description of membrane current and its application to conduction and excitation in nerve, J. Physiol., 117 (1952), 500–544. https://doi.org/10.1113/jphysiol.1952.sp004764 doi: 10.1113/jphysiol.1952.sp004764
[3]
R. FitzHugh, Impulses and physiological states in theoretical models of nerve membrane, Biophys. J., 1 (1961), 445466. https://doi.org/10.1016/s0006-3495(61)86902-6 doi: 10.1016/s0006-3495(61)86902-6
[4]
C. Morris, H. Lecar, Voltage oscillations in the barnacle giant muscle fiber, Biophys. J., 35 (1981), 193–213. https://doi.org/:10.1016/s0006-3495(81)84782-0 doi: 10.1016/s0006-3495(81)84782-0
[5]
J. L. Hindmarsh, R. M. Rose, A model of neuronal bursting using three coupled first order differential equations, Proc. R. Soc. Lond. Ser. B, 221 (1984), 87–102. https://doi.org/10.2307/35900 doi: 10.2307/35900
[6]
Q. Zheng, J. W. Shen, Turing instability induced by random network in FitzHugh-Nagumo model, Appl. Math. Comput., 381 (2020), 125304. https://doi.org/10.1016/j.amc.2020.125304 doi: 10.1016/j.amc.2020.125304
[7]
A. Mondal, R. K. Upadhyay, J. Ma, B. K. Yadav, S. K. Sharma, A. Mondal, Bifurcation analysis and diverse firing activities of a modified excitable neuron model, Cogn. Neurodyn., 13 (2019), 393–407. https://doi.org/10.1007/s11571-019-09526-z doi: 10.1007/s11571-019-09526-z
[8]
M. K. Wouapi, B. H. Fotsin, E. B. M. Ngouonkadi, F. F. Kemwoue, Z. T. Njitacke, Complex bifurcation analysis and synchronization optimal control for Hindmarsh-Rose neuron model under magnetic flow effect, Cogn. Neurodyn., 15 (2021), 315–347. https://doi.org/10.1007/s11571-020-09606-5 doi: 10.1007/s11571-020-09606-5
[9]
J. Milton, P. Jung, Brain defibrillators: synopsis, problems and future directions. In: Epilepsy as a Dynamic Disease, Springer, Berlin Heidelberg, 2003.
[10]
Y. Xie, L. Chen, Y. Kang, K. Aihara, Controlling the onset of Hopf bifurcation in the Hodgkin-Huxley model, Phys. Rev. E., 77 (2008), 061921. https://doi.org/10.1103/PhysRevE.77.061921 doi: 10.1103/PhysRevE.77.061921
[11]
L. Nguyen and K. Hong, Hopf bifurcation control via a dynamic state-feedback control, Phys. Lett. A, 376 (2012), 442–446. https://doi.org/10.1016/j.physleta.2011.11.057 doi: 10.1016/j.physleta.2011.11.057
[12]
Y. Xie, Change in types of neuronal excitebility via bifurcation control, Chin. Q. Mech., 31 (2010), 58–63. https://doi.org/10.1103/PhysRevE.77.021917 doi: 10.1103/PhysRevE.77.021917
[13]
D. J. Schulz, Plasticity and stability in neuronal output via changes in intrinsic excitability: it's what's inside that counts, J. Exp. Biol., 209 (2006), 4821–4827. https://doi.org/10.1242/jeb.02567 doi: 10.1242/jeb.02567
[14]
G. R. Chen, J. L. Moiola, H. O. Wang, Bifurcation Control: Theories, Methods, and Applications, Int. J. Bifurc. Chaos, 10 (2000), 511–548. https://doi.org/10.1142/S0218127400000360 doi: 10.1142/S0218127400000360
[15]
X. F. Liao, S. W. Li, K. W. Wong, Hopf bifurcation on a twoneuron system with distributed delays: a frequency domain approach, Nonlinear Dyn., 31 (2003), 299–326. https://doi.org/10.1023/A:1022928118143 doi: 10.1023/A:1022928118143
[16]
P. Yu, G. R. Chen, Hopf bifurcation control using nonlinear feedback with polynomial functions, Int. J. Bifurc. Chaos, 14 (2004), 1683–1704. https://doi.org/10.1142/S0218127404010291 doi: 10.1142/S0218127404010291
[17]
J. Jiang, Y. L. Song, Delay-induced Bogdanov-Takens bifurcation in a Leslie-Gower predator-prey model with nonmonotonic functional response, Commun. Nonlinear. Sci. Numer. Simul., 19 (2014), 2454–2465. https://doi.org/10.1016/j.cnsns.2013.11.020 doi: 10.1016/j.cnsns.2013.11.020
[18]
M. Xiao, D. Ho, J. Cao, Time-delayed feedback control of dynamical small-world networks at Hopf bifurcation, Nonlinear Dyn., 58 (2009), 319–344. https://doi.org/10.1007/s11071-009-9485-0 doi: 10.1007/s11071-009-9485-0
[19]
C. L. Huang, W. Sun, Z. G. Zheng, J. H. Lu, S. H. Chen, Hopf bifurcation control of the M-L neuron model with type I, Nonlinear Dyn., 87 (2017), 755–766. https://doi.org/10.1007/s11071-016-3073-x doi: 10.1007/s11071-016-3073-x
[20]
L. H. Nguyen, K. S. Hong, Hopf bifurcation control via a dynamic state-feedback control, Phys. Lett. A, 376 (2012), 442–446. https://doi.org/10.1016/j.physleta.2011.11.057 doi: 10.1016/j.physleta.2011.11.057
[21]
J. Rinzel, G. B. Ermentrout, Analysis of Neural Excitability and Oscillations, In: Koch, C. and Segev, I., Eds., Methods in Neuronal Modeling: From Synapses to Networks, MIT Press, Cambridge, 1989,135–169.
[22]
E. M. Izhikevich, Dynamical Systems in Neuroscience, Cambridge, Massachusetts: MIT Press, 2007.
[23]
W. Liu, Criterion of Hopf bifurcations without using eigenvalues, J. Math. Anal. Appl., 182 (1994), 250–256. https://doi.org/10.1006/jmaa.1994.1079 doi: 10.1006/jmaa.1994.1079
This article has been cited by:
1.
Saureesh Das,
Recurrence quantification and bifurcation analysis of electrical activity in resistive/memristive synapse coupled Fitzhugh–Nagumo type neurons,
2022,
165,
09600779,
112772,
10.1016/j.chaos.2022.112772
2.
Jiajin He, Min Xiao, Jing Zhao, Zhengxin Wang, Yi Yao, Jinde Cao,
Tree-structured neural networks: Spatiotemporal dynamics and optimal control,
2023,
164,
08936080,
395,
10.1016/j.neunet.2023.04.039
Zhihui Wang, Xindan Wei, Lixia Duan,
The Critical Role of Self-Inhibition in Epileptic Seizures Based on Dynamic Analysis of the Corticothalamic Model,
2025,
35,
0218-1274,
10.1142/S0218127425500130
Qinghua Zhu, Meng Li, Fang Han. Hopf bifurcation control of the ML neuron model with Hc bifurcation type[J]. Electronic Research Archive, 2022, 30(2): 615-632. doi: 10.3934/era.2022032
Qinghua Zhu, Meng Li, Fang Han. Hopf bifurcation control of the ML neuron model with Hc bifurcation type[J]. Electronic Research Archive, 2022, 30(2): 615-632. doi: 10.3934/era.2022032
Figure 1. (a) Bifurcation diagram of the ML model with Hc bifurcation type. (b) The enlarged representation of (a). Here, the red and black lines represent stable equilibrium point and unstable equilibrium point, respectively. The green and blue circles represent stable limit cycle and unstable limit cycle, respectively. SNi(i=1,2) represent saddle node bifurcation points, H represents Hopf bifurcation point, and Hc represents homoclinic bifurcation point. I−VI represent six types of dynamical behavior regions, respectively. SN1: I=−9.949μA/cm2, Hc: I=34.94μA/cm2, H: I=36.32μA/cm2, SN2: I=39.96μA/cm2
Figure 2. Six types of phase diagrams of the ML model with Hc bifurcation type. Here, the blue and red curves represent N-nullcline and V-nullcline, respectively. The black solid line and the black dashed line represent stable limit cycle (SLC) and unstable limit cycle (ULC), respectively. The red circle and the black circle represent the stable equilibrium point (SEP) and the unstable equilibrium point (UEP), respectively. (a) I=−30μA/cm2, (b) I=30μA/cm2, (c) I=35.5μA/cm2, (d) I=37.5μA/cm2, (e) I=40μA/cm2, (f) I=45μA/cm2
Figure 3. Homoclinic orbit of the ML model with Hc bifurcation type. I=34.94μA/cm2
Figure 4. (a) (c) (e) Time evolution diagrams of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively. (b) (d) (f) Frequency-I curves of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively
Figure 5. (a) Bifurcation diagram of the controlled system when k1=0.0893,k3=0,I0=70. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled model with different bifurcation parameter I when k1=0.0893,k3=0. (d) Frequency-I curves of the controlled system when k1=0.0893,k3=0
Figure 6. (a) Bifurcation diagram of the controlled system with k1=0.34966383 and k3=−0.00205656. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled system with different bifurcation parameter I when k1=0.34966383 and k3=−0.00205656. (d) Frequency-I curves of the controlled system when k1=0.34966383 and k3=−0.00205656
Figure 7. Time evolutions and phase diagrams of the controlled system (a-b) and original system (c-d) when I1 = 200μA/cm2