Hopf bifurcation analysis and control of three-dimensional Prescott neuron model

Neurons generate different firing patterns because of different bifurcations in the dynamical viewpoint. Various nerve diseases are relevant to the bifurcation of nervous system. Therefore, it is vital to control bifurcation since it may be potential ways of treating nerve diseases. This paper focuses on the critical Hopf bifurcation analysis and the problem of Hopf bifurcation control. We investigate the effects of key parameters on critical Hopf bifurcation and obtain the Hopf bifurcation occurrence region on parameter plane. With the theory of high-dimensional Hopf bifurcation, we analytically deduce the judgement criteria of Hopf bifurcation type for the three-dimensional models and judge the Hopf bifurcation type of Prescott model by using it. With application of the Washout filter, the subcritical Hopf bifurcation of Prescott model is controlled and converted to supercritical Hopf bifurcation. In addition, we make some discussions on Hopf bifurcation analysis of a coupled neural network. The results provided in this paper could bring new ways to controlling neurological diseases.


Introduction
The nervous system contains a variety of neurons displaying diverse firing patterns that coordinate and control the physiological processes in our body [1][2][3].Irregular firing patterns could lead to disorders in the nervous system, thus resulting in the occurrence of some nerve diseases, such as Parkinson's disease, Alzheimer's disease, epilepsy, etc. [4].In the view of dynamical systems, the appearance of different firing patterns is caused by bifurcations of neurons [5].For example, the excessive synchronization in basal ganglia causing the Parkinson's disease is considered to be the generation of a Hopf bifurcation of the firing potential [6].The seizurelike oscillations of epilepsy occur when the input to cortical neurons undergos a Hopf bifurcation of the averaged soma membrane potentials [7].In the case of schizophrenia, the symptoms appear resulting from the parameter-dependent bifurcation of neurotransmitter dopamine [8].Consequently, understanding the underlying mechanism of bifurcation and controlling bifurcation may throw new light on the diagnosis and treatment of neurological diseases.
Bifurcation refers to the phenomenon in dynamical systems when parameter values across the threshold (bifurcation value), the solution structure changes qualitatively.Bifurcation control means designing a controller to alter the bifurcation characteristics that can avoid undesired dynamical behaviors or obtain desired dynamical behaviors.Currently, various representative approaches for bifurcation control have been mentioned in literatures, such as linear state feedback [9], nonlinear state feedback [10], time-delayed feedback [11], harmonic balance approximation [12], quadratic invariants in normal forms [13], and hybrid control for discrete time model [14].However, the traditional ways of controlling bifurcation will unavoidably change the open-loop status of the system and result in wasting energy.While washout filter can avoid these disadvantages by preserving all the open-loop equilibrium points in the closed-loop system and will not fall into the degradation of system.
Model-based research is the principal method in bifurcation control.Y. Xie controlled the onset of Hopf bifurcation in Hodgkin-Huxley model [15].But taking many factors into account brings the model high dimensionality and makes the research more complexity [16].Doruk employed the FitzHugh-Nagumo model to control the firing patterns by nerve fiber cell membrane potential feedback [17].However, these simple models are hard to imitate abundant firing patterns generated in realistic neurons.Here, we adopt a three-dimensional Prescott model not only incorporating sufficient dynamical details to account for but also making the computation convenience.L. Ding did the stabilizing control of Hopf bifurcation in Hodgkin-Huxley model using a washout filter with linear control term [18], J. Wang controlled Hopf bifurcation in Hodgkin-Huxley model by a unified state-feedback method [19], but they did not get the analytical criteria to judge the type of Hopf bifurcation.M. A. Herrera-Valdez used bifurcation analysis to describe a connection between ion channel expression, bifurcation structure, and firing patterns [20], C. M. Liu studied the dynamical behaviors in Morris-Lecar model via bifurcation methods [21], however, they didn't investigate the critical Hopf bifurcation characteristics of the models.In this paper, we focus on the critical analysis and control of Hopf bifurcation in the three-dimensional Prescott model.
The paper proceeds as follows.In Section 2, we describe the three-dimensional Prescott model and analyze the critical Hopf bifurcation characteristics of the model.The impacts of two parameters on critical Hopf bifurcation are shown in this section.We deduce the analytical criteria to judge the type of Hopf bifurcation in three-dimensional model in Section 3. In Section 4, we control the type of Hopf bifurcation in Prescott model and change the subcritical Hopf bifurcation to supercritical Hopf bifurcation through a washout filter.Discussions on Hopf bifurcation analysis of a coupled network are presented in Section 5. Conclusions are given in Section 6.

Three-dimensional Prescott model and its Hopf bifurcation analysis
Two-dimensional Prescott neuron model is an improvement of Morris-Lecar neuron model, which consists of a fast activation variable and a slow recovery variable [1].Three-dimensional Prescott model introduces an adaptive ion current to make the model more realistic since adaptability is a prominent characteristic of neural dynamics.The dynamical equations are as follows: with: where represents the membrane potential, is the recovery variable of slow ionic currents.= ( )( − ) represents the fast ionic current, and = ( − ) represents the slow ionic current.= ( − ) is the adaptive current, and is the gating variable of the adaptive current.The conductances of the leak, sodium, potassium and adaptive currents are denoted by , , and .The corresponding reversal potentials are , , and , respectively.
For n-dimensional nonlinear system: where is the bifurcation parameter, ( ) is equilibrium point of the system, the eigenvalues corresponding to Jacobian matrix are one pair of complex characteristic roots , = ( ) ( ) and negative real eigenvalues ( ) ( = 3,4, ⋯ ).If = satisfies the conditions [22]: then = is the Hopf bifurcation value of Eq. (3).A set of limit cycles around the equilibrium point arise in the neighborhood of = .In the three-dimensional system, the equilibrium point becomes center -focus when Hopf bifurcation occurs, and its conjugate eigenvalues become purely imaginary.

Judging the Hopf bifurcation type of three-dimensional Prescott model
Consider a three-dimensional nonlinear system: where is the state vector, ∈ is bifurcation parameter.Let = ( , ) = 0 , find the equilibrium point ( ) = ( ), ( ), ( ) .Extract the linear portion, and the equation can be rewritten as: where is the Jacobian matrix of the equilibrium point, ( ) is nonlinear term: .The first two columns of are the real and imaginary parts of the eigenvector corresponding to eigenvalue = ( ) of .The third column of is the eigenvector corresponding to the real root .Therefore: where = is the Jordan matrix of the system, = ( + ( )) − is nonlinear term: The type of Hopf bifurcation can be determined according to the high dimensional Hopf bifurcation theory [22]: It can be simplified to: with: ) . ( As for three-dimensional system, = 3, so: with: where is the real eigenvalue of the equilibrium point, > 0 is the imaginary part of the purely imaginary characteristic roots.Therefore, the judgement formula of the Hopf bifurcation type of three-dimensional model is: When > 0, subcritical Hopf bifurcation occurs, the new equilibrium branch is unstable when the bifurcation parameter is greater than the bifurcation value.When < 0, supercritical Hopf bifurcation occurs, a stable equilibrium point becomes an unstable equilibrium point and produces a stable limit cycle when the bifurcation parameter changes to the bifurcation value.The new equilibrium branch is stable when the bifurcation parameter is greater than the bifurcation value. As for the three-dimensional Prescott model, = ( , , ), = , we can deduce the nonlinear term ( + ( )) in : The external current in Hopf bifurcation point is = 42.94μA/cm 2 , Hopf bifurcation point is ( , , ) = (−39.2767,0.0029,0.0804).The matrix in Hopf bifurcation point is: = 0.9999917 0.0 0.9998145 0.0002432193 0.0002817666 0.002762562 0.001670333 0.003698198 −0.01905972 .
Then we can calculate , and = ( + ( )) can be obtained.Use Eq. ( 20), = 0.001765 > 0 can be calculated.Thus subcritical Hopf bifurcation occurs in the model.As the bifurcation parameter increases, an unstable limit cycle shrinks to an equilibrium and makes it lose stability.The subcritical Hopf bifurcation of the three-dimensional Prescott model before control is shown in Fig. 6.It can be seen from the figure that the stable equilibrium loses its stability when the external current = 42.94μA/cm 2 .Subcritical Hopf bifurcation causes the new equilibrium branch of the system to lose stability.To avoid this happening, this paper design the controller based on Washout filter for subcritical Hopf bifurcation caused by the bifurcation parameter .Since adding the controller does not change the equilibrium position of the system, so as not to change the operating state of the system, thus to avoid waste of energy, we choose Washout filter to control the Hopf bifurcation type of neuron model.The dynamic equations after joining the controller are as follows: where is the nonlinear Washout controller, is the state variable of the controller, is the gain of the Washout controller, and is the inverse of the time constant, which is set to 0.5 in this paper.From Eq. ( 14) we know, as for four-dimensional system, = 4 , Re( 2) and Re( ) are needed to calculate Hopf bifurcation type judgement formula: When the external current is = 42.94μA/cm 2 , the Hopf bifurcation point is ( , , , ) = (−39.2767,0.0029,0.0804,−78.5533) .We can calculate and obtain the formula = −0.074957734− 0.00051748551 from Eqs. ( 20), (25-26).When < -0.0069, < 0. Therefore, the Hopf bifurcation of the model is supercritical Hopf bifurcation at = 42.94μA/cm 2 which is shown in Fig. 7.

Discussions
In this part, we will discuss how the presented method for controlling Hopf bifurcation of a single neuron could be used for a network of neurons.Here we consider a neural network of two coupled Prescott neurons, which can be described by the following equations: The critical step for Hopf bifurcation control of the coupled neurons is to obtain the characteristic polynomial coefficients and then the Hopf bifurcation points.The coupled Prescott model is a six-order transcendental equation.The difficulty of bifurcation analysis for high-dimensional nonlinear equations is the dimensionality reduction.In this paper, we reduce the six-order characteristic polynomial to two three-order characteristic polynomials, and therefore determine the Hopf bifurcation points of the model.Next step, we can apply the technique proposed previously to judge Hopf bifurcation type of the network and control the Hopf bifurcation if it is subcritical Hopf bifurcation.

Conclusions
In this paper, we study the Hopf bifurcation characteristics of the three-dimensional Prescott neuron model, analyze the effect of adaptation conductance and ion channel key parameter on the critical Hopf bifurcation point, and obtain the region on -plane which Hopf bifurcation can occur.Using high-dimensional Hopf bifurcation theory, we derive the analytical formula for judging Hopf bifurcation type of three-dimensional nonlinear model, and judge the Hopf bifurcation type of the Prescott model via the formula.By applying a Washout filter to control the three-dimensional model, the Hopf bifurcation type of Prescott model is converted from subcritical Hopf bifurcation to supercritical Hopf bifurcation.The method for judging and controlling Hopf bifurcation type proposed here also applies to a network of two coupled Prescott neurons.The results in this paper provide theoretical foundations for biophysical dynamics of neurons and lighten a new way of treating neurological diseases such as Parkinson and epilepsy.

Fig. 1 .Fig. 2 .
Fig. 1.The eigenvalue and external input curves of the three-dimensional Prescott model.The real part of , , and are indicated by blue solid, red dashed, and green dotted lines, respectively Fig. 2. The eigenvalue and external input curves in the critical Hopf bifurcation point.The real part of , , and are indicated by blue solid, red dashed, and green dotted lines, respectively Fig. 1 shows the eigenvalue curves and external current of the three-dimensional Prescott model.It can be seen from the figure that the real eigenvalue < 0 when the real part of the conjugate eigenvalues Re( ) = 0. Therefore, the two intersections of Re( ) curve and lateral axis are Hopf Bifurcation points of the model.The values of the first intersection are = -30.16mV, = 66.12 μA/cm 2 .The key parameter has a crucial impact on Hopf bifurcation of the model.Hopf bifurcation points change when has different values.When the Re( ) curve is tangent to lateral axis, which is shown in Fig. 2, there is only one Hopf bifurcation point.If Hopf bifurcation no longer happens in the model when is smaller than a certain value, we call the certain value critical Hopf bifurcation point in this paper.The values in the critical Hopf bifurcation point are = -17.3053mV, = -22.84mV.The nullcline, nullcline, phase trajectories of the three-dimensional Prescott model are shown in Fig. 3, we can see from the figure that the initial value starts from the Hopf bifurcation point,

Fig. 3 .
Fig. 3.The phase space in the critical Hopf bifurcation point.The nullcline, nullcline, and phase trajectories are indicated by blue solid, red dotted, and black dashed lines, respectively Adaptation plays a vital role in the information processing of neurons.Now we analyze the influence of the adaptive current and key parameter on Hopf bifurcation.When the conductance of adaptive current has different values, the changes of Hopf bifurcation points with are shown in Fig. 4(a).With increasing, the value of in the critical Hopf bifurcation point increases as well.On the other hand, when the values of are different, Hopf bifurcation points of the model varying with are depicted in Fig. 4(b).With growing, the value of in the critical Hopf bifurcation point also increases.Therefore, and together affect the generation of the critical Hopf bifurcation.In the -plane, the critical Hopf bifurcation curve are drawn in Fig.5.The left side of the dividing line does not exist Hopf bifurcation, that is, no matter how much the external current is, Hopf bifurcation will not occur.On the right side of the dividing line, Hopf bifurcation can occur in the three-dimensional Prescott model.

Fig. 5 .
Fig. 5.The critical Hopf bifurcation curve in the -plane Take linear transformation for the three-dimensional system: = + ( ) , where = Re( ) Im( ).The first two columns of are the real and imaginary parts of the eigenvector corresponding to eigenvalue = ( ) of .The third column of is the eigenvector corresponding to the real root .Therefore: 2182.HOPF BIFURCATION ANALYSIS AND CONTROL OF THREE-DIMENSIONAL PRESCOTT NEURON MODEL.CHUNHUA YUAN, JIANG WANG = Re + ,