Simple models of neurons!! Lecture 4!

Size: px
Start display at page:

Download "Simple models of neurons!! Lecture 4!"

Transcription

1 Simple models of neurons!! Lecture 4! 1!

2 Recap: Phase Plane Analysis! 2!

3 FitzHugh-Nagumo Model! Membrane potential K activation variable Notes from Zillmer, INFN 3!

4 FitzHugh Nagumo Model Demo! 4!

5 Phase Plane Analysis - Summary! Scholarpedia, FitzHugh Nagumo Model 5!

6 Recap HH model! 6!

7 Reducing HH models(2)! Original Reduced 7!

8 Reducing HH models(3)! 8! Dynamical Systems Neuroscience, Izhikevich

9 Reducing HH to a 2-D equation! g For all three variables! g First approximation: Replace m by its asymptoic value! Abbott and Kepler, Model Neurons: From HH to Hopfield 9!

10 Reducing HH to a 2-D equation! g 2-d model! g We want the time-dependence of U in f in reduced model to approximate the timedependence of F in the full model by changing h and n! 10!

11 Two-dimensional version of HH! 11!

12 0 Two-dimensional version of HH! 30 Limit cycle due to constant current injection (train of action potentials) 40 U V - nullcline V 20 12!

13 General HH Model Reduction Strategy! g Start with one equation for V and one for recovery variable (lets call it R)! g To match the models in phase space, the first equation has to be a cubic polynomial in V! [Wilson, Spikes, decisions, and actions, 1998]! 13!

14 General HH Model Reduction Strategy! g Start with one equation for V and one for recovery variable (lets call it R)!! g Null clines of this system of equations are:! [Wilson, Spikes, decisions, and actions, 1999]! 14!

15 General HH Model Reduction Strategy! Direct! reduction! Phase plane! reduction (eqn below)!! g Fitting to nullclines we get:! [Wilson, Spikes, decisions, and actions, 1999]! 15!

16 General HH Model Reduction Strategy! [Wilson, Spikes, decisions, and actions, 1999]! 16!

17 What is the code here?! g Intracellular recording from a locust projection neuron to 1s odor puffs! 60 mv 17!

18 One-Dimensional Reductions! g Perfect Integrate and Fire Model! I(t) C V $"(t # t i ) Whatʼs missing! In this model! compared to HH?! t ref dv (t) C dt V t = I(t) ( ) = V Thr " Fire+reset linear threshold 18!

19 One-dimensional Reductions! g Perfect Integrate and Fire Model! I(t) V $"(t # t i ) V Spike emission V Thr C reset t ref dv (t) C dt V t = I(t) ( ) = V Thr " Fire+reset linear threshold 19!

20 One dimensional reductions! g The successive times, t i, of spike occurrence:!! t i+1 " t i I(t)dt = CV th g Firing rate vs. input current of the perfect integrator:! I f = CV! Thr g If you force a refractory period T ref following a spike, such that V = 0mV for T ref period following a spike, then:! I f = CV Thr + t ref I 20!

21 One-dimensional Reductions! g Leaky Integrate and Fire Model! I(t) V $"(t # t i ) V Spike emission V Thr C R reset t ref dv (t) C dt V t + V (t) R = I(t) ( ) = V Thr " Fire+reset linear threshold 21!

22 Comparison of the two models! Perfect Integrate & Fire Leaky Integrate & Fire Model Input 7 V(mV) C = 1nF R = 10MOhm VSPK = 70mV VTHR = 15 mv I = 1 na What is the main! difference?! Time (ms) 22!

23 C I(t) R One-dimensional Reductions! g Adapting Leaky Integrate and Fire Model! V $"(t # t i ) g adapt t ref C dv (t) dt " adapt dg adapt (t) dt + V (t) R + g adaptv (t) = I(t) linear = #g adapt (t) adaptation V ( t) = V Thr " Fire+reset threshold g adapt (t) = g adapt (t) + G inc 23!

24 Comparison of the three models! Perfect Integrate & Fire Leaky Integrate & Fire Input Adapting Integrate & Fire C = 1nF R = 10MOhm VSPK = 70mV VTHR = 15 mv τ G = 50 ms G inc = 0.2 ns I = 1nA !

25 Comparison of the three models! Perfect I&F Leaky I&F Adapt I&F Firing rate (hertz) Input Amplitude (na) 25!

26 Low pass filter! 26!

27 One Dimensional Models! g A more principled approach (again based on 2-d phase plane dynamics!!!)! Input current moves! V nullcline! and! therefore alters stability of the! equilibrium point! y Slow Variable (W) FitzHugh-Nagumo Model W-nullcline meets! left or right segments of! V-nullcline then the! equilibrium point is stable! W-nullcline meets! Center segments of! V-nullcline then the! equilibrium point is unstable! x Fast Variable (V) 27!

28 One Dimensional Models! g A more principled approach (again based on 2-d phase plane dynamics!!!)! Input current moves! V nullcline! and! therefore alters stability of the! equilibrium point! y Slow Variable (W) Threshold behavior once the current alters! the stability of the fixed-point (causes spike)! FitzHugh-Nagumo Model W-nullcline meets! left or right segments of! V-nullcline then the! equilibrium point is stable! W-nullcline meets! Center segments of! V-nullcline then the! equilibrium point is unstable! x Fast Variable (V) 28!

29 One Dimensional Models! g A more principled approach (again based on 2-d phase plane dynamics!!!)! Sub-threshold dynamics captured by this highlighted region If I want to construct a integrate and fire model only this region is important!!! 29! Dynamical Systems Neuroscience, Izhikevich

30 Quadratic Integrate and Fire Model! g A more principled approach (again based on 2-d phase plane dynamics!!!)! Notice: in this highlighted region V-nullcline is a parbola U-nullcline is still a line Sub-threshold dynamics captured by this highlighted region 30! Dynamical Systems Neuroscience, Izhikevich

31 Izhikevich Model! g A simple model that captures the subthreshold behavior in a small neighborhood of the left knee (confined to the shaded square) and the initial segment of the upstroke of an action potential is given by:! where a, b, c, d are dimensionless parameters! 31!

32 Izhikevich Model! 32!

33 Izhikevich Model! 33!

34 I(t) Firing rate models! g The potential in a continuous firing rate unit has same dynamics as in a LIF neuron! V f= g(v) C R C dv (t) + V (t) dt R f = g(v ) = I(t) Sigmoidal function 34!

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

Chapter 24 BIFURCATIONS

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

Nonlinear Dynamics of Neural Firing

Nonlinear 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 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

Computational Neuroscience. Session 4-2

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

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

6.3.4 Action potential

6.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 information

Single-Cell and Mean Field Neural Models

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

3 Action Potentials - Brutal Approximations

3 Action Potentials - Brutal Approximations Physics 172/278 - David Kleinfeld - Fall 2004; Revised Winter 2015 3 Action Potentials - Brutal Approximations The Hodgkin-Huxley equations for the behavior of the action potential in squid, and similar

More information

Lecture IV: LTI models of physical systems

Lecture IV: LTI models of physical systems BME 171: Signals and Systems Duke University September 5, 2008 This lecture Plan for the lecture: 1 Interconnections of linear systems 2 Differential equation models of LTI systems 3 eview of linear circuit

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

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

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

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

Presented by Sarah Hedayat. Supervised by Pr.Cappy and Dr.Hoel

Presented by Sarah Hedayat. Supervised by Pr.Cappy and Dr.Hoel 1 Presented by Sarah Hedayat Supervised by Pr.Cappy and Dr.Hoel Outline 2 Project objectives Key elements Membrane models As simple as possible Phase plane analysis Review of important Concepts Conclusion

More information

Dynamical Systems for Biology - Excitability

Dynamical Systems for Biology - Excitability Dynamical Systems for Biology - Excitability J. P. Keener Mathematics Department Dynamical Systems for Biology p.1/25 Examples of Excitable Media B-Z reagent CICR (Calcium Induced Calcium Release) Nerve

More information

Neuronal Dynamics: Computational Neuroscience of Single Neurons

Neuronal Dynamics: Computational Neuroscience of Single Neurons Nonlinear Integrate-and-Fire Model Neronal Dynamics: Comptational Neroscience of Single Nerons Nonlinear Integrate-and-fire (NLIF) - Definition - qadratic and expon. IF - Extracting NLIF model from data

More information

NE 204 mini-syllabus (weeks 4 8)

NE 204 mini-syllabus (weeks 4 8) NE 24 mini-syllabus (weeks 4 8) Instructor: John Burke O ce: MCS 238 e-mail: jb@math.bu.edu o ce hours: by appointment Overview: For the next few weeks, we will focus on mathematical models of single neurons.

More information

Single neuron models. L. Pezard Aix-Marseille University

Single 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 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

1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model

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

Neuronal Dynamics: Computational Neuroscience of Single Neurons

Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 4 part 5: Nonlinear Integrate-and-Fire Model 4.1 From Hodgkin-Huxley to 2D Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 4 Recing detail: Two-dimensional neuron models Wulfram

More information

9.01 Introduction to Neuroscience Fall 2007

9.01 Introduction to Neuroscience Fall 2007 MIT OpenCourseWare http://ocw.mit.edu 9.01 Introduction to Neuroscience Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 9.01 Recitation (R02)

More information

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wires -- signal propagation -- processing

More information

Mathematical Foundations of Neuroscience - Lecture 9. Simple Models of Neurons and Synapses.

Mathematical Foundations of Neuroscience - Lecture 9. Simple Models of Neurons and Synapses. Mathematical Foundations of Neuroscience - Lecture 9. Simple Models of and. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010

More information

AT2 Neuromodeling: Problem set #3 SPIKE TRAINS

AT2 Neuromodeling: Problem set #3 SPIKE TRAINS AT2 Neuromodeling: Problem set #3 SPIKE TRAINS Younesse Kaddar PROBLEM 1: Poisson spike trains Link of the ipython notebook for the code Brain neuron emit spikes seemingly randomly: we will aim to model

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

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

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 Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913

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

Single-Compartment Neural Models

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

We observe the model neuron s response to constant input current, studying the dependence of:

We observe the model neuron s response to constant input current, studying the dependence of: BioE332A Lab 2, 21 1 Lab 2 December 23, 29 A Spiking Neuron Like biological neurons, the model neuron we characterize in this lab has a repertoire of (model) ion-channel populations. Here we focus on the

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

Chapter 2 The Hodgkin Huxley Theory of Neuronal Excitation

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

LIMIT CYCLE OSCILLATORS

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

REBECCA J GAFF (Under the direction of Andrew Sornborger) ABSTRACT. Imaging research of seizures in Zebrafish brains involves tracing calcium levels

REBECCA J GAFF (Under the direction of Andrew Sornborger) ABSTRACT. Imaging research of seizures in Zebrafish brains involves tracing calcium levels MODELING SEIZURES WITH A CONDUCTANCE-BASED INTEGRATE-AND-FIRE NEURONAL NETWORK MODEL by REBECCA J GAFF (Under the direction of Andrew Sornborger) ABSTRACT Imaging research of seizures in Zebrafish brains

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

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

Biomedical Instrumentation

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

Decoding. How well can we learn what the stimulus is by looking at the neural responses?

Decoding. How well can we learn what the stimulus is by looking at the neural responses? Decoding How well can we learn what the stimulus is by looking at the neural responses? Two approaches: devise explicit algorithms for extracting a stimulus estimate directly quantify the relationship

More information

Biological Modeling of Neural Networks

Biological Modeling of Neural Networks Week 3 part 1 : Rection of the Hodgkin-Huxley Model 3.1 From Hodgkin-Huxley to 2D Biological Modeling of Neural Netorks - Overvie: From 4 to 2 dimensions - MathDetour 1: Exploiting similarities - MathDetour

More information

Data-true Characterization Of Neuronal Models

Data-true Characterization Of Neuronal Models University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Data-true Characterization Of Neuronal Models 211 Jose Suarez University of Central Florida Find similar works

More information

Computational Explorations in Cognitive Neuroscience Chapter 2

Computational Explorations in Cognitive Neuroscience Chapter 2 Computational Explorations in Cognitive Neuroscience Chapter 2 2.4 The Electrophysiology of the Neuron Some basic principles of electricity are useful for understanding the function of neurons. This is

More information

arxiv: v2 [q-bio.nc] 13 Jan 2016

arxiv: v2 [q-bio.nc] 13 Jan 2016 Spike-Threshold Variation Originated from Separatrix-Crossing in Neuronal Dynamics Longfei Wang 1, Hengtong Wang 2, Lianchun Yu 1, and Yong Chen 3,4,* arxiv:159.2781v2 [q-bio.nc] 13 Jan 216 1 Institute

More information

Frequency Adaptation and Bursting

Frequency Adaptation and Bursting BioE332A Lab 3, 2010 1 Lab 3 January 5, 2010 Frequency Adaptation and Bursting In the last lab, we explored spiking due to sodium channels. In this lab, we explore adaptation and bursting due to potassium

More information

MATH 3104: THE HODGKIN-HUXLEY EQUATIONS

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

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells

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

Identification of Odors by the Spatiotemporal Dynamics of the Olfactory Bulb. Outline

Identification of Odors by the Spatiotemporal Dynamics of the Olfactory Bulb. Outline Identification of Odors by the Spatiotemporal Dynamics of the Olfactory Bulb Henry Greenside Department of Physics Duke University Outline Why think about olfaction? Crash course on neurobiology. Some

More information

Threshold Curve for the Excitability of Bidimensional Spiking Neurons

Threshold Curve for the Excitability of Bidimensional Spiking Neurons Threshold Curve for the Excitability of Bidimensional Spiking Neurons Arnaud Tonnelier INRIA-Grenoble 655 avenue de l Europe, Montbonnot 38334 Saint Ismier, France (Dated: July 29, 214) We shed light on

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

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

Comparing integrate-and-fire models estimated using intracellular and extracellular data 1

Comparing integrate-and-fire models estimated using intracellular and extracellular data 1 Comparing integrate-and-fire models estimated using intracellular and extracellular data 1 Liam Paninski a,b,2 Jonathan Pillow b Eero Simoncelli b a Gatsby Computational Neuroscience Unit, University College

More information

From neuronal oscillations to complexity

From 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 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

Sean Escola. Center for Theoretical Neuroscience

Sean Escola. Center for Theoretical Neuroscience Employing hidden Markov models of neural spike-trains toward the improved estimation of linear receptive fields and the decoding of multiple firing regimes Sean Escola Center for Theoretical Neuroscience

More information

Systems Biology: Theoretical Biology. Kirsten ten Tusscher, Theoretical Biology, UU

Systems Biology: Theoretical Biology. Kirsten ten Tusscher, Theoretical Biology, UU Systems Biology: Theoretical Biology Kirsten ten Tusscher, Theoretical Biology, UU 1 Chapter 9 Hodgkin-Huxley model 2 Introduction Hodgkin-Huxley model : detailed model for nerve action potential Perfect

More information

Determining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model

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

Dynamical phase transitions in periodically driven model neurons

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

This script will produce a series of pulses of amplitude 40 na, duration 1ms, recurring every 50 ms.

This script will produce a series of pulses of amplitude 40 na, duration 1ms, recurring every 50 ms. 9.16 Problem Set #4 In the final problem set you will combine the pieces of knowledge gained in the previous assignments to build a full-blown model of a plastic synapse. You will investigate the effects

More information

When is an Integrate-and-fire Neuron like a Poisson Neuron?

When is an Integrate-and-fire Neuron like a Poisson Neuron? When is an Integrate-and-fire Neuron like a Poisson Neuron? Charles F. Stevens Salk Institute MNL/S La Jolla, CA 92037 cfs@salk.edu Anthony Zador Salk Institute MNL/S La Jolla, CA 92037 zador@salk.edu

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

Hybrid Integrate-and-Fire Model of a Bursting Neuron

Hybrid Integrate-and-Fire Model of a Bursting Neuron LETTER Communicated by John Rinzel Hybrid Integrate-and-Fire Model of a Bursting Neuron Barbara J. Breen bbreen@ece.gatech.edu William C. Gerken wgerken@ece.gatech.edu Robert J. Butera, Jr. rbutera@ece.gatech.edu

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

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

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

Visual Selection and Attention Shifting Based on FitzHugh-Nagumo Equations

Visual Selection and Attention Shifting Based on FitzHugh-Nagumo Equations Visual Selection and Attention Shifting Based on FitzHugh-Nagumo Equations Haili Wang, Yuanhua Qiao, Lijuan Duan, Faming Fang, Jun Miao 3, and Bingpeng Ma 3 College of Applied Science, Beijing University

More information

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period

More information

Neuroscience applications: isochrons and isostables. Alexandre Mauroy (joint work with I. Mezic)

Neuroscience 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 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

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

arxiv: v1 [physics.bio-ph] 2 Jul 2008

arxiv: v1 [physics.bio-ph] 2 Jul 2008 Modeling Excitable Systems Jarrett L. Lancaster and Edward H. Hellen University of North Carolina Greensboro, Department of Physics and Astronomy, Greensboro, NC 27402 arxiv:0807.0451v1 [physics.bio-ph]

More information

Event-driven simulations of nonlinear integrate-and-fire neurons

Event-driven simulations of nonlinear integrate-and-fire neurons Event-driven simulations of nonlinear integrate-and-fire neurons A. Tonnelier, H. Belmabrouk, D. Martinez Cortex Project, LORIA, Campus Scientifique, B.P. 239, 54 56 Vandoeuvre-lès-Nancy, France Abstract

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

Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons

Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons PHYSICAL REVIEW E 69, 051918 (2004) Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons Magnus J. E. Richardson* Laboratory of Computational Neuroscience, Brain

More information

Fast neural network simulations with population density methods

Fast neural network simulations with population density methods Fast neural network simulations with population density methods Duane Q. Nykamp a,1 Daniel Tranchina b,a,c,2 a Courant Institute of Mathematical Science b Department of Biology c Center for Neural Science

More information

2013 NSF-CMACS Workshop on Atrial Fibrillation

2013 NSF-CMACS Workshop on Atrial Fibrillation 2013 NSF-CMACS Workshop on A Atrial Fibrillation Flavio H. Fenton School of Physics Georgia Institute of Technology, Atlanta, GA and Max Planck Institute for Dynamics and Self-organization, Goettingen,

More information

Spike-Frequency Adaptation: Phenomenological Model and Experimental Tests

Spike-Frequency Adaptation: Phenomenological Model and Experimental Tests Spike-Frequency Adaptation: Phenomenological Model and Experimental Tests J. Benda, M. Bethge, M. Hennig, K. Pawelzik & A.V.M. Herz February, 7 Abstract Spike-frequency adaptation is a common feature of

More information

Ionic basis of the resting membrane potential. Foundations in Neuroscience I, Oct

Ionic basis of the resting membrane potential. Foundations in Neuroscience I, Oct Ionic basis of the resting membrane potential Foundations in Neuroscience I, Oct 3 2017 The next 4 lectures... - The resting membrane potential (today) - The action potential - The neural mechanisms behind

More information

The Resonate-and-fire Neuron: Time Dependent and Frequency Selective Neurons in Neural Networks

The Resonate-and-fire Neuron: Time Dependent and Frequency Selective Neurons in Neural Networks Bucknell University Bucknell Digital Commons Master s Theses Student Theses 2010 The Resonate-and-fire Neuron: Time Dependent and Frequency Selective Neurons in Neural Networks Himadri Mukhopadhyay Bucknell

More information

MATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev

MATH 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 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

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

Physiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES

Physiology Unit 2. MEMBRANE POTENTIALS and SYNAPSES Physiology Unit 2 MEMBRANE POTENTIALS and SYNAPSES Neuron Communication Neurons are stimulated by receptors on dendrites and cell bodies (soma) Ligand gated ion channels GPCR s Neurons stimulate cells

More information

Activity Driven Adaptive Stochastic. Resonance. Gregor Wenning and Klaus Obermayer. Technical University of Berlin.

Activity Driven Adaptive Stochastic. Resonance. Gregor Wenning and Klaus Obermayer. Technical University of Berlin. Activity Driven Adaptive Stochastic Resonance Gregor Wenning and Klaus Obermayer Department of Electrical Engineering and Computer Science Technical University of Berlin Franklinstr. 8/9, 187 Berlin fgrewe,obyg@cs.tu-berlin.de

More information

Model Neurons I: Neuroelectronics

Model Neurons I: Neuroelectronics Chapter 5 Model Neurons I: Neuroelectronics 5.1 Introduction A great deal is known about the biophysical mechanisms responsible for generating neuronal activity, and these provide a basis for constructing

More information

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

Deconstructing Actual Neurons

Deconstructing Actual Neurons 1 Deconstructing Actual Neurons Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference: The many ionic

More information

Wave Pinning, Actin Waves, and LPA

Wave Pinning, Actin Waves, and LPA Wave Pinning, Actin Waves, and LPA MCB 2012 William R. Holmes Intercellular Waves Weiner et al., 2007, PLoS Biology Dynamic Hem1 waves in neutrophils Questions? How can such waves / pulses form? What molecular

More information

Lumping Izhikevich neurons

Lumping Izhikevich neurons Visser and Van Gils EPJ Nonlinear Biomedical Physics 2014, 2:6 RESEARCH Open Access Lumping Izhikevich neurons Sid Visser 1,2,3* and Stephan A Van Gils 2,3 *Correspondence: sid.visser@nottingham.ac.uk

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

Electronics 101 Solving a differential equation Incorporating space. Numerical Methods. Accuracy, stability, speed. Robert A.

Electronics 101 Solving a differential equation Incorporating space. Numerical Methods. Accuracy, stability, speed. Robert A. Numerical Methods Accuracy, stability, speed Robert A. McDougal Yale School of Medicine 21 June 2016 Hodgkin and Huxley: squid giant axon experiments Top: Alan Lloyd Hodgkin; Bottom: Andrew Fielding Huxley.

More information

BIFURCATION ANALYSIS OF A GENERAL CLASS OF NONLINEAR INTEGRATE-AND-FIRE NEURONS

BIFURCATION ANALYSIS OF A GENERAL CLASS OF NONLINEAR INTEGRATE-AND-FIRE NEURONS SIAM J. APPL. MATH. Vol. 68, No. 4, pp. 45 79 c 28 Society for Industrial and Applied Mathematics BIFURCATION ANALYSIS OF A GENERAL CLASS OF NONLINEAR INTEGRATE-AND-FIRE NEURONS JONATHAN TOUBOUL Abstract.

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

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

Neuronal Dynamics: Computational Neuroscience of Single Neurons

Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 5 part 3a :Three definitions of rate code Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 5 Variability and Noise: The question of the neural code Wulfram Gerstner EPFL, Lausanne,

More information

Biological Modeling of Neural Networks

Biological Modeling of Neural Networks Week 4 part 2: More Detail compartmental models Biological Modeling of Neural Networks Week 4 Reducing detail - Adding detail 4.2. Adding detail - apse -cable equat Wulfram Gerstner EPFL, Lausanne, Switzerland

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

Naseem Demeri. Mohammad Alfarra. Mohammad Khatatbeh

Naseem Demeri. Mohammad Alfarra. Mohammad Khatatbeh 7 Naseem Demeri Mohammad Alfarra Mohammad Khatatbeh In the previous lectures, we have talked about how the difference in permeability for ions across the cell membrane can generate a potential. The potential

More information

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Figure 1: Equivalent electric (RC) circuit of a neurons membrane Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of

More information

arxiv: v2 [math.ds] 26 Jun 2018

arxiv: v2 [math.ds] 26 Jun 2018 Stochastic FitzHugh-Nagumo neuron model in excitable regime embeds a leaky integrate-and-fire model MARIUS E. YAMAKOU, TAT DAT TRAN, LUU HOANG DUC, AND JÜRGEN JOST arxiv:806.0749v2 math.ds] 26 Jun 208

More information

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring. Elementary neuron models -- conductance based -- modelers alternatives

Basic elements of neuroelectronics -- membranes -- ion channels -- wiring. Elementary neuron models -- conductance based -- modelers alternatives Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wiring neurons together -- synapses

More information