Skip to main content

Advertisement

Low-dimensional spike rate dynamics of coupled adaptive model neurons

The spiking activity of single neurons can be well described by a two-dimensional integrate-and-fire model that includes neuronal adaptation [1] caused by slowly decaying potassium currents [2]. For fluctuating inputs sparsely coupled spiking model neurons exhibit stochastic population dynamics which can be effectively characterized using the Fokker-Planck equation (see, e.g., [35]). This approach leads to a model with an infinite-dimensional state space and non-standard boundary conditions. However, the spike rate dynamics can be approximated by a low-dimensional ordinary differential equation in different ways [4, 6, 7]. Although these approximation techniques are interrelated it is not clear which reduced model best reproduces the spike rate of the original spiking network, depending on the statistics of the input. Here we first extend each of these reduction methods to account for neuronal adaptation and then evaluate the reduced models in terms of spike rate reproduction accuracy for a range of biologically plausible input statistics, computational demand and implementation complexity (see, e.g., Figure 1). These reduced descriptions are well suited for (i) application in neural mass/mean-field based brain network models, having a link to single neuron properties retained and being computationally efficient, and (ii) mathematical analyses of, e.g., asynchronous and rhythmic network states.

Figure 1
figure1

Simulation of a large population of adaptive exponential integrate-and-fire (aEIF) neurons driven by a stochastic current with time-varying moments. Instantaneous spike rate and adaptation current averaged over 200,000 neurons are shown in gray. Overlaid are mean spike rate and adaptation current of two derived low-dimensional models receiving input with the same time-dependent moments as the population of aEIF neurons.

References

  1. 1.

    Brette R, Gerstner W: Adaptive exponential integrate-and-fire model as an effective description of neuronal activity. J Neurophysiol. 2005, 94: 3637-3642.

  2. 2.

    Ladenbauer J, Augustin M, Obermayer K: How adaptation currents change threshold, gain, and variability of neuronal spiking. J Neurophysiol. 2014, 111: 939-953.

  3. 3.

    Brunel N: Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. J Comput Neurosci. 2000, 8: 183-208.

  4. 4.

    Mattia M, Del Giudice P: Population dynamics of interacting spiking neurons. Phys Rev E. 2002, 66: 051917-

  5. 5.

    Augustin M, Ladenbauer J, Obermayer K: How adaptation shapes spike rate oscillations in recurrent neuronal networks. Front Comput Neurosci. 2013, 7: 9-

  6. 6.

    Schaffer E, Ostojic S, Abbott L: A complex-valued firing-rate model that approximates the dynamics of spiking networks. PLOS Comput Biol. 2013, 9: e1003301-

  7. 7.

    Ostojic S, Brunel N: From spiking neuron models to linear-nonlinear models. PLOS Comput Biol. 2011, 7: e1001056-

Download references

Acknowledgements

This work was supported by the DFG Collaborative Research Center SFB910.

Author information

Correspondence to Moritz Augustin.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Model Neuron
  • Single Neuron
  • Spike Activity
  • Spike Rate
  • Implementation Complexity