Patterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops

Size: px
Start display at page:

Download "Patterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops"

Transcription

1 Math. Model. Nat. Phenom. Vol. 5, No. 2, 2010, pp DOI: /mmnp/ Patterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops J. Ma 1 and J. Wu 2 1 Department of Mathematics, University of Houston, Houston TX , USA 2 Center for Disease Modeling; Laboratory for Industrial and Applied Mathematics, Department of Mathematics and Statistics, York University, Toronto, Ontario, M3J 1P3, Canada Abstract. We study the coexistence of multiple periodic solutions for an analogue of the integrateand-fire neuron model of two-neuron recurrent inhibitory loops with delayed feedback, which incorporates the firing process and absolute refractory period. Upon receiving an excitatory signal from the excitatory neuron, the inhibitory neuron emits a spike with a pattern-related delay, in addition to the synaptic delay. We present a theoretical framework to view the inhibitory signal from the inhibitory neuron as a self-feedback of the excitatory neuron with this additional delay. Our analysis shows that the inhibitory feedbacks with firing and the absolute refractory period can generate four basic types of oscillations, and the complicated interaction among these basic oscillations leads to a large class of periodic patterns and the occurrence of multistability in the recurrent inhibitory loop. We also introduce the average time of convergence to a periodic pattern to determine which periodic patterns have the potential to be used for neural information transmission and cognition processing in the nervous system. Key words: multistability, periodic pattern, neural network, time delay, pattern formation, recurrent inhibitory loops AMS subject classification: 34K13, 34K18, 92B20 Corresponding author. wujh@mathstat.yorku.ca 67 Article published by EDP Sciences and available at or

2 1. Introduction In a living nervous system, recurrent loops involving two or more neurons are ubiquitous and are particularly prevalent in cortical regions for memory such as the hippocampal-mesial temporal lobe complex 26. The recurrent neural loops include the multiple pathways involved in the control of movement 2, reciprocal thalamocortical loops involved in epileptic seizures 11 and regulation of states of arousal 6, the limbic nervous system loops related to memory 1 and epileptic seizures 24, and the cortical-basal ganglia-thalamus-brain stem-cortical loops, which participate in the control of movement and act as gate keepers for the propagation of epileptic seizures 23. In this paper, we consider a simple recurrent inhibitory loop consisting of an excitatory neuron E and an inhibitory neuron I, where neuron E gives off collateral branches and excites the inhibitory neuron I, which in turns inhibits the firing of neuron E, in a delay time. Such twoneuron recurrent inhibitory loops with delay display similar complex dynamic behaviors as larger networks and many techniques developed to deal with two-neuron networks can carry over to networks of large size. Moreover, two-neuron networks are sometimes thought of as systems of two modules, where each module represents the mean activity of a spatially localized neural population 3, 25. The focus of this paper is on the capacity of the recurrent inhibitory loop to generate multiple coexisting periodic patterns (multistability). The coexistence of multiple stable patterns in neural networks is the basis of the mechanism for (associative) content-addressable memory storage and retrieval 8, 13, 14, 20, 22 where each stable equilibrium is identified with a static memory, while stable periodic orbits are associated with temporally patterned spike trains 4, 8, 9. Periodic patterns exhibited in neural networks have been linked to a variety of rhythms, which are associated with important behavioral and cognitive states in the nervous system, including attention, working memory, associative memory, object recognition, sensory motor integration and perception processing 5, 7, 15. Time delays, a powerful mechanism for multistability, are intrinsic properties of the nervous systems and are unavoidable in electronic implementation due to axonal conduction times, distances of interneurons and the finite switching speeds of amplifiers. Multistability in a delayed neural network has been extensively studied in the literature, in particular for delayed neural recurrent loops 8, 9, and experimentally in electrical circuits 9 and in recurrently clamped neurons 10. Foss et al. 8 studied neural recurrent inhibitory loops using the well-known Hodgkin- Huxley model and found three coexisting attracting periodic solutions by computer simulation. Ma and Wu 17, 18 showed that the phenomenological spiking neuron model incorporating the firing process and the absolute refractory period can generate a large number of asymptotically stable periodic solutions with predictable patterns of oscillations. However, the aforementioned work 8, 17, 18 on inhibitory loops/feedbacks is based on the consideration of the dynamical behaviors of the single excitatory neuron in the two-neuron inhibitory loop by simplifying the effect of the inhibitory neuron on the excitatory neuron as an inhibitory self-feedback in the recurrent loops. In this paper, we study both the excitatory and inhibitory neurons and consider a system of coupled delay differential equations. We develop a systematical approach to rigorously analyze the mechanism for the coexistence of multiple periodic patterns in the recurrent inhibitory loop. 68

3 In particular, the interaction among the excitatory and inhibitory neurons, time lag, firing process and the absolute refractory period can generate some basic components of periodic patterns. This allows us to link the solution semiflow defined by the model to symbolic dynamics in which these basic components can be pinned together to form a large class of periodic patterns and the model exhibits rich dynamics in the form of the coexisting periodic patterns. Moveover, an important quantity, the average time of convergence to a periodic pattern, can help us to determine what kind of patterns can be potentially used for neural information transmission, object recognition, sensory and perception processing in the nervous system. The rest of this paper is organized as follows: we first formulate the two-dimensional integrateand-fire neuron model for the recurrent inhibitory loop which incorporates the firing process and the absolute refractory period in Section 2. In Section 3, we show that the interaction of time lag, inhibitory feedback, firing and absolute refractory period can generate four types of basic oscillations (V, W u, W d and W c ) for the excitatory neuron, and we describe in details the corresponding basic types of the oscillations for the inhibitory neuron. Section 4 is devoted to developing a theoretical framework to view the inhibitory feedback from the inhibitory neuron as a self-feedback of the excitatory neuron: we note that the inhibitory neuron, upon receiving an excitatory signal from the excitatory neuron, emits a spike with a delay. This delay is in addition to the synaptic delay. The size of this additional delay is associated with the specific type of oscillation of the inhibitory neuron, and is determined fully by the period of the corresponding type of oscillation of the excitatory neuron. In summary, we show that the inhibitory signal from the inhibitory neuron can be regarded as a self-feedback of the excitatory neuron with this additional delay. In Section 5, we develop general principles that determine how these building blocks of oscillations interact each other to generate multiple periodic patterns, and apply these general principles to an example in Section 6 to illustrate the periodic patterns we theoretically derive. Section 7 provides some numerical simulations to illustrate similar periodic patterns exhibited by the more realistic model of the recurrent loop with an alpha synaptic function. Some remarks and discussions are presented in the final section. 2. Model of Recurrent Inhibitory Loops and Simulation Results We consider the following normalized integrate-and-fire model of a recurrent inhibitory loop: { V E (t) = V E(t) F I (t) + I 0, V I (t) = V (2.1) I(t) + F E (t, τ), where I 0 is the external input (assumed to be a constant), F E (t, τ) describes the excitatory input from the excitatory neuron E to the inhibitory neuron I and F I (t) describes the inhibitory input from neuron I to neuron E. The dynamical behaviors of both excitatory and inhibitory neurons are subject to a resetting mechanism. Namely, once these potentials reach the firing threshold ϑ at a firing time t f, there is a spike followed by an absolute refractory period (t f, t f + T F R ) during which both excitatory and inhibitory neurons are unresponsive to any input, and their potentials are 69

4 reset to the after-potential V A. Here T F R is the sum of the firing process and the absolute refractory period. The firing time t f is defined by the threshold condition: V i (t f ) = ϑ and V i (t f ε) < ϑ, where i=(e,i) for sufficiently small ε > 0. Once neuron E fires, an excitatory postsynaptic potential (EPSP) is delivered to neuron I with a time lag τ. The excitatory input is given by { b if t t τ F E (t, τ) = f,e + τ, t τ f,e + τ + T EF ; 0 otherwise, (2.2) where T EF is the duration of the excitatory input and t τ f,e is the last firing time of neuron E prior to the time t τ. More precisely, t τ f,e = sup{s; s t τ, V E(s) = ϑ and V (s ε) < ϑ for sufficiently small ε > 0}. We assume that the excitatory input F E (t, τ) will cause neuron I to fire and trigger only one spike. Hence, we take b > ϑ. In order to simplify our analysis, we assume T EF is set appropriately so that once neuron I emits a spike, there will be no further impact of the EPSP on neuron I (If T EF is greater than such a value, the computation is slightly more complicated, but it does not change our final results). In turn, neuron I instantaneously delivers an inhibitory postsynaptic potential (IPSP) to inhibit the firing of neuron E. The inhibitory input is given by { a if t tf,i, t F I (t) = f,i + T IF ; 0 otherwise, (2.3) where T IF is the duration of the inhibitory input and t f,i is the firing time of neuron I. It is reasonable to assume that T F R < T IF < T, where T is the intrinsic period of neuron I, to be defined in Subsection 3.1. We also assume that a > I 0 V A so that the inhibitory input will cause the potential of neuron E to decline. However, if the inhibitory input is delivered during the firing period or the absolute refractory period of neuron E, the input has no effect on the membrane potential during these two periods. Note that we assume that neuron I delivers an IPSP to neuron E instantaneously (i.e. without a synaptic delay). Such a synaptic delay can be easily incorporated into the recurrent loop model, but we can also omit this synaptic delay by a simple change of variables Simulation Results To illustrate basic oscillations and the formation of periodic patterns exhibited by neuron E in the recurrent inhibitory loop, we briefly present some simulation results. The firing process is represented by a piecewise linear function (the membrane potential first increases from the firing threshold ϑ to the amplitude c of a spike, then decays to a reset potential V r ), and the followed absolute refractory period is represented by an exponential function (the potential to increase V r to the after-potential V A ). To obtain all possible periodic patterns, we have performed experiments using a variety of initial conditions. 70

5 4 (a) (a) MEMBRANE POTENTIAL (mv) (b) 2 (b) TIME (ms) TIME (ms) Figure 1: Two coexisting attracting periodic solutions for τ = 2.406ms when ϑ = 1mV: (a) (W u W d ) and (b) (W uu V ). Other parameters are I 0 = 1.45µA, a = 2.7µA, b = 6µA, T IF = 0.5ms, T F R = 0.35ms, V r = 0.3mV, V A = mV and amplitude of spikes is 3.5mV. The righthand side is the blow-up of solutions in a given period (not delay τ) to clearly illustrate the patterns of solutions (the dotted line represents the after-potential V A. 4 (a) (a) 2 MEMBRANE POTENTIAL (mv) 0 4 (b) (b) (c) (c) (d) (d) TIME (ms) TIME (ms) Figure 2: Four coexisting attracting periodic solutions when τ = 3.917ms and other parameters identical to those used in Figure 1: (a) (W u ), (b) (W uu W d V ), (c) (W uuu 2V ) and (d) (W uu W u V W u ). The dotted line in the right side represents the after-potential V A ). Figure 1 shows two coexisting attracting periodic solutions when τ = 2.406ms, ϑ = 1mV, I 0 = 1.45µA, a = 2.7µA, b = 6µA, T IF = 0.5ms, T F R = 0.35ms, V r = 0.3mV and V A = mV. The right-hand side is the blow-up of solutions in a given period (not delay τ) to clearly illustrate the corresponding patterns. Focusing on the downstream oscillatory part of action potentials of these patterns, we note four types of basic oscillations: two types of W -oscillations as shown in Figure 1 (a), denoted by W u -oscillation and W d -oscillation, a more complicated oscillation, denoted by W uu and a simple oscillation, denoted by V -oscillation in Figure 1 (b). Intuitively, u stands for up indicating that the membrane potential increases before 71

6 an inhibitory signal is delivered; d stands for down indicating that the membrane potential decreases due to the effect of an inhibitory signal; V is a natural pattern where the action potential is not affected by an inhibitory signal. These periodic patterns in Figure 1 can be expressed in the order of these oscillations: (a) (W u W d ) and (b) (W uu V ). Figure 2 shows four coexisting periodic patterns when τ = 3.917ms and other parameters are identical to those used in Figure 1: (a) (W u ), (b) (W uu W d V ), (c) (W uuu 2V ) and (d) (W uu W u V W u ). In the following sections, we shall define these basic oscillations and derive general principles of how periodic pattern formations develop. 3. Basic Oscillations In this section, we define basic oscillations of neuron E observed in our simulation, discuss the mechanisms for generating these oscillations and describe their implication for the oscillations exhibited by neuron I The Excitatory Neuron In the absence of recurrent inhibition (F I 0), a strong enough stimulus I 0 (I 0 > ϑ) will cause neuron E to emit a sequence of spikes, called a spike train. The period of such a spike train subject to firing and the absolute refractory period is called the intrinsic spiking period of neuron E, denoted by T. This period is the duration of two consecutive firing times and can be divided into two parts, T F R and T Aϑ, I0 V A T = T F R + T Aϑ, T Aϑ = log, (3.1) I 0 ϑ where T Aϑ is the time that the membrane potential of neuron E increases from the after-potential V A to the threshold ϑ. Then we consider the impact of recurrent inhibitory input F I (t) on the dynamical behavior of neuron E. We note that each time neuron I fires a spike, an inhibitory input instantaneously acts on neuron E. Both the number of delivered inhibitory signals and the arrival time of these signals in the cycle of an action potential (the duration between two consecutive firing times), play an important role in formatting an oscillation. We distinguish the number of the inhibitory signals delivered between two consecutive firing times in three cases: no inhibitory signal, only one signal and multiple signals, which will generate a V -oscillation, a W d -oscillation or W u -oscillation, a W c - oscillation, respectively. We use Figure 3 (a) to illustrate these basic oscillations (V, W d, W u, W c ) corresponding to three cases. No Inhibitory Signal If there is no inhibitory signal delivered between two consecutive firing times, neuron E emits a natural spike, called a V -oscillation, illustrated by the graph from A to C in Figure 3 (a). The period of a V -oscillation is equal to the intrinsic spiking period T of neuron E, where T F R is the interval between A and B, T Aϑ is the interval between B and C in Figure 3 (a). Only One Inhibitory Input 72

7 V E V I 3 (a) Feedback Feedbacks Feedbacks Feedback 2 A C F J O U 1 M S 0 D H K N Q T L R B E G I P V (b) A 2 B 2 C 2 D 2 E 2 F 2 G 2 H 2 J 2 K TIME (ms) O 2 P 2 U 2 V 2 Figure 3: (a) Basic types of oscillations of the excitatory neuron and (b) the corresponding fundamental oscillations of the inhibitory neuron. In Figure 3 (a), from A to C is a V -oscillation; from C to F is a W d -oscillation; from F to J is a W u -oscillation; from J to O is a W du -oscillation; from O to U is a W uu -oscillation. Each time the excitatory neuron E fires a spike, an excitatory signal (EFB) is delivered to the inhibitory neuron I in a time lag τ. In Figure 3 (b), the excitatory signals delivered at A 2, C 2, F 2, J 2, O 2 and U 2 correspond to the spikes at A, C, F, J, O and U in (a), respectively. In order to associate with each spike, we moved the position of V I backwards by an amount of time delay τ. The top dotted line indicates the threshold ϑ and the bottom line indicates the after-potential V A. If there is only one inhibitory input delivered between two firing times, and if this inhibitory signal is delivered during the firing process or the absolute refractory period (between C and D in Figure 3 (a)), since T IF > T F R, this inhibitory signal partially impacts on the membrane potential: after the absolute refractory period, the membrane potential first decreases in the effective inhibitory time (from D to E), denoted by t down (ms), and then increases until the neuron fires. We call such an oscillation a W d -oscillation, depicted from C to F in Figure 3 (a) ( d stands for down ). On the other hand, if this inhibitory signal is delivered after the absolute refractory period (at H in Figure 3 (a)), this signal fully impacts on the membrane potential: after the absolute refractory period, the membrane potential first increases in a certain amount of time before the signal arrives, denoted by t up (ms) (from G to H), then decreases due to the effect of the inhibitory signal and finally increases until the neuron fires. We call such an oscillation a W u -oscillation, depicted from F to J in Figure 3 (a) ( u stands for up ). The first quantity characterizing these two oscillations (W d and W u ) is the effective inhibitory time, t down, for a W d -oscillation and the rising time, t up, prior to arrival of the signal for a W u - oscillation. The second quantity to describe the inhibition of the signal is a time shift to the original cycle of the action potential. We denote by t Aϑ the duration from the time when the inhibition of 73

8 the signal wears off to the next firing time, depicted by the time from E to F for the W d -oscillation and from I to J for the W u -oscillation in Figure 3 (a). The time shift is t down + t Aϑ T Aϑ for a W d -oscillation and t up + T IF +t Aϑ T Aϑ for a W u -oscillation. To simplify the notation, we denote t Aϑ T Aϑ by d or t u corresponding to a W d or W u oscillation. Then each W d -oscillation can be characterized by a pair of variables (t down, t d ) with T IF T F R < t down T IF and t d > 0, and its period is T d = T + t down + t d. Each W u -oscillation can be characterized by a triple of variables (t up, T IF, t u ) with 0 < t up < T Aϑ, and its period is T u = T + t up + T IF + t u. The sign of t u depends on the value t up : if the signal inhibition brings the membrane potential below V A, t u > 0; otherwise, t < 0. Multiple Inhibitory Signals There can be multiple inhibitory signals delivered between two consecutive firing times, which will generate more complicated oscillations, such as a W du -oscillation from J to O and a W uu - oscillation from O to U in Figure 3 (a). Similarly to W d and W u notations, the arrival time of the first signal determines either a d-stream characterized by t c down or a u-stream characterized by (t c up, T IF ). The followed signals are all of u-streams characterized by (t c up, T IF ). In Figure 3 (a), the W du -oscillation results from two signals: the first signal delivered between J and K decreases the potential in an effective inhibitory duration t c down (from K to L) and the second signal is delivered at M characterized by (t up, T IF ), corresponding to the segments from L to M and from M to N. The W uu -oscillation is the effect of two signals delivered at Q and S, characterized by (t c,1 up, T IF ) and (t c,2 up, T IF ), corresponding to each segment between P and T. We call this type of oscillation a W c -oscillation or Wc m -oscillation ( c for complicated oscillation ) with c = du...u or c = u...u, where m is the total number of d and u s in c and it is the number of inhibitory signals. In general, a Wdu...u m -oscillation can be characterized by (tc down, tc,1 up, T IF,..., t c,m 1 up, T IF, t c ), where T IF T F R < t c down T IF, and its period is T +t c down + m 1 i=1 (tc,i up+t IF )+ t c ; a Wu...u-oscillation m can be characterized by (t c,1 up, T IF,..., t c,m up, T IF, t c ) and its period is T + m i=1 (tc,i up + T IF ) + t c, where t c is defined in the same way as d and u described above The Inhibitory Neuron Each time neuron E fires a spike, an excitatory signal is delivered to neuron I with a time lag τ. In contrast to neuron E, we define a cycle of the membrane potential of neuron I between two consecutive excitatory signals as a fundamental oscillation of neuron I. The excitatory signal delivered by a V -oscillation of neuron E generates a fundamental oscillation of neuron I, denoted by V I -oscillation, which is the graph from A 2 to C 2 in Figure 3 (b). Similarly, the spike of W d - oscillation generates a fundamental W I d -oscillation (from C 2 to F 2 in Figure 3 (b)); the spike of W u -oscillation generates a fundamental W I u -oscillation (from F 2 to J 2 in Figure 3 (b)); the spike of W c -oscillation generates a fundamental W I c -oscillation, such as a W I du -oscillation (from J 2 to O 2 ) and a W I uu-oscillation (from O 2 to U 2 ). It follows easily from the above discussion that the period of each fundamental oscillation of neuron I is equal to the period of its corresponding basic oscillation of neuron E. 74

9 4. Pattern Formation From Four Basic Oscillations The interaction of inhibitory signals, the time lag, firing and the absolute refractory period can generate four basic oscillations (V, W d, W u, W c ), which are the basic building blocks of periodic patterns exhibited by neuron E. In this section, we discuss the general principles of how these basic oscillations interact to generate a large class of periodic patterns of neuron E. If neuron E fires a spike at t f, neuron I receives an excitatory signal in a time lag τ due to such a spike, which excites neuron I to fire. Once neuron I emits a spike, an inhibitory signal is immediately delivered to neuron E. Such an inhibitory signal can be considered as a self-feedback (SFB) of neuron E due to the spike at t f. For the inhibitory neuron, we denote the duration from the time that an excitatory signal arrives until neuron I fires by τ A, which is the segment, A 2 B 2, or C 2 D 2, or F 2 G 2, or J 2 K 2, or O 2 P 2, or U 2 V 2 in Figure 3 (b). In contrast to the synaptic delay τ, we call τ A an additional delay. Hence, corresponding to each spike of the excitatory neuron, there is always an inhibitory self-feedback (SFB) to act on the excitatory neuron in a time delay τ = τ + τ A. We shall show that τ A only depends on the aforementioned basic types of oscillations of the excitatory neuron when a periodic solution is generated. The advantage of introducing this additional delay is to help us to regard the IPSP as a self-feedback of the excitatory neuron so that we can focus on analyzing the dynamics of the excitatory neuron only. For a given initial condition, the effect of such self-feedbacks eventually stabilizes the corresponding trajectory of neuron E to a periodic solution with a predicted pattern composed of the aforementioned basic oscillations. Then a periodic pattern can be expressed by σ = (π 1, π 2,..., π N 1, π N ) where π i {V, W d, W u, W c }. (4.1) We denote by T i := T (π i ) the period of the basic oscillation π i and by p := N i=1 T i the period of such a pattern. We note that a pattern which is the permutation of these oscillations π 1,..., π N is equivalent to the above pattern, i.e., σ 1 = (π 2, π 3,..., π N, π 1 ) is equivalent to σ The Inhibitory Neuron We denote the fundamental oscillations of the inhibitory neuron generated by spikes of π 1, π 2,..., π N by π1, I π2, I..., πn I, respectively. Additional delays corresponding to these fundamental oscillations π1, I π2, I..., πn I are denoted by τ 1, τ 2,..., τ N, respectively. In what follows, we show that τ 1,..., τ N can be determined by T 1,..., T N via a single scalar function. Theorem 1. For a given periodic pattern σ = (π 1, π 2,..., π N 1, π N ) with π i {V, W d, W u, W c }, there exists a continuous function g such that τ 1 = g(t 1, T 2,..., T N ), τ 2 = g(t 2, T 3,..., T 1 ),..., τ N = g(t N, T 1,..., T N 1 ). Proof: We first illustrate our argument by a simple case where σ = (π 1, π 2 ). We denote the membrane potential of the inhibitory neuron by V1 I at the time when the excitatory signal due to the spike π 1 is delivered, and by V2 I at the time when the excitatory signal due to the spike π 2 is delivered. 75

10 Now we discuss the time course of membrane potential of the fundamental oscillation π1. I The potential evolves in three steps: (i) increasing from V1 I to the threshold; (ii) the firing process and the absolute refractory period; (iii) increasing from V A to V2 I. We follow these steps to calculate V1 I and V2 I. First Step: Starting from the time when the excitatory signal due to the spike π 1 is delivered, it takes τ 1 (ms) for the potential to increase from V1 I to the firing threshold ϑ. Integrating the second equation in (2.1) yields ϑ = V I 1 e τ 1 + b(1 e τ 1 ). (4.2) Second Step: After neuron I fires a spike and passes the absolute refractory period, the potential reaches the after-potential V A. Third Step: After firing and the absolute refractory period, it takes the duration of T 1 τ 1 T F R for the potential to increase from the after-potential V A to V2 I. System (2.1) gives rise to V I 2 = V A e (T 1 τ 1 T F R ). (4.3) Similarly, we consider the fundamental oscillation π I 2 and obtain two equations analogous to equations (4.2) and (4.3). Combining these equations together, we obtain ϑ = V1 I e τ 1 + b(1 e τ 1 ), V1 I = V A e (T 2 τ 2 T F R ), ϑ = V2 I e τ 2 + b(1 e τ 2 ), V2 I = V A e (T 1 τ 1 T F R ). We solve τ 1, τ 2 to obtain τ 1 = log τ 2 = log b(b ϑ V A e T 2 +T F R ) (b ϑ) 2 VA 2e T 1 T 2 +2T F R b(b ϑ V A e T 1 +T F R ) (b ϑ) 2 V 2 A e T 1 T 2 +2T F R,. (4.4) Therefore, the function g is given by b(b ϑ VA e x 2+T F R ) g(x 1, x 2 ) = log (b ϑ) 2 VA 2e x 1 x 2 +T F R. Clearly, τ 1 = g(t 1, T 2 ) and τ 2 = g(t 2, T 1 ). Applying the same argument to the periodic pattern σ = (π 1, π 2, π 3 ), we obtain { } b(b ϑ)(b ϑ V τ 1 = log A e T 3 +T F R )+VA 2e T 2 T 3 +2T F R, (b ϑ) 3 +VA 3e T 1 T 2 T 3 +3T F R { } b(b ϑ)(b ϑ V τ 2 = log A e T 1 +T F R )+VA 2e T 1 T 3 +2T F R, (b ϑ) 3 +VA 3e T 1 T 2 T 3 +3T F R { } τ 3 = log. b(b ϑ)(b ϑ V A e T 2 +T F R )+V 2 A e T 1 T 2 +2T F R (b ϑ) 3 +V 3 A e T 1 T 2 T 3 +3T F R (4.5) 76

11 The function g is given by { b(b ϑ)(b ϑ VA e x 3+T F R ) + VA 2 g(x 1, x 2, x 3 ) = log e x 2 x 3 +2T F R } (b ϑ) 3 + VA 3e x 1 x 2 x 3 +3T F R and τ 1 = g(t 1, T 2, T 3 ), τ 2 = g(t 2, T 3, T 1 ) and τ 3 = g(t 3, T 1, T 2 ). This argument can apply for the general periodic pattern σ = (π 1, π 2,..., π N 1, π N ). This completes our proof. Therefore, there exists a single continuous function to describe the relationship of each τ i and T 1,..., T N. The observed cyclic property can be interpreted as follows: if τ 1 is considered as the function g with respect to the periodic pattern σ = (π 1, π 2,..., π N ), then τ i is the function g with respect to the equivalent periodic pattern σ i = (π i, π i+1,..., π i 1 ). In other words, the additional delay τ i can be expressed in terms of the function g with respect to the equivalent periodic pattern. In what follows, we use notations of τ v, τ u, τ d, τ c (such as τ du, τ uu ) to denote the additional delays corresponding to the V -oscillation, W u -oscillation, W d -oscillation, W c -oscillation (W du, W uu ), respectively The Excitatory Neuron After introducing the self-feedback (SFB) and establishing the relationship of τ 1,..., τ N and T 1,..., T N, we can now concentrate on the excitatory neuron to study periodic patterns. We present the general principles to determine how the aforementioned basic oscillations are interacted to form a periodic pattern. For a given periodic pattern σ = (π 1, π 2,..., π N 1, π N ) where π i {V, W d, W u, W c }, we consider an action potential π i which delivers a self-feedback (SFB) to the action potential π i+r with 0 R N 1. The action of the self-feedback (SFB) from π i on π i+r yields i+r 1 k=i T k + T F R τ + τ i + T IF and π i = W c iff π i+r = V ; i+r 1 k=i T k + T F R + t down = τ + τ i + T IF iff π i+r = W d ; i+r 1 k=i T k + T F R + t up = τ + τ i iff π i+r = W u ;. i+r 1 k=j T k + T F R + t c down = τ + τ i + T IF iff d of π i+r = Wdu...u m ; i+r 1 k=i T k + T F R + t c down + h j=1 (tc,j up + T F D ) = τ + τ i + T IF iff the h-th u of π i+r = Wdu...u m ; T k + T F R + h j=1 tc,j up = τ + τ i iff the h-th u of π i+r = Wu...u, m i+r 1 k=i where τ i is the additional delay corresponding to the basic oscillation π i, given by τ i = g(t i, T i+1..., T i 1 ). Analogously, N 1 other patterns equivalent to σ give N 1 relationships, which determine the types of periodic patterns generated by the recurrent inhibitory loop for a given time delay. Furthermore, these relationships determine the minimum and maximum values of τ for the existence of a given periodic pattern. Now we consider the possible value of R to determine all possible periodic patterns exhibited by the excitatory neuron. According to the value of R, we distinguish periodic patterns in four cases, as shown in Figure 4: 77 (4.6)

12 MEMBRANE POTENTIAL (mv) A B far range inhibition nearest neighborhood inhibition self inhibition G E C H F D I J K L M N O threshold after potential V A TIME Figure 4: Schematic of three types of periodic patterns composed of a W c -oscillation (from B to G). When the excitatory neuron fires a spike at B, an inhibitory self-feedback is delivered at a time lag τ = τ + τ A. If this feedback is delivered to the W c -oscillation itself (at E), it is called a self-inhibitory pattern; if the feedback is delivered to its nearest neighbor of W c -oscillation (at H), it is called a nearest-neighbor-inhibitory pattern; otherwise, a far-range-inhibitory pattern is generated (for example, the feedback arrives at K). R = N 1: the first item π 1 always acts on the last item π N within the pattern σ and this periodic pattern is composed of only W u and W d oscillations; R = 0 for π i = W c : the spike of W c -oscillation delivers a self-feedback to itself and we call it a self-inhibitory pattern; R = 1 for π j = W c : the spike of W c -oscillation delivers a self-feedback to its nearest neighbor and we call it a nearest-neighbor-inhibitory pattern; R > 1 for π i = W c : the spike of W c -oscillation delivers a self-feedback to an oscillation far away from its nearest neighbor and we call it a far-range-inhibitory pattern. 5. Periodic Patterns Exhibited by the Excitatory Neuron In this section, we discuss four types of periodic patterns in details: periodic patterns composed of only W u -oscillations and W d -oscillations, self-inhibitory periodic patterns, nearest-neighborinhibitory periodic patterns and far-range-inhibitory patterns. 78

13 5.1. Periodic Patterns Composed of W u -oscillations and W d -oscillations If R = N 1 for all equivalent periodic patterns of σ = (π 1, π 2,..., π N ), then the first item π 1 always acts on the last item π N within the pattern σ. It follows easily from the definition in Subsection 3.1 that neither the V -oscillation nor the W c -oscillation can appear in such a periodic pattern. Lemma 2. For a given periodic pattern σ = (π 1, π 2,..., π N ) where each π i is one of the basic oscillation, if R = N 1 for all equivalent periodic patterns of σ, then π i is either a W u -oscillation or a W d -oscillation. Examples of periodic patterns composed of only W u -oscillation and W d -oscillation are (W u W d ) in Figure 1 (a) and (1W u ) in Figure 2 (a). Lemma 3. For a given periodic pattern σ = (π 1, π 2,..., π N ) where π i {W u, W d }, all W u - oscillations must be the same in terms of (t up, T IF, t u ) and all W d -oscillations must be the same in terms of (t down, t d ). Proof: We illustrate our argument by a simple case, where σ = (π 1, π 2 ) = (W u (1), W u (2) ). π 1 is characterized by (t (1) up, T IF, t (1) u ) and π 2 is characterized by (t (2) up, T IF, t (2) u ). Condition (4.6) gives rise to { T 2 + T F R + t (1) up = τ + τ 2, T 1 + T F R + t (2) up = τ + τ 1. The above system yields T 1 t (1) up τ 1 = T 2 t (2) up τ 2. Substituting T 1, T 2 into the above equation, we obtain that t (1) u τ 1 = t (2) u τ 2. The potential of a W u -oscillation evolves in three steps after the firing and the absolute refractory period: (i) increases from the after-potential V A to V 1 in the duration of t up (ms); (ii) then decreases from V 1 to V 2 in the duration of T IF (ms) due to an inhibitory feedback; (iii) after the feedback wears off, the potential increases from V 2 to the threshold ϑ in the duration of t Aϑ = T Aϑ + t u. These three steps give rise to V 1 = I 0 (I 0 V A )e t up, V 2 = I 0 a (I 0 a V 1 )e t up, I0 V 2 T Aϑ + t u = log. I 0 ϑ Combining the above three equations and using equation (3.1) of T Aϑ, we obtain (I 0 V A )e t u e T IF t up = a(1 e T IF ). 79

14 For both the W u (1) -oscillation and W u (2) -oscillation, we have (I 0 V A )e t(1) u e T IF t (1) up = a(1 e T IF ), (I 0 V A )e t(2) u e T IF t (2) up = a(1 e T IF ). Combining the above two equations with t (1) u τ 1 = t (2) u τ 2, we obtain e τ 2 (e T F R e T 1 ) = e τ 1 (e T F R e T 2 ). Substituting τ 1 = g(t 1, T 2 ) and τ 2 = g(t 2, T 1 ) of equation (4.4), we obtain b(b ϑ VA e T 1+T F R ) (b ϑ) 2 VA 2e T 1 T 2 (e T F R e T 1 b(b ϑ VA e T 2+T F R ) ) = +2T F R (b ϑ) 2 VA 2e T 1 T 2 +2T F R (e T F R e T 2 ). It follows easily from the above equation that T 1 = T 2. Hence, t (1) up = t (2) up and t (1) u = t (2) u. The same argument with τ i = g(t i, T i+1,..., T i 1 ) can apply to show that all W u -oscillations in a periodic pattern must be the same in terms of (t up, T IF, t u ) and all W d -oscillations must be the same in terms of (t down, t d ). Lemma 4. For a given periodic pattern σ = (π 1, π 2,..., π N ) where π i {W u, W d }, all additional delays τ i related to the W u -oscillations must be the same and all additional delays τ i related to the W d -oscillations must be the same. Lemma 5. For a given periodic pattern σ = (π 1, π 2,..., π N ) where π i {W u, W d }, if a W u - oscillation and a W d -oscillation coexist in this periodic pattern, then t u τ u = t d τ d, where τ u, τ d are the additional delays corresponding to the W u -oscillation and the W d -oscillation, respectively. Proof: Again, we illustrate our argument by a simple case where σ = (π 1, π 2 ) = (W u, W d ). π 1 is characterized by (t up, T IF, t u ) and π 2 is characterized by (t down, t d ). Condition (4.6) gives rise to { T1 + T F R + t down = τ + τ u + T IF, The above system yields T 2 + T F R + t up = τ + τ d. T 1 t up τ u = T 2 t down + T IF τ d. Substituting T 1, T 2 into the above equation, we obtain that t u τ u = t d τ d.. The relationship of t down and t d for a W d -oscillation is given by I 0 V A a t down = log. (5.1) (I 0 V A )e t d a 80

15 The relationship of t up and t u for a W u -oscillation is given by I 0 V A t up = log (I 0 V A )e T IF + t u. (5.2) ae T IF + a When t down = T IF for a W d -oscillation and t up = 0 for a W u -oscillation, the W d -oscillation coincides with the W u -oscillation ( t u = t d, and correspondingly τ u = τ d ) and two action potentials are the same. In such a case, t d and t u reach their maximum values, denoted by t max = log 1 + a(et IF + 1) T IF > 0. (5.3) I 0 V A For the W d -oscillation, when t down = T IF T F R, t d reaches its minimum value t min d (I0 V A a)e (T IF T F R ) + a t min d = log I 0 V A given by > 0. (5.4) Theorem 6. The periodic pattern (1W u ) for σ = (kw u ) with k 2 can be generated if and only if { tup = τ+τ u (k 1)T +T F R +(k 1)T IF +(k 1) t u, k (5.5) 0 < t up < T Aϑ, where τ u is the additional delay of the W u -oscillation. The minimum and maximum values of τ for the existence of such a periodic pattern are given by { τmin = (k 1)T + T F R + (k 1)(T IF + t max ) τ u tmax, (5.6) τ max = (k 1)T + T F R + (k 1)T IF + kt Aϑ + (k 1) t min u τ u t min u, where τ u tmax is the value of τ u subject to the condition where t up = 0 and t u = t max, and τ u t min is the value of τ u u subject to the condition where t up = T Aϑ and t u = t min u. Here t min I0 ϑ + a(e T IF 1) u = log T IF < 0. (5.7) I 0 V A Proof: We consider σ = (π 1,..., π k ) = (W u,..., W u ) with the W u -oscillation being characterized by (t up, T IF, t u ). Condition (4.6) gives rise to (k 1)T + (k 1)(t up + T IF + t u ) + T F R + t up = τ + τ u. Re-arranging the above formula yields equation (5.5) and t up must satisfy 0 < t up < T Aϑ. When t up = 0, t u reaches its maximum value t max given by equation (5.3). The corresponding additional delay τ u is denoted by τ u tmax. This situation gives rise to the minimum value τ min of equation (5.5). of equation. This situation yields the On the other hand, as t up approaches T Aϑ, equation (5.2) yields the value t min u (5.7). The corresponding additional delay τ u is denoted by τ u t min u maximum value τ max of equation (5.5). 81

16 Now we apply the above theorem to the simple pattern σ = (2W u ), as shown in Figure 2 (a). The period of a W u -oscillation, condition (4.6), the additional time delay (4.4), and equation (5.2) yield T u = T + t up + T IF + t u, T u + T F R + t up = τ + τ u, b(b ϑ VA e T u+t F R ) τ u = log (b ϑ) 2 VA 2e 2T u+2t F R t up = log, I 0 V A (I 0 V A )e T IF + t u ae T IF + a For any given τ, the solution of variables T u, t up, t u and τ u uniquely determines the periodic pattern. On the other hand, we can calculate the minimum value τ min subject to the condition where t up = 0 and the maximum value τ max subject to the condition where t up = T Aϑ. The segment between A and B in Figure 5 represents the interval for the existence of the periodic pattern (1W u ) with σ = (2W u ). The W d -oscillation always transits to the W u -oscillation at certain values of time delay τ. This transition implies that any periodic pattern in the case where R = N 1 contains at most one W d -oscillation. We now discuss the periodic patterns (1W d, hw u ) with h 1. Condition (4.6) immediately gives rise to Theorem 7. (i) The periodic patterns (1W d, hw u ) with h 1 can be generated if and only if { τ = ht + TF R + (h 1)T IF + h(t up + t u ) + t down τ u ; (5.8) t u τ u = t d τ d. The minimum value of τ for the existence of such a periodic pattern is given by τ min = NT + T F R + (h 1)T IF + min t min d t d t max h(t up + t u ) + t down τ u, subject to the condition t u τ u = t d τ d. The maximum value of τ for the existence of such a periodic pattern is given by τ max = ht + T F R + h(t IF + t max ) τ u tmax, where τ u tmax is defined in Theorem 6. (ii) If τ reaches its minimum value at t = t max, such periodic patterns can not be generated. (iii) If τ reaches its minimum at t < t max, then a pattern transition occurs at τ max. Now we apply the above theorem to the periodic pattern σ = (W u, W d ), as shown in Figure 1 (a). The periods of the W u -oscillation and the W d -oscillation, condition (4.6), the additional time. 82

17 p/t H I J K 4 G F 3 E C D 2 B A τ/t Figure 5: Plot of p/t versus τ/t of six periodic patterns where p is the period of a pattern: (W u ) with σ = (2W u ) from A to B; (W u W d ) from C to D; (W du V ) from E to F ; (W uu V ) from F to G; (W du W d V ) from H to I; and (W uu W d V ) from J to K. The periodic pattern (W du V ) transits to the periodic pattern (W uu V ) at F and the periodic pattern (W u W d ) transits to the periodic pattern (W u ) with σ = (2W u ) at the dotted line. delay (4.4), equations (5.1) and (5.2) yield T u = T + t up + T IF + t u, T d = T + t down + t d, T d + T F R + t up = τ + τ d, t u τ u = t d τ d, b(b ϑ VA e T d+t F R ) τ u = log (b ϑ) 2 VA 2e T d T u, +2T F R b(b ϑ VA e Tu+T F R ) τ d = log (b ϑ) 2 VA 2e T d T u, +2T F R I 0 V A t up = log (I 0 V A )e T IF + t u ae T IF + a, t down = log I 0 V A a (I 0 V A )e t d a For any given τ, the solution of the above system uniquely determines the periodic pattern. On the other hand, we can calculate the minimum value τ min subject to t down = T IF T F R and the maximum value τ max subject to t down = T IF. The segment between C and D in Figure 5 represents the interval for the existence of the periodic pattern (W u, W d ). The periodic pattern (W u, W d ) transits to the periodic pattern (W u ) with σ = (2W u ) at the dotted line in Figure Self-Inhibitory Periodic Patterns Self-inhibitory periodic patterns are composed of a W m c -oscillation and (m 1) V -oscillations, so these periodic patterns can be expressed as (W m du..u, (m 1)V ) or (W m u..u, (m 1)V ). The spike emitted by the W c -oscillation along with spikes emitted by the V -oscillations generate the. 83

18 W c -oscillation. Examples of self-inhibitory periodic patterns are (W uu V ) in Figure 1 (b) and (W uuu 2V ) in Figure 2 (c). Lemma 8. For the self-inhibitory pattern (W m c, (m 1)V ), the following holds: (i) t c,1 up =... = t c,m 2 up = T T IF, t c,m 1 up = T T IF + τ v τ c for a Wdu...u m -oscillation, and the period of Wdu...u m (m) -oscillation is T du...u = mt + tc down + t c + τ v τ c ; (ii) t c,2 up =... = t c,m 1 up = T T IF, t c,m up = T T IF + τ v τ c for a W m u...u-oscillation, and the period of W m u...u-oscillation is T (m) u...u = mt + t c,1 up + t c + τ v τ c. Proof: We consider the pattern σ = (π 1, π 2,..., π m ) = (W c, (m 1)V ). For the W c -oscillation, the evolution of the potential is characterized by (t c,i up, T IF ) between two consecutive self-feedbacks due to π i and π i+1. Without the consideration of the additional delay, we have t c,i up +T IF = T due to the oscillation π i = V. The effect of the additional delay leads to t c,i up + T IF + τ i+1 τ i = T. In the situation where π i = π i+1 = V, we still have t c,i up + T IF = T since τ i = τ i+1. On the other hand, in the situation where π i = V and π i+1 = W c (it happens when i = m 1 for the Wdu...u m -oscillation and i = m for the Wu...u-oscillation), m we have t c,i up + T IF + τ c τ v = T. This completes the proof. We now consider the periodic pattern (W du V ). The W du -oscillation is characterized by (t c down, T T IF + τ du τ v, T IF, t c ) and the corresponding additional delay is denoted by τ du. Following the evolution of the potential of the W du -oscillation (similarly to the calculation in Lemma 3), we establish the relationship of t c down and t c. t c = log e T tc down τ du+τ v + a 1 e T IF + e T τ du+τ v (1 e tc down ). I 0 V A When the self-feedback of the W du -oscillation is delivered, the potential is given by V 2 = I 0 (I 0 V A )e tc down + a(1 e t c down )e T +T IF τ du +τ v. (5.9) The period of W du -oscillation, condition (4.6) and additional delays (4.4) yield the following system T du = 2T + t c down + t c + τ v τ du, (5.10) T F R + t c down + T T IF = τ + τ du, (5.11) b(b ϑ V A e T +T F R ) τ du = log, (5.12) (b ϑ) 2 VA 2e T T du+2t F R b(b ϑ VA e T du+t F R ) τ v = log, (5.13) (b ϑ) 2 VA 2e T T du+2t F R t c = log e T tc down τ du+τ v + a 1 e T IF + e T τ du+τ v (1 e tc down ). (5.14) I 0 V A 84

19 Lemma 9. For the self-inhibitory periodic pattern (W du, V ), the minimum value of τ for the existence of such a periodic pattern is τ min = T +T F R T IF +(t c down τ du) min, where (t c down τ du) min is the minimum value of t c down τ du of system ( ) subject to t c down = T IF T F R and the condition where V 2 = ϑ of equation (5.9); the maximum value is τ max = T +T F R (τ du ) t c down =T IF where (τ du ) t c down =T IF is the solution of τ du of system ( ) subject to t c down = T IF. On the other hand, for any given value τ (τ min, τ max ), the solution of system ( ) uniquely determines the periodic pattern. Furthermore, the periodic pattern (W du, V ) transits to the periodic pattern (W uu, V ) at τ max. We then consider the periodic pattern (W uu V ). The W uu -oscillation is characterized by (t c,1 up, T IF, T T IF + τ uu τ v, T IF, t c ) and the corresponding additional delay is denoted by τ uu. Following the evolution of the potential of the W uu -oscillation, we establish the relationship of t c,1 up and t c. obtain t c = log e T T IF t c,1 up τ uu+τ v + a(1 e T IF )(1 + e T τ uu+τ v ) I 0 V A When the self-feedback of the W du -oscillation is delivered, the potential is given by V 2 = I 0 (I 0 V A )e tc,1 up + a(e T IF 1)e T τ uu+τ v. (5.15) The period of W uu -oscillation, condition (4.6), additional time delays (4.4) yield the following system T uu = 2T + t c,1 up + T IF + t c + τ v τ uu, (5.16) T F R + t c,1 up + T = τ + τ uu, (5.17) b(b ϑ V A e T +T F R ) τ uu = log, (5.18) (b ϑ) 2 VA 2e T T uu+2t F R b(b ϑ VA e Tuu+T F R ) τ v = log, (5.19) t c = log (b ϑ) 2 VA 2e T T uu+2t F R e T T IF t c,1. up τ uu +τ v + a(1 e T IF )(1 + e T τuu+τv ) I 0 V A. (5.20) Lemma 10. For the self-inhibitory periodic pattern (W uu, V ), the minimum value of τ for the existence of such a periodic pattern is τ min = T + T F R (τ uu ) t c,1 up =0, where (τ uu) t c,1 up =0 is the solution of τ uu of system (5.16) (5.20) subject to t c,1 up = 0; the maximum value of τ is τ max = T + T F R + (t c,1 up τ uu ) V2, where (t c,1 up τ uu ) V2 is the value of t c,1 up τ uu of system (5.16) (5.20) subject to V 2 = ϑ in equation (5.15). On the other hand, for any given value τ (τ min, τ max ), the solution of system ( ) uniquely determines the periodic pattern. The segment between E and F in Figure 5 represents the interval for the existence of the periodic pattern (W du, V ). The segment between F and G in Figure 5 represents the interval for the existence of the periodic pattern (W uu, V ). The periodic pattern (W du, V ) transits to the periodic pattern (W uu, V ) at F. Some other self-inhibitory periodic patterns such as (W duu, 2V ) and (W uuu, 2V ) will be presented in our case study. 85

20 5.3. Nearest-Neighbor-Inhibitory Periodic Patterns The simple nearest-neighbor-inhibitory periodic patterns are (Wc m, W d, (m 1)V ) with m 2 such as (W du W d V ) and (W uu W d V ). For the periodic pattern (Wc m, W d, (m 1)V ), the inhibitory self-feedback due to the Wc m -oscillation generates the W d -oscillation that in turn, along with the V -oscillations, generates the Wc m -oscillation. Examples of nearest-neighbor-inhibitory periodic patterns are (W uu W d V ) and (W uu W u V W u ) in Figure 2. We first consider the periodic pattern (W du W d V ). The W d -oscillation is characterized by (t down, t d ) and the W du -oscillation is characterized by (t c down, tc,1 up, T IF, t c ). For the W du -oscillation, when the self-feedback due to the V -oscillation arrives, the potential is given by The W du -oscillation yields V 2 = I 0 (I 0 V A )e tc,1 up t c down a(1 e t c down )e t c,1 up. (5.21) T du = T + t c down + t c,1 up + T IF + t c, (5.22) t c = log e T IF t c down tc,1 up + a 1 e T IF + (1 e tc down )e t c,1 up T IF. (5.23) I 0 V A The W d -oscillation leads to Condition (4.6) gives rise to T d = T + t down + t d, (5.24) I 0 V A a t down = log. (5.25) (I 0 V A )e t d a 2T + t down + t d + T F R + t c down = τ + τ d + T IF, (5.26) t c,1 up = T + t down + t d T IF τ d + τ v, (5.27) t c + t down τ du + τ v = 0. (5.28) Additional delays are given by { b(b ϑ)(b ϑ VA e T +T F R ) + VA 2 τ du = log e T d T +2T F R }, (5.29) (b ϑ) 3 + VA 3e T du T d T +3T F R { b(b ϑ)(b ϑ VA e T du+t F R ) + VA 2 τ d = log e T du T +2T F R }, (5.30) (b ϑ) 3 + VA 3e T du T d T +3T F R { b(b ϑ)(b ϑ VA e T d+t F R ) + VA 2 τ v = log e T du T d +2T F R }. (5.31) (b ϑ) 3 + VA 3e T du T d T +3T F R Lemma 11. For the nearest-neighbor-inhibitory periodic pattern (W du W d V ), the minimum value of τ for the existence of such a periodic pattern is τ min = 2T +T F R T IF +(t down +t c down + t d τ d ) max, where (t down + t c down + t d τ d ) max is the maximum value of t down + t c down + t d τ u 86

21 of system (5.22) (5.31) subject to t c down = T IF T F R or t c down = T IF or V 2 = ϑ of equation (5.21); the maximum value of τ is τ max = T + T F R (t down + t c down + t d τ d ) min, where (t down + t c down + t d τ d ) min is the minimum value of t down + t c down + t d τ d of system (5.22) (5.31) subject to t c down = T IF or t down = T IF T F R. On the other hand, for any given value τ (τ min, τ max ), the solution of system (5.22) (5.31) uniquely determines the periodic pattern. We then consider the periodic pattern (W uu W d V ). The W d -oscillation is characterized by (t down, t d ) and the W uu -oscillation is characterized by (t c,1 up, T IF, t c,2 up, T IF, t c ). For the W uu - oscillation, when the self-feedback due to the V -oscillation arrives, the potential is given by The W uu -oscillation yields The W d -oscillation leads to Condition (4.6) gives rise to V 2 = I 0 (I 0 V A )e tc,1 up t c,2 up T IF a(1 e T IF )e tc,2 up. (5.32) T uu = T + t c,1 up + t c,2 up + 2T IF + t c, (5.33) t c = log e 2T IF t c,1 up tup c,2 + a (1 e T IF )(1 + e tc,2 up T IF ). (5.34) I 0 V A T d = T + t down + t d, (5.35) I 0 V A a t down = log. (5.36) (I 0 V A )e t d a 2T + t down + t d + T F R + t c,1 up = τ + τ d, (5.37) t c,2 up = T + t down + t d T IF τ d + τ v, (5.38) t c + t down τ uu + τ v = 0. (5.39) Additional delays are given by { b(b ϑ)(b ϑ VA e T +T F R ) + VA 2 τ uu = log e T d T +2T F R }, (5.40) (b ϑ) 3 + VA 3e Tuu T d T +3T F R { b(b ϑ)(b ϑ VA e T uu+t F R ) + VA 2 τ d = log e T uu T +2T F R }, (5.41) (b ϑ) 3 + VA 3e Tuu T d T +3T F R { b(b ϑ)(b ϑ VA e T d+t F R ) + VA 2 τ v = log e T uu T d +2T F R }. (5.42) (b ϑ) 3 + VA 3e Tuu T d T +3T F R Lemma 12. For the the nearest-neighbor-inhibitory periodic pattern (W uu W d V ), the minimum value of τ for the existence of such a periodic pattern is τ min = 2T +T F R +(t down + t d τ d ) max, where (t down + t d τ d ) max is the maximum value of t down + t d τ d of system (5.33) (5.42) 87

Consider the following spike trains from two different neurons N1 and N2:

Consider the following spike trains from two different neurons N1 and N2: About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in

More information

An Introductory Course in Computational Neuroscience

An Introductory Course in Computational Neuroscience An Introductory Course in Computational Neuroscience Contents Series Foreword Acknowledgments Preface 1 Preliminary Material 1.1. Introduction 1.1.1 The Cell, the Circuit, and the Brain 1.1.2 Physics of

More information

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/319/5869/1543/dc1 Supporting Online Material for Synaptic Theory of Working Memory Gianluigi Mongillo, Omri Barak, Misha Tsodyks* *To whom correspondence should be addressed.

More information

Chapter 48 Neurons, Synapses, and Signaling

Chapter 48 Neurons, Synapses, and Signaling Chapter 48 Neurons, Synapses, and Signaling Concept 48.1 Neuron organization and structure reflect function in information transfer Neurons are nerve cells that transfer information within the body Neurons

More information

Phase Response Properties and Phase-Locking in Neural Systems with Delayed Negative-Feedback. Carter L. Johnson

Phase Response Properties and Phase-Locking in Neural Systems with Delayed Negative-Feedback. Carter L. Johnson Phase Response Properties and Phase-Locking in Neural Systems with Delayed Negative-Feedback Carter L. Johnson Faculty Mentor: Professor Timothy J. Lewis University of California, Davis Abstract Oscillatory

More information

Neural Modeling and Computational Neuroscience. Claudio Gallicchio

Neural 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 information

Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity

Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity Jorge F. Mejias 1,2 and Joaquín J. Torres 2 1 Department of Physics and Center for

More information

Chapter 9. Nerve Signals and Homeostasis

Chapter 9. Nerve Signals and Homeostasis Chapter 9 Nerve Signals and Homeostasis A neuron is a specialized nerve cell that is the functional unit of the nervous system. Neural signaling communication by neurons is the process by which an animal

More information

Neurons, Synapses, and Signaling

Neurons, Synapses, and Signaling LECTURE PRESENTATIONS For CAMPBELL BIOLOGY, NINTH EDITION Jane B. Reece, Lisa A. Urry, Michael L. Cain, Steven A. Wasserman, Peter V. Minorsky, Robert B. Jackson Chapter 48 Neurons, Synapses, and Signaling

More information

Introduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten. Lecture 2a. The Neuron - overview of structure. From Anderson (1995)

Introduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten. Lecture 2a. The Neuron - overview of structure. From Anderson (1995) Introduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten Lecture 2a The Neuron - overview of structure From Anderson (1995) 2 Lect_2a_Mathematica.nb Basic Structure Information flow:

More information

Dendrites - receives information from other neuron cells - input receivers.

Dendrites - receives information from other neuron cells - input receivers. The Nerve Tissue Neuron - the nerve cell Dendrites - receives information from other neuron cells - input receivers. Cell body - includes usual parts of the organelles of a cell (nucleus, mitochondria)

More information

BIOLOGY. 1. Overview of Neurons 11/3/2014. Neurons, Synapses, and Signaling. Communication in Neurons

BIOLOGY. 1. Overview of Neurons 11/3/2014. Neurons, Synapses, and Signaling. Communication in Neurons CAMPBELL BIOLOGY TENTH EDITION 48 Reece Urry Cain Wasserman Minorsky Jackson Neurons, Synapses, and Signaling Lecture Presentation by Nicole Tunbridge and Kathleen Fitzpatrick 1. Overview of Neurons Communication

More information

Topics in Neurophysics

Topics in Neurophysics Topics in Neurophysics Alex Loebel, Martin Stemmler and Anderas Herz Exercise 2 Solution (1) The Hodgkin Huxley Model The goal of this exercise is to simulate the action potential according to the model

More information

Nervous Tissue. Neurons Electrochemical Gradient Propagation & Transduction Neurotransmitters Temporal & Spatial Summation

Nervous Tissue. Neurons Electrochemical Gradient Propagation & Transduction Neurotransmitters Temporal & Spatial Summation Nervous Tissue Neurons Electrochemical Gradient Propagation & Transduction Neurotransmitters Temporal & Spatial Summation What is the function of nervous tissue? Maintain homeostasis & respond to stimuli

More information

Information processing. Divisions of nervous system. Neuron structure and function Synapse. Neurons, synapses, and signaling 11/3/2017

Information processing. Divisions of nervous system. Neuron structure and function Synapse. Neurons, synapses, and signaling 11/3/2017 Neurons, synapses, and signaling Chapter 48 Information processing Divisions of nervous system Central nervous system (CNS) Brain and a nerve cord Integration center Peripheral nervous system (PNS) Nerves

More information

Hopfield Neural Network and Associative Memory. Typical Myelinated Vertebrate Motoneuron (Wikipedia) Topic 3 Polymers and Neurons Lecture 5

Hopfield Neural Network and Associative Memory. Typical Myelinated Vertebrate Motoneuron (Wikipedia) Topic 3 Polymers and Neurons Lecture 5 Hopfield Neural Network and Associative Memory Typical Myelinated Vertebrate Motoneuron (Wikipedia) PHY 411-506 Computational Physics 2 1 Wednesday, March 5 1906 Nobel Prize in Physiology or Medicine.

More information

Reducing neuronal networks to discrete dynamics

Reducing neuronal networks to discrete dynamics Physica D 237 (2008) 324 338 www.elsevier.com/locate/physd Reducing neuronal networks to discrete dynamics David Terman a,b,, Sungwoo Ahn a, Xueying Wang a, Winfried Just c a Department of Mathematics,

More information

Nervous Systems: Neuron Structure and Function

Nervous Systems: Neuron Structure and Function Nervous Systems: Neuron Structure and Function Integration An animal needs to function like a coherent organism, not like a loose collection of cells. Integration = refers to processes such as summation

More information

Neurons, Synapses, and Signaling

Neurons, Synapses, and Signaling Chapter 48 Neurons, Synapses, and Signaling PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions

More information

Tuning tuning curves. So far: Receptive fields Representation of stimuli Population vectors. Today: Contrast enhancment, cortical processing

Tuning tuning curves. So far: Receptive fields Representation of stimuli Population vectors. Today: Contrast enhancment, cortical processing Tuning tuning curves So far: Receptive fields Representation of stimuli Population vectors Today: Contrast enhancment, cortical processing Firing frequency N 3 s max (N 1 ) = 40 o N4 N 1 N N 5 2 s max

More information

MEMBRANE POTENTIALS AND ACTION POTENTIALS:

MEMBRANE POTENTIALS AND ACTION POTENTIALS: University of Jordan Faculty of Medicine Department of Physiology & Biochemistry Medical students, 2017/2018 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Review: Membrane physiology

More information

Information Theory and Neuroscience II

Information Theory and Neuroscience II John Z. Sun and Da Wang Massachusetts Institute of Technology October 14, 2009 Outline System Model & Problem Formulation Information Rate Analysis Recap 2 / 23 Neurons Neuron (denoted by j) I/O: via synapses

More information

CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE. by Pamela Reitsma. B.S., University of Maine, 2007

CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE. by Pamela Reitsma. B.S., University of Maine, 2007 CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE by Pamela Reitsma B.S., University of Maine, 27 Submitted to the Graduate Faculty of the Department of Mathematics in partial

More information

Dynamical Constraints on Computing with Spike Timing in the Cortex

Dynamical Constraints on Computing with Spike Timing in the Cortex Appears in Advances in Neural Information Processing Systems, 15 (NIPS 00) Dynamical Constraints on Computing with Spike Timing in the Cortex Arunava Banerjee and Alexandre Pouget Department of Brain and

More information

BIOLOGY 11/10/2016. Neurons, Synapses, and Signaling. Concept 48.1: Neuron organization and structure reflect function in information transfer

BIOLOGY 11/10/2016. Neurons, Synapses, and Signaling. Concept 48.1: Neuron organization and structure reflect function in information transfer 48 Neurons, Synapses, and Signaling CAMPBELL BIOLOGY TENTH EDITION Reece Urry Cain Wasserman Minorsky Jackson Lecture Presentation by Nicole Tunbridge and Kathleen Fitzpatrick Concept 48.1: Neuron organization

More information

Control and Integration. Nervous System Organization: Bilateral Symmetric Animals. Nervous System Organization: Radial Symmetric Animals

Control and Integration. Nervous System Organization: Bilateral Symmetric Animals. Nervous System Organization: Radial Symmetric Animals Control and Integration Neurophysiology Chapters 10-12 Nervous system composed of nervous tissue cells designed to conduct electrical impulses rapid communication to specific cells or groups of cells Endocrine

More information

Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting

Dynamical 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 information

Synchronization in two neurons: Results for a two-component dynamical model with time-delayed inhibition

Synchronization in two neurons: Results for a two-component dynamical model with time-delayed inhibition Physica D 4 (998 47 7 Synchronization in two neurons: Results for a two-component dynamical model with time-delayed inhibition Gilles Renversez Centre de Physique Théorique, UPR 76 du Centre National de

More information

2 1. Introduction. Neuronal networks often exhibit a rich variety of oscillatory behavior. The dynamics of even a single cell may be quite complicated

2 1. Introduction. Neuronal networks often exhibit a rich variety of oscillatory behavior. The dynamics of even a single cell may be quite complicated GEOMETRIC ANALYSIS OF POPULATION RHYTHMS IN SYNAPTICALLY COUPLED NEURONAL NETWORKS J. Rubin and D. Terman Dept. of Mathematics; Ohio State University; Columbus, Ohio 43210 Abstract We develop geometric

More information

Compartmental Modelling

Compartmental Modelling Modelling Neurons Computing and the Brain Compartmental Modelling Spring 2010 2 1 Equivalent Electrical Circuit A patch of membrane is equivalent to an electrical circuit This circuit can be described

More information

Synaptic dynamics. John D. Murray. Synaptic currents. Simple model of the synaptic gating variable. First-order kinetics

Synaptic dynamics. John D. Murray. Synaptic currents. Simple model of the synaptic gating variable. First-order kinetics Synaptic dynamics John D. Murray A dynamical model for synaptic gating variables is presented. We use this to study the saturation of synaptic gating at high firing rate. Shunting inhibition and the voltage

More information

NON-MARKOVIAN SPIKING STATISTICS OF A NEURON WITH DELAYED FEEDBACK IN PRESENCE OF REFRACTORINESS. Kseniia Kravchuk and Alexander Vidybida

NON-MARKOVIAN SPIKING STATISTICS OF A NEURON WITH DELAYED FEEDBACK IN PRESENCE OF REFRACTORINESS. Kseniia Kravchuk and Alexander Vidybida MATHEMATICAL BIOSCIENCES AND ENGINEERING Volume 11, Number 1, February 214 doi:1.3934/mbe.214.11.xx pp. 1 xx NON-MARKOVIAN SPIKING STATISTICS OF A NEURON WITH DELAYED FEEDBACK IN PRESENCE OF REFRACTORINESS

More information

DISCRETE EVENT SIMULATION IN THE NEURON ENVIRONMENT

DISCRETE EVENT SIMULATION IN THE NEURON ENVIRONMENT Hines and Carnevale: Discrete event simulation in the NEURON environment Page 1 Preprint of a manuscript that will be published in Neurocomputing. DISCRETE EVENT SIMULATION IN THE NEURON ENVIRONMENT Abstract

More information

Modelling stochastic neural learning

Modelling stochastic neural learning Modelling stochastic neural learning Computational Neuroscience András Telcs telcs.andras@wigner.mta.hu www.cs.bme.hu/~telcs http://pattern.wigner.mta.hu/participants/andras-telcs Compiled from lectures

More information

Neural Conduction. biologyaspoetry.com

Neural Conduction. biologyaspoetry.com Neural Conduction biologyaspoetry.com Resting Membrane Potential -70mV A cell s membrane potential is the difference in the electrical potential ( charge) between the inside and outside of the cell. The

More information

Lecture 11 : Simple Neuron Models. Dr Eileen Nugent

Lecture 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 information

Localized activity patterns in excitatory neuronal networks

Localized activity patterns in excitatory neuronal networks Localized activity patterns in excitatory neuronal networks Jonathan Rubin Amitabha Bose February 3, 2004 Abstract. The existence of localized activity patterns, or bumps, has been investigated in a variety

More information

Ch 7. The Nervous System 7.1 & 7.2

Ch 7. The Nervous System 7.1 & 7.2 Ch 7 The Nervous System 7.1 & 7.2 SLOs Describe the different types of neurons and supporting cells, and identify their functions. Identify the myelin sheath and describe how it is formed in the CNS and

More information

9 Generation of Action Potential Hodgkin-Huxley Model

9 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 information

Neurons, Synapses, and Signaling

Neurons, Synapses, and Signaling CAMPBELL BIOLOGY IN FOCUS URRY CAIN WASSERMAN MINORSKY REECE 37 Neurons, Synapses, and Signaling Lecture Presentations by Kathleen Fitzpatrick and Nicole Tunbridge, Simon Fraser University SECOND EDITION

More information

DEVS Simulation of Spiking Neural Networks

DEVS Simulation of Spiking Neural Networks DEVS Simulation of Spiking Neural Networks Rene Mayrhofer, Michael Affenzeller, Herbert Prähofer, Gerhard Höfer, Alexander Fried Institute of Systems Science Systems Theory and Information Technology Johannes

More information

arxiv:physics/ v1 [physics.bio-ph] 19 Feb 1999

arxiv:physics/ v1 [physics.bio-ph] 19 Feb 1999 Odor recognition and segmentation by coupled olfactory bulb and cortical networks arxiv:physics/9902052v1 [physics.bioph] 19 Feb 1999 Abstract Zhaoping Li a,1 John Hertz b a CBCL, MIT, Cambridge MA 02139

More information

CSE/NB 528 Final Lecture: All Good Things Must. CSE/NB 528: Final Lecture

CSE/NB 528 Final Lecture: All Good Things Must. CSE/NB 528: Final Lecture CSE/NB 528 Final Lecture: All Good Things Must 1 Course Summary Where have we been? Course Highlights Where do we go from here? Challenges and Open Problems Further Reading 2 What is the neural code? What

More information

Executable Symbolic Modeling of Neural Processes

Executable Symbolic Modeling of Neural Processes Executable Symbolic Modeling of Neural Processes M Sriram Iyengar 1, Carolyn Talcott 2, Riccardo Mozzachiodi 3, Douglas Baxter 3 1 School of Health Information Sciences, Univ. of Texas Health Science Center

More information

Chapter 37 Active Reading Guide Neurons, Synapses, and Signaling

Chapter 37 Active Reading Guide Neurons, Synapses, and Signaling Name: AP Biology Mr. Croft Section 1 1. What is a neuron? Chapter 37 Active Reading Guide Neurons, Synapses, and Signaling 2. Neurons can be placed into three groups, based on their location and function.

More information

FRTF01 L8 Electrophysiology

FRTF01 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 information

Conductance-Based Integrate-and-Fire Models

Conductance-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 information

Neurons, Synapses, and Signaling

Neurons, Synapses, and Signaling Chapter 48 Neurons, Synapses, and Signaling PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions

More information

Neural Networks 1 Synchronization in Spiking Neural Networks

Neural Networks 1 Synchronization in Spiking Neural Networks CS 790R Seminar Modeling & Simulation Neural Networks 1 Synchronization in Spiking Neural Networks René Doursat Department of Computer Science & Engineering University of Nevada, Reno Spring 2006 Synchronization

More information

Dynamical Systems in Neuroscience: Elementary Bifurcations

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 information

THE LOCUST OLFACTORY SYSTEM AS A CASE STUDY FOR MODELING DYNAMICS OF NEUROBIOLOGICAL NETWORKS: FROM DISCRETE TIME NEURONS TO CONTINUOUS TIME NEURONS

THE LOCUST OLFACTORY SYSTEM AS A CASE STUDY FOR MODELING DYNAMICS OF NEUROBIOLOGICAL NETWORKS: FROM DISCRETE TIME NEURONS TO CONTINUOUS TIME NEURONS 1 THE LOCUST OLFACTORY SYSTEM AS A CASE STUDY FOR MODELING DYNAMICS OF NEUROBIOLOGICAL NETWORKS: FROM DISCRETE TIME NEURONS TO CONTINUOUS TIME NEURONS B. QUENET 1 AND G. HORCHOLLE-BOSSAVIT 2 1 Equipe de

More information

How do synapses transform inputs?

How do synapses transform inputs? Neurons to networks How do synapses transform inputs? Excitatory synapse Input spike! Neurotransmitter release binds to/opens Na channels Change in synaptic conductance! Na+ influx E.g. AMA synapse! Depolarization

More information

Neurons, Synapses, and Signaling

Neurons, Synapses, and Signaling Chapter 48 Neurons, Synapses, and Signaling PowerPoint Lectures for Biology, Eighth Edition Lectures by Chris Romero, updated by Erin Barley with contributions from Joan Sharp and Janette Lewis Copyright

More information

Nerve Signal Conduction. Resting Potential Action Potential Conduction of Action Potentials

Nerve Signal Conduction. Resting Potential Action Potential Conduction of Action Potentials Nerve Signal Conduction Resting Potential Action Potential Conduction of Action Potentials Resting Potential Resting neurons are always prepared to send a nerve signal. Neuron possesses potential energy

More information

Phase Response. 1 of of 11. Synaptic input advances (excitatory) or delays (inhibitory) spiking

Phase Response. 1 of of 11. Synaptic input advances (excitatory) or delays (inhibitory) spiking Printed from the Mathematica Help Browser 1 1 of 11 Phase Response Inward current-pulses decrease a cortical neuron's period (Cat, Layer V). [Fetz93] Synaptic input advances (excitatory) or delays (inhibitory)

More information

A Three-dimensional Physiologically Realistic Model of the Retina

A Three-dimensional Physiologically Realistic Model of the Retina A Three-dimensional Physiologically Realistic Model of the Retina Michael Tadross, Cameron Whitehouse, Melissa Hornstein, Vicky Eng and Evangelia Micheli-Tzanakou Department of Biomedical Engineering 617

More information

Probabilistic Models in Theoretical Neuroscience

Probabilistic Models in Theoretical Neuroscience Probabilistic Models in Theoretical Neuroscience visible unit Boltzmann machine semi-restricted Boltzmann machine restricted Boltzmann machine hidden unit Neural models of probabilistic sampling: introduction

More information

On oscillating systems of interacting neurons

On oscillating systems of interacting neurons Susanne Ditlevsen Eva Löcherbach Nice, September 2016 Oscillations Most biological systems are characterized by an intrinsic periodicity : 24h for the day-night rhythm of most animals, etc... In other

More information

Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System

Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Dimitra Despoina Pagania 1, Adam Adamopoulos 1,2 and Spiridon D. Likothanassis 1 1 Pattern Recognition Laboratory,

More information

The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs

The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs Jonathan Rubin Department of Mathematics University of Pittsburgh Pittsburgh, PA 15260 USA Krešimir Josić

More information

Electrophysiology of the neuron

Electrophysiology 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 information

Approximate, not perfect synchrony maximizes the downstream effectiveness of excitatory neuronal ensembles

Approximate, not perfect synchrony maximizes the downstream effectiveness of excitatory neuronal ensembles Börgers et al. RESEARCH Approximate, not perfect synchrony maximizes the downstream effectiveness of excitatory neuronal ensembles Christoph Börgers *, Jie Li and Nancy Kopell 2 * Correspondence: cborgers@tufts.edu

More information

Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued

Mathematical 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 information

Neurons and Nervous Systems

Neurons and Nervous Systems 34 Neurons and Nervous Systems Concept 34.1 Nervous Systems Consist of Neurons and Glia Nervous systems have two categories of cells: Neurons, or nerve cells, are excitable they generate and transmit electrical

More information

Coupling in Networks of Neuronal Oscillators. Carter Johnson

Coupling in Networks of Neuronal Oscillators. Carter Johnson Coupling in Networks of Neuronal Oscillators Carter Johnson June 15, 2015 1 Introduction Oscillators are ubiquitous in nature. From the pacemaker cells that keep our hearts beating to the predator-prey

More information

Temporal Pattern Analysis

Temporal Pattern Analysis LIACS Leiden Institute of Advanced Computer Science Master s Thesis June 17, 29 Temporal Pattern Analysis Using Reservoir Computing Author: Ron Vink Supervisor: Dr. Walter Kosters 1 Contents 1 Introduction

More information

Structure and Measurement of the brain lecture notes

Structure and Measurement of the brain lecture notes Structure and Measurement of the brain lecture notes Marty Sereno 2009/2010!"#$%&'(&#)*%$#&+,'-&.)"/*"&.*)*-'(0&1223 Neurons and Models Lecture 1 Topics Membrane (Nernst) Potential Action potential/voltage-gated

More information

Neurophysiology. Danil Hammoudi.MD

Neurophysiology. Danil Hammoudi.MD Neurophysiology Danil Hammoudi.MD ACTION POTENTIAL An action potential is a wave of electrical discharge that travels along the membrane of a cell. Action potentials are an essential feature of animal

More information

CHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING

CHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING CHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING 4.1 Introduction Thus far we have studied synchronization in networks of locally coupled neurobiologically based oscillator models. One additional

More information

Nervous Tissue. Neurons Neural communication Nervous Systems

Nervous Tissue. Neurons Neural communication Nervous Systems Nervous Tissue Neurons Neural communication Nervous Systems What is the function of nervous tissue? Maintain homeostasis & respond to stimuli Sense & transmit information rapidly, to specific cells and

More information

Liquid Computing in a Simplified Model of Cortical Layer IV: Learning to Balance a Ball

Liquid Computing in a Simplified Model of Cortical Layer IV: Learning to Balance a Ball Liquid Computing in a Simplified Model of Cortical Layer IV: Learning to Balance a Ball Dimitri Probst 1,3, Wolfgang Maass 2, Henry Markram 1, and Marc-Oliver Gewaltig 1 1 Blue Brain Project, École Polytechnique

More information

Bio 449 Fall Exam points total Multiple choice. As with any test, choose the best answer in each case. Each question is 3 points.

Bio 449 Fall Exam points total Multiple choice. As with any test, choose the best answer in each case. Each question is 3 points. Name: Exam 1 100 points total Multiple choice. As with any test, choose the best answer in each case. Each question is 3 points. 1. The term internal environment, as coined by Clause Bernard, is best defined

More information

1. Neurons & Action Potentials

1. Neurons & Action Potentials Lecture 6, 30 Jan 2008 Vertebrate Physiology ECOL 437 (MCB/VetSci 437) Univ. of Arizona, spring 2008 Kevin Bonine & Kevin Oh 1. Intro Nervous System Fxn (slides 32-60 from Mon 28 Jan; Ch10) 2. Neurons

More information

37 Neurons, Synapses, and Signaling

37 Neurons, Synapses, and Signaling CAMPBELL BIOLOGY IN FOCUS Urry Cain Wasserman Minorsky Jackson Reece 37 Neurons, Synapses, and Signaling Lecture Presentations by Kathleen Fitzpatrick and Nicole Tunbridge Overview: Lines of Communication

More information

Roles of Fluctuations. in Pulsed Neural Networks

Roles of Fluctuations. in Pulsed Neural Networks Ph.D. Thesis Roles of Fluctuations in Pulsed Neural Networks Supervisor Professor Yoichi Okabe Takashi Kanamaru Department of Advanced Interdisciplinary Studies, Faculty of Engineering, The University

More information

Infinite systems of interacting chains with memory of variable length - a stochastic model for biological neural nets

Infinite systems of interacting chains with memory of variable length - a stochastic model for biological neural nets Infinite systems of interacting chains with memory of variable length - a stochastic model for biological neural nets Antonio Galves Universidade de S.Paulo Fapesp Center for Neuromathematics Eurandom,

More information

9 Generation of Action Potential Hodgkin-Huxley Model

9 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 information

Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System

Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Dimitra-Despoina Pagania,*, Adam Adamopoulos,2, and Spiridon D. Likothanassis Pattern Recognition Laboratory,

More information

Dendritic computation

Dendritic computation Dendritic computation Dendrites as computational elements: Passive contributions to computation Active contributions to computation Examples Geometry matters: the isopotential cell Injecting current I

More information

IN THIS turorial paper we exploit the relationship between

IN THIS turorial paper we exploit the relationship between 508 IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 3, MAY 1999 Weakly Pulse-Coupled Oscillators, FM Interactions, Synchronization, Oscillatory Associative Memory Eugene M. Izhikevich Abstract We study

More information

The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs

The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs The firing of an excitable neuron in the presence of stochastic trains of strong synaptic inputs Jonathan Rubin Department of Mathematics University of Pittsburgh Pittsburgh, PA 15260 USA Krešimir Josić

More information

Voltage-clamp and Hodgkin-Huxley models

Voltage-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 information

Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback

Neural 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 information

Approximate, not Perfect Synchrony Maximizes the Downstream Effectiveness of Excitatory Neuronal Ensembles

Approximate, not Perfect Synchrony Maximizes the Downstream Effectiveness of Excitatory Neuronal Ensembles Journal of Mathematical Neuroscience (2014) 4:10 DOI 10.1186/2190-8567-4-10 RESEARCH OpenAccess Approximate, not Perfect Synchrony Maximizes the Downstream Effectiveness of Excitatory Neuronal Ensembles

More information

Neurons. The Molecular Basis of their Electrical Excitability

Neurons. The Molecular Basis of their Electrical Excitability Neurons The Molecular Basis of their Electrical Excitability Viva La Complexity! Consider, The human brain contains >10 11 neurons! Each neuron makes 10 3 (average) synaptic contacts on up to 10 3 other

More information

Artificial Neural Network and Fuzzy Logic

Artificial Neural Network and Fuzzy Logic Artificial Neural Network and Fuzzy Logic 1 Syllabus 2 Syllabus 3 Books 1. Artificial Neural Networks by B. Yagnanarayan, PHI - (Cover Topologies part of unit 1 and All part of Unit 2) 2. Neural Networks

More information

Dynamic Modeling of Brain Activity

Dynamic Modeling of Brain Activity 0a Dynamic Modeling of Brain Activity EIN IIN PC Thomas R. Knösche, Leipzig Generative Models for M/EEG 4a Generative Models for M/EEG states x (e.g. dipole strengths) parameters parameters (source positions,

More information

Fast and exact simulation methods applied on a broad range of neuron models

Fast and exact simulation methods applied on a broad range of neuron models Fast and exact simulation methods applied on a broad range of neuron models Michiel D Haene michiel.dhaene@ugent.be Benjamin Schrauwen benjamin.schrauwen@ugent.be Ghent University, Electronics and Information

More information

Nervous System: Part II How A Neuron Works

Nervous System: Part II How A Neuron Works Nervous System: Part II How A Neuron Works Essential Knowledge Statement 3.E.2 Continued Animals have nervous systems that detect external and internal signals, transmit and integrate information, and

More information

Overview Organization: Central Nervous System (CNS) Peripheral Nervous System (PNS) innervate Divisions: a. Afferent

Overview Organization: Central Nervous System (CNS) Peripheral Nervous System (PNS) innervate Divisions: a. Afferent Overview Organization: Central Nervous System (CNS) Brain and spinal cord receives and processes information. Peripheral Nervous System (PNS) Nerve cells that link CNS with organs throughout the body.

More information

Memories Associated with Single Neurons and Proximity Matrices

Memories Associated with Single Neurons and Proximity Matrices Memories Associated with Single Neurons and Proximity Matrices Subhash Kak Oklahoma State University, Stillwater Abstract: This paper extends the treatment of single-neuron memories obtained by the use

More information

Introduction and the Hodgkin-Huxley Model

Introduction 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 information

Neurophysiology of a VLSI spiking neural network: LANN21

Neurophysiology of a VLSI spiking neural network: LANN21 Neurophysiology of a VLSI spiking neural network: LANN21 Stefano Fusi INFN, Sezione Roma I Università di Roma La Sapienza Pza Aldo Moro 2, I-185, Roma fusi@jupiter.roma1.infn.it Paolo Del Giudice Physics

More information

LESSON 2.2 WORKBOOK How do our axons transmit electrical signals?

LESSON 2.2 WORKBOOK How do our axons transmit electrical signals? LESSON 2.2 WORKBOOK How do our axons transmit electrical signals? This lesson introduces you to the action potential, which is the process by which axons signal electrically. In this lesson you will learn

More information

Balance of Electric and Diffusion Forces

Balance of Electric and Diffusion Forces Balance of Electric and Diffusion Forces Ions flow into and out of the neuron under the forces of electricity and concentration gradients (diffusion). The net result is a electric potential difference

More information

Computing with Inter-spike Interval Codes in Networks of Integrate and Fire Neurons

Computing with Inter-spike Interval Codes in Networks of Integrate and Fire Neurons Computing with Inter-spike Interval Codes in Networks of Integrate and Fire Neurons Dileep George a,b Friedrich T. Sommer b a Dept. of Electrical Engineering, Stanford University 350 Serra Mall, Stanford,

More information

Modeling Action Potentials in Cell Processes

Modeling 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 information

Biological Modeling of Neural Networks:

Biological Modeling of Neural Networks: Week 14 Dynamics and Plasticity 14.1 Reservoir computing - Review:Random Networks - Computing with rich dynamics Biological Modeling of Neural Networks: 14.2 Random Networks - stationary state - chaos

More information

The Spike Response Model: A Framework to Predict Neuronal Spike Trains

The Spike Response Model: A Framework to Predict Neuronal Spike Trains The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology

More information

An algorithm for detecting oscillatory behavior in discretized data: the damped-oscillator oscillator detector

An algorithm for detecting oscillatory behavior in discretized data: the damped-oscillator oscillator detector An algorithm for detecting oscillatory behavior in discretized data: the damped-oscillator oscillator detector David Hsu, Murielle Hsu, He Huang and Erwin B. Montgomery, Jr Department of Neurology University

More information