Synchrony in Stochastic Pulse-coupled Neuronal Network Models
|
|
- Ruth Henderson
- 5 years ago
- Views:
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 Katherine A. Newhall, Gregor Kovačič, Peter R. Kramer, and David Cai 2,3, Mathematical Sciences Department, Rensselaer Polytechnic Institute,
More informationDYNAMICS 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 informationDedicated 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 informationDYNAMICS 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 informationFast 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 informationConsider 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 informationMathematical 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 informationTHE 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 information1. 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 informationPseudo-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 informationNeuronal 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 informationPopulation 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 information16. 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 informationEfficient 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 informationDynamics 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 informationLINEAR 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 informationPhase 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 informationActivity 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 informationEntrainment and Chaos in the Hodgkin-Huxley Oscillator
Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42 Overview (1) Goal: Show that the
More informationDiffusion 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 information2. 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 informationMost 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 information1. 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 informationIs 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 informationOn 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 informationLibrary-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 informationHandbook 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 informationSynaptic 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 informationSynchronization 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 informationStochastic 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 informationPhase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable
Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a
More informationKatherine 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 informationCSE/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 informationStructure 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 informationStochastic 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 informationNon-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 informationSupporting 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 informationSynchrony 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 informationIntroduction. 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 informationNE 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 informationCISC 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 informationLecture 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 informationTopics 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 informationDerivation 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 informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationThe 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 informationNeural 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 informationPath 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 informationShort-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 informationDendritic 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 informationApproximate, 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 informationChapter 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 informationThe 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 informationCoarse-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 information6 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 informationPhase 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 informationFluctuation 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 informationHandbook 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 informationWe 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 informationTheoretical 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 informationStability 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 informationOrdinary 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 informationNeuronal 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 informationLecture 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 informationMarkov 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 informationStochastic 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 informationSimulation 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 informationBeyond 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 informationData 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 informationWhen 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 information1 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 informationLocalized 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 informationMemory 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 informationEffects 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 information1.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 informationLangevin 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 informationComparing 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 informationComputing 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 informationFirst 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 informationControlled 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 informationA 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 informationDATA-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 informationSynchrony 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 informationProbabilistic 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 informationQuantifying 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 informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More informationA 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 informationSongting 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 informationA Mathematical Study of Electrical Bursting in Pituitary Cells
A Mathematical Study of Electrical Bursting in Pituitary Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Collaborators on
More informationInformation 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 informationPatterns 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 informationDynamical 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 informationTowards 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 informationLecture 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 informationAdvanced 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 informationCoupling 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 informationEEE105 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 informationLecture 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 informationON 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