# Spike-response model

Wulfram Gerstner (2008), Scholarpedia, 3(12):1343. | doi:10.4249/scholarpedia.1343 | revision #91800 [link to/cite this article] |

The ** Spike Response Model **
is a generalization of the leaky integrate-and-fire model
and gives a simple description of action potential generation in neurons.
Just as in the integrate-and-fire model, action potentials are generated
when the voltage passes a threshold from below. In contrast to the leaky integrate-and-fire model, the spike response model includes *refractoriness*.
A notational difference is that integrate-and-fire models
are formulated using differential equations for the
voltage, whereas the Spike Response Model is formulated using filters.

## Mathematical formulation

The membrane potential in the spike response model is given by \[\tag{1} u(t) = \eta(t-\hat{t}) + \int_0^\infty \kappa(t-\hat{t},s) I(t-s) ds \]

where \(\hat{t}\) is the firing time of the last spike of the neuron, \(\eta\) describes the form of the action potential and its spike after-potential, \(\kappa\) the linear response to an input pulse and \(I(t)\) a stimulating current. The next spike occurs if the membrane potential \(u\) hits a threshold \(\theta(t-\hat{t})\) from below in which case \(\hat{t}\) is updated.

- The name spike response model stems from the fact that \(\kappa\) describes the
*response*of the neuron to an incoming short pulse (such as a spike arriving from another neuron) and \(\eta\) describes the response of the membrane to its own spike. The functions \(\eta\) and \(\kappa\) are also called kernels, \(\kappa\) is also called the linear filter of the membrane.

- From the mathematical point of view it is not necessary to keep both a dynamic threshold and the spike shape \(\eta\ ,\) since only the difference \(x=u - \theta\) between the membrane potential and the threshold matters for the spike dynamics.

## Main features and examples

- The threshold \(\theta\) is not fixed but depends on the time since the last spike. Typically the threshold is higher immediately after a spike and decays then back to its resting value.

- The spike shape \(\eta\) is a function of the time since the last spike. It can describe a depolarizing, hyperpolarizing, or resonating spike-after potential (see Figure 3 and Figure 2).

- The responsiveness \(\kappa\) to an input pulse depends on the time since the last spike, because typically the effective membrane time constant after a spike is shorter, since many ion channels are open.

- The time course of the response \(\kappa\) can include a single exponential, combinations of exponentials with different time constants, or resonating behavior in form of a delayed oscillation. This is the case if the standard Hodgkin-Huxley model is approximated by the Spike Response Model.

- Refractoriness can be modelled as a combination of increased threshold, hyperpolarizing afterpotential, and reduced responsiveness after a spike, as observed in real neurons (Badel et al., 2008).

## Fits to experimental data

The Spike Response Model can be fitted to experimental data where a neuron is stimulated by a rapidly varying time dependent current or conductance, see Figure 6 To do so, use the following steps.

- subtract from the experimental data the resting potential. This gives a normalized voltage trace \(u(t)\ .\)

- align spikes and determine the mean shape of the spike and spike-afterpotential. This gives \(\eta(t-\hat{t})\ .\)

- remove spikes and calculate the subthreshold membrane potential \( y(t) = u(t) - \eta(t-\hat{t})\)

- determine the best linear filter to approximate the subthreshold potential \(y(t)\) by \(\int_0^\infty \kappa(t-\hat{t},s) I^{\rm ext}(t-s) ds\ .\)

- optimize the threshold \(\theta\) so as to get the correct mean firing rate of the neuron.

The above parameter extraction procedure yields a Spike Response Model that fits experimental data to a high degree of accuracy and predicts a large fraction of spikes with a precision of +/-2ms (Jolivet et al., 2006). The same procedure has also been used to approximate detailed neuron models of the Hodgkin-Huxley type by the Spike Response Model (Kistler et al., 1997, Jolivet et al., 2004).

In a public competition of spike-time prediction
under random conductance injection, the Spike Response
Model was in the group of winning models
(Jolivet et al., 2008) whereas a standard leaky
integrate-and-fire model performed significantly worse.
These and other results show that inclusion of refractory properties
is important. Refractory properties can be measured
directly and manifest themselves in the time course of the spike afterpotenial
\(\eta(t-\hat{t})\ ,\) by a reduction in the effective membrane time constant
influencing \(\kappa(t-\hat{t})\ ,\) and by an increase of the threshold \( \theta \)
after a spike (Badel et al., 2008).

## Special cases and variants of the Spike Response Model

### Leaky Integrate-and-fire model

The leaky integrate and fire model is a special case of the Spike Response Model. To see this, take the differential equation of the leaky integrate-and-fire model \[ {dv \over dt} = -{{v - v_{eq}} \over {\tau}} + I(t) \] where \(v\) is the membrane potential, \(v_{eq}\) is the equilibrium potential and \(\tau = RC\) is the membrane time constant. Integration of the differential equation for arbitrary input \(I(t)\ ,\) yields \[\tag{2} u(t) = \eta(t-\hat{t}) + \int_0^\infty \Theta(t-\hat{t}-s) \exp(-s/\tau)\, I(t-s) ds \]

where \(u = v-v_{eq}\) and \(\eta(t-\hat{t})= (v_{reset} - v_{eq}) \, \exp[-(t-\hat{t})/\tau]\) comes from the reset after each spike. This is a special case of the Spike Response Model defined in Eq. (1) with \(\kappa(t-\hat{t},s) = \Theta(t-\hat{t}-s) \exp(-s/\tau)\) where \(\Theta\) denotes the Heaviside step function.

### Leaky Integrate-and-fire model with time-dependent time constant

The effective time constant of a neuron is different immediately after an action potential. Similarly, the momentary equilibrium potential is different and varies as a function of the time since the last spike. Hence we can write an integrate-and-fire model with time-dependent parameters \[ {dv\over dt} = -{{v-v_{eq}(t-\hat{t})}\over{\tau(t-\hat{t})}} + I(t)\] Integration of this equation gives a special case of the Spike Response Model (1).

### Spike input and synaptic response kernel \(\epsilon\)

In the case of input generated by synaptic current pulses caused by presynaptic spike arrival, the Spike Response Model can be written as \[\tag{3} u(t) = \eta(t-\hat{t}) + \sum_j \sum_f w_j \epsilon(t-\hat{t},t-t_j^f) \]

where \(w_j\) is the weight of synapse \(j\ ,\) \(t_j^f\) the arrival time of the \(f\) th spike at synapse \(j\ ,\) and \(\epsilon\) the time course of a postsynaptic potential caused by spike arrival.

To see the connection with Eq. (1), suppose that the input to the Spike Response Model consists not of an imposed current, but of synaptic input currents of amplitude \(w_j\) and time course \(\alpha(t-t_j^f)\) where \(t_j^f\) is the spike arrival time at synapse \(j\ ,\) that is, \(I(t) = \sum_j \sum_f w_j \, \alpha(t-t_j^f)\ .\) Convolution of the kernel \(\kappa\) with the current \(\alpha\) yields the postsynaptic potential \(\epsilon\ .\)

### Spike Response Model SRM\(_0\)

SRM\(_0\) is a simplified version of the Spike Response Model. It does not include a dependence of the response kernel \(\kappa\) upon the time since the last spike. Hence we have \[\tag{4} u(t) = \eta(t-\hat{t}) + \int \kappa(s) I^{\rm ext}(t-s) ds \]

instead of (1). The threshold can be dynamic as before. The variant SRM\(_0\) is easier to fit to experimental data than the full Spike Response Model, since it needs less data (Jolivet et al. 2006).

### Cumulative Spike Response Model: bursting and adaptation

In some papers (Gerstner et al. 1996) the term Spike Response Model was used
for a model where refractoriness and adaptation were modeled
by the combined effects of the
spike after potentials of * several* previous spikes,
rather than only the most recent spike.
Hence equation (1) is replaced by
\[\tag{5}
u(t) = \sum_{t^k} \eta(t-t^k) + \int \kappa(t-\hat{t},s) I^{\rm ext}(t-s) ds
\]

where \(t^k<t\) denotes previous moments of spike firing. Similarly, the threshold increases .

The reason for keeping in the normal Spike Response Model only the effect of the most recent spike is that a closed-form mathematical analysis of large networks in the form of mean-field equations (Gerstner 2000, Gerstner and van Hemmen 1992) is much easier in the standard Spike Response Model than in the cumulative Spike Response Model. The advantage of the cumulative model is that it accounts for adaptation and bursting (Gerstner and van Hemmen 1992, Gerstner et al. 1996).

### Resonate-and-fire model

The resonate-and-fire model (Izhikevich 2001, Richardson et al. 2003) is a special case of the Cumulative Spike Response Model. The resonate and fire model consists of two linear equations: \[ {du \over dt} = -{{v - v_{eq}} \over {\tau}} - c w + I(t) \] \[ \tau_w {dw \over dt} = a (v - v_{eq}) - w \] where \( w \) is a second variable that summarizes the effect of subthreshold membrane current. The voltage variable \( v \) is reset whenever it reaches a firing threshold. Integration of these linear equations gives a solution of the form (1) with a specific choice for the filters \(\eta\) and \(\kappa\ .\)

### Spike Response Model with a cumulative dynamic threshold

The threshold \(\theta\) is calculated as \[\tag{6} \theta(t) = \theta_0 + \sum_{t^k} \vartheta(t-t^k) \]

where \(t^k<t\) denotes previous moments of spike firing, \(\theta_0\) is the value of the threshold at rest, and \(\vartheta(t-t^k)\) describes the effect of a spike at time \(t^k\) on the value of the threshold at time \(t\ .\) The difference to the standard form of the spike response model is that now the value of the threshold depends on all previous spikes, not only the most recent one. Having a cumulative threshold is one possible way to incorporate adaptation. Furthermore, a spike response model with cumulative threshold makes it possible to have a single model with a fixed set of parameters that fits experimental data across a broad range of firing rates (Jolivet et al. 2006).

## Noise in the Spike Response Model

Noise can be included into the Spike Response Model by replacing the strict threshold criterion \(u(t) = \theta\) by a stochastic process. The probability \(P\) of firing a spike within a very short time \(\Delta t\) is \(P=\rho(t)\,\Delta t\) where the instantaneous firing rate or firing intensity \(\rho(t)\) is a function of the momentary difference between the membrane potential \(u(t)\) and the threshold \(\theta(t),\) \[\tag{SRMcum-eq:label exists!} \rho = f(u-\theta). \]

This noise model has been called stochastic threshold or escape noise (Gerstner and Kistler 2002). Typical functional forms for \(f\) could be an exponential function (Jolivet et al. 2006) or a rectified linear function. A saturating form for \(f\) is not a good choice, because it introduces a finite response time even for strong input well above threshold (Gerstner and Kistler 2002).

The advantage of the escape noise model in combination
with a Spike Response Model is that interspike interval
distributions can be expressed analytically for arbitrary
time-dependent input. Given a spike at time \(\hat{t}\)
and an input \(I(t)\ ,\) the next spike will occur at time \(t\)
with probability density
\[\tag{SRMcum-eq:label exists!}
P(t|\hat{t}) = \rho(t) \, \exp[ \int_\hat{t}^t -\rho(t')\, dt']
\]

From the point of view of stochastic processes, this is a time-dependent (inhomogeneous) version of a renewal model.

Because stochastic spike arrival with time-dependent input rates can be well represented by an equivalent escape noise model (Plesser and Gerstner 2000), one can calculate interspike interval distributions even for cases where the first-passage time problem of the leaky integrate-and-fire model with diffusive noise cannot be solved analytically. The Spike Response Model with escape noise has turned out to be useful for description of experimental data (Jolivet et al. 2006, Pillow et al. 2008).

## Effects not captured by a Spike Response Model

The Spike Response Model is not suited to describe the following effects:

- Pharmacological blocking of ion channels.

The biophysics of the neuronal membrane is not described explicitly in the Spike Response Model. Instead, the combined effects of several ion channel are captured phenomenologically in the spike shape function \(\eta\) and the filter \(\kappa\ .\) Therefore the model cannot make predictions about blocking of individual ion channels. A Hodgkin-Huxley model is better suited to describe the effects of individual channels.

- Delayed spike initiation.

For isolated input pulses, type I neuron models exhibit action potentials of standard shape, but different delay depending on the amplitude of the input pulse. The Spike Response Model cannot capture these effects, because of the strict threshold criterion. A quadratic (Latham et al. 2000) or an exponential integrate-and-fire model (Fourcaud-Trocme et al. 2003) is better suited to describe these effects.

- Dependence of the threshold upon the input.

The critical voltage for spike initiation with very slow ramp currents or constant currents (rheobase threshold) is different in real neurons from the critical voltage for spike initiation with short current pulses. The Spike Response Model cannot capture these effects because of the strict threshold criterion. An adaptive exponential integrate-and-fire (Brette and Gerstner, 2006) model is better suited to describe these effects.

## History

Simple spiking neurons with a formal firing threshold that potentially increases after each spike have a long history (e.g., Lapique 1907, Hill 1936, Weiss 1966). Models closely related to the Spike Response Model have been used for quantitative prediction of neural spike trains by several groups (Brillinger 1992, Keat et al. 2001, Jolivet et al. 2004, Carandini 2007, Pillow 2008). The term Spike Response Model was introduced around 1993 (Gerstner et al., 1993; Gerstner 1995).

## Summary

A major advantage of formal spiking neurons
such as the Spike Response Model
is their * simplicity* which has
several
important consequences.

- It is possible to simulate neural networks with a large number of neurons at a reasonable numerical cost.

- Network properties such as the mean firing rate of neurons in a network of connected spiking units can be studied analytically using tools from mathematical probability theory, statistical physics, and bifurcation theory.

- Questions of neural coding can be addressed in a transparent fashion.

- Simple spiking neuron models can be fitted to experimental data.

- Compared to standard leaky integrate-and-fire models the Spike Response Model allows to cover refractoriness.

- The spike response model is the most general model that combines linear filtering with a strict threshold. Hence it incorporates several spiking neuron models as special cases.

## References

- Badel L., Lefort S., Brette R., Petersen C., Gerstner W. and Richardson M.J.E. (2008), Dynamic I-V Curves Are Reliable Predictors of Naturalistic Pyramidal-Neuron Voltage Traces, J Neurophysiol 99: 656 - 666, 2008.

- Brillinger, D. R. (1992), Nerve cell spike train data analysis: a progression of technique. J. Am. Statist. Assoc., 87:260-271.

- Brette R. and Gerstner W. (2005), Adaptive Exponential Integrate-and-Fire Model as an Effective Description of Neuronal Activity, J. Neurophysiol., Vol. 94, pp. 3637 - 3642, 2005.

- Carandini M, Horton J.C. and Sincich L.M. (2007), Thalamic filtering of retinal spike trains by postsynaptic summation, Journal of Vision (2007) 7:20, 1–11

- Fourcaud-Trocme N., Hansel D., van Vreeswijk C., and Brunel N. (2003), How spike generation mechanisms determine the neuronal response to fluctuating inputs, J. Neuroscience 23:11628-11640

- Gerstner W. and van Hemmen J.L. (1992), Associative memory in a network of 'spiking' neurons. Network 3 :139-164

- Gerstner W, Ritz R, and van Hemmen JL (1993), Why spikes? Hebbian learning and retrieval of time-resolved excitation patterns. Biol. Cybern. 69:503-515

- Gerstner W., (1995), Time structure of the activity in neural network models. Phys. Rev. E, 51 :738-758

- Gerstner W., vanHemmen J.L., Cowan JD (1996), What matters in neuronal locking,Neural Computation 8:1653-1676

- Gerstner W. (2000), Population Dynamics of Spiking Neurons: Fast Transients, Asynchronous States, and Locking. Neural Computation 12:43-89.

- Gerstner and Kistler W.M. (2002),
*Spiking Neuron Models - Single Neurons, Populations, Plasticity*, Cambridge Univ. Press

- Hill, A. V. (1936), Excitation and accommodation in nerve. Proc. R. Soc. B, 119:305-355.

- Izhikevich E.M. (2003), Simple model of spiking neurons, IEEE Transactions On Neural Networks, 14:1569-1572

- Izhikevich, E.M. (2001), Resonate-and-fire neurons, Neural Networks, 14:883-894

- Jolivet R., Lewis T.J., and Gerstner W. (2004), Generalized Integrate-and-Fire Models of Neuronal Activity Approximate Spike Trains of a Detailed Model to a High Degree of Accuracy. J. Neurophysiology 92: 959-976

- Jolivet R., Rauch A., Luscher H.R., and Gerstner, W. (2006), Predicting spike timing of neocortical pyramidal neurons by simple threshold models. Journal of Computational Neuroscience 21:35-49

- Jolivet R., Kobayashi R., Rauch A., Naud R. , Shinomoto S., and Gerstner W (2008), A benchmark test for a quantitative assessment of simple neuron models. Journal of Neuroscience Methods 169: 417-424

- Keat J, Reinagel P., Reid R.C., and Meister M. (2001), Predicting every spike: A model for the responses of visual neurons, Neuron 30:803-817

- Lapicque, L. (1907), Recherches quantitatives sur l'excitation electrique des nerfs traitée comme une polarisation. J. Physiol. Pathol. Gen., 9:620-635.

- Latham, P.E., Richmond, B.J., Nelson, P.G., and Nirenberg, S. (2000), Intrinsic dynamics in neuronal networks. I. Theory. J. Neurophysiol., 83:808-827.

- Pillow, J.W., Shlens, J., Paninski, L., Sher, A., Litke, A.M., Chichilnisky, E.J., and Simoncelli, E.P. (2008) Spatio-temporal correlations and visual signaling in a complete neuronal population. Nature 454: 995-999

- Kistler W.M., Gerstner W., and vanHemmen J.L. (1997), Reduction of Hodgkin-Huxley equations to a threshold model. Neural Comput. 9:1069-1100

- Richardson M.J.E., Brunel N., and Hakim V (2003), From Subthreshold to Firing-Rate Resonance. Journal of Neurophysiology 89:2538-2554

- Weiss, T. F. (1966). A model of the peripheral auditory system. Kybernetik, 3:153-175.

**Internal references**

- Rodolfo Llinas (2008) Neuron. Scholarpedia, 3(8):1490.

## Recommended reading

- Gerstner W. and Kistler W.K. (2002) Spiking Neuron Models - single neurons, populations, plasticity. Cambridge University Press ISBN 0521890799

## External links

## See also

Spiking Neurons, Leaky Integrate-and-Fire Neuron, Integrate-and-fire model, Neuron, Renewal Model, Refractoriness