One-Dimensional Population Density Approaches to Recurrently Coupled Networks of Neurons with Noise

Mean-field systems have been previously derived for networks of coupled, two-dimensional, integrate-and-fire neurons such as the Izhikevich, adapting exponential, and quartic integrate-and-fire, among others. Unfortunately, the mean-field systems have a degree of frequency error, and the networks analyzed often do not include noise when there is adaptation. Here, we derive a one-dimensional partial differential equation (PDE) approximation for the marginal voltage density under a first order moment closure for coupled networks of integrate-and-fire neurons with white noise inputs. The PDE has substantially less frequency error than the mean-field system and provides a great deal more information, at the cost of analytical tractability. The convergence properties of the mean-field system in the low noise limit are elucidated. A novel method for the analysis of the stability of the asynchronous tonic firing solution is also presented and implemented. Unlike in previous attempts at stability analysis with these network types, information about the marginal densities of the adaptation variables is used. This method can in principle be applied to other systems with nonlinear PDEs.

[1]  M. J. Richardson,et al.  Firing-rate response of linear and nonlinear integrate-and-fire neurons to modulated current-based and conductance-based synaptic drive. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  D. Hansel,et al.  Existence and stability of persistent states in large neuronal networks. , 2001, Physical review letters.

[3]  Germán Mato,et al.  Asynchronous States and the Emergence of Synchrony in Large Networks of Interacting Excitatory and Inhibitory Neurons , 2003, Neural Computation.

[4]  Bruce W. Knight,et al.  Dynamics of Encoding in Neuron Populations: Some General Mathematical Features , 2000, Neural Computation.

[5]  Klaus Obermayer,et al.  How adaptation shapes spike rate oscillations in recurrent neuronal networks , 2012, Front. Comput. Neurosci..

[6]  Nicolas Brunel,et al.  Dynamics of the Firing Probability of Noisy Integrate-and-Fire Neurons , 2002, Neural Computation.

[7]  Nicolas Brunel,et al.  Analytical approximations of the firing rate of an adaptive exponential integrate-and-fire neuron in the presence of synaptic noise , 2014, Front. Comput. Neurosci..

[8]  Eugene M. Izhikevich,et al.  Simple model of spiking neurons , 2003, IEEE Trans. Neural Networks.

[9]  Eugene M. Izhikevich,et al.  Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting , 2006 .

[10]  Paul C. Bressloff,et al.  Stochastic Neural Field Theory and the System-Size Expansion , 2009, SIAM J. Appl. Math..

[11]  Frances K. Skinner,et al.  Simple, biologically-constrained CA1 pyramidal cell models using an intact, whole hippocampus context , 2014, F1000Research.

[12]  A. Föhrenbach,et al.  SIMPLE++ , 2000, OR Spectr..

[13]  Abbott,et al.  Asynchronous states in networks of pulse-coupled oscillators. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  Duane Q. Nykamp,et al.  A kinetic theory approach to capturing interneuronal correlation: the feed-forward case , 2008, Journal of Computational Neuroscience.

[15]  Alla Borisyuk,et al.  Fluctuation-driven rhythmogenesis in an excitatory neuronal network with slow adaptation , 2008, Journal of Computational Neuroscience.

[16]  Nicolas Brunel,et al.  Fast Global Oscillations in Networks of Integrate-and-Fire Neurons with Low Firing Rates , 1999, Neural Computation.

[17]  Magnus J. E. Richardson,et al.  Spike-train spectra and network response functions for non-linear integrate-and-fire neurons , 2008, Biological Cybernetics.

[18]  Willy Govaerts,et al.  MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs , 2003, TOMS.

[19]  Wilten Nicola,et al.  Mean-field models for heterogeneous networks of two-dimensional integrate and fire neurons , 2013, Front. Comput. Neurosci..

[20]  Wilten Nicola,et al.  Network bursting using experimentally constrained single compartment CA3 hippocampal neuron models with adaptation , 2011, Journal of Computational Neuroscience.

[21]  Nicolas Brunel,et al.  Firing Rate of the Noisy Quadratic Integrate-and-Fire Neuron , 2003, Neural Computation.

[22]  Wulfram Gerstner,et al.  Adaptive exponential integrate-and-fire model as an effective description of neuronal activity. , 2005, Journal of neurophysiology.

[23]  Duane Q. Nykamp,et al.  A Population Density Approach That Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning , 2004, Journal of Computational Neuroscience.

[24]  Wilten Nicola,et al.  Bifurcations of large networks of two-dimensional integrate and fire neurons , 2013, Journal of Computational Neuroscience.

[25]  Wilten Nicola,et al.  Non-smooth Bifurcations of Mean Field Systems of Two-Dimensional Integrate and Fire Neurons , 2014, 1408.4767.

[26]  Michael A. Buice,et al.  Systematic Fluctuation Expansion for Neural Network Activity Equations , 2009, Neural Computation.

[27]  N. Parga,et al.  Role of synaptic filtering on the firing response of simple model neurons. , 2004, Physical review letters.

[28]  Néstor Parga,et al.  Theory of Input Spike Auto- and Cross-Correlations and Their Effect on the Response of Spiking Neurons , 2007, Neural Computation.

[29]  Aaditya V Rangan.  Diagrammatic expansion of pulse-coupled network dynamics. , 2009, Physical review letters.

[30]  Cheng Ly,et al.  A Principled Dimension-Reduction Method for the Population Density Approach to Modeling Networks of Neurons with Synaptic Dynamics , 2013, Neural Computation.

[31]  Julien Modolo,et al.  Development and validation of a neural population model based on the dynamics of a discontinuous membrane potential neuron model. , 2007, Journal of integrative neuroscience.

[32]  S. Orszag,et al.  Advanced Mathematical Methods For Scientists And Engineers , 1979 .

[33]  Cheng Ly,et al.  Population density methods for stochastic neurons with realistic synaptic kinetics: Firing rate dynamics and fast computational methods , 2006, Network.

[34]  M. J. Richardson,et al.  Dynamics of populations and networks of neurons with voltage-activated and calcium-activated currents. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Jonathan Touboul,et al.  Bifurcation Analysis of a General Class of Nonlinear Integrate-and-Fire Neurons , 2008, SIAM J. Appl. Math..

[36]  S. Orszag,et al.  Advanced mathematical methods for scientists and engineers I: asymptotic methods and perturbation theory. , 1999 .

[37]  S. Strogatz,et al.  Stability of incoherence in a population of coupled oscillators , 1991 .

[38]  Cheng Ly,et al.  Critical Analysis of Dimension Reduction by a Moment Closure Method in a Population Density Approach to Neural Network Modeling , 2007, Neural Computation.