Neurocomputing 81 (2012) Contents lists available at SciVerse ScienceDirect. Neurocomputing. journal homepage:
|
|
- Cordelia Austin
- 6 years ago
- Views:
Transcription
1 Neurocomputing 81 (212) Contents lists available at SciVerse ScienceDirect Neurocomputing journal homepage: Bifurcations in the Hodgkin Huxley model exposed to DC electric fields Yanqiu Che a,n, Jiang Wang b, Bin Deng b, Xile Wei b, Chunxiao Han a a Tianjin Key Laboratory of Information Sensing & Intelligent Control, School of Automation and Electrical Engineering, Tianjin University of Technology and Education, Tianjin 3222, China b School of Electrical Engineering and Automation, Tianjin University, Tianjin 372, China article info Article history: Received 1 March 211 Received in revised form 23 September 211 Accepted 1 November 211 Communicated by J. Torres Available online 27 December 211 Keywords: Hodgkin Huxley model Hopf bifurcation Electric fields abstract Diverse behaviors of the original Hodgkin Huxley (HH) model, depending on the parameter values, have been studied extensively. This paper proposes modified HH equations exposed to externally applied extremely low frequency (ELF) electric fields. We investigate the effect of the DC electric fields on the dynamics of the modified HH model using bifurcation analysis. The obtained bifurcation sets partition the two dimensional parameter space, representing intensity of externally applied DC current and trans-membrane voltage induced by external DC electric fields, in terms of the qualitatively different behaviors of the HH model. Thus the neuronal information encodes the stimulus information, and vice versa. We also illustrate that the multi-stability phenomena in the HH model are associated with Hopf and double cycle bifurcations. & 211 Elsevier B.V. All rights reserved. 1. Introduction A neuron s response to external electrical stimulation varies depending on the neuron s characteristics as well as the stimulation. With a constant current, a single spike followed by resting potential (i.e. excitability) and periodic tonic spiking are two common types of responses. For the dynamical system [1], considering the steady states for these two behaviors, they correspond to a stable equilibrium and a stable limit cycle, respectively. The Hodgkin Huxley (HH) model [2] which quantitatively describes electrical excitations of squid giant axons is a prototype for excitable cell models. Under the HH formalism (the method used in deriving the HH model), a large family of mathematical models (HH-type models) for diverse neurons and excitable cells are established [3 7]. Since these HH-type models share many common dynamical characteristics with HH model, a detailed analysis of the dynamics in HH model under electric fields will also provide a basic understanding of that in HH-type models under electric fields. A bifurcation is a qualitative change in the behavior of a nonlinear dynamical system as its parameters pass through critical values [1]. The study of bifurcations in neural models is a keystone for understanding the dynamical origin of many single-neuron and circuit phenomena involved in neural information processing and the organization of behavior [8]. So far, many studies have been done on the bifurcation analysis of the original n Corresponding author. Tel./fax: þ address: yqche@tju.edu.cn (Y. Che). HH model [9 21]. With a change in the externally applied DC current I ext, Hopf bifurcation would occur at the equilibrium of HH model, generating limit cycles [9 11]. The corresponding stable and unstable periodic solutions are numerically calculated [12]. With appropriate I ext and V K, representing the equilibrium potential of potassium, two stable equilibrium potentials coexist in the HH model [13], which coincides with the experimental result [14]. Guckenheimer and Labouriau [15] give the detailed bifurcation diagrams of HH model in two-parameter space of I ext and V K. Bedrov et al. reveal the possible bifurcations with changes of g Na and g K, representing the maximal conductance of sodium and potassium, respectively [16,17]. The global structure of bifurcations in multiple-parameter space of the HH model is examined [18], and the details of the degenerate Hopf bifurcations are analyzed using the singularity theoretic approach [19]. It has also been shown that the chaotic solutions can exist in the HH equations by varying the time constant [2] or the original parameters [21]. On the other hand, the interaction between electric fields and biological tissues has been of longstanding interest. The electric fields can alter neuronal excitability in normal physiological conditions [22 27], change the neuronal ensemble activity of epileptiform [28 32]. The electric field stimulation becomes a powerful physiotherapy [33]. In spite of the importance of HH model, the detailed bifurcation analysis of the HH model exposed to external electric fields has not been carried out to the best of our knowledge. The aim of this paper is to illustrate how the DC electric fields affect the dynamics of neural activity in the context of HH equations. Based on the results we obtained in the literature /$ - see front matter & 211 Elsevier B.V. All rights reserved. doi:1.116/j.neucom
2 42 Y. Che et al. / Neurocomputing 81 (212) [34], the modified HH model exposed to extremely low frequency (ELF) external electric field is given. In the modified HH model a new parameter, induced trans-membrane potential, is introduced to denote the effect of external electric fields. We are concerned with the bifurcations of the modified HH model in two-parameter space of the externally applied DC current and the induced transmembrane voltage, while the other parameters are kept as the original values in HH. Computational tools based upon the theory of nonlinear dynamical systems are used. The researches will throw some light on the interference mechanism between biological systems and electric field exposure. 2. Neuron models 2.1. The original HH model The original HH neuron model consists of the following ordinary differential equations: dv C m dt ¼ I ext ½g Na m 3 hðv V Na Þþg K n 4 ðv V K Þþg L ðv V L ÞŠ dm dt ¼½a mðvþð1 mþ b m ðvþmšf dh dt ¼½a hðvþð1 hþ b h ðvþhšf dn dt ¼½a nðvþð1 nþ b n ðvþnšf ð1þ where V represents membrane potential in mv which is produced by the accumulation and transportation of ions such as sodium, potassium and chlorine across the membrane. The time t is measured in ms, the variables m, h and n represent the sodium activation, the sodium inactivation and the potassium activation, respectively. Parameters V Na, V K and V L are the equilibrium potentials of sodium, potassium and leak current, respectively. They are determined uniquely by Nernst s equations. g Na, g K and g L are the maximal conductance of the corresponding ionic currents. They reflect the ionic channel density distributed over the membrane. C m is the membrane capacitance, I ext is the externally applied current. The factor f ¼ 3 ðt 6:3Þ=1 alters the reaction rate constants a y ðvþ and b y ðvþðy ¼ m,h,nþ. a y ðvþ and b y ðvþðy ¼ m,h,nþ are nonlinear functions of V, given by the following equations: a m ðvþ¼:1 ð25 VÞ=½expðð25 VÞ=1Þ 1Š b m ðvþ¼4 expð V=18Þ a h ðvþ¼:7 expð V=2Þ b h ðvþ¼1:=½expðð V þ3þ=1þþ1š a n ðvþ¼:1 ð1 VÞ=½expðð1 VÞ=1Þ 1Š b n ðvþ¼:125 expð V=8Þ 2.2. The modified HH model exposed to external electric field For biological cells under external ELF electric fields, the capacitance effect of both cytoplasm and outside medium can be neglected because of their large reactance, while the capacitance effect of membrane must be taken into account due to its high resistance. Thus at cellular level, an exposure of a cell to an external electric field results in the induced membrane potential V E that superimposes to the ionic membrane potential V [34 36], ð2þ i.e. the overall membrane potential becomes V þv E. Then the original HH model is modified as follows: dðv þv E Þ C m ¼ I ext ½g dt Na m 3 hðv þv E V Na Þþg K n 4 ðv þv E V K Þ þg L ðv þv E V L ÞŠ dy dt ¼ a yðvþð1 yþ b y ðvþy ðy ¼ m,h,nþ ð3þ The induced membrane potential V E can be calculated as [34,35] V E ¼ KEðtÞR cos y½1 expð t=tþš ð4þ where E is the strength of ELF electric fields, K is a constant determined by the geometric and electrical characteristics of the cell, t is the time constant of membrane, R is the radius of the cell, and y is the angle between V E ðtþ and E(t). Since the effect of AC electric fields on the dynamics of HH model has been studied elsewhere [37], we only study the effect of DC electric fields in this paper, then in the left of Eq. (3), the term C m ðdv E =dtþ¼. The ELF electric fields directly act over membrane potential through an additive perturbation V E, which in turn influence the activity of voltage-dependent ionic channels. V E does not change the basic structure of HH model, but simultaneously shifts the three ionic reverse potentials, which is different from the situation studied in [38]. According to control theory, the introduction of V E could be regarded as a disturbance applied to the original system and its dynamic performance under disturbance should be investigated. Fixed parameter values used in this paper are listed in Table 1, with which the neuron keeps silent without externally applied electric currents and electric fields. 3. Definition and detection of bifurcations A bifurcation is a change of qualitative behavior of a dynamical system at special values of the parameters. The codimension of a bifurcation is the number of parameters which must be varied for the bifurcation to occur [15]. In this study, we consider the following bifurcations [15,18] Codimension one bifurcations Hopf bifurcation (H): A continuous fundamental path of an equilibrium point loses its stability as it intersects a secondary path of a periodic solution. The location of a Hopf bifurcation on the equilibrium point is characterized by a complex conjugate pair of linear eigenvalues of the Jacobian matrix whose real part passes through zero. When the secondary path is stable, it is the supercritical Hopf bifurcation (sh). Conversely, when the secondary path is unstable, it is the subcritical Hopf bifurcation (uh). Table 1 Fixed parameters for the HH model. Parameters Values Units Variables C m 1. mf=cm 2 Membrane capacitance V Na 115 mv Sodium equilibrium potential V K 12 mv Potassium equilibrium potential V L mv Leak equilibrium potential g Na 12 ms/cm 2 Maximal sodium conductance g K 36 ms/cm 2 Maximal potassium conductance g L.3 ms/cm 2 Maximal leak conductance T 6.3 1C Temperature
3 Y. Che et al. / Neurocomputing 81 (212) Saddle-node bifurcation (sn): Two equilibrium points coalesce and disappear. At this bifurcation point, the Jacobian matrix of the equations at the equilibrium point has a zero eigenvalue. Double cycle bifurcation (dc, saddle-node of cycles): Two periodic solutions with finite amplitude coalesce and disappear Codimension two bifurcations Cusp (c): Three equilibrium points coalesce into one. Cusp of periodic (cp): Three periodic solutions coalesce into one. On an appropriate Poincare section, this bifurcation corresponds to the cusp of the equilibrium points of the Poincare map. Takens Bogdanov bifurcation (TB): The Jacobian matrix of the equations at the equilibrium point has two zero eigenvalues. On a two-dimensional bifurcation diagram (2BD), TB locates on the sn curve and the Hopf bifurcation curve is tangent to the sn curve at this point. Degenerate Hopf bifurcation (dh): The stability of the periodic solution which bifurcates at the Hopf bifurcation point changes. A dc curve is terminated at this point on a 2BD. In this study, numerical computations for the bifurcation diagrams are conducted using MATCONT, a new MATLAB continuation package for the study of parameterized ODE systems. Matcont can detect several bifurcation points automatically and can trace both stable and unstable branches of equilibria and periodic solutions [39]. Numerical integration of the system equations for obtaining the trajectories is conducted in MATLAB using ode Bifurcation analysis In this section, we first present a global bifurcation diagram for the modified HH model (Eq. (3)) in I ext,v E two-dimensional parameter space. Then, we illustrate how the parameter V E affects the bifurcation process by analyzing a series of one parameter bifurcation diagrams for I ext with different fixed values of V E Bifurcations in I ext,v E parameter space Fig. 1 shows the two-parameter bifurcation diagram of the modified HH model. The abscissa and the ordinate are I ext and V E, 5 5 A C Fig. 1. Bifurcation diagram of HH model in (I ext, V E ) parameter space. The abscissa is constant DC externally applied current I ext and the ordinate is the electric fields parameter V E. The supercritical Hopf, subcritical Hopf, double-cycle, and saddle-node bifurcation curves are labeled as sh, uh, dc and sn, respectively. The parameter planes are divided into regions A I (regions F I are shown in Fig. 2) bythesebifurcation curves. In each region, a schematic phase portrait is shown. Solid dots and crosses represent stable and unstable equilibrium points, respectively. The curved arrows represent schematic trajectories of the equations. Stable and unstable periodic solutions are represented by solid and dashed circles, respectively. B A respectively. Fig. 2 gives detailed structure for some regions in Fig. 1. Three types of curves of codimension one bifurcations partition the whole parameter plane into several regions (A I). The solid one at the center is the Hopf bifurcation curve (uh and sh). Two dashed curves, lying in upper right, are double cycle bifurcation curves (dc). And the dash-dotted curve in the downleft corner is saddle-node bifurcation curve (sn). Hopf bifurcation changes stability at dh points and terminates at TB points. Two branches of sn join at a cusp (c) point, creating a region of coexistence of three equilibria. Two cp points are joint points of dc curves. Schematic phase portraits for each region are given. Stable equilibrium points are shown as solid dots, unstable ones are crosses; stable limit cycles are closed curves with solid lines, and unstable periodic orbits are dashed lines. Different phase portraits indicate qualitatively different behaviors of HH model. Fig. 3 shows some examples of waveforms of the membrane potentials in HH model for parameters located in the different regions of bistability (Fig. 3a e, h) or tri-stability (Fig. 3f, g). Properly added perturbations (arrows in figures) to the membrane potentials at certain time induce transitions from one stable state to another. In region A, the unique steady state is an stable equilibrium point. Parameter set ði ext,v E Þ¼ð; Þ locates in this region. The membrane is excitable, i.e. the membrane potential returns to the steady state after an excitation in respond to pulse stimulation. In region B, the periodic solution is the unique stable steady state and an unstable equilibrium point exist within it. The asymptotic dynamics of the modified HH is the periodic oscillation. In C, three equilibrium points coexist and two of them are stable, i.e. bistability of equilibrium points. Asymptotic dynamics is either one of them, depending on its initial conditions. Fig. 3(a) and (b) gives two examples of membrane potential waveforms when the parameter values are set in C ði ext,v E Þ¼ð 5, 5Þ and ði ext,v E Þ¼ð 1:755, 15:87Þ, respectively. Perturbations to the membrane potential induce transitions between the two steady states. In D, two stable attractors, an equilibrium point and a periodic solution, coexist with an unstable periodic solution in-between, i.e. bistability of an equilibrium point and a limit cycle. In E (Fig. 1), there is another bistability of an equilibrium point and a limit cycle but it is qualitatively different from that in region D. As illustrated by examples in Fig. 3(c) for region D, ði ext,v E Þ¼ð5; 25Þ and in Fig. 3(d) for region E, ði ext,v E Þ¼ð 3, 26Þ, respectively, the limit cycle encircles the equilibrium in D while it does not in E. In F (Fig. 2(a)), bistability of limit cycles is observed. In Fig. 3(e) for region F, ði ext,v E Þ¼ð115,22:9Þ, perturbations induce transitions between the small amplitude oscillation and the larger one. In G (Fig. 2(b)), four limit cycles with different periods and amplitudes coexist with one equilibrium point. Two of the limit cycles and the equilibrium point are stable (tri-stability) and the others are unstable. In Fig. 3(f) for region G, ði ext,v E Þ¼ð98:735,23:1993Þ, properly delivered perturbations transform membrane potential from resting state to with very small amplitude to oscillation with large amplitude. Note that the attractor regions of the equilibrium and the small limit cycle are very small. Region H(Fig. 2(c)) is another tri-stability region. It s boundary consists of dc 1, sn and uh. One stable and one unstable equilibrium points locate outside the stable limit cycle which encloses another stable equilibrium point and an unstable limit cycle. Fig. 3(g) for region H, ði ext,v E Þ¼ð 1:7497, 15:8776Þ, illustrates transitions between the three stable states. The narrow zone enclosed by curves of dc 6 and dc 7 between dc 1 and uh is denoted as Region I (Fig. 2(d)). In this region, between a stable equilibrium point and a stable limit cycle, there exist three even more unstable limit cycles [12,2]. Fig. 3(h) exemplifies a membrane potential waveform when the parameter values are set in I ði ext,v E Þ¼ð 1:, 5:43Þ.
4 44 Y. Che et al. / Neurocomputing 81 (212) Fig. 2. Detail bifurcation structures of regions F I shown in Fig. 1. Perturbations to the membrane potential induce transitions between the two steady states One parameter bifurcations for I ext To illustrate how V E affects the dynamics of HH model and how the multi-stabilities emerge, we investigate the changes of the one parameter bifurcation diagrams for I ext with different values of V E. The additive term V E is fixed to certain values and the bifurcation parameter I ext is continuously changed. Fig. 4(a i) shows the results for V E ¼, 2, 23, , , 23.25, 23.3, and 3 mv, respectively. In each diagram, extrema of membrane potential V extrema of the steady solutions of HH are plotted as functions of I ext. The central curve represents the equilibrium points and the upper and lower ones represent the maximum and minimum values of periodic solutions, respectively. Stable and unstable solutions are shown by solid and dotted curves, respectively. The I ext axis is divided into sub-intervals which are associated with regions in Figs. 1 and 2. Without or with very weak electric fields, e.g. V E ¼ mv (Fig. 4(a)), as I ext increases from small value, first reaches a dc bifurcation dc 1,andapairofstableandunstableperiodicsolutions emerge. Since dc does not affect the existing stable equilibrium point, the bistability of limit cycle and equilibrium point occurs. As I ext increases further, in a small interval of I ext, three branches of unstable limit cycles show azshapecurvejoinedbydc 6 and dc 7 bifurcation. When I ext reaches the left Hopf point uh, the stability of the equilibrium changes from stable to unstable, the bistability disappears. With further increase of I ext, amplitude of the stable limit cycle decreases and vanishes at the right Hopf point sh point, and the equilibrium becomes stable again. The neuron becomes to be quiescent. For a relative larger V E (e.g. V E ¼ 2 mv in Fig. 4(b)), the only different feature of the bifurcation diagram from that in Fig. 4(a) is the disappearance of the Z shape curve, leaving a single unstable periodic solution branch, due to collision and disappearance of dc 6 and dc 7.ForlargerV E, the bifurcation diagrams change dramatically in a small region, due to the coexistence of several codimension two bifurcations (see Fig. 2(b)). As shown in Fig. 4(c) for V E ¼ 23 mv, the right half of the stable limit cycle branch changes to be a Z shape, caused by the occurrence of bifurcations dc 2 and dc 3 between which two stable and one unstable limit cycles coexist with an unstable equilibrium point. When I ext is changed gradually passing this bistability region, the modified HH shows a hysteresis phenomenon, i.e. its steady state jumps from one type of periodic firing to the other at points dc 6 and dc 7. Hopf points uh and sh become closer. For V E ¼ 23:199 mv and V E ¼ 23:1993 mv (Fig. 4(d, e)), the left half of the unstable limit cycle branch creates a Z shape in a narrow region close to uh as shown in the enlargements. dc 4 and dc 5 bifurcations account for that changes. The middle limit cycle is stable rather than unstable as that occurs in Fig. 4(a), thus in the narrow parameter range, the modified HH exhibits tri-stability. The difference between (d) and (e) is the relative position of dc 4 and uh. For V E ¼ 23:25 mv (Fig. 4, dc 5 meets uh point, dc 5 disappears and uh changes into sh. The periodic solution bifurcates form this sh point is stable rather than unstable and the bifurcation direction changes. The dc 3 and dc 4 get closer and nearly collide. For V E ¼ 23:3mV (Fig. 4, afterdc 3 and dc 4 coalesce and disappear, an isolate periodic solution branch which is not associated with Hopf point (named isola [1]) exists in a wide range of I ext. The periodic solutions bifurcate from the two Hopf points join at certain value of I ext.forv E ¼ 23:455 mv (Fig. 4, the two sh points both occur before dc 2 and get much closer. For large V E (in Fig. 4, V E ¼ 3 mv), the two Hopf points coalesce and disappear after a dh bifurcation dh 2, leaving an isola and a stable equilibrium point. In a summary, as V E increases, the dc bifurcations cause the occurrence and disappearance of Z-shaped branch of the periodic
5 Y. Che et al. / Neurocomputing 81 (212) C 5 C D 1 E F 8 G H 15 I Fig. 3. Examples of membrane potential waveforms in bistability or tri-stability regions C I in Figs. 1 and 2. The abscissa and the ordinate are time (ms) and the membrane potential (mv), respectively. In each case, all the parameters of the HH are fixed, and appropriate impulsive perturbations delivered to the membrane potential V lead to transitions from one steady state to another (see text). solution, while the Hopf bifurcation points uh and sh get closer, and the left uh becomes sh. Then, the periodic solution branch is disconnected with equilibrium, leading to the appearance of the isola. As V E increases further, the left and right sh collide and disappear by dh 2. The isola remains for higher values of V E.SothedcandHopf bifurcations mainly contribute to the emergence of multi-stabilities.
6 46 Y. Che et al. / Neurocomputing 81 (212) Fig. 4. One parameter bifurcation diagrams for I ext with different values of V E. The abscissa and the ordinate are I ext and the extrema of membrane potential V extrema, respectively. In each diagram, the central curve represents the equilibrium potential. The upper and lower curves represent the maximum and minimum values, respectively, of the membrane potential of periodic solutions. Solid and dashed curves represent stable and unstable solutions, respectively. sh, uh and dc are supercritical, subcritical Hopf bifurcation, and double-cycle bifurcation, respectively.
7 Y. Che et al. / Neurocomputing 81 (212) Conclusion Many studies have been done on the diverse behaviors and bifurcations in the original HH model. This paper have investigated the effect of DC electric fields on the bifurcation structure of the HH model. We summarize the results on a two-dimensional I ext, V E bifurcation diagram. The loci of codimension one bifurcations, Hopf, saddle-node and double cycle considered in this paper, partition the parameter space into several qualitatively different regions. The stimulus parameters can be quantified by the various neuronal firing patterns. On the one hand, one can interpret the evoked spike sequences with various patterns as a representation of the parameter values. On the other hand, one can determine the parameter values to make the neuron elicit certain specific spike sequences, which is significant in neural control. By analyzing the systematic changes of bifurcation diagrams for I ext with different values of V E,weillustrate that the Hopf and double cycle bifurcations account for the bistabilities and multi-stabilities phenomena. Acknowledgement This work is supported by the Key Program of the National Natural Science Foundation of China (Grant No ), the Young Scientists Fund of the National Natural Science Foundation of China (Grants No and No ), the National Natural Science Foundation of China (Grant No ) and the National Science Foundation for Post-doctoral Scientists of China (Grant No ). References [1] J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcation of Vector Fields, Springer, [2] A.L. Hodgkin, A.F. Huxley, A quantitative description of membrane and its application to conduction and excitation in nerve, J. Physiol. 117 (1952) [3] J. Cronin, Mathematical Aspects of Hodgkin Huxley Neural Theory, Cambridge University Press, [4] T.R. Chay, J. Keizer, Minimal model for membrane oscillations in the pancreatic b-cell, Biophys. J. 42 (1983) [5] R.D. Traub, R.K.S. Wong, R. Miles, H. Michelson, A model of a CA3 hippocampal pyramidal neuron incorporating voltage-clamp data on intrinsic conductances, J. Neurophysiol. 66 (1991) [6] D. Noble, The development of mathematical models of the heart, Chaos Solitons Fractals 5 (1995) [7] J. Keener, J. Sneyd, Mathematical Physiology. Interdisciplinary Applied Mathematics, vol. 8, Springer-Verlag, New York, [8] M.I. Rabinovich, P. Varona, A.I. Selverston, H.D.I. Abarbanel, Dynamical principles in neuroscience, Rev. Mod. Phys. 78 (26) [9] B. Hassard, Bifurcation of periodic solutions of the Hodgkin Huxley model for the squid giant axon, J. Theory Biol. 71 (1978) [1] J. Rinzel, One repetitive activity in nerve, Fed. Proc. 37 (1978) [11] W.C. Troy, Bifurcations of periodic-solutions in Hodgkin Huxley equations, Q. Appl. Math. 36 (1978) [12] J. Rinzel, R.N. Miller, Numerical calculation of stable and unstable solutions to the Hodgkin Huxley equations, Math. Biosci. 49 (198) [13] K. Aihara, G. Matsumoto, Two stable steady states in the Hodgkin Huxley axons, J. Biophys. 41 (1983) [14] I. Tasaki, Demonstration of two stable states of the nerve membrane in potassium-rich media, J. Physiol. (London) 148 (1959) [15] J. Guckenheimer, I.S. Labouriau, Bifurcation of the Hodgkin and Huxley equations: a new twist, Bull. Math. Biol. 55 (1993) [16] Y.A. Bedrov, G.N. Akoev, O.E. Dick, Partition of the Hodgkin Huxley type model parameter space in to the region of qualitatively different solutions, Biol. Cybern. 66 (1992) [17] Y.A. Bedrov, G.N. Akoev, O.E. Dick, On the relationship between the number of negative slope regions in the voltage current curve of the Hodgkin-Huxley model and its parameter values, Biol. Cybern. 73 (1995) [18] H. Fukai, S. Doi, T. Nomura, S. Sato, Hopf bifurcations in multiple parameter space of the Hodgkin Huxley equations. I. Global organization of bistable periodic solutions, Biol. Cybern. 82 (2) [19] H. Fukai, S. Doi, T. Nomura, S. Sato, Hopf bifurcations in multiple parameter space of the Hodgkin Huxley equations. II. Singularity theoretic approach and highly degenerate bifurcations, Biol. Cybern. 82 (2) [2] S. Doi, S. Nabetani, S. Kumagai, Complex nonlinear dynamics of the Hodgkin Huxley equations induced by time scale changes, Biol. Cybern. 85 (21) [21] J. Guckenheimer, A. Oliva, Chaos in the Hodgkin Huxley model, SIAM J. Appl. Dyn. Syst. 1 (22) [22] C.Y. Chan, C. Nicholson, Modulation by applied electric fields of Purkinje and stellate cell activity in the isolated turtle cerebellum, J. Physiol. 371 (1986) [23] C.Y. Chan, J. Hounsgaard, C. Nicholson, Effects of electric fields on transmembrane potential and excitability of turtle cerebellar Purkinje cells in vitro, J. Physiol. 42 (1988) [24] J.G.R. Jefferys, Influence of electric fields on the excitability of granule cells in guinea-pig hippocampal slices, J. Physiol. 319 (1981) [25] S.M. Bawin, A.R. Sheppard, M.D. Mahoney, W.R. Adey, Influences of sinusoidal electric fields on excitability in the rat hippocampal slice, Brain Res. 323 (1984) [26] J. Berzhanskaya, A. Gorchetchnikov, S.J. Schiff, Switching between gamma and theta: dynamic network control using subthreshold electric fields, Neurocomputing 7 (27) [27] M. Bikson, M. Inoue, H. Akiyama, J.K. Deans, J.E. Fox, H. Miyakawa, G.R. Jefferys, Effects of uniform extracellular DC electric fields on excitability in rat hippocampal slices in vitro, J. Physiol. 557 (24) [28] B.J. Gluckman, E.J. Neel, T.I. Netoff, W.L. Ditto, M.L. Spano, S.J. Schiff, Electric field suppression of epileptiform activity in hippocampal slices, J. Neurophysiol. 76 (1996) [29] D.H. Duong, T. Chang, The influence of electric fields on the epileptiform bursts induced by high potassium in CA3 region of rat hippocampal slices, Neurol. Res. 2 (1998) [3] R. Ghai, M. Bikson, D.M. Durand, Effects of applied electric field on lowcalcium epileptiform activity in the CAI region of rat hippocampal slices, J. Neurophysiol. 84 (2) [31] B.J. Gluckman, H. Nguyen, S.L. Weinstein, S.J. Schiff, Adaptive electric field suppression of epileptic seizures, J. Neurosci. 21 (21) [32] E.-H. Park, P. So, E. Barreto, B.J. Gluckman, S.J. Schiff, Electric field modulation of synchronization in neuronal networks, Neurocomputing (22) [33] J.P. Reilly, Electrical Stimulation and Electropathology, Cambridge University Press, Cambridge, [34] J. Wang, K.M. Tsang, H. Zhang, Hopf bifurcation in the Hodgkin Huxley model exposed to ELF electrical field, Chaos Solitons Fractals 2 (24) [35] T. Kotnik, G. Pucihar, D. Miklavc ic, The cell in the electric field, in: T.K. Stephen, G. Julie, W.L. Edward (Eds.), Clinical Aspects of Electroporation, 1st ed.,springer, 211, pp [36] C. Polk, E. Postow (Eds.), CRC Handbook of Biological Effects of Electromagnetic Fields, second ed., CRC Press, Boca Raton, [37] Y.Q. Che, J. Wang, W.J. Si, X.Y. Fei, Phase-locking and chaos in a silent Hodgkin Huxley neuron exposed to sinusoidal electric field, Chaos Solitons Fractals 39 (29) [38] Y.A. Bedrov, G.N. Akoev, O.E. Dick, Partition of the Hodgkin Huxley type model parameter space into the regions of qualitatively different solutions, Biol. Cybern. 66 (5) (1991) [39] A. Dhooge, W. Govaerts, Y.A. Kuznetsov, MATCONT: a MATLAB package for numerical bifurcation analysis of ODEs, ACM Trans. Math. Software 29 (23) Yanqiu Che received BS, M.Eng and the Ph.D. degrees in Control Theory and Engineering from Tianjin University, Tianjin, China in 23, 25 and 28, respectively. During28 21hewasapostdocfellowatTianjin University. Since 21, He is an associate professor at the School of Automation and Electrical Engineering, Tianjin University of Technology and Education. His research interests focus on Neurodynamics, Computational Neuroscience, Neural Control Engineering, Bioelectromagnetics and Brain Computer Interface. Jiang Wang received the Ph.D. degree in Control Theory and Engineering from Tianjin University in He is currently a professor at the School of Electrical Engineering and Automation, Tianjin University. His research interests are nonlinear dynamical systems, computational neuroscience, and information processing and detecting.
8 48 Y. Che et al. / Neurocomputing 81 (212) Bin Deng received BS, M.Eng and the Ph.D. degrees in Control Theory and Engineering from Tianjin University, Tianjin, China in 21, 24 and 27, respectively. During he was a postdoc fellow at Tianjin University. Since 29, He is an associate professor of School of Electrical Engineering and Automation, Tianjin University. His research interests are analysis and control of nonlinear system, the dynamics of neural system, and the analysis of neural signals. Chunxiao Han received the BS and Ph.D. degrees in Control Theory and Engineering from Tianjin University, Tianjin, China in 23 and 28, respectively. She is currently a lecturer at the School of Automation and Electrical Engineering, Tianjin University of Technology and Education. Her research interests focus on biosignal processing and neurodynamics. Xile Wei received BS, M.Eng and the Ph.D. degrees in Control Theory and Engineering from Tianjin University, Tianjin, China in 1999, 24 and 27, respectively. During he was a postdoc fellow at Tianjin University. Since 29, He is an associate professor of School of Electrical Engineering and Automation, Tianjin University. His current research interests are in the areas of Bioelectromagnetics, Neuromorphic Modeling, and Neural Control Engineering.
Dynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More informationChaos in the Hodgkin Huxley Model
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 1, No. 1, pp. 105 114 c 2002 Society for Industrial and Applied Mathematics Chaos in the Hodgkin Huxley Model John Guckenheimer and Ricardo A. Oliva Abstract. The
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationBifurcation Analysis of Regulation of Bursting Discharge in DRG Neurons
Chaotic Modeling and Simulation (CMSIM) 2: 181-191, 217 Bifurcation Analysis of Regulation of Bursting Discharge in DRG Neurons Olga E. Dick Pavlov Institute of Physiology of Russian Academy of Science,
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More information2182. Hopf bifurcation analysis and control of three-dimensional Prescott neuron model
2182. Hopf bifurcation analysis and control of three-dimensional Prescott neuron model Chunhua Yuan 1, Jiang Wang 2 School of Electrical Engineering and Automation, Tianjin University, Tianjin, China 2
More informationDynamical Systems in Neuroscience: The Geometry of Excitability and Bursting
Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting Eugene M. Izhikevich The MIT Press Cambridge, Massachusetts London, England Contents Preface xv 1 Introduction 1 1.1 Neurons
More informationChapter 2 The Hodgkin Huxley Theory of Neuronal Excitation
Chapter 2 The Hodgkin Huxley Theory of Neuronal Excitation Hodgkin and Huxley (1952) proposed the famous Hodgkin Huxley (hereinafter referred to as HH) equations which quantitatively describe the generation
More informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationComputational Neuroscience. Session 4-2
Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.
More informationAll-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model
All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model S. A. Sadegh Zadeh, C. Kambhampati International Science Index, Mathematical and Computational Sciences waset.org/publication/10008281
More informationA Mathematical Study of Electrical Bursting in Pituitary Cells
A Mathematical Study of Electrical Bursting in Pituitary Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Collaborators on
More informationLecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers
Lecture Notes 8C2 Inleiding Meten en Modelleren Cellular electrophysiology: modeling and simulation Nico Kuijpers nico.kuijpers@bf.unimaas.nl February 9, 2 2 8C2 Inleiding Meten en Modelleren Extracellular
More informationNumerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a Mathematical Model of Snail Neurons
International Journal of Theoretical and Mathematical Physics 2015, 5(5): 145-150 DOI: 10.5923/j.ijtmp.20150505.08 Numerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a
More informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best Koch, Chapters 6, 8, 9 See also Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913
More informationStrange Nonchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley Neuron
Journal of the Korean Physical Society, Vol. 57, o. 1, July 2010, pp. 23 29 Strange onchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley euron Woochang Lim and Sang-Yoon Kim Department of Physics,
More informationFrom neuronal oscillations to complexity
1/39 The Fourth International Workshop on Advanced Computation for Engineering Applications (ACEA 2008) MACIS 2 Al-Balqa Applied University, Salt, Jordan Corson Nathalie, Aziz Alaoui M.A. University of
More informationCHAPTER 1. Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting
CHAPTER 1 Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting John Guckenheimer Mathematics Department, Cornell University Ithaca, NY 14853 Joseph H. Tien Center for Applied Mathematics,
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More information1 Introduction and neurophysiology
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 (2009) 535-549 Copyright c 2009 Watam Press http://www.watam.org ASYMPTOTIC DYNAMICS OF THE SLOW-FAST HINDMARSH-ROSE
More informationSingle-Compartment Neural Models
Single-Compartment Neural Models BENG/BGGN 260 Neurodynamics University of California, San Diego Week 2 BENG/BGGN 260 Neurodynamics (UCSD) Single-Compartment Neural Models Week 2 1 / 18 Reading Materials
More informationQUARTERLY OF APPLIED MATHEMATICS
QUARTERLY OF APPLIED MATHEMATICS Volume LIV December 1996 Number 4 DECEMBER 1996, PAGES 601-607 ON EXISTENCE OF PERIODIC ORBITS FOR THE FITZHUGH NERVE SYSTEM By S. A. TRESKOV and E. P. VOLOKITIN Institute
More informationMathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued
Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland
More informationMathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells
Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University
More informationBo Deng University of Nebraska-Lincoln UNL Math Biology Seminar
Mathematical Model of Neuron Bo Deng University of Nebraska-Lincoln UNL Math Biology Seminar 09-10-2015 Review -- One Basic Circuit By Kirchhoff's Current Law 0 = I C + I R + I L I ext By Kirchhoff s Voltage
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationPhase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable
Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary Figure S1. Pulses >3mJ reduce membrane resistance in HEK cells. Reversal potentials in a representative cell for IR-induced currents with laser pulses of 0.74 to
More information6.3.4 Action potential
I ion C m C m dφ dt Figure 6.8: Electrical circuit model of the cell membrane. Normally, cells are net negative inside the cell which results in a non-zero resting membrane potential. The membrane potential
More informationReduction of Conductance Based Models with Slow Synapses to Neural Nets
Reduction of Conductance Based Models with Slow Synapses to Neural Nets Bard Ermentrout October 1, 28 Abstract The method of averaging and a detailed bifurcation calculation are used to reduce a system
More informationBursting Types and Bifurcation Analysis in the Pre-Bötzinger Complex Respiratory Rhythm Neuron
International Journal of Bifurcation and Chaos Vol. 27 No. 1 (2017) 1750010 (17 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127417500109 Bursting Types and Bifurcation Analysis in the
More informationPhase Response Curves, Delays and Synchronization in Matlab
Phase Response Curves, Delays and Synchronization in Matlab W. Govaerts and B. Sautois Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 281-S9, B-9000 Ghent, Belgium
More informationLecture 11 : Simple Neuron Models. Dr Eileen Nugent
Lecture 11 : Simple Neuron Models Dr Eileen Nugent Reading List Nelson, Biological Physics, Chapter 12 Phillips, PBoC, Chapter 17 Gerstner, Neuronal Dynamics: from single neurons to networks and models
More informationTopological target patterns and population oscillations in a network with random gap junctional coupling
0 0 0 0 Neurocomputing }0 (00) } Topological target patterns and population oscillations in a network with random gap junctional coupling Timothy J. Lewis*, John Rinzel Center for Neural Science and Courant
More informationStochastic Oscillator Death in Globally Coupled Neural Systems
Journal of the Korean Physical Society, Vol. 52, No. 6, June 2008, pp. 19131917 Stochastic Oscillator Death in Globally Coupled Neural Systems Woochang Lim and Sang-Yoon Kim y Department of Physics, Kangwon
More informationModeling Action Potentials in Cell Processes
Modeling Action Potentials in Cell Processes Chelsi Pinkett, Jackie Chism, Kenneth Anderson, Paul Klockenkemper, Christopher Smith, Quarail Hale Tennessee State University Action Potential Models Chelsi
More informationNonlinear Dynamics of Neural Firing
Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials
More informationBIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1
KYBERNETIKA VOLUME 3 0 (1994), NUMBER 2, PAG ES 121-128 BIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1 YURI A. KUZNETSOV AND CARLO PICCARDI An adaptive control system
More informationIntroduction and the Hodgkin-Huxley Model
1 Introduction and the Hodgkin-Huxley Model Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference:
More information1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model
1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model Alan Hodgkin and Andrew Huxley developed the first quantitative model of the propagation of an electrical signal (the action potential)
More informationEntrainment and Chaos in the Hodgkin-Huxley Oscillator
Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42 Overview (1) Goal: Show that the
More informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More informationTHREE SCENARIOS FOR CHANGING OF STABILITY IN THE DYNAMIC MODEL OF NERVE CONDUCTION
THREE SCENARIOS FOR CHANGING OF STABILITY IN THE DYNAMIC MODEL OF NERVE CONDUCTION E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper deals with the specific cases
More informationDifference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay
Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,
More informationBiomedical Instrumentation
ELEC ENG 4BD4: Biomedical Instrumentation Lecture 5 Bioelectricity 1. INTRODUCTION TO BIOELECTRICITY AND EXCITABLE CELLS Historical perspective: Bioelectricity first discovered by Luigi Galvani in 1780s
More informationSTUDENT PAPER. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters Education Program 736 S. Lombard Oak Park IL, 60304
STUDENT PAPER Differences between Stochastic and Deterministic Modeling in Real World Systems using the Action Potential of Nerves. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters
More informationAN ELECTRIC circuit containing a switch controlled by
878 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 46, NO. 7, JULY 1999 Bifurcation of Switched Nonlinear Dynamical Systems Takuji Kousaka, Member, IEEE, Tetsushi
More informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationEvaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons
International Journal of Theoretical and Mathematical Physics 2016, 6(3): 99-103 DOI: 10.5923/j.ijtmp.20160603.02 Evaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons Takaaki Shirahata
More informationBifurcation structure of a model of bursting pancreatic cells
BioSystems 63 (2001) 3 13 www.elsevier.com/locate/biosystems Bifurcation structure of a model of bursting pancreatic cells Erik Mosekilde a, *, Brian Lading a, Sergiy Yanchuk b, Yuri Maistrenko b a Department
More informationThe selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales
CHAOS 18, 015105 2008 The selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales Jonathan Rubin Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania
More informationBursting Oscillations of Neurons and Synchronization
Bursting Oscillations of Neurons and Synchronization Milan Stork Applied Electronics and Telecommunications, Faculty of Electrical Engineering/RICE University of West Bohemia, CZ Univerzitni 8, 3064 Plzen
More informationLIMIT CYCLE OSCILLATORS
MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin
More informationSystem Identification for the Hodgkin-Huxley Model using Artificial Neural Networks
Proceedings of International Joint Conference on Neural Networks, Orlando, Florida, USA, August 12-17, 2007 System Identification for the Hodgkin-Huxley Model using Artificial Neural Networks Manish Saggar,
More informationBursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model
Bursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model Abhishek Yadav *#, Anurag Kumar Swami *, Ajay Srivastava * * Department of Electrical Engineering, College of Technology,
More informationParameters for Minimal Model of Cardiac Cell from Two Different Methods: Voltage-Clamp and MSE Method
Parameters for Minimal Model of Cardiac Cell from Two Different Methods: oltage-clamp and MSE Method Soheila Esmaeili 1, * and Bahareh beheshti 1 Department of Biomedical engineering, ran University of
More informationMath 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction Junping Shi College of William and Mary November 8, 2018 Neuron Neurons Neurons are cells in the brain and other subsystems
More informationMixed mode oscillations as a mechanism for pseudo-plateau bursting
J Comput Neurosci (2010) 28:443 458 DOI 10.1007/s10827-010-0226-7 Mixed mode oscillations as a mechanism for pseudo-plateau bursting Theodore Vo Richard Bertram Joel Tabak Martin Wechselberger Received:
More informationDynamic Systems: Ordinary Differential Equations. Ordinary Differential Equations
Dynamic Systems: Ordinary Differential Equations Adapted From: Numerical Methods in Biomedical Engineering Stanley M. Dunn, Alkis Constantinides, Prabhas V. Moghe Chapter 7 Kim Ferlin and John Fisher Ordinary
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Hodgkin-Huxley Model for Nerve Cell Action Potential Part 1 Dr. Zvi Roth (FAU) 1 References Hoppensteadt-Peskin Ch. 3 for all the mathematics. Cooper s The Cell
More informationMANY scientists believe that pulse-coupled neural networks
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 3, MAY 1999 499 Class 1 Neural Excitability, Conventional Synapses, Weakly Connected Networks, and Mathematical Foundations of Pulse-Coupled Models Eugene
More informationarxiv: v2 [physics.optics] 15 Dec 2016
Bifurcation analysis of the Yamada model for a pulsing semiconductor laser with saturable absorber and delayed optical feedback Abstract Soizic Terrien 1, Bernd Krauskopf 1, Neil G.R. Broderick 2 arxiv:1610.06794v2
More informationFRTF01 L8 Electrophysiology
FRTF01 L8 Electrophysiology Lecture Electrophysiology in general Recap: Linear Time Invariant systems (LTI) Examples of 1 and 2-dimensional systems Stability analysis The need for non-linear descriptions
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationNeuroscience applications: isochrons and isostables. Alexandre Mauroy (joint work with I. Mezic)
Neuroscience applications: isochrons and isostables Alexandre Mauroy (joint work with I. Mezic) Outline Isochrons and phase reduction of neurons Koopman operator and isochrons Isostables of excitable systems
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 2, W.W. Lytton, Hodgkin-Huxley Model) 9. Passive and active membrane models In the previous lecture we have considered a passive
More informationLinearization of F-I Curves by Adaptation
LETTER Communicated by Laurence Abbott Linearization of F-I Curves by Adaptation Bard Ermentrout Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A. We show that negative
More informationConductance-Based Integrate-and-Fire Models
NOTE Communicated by Michael Hines Conductance-Based Integrate-and-Fire Models Alain Destexhe Department of Physiology, Laval University School of Medicine, Québec, G1K 7P4, Canada A conductance-based
More informationEffects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell
Effects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell 1. Abstract Matthew Dunlevie Clement Lee Indrani Mikkilineni mdunlevi@ucsd.edu cll008@ucsd.edu imikkili@ucsd.edu Isolated
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 12, W.W. Lytton, Hodgkin-Huxley Model) 9.1 Passive and active membrane models In the previous lecture we have considered a passive
More informationDynamics of four coupled phase-only oscillators
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phase-only oscillators Richard Rand *, Jeffrey
More informationNonlinear Observer Design and Synchronization Analysis for Classical Models of Neural Oscillators
Nonlinear Observer Design and Synchronization Analysis for Classical Models of Neural Oscillators Ranjeetha Bharath and Jean-Jacques Slotine Massachusetts Institute of Technology ABSTRACT This work explores
More informationSimulation of Cardiac Action Potentials Background Information
Simulation of Cardiac Action Potentials Background Information Rob MacLeod and Quan Ni February 7, 2 Introduction The goal of assignments related to this document is to experiment with a numerical simulation
More informationDynamics from Seconds to Hours in Hodgkin Huxley Model
Dynamics from Seconds to Hours in Hodgkin Huxley Model with Time Dependent Ion Concentrations and Buffer Reservoirs Niklas Hübel Technische Universität Berlin July 9, 214 OUTLINE Closed Models model review
More informationNonlinear dynamics vs Neuronal Dynamics
Nonlinear dynamics vs Neuronal Dynamics Cours de Methodes Mathematiques en Neurosciences Brain and neuronal biological features Neuro-Computationnal properties of neurons Neuron Models and Dynamical Systems
More informationIntegral Equation and PDE Neuronal Models:
Integral Equation and PDE Neuronal Models: 1952-2015 W. C. Troy Integral Equation and PDE Neuronal Models: 1952-2015 p.1/51 Part 1. 1D Single Neuron Models. (i) Hodgkin-Huxley (1952), Fitzhugh-Nagumo (1961)
More informationControlling chaos in Colpitts oscillator
Chaos, Solitons and Fractals 33 (2007) 582 587 www.elsevier.com/locate/chaos Controlling chaos in Colpitts oscillator Guo Hui Li a, *, Shi Ping Zhou b, Kui Yang b a Department of Communication Engineering,
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationSingle-Cell and Mean Field Neural Models
1 Single-Cell and Mean Field Neural Models Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 The neuron
More informationLecture 10 : Neuronal Dynamics. Eileen Nugent
Lecture 10 : Neuronal Dynamics Eileen Nugent Origin of the Cells Resting Membrane Potential: Nernst Equation, Donnan Equilbrium Action Potentials in the Nervous System Equivalent Electrical Circuits and
More informationHopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators
The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of
More informationDynamics and complexity of Hindmarsh-Rose neuronal systems
Dynamics and complexity of Hindmarsh-Rose neuronal systems Nathalie Corson and M.A. Aziz-Alaoui Laboratoire de Mathématiques Appliquées du Havre, 5 rue Philippe Lebon, 766 Le Havre, FRANCE nathalie.corson@univ-lehavre.fr
More informationDetermining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model
JComputNeurosci(2014)37:403 415 DOI 10.1007/s10827-014-0511-y Determining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model Sevgi Şengül & Robert Clewley & Richard Bertram
More informationDynamical phase transitions in periodically driven model neurons
Dynamical phase transitions in periodically driven model neurons Jan R. Engelbrecht 1 and Renato Mirollo 2 1 Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA 2 Department
More informationSlow Manifold of a Neuronal Bursting Model
Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,
More informationSingle neuron models. L. Pezard Aix-Marseille University
Single neuron models L. Pezard Aix-Marseille University Biophysics Biological neuron Biophysics Ionic currents Passive properties Active properties Typology of models Compartmental models Differential
More informationGenerating hyperchaotic Lu attractor via state feedback control
Physica A 364 (06) 3 1 www.elsevier.com/locate/physa Generating hyperchaotic Lu attractor via state feedback control Aimin Chen a, Junan Lu a, Jinhu Lu b,, Simin Yu c a College of Mathematics and Statistics,
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Small-signal neural models and their applications Author(s) Basu, Arindam Citation Basu, A. (01). Small-signal
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationA model of slow plateau-like oscillations based upon the fast Na current in a window mode
Neurocomputing 38}40 (2001) 159}166 A model of slow plateau-like oscillations based upon the fast Na current in a window mode G.S. Cymbalyuk*, R.L. Calabrese Biology Department, Emory University, 1510
More informationChaos in the Hodgkin-Huxley Equations: The Takens-Bodganov Cusp Bifurcation. Rami Mohieddine
CORNELL UNIVERSITY MATHEMATICS DEPARTMENT SENIOR THESIS Chaos in the Hodgkin-Huxley Equations: The Takens-Bodganov Cusp Bifurcation A THESIS PRESENTED IN PARTIAL FULFILLMENT OF CRITERIA FOR HONORS IN MATHEMATICS
More informationTitle. Author(s)Fujii, Hiroshi; Tsuda, Ichiro. CitationNeurocomputing, 58-60: Issue Date Doc URL. Type.
Title Neocortical gap junction-coupled interneuron systems exhibiting transient synchrony Author(s)Fujii, Hiroshi; Tsuda, Ichiro CitationNeurocomputing, 58-60: 151-157 Issue Date 2004-06 Doc URL http://hdl.handle.net/2115/8488
More informationTowards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University
Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Dynamical systems with multiple time scales arise naturally in many domains. Models of neural systems
More informationMATH 3104: THE HODGKIN-HUXLEY EQUATIONS
MATH 3104: THE HODGKIN-HUXLEY EQUATIONS Parallel conductance model A/Prof Geoffrey Goodhill, Semester 1, 2009 So far we have modelled neuronal membranes by just one resistance (conductance) variable. We
More informationA note on discontinuous rate functions for the gate variables in mathematical models of cardiac cells
Procedia Computer Science (2) (22) 6 945 95 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science ICCS 2 A note on discontinuous rate functions for
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More informationMATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev
MATH 723 Mathematical Neuroscience Spring 28 Instructor: Georgi Medvedev 2 Lecture 2. Approximate systems. 2. Reduction of the HH model to a 2D system. The original HH system consists of 4 differential
More informationarxiv:q-bio/ v1 [q-bio.sc] 11 Jan 2006
arxiv:q-bio/0601013v1 [q-bio.sc] 11 Jan 2006 Effect of channel block on the spiking activity of excitable membranes in a stochastic Hodgkin-Huxley model G. Schmid, I. Goychuk and P. Hänggi Institut für
More information