Non-smooth Bifurcations of Mean Field Systems of Two-Dimensional Integrate and Fire Neurons

Mean-field systems have been recently derived that adequately predict the behaviors of large networks of coupled integrate-and-fire neurons [14]. The mean-field system for a network of neurons with spike frequency adaptation is typically a pair of differential equations for the mean adaptation and mean synaptic gating variable of the network. These differential equations are non-smooth, and in particular are piecewise smooth continuous (PWSC). Here, we analyze the smooth and non-smooth bifurcation structure of these equations and show that the system is organized around a pair of co-dimension two bifurcations that involve, respectively, the collision between a Hopf equilibrium point and a switching manifold, and a saddle-node equilibrium point and a switching manifold. These two co-dimension 2 bifurcations can coalesce into a co-dimension 3 non-smooth bifurcation. As the mean-field system we study is a non-generic piecewise smooth continuous system, we discuss possible regularizations of this system and how the bifurcations which occur are related to non-smooth bifurcations displayed by generic PWSC systems.

[1]  Alessandro Colombo,et al.  Bifurcations of piecewise smooth flows: Perspectives, methodologies and open problems , 2012 .

[2]  Mario di Bernardo,et al.  Bifurcations in Nonsmooth Dynamical Systems , 2008, SIAM Rev..

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

[4]  Carl van Vreeswijk,et al.  Patterns of Synchrony in Neural Networks with Spike Adaptation , 2001, Neural Computation.

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

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

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

[8]  Marco Antonio Teixeira,et al.  A singular approach to discontinuous vector fields on the plane , 2006 .

[9]  Wulfram Gerstner,et al.  Firing patterns in the adaptive exponential integrate-and-fire model , 2008, Biological Cybernetics.

[10]  Cheng Ly,et al.  One-Dimensional Population Density Approaches to Recurrently Coupled Networks of Neurons with Noise , 2015, SIAM J. Appl. Math..

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

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

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

[14]  Bard Ermentrout,et al.  Reduction of Conductance-Based Models with Slow Synapses to Neural Nets , 1994, Neural Computation.

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

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