Qualitative Analysis of a Quadratic Integrate-and-Fire Neuron Model with State-Dependent Feedback Control

Spiking neuron models which exhibit rich dynamics are usually defined by hybrid dynamical systems. It is revealed that mathematical analysis of thesemodels has important significance.Therefore, in this work, we provide a comprehensively qualitative analysis for a quadratic integrate-and-fire model by using the theories of hybrid dynamical system. Firstly, the exact impulsive and phase sets are defined according to the phase portraits of the proposed model, and then the Poincaré map is constructed. Furthermore, the conditions for the existence and stability of an order 1 periodic solution are provided. Moreover, the existence and nonexistence of an order k (k ≥ 2) periodic solution have been studied theoretically and numerically, and the results show that the system has periodic solutions with any period. Finally, some biological implications of the mathematical results are discussed.


Introduction
To understand the working mechanism of the brain, it is necessary to combine experimental studies of nervous systems with numerical simulations of neuron models [1].Therefore, a large number of differential equations have been proposed to model neurons in the neurocomputing community.However, there are two major impediments in computational neuroscience including the need for a computationally simple and biologically realistic model of the neuron [2].To conquer these issues, FitzHugh and Nagumo introduced a new neuron model with cubic nonlinearity and two variables [3,4].Then this novel model was reduced to a class of nonlinear bidimensional spiking models with adaptation by many authors which are mathematically tractable, are efficiently implemented, and are able to reproduce many electrophysiological signatures [5][6][7].These models often consist of two variables: the membrane potential of the neuron V() and the adaptation variable ().In this paper, the leak and spike initiation currents function (V) is assumed to be V 2 ; then the quadratic adaptive model can be described by the following ordinary differential equations: where  represents the input intensity of the neuron, 1/ ( > 0) is the characteristic time of the adaptation variable, and  accounts for the interaction between the membrane potential and the adaptation variable.System (1) has been recently used by Izhikevich and Edelman in very large-scale simulations of neural networks [7].In addition, the rich dynamics that system (1) presented have been investigated in detail by many authors [2,8,9].In view of the spike and reset condition, a spike is emitted at time  when the membrane potential V() reaches a threshold (or a cutoff value) V  .Then the membrane potential is reset to a constant value  and the adaptation variable is updated to () + .Therefore, system (1) combining with the reset process can be written as V ()  =  + V 2 () −  () ,  ()  =  (V () −  ()) , V () < V  , V ( + ) = ,  ( + ) =  () + , where  ( > 0) describes the total amount of outward minus inward currents activated during the spike and  is the reset membrane potential.System (2) is a hybrid dynamical system which is also known as a planar impulsive semidynamical system (ISDS) [10,11].Recently, Shlizerman and Holmes have studied system (2) and obtained explicit approximations of instantaneous firing rates for fixed values of the recovery variable, and then they used the averaging theorem to obtain asymptotic firing rates as a function of current and reset parameters, which provided explicit tools for the interpretation of different spiking patterns [8], whereas Touboul and Brette constructed an adaptation map to study the patterns of spikes numerically when the adaptation variable blows up [2].In this paper, the global dynamical behaviors of system (2) will be studied theoretically.To this end, some useful definitions and lemmas of the planar ISDS will be presented in Section 2; then we will briefly review the rich dynamics of system (1).In Section 3, we present our main results.In particular, the impulsive and phase sets will be defined exactly in consideration of the phase portraits of model (1) in Section 3.1.In Section 3.2, the Poincaré map for impulsive point series defined in the exact phase set is obtained.In Section 3.3, the existence and stability of order 1 periodic solution will be addressed in detail.In Section 3.4, the existence of order  ( ≥ 2) periodic solutions will be studied theoretically and numerically.In the last section, Section 4, the related biological implications are discussed from a neurocomputational point of view.

Preliminaries and Main Properties of System (1)
The generalized planar ISDS with state-dependent feedback control can be described as follows: (, ) ∉ I, where (, ) ∈  2 .We denote  + = ( + ) and  + = ( + ) for simplicity, and , ,  1 , and  2 are continuous functions from  2 into , and I ⊂  2 denotes the impulsive set.For each point (, ) ∈ I, the map  :  2 →  2 is defined: and  + is called an impulsive point of .
Let P = (I) be the phase set (i.e., for any  ∈ I, () =  + ∈ P), and P ∩ I = Ø.System (3) is generally known as a planar ISDS.We note that system (2) is an ISDS, where impulsive set 2 .It follows that the phase set Without loss of generality, unless otherwise specified we assume the initial point (V + 0 ,  + 0 ) ∈ P. In the following we briefly list some definitions related to ISDS, which are useful in this work.
Let (, Π, ) or (, Π) be a semidynamical system [10,11], where  is a metric space and  is the set of all nonnegative reals.For any  ∈ , the function Π :  →  is a continuous function such that Π(, 0) =  for all  ∈ , and Π(Π(, ), ) = Π(,  + ) for all  ∈  and ,  ∈  (denoted by Π  () = Π(, )).Throughout the paper, we denote the points of discontinuity of Π  by { +  } and call  +  an impulsive point of   .We define a function Φ from  into the extended positive reals  ∪ {∞} as follows: let  ∈ ; if I + () = Ø we set Φ() = ∞; otherwise I + () ̸ = Ø and we set Φ() = , where Π(, ) ∉ I for 0 <  <  but Π(, ) ∈ I. Definition 2. A trajectory Π  in (, Π, I, ) is said to be periodic of period   and order  if there exist nonnegative integers  ≥ 0 and  ≥ 1 such that  is the smallest integer for which For simplification, we denote a periodic trajectory of period   and order  by an order  periodic solution.An order  periodic solution is called an order  limit cycle if it is isolated.For more details of the concepts and properties of continuous dynamical systems and impulsive dynamical systems, see [12,15,17].
Lemma 3 (Analogue of Poincaré Criterion [10,11]).The periodic solution (, ) = ((), ()) of system Since system (1) has been investigated by many scholars, bifurcations and phase portraits were also addressed in detail [2,8,9].Here we briefly review the dynamics of system (1) which are useful in this study.The two isoclines of system (1) are denoted by  1 and  2 , where  1 :  =  + V 2 ,  2 :  = V.If we fix parameters  and , choosing  as a variable, then the following results can be obtained easily.
In order to address the global dynamics of system (2), we consider case ( =  hom ) at which there exists a homoclinic orbit based on the qualitative behaviors of ODE model (1).In this case the system has two fixed points  1 and  2 , and the solution initiating from the inside of the homoclinic orbit may not reach the threshold V  under certain conditions (see more details later), which means that the impulsive effects may not happen.This indicates that system (2) may exhibit rich dynamics due to different values of the threshold and the dynamics presented in this case could be more complicated than other cases.Therefore, if we can investigate the global dynamics of this case in detail by using the theories of the ISDS and analytic techniques, then by using similar methods the dynamics of the remaining cases could be studied.

Impulsive Set and Phase Set.
In order to investigate the existence of order  periodic solutions, the Poincaré map is constructed first.Nevertheless, it is essential to know that the exact conditions under which the solution of system (2) starting from (V + 0 ,  + 0 ) ∈ P is free from impulsive effects, that is to say, the more exact phase set P, should be provided.Moreover, the impulsive set I defined in Section 2 is the maximum interval for the vertical coordinates.Therefore, the part of I which the solution of system (2) can not reach will be removed and then the exact domains of impulsive set can be obtained.From now on, unless otherwise specified we assume that  =  hom ( > ) throughout the paper.
Based on different positions of V  we consider the following three cases: The homoclinic cycle is denoted by Γ ℎ , and we let V 4 be the horizontal component of the small intersection point of the homoclinic cycle Γ ℎ with the line  1 (denoted by  4 = (V 4 ,  4 ); see Figure 1).For case ( 1 ), V  intersects with the right branch of the homoclinic cycle at two points and we denote the lower point by (V  ,   ).According to the relations among , V 1 , and V 4 , there may be three cases for : then there exists a trajectory Γ 1 such that the line V =  tangents to the curve Γ 1 at the point (,   ); here   =  +  2 .Clearly, the curve Γ 1 also intersects the line V = V  at two points and we denote the lower point by (V  ,   ).Therefore, if  < V 4 , then the impulsive set is defined by If so, any trajectory of system (2) initiating from the interior of the curve Γ 1 either cannot reach Table 1: Exact domains of impulsive set and phase set of system (2).
the impulsive set or reaches the interval (  , V 2  + ] and then maps to the phase set P 1 after one time impulsive effect.If V 4 ≤  ≤ V 1 , then the line V =  intersects the left branch of homoclinic cycle Γ ℎ at two points, denoted by  1 = (,   1 ) and  2 = (,   2 ).Thus, the impulsive set can be defined as Moreover, the trajectories initiating from (V + ,  + ) ∈ P with   1 ≤  + ≤   2 will be free from the impulsive effects and will tend to the stable equilibrium  2 .If  > V 1 , then by using the same methods as subcase  < V 4 the impulsive set is defined by I 1 , and the phase set is defined by P 1 .
For case ( 2 ) (as shown in Figure 2 3 + ) such that Γ 2 tangents to the line  1 at this point.If  ≤ V 3 , then the curve Γ 1 which tangents to the line V =  at the point  must intersect the line V = V  at two points; the lower intersect point is denoted by .Therefore, the impulsive set is defined by I 1 , and the phase set is defined by P 1 .If  > V 3 , then the curve Γ 2 must intersect the line V =  at two points, denoted by  1 = (,   1 ) and  2 = (,   2 ).Thus, the impulsive set can be defined as I and the phase set is defined by Moreover, any trajectories initiating from (V + ,  + ) ∈ P with   1 <  + <   2 will be free from the impulsive effects and will tend to the stable equilibrium  2 .
For case ( 3 ) (as shown in Figure 2(b)), by using the same methods as cases ( 1 ) and ( 2 ), the impulsive set and the phase set can be defined as I 1 and P 1 , respectively.
In conclusion, we list all possible cases for the domains of impulsive sets and phase sets of model ( 2) in Table 1 based on the above analyses.It is found that the basic phase set P can not be used to define the real phase set of model (2) for any case, which indicates that the exact domains of phase set of model ( 2) should be carefully discussed before carrying out qualitative analyses.
In order to provide the exact domains of the Poincaré map, we first need to know under what conditions the trajectory initiating from  +  ∈ P  can not reach the impulsive set.From the analysis of Section 3.1, there may be two possible cases.
Case 1.The first case is For Case 1, it follows from Section 3.1 that if the initial point  +  lies on the homoclinic cycle Γ ℎ or its interior, then the trajectory can not reach the line V = V  forever, which indicates that this trajectory is free from impulsive effects, as shown in Figure 1, so we need  +  ∈ P 2 .For Case 2, if the initial point  +  lies in the interior of the curve Γ 2 , then the trajectory can not reach the line V = V  and will tend to the stable equilibrium  2 , which means that this trajectory is also free from impulsive effects, and then we need  +  ∈ P 3 , as shown in Figure 2(b).
In addition to the above two cases, any trajectory of system (2) initiating from  +  ∈ P  can reach at the point  +1 ∈ I. Due to the impulsive effects, we have To explore the dynamical behaviors of system (2), the Poincaré map is constructed in the exact domains of the phase set.
Lemma 6 (see [18]).The Poincaré map for the impulsive points of model ( 2) defined in phase set can be defined as where the Poincaré map ( +  ) is defined by (9).
Since the function  and consequently the function  are continuously differentiable with respect to  +  , for each fixed point of Poincaré map , there exists an associated periodic solution for system (2).

Existence and Stability of Order 1 Periodic Solution.
In this subsection, we mainly focus on investigating the existence and stability of an order 1 periodic solution of system (2) based on cases ( 1 ), ( 2 ), and ( 3 ), which could be realized by studying the existence and stability of the fixed point of the Poincaré map.To this end, we have following main results.
Figure 3: Location of the periodic solution of system (2) for case The case when the pulsed point  + is above the point ; (b) the case when the pulsed point  + lies below the point .
Poincaré map  exists and consequently there exists an order 1 periodic solution for system (2).
Proof.Firstly, we prove that there exists a fixed point for the Poincaré map  if  < V 4 or  > V 1 for case ( 1 ).Then by using the same methods, it can be demonstrated that the results are also true for case ( 2 ) when  ≤ V 3 or for case ( 3 ) when  < V 2 .
For the first case, the curve Γ 1 initiating from the point (,   ) tangents to the line V =  and reaches the line V = V  at the point (V  ,   ) (clearly,   = (  )), and then the point  maps to the point  + (,   + ) (or  + (,   + )) due to the impulsive effect.To determine the location of the point  + , we assume that there exists  0 such that   +  0 =  +  2 ; then the point  + coincides with the point  for  =  0 .Obviously, the Poincaré map  has a fixed point and consequently the curve P forms an order 1 periodic solution for system (2).
If  >  0 (see Figure 3(a)), then the pulsed point  + is above the point .Thus, the following inequality holds: Furthermore, the trajectory initiating from the point  + will meet the line V  at a point, denoted by  1 (V  , (  + )).In view of disjointness of any two trajectories, the point  1 lies below the point .Then the point  1 maps to the point  + 1 (,   + 1 ) after one time impulsive effect.Besides, the two line segments  + and  1  + 1 satisfy  + ‖  1  + 1 , which indicates that the point  + 1 must lie below the point  + .So another inequality holds: It follows from ( 13) and ( 14) that the Poincaré map has a fixed point, which corresponds to an order 1 periodic solution of system (2).If  <  0 , then the pulsed point  + lies below the point  (see Figure 3(b)).Thus, the trajectory starting from the point  + will reach the line V = V  at the point  1 , and it is easy to see that the point  1 lies below the point  because of the uniqueness of solutions.Then the point  1 maps to the point  + 1 after one time impulsive effect, where the point  + 1 lies below the point  + .By induction, it is noted that the point  + +1 lies below the point  +  ,  = 1, 2, . . ., , which indicates that the impulsive point series  +  is monotonically decreasing.That is, Therefore, what we want to show in the following is that the impulsive point series   +  is bounded.First of all, V/ > 0 for V ∈ [, V  ]; that is to say, V is monotonically increasing with respect to .This indicates that any solution initiating from From the second equation of system (1) we have According to Comparison Theorem of ODE [19] we can get Figure 4: Location of the periodic solution of system (2) for case According to (17), it is known that the pulsed sequence   +  (or Poincaré map) is bounded and monotonically decreasing.This indicates that the sequence   +  has a limit point.Therefore, the Poincaré map has a fixed point when  < V 4 or  > V 1 for case ( 1 ), which corresponds to an order 1 periodic solution of system (2).This completes the proof.Remark 8.For case ( 1 ) when  > V 1 , there exists a unique homoclinic cycle in system (2).In particular, there exists  * such that the point  maps to the point   which lies on the stable manifold of the right branch of homoclinic cycle Γ ℎ in the phase set; it is easy to see that the curve Ĥ  1  forms a unique homoclinic cycle with impulsive effects.Moreover, if  >  * or  <  * , then the homoclinic cycle is broken and an order 1 periodic solution is generated for system (2).
where   + 1 is the vertical coordinate of the initial point (,   2 ) from which the solution experiences two time impulsive effects, then there exists a fixed point for the Poincaré map  which corresponds to an order 1 periodic solution for system (2).
Proof.For this case (as shown in Figure 4), the line V =  intersects the left branch of homoclinic cycle Γ ℎ at two points  1 = (,   1 ) and  2 = (,   2 ).In the region above the line  2 , we have ()/ < 0; that is to say, the adaptation variable () is decreasing in this area.While in the region lying below the line  2 one has ()/ > 0, this indicates that the adaptation variable () is increasing.Therefore, the following inequality must hold true: Due to impulsive effects, the point  maps to the point  + which is above the point  2 due to  > 0 and   >   2 .Thus, the trajectory initiating from the point  + will meet the line V  at two points, the lower point denoted by  1 (V  ,   1 ).The point  1 must lie below the point  in view of disjointness of any two trajectories.Then the point  1 maps to the point  + 1 (,   + 1 ) after one time impulsive effect.For the location of the point  + 1 , it follows from the assumption that the point  + 1 must lie above the point  2 and lies below the point  + .Thus, we have following inequality: Similarly, the solution starting from the point  + 1 meets the line V = V  at the point  2 which lies between  and  1 , and then  2 maps to the point  + 2 which lies between  + and  + 1 .Thus, the following inequality holds: It follows from ( 19) and ( 20) that there exists a fixed point for the Poincaré map, which is also an order 1 periodic solution of system (2).This completes the proof.) from which the solution experiences two time impulsive effects, then the fixed point of the Poincaré map  exists and consequently there exists an order 1 periodic solution for system (2).
Proof.Assuming that the first intersection point of the line  2 with the curve Γ 2 (initiating from the point  2 ) is denoted as  5 (V 5 ,  5 ), see Figure 2(a).Hence, if  ≥ V 5 , then the point  is above the point  1 because the adaptation variable () is increasing in the region that lies below the line  2 .Further, the pulsed point  + is above the point  1 after one time impulsive effect due to  > 0. Therefore, if   + ≥   2 and   + 1 ≥   2 , then, by using the same method as the proof of Theorem 9, it can be proved that there exists a fixed point for the Poincaré map .
If V 3 <  < V 5 , then there are two possible cases for the relations of   1 and   :   >   1 or   ≤   1 .For the former, it is clear that the point  is still above the point  1 .From the above analysis, there exists a fixed point for the Poincaré map.For the latter, if   + ≥   2 and   + 1 ≥   2 , then there exists a fixed point for the Poincaré map; if   + ≤   1 , then, by using the same analytic techniques as shown in Theorem 7 for the case  ≤  0 , it can be shown that there exists a fixed point for the Poincaré map , which means that system (2) has an order 1 periodic solution.This completes the proof.
To sum up, the conditions for the existence of an order 1 periodic solution of system (2) are provided for cases ( 1 ), ( 2 ), and ( 3 ).In addition, parameter sets which guarantee the existence of a unique homoclinic cycle and phasic spiking are also obtained.Furthermore, in the following the stability of the order 1 periodic solution will be investigated by using the Analogue of Poincaré Criterion which is introduced in Section 2.
Proof.We first show that the order 1 periodic solution which is obtained in Theorem 7 is orbitally asymptotically stable.
To show this, assume that the order 1 periodic solution with period  passes through the points (V  ,  0 ) and  + (,  0 +).
Note that the explicit expression and period  of the order 1 periodic solution are unknown.Thus, Lemma 3 provides a necessary tool to discuss stability of this periodic solution.To do this, we denote Then Moreover, we have Therefore, From the proof of Theorem 7, there exists  0 such that the point  + coincides with the point  when  0 satisfies   +  0 =  +  2 .Hence, | 2 | = 0 < 1 leads to the fact that the order 1 periodic solution is orbitally asymptotically stable.
Therefore, if | 2 | < 1, that is, inequality (21) holds true, then the order 1 periodic solution ((), ()) is orbitally asymptotically stable and enjoys the property of asymptotic phase.Similarly, we can prove that the order 1 periodic solution ((), ()) obtained from Theorems 9 and 10 is orbitally asymptotically stable as well.This completes the proof.

Existence and Nonexistence of Order 𝑘 Periodic Solutions.
Our main purpose in this subsection is to study the existence and nonexistence of order  ( ≥ 2) periodic solutions of system (2).Theorem 13.Assume that  < V 4 or  > V 1 for case ( 1 ) (or  ≤ V 3 for case ( 2 ) or  < V 2 for case ( 3 )); if  +  < + 2 , then system (2) does not have an order  ( ≥ 2) periodic solution other than an order 1 periodic solution.
Proof.If  +  <  +  2 , then all the impulsive points initiating from phase set P 1 lie below the point .From the proof of Theorem 7, the impulsive point series  +  satisfies (15).Besides, the pulsed segments  + ,  1  + 1 , . . .,    +  are parallel, that is, It follows from ( 15) that an order  ( ≥ 2) periodic solution does not exist for system (2).This completes the proof. we have the following two possibilities ( = 1, 2, . . ., ; without loss of generality, we assume that there exists a positive integer  such that  = 2 + 1): From (a) and (b), we conclude that the impulsive point series { +  2 } is decreasing in the interval (  2 ,   + ] and { +  2+1 } is increasing in the interval (  2 ,   + ], and they both converge to a fixed point in this bounded region.Moreover, note that all the impulsive lines are parallel.Therefore, the Poincaré map  has either a fixed point or a period two-point cycle; that is to say, system (2) does not have an order  ( ≥ 3) periodic solution if   + 1 >   2 .This completes the proof.
From Theorem 14, the conditions for the nonexistence of an order  ( ≥ 3) periodic solution are provided for case >   2 .By using the same methods, it is also confirmed that system (2) does not have an order  ( ≥ 3) periodic solution for another three cases: case ( 1 ) with  < V 4 (or  > V 1 ) and  +  1 ≥  + Further, if  → ∞, then, by using the same methods as those shown in Theorem 13, there is no order  ( ≥ 2) periodic solution other than an order 1 periodic solution in system (2).However, if there exists  <  such that  + −1 ≤  +  2 and  +  >  +  2 , then we cannot determine the accurate position of  + +1 , which indicates that we can not determine the nonexistence of order  ( ≥ 3) periodic solution in this case.In fact, the kinetic behaviors of system (2) displayed in this case are very complicated and it is difficult to address the qualitative behaviors analytically.Nevertheless, a traditional approach to gain preliminary insight into the properties of a dynamic system is to carry out a one-dimensional bifurcation analysis.One-dimensional bifurcation diagrams show the dependence of the dynamic behaviors on a certain parameter.Therefore, we resort to numerical investigations to show what dynamics will be presented in system (2).To show this, we set the threshold V  = 3 and fix all parameters from [8] as follows:  hom = 0.1485,  = 0.5,  = 1, and  = 2.5 and then choosing  as a bifurcation parameter.
Figure 5 illustrates the bifurcation diagrams for system (2) with  as the control parameter.Figure 5(b) is the magnified part of Figure 5(a).For parameter  > 5.93, a stable order 1 periodic solution is observed.As parameter  is reduced, the order 1 periodic solution loses its stability and then there exists an order 2 periodic solution for system (2) via a period-doubling bifurcation.As  further decreased about  ≈ 3, it can be seen that system (2) exhibits transitions from order  periodic solutions to order ( + 1) periodic solutions ( = 2, 3, 4, 5) via period-adding bifurcations involving chaotic bands.When parameter  is slightly less than 1.182, the periodic attractor suddenly disappears and the chaotic attractor abruptly appears, thus constituting a type of attractor crisis (the phenomenon of crisis shows that chaotic attractors can suddenly appear or disappear or change size discontinuously as a parameter smoothly varies [20]) (Figure 6).When 1.145 <  < 1.182, the dynamics presented in model ( 2) is very complicated, including chaotic bands, periodic windows, period-double bifurcations, and periodhalving bifurcations (for details see Figure 5(b)).With further reduction of , period-halving bifurcations lead system (2) into a stable state of an order 1 periodic solution.

Conclusion
The quadratic integrate-and-fire model with a recovery variable has been recently used in very large-scale simulations of neural networks; it is stressed that this simple model is capable of reproducing the rich behavior of biological neurons, such as spiking, bursting, and mixed mode firing patterns and continuous spiking with frequency adaptation [1].Consequently, mathematical analysis of the model not only has important theoretical significance, but also has direct biological relevance.Recently, the firing rates for fixed values of the recovery variable were studied by Shlizerman and Holmes [8].Then the subthreshold dynamics were discussed in detail by Touboul [21].Subsequently, the patterns of spikes were addressed numerically when the adaptation variable blows up [2,22].Here the main purpose is to provide a comprehensively qualitative analysis for all possible dynamics on whole parameter space for system (2) by using the theories of impulsive semidynamic system, and those theories have been applied broadly in many fields; for details see [23][24][25][26][27].
This paper mainly focuses on the most complex case  =  hom since other cases are easier and can be studied analogously.First of all, the impulsive and phase sets are defined exactly in line with the phase portraits of model (1), and then the Poincaré map is constructed in the exact phase sets.Further, the different parameter spaces for the existence and stability of order 1 periodic solution are investigated in detail by using the Poincaré map and the Analogue of Poincaré Criterion.In addition, the conditions for the homoclinic cycle with impulsive effects and phasic spiking are provided.Moreover, the existence and nonexistence of order  ( ≥ 2) periodic solution are studied theoretically and numerically.
The Poincaré map has been constructed to investigate the existence of periodic solutions in neuron model, which allow us to better understand the working mechanism of the brain.The conditions for the regular spiking, bursting, and phasic spiking are provided; as regular spiking corresponds to an order 1 periodic solution and bursting corresponds to an order  ( ≥ 2) periodic solution, this indicates that the parameters should be chosen carefully in order to control the number of spikes per burst.When an order 1 periodic solution loses its stability, model (2) exhibits very intricate dynamics, such as periodic windows, chaotic bands, period-doubling bifurcations, period-adding bifurcations, and period-halving bifurcations (for details see Figure 5).Particularly, it is noted that chaos whose electrophysiological signature is bursting has been observed in the simple neuron model as well as real neurons in vitro [28][29][30][31][32][33].The results imply that the analytical methods can be used as the basis for understanding dynamical behaviors of neuron models.
To provide a more comprehensive theoretical analysis for system (2), the dynamics presented in the remaining cases will be studied in the future.Recently, with the aim of improving the effective simulation studies, an alternative has been made to develop minimal models of spiking neurons which reduced the dimensionality of variable space [34].To this end, piecewise linear planar neuron models, which are described by a nonsmooth dynamical system, are a good choice and have become one of the most popular topics [34][35][36][37].Therefore, these will be the main focus of our future research.It is hoped that such research, planned for the near future and to be reported elsewhere, will be useful for computational neuroscience.

Figure 1 :
Figure 1: Illustrations of the domains of phase set and impulsive set for case ( 1 ).The blue lines represent the phase sets and the yellow line denotes the impulsive set.

Figure 2 :
Figure 2: Illustrations of the domains of phase set and impulsive set: (a) for case ( 2 ); (b) for case ( 3 ).