Synchrony in Stochastic Pulse-coupled Neuronal Network Models

Size: px
Start display at page:

Download "Synchrony in Stochastic Pulse-coupled Neuronal Network Models"

Transcription

1 Synchrony in Stochastic Pulse-coupled Neuronal Network Models Katie Newhall Gregor Kovačič and Peter Kramer Aaditya Rangan and David Cai 2 Rensselaer Polytechnic Institute, Troy, New York 2 Courant Institute, New York, New York January 9, 29 Stochastic Models in Neuroscience Marseille, France

2 Motivation and Goals Motivation: Neuronal Model Behavior Study the most basic point neuron models, network models, and coarse-grained models for use in computations Develop new analytical techniques Look for simple, universal network mechanisms Goals: Understand Oscillatory Dynamics of IF Model Generalize classical problem of perfectly synchronous Integrate-and-Fire (IF) network oscillations to stochastic setting Quantitative analysis of mechanism sustaining synchrony Lessons for understanding more complicated models and types of oscillations?

3 Outline Network dynamics: model and simulations Characterization of mechanism for sustaining stochastically driven synchronized dynamics Computation of probability of repeated total firing events Computation of firing rates based on first passage time First passage time calculation Approximation of first passage time

4 Excitatory Neuronal Network All-to-all connected, current based, integrate-and-fire (IF) network dv i (t) dt = g L (v i (t) V R ) + I i (t), i =,...,N I i (t) = f k δ(t t ik ) + S δ(t t jk ) N j i k v v3 v2 v4 V R v i (t) < V T V R : reset potential (=) V T : firing threshold (=) g L : leakage rate (=) f and S: coupling strengths (> )

5 Integrate & Fire Dynamics v i (t) evolves according to the voltage ODE until t ik, when v i ( t ik ) V T At t ik, v i (t) is reset to V R : v i ( t + ik ) = V R S/N is added to the voltage, v j, for j =,...,N,j i Reset: Membrane Potential: V R V T Voltage Time

6 Classification by Mean External Driving Strength Mean voltage driven to V R + fν/g L without network coupling Superthreshold V R + fν/g L > V T Subthreshold V R + fν/g L < V T voltage time f =.,ν = 2 fν =.2 > voltage f =.,ν = 9 fν =.9 < time

7 Classification by Fluctuation Size Approximate Poisson point process with rate ν and amplitude f by drift-diffusion process with drift coef fν and diffusion coef f 2 ν/2 Diffusion Approx. Valid Diffusion Approx. Not Valid Poisson Process Driven Poisson Process Driven voltage.5 voltage Drift Diffusion Process Driven Drift Diffusion Process Driven voltage time f =.5 (V T V R )/g L ν = 24,fν =.2,N = voltage time f =. = O((V T V R )/g L ) ν = 2,fν =.2,N =

8 Types of Synchronous Firing: Partial Synchrony All-to-all connected network of N = neurons number of firing events neuron number time f =.,fν =.2,S = time

9 Types of Synchronous Firing: Imperfect Synchrony All-to-all connected network of N = neurons 2 number of firing events neuron number time f =.,fν =.2,S = time

10 Types of Synchronous Firing: Perfect Synchrony All-to-all connected network of N = neurons Consistent total firing events; synchronizable network 2 number of firing events neuron number time f =.,fν =.2,S = time

11 Model for Self-Sustaining Synchrony Fluctuations in input desynchronize the network Between total firing events N neurons behave independently Pulse coupling synchronizes the network First neuron firing causes cascading event, pulling all other neurons over threshold due to synaptic coupling (S) After each total firing event, all neurons reset, process repeats Which systems are synchronous? What is the mean firing rate of the network?

12 () What is the probability to see repeated total firing events?

13 Model for Synchronous Firing Voltages of the other N neurons when the first neuron fires: Cascade-Susceptible Not Cascade-Susceptible number of neurons number of neurons voltage f =.8,fν = voltage f =.8,fν =. S = 2,N = (Bin size = S/N)

14 For What Parameters is Network Synchronizable? number of neurons < 5 S/N voltage Cascade-susceptibility described by event: C = N k= Neuron will fire in a cascade if instantaneous coupling from other firing neurons pushes its voltage over threshold Probability all neurons included in cascading event: P(C) C k where C k : V T V (k) (N k) S N C,C k [V R,V T ] N and V () V (2) V (N )

15 Computation of Cascade-Susceptibility Probability Condition upon time of first threshold crossing: P(C) = P(C T () = t)p () T (t)dt. Approximate condition: max neuron at threshold at time t P(C T () = t) P(C V (N) (t) = V T ) Manipulate using elementary probability, P(C V (N) (t) = V T ) = P(C c V (N) (t) = V T ) where A k C c k N j=k+ (C j ) N = P(A j V (N) (t) = V T ) j= denotes event that cascade fails first at kth neuron

16 Computation of P(A j V (N) (t) = V T ) A 4 : Reduce to elementary combinatorial calculation of how neurons are distributed in bins (width S/N), each independently distributed with p v (x,t) VT for x [V R,V T ], p v V (N)(x,t) V R p v (x,t)dx otherwise, Approximate single-neuron freely evolving voltage pdf p v (x,t) as Gaussian Sum probabilities of each configuration consistent with A j V T

17 Probability that Network is Cascade-Susceptible Larger fluctuations need larger coupling to maintain synchrony.8 N= N=5 N=25 N=.8 P(C).6.4 P(C) S S/N fν =.2,f =. P(C).6.4 f=. f=.5.2 f=. f=.2 f= S/N fν =.2,N =

18 Probability that Network is Cascade-Susceptible What about for networks where connections are not all-to-all? Synaptic failure: p f Sparsity: p c.8.8 P(C) p f =. p f =.25 p f =.5 p f = S N =,fν =.2,f =. P(C) p c =. p c =.25 p c =.5 p c = S N =,fν =.2,N = Adjust P(C) by taking p v V (N)(x,t) { ( pf ) ( p c ) p v V (N)(x,t)

19 Probability that Network is Cascade-Susceptible What about scale-free networks? Effect of local topology Distribution of connections: P(K = k) k 3 Account for topology in P(A 2 ): L K A=4 A B K =3 B L L=2 red only neurons that fire blue neurons result of clustering neuron A connects to B P(C) Consider higher order terms: P(k B ) P(k B k A,A B) S Thr.: m=3 Sim.: m=3 Thr.: m=5 Sim.: m=5 Thr.: m= Sim.: m= N = 4,fν =.2,f =. (also with M. Shkarayev)

20 Probability that Network is Cascade-Susceptible What about non instantaneous synaptic coupling? P(C) = Synaptic Delay: T D Synaptic Time Course: H(t) t e t/τ τe 2 E 8 8 neuron number 6 4 neuron number time N =,S =., T D =.2 fν =.2,f = time N =,S =.,τ E =.2 fν =.2,f =.,

21 (2) What is the mean time between total firing events?

22 Synchronous First Exit Problem N neurons - don t want to solve N-dimensional mean exit time problem! Obtain PDF of one neuron s exit time, p T (t), and find PDF of minimum exit time, p () T, of N independent neurons neuron number Mean first passage time of N neurons t = tp () T (t)dt p () T (t) = Np T(t)( F T (t)) N time

23 Synchronous First Exit Problem pv(x, t) t=.6 t=3.3 t=4.9 t=6.6 FT(t) voltage time Neuron reaches threshold - removed from system (absorbed) Probability neuron not fired is probability still in V R V T P(T t) = F T (t) = VT V R p v (x,t)dx p v (x,t) solves Kolmogorov Forward Equation (KFE)

24 Single Neuron Kolmogorov Forward Equation KFE: t p v(x,t) = x [g L(x V R )p v (x,t)] + ν [p v (x f,t) p v (x,t)] Taylor Expand KFE - Diffusion Approximation: t p v(x,t) = x [(g L(x V R ) fν)p v (x,t)] + f2 ν 2 2 x 2p v(x,t) Boundary Conditions: absorbing at V T : p v (V T,t) = reflecting at V R : J[p v ](V R,t) = Flux: J[p v ](x,t) = [(g L (x V R ) fν)p v (x,t)] f2 ν 2 x p v(x,t)

25 Solution to Single Neuron KFE Eigenfunction expansion: p v (x,t) = p s (x) A n Q n (x)e λnt n= Define: z(x) = g L(x V R ) fν f νg L Eigenfunctions are Confluent Hypergeometric Functions: ( ) ( ) c M λn 2g Q n (x) = L, 2,z2 (x) + c 2 U λn 2g L, 2,z2 (x) ( ) ( ) c 3 M λn 2g L, 2,z2 (x) + c 4 U λn 2g L, 2,z2 (x) for z(x) < for z(x) Stationary Solution: p s (x) = Ne z2 (x) Initial Conditions: p v (x,) = δ(x V R ) A n = Q n(v R ) R VT V p R s(x)q 2 n (x)dx

26 Synchronous Gain Curves dashed - Gaussian approx. solid - first passage time circles - simulations black - deterministic rate Good agreement where diffusion approximation valid For fixed (superthreshold) fν, dependence f (a) firing rate 2.5 m /ˆτ.5..2 f.5 f=.5 f=. f=.2 f= fν N=, S=5.

27 Synchronous Gain Curves dashed - Gaussian approx. solid - first passage time circles - simulations For fixed (superthreshold) fν, dependence ln N (b) firing rate N= N=5 N= m /ˆτ ln(n) fν.4 f =.,S = 2.

28 Summary Synchronous firing over a wide range of parameters Balance between fluctuations and coupling strength Synchronous firing rate in terms of single neuron properties Driven by time first neuron crosses threshold Network properties are calculated accurately through exit time problems where diffusion approximation is valid This work was supported by a NSF graduate fellowship (Newhall & al. (2) Comm in Math Sci, Vol. 8, No. 2, pp. 54-6)

Cascade-Induced Synchrony in Stochastically-Driven Neuronal Networks

Cascade-Induced Synchrony in Stochastically-Driven Neuronal Networks Cascade-Induced Synchrony in Stochastically-Driven Neuronal Networks Katherine A. Newhall, Gregor Kovačič, Peter R. Kramer, and David Cai 2,3, Mathematical Sciences Department, Rensselaer Polytechnic Institute,

More information

DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS

DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS KATHERINE A. NEWHALL, GREGOR KOVAČIČ, PETER R. KRAMER, DOUGLAS ZHOU, AADITYA V. RANGAN, AND DAVID CAI Abstract. Synchronous

More information

Dedicated to the sixtieth birthday of Professor Andrew Majda

Dedicated to the sixtieth birthday of Professor Andrew Majda COMMUN. MATH. SCI. Vol. 8, No. 2, pp. 541 6 c 21 International Press DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS KATHERINE A. NEWHALL, GREGOR KOVAČIČ, PETER R. KRAMER,

More information

DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS

DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS KATHERINE A. NEWHALL, GREGOR KOVAČIČ, PETER R. KRAMER, DOUGLAS ZHOU, AADITYA V. RANGAN, AND DAVID CAI Abstract. Synchronous

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

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

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( ) Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical

More information

THE ROLE OF FLUCTUATIONS IN COARSE-GRAINED DESCRIPTIONS OF NEURONAL NETWORKS. To Dave Levermore on his 60th birthday with friendship and appreciation

THE ROLE OF FLUCTUATIONS IN COARSE-GRAINED DESCRIPTIONS OF NEURONAL NETWORKS. To Dave Levermore on his 60th birthday with friendship and appreciation COMMUN. MATH. SCI. Vol. 1, No. 1, pp. 37 354 c 212 International Press THE ROLE OF FLUCTUATIONS IN COARSE-GRAINED DESCRIPTIONS OF NEURONAL NETWORKS DAVID CAI, LOUIS TAO, MAXIM S. SHKARAYEV, AADITYA V.

More information

1. Introduction - Reproducibility of a Neuron. 3. Introduction Phase Response Curve. 2. Introduction - Stochastic synchronization. θ 1. c θ.

1. Introduction - Reproducibility of a Neuron. 3. Introduction Phase Response Curve. 2. Introduction - Stochastic synchronization. θ 1. c θ. . Introduction - Reproducibility of a euron Science (995 Constant stimuli led to imprecise spike trains, whereas stimuli with fluctuations produced spike trains with timing reproducible to less than millisecond.

More information

Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics

Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics J Comput Neurosci (1) 8:7 66 DOI 1.17/s187-9 Pseudo-s and predictability of Hodgkin-Huxley neuronal network dynamics Yi Sun Douglas Zhou Aaditya V. Rangan David Cai Received: 16 May 9 / Revised: 3 November

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

Population dynamics of interacting spiking neurons

Population dynamics of interacting spiking neurons Population dynamics of interacting spiking neurons Maurizio Mattia* Physics Laboratory, Istituto Superiore di Sanità, INFN - Gr. Coll. Roma I, V.le Regina Elena 299, 00161 Roma, Italy Paolo Del Giudice

More information

16. Working with the Langevin and Fokker-Planck equations

16. Working with the Langevin and Fokker-Planck equations 16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation

More information

Efficient Leaping Methods for Stochastic Chemical Systems

Efficient Leaping Methods for Stochastic Chemical Systems Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some

More information

Dynamics of the Exponential Integrate-and-Fire Model with Slow Currents and Adaptation

Dynamics of the Exponential Integrate-and-Fire Model with Slow Currents and Adaptation Noname manuscript No. (will be inserted by the editor) Dynamics of the Exponential Integrate-and-Fire Model with Slow Currents and Adaptation Victor J. Barranca Daniel C. Johnson Jennifer L. Moyher Joshua

More information

LINEAR RESPONSE THEORY

LINEAR RESPONSE THEORY MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior

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

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

Entrainment and Chaos in the Hodgkin-Huxley Oscillator

Entrainment and Chaos in the Hodgkin-Huxley Oscillator Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42 Overview (1) Goal: Show that the

More information

Diffusion in the cell

Diffusion in the cell Diffusion in the cell Single particle (random walk) Microscopic view Macroscopic view Measuring diffusion Diffusion occurs via Brownian motion (passive) Ex.: D = 100 μm 2 /s for typical protein in water

More information

2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES

2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES 2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES 2.0 THEOREM OF WIENER- KHINTCHINE An important technique in the study of deterministic signals consists in using harmonic functions to gain the spectral

More information

Most Probable Escape Path Method and its Application to Leaky Integrate And Fire Neurons

Most Probable Escape Path Method and its Application to Leaky Integrate And Fire Neurons Most Probable Escape Path Method and its Application to Leaky Integrate And Fire Neurons Simon Fugmann Humboldt University Berlin 13/02/06 Outline The Most Probable Escape Path (MPEP) Motivation General

More information

1. Synchronization Phenomena

1. Synchronization Phenomena 1. Synchronization Phenomena In nature synchronization happens all the time. In mechanical systems, in biological systems, in epidemiology, basically everywhere. When we talk about synchronization we usually

More information

Is the superposition of many random spike trains a Poisson process?

Is the superposition of many random spike trains a Poisson process? Is the superposition of many random spike trains a Poisson process? Benjamin Lindner Max-Planck-Institut für Physik komplexer Systeme, Dresden Reference: Phys. Rev. E 73, 2291 (26) Outline Point processes

More information

On the Dynamics of Delayed Neural Feedback Loops. Sebastian Brandt Department of Physics, Washington University in St. Louis

On the Dynamics of Delayed Neural Feedback Loops. Sebastian Brandt Department of Physics, Washington University in St. Louis On the Dynamics of Delayed Neural Feedback Loops Sebastian Brandt Department of Physics, Washington University in St. Louis Overview of Dissertation Chapter 2: S. F. Brandt, A. Pelster, and R. Wessel,

More information

Library-based numerical reduction of the Hodgkin Huxley neuron for network simulation

Library-based numerical reduction of the Hodgkin Huxley neuron for network simulation DOI 1.17/s1827-9-151-9 Library-based numerical reduction of the Hodgkin Huxley neuron for network simulation Yi un Douglas Zhou Aaditya V. Rangan David Cai Received: 8 August 28 / Revised: 13 February

More information

Handbook of Stochastic Methods

Handbook of Stochastic Methods C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical

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

Synchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014

Synchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014 Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,

More information

Stochastic modelling of epidemic spread

Stochastic modelling of epidemic spread Stochastic modelling of epidemic spread Julien Arino Department of Mathematics University of Manitoba Winnipeg Julien Arino@umanitoba.ca 19 May 2012 1 Introduction 2 Stochastic processes 3 The SIS model

More information

Phase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable

Phase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a

More information

Katherine Newhall

Katherine Newhall Katherine Newhall knewhall@unc.edu http://knewhall.web.unc.edu/ Last modified: January 20, 2016 Institute Affiliation University of North Carolina at Chapel Hill Department of Mathematics CB# 3250 Chapel

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

Structure of Correlations in Neuronal Networks

Structure of Correlations in Neuronal Networks Structure of Correlations in Neuronal Networks Krešimir Josić University of Houston James Trousdale (UH), Yu Hu (UW) Eric Shea-Brown (UW) The Connectome Van J. Wedeen, MGH/Harvard The Connectome Van J.

More information

Stochastic modelling of epidemic spread

Stochastic modelling of epidemic spread Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca

More information

Non-equilibrium phase transitions

Non-equilibrium phase transitions Non-equilibrium phase transitions An Introduction Lecture III Haye Hinrichsen University of Würzburg, Germany March 2006 Third Lecture: Outline 1 Directed Percolation Scaling Theory Langevin Equation 2

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

Synchrony in Stochastically Driven Neuronal Networks with Complex Topologies

Synchrony in Stochastically Driven Neuronal Networks with Complex Topologies Synchrony in Stochastically Driven Neuronal Networks with Coplex Topologies Katherine A. Newhall, Maxi S. Shkarayev, 2 Peter R. Kraer, 3 Gregor Kovačič, 3 and David Cai 4,5,6, Departent of Matheatics,

More information

Introduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations.

Introduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations. 19th May 2011 Chain Introduction We will Show the relation between stochastic differential equations, Gaussian processes and methods This gives us a formal way of deriving equations for the activity of

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

CISC 3250 Systems Neuroscience

CISC 3250 Systems Neuroscience CISC 3250 Systems Neuroscience Systems Neuroscience How the nervous system performs computations How groups of neurons work together to achieve intelligence Professor Daniel Leeds dleeds@fordham.edu JMH

More information

Lecture 12: Detailed balance and Eigenfunction methods

Lecture 12: Detailed balance and Eigenfunction methods Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),

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

Derivation of Itô SDE and Relationship to ODE and CTMC Models

Derivation of Itô SDE and Relationship to ODE and CTMC Models Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)

More information

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

The homogeneous Poisson process

The homogeneous Poisson process The homogeneous Poisson process during very short time interval Δt there is a fixed probability of an event (spike) occurring independent of what happened previously if r is the rate of the Poisson process,

More information

Neural Encoding: Firing Rates and Spike Statistics

Neural Encoding: Firing Rates and Spike Statistics Neural Encoding: Firing Rates and Spike Statistics Dayan and Abbott (21) Chapter 1 Instructor: Yoonsuck Choe; CPSC 644 Cortical Networks Background: Dirac δ Function Dirac δ function has the following

More information

Path integrals for classical Markov processes

Path integrals for classical Markov processes Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field

More information

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,

More information

Dendritic cable with active spines: a modelling study in the spike-diffuse-spike framework

Dendritic cable with active spines: a modelling study in the spike-diffuse-spike framework Dendritic cable with active spines: a modelling study in the spike-diffuse-spike framework Yulia Timofeeva a, Gabriel Lord a and Stephen Coombes b a Department of Mathematics, Heriot-Watt University, Edinburgh,

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

Chapter 6. Random Processes

Chapter 6. Random Processes Chapter 6 Random Processes Random Process A random process is a time-varying function that assigns the outcome of a random experiment to each time instant: X(t). For a fixed (sample path): a random process

More information

The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is

The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is 1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles

More information

Coarse-grained event tree analysis for quantifying Hodgkin-Huxley neuronal network dynamics

Coarse-grained event tree analysis for quantifying Hodgkin-Huxley neuronal network dynamics J Comput Neurosci (212) 32:55 72 DOI 1.17/s1827-11-339-7 Coarse-grained event tree analysis for quantifying Hodgkin-Huxley neuronal network dynamics Yi Sun Aaditya V. Rangan Douglas Zhou David Cai Received:

More information

6 Continuous-Time Birth and Death Chains

6 Continuous-Time Birth and Death Chains 6 Continuous-Time Birth and Death Chains Angela Peace Biomathematics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology.

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

Fluctuation theorem in systems in contact with different heath baths: theory and experiments.

Fluctuation theorem in systems in contact with different heath baths: theory and experiments. Fluctuation theorem in systems in contact with different heath baths: theory and experiments. Alberto Imparato Institut for Fysik og Astronomi Aarhus Universitet Denmark Workshop Advances in Nonequilibrium

More information

Handbook of Stochastic Methods

Handbook of Stochastic Methods Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei

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

Theoretical connections between mathematical neuronal models corresponding to different expressions of noise

Theoretical connections between mathematical neuronal models corresponding to different expressions of noise Theoretical connections between mathematical neuronal models corresponding to different expressions of noise Grégory Dumont, Jacques Henry, Carmen Oana Tarniceriu To cite this version: Grégory Dumont,

More information

Stability of the splay state in pulse-coupled networks

Stability of the splay state in pulse-coupled networks Krakow August 2013 p. Stability of the splay state in pulse-coupled networks S. Olmi, A. Politi, and A. Torcini http://neuro.fi.isc.cnr.it/ Istituto dei Sistemi Complessi - CNR - Firenze Istituto Nazionale

More information

Ordinary differential equations. Phys 420/580 Lecture 8

Ordinary differential equations. Phys 420/580 Lecture 8 Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue

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

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

Markov Reliability and Availability Analysis. Markov Processes

Markov Reliability and Availability Analysis. Markov Processes Markov Reliability and Availability Analysis Firma convenzione Politecnico Part II: Continuous di Milano e Time Veneranda Discrete Fabbrica State del Duomo di Milano Markov Processes Aula Magna Rettorato

More information

Stochastic process. X, a series of random variables indexed by t

Stochastic process. X, a series of random variables indexed by t Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,

More information

Simulation methods for stochastic models in chemistry

Simulation methods for stochastic models in chemistry Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation

More information

Beyond blow-up in excitatory integrate and fire neuronal networks: refractory period and spontaneous activity

Beyond blow-up in excitatory integrate and fire neuronal networks: refractory period and spontaneous activity Beyond blow-up in excitatory integrate and fire neuronal networks: refractory period and spontaneous activity Maria J. Caceres, Benoît Perthame To cite this version: Maria J. Caceres, Benoît Perthame.

More information

Data assimilation with and without a model

Data assimilation with and without a model Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,

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

1 Balanced networks: Trading speed for noise

1 Balanced networks: Trading speed for noise Physics 178/278 - David leinfeld - Winter 2017 (Corrected yet incomplete notes) 1 Balanced networks: Trading speed for noise 1.1 Scaling of neuronal inputs An interesting observation is that the subthresold

More information

Localized Excitations in Networks of Spiking Neurons

Localized Excitations in Networks of Spiking Neurons Localized Excitations in Networks of Spiking Neurons Hecke Schrobsdorff Bernstein Center for Computational Neuroscience Göttingen Max Planck Institute for Dynamics and Self-Organization Seminar: Irreversible

More information

Memory and hypoellipticity in neuronal models

Memory and hypoellipticity in neuronal models Memory and hypoellipticity in neuronal models S. Ditlevsen R. Höpfner E. Löcherbach M. Thieullen Banff, 2017 What this talk is about : What is the effect of memory in probabilistic models for neurons?

More information

Effects of Noisy Drive on Rhythms in Networks of Excitatory and Inhibitory Neurons

Effects of Noisy Drive on Rhythms in Networks of Excitatory and Inhibitory Neurons Effects of Noisy Drive on Rhythms in Networks of Excitatory and Inhibitory Neurons Christoph Börgers 1 and Nancy Kopell 2 1 Department of Mathematics, Tufts University, Medford, MA 2155 2 Department of

More information

1.1. BASIC ANTI-DIFFERENTIATION 21 + C.

1.1. BASIC ANTI-DIFFERENTIATION 21 + C. .. BASIC ANTI-DIFFERENTIATION and so e x cos xdx = ex sin x + e x cos x + C. We end this section with a possibly surprising complication that exists for anti-di erentiation; a type of complication which

More information

Langevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany

Langevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace

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

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

First Passage Time Calculations

First Passage Time Calculations First Passage Time Calculations Friday, April 24, 2015 2:01 PM Homework 4 will be posted over the weekend; due Wednesday, May 13 at 5 PM. We'll now develop some framework for calculating properties of

More information

Controlled Diffusions and Hamilton-Jacobi Bellman Equations

Controlled Diffusions and Hamilton-Jacobi Bellman Equations Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter

More information

A realistic neocortical axonal plexus model has implications for neocortical processing and temporal lobe epilepsy

A realistic neocortical axonal plexus model has implications for neocortical processing and temporal lobe epilepsy A realistic neocortical axonal plexus model has implications for neocortical processing and temporal lobe epilepsy Neocortical Pyramidal Cells Can Send Signals to Post-Synaptic Cells Without Firing Erin

More information

DATA-DRIVEN TECHNIQUES FOR ESTIMATION AND STOCHASTIC REDUCTION OF MULTI-SCALE SYSTEMS

DATA-DRIVEN TECHNIQUES FOR ESTIMATION AND STOCHASTIC REDUCTION OF MULTI-SCALE SYSTEMS DATA-DRIVEN TECHNIQUES FOR ESTIMATION AND STOCHASTIC REDUCTION OF MULTI-SCALE SYSTEMS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial Fulfillment

More information

Synchrony in Neural Systems: a very brief, biased, basic view

Synchrony in Neural Systems: a very brief, biased, basic view Synchrony in Neural Systems: a very brief, biased, basic view Tim Lewis UC Davis NIMBIOS Workshop on Synchrony April 11, 2011 components of neuronal networks neurons synapses connectivity cell type - intrinsic

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

Quantifying Intermittent Transport in Cell Cytoplasm

Quantifying Intermittent Transport in Cell Cytoplasm Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent

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

A Population Density Approach that Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning

A Population Density Approach that Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning Journal of Computational Neuroscience, 8, 19 5 (2) c 2 Kluwer Academic Publishers, Boston. Manufactured in The Netherlands. A Population Density Approach that Facilitates Large-Scale Modeling of Neural

More information

Songting Li. Applied Mathematics, Mathematical and Computational Neuroscience, Biophysics

Songting Li. Applied Mathematics, Mathematical and Computational Neuroscience, Biophysics Songting Li Contact Information Phone: +1-917-930-3505 email: songting@cims.nyu.edu homepage: http://www.cims.nyu.edu/ songting/ Address: Courant Institute, 251 Mercer Street, New York, NY, United States,

More information

A Mathematical Study of Electrical Bursting in Pituitary Cells

A Mathematical Study of Electrical Bursting in Pituitary Cells A Mathematical Study of Electrical Bursting in Pituitary Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Collaborators on

More 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

Patterns of Synchrony in Neural Networks with Spike Adaptation

Patterns of Synchrony in Neural Networks with Spike Adaptation Patterns of Synchrony in Neural Networks with Spike Adaptation C. van Vreeswijky and D. Hanselz y yracah Institute of Physics and Center for Neural Computation, Hebrew University, Jerusalem, 9194 Israel

More information

Dynamical modelling of systems of coupled oscillators

Dynamical modelling of systems of coupled oscillators Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)

More information

Towards a theory of cortical columns: From spiking neurons to interacting neural populations of finite size. Abstract.

Towards a theory of cortical columns: From spiking neurons to interacting neural populations of finite size. Abstract. Schwalger T, Deger M, Gerstner W. PLoS Comput Biol. 2017;13(4):e1005507 Towards a theory of cortical columns: From spiking neurons to interacting neural populations of finite size Tilo Schwalger 1,*, Moritz

More information

Lecture Notes 7 Random Processes. Markov Processes Markov Chains. Random Processes

Lecture Notes 7 Random Processes. Markov Processes Markov Chains. Random Processes Lecture Notes 7 Random Processes Definition IID Processes Bernoulli Process Binomial Counting Process Interarrival Time Process Markov Processes Markov Chains Classification of States Steady State Probabilities

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechanics Rajdeep Sensarma! sensarma@theory.tifr.res.in Lecture #22 Path Integrals and QM Recap of Last Class Statistical Mechanics and path integrals in imaginary time Imaginary time

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

EEE105 Teori Litar I Chapter 7 Lecture #3. Dr. Shahrel Azmin Suandi Emel:

EEE105 Teori Litar I Chapter 7 Lecture #3. Dr. Shahrel Azmin Suandi Emel: EEE105 Teori Litar I Chapter 7 Lecture #3 Dr. Shahrel Azmin Suandi Emel: shahrel@eng.usm.my What we have learnt so far? Chapter 7 introduced us to first-order circuit From the last lecture, we have learnt

More information

Lecture 12: Detailed balance and Eigenfunction methods

Lecture 12: Detailed balance and Eigenfunction methods Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner

More information

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose

More information