^{1}

^{*}

^{1}

In this paper,we study the complicated dynamics of general Morris-Lecar model with the impact of Cl
^{-} fluctuations on firing patterns of this neuron model. After adding Cl
^{-} channel in the original Morris-Lecar model, the dynamics of the original model such as its bifurcations of equilibrium points would be changed and they occurred at different values compared to the primary model. We discover these qualitative changes in the point of dynamical systems and neuroscience. We will conduct the co-dimension two bifurcations analysis with respect to different control parameters to explore the complicated behaviors for this new neuron model.

Computational neuroscience uses knowledge of biology and combines with mathematical modeling to simulate the fluctuations of neuron cells and biological physiological characteristics of them [

The Morris-Lecar model is a two-dimensional biological neuron model which produces action potentials and neuronal fluctuations; in dynamical system point of view, oscillatory behaviors related to its ionic channels in giant barnacle muscle fibers [

{ C M d V d t = I a p p − g L ( V − E L ) − g K n ( V − E K ) − g C a m ∞ ( V ) ( V − E C a ) = I a p p − I i o n ( V , n ) , d n d t = ϕ ( n ∞ ( V ) − n ) / τ n ( V ) ,

where

m ∞ ( V ) = 1 2 [ 1 + tanh ( V − V 1 ) / V 2 ] ,

τ n ( V ) = 1 / cosh ( ( V − V 3 ) / ( 2 V 4 ) ) ,

n ∞ ( V ) = 1 2 [ 1 + tanh ( ( V − V 3 ) / V 4 ) ] .

and

I i o n ( V , n ) = g L ( V − E L ) + g K n ( V − E K ) + g C a m ∞ ( V ) ( V − E C a )

where V demonstrates membrane potential, and n the activation variable of the persistent K^{+} current, so it is a two-dimensional vector ( V , n ) . E K , E C a , and E L denote the Nernst equilibrium potentials. I a p p demonstrates the injected current and I i o n the ionic current. Parameter ϕ is a temperature factor. g L is leak membrane conductance, g K is potassium membrane conductance and g C a is calcium membrane conductance. Moreover, C M is the total membrane capacitance. Also, the voltage-sensitive steady-state activation function m ∞ ( V ) and n ∞ ( V ) , and the time constant τ n ( V ) can be measured experimentally.

In this study, we will discover the dynamics of Morris-Lecar model with the impact of Cl^{−} channel that has the following form:

{ C M d V d t = I a p p − g L ( V − E L ) − g K n ( V − E K ) − g N a m ∞ ( V ) ( V − E N a ) − g e x c ( V − E e x c ) − g i n h ( V − E C l ) , d n d t = Φ ( n ∞ ( V ) − n ) / τ n ( V ) , (1)

where

m ∞ ( V ) = 1 2 [ 1 + tanh ( ( V − V 1 ) / V 2 ) ] ,

τ n ( V ) = 1 / cosh ( ( V − V 3 ) / ( 2 V 4 ) ) ,

n ∞ ( V ) = 1 2 [ 1 + tanh ( ( V − V 3 ) / V 4 ) ] .

Here, E C l , E e x c , E L , E K , and E N a the Nernst equilibrium potentials. I a p p the injected current, g i n h , g e x c , g L , g K , g N a , inhibitory, excitatory, leak, potassium and sodium membrane conductance, C M the total membrane capacitance. In Neuroscience point of view, we can describe the electrical properties of a neuron with the help of an Equivalent Circuit. Here, we can write the total current I as the following form

I = C V ˙ + I N a + I K + I C l Kirchhoff’s Law

And Then

C V ˙ = I − I N a − I K − I C l

And also we know

E K < E C l < V < E N a

where I N a (inward currents) is negative and also I K , I C l are positive. Moreover, The inward currents depolarize the neuron and outward currents hyperpolarize it. Here, we have designed an Equivalent Circuit with the help of Simulink for model (1) (see

In this paper, we focus on this novel Morris-Lecar model for a certain range of parameters value which demonstrates different types of local bifurcations such as Hopf bifurcation and homoclinic bifurcation. We use rigorous analytic dynamical systems tools to study the complicated neuronal behaviors of this model such as normal form theory. We conduct bifurcation analysis with respect to different parameters to display the effect of various biological parameters on spiking dynamics of the system. Moreover, we are interested to discover the co-dimension two bifurcations such as Bautin or generalized Hopf and Bogdanov-Takens and cusp bifurcations and we present the normal form of these bifurcations.

In dynamical system point of view, for the Hopf bifurcation, a stable focus loses

its stability and becomes an unstable focus as we change the parameter for bifurcation, and the stable focus which is an attractor becomes a limit cycle that this limit cycle in phase space is a closed curve [

Using the parameters values corresponding to the Hopf case in

We have demonstrated the time series for Hopf case in

And for other values of I a p p near the Hopf bifurcation, we can see the changing in the behaviors of solutions. The neurons that show the subcritical Hopf bifurcation are resonators that exhibit bistability and the neurons that show supercritical Hopf bifurcations are resonators that exhibit monostability [

Continuing further, we compare the trajectories of two types of Hopf bifurcation. For model (1), we have both types of Hopf bifurcation. For supercritical Hope bifurcation that occurs for I a p p ≈ 208 a small limit cycle is born at the bifurcation point. This limit cycle will grow as we increase further the injected current (the amplitude of limit cycle increases with further increase of parameter). In

Parameter | Hopf | SNLC | Homoclinic |
---|---|---|---|

∅ | 0.04 | 0.067 | 0.23 |

g_{Ca} | 4.4 | 4 | 4 |

V_{3} | 2 | 12 | 12 |

V_{4} | 30 | 17.4 | 17.4 |

E_{Ca} | 120 | 120 | 120 |

E_{K} | −84 | −84 | −84 |

E_{L} | −60 | −60 | −60 |

g_{K} | 8 | 8 | 8 |

g_{L} | 2 | 2 | 2 |

V_{1} | −1.2 | −1.2 | −1.2 |

V_{2} | 18 | 18 | 18 |

C_{M} | 20 | 20 | 20 |

In a supercritical Hopf bifurcation, we have an equilibrium point that one limit cycle bifurcates from it. It means, for the parameter values before the parameter bifurcation, the limit cycle is degenerate and its amplitude is zero, and the amplitude extends as the parameter enhances further. But, in subcritical Hopf bifurcation, we have an unstable limit cycle around the equilibrium point, and a stable limit cycle surrounds it. This unstable limit cycle is as separatrix cycle

that separates two stable states and we could see the bistability that is not available in supercritical Hopf bifurcation and as a result only mono stability occurs in this case. The unstable limit cycle in subcritical Hopf bifurcation, degenerates to the equilibrium point, and the equilibrium points becomes unstable during this bifurcation. This type of bifurcation creates catastrophic changes in the behavior as we have suddenly changed the behavior from stable focus to oscillations with large amplitude [

The bifurcation diagram for model (1) has been displayed in

To rigorous analysis of Hopf bifurcation, as it has been demonstrated in

label = H, x = ( 1.318376 0.488642 207.848139 )

First Lyapunov coefficient = −2.972401e−04

label = H, x = ( − 12.596094 0.274268 160.995457 )

First Lyapunov coefficient = 8.624905e−05

As a result, for I a p p = 160.995457 and ( V , n ) = ( − 12.596094,0.274268 ) we have the first Hopf point with positive First Lyapunov coefficient, and for I a p p = 207.848139 and ( V , n ) = ( 1.318376,0.488642 ) we can see the second Hopf point with negative First Lyapunov coefficient.

Thus, there should exist an unstable limit cycle, bifurcating from the equilibrium and it indicates a subcritical Hopf Bifurcation for the first point and for the second point we have supercritical Hopf bifurcation a very small limit cycle is bifurcated at the bifurcation point. This limit cycle will grow with the further increasing of the injected current.

For the first point, the eigenvalues are ( λ 1 , λ 2 ) = ( 2.70244 e − 08 + i ( 0.176907 ) ,2.70244 e − 08 − i ( 0.176907 ) ) . As we see, we have 2 complex conjugate eigenvalues with positive real part.

We can write the normal form for first Hopf point ( V , n ) = ( − 12.596094,0.274268 ) with the help of eigenvalues and the first Lyapunov coefficient as the form below

r ˙ = ( 2.70244 e − 08 ) r + ( 8.624905 e − 05 ) r 3

θ ˙ = 0.176907

Here, θ ˙ > 0 as we know θ is the angle of oscillations that here is positive and increasing because the frequency of damped or sustained oscillations around this point ω 0 , is positive. But for the following normal form

r ˙ = ( 2.70244 e − 08 ) r + ( 8.624905 e − 05 ) r 3

θ ˙ = − 0.176907

Here, θ ˙ < 0 as we know θ is the angle of oscillations that here is negative and decreasing because the frequency of damped or sustained oscillations around this point ω 0 , is negative.

However, the analysis of normal form is just limited to the first equation of normal form.

( 2.70244 e − 08 ) r + ( 8.624905 e − 05 ) r 3 = 0 → r ( ( 2.70244 e − 08 ) + ( 8.624905 e − 05 ) r 2 ) = 0

So

r = 0 , ( 2.70244 e − 08 ) + ( 8.624905 e − 05 ) r 2 = 0

Here, r = 0 is an equilibrium and because for r = 0

d d r [ ( 2.70244 e − 08 ) r + ( 8.624905 e − 05 ) r 3 ] = ( 2.70244 e − 08 ) > 0

As a result, this equilibrium is unstable. The equation ( 2.70244 e − 08 ) + ( 8.624905 e − 05 ) r 2 = 0 does not give us any periodic solutions or oscillatory behaviors.

For other Hopf point with the eigenvalues ( λ 1 , λ 2 ) = ( − 2.39388 e − 07 + i ( 0.282179 ) , − 2.39388 e − 07 − i ( 0.282179 ) ) we have the following normal form (for ( V , n ) = ( 1.318376,0.488642 ) )

r ˙ = ( − 2.39388 e − 07 ) r − ( − 2.972401 e − 04 ) r 3

θ ˙ = 0.282179

Here, θ ˙ > 0 like above we know θ is the phase of oscillations that here is positive and increasing because the frequency of damped or sustained oscillations around this point, ω 0 , is positive. But for the following normal form

r ˙ = ( − 2.39388 e − 07 ) r − ( − 2.972401 e − 04 ) r 3

θ ˙ = − 0.282179

θ ˙ < 0 as we know θ is the angle of oscillations that here it is negative and decreasing because the frequency of damped or sustained oscillations around this point ω 0 , is negative.

But, we know that the analysis of normal form depends on the first equation of normal form

( − 2.39388 − 07 ) r − ( − 2.972401 e − 04 ) r 3 = 0 → r ( ( − 2.39388 e − 07 ) − ( − 2.972401 e − 04 ) r 2 ) = 0

Therefore

r = 0 , ( − 2.39388 e − 07 ) − ( − 2.972401 e − 04 ) r 2 = 0

Here, r = 0 is an equilibrium and because for r = 0

d d r [ ( − 2.39388 e − 07 ) r − ( − 2.972401 e − 04 ) r 3 ] = ( − 2.39388 e − 07 ) < 0

As a result, this equilibrium is stable. The equation ( − 2.39388 e − 07 ) − ( − 2.972401 e − 04 ) r 2 = 0 gives us the stable periodic solution with amplitude

r = 2.39388 e − 07 2.972401 e − 04

In this section, we consider g i n h and I a p p as free parameters.

If we continue to increase I a p p , the period of limit cycle becomes 2-times and then 2^{2} and … Here we say that the system undergoes the period-doubling bifurcation. As we do continuation of Hopf point, the trajectory repeats itself while we passed through PD point and we can see the trajectory is going around the loop twice before coming back to PD point. For period-doubling bifurcation a new limit cycle creates from the primary limit cycle, and the period of this new limit cycle is two times the previous limit cycle (see

If we increase I a p p further, the limit cycle that is born from supercritical Hopf bifurcation undergoes another period-doubling bifurcation. The Period-doubling can continue to infinity. But, in this procedure, the values of I a p p at which this bifurcation occurs become closer than before.

For localization of Branch Points of cycles, a special algorithm which is of Newton type Continuation Method has been used to show these singularities. Also, the equilibrium curve can have two generic codimension one bifurcation limit point and Hopf point and plus continuation of equilibrium curve gives us Branch Point that in

The branch point Z 0 for the complex valued function f ( z ) , has this property that the value of f ( z ) does not return to its initial value as a closed curve around a point is traced, in such a way that f varies continuously as the pass is traced.

Moreover, In the contex of complex analysis, a branch point of a multi-valued function is a point at which the function is discontinuous when it is going in a neighborhood of this point.

If we only consider g i n h as bifurcation parameter, we have

label = H, x = ( − 19.573692 0.191817 0.071067 )

Neutral saddle

label = LP, x = ( − 4.616109 0.391485 0.187732 )

a = 8.199830 e − 03

label = H, x = ( 6.041592 0.566955 − 0.192795 )

First Lyapunov coefficient = 7.969763e−04

In this case, we have a natural saddle in ( V , n ) = ( − 19.573692,0.191817 ) and g i n h = 0.071067 for which λ 1 + λ 2 = 0 .

Also for ( V , n ) = ( 6.041592,0.566955 ) and g i n h = − 0.192795 and for the eigenvalues ( λ 1 , λ 2 ) = ( − 1.75069 e − 07 + i ( 0.291426 ) , − 1.75069 e − 07 − i ( 0.291426 ) ) we have a subcritical Hopf bifurcation with the following normal form

r ˙ = ( − 1.75069 e − 07 ) r + ( 7.969763 e − 04 ) r 3

θ ˙ = 0.291426

Here, θ ˙ > 0 as we know θ is the angle of oscillations that here is positive and increasing because the frequency of damped or sustained oscillations around this point ω 0 , is positive. But for the following normal form

r ˙ = ( − 1.75069 e − 07 ) r + ( 7.969763 e − 04 ) r 3

θ ˙ = − 0.291426

Here, θ ˙ < 0 as we know θ is the angle of oscillations that here is negative and decreasing because the frequency of damped or sustained oscillations around this point ω 0 , is negative.

However, the analysis of normal form is just limited to the first equation of normal form.

( − 1.75069 e − 07 ) r + ( 7.969763 e − 04 ) r 3 = 0 → r ( ( − 1.75069 e − 07 ) + ( 7.969763 e − 04 ) r 2 ) = 0

So

r = 0 , ( − 1.75069 e − 07 ) + ( 7.969763 e − 04 ) r 2 = 0

Here, r = 0 is an equilibrium and because for r = 0

d d r [ ( − 1.75069 e − 07 ) r + ( 7.969763 e − 04 ) r 3 ] = ( − 1.75069 e − 07 ) < 0

As a result, this equilibrium is stable. The equation ( − 1.75069 e − 07 ) + ( 7.969763 e − 04 ) r 2 = 0 gives us the stable periodic solution with amplitude

r = 1.75069 e − 07 7.969763 e − 04

The complete bifurcation diagram with respect to g i n h has been demonstrated in

In

u ˙ = v

v ˙ = a + b u + u 2 + σ u v

where, a , b are the normal form coefficients, and the parameter σ = 1 , − 1 depends on different types of supercritical and subcritical.

In another word, when cusp bifurcation happens, we have a Saddle-Node bifurcation at an equilibrium point, when for the system x ˙ = f ( x , γ ) we have f x = 0 and f x x ≠ 0 and this system is equivalent with x ˙ = a ( γ ) + x 2 .

When we have f x x = 0 and f x x x ≠ 0 , then the equilibrium point undergoes cusp bifurcation that is a codimension two bifurcation with respect to I a p p , g i n h as bifurcation parameters. The system near the cusp point obeys the following normal form

x ˙ = a 1 ( γ ) + a 2 ( γ ) x + c x 3

So that

a 1 ( γ ) = f ( x γ , γ ) , a 2 ( γ ) = f x ( x γ , γ ) , c = f x x x / 6 ≠ 0

If c > 0 we have subcritical cusp bifurcation and if c < 0 , we have supercritical cusp bifurcation. To find the bifurcation set:

d d x ( a 1 ( γ ) + a 2 ( γ ) x + c x 3 ) = a 2 ( γ ) + 3 c x 2 = 0

which gives the following saddle-node curves

a 1 ( γ ) = + 2 a ( a 2 ( γ ) 3 ) 3 / 2 , a 2 ( γ ) = − 2 a ( a 2 ( γ ) 3 ) 3 / 2

Also, we know that if a 1 ( γ ) = 0 and a 2 ( γ ) = γ , we have the following normal form

x ˙ = γ x + c x 3

which is the normal form for Pitchfork bifurcation.

Here, for the first point

label = CP, x = ( − 10.008656 0.309902 221.927487 0.996196 )

c = 5.411489 e − 04

( a 1 ( γ ) , a 2 ( γ ) ) = ( I a p p , g i n h ) = ( 221.927487 , 0.996196 )

So that we can easily write the normal form for this equilibrium point and we can find the bifurcation set and saddle-node curves.

In this section, we consider g e x c and I a p p as free parameters.

At first, we assume g e x c as bifurcation parameter (see

In

label = H, x = ( − 7.210500 0.351140 1.975854 )

First Lyapunov coefficient = −9.396562e−04

label = H, x = ( − 12.596098 0.274268 0.209999 )

First Lyapunov coefficient = 8.625008e−05

For the first Hopf point ( V , n ) = ( − 7.210500,0.351140 ) and g e x c = 1.975854 and the eigenvalues ( λ 1 , λ 2 ) = ( − 3.66907 e − 08 + i ( 0.232446 ) , − 3.66907 e − 08 − i ( 0.232446 ) ) we have a supercritical Hopf bifurcation with the following normal form

r ˙ = ( − 3.66907 e − 08 ) r − ( − 9.396562 e − 04 ) r 3

θ ˙ = 0.232446

Here, θ ˙ > 0 as we know θ is the angle of oscillations that here is positive and increasing because the frequency of damped or sustained oscillations around this point ω 0 , is positive. But for the following normal form

r ˙ = ( − 3.66907 e − 08 ) r − ( − 9.396562 e − 04 ) r 3

θ ˙ = − 0.232446

Here, θ ˙ < 0 as we know θ is the angle of oscillations that here is negative and decreasing because the frequency of damped or sustained oscillations around this point ω 0 , is negative.

However, the analysis of normal form is just limited to the first equation of normal form.

( − 3.66907 e − 08 ) r − ( − 9.396562 e − 04 ) r 3 = 0 → r ( ( − 3.66907 e − 08 ) − ( − 9.396562 e − 04 ) r 2 ) = 0

So

r = 0 , ( − 3.66907 e − 08 ) − ( − 9.396562 e − 04 ) r 2 = 0

Here, r = 0 is an equilibrium and because for r = 0

d d r [ ( − 3.66907 e − 08 ) r − ( − 9.396562 e − 04 ) r 3 ] = ( − 3.66907 e − 08 ) < 0

As a result, this equilibrium is stable. The equation ( − 3.66907 e − 08 ) − ( − 9.396562 e − 04 ) r 2 = 0 gives us the stable periodic solution with amplitude

r = 3.66907 e − 08 9.396562 e − 04

For other Hopf point with the eigenvalues ( λ 1 , λ 2 ) = ( 2.6594 e − 11 + i ( 0.176907 ) ,2.6594 e − 11 + i ( 0.176907 ) ) we have the following normal form (for ( V , n ) = ( − 12.596098,0.274268 ) and g e x c = 0.209999 )

r ˙ = ( 2.6594 e − 11 ) r + ( 8.625008 e − 05 ) r 3

θ ˙ = 0.176907

Here, θ ˙ > 0 like above we know θ is the phase of oscillations that here is positive and increasing because the frequency of damped or sustained oscillations around this point, ω 0 , is positive. But for the following normal form

r ˙ = ( 2.6594 e − 11 ) r + ( 8.625008 e − 05 ) r 3

θ ˙ = − 0.176907

θ ˙ < 0 as we know θ is the angle of oscillations that here it is negative and decreasing because the frequency of damped or sustained oscillations around this point ω 0 , is negative.

But, we know that the analysis of normal form depends on the first equation of normal form

( 2.6594 e − 11 ) r + ( 8.625008 e − 05 ) r 3 = 0 → r ( ( 2.6594 e − 11 ) + ( 8.625008 e − 05 ) r 2 ) = 0

Therefore

r = 0, ( 2.6594 e − 11 ) + ( 8.625008 e − 05 ) r 2 = 0

Here, r = 0 is an equilibrium and because for r = 0

d d r [ ( 2.6594 e − 11 ) r + ( 8.625008 e − 05 ) r 3 ] = ( 2.6594 e − 11 ) > 0

As a result, this equilibrium is unstable. The equation ( 2.6594 e − 11 ) + ( 8.625008 e − 05 ) r 2 = 0 does not give us the periodic solutions.

Also, at g i n h = 0.071067341 and considering g e x c as bifurcation, parameter with changing the period we have the result in

In this section, we study the homoclinic bifurcation using the parameter values from

When a saddle homoclinic bifurcation occurs a limit cycle disappears or born. A saddle point has two stable and unstable submanifolds. When we change the parameters for bifurcation, these two submanifolds collide and create one homoclinic orbit as we see clearly in

In neuroscience point of view, before the beginning of bifurcation, all trajectories go back to the resting state after a counter-clockwise excursion (firing spike) and we have limit cycle attractors corresponding to periodic spiking. At the moment of bifurcation, we have periodic spiking activity with the period that goes to infinity (the period of limit cycle goes to infinity). After the bifurcation the state of the system again goes back to the rest and we do not have tonic spiking.

There is another bifurcation similar to saddle Homoclinic bifurcation and in both of them we have a Homoclinic orbit that begins and finishes in the same equilibrium point. But saddle-Homoclinic is of codimension 1 and saddle-node Homoclinic is of codimension two. For saddle-node Homoclinic we must have the condition for occurrence of saddle-node bifurcation λ 1 = 0 and also the saddle-node quantity λ 1 + λ 2 < 0 it means that λ 2 < 0 . As a result, in saddle-node Homoclinic, we observe the appearance or disappearance of a stable limit cycle.

In this study we discovered the biological and spiking behaviors of Morris-Lecal neuron model including a new more channel, chloride channel, and we compared it to the original model. Morris-Lecal neuron model with two ordinary differential equations is a reduction version of Hodgkin-Huxley model. We have shown that for different biological parameters values, the model displays quiescent and spiking behaviors. To explore the exciting behaviors of neuron, we conducted different numerical simulations and bifurcation analysis. We obtained numerically the Hopf bifurcation and homocinic bifurcation at larger values for injected current compared to the original model and we derived them analytically by finding their normal forms. Specifically, we conducted rigorous bifurcation analysis and continuation of equilibrium point with respect to the injected current or changing other biological parameters and we explored the complicated and interesting dynamics of new Morris-Lecar model such as different fluctuations neuronal patterns.

In supercritical Hopf bifurcation, when we changed the control parameters (bifurcation parameters), we could see a stable equilibrium point lost its stability and a stable limit cycle bifurcated from it. In this case, the state of the system stayed near the equilibrium and we could see only small amplitude oscillations around the equilibrium point. When we changed the parameter in the opposite direction the limit cycle degenerated to a point and we reached the previous equilibrium point. We noticed that as it is expected this neuron that has undertaken supercritical Hopf bifurcations did not show an immediate spike and it remained quiescent (it displayed small amplitude sustained oscillations).

In other cases, for subcritical Hopf bifurcation, when we changed the bifurcation parameter, we noticed the immediate spike or a large amplitude jump in time series of system equations. In this case, when we changed the parameter in opposite direction, the equilibrium became stable again but the state of the system went back to its similar saddle-node bifurcation or saddle-Homoclinic bifurcations.

As we have seen for model (1), there are different types of codimension-two bifurcations which occurred when we had two bifurcation parameters. One of them was the Bogdanov-Takens bifurcation which is a codimension-two bifurcation of an equilibrium point in a two-dimensional system for which the equilibrium point has two zero eigenvalues. Near the bifurcation parameter, we have two equilibrium points, a saddle point and others that are not a saddle point which collides and disappears via a saddle-node bifurcation. For the equilibrium that was not saddle, we have seen Hopf bifurcation that produced a limit cycle. This single limit cycle degenerated to a Homoclinic orbit to the saddle and disappeared via a saddle homoclinic bifurcation. Bogdanov-Takens bifurcation can occur when an equilibrium undergoes Hopf bifurcation and saddle-node bifurcation simultaneously and also it occurs when we have two-dimensional systems. Therefore, for this model we showed that we have a transition between integrator and resonator as the neuron that undergoes saddle-node bifurcation is an integrator (without damped subthreshold oscillations) and the neuron that undergoes Hopf bifurcation is a resonator (with damped subthreshold oscillations).

Continuation of limit cycles from the Hopf point, gave us the limit point (Fold/Saddle Node) bifurcation. Limit point cycles are a set of limit cycles that bifurcate from the Hopf point, and moreover, limit point cycle (LPC) is a fold bifurcation for which two limit cycles with different periods are near the LPC point. Likewise, as we saw, the continuation of cycles from a Hopf point also gave period doubling (flip) bifurcations. Period doubling bifurcation is a bifurcation for which a new limit cycle creates from the limit cycle that exists and the period of the new limit cycle is twice the first limit cycle. Moreover, we have shown changing the type of Hopf bifurcation from supercritical to subcritical which means that the type of stability changes from monostability in supercritical to bistability in subcritical Hopf bifurcation.

The authors declare no conflicts of interest regarding the publication of this paper.

Azizi, T. and Alali, B. (2020) Impact of Chloride Channel on Spiking Patterns of Morris-Lecar Model. Applied Mathematics, 11, 650-669. https://doi.org/10.4236/am.2020.117044