Population dynamics of interacting spiking neurons

A dynamical equation is derived for the spike emission rate nu(t) of a homogeneous network of integrate-and-fire (IF) neurons in a mean-field theoretical framework, where the activity of the single cell depends both on the mean afferent current (the field) and on its fluctuations. Finite-size effects are taken into account, by a stochastic extension of the dynamical equation for the nu; their effect on the collective activity is studied in detail. Conditions for the local stability of the collective activity are shown to be naturally and simply expressed in terms of (the slope of) the single neuron, static, current-to-rate transfer function. In the framework of the local analysis, we studied the spectral properties of the time-dependent collective activity of the finite network in an asynchronous state; finite-size fluctuations act as an ongoing self-stimulation, which probes the spectral structure of the system on a wide frequency range. The power spectrum of nu exhibits modes ranging from very high frequency (depending on spike transmission delays), which are responsible for instability, to oscillations at a few Hz, direct expression of the diffusion process describing the population dynamics. The latter diffusion slow modes do not contribute to the stability conditions. Their characteristic times govern the transient response of the network; these reaction times also exhibit a simple dependence on the slope of the neuron transfer function. We speculate on the possible relevance of our results for the change in the characteristic response time of a neural population during the learning process which shapes the synaptic couplings, thereby affecting the slope of the transfer function. There is remarkable agreement of the theoretical predictions with simulations of a network of IF neurons with a constant leakage term for the membrane potential.