Lumping Izhikevich neurons

We present the construction of a planar vector field that yields the firing rate of a bursting Izhikevich neuron can be read out, while leaving the sub-threshold behavior intact. This planar vector field is used to derive lumped formulations of two complex heterogeneous networks of bursting Izhikevich neurons. In both cases, the lumped model is compared with the spiking network. There is excellent agreement in terms of duration and number of action potentials within the bursts, but there is a slight mismatch of the burst frequency. The lumped model accurately accounts for both intrinsic bursting and post inhibitory rebound potentials in the neuron model, features which are absent in prevalent neural mass models.

[1]  R. FitzHugh.  Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[2]  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.

[3]  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.

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

[5]  C. Morris,et al.  Voltage oscillations in the barnacle giant muscle fiber. , 1981, Biophysical journal.

[6]  J. Cowan,et al.  Excitatory and inhibitory interactions in localized populations of model neurons. , 1972, Biophysical journal.

[7]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[8]  Jonathan Touboul,et al.  Spiking Dynamics of Bidimensional Integrate-and-Fire Neurons , 2009, SIAM J. Appl. Dyn. Syst..

[9]  Philip Holmes,et al.  Neural Dynamics, Bifurcations, and Firing Rates in a Quadratic Integrate-and-Fire Model with a Recovery Variable. I: Deterministic Behavior , 2012, Neural Computation.

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

[11]  P. Bressloff.  Spatiotemporal dynamics of continuum neural fields , 2012 .

[12]  Lawrence Sirovich,et al.  On the Simulation of Large Populations of Neurons , 2004, Journal of Computational Neuroscience.

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

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

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

[16]  S. Amari.  Dynamics of pattern formation in lateral-inhibition type neural fields , 1977, Biological Cybernetics.

[17]  Donald O. Walter,et al.  Mass action in the nervous system , 1975 .

[18]  R. Jindra.  Mass action in the nervous system W. J. Freeman, Academic Press, New York (1975), 489 pp., (hard covers). $34.50 , 1976, Neuroscience.

[19]  F. H. Lopes da Silva,et al.  Model of brain rhythmic activity , 1974, Kybernetik.

[20]  J. Cowan,et al.  A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue , 1973, Kybernetik.