Low-dimensional firing-rate dynamics for populations of renewal-type spiking neurons

Abstract

The macroscopic dynamics of large populations of neurons can be mathematically analyzed using low-dimensional firing-rate or neural-mass models. However, these models fail to capture spike synchronization effects of stochastic spiking neurons such as the non-stationary population response to rapidly changing stimuli. Here, we derive low-dimensional firing-rate models for homogeneous populations of general renewal-type neurons, including integrate-and-fire models driven by white noise. Renewal models account for neuronal refractoriness and spike synchronization dynamics. The derivation is based on an eigenmode expansion of the associated refractory density equation, which generalizes previous spectral methods for Fokker-Planck equations to arbitrary renewal models. We find a simple relation between the eigenvalues, which determine the characteristic time scales of the firing rate dynamics, and the Laplace transform of the interspike interval density or the survival function of the renewal process. Analytical expressions for the Laplace transforms are readily available for many renewal models including the leaky integrate-and-fire model. Retaining only the first eigenmode yields already an adequate low-dimensional approximation of the firing-rate dynamics that captures spike synchronization effects and fast transient dynamics at stimulus onset. We explicitly demonstrate the validity of our model for a large homogeneous population of Poisson neurons with absolute refractoriness, and other renewal models that admit an explicit analytical calculation of the eigenvalues. The here presented eigenmode expansion provides a systematic framework for novel firing-rate models in computational neuroscience based on spiking neuron dynamics with refractoriness.

arXiv comment24 pages, 7 figures

References in corpus (13)

• Towards a theory of cortical columns: From spiking neurons to interacting neural populations of finite size
• Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks
• Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: Comparison and implementation
• Network events on multiple space and time scales in cultured neural networks and in a stochastic rate model
• An exact firing rate model reveals the differential effects of chemical versus electrical synapses in spiking networks
• Patterns of interval correlations in neural oscillators with adaptation
• How single neuron properties shape chaotic dynamics and signal transmission in random neural networks
• Nonlinear excitations in DNA: Aperiodic models vs actual genome sequences
• Automated point-neuron simplification of data-driven microcircuit models
• Equivalence between synaptic current dynamics and heterogeneous propagation delays in spiking neuron networks
• Low-dimensional firing rate dynamics of spiking neuron networks
• Renewal Approach to the Analysis of the Asynchronous State for Coupled Noisy Oscillators
• Mesoscopic population equations for spiking neural networks with synaptic short-term plasticity

Cited by in corpus (5)

• A reduction methodology for fluctuation driven population dynamics
• A multiple time renewal equation for neural assemblies with elapsed time model
• Analytic Investigation for Spatio-temporal Patterns Propagation in Spiking Neural Networks
• Adjoint Method for Macroscopic Phase-Resetting Curves of Generic Spiking Neural Networks
• Damped oscillations of the probability of random events followed by absolute refractory period: exact analytical results